Tracing the Impact of Public Health Interventions on HIV-1 Transmission in Portugal Using Molecular Epidemiology.

BACKGROUND: Estimation of temporal changes in human immunodeficiency virus (HIV) transmission patterns can help to elucidate the impact of preventive strategies and public health policies. METHODS: Portuguese HIV-1 subtype B and G pol genetic sequences were appended to global reference data sets to identify country-specific transmission clades. Bayesian birth-death models were used to estimate subtype-specific effective reproductive numbers (Re). Discrete trait analysis (DTA) was used to quantify mixing among transmission groups. RESULTS: We identified 5 subtype B Portuguese clades (26-79 sequences) and a large monophyletic subtype G Portuguese clade (236 sequences). We estimated that major shifts in HIV-1 transmission occurred around 1999 (95% Bayesian credible interval [BCI], 1998-2000) and 2000 (95% BCI, 1998-2001) for subtypes B and G, respectively. For subtype B, Re dropped from 1.91 (95% BCI, 1.73-2.09) to 0.62 (95% BCI,.52-.72). For subtype G, Re decreased from 1.49 (95% BCI, 1.39-1.59) to 0.72 (95% BCI, .63-.8). The DTA suggests that people who inject drugs (PWID) and heterosexuals were the source of most (>80%) virus lineage transitions for subtypes G and B, respectively. CONCLUSIONS: The estimated declines in Re coincide with the introduction of highly active antiretroviral therapy and the scale-up of harm reduction for PWID. Inferred transmission events across transmission groups emphasize the importance of prevention efforts for bridging populations.

The global spread of HIV-1 subtype B epidemic.

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

Sub-Epidemics Explain Localized High Prevalence of Reduced Susceptibility to Rilpivirine in Treatment-Naive HIV-1-Infected Patients: Subtype and Geographic Compartmentalization of Baseline Resistance Mutations.

OBJECTIVE: The latest nonnucleoside reverse transcriptase inhibitor (NNRTI) rilpivirine (RPV) is indicated for human immunodeficiency virus type-1 (HIV-1) patients initiating antiretroviral treatment, but the extent of genotypic RPV resistance in treatment-naive patients outside clinical trials is poorly defined. STUDY DESIGN: This retrospective observational study of clinical data from Belgium and Portugal evaluates genotypic information from HIV-1 drug-naive patients obtained for the purpose of drug resistance testing. Rilpivirine resistance-associated mutations (RPV-RAMs) were defined based on clinical trials, phenotypic studies, and expert-based resistance algorithms. Viral susceptibility to RPV alone and to the single-tablet regimen was estimated using expert-based resistance algorithms. RESULTS: In 4,631 HIV-1 treatment-naive patients infected with diverse HIV-1 subtypes, major RPV-RAMs were detected in 4.6%, while complete viral susceptibility to RPV was estimated in 95% of patients. Subtype C- and F1-infected patients displayed the highest levels of reduced viral susceptibility at baseline, respectively 13.2% and 9.3%, mainly due to subtype- and geographic-dependent occurrence of RPV-RAMs E138A and A98G as natural polymorphisms. Strikingly, a founder effect in Portugal resulted in a 138A prevalence of 13.2% in local subtype C-infected treatment-naive patients. The presence of transmitted drug resistance did not impact our estimates. CONCLUSION: RPV is the first HIV-1 inhibitor for which, in the absence of transmitted drug resistance, intermediate or high-level genotypic resistance can be detected in treatment-naive patients. The extent of RPV susceptibility in treatment-naive patients differs depending on the HIV-1 subtype and dynamics of local compartmentalized epidemics. The highest prevalence of reduced susceptibility was found to be 15.7% in Portuguese subtype C-infected treatment-naive patients. In this context, even in the absence of transmitted HIV-1 drug resistance (TDR), drug resistance testing at baseline should be considered extremely important before starting treatment with this NNRTI.

Assessing transmissibility of HIV-1 drug resistance mutations from treated and from drug-naive individuals.

