HIV-1 transmission cluster with T215D revertant mutation among newly diagnosed patients from the Basque Country, Spain

Transmission Clusters, Predominantly Associated With Men Who Have Sex With Men, Play a Main Role in the Propagation of HIV-1 in Northern Spain (2013–2018)

Viruses of HIV-1-infected individuals whose transmission is related group phylogenetically in transmission clusters (TCs). The study of the phylogenetic relations of these viruses and the factors associated with these individuals is essential to analyze the HIV-1 epidemic. In this study, we examine the role of TCs in the epidemiology of HIV-1 infection in Galicia and the Basque County, two regions of northern Spain. A total of 1,158 HIV-1-infected patients from both regions with new diagnoses (NDs) in 2013–2018 were included in the study. Partial HIV-1 pol sequences were analyzed phylogenetically by approximately maximum-likelihood with FastTree 2. In this analysis, 10,687 additional sequences from samples from HIV-1-infected individuals collected in Spain in 1999–2019 were also included to assign TC membership and to determine TCs’ sizes. TCs were defined as those which included viruses from ≥4 individuals, at least 50% of them Spaniards, and with ≥0.95 Shimodaira-Hasegawa-like node support in the phylogenetic tree. Factors associated to TCs were evaluated using odds ratios (OR) and their 95% CI. Fifty-one percent of NDs grouped in 162 TCs. Male patients (OR: 2.6; 95% CI: 1.5–4.7) and men having sex with men (MSM; OR: 2.1; 95% CI: 1.4–3.2) had higher odds of belonging to a TC compared to female and heterosexual patients, respectively. Individuals from Latin America (OR: 0.3; 95% CI: 0.2–0.4), North Africa (OR: 0.4; 95% CI: 0.2–1.0), and especially Sub-Saharan Africa (OR: 0.02; 95% CI: 0.003–0.2) were inversely associated to belonging to TCs compared to native Spaniards. Our results show that TCs are important components of the HIV-1 epidemics in the two Spanish regions studied, where transmission between MSM is predominant. The majority of migrants were infected with viruses not belonging to TCs that expand in Spain. Molecular epidemiology is essential to identify local peculiarities of HIV-1 propagation. The early detection of TCs and prevention of their expansion, implementing effective control measures, could reduce HIV-1 infections.

A cross-sectional study to characterize local HIV-1 dynamics in Washington, DC using next-generation sequencing

Washington, DC continues to experience a generalized HIV-1 epidemic. We characterized the local phylodynamics of HIV-1 in DC using next-generation sequencing (NGS) data. Viral samples from 68 participants from 2016 through 2017 were sequenced and paired with epidemiological data. Phylogenetic and network inferences, drug resistant mutations (DRMs), subtypes and HIV-1 diversity estimations were completed. Haplotypes were reconstructed to infer transmission clusters. Phylodynamic inferences based on the HIV-1 polymerase (pol) and envelope genes (env) were compared. Higher HIV-1 diversity (n.s.) was seen in men who have sex with men, heterosexual, and male participants in DC. 54.0% of the participants contained at least one DRM. The 40-49 year-olds showed the highest prevalence of DRMs (22.9%). Phylogenetic analysis of pol and env sequences grouped 31.9-33.8% of the participants into clusters. HIV-TRACE grouped 2.9-12.8% of participants when using consensus sequences and 9.0-64.2% when using haplotypes. NGS allowed us to characterize the local phylodynamics of HIV-1 in DC more broadly and accurately, given a better representation of its diversity and dynamics. Reconstructed haplotypes provided novel and deeper phylodynamic insights, which led to networks linking a higher number of participants. Our understanding of the HIV-1 epidemic was expanded with the powerful coupling of HIV-1 NGS data with epidemiological data.

Diverse Large HIV-1 Non-subtype B Clusters Are Spreading Among Men Who Have Sex With Men in Spain

In Western Europe, the HIV-1 epidemic among men who have sex with men (MSM) is dominated by subtype B. However, recently, other genetic forms have been reported to circulate in this population, as evidenced by their grouping in clusters predominantly comprising European individuals. Here we describe four large HIV-1 non-subtype B clusters spreading among MSM in Spain. Samples were collected in 9 regions. A pol fragment was amplified from plasma RNA or blood-extracted DNA. Phylogenetic analyses were performed via maximum likelihood, including database sequences of the same genetic forms as the identified clusters. Times and locations of the most recent common ancestors (MRCA) of clusters were estimated with a Bayesian method. Five large non-subtype B clusters associated with MSM were identified. The largest one, of F1 subtype, was reported previously. The other four were of CRF02_AG (CRF02_1; n = 115) and subtypes A1 (A1_1; n = 66), F1 (F1_3; n = 36), and C (C_7; n = 17). Most individuals belonging to them had been diagnosed of HIV-1 infection in the last 10 years. Each cluster comprised viruses from 3 to 8 Spanish regions and also comprised or was related to viruses from other countries: CRF02_1 comprised a Japanese subcluster and viruses from 8 other countries from Western Europe, Asia, and South America; A1_1 comprised viruses from Portugal, United Kingom, and United States, and was related to the A1 strain circulating in Greece, Albania and Cyprus; F1_3 was related to viruses from Romania; and C_7 comprised viruses from Portugal and was related to a virus from Mozambique. A subcluster within CRF02_1 was associated with heterosexual transmission. Near full-length genomes of each cluster were of uniform genetic form. Times of MRCAs of CRF02_1, A1_1, F1_3, and C_7 were estimated around 1986, 1989, 2013, and 1983, respectively. MRCA locations for CRF02_1 and A1_1 were uncertain (however initial expansions in Spain in Madrid and Vigo, respectively, were estimated) and were most probable in Bilbao, Spain, for F1_3 and Portugal for C_7. These results show that the HIV-1 epidemic among MSM in Spain is becoming increasingly diverse through the expansion of diverse non-subtype B clusters, comprising or related to viruses circulating in other countries.