OBJECTIVES: Surveillance drug resistance mutations (SDRMs) in drug-naive patients are typically used to survey HIV-1-transmitted drug resistance (TDR). We test here how SDRMs in patients failing treatment, the original source of TDR, contribute to assessing TDR, transmissibility and transmission source of SDRMs. DESIGN: This is a retrospective observational study analyzing a Portuguese cohort of HIV-1-infected patients. METHODS: The prevalence of SDRMs to protease inhibitors, nucleoside reverse transcriptase inhibitors (NRTIs) and nonnucleoside reverse transcriptase inhibitors (NNRTIs) in drug-naive and treatment-failing patients was measured for 3554 HIV-1 subtype B patients. Transmission ratio (prevalence in drug-naive/prevalence in treatment-failing patients), average viral load and robust linear regression with outlier detection (prevalence in drug-naive versus in treatment-failing patients) were analyzed and used to interpret transmissibility. RESULTS: Prevalence of SDRMs in drug-naive and treatment-failing patients were linearly correlated, but some SDRMs were classified as outliers - above (PRO: D30N, N88D/S, L90 M, RT: G190A/S/E) or below (RT: M184I/V) expectations. The normalized regression slope was 0.073 for protease inhibitors, 0.084 for NRTIs and 0.116 for NNRTIs. Differences between SDRMs transmission ratios were not associated with differences in viral loads. CONCLUSION: The significant linear correlation between prevalence of SDRMs in drug-naive and in treatment-failing patients indicates that the prevalence in treatment-failing patients can be useful to predict levels of TDR. The slope is a cohort-dependent estimate of rate of TDR per drug class and outlier detection reveals comparative persistence of SDRMs. Outlier SDRMs with higher transmissibility are more persistent and more likely to have been acquired from drug-naive patients. Those with lower transmissibility have faster reversion dynamics after transmission and are associated with acquisition from treatment-failing patients.

Discordant predictions of residual activity could impact dolutegravir prescription upon raltegravir failure.

BACKGROUND: Dolutegravir is approved for the treatment of HIV-1 patients exposed to other integrase inhibitors, but the decision to use dolutegravir in this setting should be informed by drug resistance testing. OBJECTIVES: This study determined the extent of disagreement in predicted residual dolutegravir activity after raltegravir use, and identified individual mutational patterns for which uncertainty exists among HIV-1 expert systems. STUDY DESIGN: Mutation patterns were classified in raltegravir signature pathways including positions 143, 148 and 155, and interpreted into clinically informative resistance levels using genotypic drug resistance interpretation systems ANRS v24, HIVdb v7.0 and Rega v9.1.0, and instructions of dolutegravir use as approved by the Food and Drug Administration and the European Medicines Agency. RESULTS: In 216HIV-1 patients failing raltegravir-therapy, 87% patients displayed mutations associated with resistance towards integrase inhibitors. A total of 141 unique mutational patterns were observed, with N155H (25.4%), Q148H (16.2%) and Y143R (8.3%) the most prevalent signature mutations. The Q148 pathway occurred almost exclusively in HIV-1 subtype B viruses. Concordances in predicted dolutegravir susceptibility scores among 5 systems were obtained in 57.8% of patients, and concordant intermediate resistant and concordant resistant scores were only observed in 6.5% and 0.9% of patients, respectively. However, systems individually scored higher levels of dolutegravir intermediate resistance and resistance, ranging from 4.2% to 10.2% and from 14.8% to 22.7% of patients, respectively. A consensus on interpreting the extent of residual activity was lacking in 34.7% of patients and was highly resistance pathway-specific. CONCLUSIONS: Dolutegravir may potentially be effective in the majority of HIV-1 patients failing raltegravir, but concern over the uncertainty in predicted residual activity could withhold clinicians from prescribing dolutegravir during its clinical assessment.

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.

Pharmacy refill adherence outperforms self-reported methods in predicting HIV therapy outcome in resource-limited settings.

BACKGROUND: Optimal adherence to antiretroviral therapy is critical to prevent HIV drug resistance (HIVDR) epidemic. The objective of the study was to investigate the best performing adherence assessment method for predicting virological failure in resource-limited settings (RLS). METHOD: This study was a single-centre prospective cohort, enrolling 220 HIV-infected adult patients attending an HIV/AIDS Care and Treatment Centre in Dar es Salaam, Tanzania, in 2010. Pharmacy refill, self-report (via visual analog scale [VAS] and the Swiss HIV Cohort study-adherence questionnaire), pill count, and appointment keeping adherence measurements were taken. Univariate logistic regression (LR) was done to explore a cut-off that gives a better trade-off between sensitivity and specificity, and a higher area under the curve (AUC) based on receiver operating characteristic curve in predicting virological failure. Additionally, the adherence models were evaluated by fitting multivariate LR with stepwise functions, decision trees, and random forests models, assessing 10-fold multiple cross validation (MCV). Patient factors associated with virological failure were determined using LR. RESULTS: Viral load measurements at baseline and one year after recruitment were available for 162 patients, of whom 55 (34%) had detectable viral load and 17 (10.5%) had immunological failure at one year after recruitment. The optimal cut-off points significantly predictive of virological failure were 95%, 80%, 95% and 90% for VAS, appointment keeping, pharmacy refill, and pill count adherence respectively. The AUC for these methods ranged from 0.52 to 0.61, with pharmacy refill giving the best performance at AUC 0.61. Multivariate logistic regression with boost stepwise MCV had higher AUC (0.64) compared to all univariate adherence models, except pharmacy refill adherence univariate model, which was comparable to the multivariate model (AUC = 0.64). Decision trees and random forests models were inferior to boost stepwise model. Pharmacy refill adherence (<95%) emerged as the best method for predicting virological failure. Other significant predictors in multivariate LR were having a baseline CD4 T lymphocytes count < 200 cells/µl, being unable to recall the diagnosis date, and a higher weight. CONCLUSION: Pharmacy refill has the potential to predict virological failure and to identify patients to be considered for viral load monitoring and HIVDR testing in RLS.

A near-full length genotypic assay for HCV1b.

A near-full genome genotypic assay for HCV1b was developed, which may prove useful to investigate antiviral drug resistance, given new combination therapies for HCV1 infection. The assay consists of three partially overlapping PCRs followed by Sanger population or Illumina next-generation sequencing. Seventy-seven therapy-naïve samples, spanning the entire diversity range of currently known HCV1b, were used for optimization of PCRs, of which ten were sequenced using Sanger and of these ten, four using Illumina. The median detection limits for the three regions, 5'UTR-NS2, E2-NS5A and NS4B-NS5B, were 570, 5670 and 56,670 IU/ml respectively. The number of Illumina reads mapped varied according to the software used, Segminator II being the best performing (81%). Consensus Illumina and Sanger sequencing results accord largely (0.013% major discordances). Differences were due almost exclusively to a larger number of ambiguities (presumably minority variants) scored by Illumina (1.50% minor discordances). The assay is easy to perform in an equipped laboratory; nevertheless, it was difficult to reach high sensitivity and reproducibility, due to the high genetic viral variability. This assay proved to be suitable for detecting drug resistance mutations and can also be used for epidemiological research, even though only a limited set of samples was used for validation.

Increase in transmitted resistance to non-nucleoside reverse transcriptase inhibitors among newly diagnosed HIV-1 infections in Europe.

BACKGROUND: One out of ten newly diagnosed patients in Europe was infected with a virus carrying a drug resistant mutation. We analysed the patterns over time for transmitted drug resistance mutations (TDRM) using data from the European Spread program. METHODS: Clinical, epidemiological and virological data from 4317 patients newly diagnosed with HIV-1 infection between 2002 and 2007 were analysed. Patients were enrolled using a pre-defined sampling strategy. RESULTS: The overall prevalence of TDRM in this period was 8.9% (95% CI: 8.1-9.8). Interestingly, significant changes over time in TDRM caused by the different drug classes were found. Whereas nucleoside resistance mutations remained constant at 5%, a significant decline in protease inhibitors resistance mutations was observed, from 3.9% in 2002 to 1.6% in 2007 (p = 0.001). In contrast, resistance to non-nucleoside reverse transcriptase inhibitors (NNRTIs) doubled from 2.0% in 2002 to 4.1% in 2007 (p = 0.004) with 58% of viral strains carrying a K103N mutation. Phylogenetic analysis showed that these temporal changes could not be explained by large clusters of TDRM. CONCLUSION: During the years 2002 to 2007 transmitted resistance to NNRTI has doubled to 4% in Europe. The frequent use of NNRTI in first-line regimens and the clinical impact of NNRTI mutations warrants continued monitoring.

RegaDB: community-driven data management and analysis for infectious diseases.