Changes from 2000 to 2009 in the Prevalence of HIV-1 Containing Drug Resistance-Associated Mutations from Antiretroviral Therapy-Naive, HIV-1-Infected Patients in the United States

Pre-existing HIV drug resistance can jeopardize first-line antiretroviral therapy (ART) success. Changes in the prevalence of drug resistance-associated mutations (DRMs) were analyzed from HIV-infected, ART-naive, U.S. individuals seeking ART treatment from 2000 to 2009. HIV DRM data from 3,829 ART-naive subjects were analyzed by year of sample collection using International Antiviral Society-United States (IAS-USA) and World Health Organization (WHO) “surveillance” DRM definitions; minor IAS-USA-defined DRMs were excluded. IAS-USA DRM prevalence between 2000 and 2009 was 14%, beginning with 8% in 2000 and 13% in 2009. The greatest incidence was observed in 2007 (17%). Overall, IAS-USA-defined non-nucleoside reverse transcriptase inhibitor (NNRTI) DRMs were 9.5%; nucleoside reverse transcriptase inhibitor (NRTI): 4%, and major protease inhibitor (PI): 3%. The most frequently detected IAS-USA-defined DRMs by class were NNRTI: K103N/S (4%), NRTI: M41L (1.5%), and PI: L90M (1%). Overall, WHO-defined DRM prevalence was 13% (5% in 2000; 13% in 2009). By class, NNRTI prevalence was 6%, NRTI: 6%, and PI: 3.2%. The most frequent WHO-defined DRMs were NRTI: codon T215 (3.0%), NNRTI: K103N/S (4%), and PI: L90 (1%). WHO-defined NNRTI DRMs declined significantly (p = .0412) from 2007 to 2009. The overall prevalence of HIV-1 containing major IAS-USA or WHO-defined DRMs to ≥2 or ≥3 classes was 2% and <1%, respectively. The prevalence of HIV-1 with WHO-defined dual- or triple-class resistance significantly declined (p = .0461) from 2008 (4%) to 2009 (<1%). In this U.S. cohort, the prevalence of HIV-1 DRMs increased from 2000 onward, peaked between 2005 and 2007, and then declined between 2008 and 2009; the detection of WHO-defined dual- or triple-class DRM similarly decreased from 2008 to 2009.

High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain: A phylogeographic approach

Phylogenetic studies are a valuable tool to understand viral transmission patterns and the role of immigration in HIV-1 spread. We analyzed the spatio-temporal relationship of different HIV-1 non-B subtype variants over time using phylogenetic analysis techniques. We collected 693 pol (PR+RT) sequences that were sampled from 2005 to 2012 from naïve patients in different hospitals in southern Spain. We used REGA v3.0 to classify them into subtypes and recombinant forms, which were confirmed by phylogenetic analysis through maximum likelihood (ML) using RAxML. For the main HIV-1 non-B variants, publicly available, genetically similar sequences were sought using HIV-BLAST. The presence of HIV-1 lineages circulating in our study population was established using ML and Bayesian inference (BEAST v1.7.5) and transmission networks were identified. We detected 165 (23.4%) patients infected with HIV-1 non-B variants: 104 (63%) with recombinant viruses in pol: CRF02_AG (71, 43%), CRF14_BG (8, 4.8%), CRF06_cpx (5, 3%) and nine other recombinant forms (11, 6.7%) and unique recombinants (9, 5.5%). The rest (61, 37%) were infected with non-recombinant subtypes: A1 (30, 18.2%), C (7, [4.2%]), D (3, [1.8%]), F1 (9, 5.5%) and G (12, 7.3%). Most patients infected with HIV-1 non-B variants were men (63%, p < 0.001) aged over 35 (73.5%, p < 0.001), heterosexuals (92.2%, p < 0.001), from Africa (59.5%, p < 0.001) and living in the El Ejido area (62.4%, p<0.001). We found lineages of epidemiological relevance (mainly within Subtype A1), imported primarily through female sex workers from East Europe. We detected 11 transmission clusters of HIV-1 non-B Subtypes, which included patients born in Spain in half of them. We present the phylogenetic profiles of the HIV-1 non-B variants detected in southern Spain, and explore their putative geographical origins. Our data reveals a high HIV-1 genetic diversity likely due to the import of viral lineages that circulate in other countries. The highly immigrated El Ejido area acts as a gateway through which different subtypes are introduced into other regions, hence the importance of setting up epidemiological control measures to prevent future outbreaks.

Characterization of HIV diversity, phylodynamics and drug resistance in Washington, DC

BACKGROUND

Washington DC has a high burden of HIV with a 2.0% HIV prevalence. The city is a national and international hub potentially containing a broad diversity of HIV variants; yet few sequences from DC are available on GenBank to assess the evolutionary history of HIV in the US capital. Towards this general goal, here we analyze extensive sequence data and investigate HIV diversity, phylodynamics, and drug resistant mutations (DRM) in DC.

METHODS

Molecular HIV-1 sequences were collected from participants infected through 2015 as part of the DC Cohort, a longitudinal observational study of HIV+ patients receiving care at 13 DC clinics. Sequences were paired with Cohort demographic, risk, and clinical data and analyzed using maximum likelihood, Bayesian and coalescent approaches of phylogenetic, network and population genetic inference. We analyzed 601 sequences from 223 participants for int (~864 bp) and 2,810 sequences from 1,659 participants for PR/RT (~1497 bp).

RESULTS