SUMMARY: RegaDB is a free and open source data management and analysis environment for infectious diseases. RegaDB allows clinicians to store, manage and analyse patient data, including viral genetic sequences. Moreover, RegaDB provides researchers with a mechanism to collect data in a uniform format and offers them a canvas to make newly developed bioinformatics tools available to clinicians and virologists through a user friendly interface. AVAILABILITY AND IMPLEMENTATION: Source code, binaries and documentation are available on http://rega.kuleuven.be/cev/regadb. RegaDB is written in the Java programming language, using a web-service-oriented architecture.

Automated subtyping of HIV-1 genetic sequences for clinical and surveillance purposes: performance evaluation of the new REGA version 3 and seven other tools.

BACKGROUND: To investigate differences in pathogenesis, diagnosis and resistance pathways between HIV-1 subtypes, an accurate subtyping tool for large datasets is needed. We aimed to evaluate the performance of automated subtyping tools to classify the different subtypes and circulating recombinant forms using pol, the most sequenced region in clinical practice. We also present the upgraded version 3 of the Rega HIV subtyping tool (REGAv3). METHODOLOGY: HIV-1 pol sequences (PR+RT) for 4674 patients retrieved from the Portuguese HIV Drug Resistance Database, and 1872 pol sequences trimmed from full-length genomes retrieved from the Los Alamos database were classified with statistical-based tools such as COMET, jpHMM and STAR; similarity-based tools such as NCBI and Stanford; and phylogenetic-based tools such as REGA version 2 (REGAv2), REGAv3, and SCUEAL. The performance of these tools, for pol, and for PR and RT separately, was compared in terms of reproducibility, sensitivity and specificity with respect to the gold standard which was manual phylogenetic analysis of the pol region. RESULTS: The sensitivity and specificity for subtypes B and C was more than 96% for seven tools, but was variable for other subtypes such as A, D, F and G. With regard to the most common circulating recombinant forms (CRFs), the sensitivity and specificity for CRF01_AE was ~99% with statistical-based tools, with phylogenetic-based tools and with Stanford, one of the similarity based tools. CRF02_AG was correctly identified for more than 96% by COMET, REGAv3, Stanford and STAR. All the tools reached a specificity of more than 97% for most of the subtypes and the two main CRFs (CRF01_AE and CRF02_AG). Other CRFs were identified only by COMET, REGAv2, REGAv3, and SCUEAL and with variable sensitivity. When analyzing sequences for PR and RT separately, the performance for PR was generally lower and variable between the tools. Similarity and statistical-based tools were 100% reproducible, but this was lower for phylogenetic-based tools such as REGA (~99%) and SCUEAL (~96%). CONCLUSIONS: REGAv3 had an improved performance for subtype B and CRF02_AG compared to REGAv2 and is now able to also identify all epidemiologically relevant CRFs. In general the best performing tools, in alphabetical order, were COMET, jpHMM, REGAv3, and SCUEAL when analyzing pure subtypes in the pol region, and COMET and REGAv3 when analyzing most of the CRFs. Based on this study, we recommend to confirm subtyping with 2 well performing tools, and be cautious with the interpretation of short sequences.

Limited cross-border infections in patients newly diagnosed with HIV in Europe.

BACKGROUND: International travel plays a role in the spread of HIV-1 across Europe. It is, however, not known whether international travel is more important for spread of the epidemic as compared to endogenous infections within single countries. In this study, phylogenetic associations among HIV of newly diagnosed patients were determined across Europe. RESULTS: Data came from the SPREAD programme which collects samples of newly diagnosed patients that are representative for national HIV epidemics. 4260 pol sequences from 25 European countries and Israel collected in 2002-2007 were included.We identified 457 clusters including 1330 persons (31.2% of all patients). The cluster size ranged between 2 and 28. A number of 987 patients (74.2%) were part of a cluster that consisted only of patients originating from the same country. In addition, 135 patients (10.2%) were in a cluster including only individuals from neighboring countries. Finally, 208 patients (15.6%) clustered with individuals from countries without a common border. Clustering with patients from the same country was less prevalent in patients being infected with B subtype (P-value <0.0001), in men who have sex with men (P-value <0.0001), and in recently infected patients (P-value =0.045). CONCLUSIONS: Our findings indicate that the transmission of HIV-1 in Europe is predominantly occurring between patients from the same country. This could have implications for HIV-1 transmission prevention programmes. Because infections through travelling between countries is not frequently observed it is important to have good surveillance of the national HIV-1 epidemics.

HIV-1 fitness landscape models for indinavir treatment pressure using observed evolution in longitudinal sequence data are predictive for treatment failure.

We previously modeled the in vivo evolution of human immunodeficiency virus-1 (HIV-1) under drug selective pressure from cross-sectional viral sequences. These fitness landscapes (FLs) were made by using first a Bayesian network (BN) to map epistatic substitutions, followed by scaling the fitness landscape based on an HIV evolution simulator trying to evolve the sequences from treatment naïve patients into sequences from patients failing treatment. In this study, we compared four FLs trained with different sequence populations. Epistatic interactions were learned from three different cross-sectional BNs, trained with sequence from patients experienced with indinavir (BNT), all protease inhibitors (PIs) (BNP) or all PI except indinavir (BND). Scaling the fitness landscape was done using cross-sectional data from drug naïve and indinavir experienced patients (Fcross using BNT) and using longitudinal sequences from patients failing indinavir (FlongT using BNT, FlongP using BNP, FlongD using BND). Evaluation to predict the failing sequence and therapy outcome was performed on independent sequences of patients on indinavir. Parameters included estimated fitness (LogF), the number of generations (GF) or mutations (MF) to reach the fitness threshold (average fitness when a major resistance mutation appeared), the number of generations (GR) or mutations (MR) to reach a major resistance mutation and compared to genotypic susceptibility score (GSS) from Rega and HIVdb algorithms. In pairwise FL comparisons we found significant correlation between fitness values for individual sequences, and this correlation improved after correcting for the subtype. Furthermore, FLs could predict the failing sequence under indinavir-containing combinations. At 12 and 48 weeks, all parameters from all FLs and indinavir GSS (both for Rega and HIVdb) were predictive of therapy outcome, except MR for FlongT and FlongP. The fitness landscapes have similar predictive power for treatment response under indinavir-containing regimen as standard rules-based algorithms, and additionally allow predicting genetic evolution under indinavir selective pressure.

HIV-1 subtype distribution and its demographic determinants in newly diagnosed patients in Europe suggest highly compartmentalized epidemics.

BACKGROUND: Understanding HIV-1 subtype distribution and epidemiology can assist preventive measures and clinical decisions. Sequence variation may affect antiviral drug resistance development, disease progression, evolutionary rates and transmission routes. RESULTS: We investigated the subtype distribution of HIV-1 in Europe and Israel in a representative sample of patients diagnosed between 2002 and 2005 and related it to the demographic data available. 2793 PRO-RT sequences were subtyped either with the REGA Subtyping tool or by a manual procedure that included phylogenetic tree and recombination analysis. The most prevalent subtypes/CRFs in our dataset were subtype B (66.1%), followed by sub-subtype A1 (6.9%), subtype C (6.8%) and CRF02_AG (4.7%). Substantial differences in the proportion of new diagnoses with distinct subtypes were found between European countries: the lowest proportion of subtype B was found in Israel (27.9%) and Portugal (39.2%), while the highest was observed in Poland (96.2%) and Slovenia (93.6%). Other subtypes were significantly more diagnosed in immigrant populations. Subtype B was significantly more diagnosed in men than in women and in MSM > IDUs > heterosexuals. Furthermore, the subtype distribution according to continent of origin of the patients suggests they acquired their infection there or in Europe from compatriots. CONCLUSIONS: The association of subtype with demographic parameters suggests highly compartmentalized epidemics, determined by social and behavioural characteristics of the patients.

Treatment-associated polymorphisms in protease are significantly associated with higher viral load and lower CD4 count in newly diagnosed drug-naive HIV-1 infected patients.