Ninety-nine and 94% of the int and PR/RT sequences, respectively, were identified as subtype B, with 14 non-B subtypes also detected. Phylodynamic analyses of US born infected individuals showed that HIV population size varied little over time with no significant decline in diversity. Phylogenetic analyses grouped 13.5% of the int sequences into 14 clusters of 2 or 3 sequences, and 39.0% of the PR/RT sequences into 203 clusters of 2-32 sequences. Network analyses grouped 3.6% of the int sequences into 4 clusters of 2 sequences, and 10.6% of the PR/RT sequences into 76 clusters of 2-7 sequences. All network clusters were detected in our phylogenetic analyses. Higher proportions of clustered sequences were found in zip codes where HIV prevalence is highest (r = 0.607; P<0.00001). We detected a high prevalence of DRM for both int (17.1%) and PR/RT (39.1%), but only 8 int and 12 PR/RT amino acids were identified as under adaptive selection. We observed a significant (P<0.0001) association between main risk factors (men who have sex with men and heterosexuals) and genotypes in the five well-supported clusters with sufficient sample size for testing.

DISCUSSION

Pairing molecular data with clinical and demographic data provided novel insights into HIV population dynamics in Washington, DC. Identification of populations and geographic locations where clustering occurs can inform and complement active surveillance efforts to interrupt HIV transmission.

The molecular epidemiology of HIV-1 in the Comunidad Valenciana (Spain): analysis of transmission clusters

HIV infections are still a very serious concern for public heath worldwide. We have applied molecular evolution methods to study the HIV-1 epidemics in the Comunidad Valenciana (CV, Spain) from a public health surveillance perspective. For this, we analysed 1804 HIV-1 sequences comprising protease and reverse transcriptase (PR/RT) coding regions, sampled between 2004 and 2014. These sequences were subtyped and subjected to phylogenetic analyses in order to detect transmission clusters. In addition, univariate and multinomial comparisons were performed to detect epidemiological differences between HIV-1 subtypes, and risk groups. The HIV epidemic in the CV is dominated by subtype B infections among local men who have sex with men (MSM). 270 transmission clusters were identified (>57% of the dataset), 12 of which included ≥10 patients; 11 of subtype B (9 affecting MSMs) and one (n = 21) of CRF14, affecting predominately intravenous drug users (IDUs). Dated phylogenies revealed these large clusters to have originated from the mid-80s to the early 00 s. Subtype B is more likely to form transmission clusters than non-B variants and MSMs to cluster than other risk groups. Multinomial analyses revealed an association between non-B variants, which are not established in the local population yet, and different foreign groups.

Decreasing prevalence of transmitted drug resistance among ART-naive HIV-1-infected patients in Iceland, 1996-2012

Introduction: Resistance to antiretroviral drugs can complicate the management of HIV-1 infection and impair control of its spread. The aim of the current study was to investigate the prevalence and transmission of HIV-1 drug resistance among 106 antiretroviral therapy (ART)-naïve patients diagnosed in Iceland (1996–2012).

Methods: HIV-1 polymerase sequences were analysed using the Calibrated Population Resistance tool. Domestic spread of transmitted drug resistance (TDR) was investigated through maximum likelihood and Bayesian approaches.

Results: Among ART-naïve patients, the prevalence of TDR to any of the following classes (NRTIs, NNRTIs and PIs) was 8.5% (95% CI: 4.5%- 15.4%): 6.6% to NRTIs, 0.9% to NNRTIs, and 1.9% to PIs. The most frequent NRTI mutation detected was T215C/D (n=7, 5.7%). The only NNRTI mutation detected was K103N (n=1, 0.9%). PI mutations detected were M46I (n=1, 0.9%) and L90M (n=1, 0.9%). Six patients harbouring T215C/D, were linked in a supported phylogenetic cluster. No significant association was found between TDR and demographic or risk groups. Trend analysis showed a decrease in the prevalence of TDR (1996–2012, p=0.003).

Conclusions: TDR prevalence in Iceland was at a moderate level and decreased during 1996-2012. Screening for TDR is recommended to limit its local spread and to optimize HIV-1 therapy.

Abbreviations: ART: Anti-retroviral therapy; ARV: antiretroviral; ATV/r: atazanavir/ritonavir; AZT: azidothymidine; BEAST: Bayesian evolutionary analysis by sampling trees; CI: confidence interval; CPR: calibrated population resistance; CRF: circulating recombinant form; d4T: stavudine; EFV: efavirenz; FET: Fishers’ exact test; FPV/r: fosamprenavir/ritonavir; HET: heterosexual; IDU: injection drug use; IDV/r: indinavir/ritonavir; LPV/r: lopinavir/ritonavir; MSM: men who have sex with men; M-W: Mann–Whitney U test; NFV: nelfinavir; NNRTIs: non-nucleoside reverse transcriptase inhibitors; NRTIs: nucleoside reverse transcriptase inhibitors; NVP: nevirapine; PIs: protease inhibitors; pol: polymerase gene; SDRM: surveillance drug resistance mutation; SQV/r: saquinavir/ritonavir; TDR: transmitted drug resistance

Defining HIV-1 transmission clusters based on sequence data

Understanding HIV-1 transmission dynamics is relevant to both screening and intervention strategies of HIV-1 infection. Commonly, HIV-1 transmission chains are determined based on sequence similarity assessed either directly from a sequence alignment or by inferring a phylogenetic tree. This review is aimed at both nonexperts interested in understanding and interpreting studies of HIV-1 transmission, and experts interested in finding the most appropriate cluster definition for a specific dataset and research question. We start by introducing the concepts and methodologies of how HIV-1 transmission clusters usually have been defined. We then present the results of a systematic review of 105 HIV-1 molecular epidemiology studies summarizing the most common methods and definitions in the literature. Finally, we offer our perspectives on how HIV-1 transmission clusters can be defined and provide some guidance based on examples from real life datasets.

Phylodynamic and Phylogeographic Profiles of Subtype B HIV-1 Epidemics in South Spain

BACKGROUND

Since 1982, HIV-1 epidemics have evolved to different scenarios in terms of transmission routes, subtype distribution and characteristics of transmission clusters. We investigated the evolutionary history of HIV-1 subtype B in south Spain.

PATIENTS & METHODS

We studied all newly diagnosed HIV-1 subtype B patients in East Andalusia during the 2005-2012 period. For the analysis, we used the reverse transcriptase and protease sequences from baseline resistance, and the Trugene® HIV Genotyping kit (Siemens, Barcelona, Spain). Subtyping was done with REGA v3.0. The maximum likelihood trees constructed with RAxML were used to study HIV-1 clustering. Phylogeographic and phylodynamic profiles were studied by Bayesian inference methods with BEAST v1.7.5 and SPREAD v1.0.6.

RESULTS

Of the 493 patients infected with HIV-1 subtype B, 234 grouped into 55 clusters, most of which were small (44 clusters ≤ 5 patients, 31 with 2 patients, 13 with 3). The rest (133/234) were grouped into 11 clusters with ≥ 5 patients, and most (82%, 109/133) were men who have sex with men (MSM) grouped into 8 clusters. The association with clusters was more frequent in Spanish (p = 0.02) men (p< 0.001), MSM (p<0.001) younger than 35 years (p = 0.001) and with a CD4+ T-cell count above 350 cells/ul (p<0.001). We estimated the date of HIV-1 subtype B regional epidemic diversification around 1970 (95% CI: 1965-1987), with an evolutionary rate of 2.4 (95%CI: 1.7-3.1) x 10-3 substitutions/site/year. Most clusters originated in the 1990s in MSMs. We observed exponential subtype B HIV-1 growth in 1980-1990 and 2005-2008. The most significant migration routes for subtype B went from inland cities to seaside locations.

CONCLUSIONS

We provide the first data on the phylodynamic and phylogeographic profiles of HIV-1 subtype B in south Spain. Our findings of transmission clustering among MSMs should alert healthcare managers to enhance preventive measures in this risk group in order to prevent future outbreaks.

Reconstructing contact network parameters from viral phylogenies

Models of the spread of disease in a population often make the simplifying assumption that the population is homogeneously mixed, or is divided into homogeneously mixed compartments. However, human populations have complex structures formed by social contacts, which can have a significant influence on the rate of epidemic spread. Contact network models capture this structure by explicitly representing each contact which could possibly lead to a transmission. We developed a method based on approximate Bayesian computation (ABC), a likelihood-free inference strategy, for estimating structural parameters of the contact network underlying an observed viral phylogeny. The method combines adaptive sequential Monte Carlo for ABC, Gillespie simulation for propagating epidemics though networks, and a kernel-based tree similarity score. We used the method to fit the Barabási-Albert network model to simulated transmission trees, and also applied it to viral phylogenies estimated from ten published HIV sequence datasets. This model incorporates a feature called preferential attachment (PA), whereby individuals with more existing contacts accumulate new contacts at a higher rate. On simulated data, we found that the strength of PA and the number of infected nodes in the network can often be accurately estimated. On the other hand, the mean degree of the network, as well as the total number of nodes, was not estimable with ABC. We observed sub-linear PA power in all datasets, as well as higher PA power in networks of injection drug users. These results underscore the importance of considering contact structures when performing phylodynamic inference. Our method offers the potential to quantitatively investigate the contact network structure underlying viral epidemics.

Transmission dynamics of HIV-1 subtype B in the Basque Country, Spain

This work was aimed to study the HIV-1 subtype B epidemics in the Basque Country, Spain. 1727 HIV-1 subtype B sequences comprising protease and reverse transcriptase (PR/RT) coding regions, sampled between 2001 and 2008, were analyzed. 156 transmission clusters were detected by means of phylogenetic analyses. Most of them comprised less than 4 individuals and, in total, they included 441 patients. Six clusters comprised 10 or more patients and were further analyzed in order to study their origin and diversification. Four clusters included men who had unprotected homosexual sex (MSM), one group was formed by intravenous drug users (IDUs), and another included both IDUs and people infected through unprotected heterosexual sex (HTs). Most of these clusters originated from the mid-1980s to the mid-1990s. Only one cluster, formed by MSM, originated after 2000. The time between infections was significantly lower in MSM groups than in those containing IDUs (P-value <0.0001). Nucleoside RT and non-nucleoside RT inhibitor (NRTI and NNRTI)-resistance mutations to antiretroviral treatment were found in these six clusters except the most recent MSM group, but only the IDU clusters presented protease inhibitor (PI)-resistance mutations. The most prevalent mutations for each inhibitor class were PI L90M, NRTI T215D/Y/F, and NNRTI K103N, which were also among the most prevalent resistant variants in the whole dataset. In conclusion, while most infections occur as isolated introductions into the population, the number of infections found to be epidemiologically related within the Basque Country is significant. Public health control measures should be reinforced to prevent the further expansion of transmission clusters and resistant mutations occurring within them.

Phylogeny and Phylogeography of a Recent HIV-1 Subtype F Outbreak among Men Who Have Sex with Men in Spain Deriving from a Cluster with a Wide Geographic Circulation in Western Europe

We recently reported the rapid expansion of an HIV-1 subtype F cluster among men who have sex with men (MSM) in the region of Galicia, Northwest Spain. Here we update this outbreak, analyze near full-length genomes, determine phylogenetic relationships, and estimate its origin. For this study, we used sequences of HIV-1 protease-reverse transcriptase and env V3 region, and for 17 samples, near full-length genome sequences were obtained. Phylogenetic analyses were performed via maximum likelihood. Locations and times of most recent common ancestors were estimated using Bayesian inference. Among samples analyzed by us, 100 HIV-1 F1 subsubtype infections of monophyletic origin were diagnosed in Spain, including 88 in Galicia and 12 in four other regions. Most viruses (n = 90) grouped in a subcluster (Galician subcluster), while 7 from Valladolid (Central Spain) grouped in another subcluster. At least 94 individuals were sexually-infected males and at least 71 were MSM. Seventeen near full-length genomes were uniformly of F1 subsubtype. Through similarity searches and phylogenetic analyses, we identified 18 viruses from four other Western European countries [Switzerland (n = 8), Belgium (n = 5), France (n = 3), and United Kingdom (n = 2)] and one from Brazil, from samples collected in 2005–2011, which branched within the subtype F cluster, outside of both Spanish subclusters, most of them corresponding to recently infected individuals. The most probable geographic origin and age of the Galician subcluster was Ferrol, Northwest Galicia, around 2007, while the Western European cluster probably emerged in Switzerland around 2002. In conclusion, a recently expanded HIV-1 subtype F cluster, the largest non-subtype B cluster reported in Western Europe, continues to spread among MSM in Spain; this cluster is part of a larger cluster with a wide geographic circulation in diverse Western European countries.