BACKGROUND: The effect of drug resistance transmission on disease progression in the newly infected patient is not well understood. Major drug resistance mutations severely impair viral fitness in a drug free environment, and therefore are expected to revert quickly. Compensatory mutations, often already polymorphic in wild-type viruses, do not tend to revert after transmission. While compensatory mutations increase fitness during treatment, their presence may also modulate viral fitness and virulence in absence of therapy and major resistance mutations. We previously designed a modeling technique that quantifies genotypic footprints of in vivo treatment selective pressure, including both drug resistance mutations and polymorphic compensatory mutations, through the quantitative description of a fitness landscape from virus genetic sequences. RESULTS: Genotypic correlates of viral load and CD4 cell count were evaluated in subtype B sequences from recently diagnosed treatment-naive patients enrolled in the SPREAD programme. The association of surveillance drug resistance mutations, reported compensatory mutations and fitness estimated from drug selective pressure fitness landscapes with baseline viral load and CD4 cell count was evaluated using regression techniques. Protease genotypic variability estimated to increase fitness during treatment was associated with higher viral load and lower CD4 cell counts also in treatment-naive patients, which could primarily be attributed to well-known compensatory mutations at highly polymorphic positions. By contrast, treatment-related mutations in reverse transcriptase could not explain viral load or CD4 cell count variability. CONCLUSIONS: These results suggest that polymorphic compensatory mutations in protease, reported to be selected during treatment, may improve the replicative capacity of HIV-1 even in absence of drug selective pressure or major resistance mutations. The presence of this polymorphic variation may either reflect a history of drug selective pressure, i.e. transmission from a treated patient, or merely be a result of diversity in wild-type virus. Our findings suggest that transmitted drug resistance has the potential to contribute to faster disease progression in the newly infected host and to shape the HIV-1 epidemic at a population level.

Effect of natural polymorphisms in the HIV-1 CRF02_AG protease on protease inhibitor hypersusceptibility.

Hypersusceptibility (HS) to inhibition by different antiretroviral drugs (ARVs) among diverse HIV-infected individuals may be a misnomer because clinical response to treatment is evaluated in relation to subtype B infections while drug susceptibility of the infecting virus, regardless of subtype, is compared to a subtype B HIV-1 laboratory strain (NL4-3 or IIIB). Mounting evidence suggests that HS to different ARVs may result in better treatment outcome just as drug resistance leads to treatment failure. We have identified key amino acid polymorphisms in the protease coding region of a non-B HIV-1 subtype linked to protease inhibitor HS, namely, 17E and 64M in CRF02_AG. These HS-linked polymorphisms were introduced in the BD6-15 CRF02_AG molecular clone and tested for inhibition using a panel of protease inhibitors. In general, suspected HS-linked polymorphisms did increase susceptibility to specific protease inhibitors such as amprenavir and atazanavir, but the combination of the 17E/64M polymorphisms showed greater HS. These two mutations were found at low frequencies but linked in a sequence database of over 700 protease sequences of CRF02_AG. In direct head-to-head virus competitions, CRF02_AG harboring the 17E/64M polymorphisms also had higher replicative fitness than did the 17E or the 64M polymorphism in the CFR02_AG clone. These findings suggest that subtype-specific, linked polymorphisms can result in hypersusceptibility to ARVs. Considering the potential benefit of HS to treatment outcome, screening for potential HS-linked polymorphisms as well as preexisting drug resistance mutations in treatment-naïve patients may guide the choice of ARVs for the best treatment outcome.

Phylodynamics of the HIV-1 CRF02_AG clade in Cameroon.

Evolutionary analyses have revealed an origin of pandemic HIV-1 group M in the Congo River basin in the first part of the XX century, but the patterns of historical viral spread in or around its epicentre remain largely unexplored. Here, we combine epidemiologic and molecular sequence data to investigate the spatiotemporal patterns of the CRF02_AG clade. By explicitly integrating prevalence counts and genetic population size estimates we date the epidemic emergence of CRF02_AG at 1973.1 (1972.1, 1975.3, 95% CI). To infer the phylogeographic signature of this clade at a regional scale, we analyze pol and env time-stamped sequence data from 10 countries using a Bayesian phylogeographic approach based on an asymmetric discretized diffusion model. Our data confirms a spatial origin of CRF02_AG in the Democratic Republic of Congo (DRC) and suggests that viral dissemination to Cameroon occurred at an early stage of the evolutionary history of CRF02_AG. We find considerable support for epidemiological linkage between neighbour countries. Compilation of ethnographic data suggested that well-supported viral migration did not reflect sustained human migratory flows. Finally, using sequence data from 15 locations in Cameroon, we use relaxed random walk models to explore the spatiotemporal dynamics of CRF02_AG at a finer geographical detail. Phylogeographic dispersal in continuous space reveals that at least two distinct CRF02_AG lineages are circulating in overlapping regions that are evolving at different evolutionary and diffusion rates. In conclusion, by combining molecular and epidemiological data, our results provide a time scale for CRF02_AG, early 70s, place its spatial root in the DRC within the putative root of group-M diversity and propose a scenario of chance-exportation events for the spatiotemporal patterns of a successful HIV-1 lineage both at a regional and country-scale.