Clinical, virological and phylogenetic characterization of a multiresistant HIV-1 strain outbreak in naive patients in southern Spain

BACKGROUND

We describe the characteristics of an HIV-1 strain with six viral reverse transcriptase mutations (D67N, T69N/D, V118I, V179D, T215S and K219Q), which we have called the Malaga strain. This strain was detected in treatment-naive patients from southern Spain.

METHODS

The study was undertaken at the Virgen de la Victoria Hospital, Malaga, a reference centre for the study of HIV-1 genotype resistance in Andalusia (the 'Costa del Sol'), Spain. Genotypic resistance testing was done in an automated sequencer. Phylogenetic analysis was performed using a 630 bp region of the reverse transcriptase with the mutations mentioned.

RESULTS

Between 2007 and 2014, we detected the Malaga strain in 30 treatment-naive patients. All were MSM, seen at five hospitals on the Costa del Sol. In all cases, the HIV-1 was subtype B with viral tropism R5. Phylogenetic analysis based on the reverse transcriptase sequence showed consistent grouping (with a bootstrap value of the common node of 100%) of the isolates that shared the mutation pattern mentioned. This strain has not been detected elsewhere or in previously treated patients. All of the patients treated with first-line combination ART responded.

CONCLUSIONS

We report a cluster of an HIV-1 strain with multiple resistance mutations that was transmitted over a period of >8 years, affecting 30 naive patients from the same geographical area. The strain was susceptible to first-line combination ART.

Epidemiological Surveillance of HIV-1 Transmitted Drug Resistance in Spain in 2004-2012: Relevance of Transmission Clusters in the Propagation of Resistance Mutations

Our objectives were to carry out an epidemiological surveillance study on transmitted drug resistance (TDR) among individuals newly diagnosed of HIV-1 infection during a nine year period in Spain and to assess the role of transmission clusters (TC) in the propagation of resistant strains. An overall of 1614 newly diagnosed individuals were included in the study from January 2004 through December 2012. Individuals come from two different Spanish regions: Galicia and the Basque Country. Resistance mutations to reverse transcriptase inhibitors (RTI) and protease inhibitors (PI) were analyzed according to mutations included in the surveillance drug-resistance mutations list updated in 2009. TC were defined as those comprising viruses from five or more individuals whose sequences clustered in maximum likelihood phylogenetic trees with a bootstrap value ≥90%. The overall prevalence of TDR to any drug was 9.9%: 4.9% to nucleoside RTIs (NRTIs), 3.6% to non-nucleoside RTIs (NNRTIs), and 2.7% to PIs. A significant decrease of TDR to NRTIs over time was observed [from 10% in 2004 to 2% in 2012 (p=0.01)]. Sixty eight (42.2%) of 161 sequences with TDR were included in 25 TC composed of 5 or more individuals. Of them, 9 clusters harbored TDR associated with high level resistance to antiretroviral drugs. T215D revertant mutation was transmitted in a large cluster comprising 25 individuals. The impact of epidemiological networks on TDR frequency may explain its persistence in newly diagnosed individuals. The knowledge of the populations involved in TC would facilitate the design of prevention programs and public health interventions.

HIV-1 transmission during early infection in men who have sex with men: a phylodynamic analysis

BACKGROUND

Conventional epidemiological surveillance of infectious diseases is focused on characterization of incident infections and estimation of the number of prevalent infections. Advances in methods for the analysis of the population-level genetic variation of viruses can potentially provide information about donors, not just recipients, of infection. Genetic sequences from many viruses are increasingly abundant, especially HIV, which is routinely sequenced for surveillance of drug resistance mutations. We conducted a phylodynamic analysis of HIV genetic sequence data and surveillance data from a US population of men who have sex with men (MSM) and estimated incidence and transmission rates by stage of infection.

METHODS AND FINDINGS

We analyzed 662 HIV-1 subtype B sequences collected between October 14, 2004, and February 24, 2012, from MSM in the Detroit metropolitan area, Michigan. These sequences were cross-referenced with a database of 30,200 patients diagnosed with HIV infection in the state of Michigan, which includes clinical information that is informative about the recency of infection at the time of diagnosis. These data were analyzed using recently developed population genetic methods that have enabled the estimation of transmission rates from the population-level genetic diversity of the virus. We found that genetic data are highly informative about HIV donors in ways that standard surveillance data are not. Genetic data are especially informative about the stage of infection of donors at the point of transmission. We estimate that 44.7% (95% CI, 42.2%-46.4%) of transmissions occur during the first year of infection.

CONCLUSIONS

In this study, almost half of transmissions occurred within the first year of HIV infection in MSM. Our conclusions may be sensitive to un-modeled intra-host evolutionary dynamics, un-modeled sexual risk behavior, and uncertainty in the stage of infected hosts at the time of sampling. The intensity of transmission during early infection may have significance for public health interventions based on early treatment of newly diagnosed individuals.

The emerging and re-emerging human immunodeficiency virus epidemics in Europe

We provide a summary of the current status of the human immunodeficiency virus (HIV) epidemic, including data on the transmission of drug-resistant virus in the European Region of the WHO. The review was conducted by searching the reports of the European Centre for Disease Prevention and Control and the United Nations General Assembly Special Session country reports to identify the number of HIV cases reported in 2002-2011, the number of HIV tests performed, and the results of the most recent HIV surveys in at-risk groups. In the West, a 5-year linear trend analysis suggests an increase in the number of reported HIV cases in Belgium, Greece, and Iceland, and a decline in The Netherlands, Ireland, Portugal, Switzerland, and the UK. In the Centre, the number of reported cases increased in Bulgaria and Turkey, and in the East, the number of reported cases increased in Armenia, Georgia, and Ukraine. Estonia and Latvia reported a significant downward trend. HIV transmission in men who have sex with men (MSM) and heterosexual transmission accounted for, respectively, 40% and 38% of newly reported HIV cases in the West in 2011, whereas the epidemic in the Centre is predominantly concentrated in MSM. Although injection drug use remains the key driving force of the epidemic in the East, there is increasing evidence of sexual transmission from the core groups of injection drug users, and increasing spread within MSM. The pattern of transmitted drug resistance (TDR) is heterogeneous accross Europe; a significant number of clusters of TDR have been recently observed in several European countries.

Molecular typing of the local HIV-1 epidemic in Serbia

Worldwide HIV-1 pandemic is becoming increasingly complex, with growing heterogeneity of subtypes and recombinant viruses. Previous studies have documented HIV-1 subtype B as the predominant one in Serbia, with limited presence and genetic diversity of non B subtypes. In recent years, MSM transmission has become the most frequently reported risk for HIV infection among newly diagnosed patients in Serbia, but very little is known of the network structure and dynamics of viral transmission in this and other risk groups. To gain insight about the HIV-1 subtypes distribution pattern as well as characteristics of HIV-1 transmission clusters in Serbia, we analyzed the genetic diversity of the pol gene segment in 221 HIV-1-infected patients sampled during 2002-2011. Subtype B was found to still be the most prevalent one in Serbia, accounting for over 90% of samples, while greater diversity of other subtypes was found than previously reported, including subtypes G, C, A, F, CRF01 and CRF02. In total, 41.3% of analyzed subtype B sequences were found associated in transmission clusters/network, that are highly related with MSM transmission route.

Simple epidemiological dynamics explain phylogenetic clustering of HIV from patients with recent infection

Phylogenies of highly genetically variable viruses such as HIV-1 are potentially informative of epidemiological dynamics. Several studies have demonstrated the presence of clusters of highly related HIV-1 sequences, particularly among recently HIV-infected individuals, which have been used to argue for a high transmission rate during acute infection. Using a large set of HIV-1 subtype B pol sequences collected from men who have sex with men, we demonstrate that virus from recent infections tend to be phylogenetically clustered at a greater rate than virus from patients with chronic infection ('excess clustering') and also tend to cluster with other recent HIV infections rather than chronic, established infections ('excess co-clustering'), consistent with previous reports. To determine the role that a higher infectivity during acute infection may play in excess clustering and co-clustering, we developed a simple model of HIV infection that incorporates an early period of intensified transmission, and explicitly considers the dynamics of phylogenetic clusters alongside the dynamics of acute and chronic infected cases. We explored the potential for clustering statistics to be used for inference of acute stage transmission rates and found that no single statistic explains very much variance in parameters controlling acute stage transmission rates. We demonstrate that high transmission rates during the acute stage is not the main cause of excess clustering of virus from patients with early/acute infection compared to chronic infection, which may simply reflect the shorter time since transmission in acute infection. Higher transmission during acute infection can result in excess co-clustering of sequences, while the extent of clustering observed is most sensitive to the fraction of infections sampled.

HIV-1 subtype D infections among Caucasians from Northwestern Poland--phylogenetic and clinical analysis

BACKGROUND

HIV-1 subtype D infections, which are associated with a faster rate of progression and lymphocyte CD4 decline, cognitive deficit and higher mortality, have rarely been found in native Europeans. In Northwestern Poland, however, infections with this subtype had been identified. This study aimed to analyze the sequence and clinical data for patients with subtype D using molecular phylogeography and identify transmission clusters and ancestry, as well as drug resistance, baseline HIV tropism and antiretroviral treatment efficacy.

METHODS

Phylogenetic analyses of local HIV-1 subtype D sequences were performed, with time to the most recent common ancestor inferred using bayesian modeling. Sequence and drug resistance data were linked with the clinical and epidemiological information.

RESULTS