Phylogeographical footprint of colonial history in the global dispersal of human immunodeficiency virus type 2 group A.

Human immunodeficiency virus type 2 (HIV-2) emerged in West Africa and has spread further to countries that share socio-historical ties with this region. However, viral origins and dispersal patterns at a global scale remain poorly understood. Here, we adopt a Bayesian phylogeographic approach to investigate the spatial dynamics of HIV-2 group A (HIV-2A) using a collection of 320 partial pol and 248 partial env sequences sampled throughout 19 countries worldwide. We extend phylogenetic diffusion models that simultaneously draw information from multiple loci to estimate location states throughout distinct phylogenies and explicitly attempt to incorporate human migratory fluxes. Our study highlights that Guinea-Bissau, together with Côte d'Ivoire and Senegal, have acted as the main viral sources in the early stages of the epidemic. We show that convenience sampling can obfuscate the estimation of the spatial root of HIV-2A. We explicitly attempt to circumvent this by incorporating rate priors that reflect the ratio of human flow from and to West Africa. We recover four main routes of HIV-2A dispersal that are laid out along colonial ties: Guinea-Bissau and Cape Verde to Portugal, Côte d'Ivoire and Senegal to France. Within Europe, we find strong support for epidemiological linkage from Portugal to Luxembourg and to the UK. We demonstrate that probabilistic models can uncover global patterns of HIV-2A dispersal providing sampling bias is taken into account and we provide a scenario for the international spread of this virus.

European recommendations for the clinical use of HIV drug resistance testing: 2011 update.

The European HIV Drug Resistance Guidelines Panel, established to make recommendations to clinicians and virologists, felt that sufficient new information has become available to warrant an update of its recommendations, explained in both pocket guidelines and this full paper. The Panel makes the following recommendations concerning the indications for resistance testing: for HIV-1 (i) test earliest sample for protease and reverse transcriptase drug resistance in drug-naive patients with acute or chronic infection; (ii) test protease and reverse transcriptase drug resistance at virologic failure, and other drug targets (integrase and envelope) if such drugs were part of the failing regimen; (iii) consider testing for CCR5 tropism at virologic failure or when a change of therapy has to be made in absence of detectable viral load, and in the latter case test DNA or last detectable plasma RNA; (iv) consider testing earliest detectable plasma RNA when a successful nonnucleoside reverse transcriptase inhibitor-containing therapy was inappropriately interrupted; (v) genotype source patient when postexposure prophylaxis is considered; for HIV-2, (vi) consider resistance testing where treatment change is needed after treatment failure. The Panel recommends genotyping in most situations, using updated and clinically evaluated interpretation systems. It is mandatory that laboratories performing HIV resistance tests take part regularly in external quality assurance programs, and that they consider storing samples in situations where resistance testing cannot be performed as recommended. Similarly, it is necessary that HIV clinicians and virologists take part in continuous education and discuss problematic clinical cases. Indeed, resistance test results should be used in the context of all other clinically relevant information for predicting therapy response.

Appearance of a single amino acid insertion at position 33 in HIV type 1 protease under a lopinavir-containing regimen, associated with reduced protease inhibitor susceptibility.

HIV drug resistance is a multifactorial phenomenon and constitutes a major concern as it results in therapy failure. The aim of this study was to assess the impact of an amino acid insertion identified at position 33 of the protease gene, derived from samples from three patients under lopinavir therapy, on viral fitness and protease inhibitor (PI) resistance. Successive samples were available from one of the patients for genotypic and phenotypic testing in order to investigate the role of this insertion. The patient had been pretreated with various antiretroviral drugs and showed poor virological response from the point of the acquisition of the mutation onward. The insertion was acquired in the context of a number of other PI mutations and was stable following acquisition. Phenotypic testing revealed reduced susceptibility to various PIs and a reduction of the replicative capacity (RC) of the virus. In the presence of the insertion alone, a decrease of the RC was observed, which seemed to be compensated by the presence of other mutations. The L33ins might have a potential role in PI resistance pathways but further investigation in a larger number of clinical samples is required in order to elucidate this resistance mechanism.