Subtype D was found in 24 non-immigrant Caucasian, heterosexually infected patients (75% of females, median age at diagnosis of 49.5 years; IQR: 29-56 years). Partial pol sequences clustered monophyletically with the clades of Ugandan origin and no evidence of transmission from other European countries was found. Time to the most common recent ancestor was 1989.24 (95% HPD: 1968.83-1994.46). Baseline drug resistance to nucleoside reverse transcriptase inhibitors was observed in 54.5% of cases (mutations: M41L, K103N, T215S/D) with evidence of clustering, no baseline integrase or protease resistance and infrequent non-R5 tropism (13.6%). Virologic failure was observed in 60% of cases and was associated with poor adherence (p<0.001) and subsequent development of drug resistance (p = 0.008, OR: 20 (95%CI: 1.7-290).

CONCLUSIONS

Local subtype D represented an independently transmitted network with probably single index case, high frequency of primary drug resistance and evidence of transmission clusters.

Molecular phylogenetics of transmitted drug resistance in newly diagnosed HIV Type 1 individuals in Denmark: a nation-wide study

Highly active antiretroviral treatment is compromised by viral resistance mutations. Transmitted drug resistance (TDR) is therefore monitored closely, but follow-up studies of these patients are limited. Virus from 1405 individuals diagnosed with HIV-1 in Denmark between 2001 and 2009 was analyzed for TDR, and molecular-epidemiological links and progression of the infection were described based on data from standardized questionnaires, the prospective Danish HIV Cohort Study, and by phylogenetic analysis. Eighty-five individuals were found to be infected with virus harboring mutations resulting in a prevalence of 6.1%, with no changes over time. The main resistance mutations were nucleoside reverse transcriptase inhibitor (NRTI) mutation 215 revertants, as well as nonnucleoside reverse transcriptase inhibitor (NNRTI) mutation 103N/S and protease inhibitor (PI) mutations 90M and 85V. Phylogenetic analysis confirmed 12 transmission chains involving 37 TDR individuals. Of these 21 were also documented epidemiologically. The virus included in the transmission chain carried similar resistance mutations to the TDR index case, whereas controls chains from index cases without TDR were generally without resistance mutations. We observed no difference in progression of the infection between individuals infected with TDR and individuals infected with wild-type HIV-1. The prevalence of TDR is low in Denmark and transmission of dual-drug-resistant HIV-1 is infrequent. The TDR isolates were shown to originate from local patients failing therapy.

Molecular and epidemiological characterization of HIV-1 infection networks involving transmitted drug resistance mutations in Northern Greece

OBJECTIVES

To determine the contribution of transmission clusters to transmitted drug resistance (TDR) in newly diagnosed antiretroviral-naive HIV-1-infected patients in Northern Greece during 2000-07.

METHODS

The prevalence of TDR was estimated in 369 individuals who were diagnosed with HIV-1 infection in the period 2000-07 at the National AIDS Reference Laboratory of Northern Greece. Phylogenetic analysis was performed using a maximum likelihood method on partial pol sequences. TDR was defined in accordance with the surveillance drug resistance mutation list (2009 update).

RESULTS

The overall prevalence of TDR in our population was 12.5% [46/369, 95% confidence interval (CI) 9.1%-15.8%], comprising 7.6% (28/369) resistant to nucleoside reverse transcriptase inhibitors, 5.4% (20/369) resistant to non-nucleoside reverse transcriptase inhibitors and 3.3% (12/369) resistant to protease inhibitors. Dual class resistance was identified in 3.8% (14/369). Infection with subtype A was the sole predictor associated with TDR in multivariate analysis (odds ratio 2.15, 95% CI 1.10-4.19, P = 0.025). Phylogenetic analyses revealed three statistically robust transmission clusters involving drug-resistant strains, including one cluster of 12 patients, 10 of whom were infected with a strain carrying both T215 revertants and Y181C mutations.

CONCLUSIONS

Our findings underline the substantial impact of transmission networks on TDR in our population.

Phylogenetic analysis of full-length pol gene from Korean hemophiliacs and plasma donors infected with Korean subclade B of HIV type 1

There was an outbreak of HIV-1 transmission among 20 out of 122 Korean hemophiliacs from 1990. We assessed the genetic relationships among HIV-1 viruses found in three cash-paid plasma donors whose preseroconversion plasma was used to produce Korean-made clotting factor, 20 hemophiliacs infected with HIV-1 in Korea, three hemophiliacs infected with HIV-1 from clotting factor manufactured outside Korea, and 71 local control patients infected with the Korean subclade of HIV-1 subtype B (KSB). Full-length pol gene sequences (2841 bp) of viruses from frozen stored serum, samples obtained 1-3 years after diagnosis, were amplified by RT-PCR and sequenced by direct DNA sequencing. Phylogenetic and signature pattern analyses were used to investigate the relationships among the sequences. Donors O and P were associated with two clusters, of 8 and 12 hemophiliacs, respectively, which were demarcated from the 71 KSB-infected local control patients and donor R. These data confirm that HIV-1 transmission to 20 hemophiliacs occurred through infusion of Korean-made clotting factor.

Phylodynamics of HIV-1 from a phase III AIDS vaccine trial in Bangkok, Thailand

BACKGROUND

In 2003, a phase III placebo-controlled trial (VAX003) was completed in Bangkok, Thailand. Of the 2,546 individuals enrolled in the trial based on high risk for infection through injection drug use (IDU), we obtained clinical samples and HIV-1 sequence data (envelope glycoprotein gene gp120) from 215 individuals who became infected during the trial. Here, we used these data in combination with other publicly available gp120 sequences to perform a molecular surveillance and phylodynamic analysis of HIV-1 in Thailand.

METHODOLOGY AND FINDINGS

Phylogenetic and population genetic estimators were used to assess HIV-1 gp120 diversity as a function of vaccination treatment, viral load (VL) and CD4(+) counts, to identify transmission clusters and to investigate the timescale and demographics of HIV-1 in Thailand. Three HIV-1 subtypes were identified: CRF01_AE (85% of the infections), subtype B (13%) and CRF15_AE (2%). The Bangkok IDU cohort showed more gp120 diversity than other Asian IDU cohorts and similar diversity to that observed in sexually infected individuals. Moreover, significant differences (P<0.02) in genetic diversity were observed in CRF01_AE IDU with different VL and CD4(+) counts. No phylogenetic structure was detected regarding any of the epidemiological and clinical factors tested, although high proportions (35% to 50%) of early infections fell into clusters, which suggests that transmission chains associated with acute infection play a key role on HIV-1 spread among IDU. CRF01_AE was estimated to have emerged in Thailand in 1984.5 (1983-1986), 3-6 years before the first recognition of symptomatic patients (1989). The relative genetic diversity of the HIV-1 population has remained high despite decreasing prevalence rates since the mid 1990s.

CONCLUSIONS

Our study and recent epidemiological reports indicate that HIV-1 is still a major threat in Thailand and suggest that HIV awareness and prevention needs to be strengthened to avoid AIDS resurgence.

Prevalence of transmitted drug resistance and impact of transmitted resistance on treatment success in the German HIV-1 Seroconverter Cohort

BACKGROUND

The aim of this study is to analyse the prevalence of transmitted drug resistance, TDR, and the impact of TDR on treatment success in the German HIV-1 Seroconverter Cohort.

METHODS

Genotypic resistance analysis was performed in treatment-naïve study patients whose sample was available 1,312/1,564 (83.9% October 2008). A genotypic resistance result was obtained for 1,276/1,312 (97.3%). The resistance associated mutations were identified according to the surveillance drug resistance mutations list recommended for drug-naïve patients. Treatment success was determined as viral suppression below 500 copies/ml.

RESULTS

Prevalence of TDR was stable at a high level between 1996 and 2007 in the German HIV-1 Seroconverter Cohort (N = 158/1,276; 12.4%; CI(wilson) 10.7-14.3; p(for trend) = 0.25). NRTI resistance was predominant (7.5%) but decreased significantly over time (CI(Wilson): 6.2-9.1, p(for trend) = 0.02). NNRTI resistance tended to increase over time (NNRTI: 3.5%; CI(Wilson): 2.6-4.6; p(for trend)= 0.07), whereas PI resistance remained stable (PI: 3.0%; CI(Wilson): 2.1-4.0; p(for trend) = 0.24). Resistance to all drug classes was frequently caused by singleton resistance mutations (NRTI 55.6%, PI 68.4%, NNRTI 99.1%). The majority of NRTI-resistant strains (79.8%) carried resistance-associated mutations selected by the thymidine analogues zidovudine and stavudine. Preferably 2NRTI/1PIr combinations were prescribed as first line regimen in patients with resistant HIV as well as in patients with susceptible strains (susceptible 45.3%; 173/382 vs. resistant 65.5%; 40/61). The majority of patients in both groups were treated successfully within the first year after ART-initiation (susceptible: 89.9%; 62/69; resistant: 7/9; 77.8%).

CONCLUSION

Overall prevalence of TDR remained stable at a high level but trends of resistance against drug classes differed over time. The significant decrease of NRTI-resistance in patients newly infected with HIV might be related to the introduction of novel antiretroviral drugs and a wider use of genotypic resistance analysis prior to treatment initiation.

Short communication: Biological and genetic characterization of HIV type 1 subtype B and nonsubtype B transmitted viruses: usefulness for vaccine candidate assessment

Due to the extraordinary degree of genetic diversity of HIV-1 and the structural complexity of its envelope glycoproteins, designing an effective vaccine is difficult, requiring the development of viral reagents to assess vaccine-elicited neutralizing antibodies. The aim of this study was to improve on our previously developed panel of HIV-1 strains of different genetic forms, focusing on strains from acute and recently acquired infections as the most representative of the transmitted viruses. HIV-1 primary isolates were expanded in peripheral blood mononuclear cells. Viral stocks of 40 ml each were produced. Syncytium-inducing (SI) phenotype, coreceptor use, and TCID(50)/ml were determined. Near full-length HIV-1 genomes were amplified by RT-nested PCR in four overlapping segments. Phylogenetic analyses were performed with neighbor-joining trees and bootscanning. Forty-four HIV-1 strains were included in the panel. Twenty-four (54.1%) strains were from early infections (16 acute and 8 recent); of them, 21 (87%) were sexually transmitted. NSI/R5 phenotype was detected in 37 (84.1%) viruses and SI/R5,X4 in another 7 (15.9%). TCID(50)/ml ranged between 10(4) and 10(6.6). Twelve different genetic forms constituted this panel: subtypes A1, B, C, F1, and G; circulating recombinant forms CRF02_AG, CRF14_BG, and CRF24_BG; and unique recombinant forms CRF02_AG/A3, BF1, CRF12_BF/B, and DF1G. In conclusion, in this study, we report the development of a comprehensive and well-characterized panel of HIV-1 isolates for assessing neutralization in HIV vaccine research. This panel is available for distribution through the Programme EVA Centre for AIDS Reagents, National Institute for Biological Standard and Control (NIBSC).

Identification of a new HIV type 1 circulating BF intersubtype recombinant form (CRF47_BF) in Spain

We report the identification of a new HIV-1 circulating recombinant form (CRF47_BF) derived from subtypes B and F. It was initially identified in protease-reverse transcriptase sequences from nine individuals from three separate regions of Spain who acquired HIV-1 infection via sexual contact. All nine sequences formed a strongly supported phylogenetic cluster, branching apart from all known CRFs, and in bootscan analyses were BF mosaics with two coincident breakpoints. Two epidemiologically unlinked viruses were sequenced in near full-length genomes, which exhibited identical mosaic structures, with 16 intersubtype breakpoints in a genome predominantly of subtype B. Subtype F segments of the new CRF failed to cluster with any of the near full-length genome subtype F sequences available in public databases. Recent dates of HIV-1 diagnoses and short genetic distances suggest a recent origin of this CRF. This is the tenth reported CRF_BF, the first apparently having originated outside of South America.

Short communication: Molecular epidemiology of HIV type 1 in the Republic of Dagestan, Russian Federation: virtually uniform circulation of subtype A, former Soviet Union variant, with predominance of the V77I(PR) subvariant

We examine the distribution of viral genetic forms and the presence of antiretroviral drug resistance mutations in HIV-1 infections in the Republic of Dagestan, in the North Caucasus area of Russia, where a recent large increase in HIV-1 infections has been documented. Samples were collected from 41 HIV-1-infected individuals from Dagestan, most of them from the cities of Derbent (n = 21) and Mahachkala (n = 18). Thirty six were injecting drug users and five were infected by heterosexual contact. None was on antiretroviral drug treatment. HIV-1 protease and a segment of reverse transcriptase were amplified by RT-PCR from plasma RNA and sequenced, and phylogenetic trees were constructed via maximum likelihood. Forty (97.6%) of 41 samples were of subtype A, former Soviet Union variant (A(FSU)), of which 27 (67.5%) clustered with the subvariant containing the V77I substitution in protease (V77I(PR)). Within this cluster, 13 viruses formed a local subcluster, 10 of which were from Derbent. Four viruses clustered with the A(SP2) subcluster, recently identified in St. Petersburg, two with a virus from Georgia and one with a virus from Azerbaijan. No mutations associated with antiretroviral drug resistance were detected. The results, therefore, show the relationship of the HIV-1 epidemic in Dagestan with that of other areas of Russia and of neighboring countries, and reveal the spread of the A(FSU) V77I(PR) variant in the North Caucasus area.