Search

Objectives: The use of tenofovir is highly associated with the emergence of mutation K65R, which confers broad resistance to nucleoside/nucleotide analogue reverse transcriptase inhibitors (NRTIs), especially when tenofovir is combined with other NRTIs also selecting for K65R. Although recent HIV-1 treatment guidelines discouraging these combinations resulted in reduced K65R selection with tenofovir, updated information on the impact of currently recommended regimens on the population selection rate of K65R is presently lacking.Methods: In this study, we evaluated changes over time in the selection rate of resistance mutation K65R in a large population of 2736 HIV-1-infected patients failing combination antiretroviral treatment between 2002 and 2010.Results: The K65R resistance mutation was detected in 144 patients, a prevalence of 5.3%. A large majority of observed K65R cases were explained by the use of tenofovir, reflecting its wide use in clinical practice. However, changing patterns over time in NRTIs accompanying tenofovir resulted in a persistent decreasing probability of K65R selection by tenofovir-based therapy. The currently recommended NRTI combination tenofovir/emtricitabine was associated with a low probability of K65R emergence. For any given dual NRTI combination including tenofovir, higher selection rates of K65R were consistently observed with a non-nucleoside reverse transcriptase inhibitor than with a protease inhibitor as the third agent.Discussion: Our finding of a stable time trend of K65R despite elevated use of tenofovir illustrates increased potency of current HIV-1 therapy including tenofovir. [ABSTRACT FROM AUTHOR]

The diagnosis of chylothorax during pregnancy is rare and is generally triggered during labor. The study of pleural fluid and imaging studies can facilitate diagnosis, even though determining its etiology can be complicated. This case involves a young pregnant woman presenting with disabling dyspnea and chest pain due to pleural effusion. The conventional study of the pleural fluid was not conclusive; however, the measurement of chylomicrons guided us to determine the diagnosis of chylothorax. Based on imaging studies and specific markers, some of the main etiologies were not determined, leading to classification as idiopathic. After drainage and dietary adjustment, a significant improvement in symptoms was achieved with the help of an interdisciplinary team, which is crucial for the prevention of complications and adequate resolution of the pregnancy., Competing Interests: The authors have declared that no competing interests exist., (Copyright © 2023, Martínez Camacho et al.)

Background: Portugal has one of the most severe HIV-1 epidemics in Western Europe. Two subtypes circulate in parallel since the beginning of the epidemic. Comparing their transmission patterns and its association with transmitted drug resistance (TDR) is important to pinpoint transmission hotspots and to develop evidence-based treatment guidelines., Methods: Demographic, clinical and genomic data were collected from 3599 HIV-1 naive patients between 2001 and 2014. Sequences obtained from drug resistance testing were used for subtyping, TDR determination and transmission clusters (TC) analyses., Results: In Portugal, transmission of subtype B was significantly associated with young males, while transmission of subtype G was associated with older heterosexuals. In Portuguese originated people, there was a decreasing trend both for prevalence of subtype G and for number of TCs in this subtype. The active TCs that were identified (i.e. clusters originated after 2008) were associated with subtype B-infected males residing in Lisbon. TDR was significantly different when comparing subtypes B (10.8% [9.5-12.2]) and G (7.6% [6.4-9.0]) (p = 0.001)., Discussion: TC analyses shows that, in Portugal, the subtype B epidemic is active and fueled by young male patients residing in Lisbon, while transmission of subtype G is decreasing. Despite similar treatment rates for both subtypes in Portugal, TDR is significantly different between subtypes., Competing Interests: ABA received funding from two commercial sources: LÓreal Portugal, through the research award L'Oréal Portugal Medals of Honor for Women in Science (2012) and Gilead Génese HIVLatePresenters. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

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 Author(s) 2019. Published by Oxford University Press for the Infectious Diseases Society of America.)

Background: CCR5-coreceptor antagonists can be used for treating HIV-2 infected individuals. Before initiating treatment with coreceptor antagonists, viral coreceptor usage should be determined to ensure that the virus can use only the CCR5 coreceptor (R5) and cannot evade the drug by using the CXCR4 coreceptor (X4-capable). However, until now, no online tool for the genotypic identification of HIV-2 coreceptor usage had been available. Furthermore, there is a lack of knowledge on the determinants of HIV-2 coreceptor usage. Therefore, we developed a data-driven web service for the prediction of HIV-2 coreceptor usage from the V3 loop of the HIV-2 glycoprotein and used the tool to identify novel discriminatory features of X4-capable variants., Results: Using 10 runs of tenfold cross validation, we selected a linear support vector machine (SVM) as the model for geno2pheno[coreceptor-hiv2], because it outperformed the other SVMs with an area under the ROC curve (AUC) of 0.95. We found that SVMs were highly accurate in identifying HIV-2 coreceptor usage, attaining sensitivities of 73.5% and specificities of 96% during tenfold nested cross validation. The predictive performance of SVMs was not significantly different (p value 0.37) from an existing rules-based approach. Moreover, geno2pheno[coreceptor-hiv2] achieved a predictive accuracy of 100% and outperformed the existing approach on an independent data set containing nine new isolates with corresponding phenotypic measurements of coreceptor usage. geno2pheno[coreceptor-hiv2] could not only reproduce the established markers of CXCR4-usage, but also revealed novel markers: the substitutions 27K, 15G, and 8S were significantly predictive of CXCR4 usage. Furthermore, SVMs trained on the amino-acid sequences of the V1 and V2 loops were also quite accurate in predicting coreceptor usage (AUCs of 0.84 and 0.65, respectively)., Conclusions: In this study, we developed geno2pheno[coreceptor-hiv2], the first online tool for the prediction of HIV-2 coreceptor usage from the V3 loop. Using our method, we identified novel amino-acid markers of X4-capable variants in the V3 loop and found that HIV-2 coreceptor usage is also influenced by the V1/V2 region. The tool can aid clinicians in deciding whether coreceptor antagonists such as maraviroc are a treatment option and enables epidemiological studies investigating HIV-2 coreceptor usage. geno2pheno[coreceptor-hiv2] is freely available at http://coreceptor-hiv2.geno2pheno.org .

Background: Epidemic HIV-2 (groups A and B) emerged in humans circa 1930-40. Its closest ancestors are SIVsmm infecting sooty mangabeys from southwestern Côte d'Ivoire. The earliest large-scale serological surveys of HIV-2 in West Africa (1985-91) show a patchy spread. Côte d'Ivoire and Guinea-Bissau had the highest prevalence rates by then, and phylogeographical analysis suggests they were the earliest epicenters. Wars and parenteral transmission have been hypothesized to have promoted HIV-2 spread. Male circumcision (MC) is known to correlate negatively with HIV-1 prevalence in Africa, but studies examining this issue for HIV-2 are lacking., Methods: We reviewed published HIV-2 serosurveys for 30 cities of all West African countries and obtained credible estimates of real prevalence through Bayesian estimation. We estimated past MC rates of 218 West African ethnic groups, based on ethnographic literature and fieldwork. We collected demographic tables specifying the ethnic partition in cities. Uncertainty was incorporated by defining plausible ranges of parameters (e.g. timing of introduction, proportion circumcised). We generated 1,000 sets of past MC rates per city using Latin Hypercube Sampling with different parameter combinations, and explored the correlation between HIV-2 prevalence and estimated MC rate (both logit-transformed) in the 1,000 replicates., Results and Conclusions: Our survey reveals that, in the early 20th century, MC was far less common and geographically more variable than nowadays. HIV-2 prevalence in 1985-91 and MC rates in 1950 were negatively correlated (Spearman rho = -0.546, IQR: -0.553--0.546, p≤0.0021). Guinea-Bissau and Côte d'Ivoire cities had markedly lower MC rates. In addition, MC was uncommon in rural southwestern Côte d'Ivoire in 1930.The differential HIV-2 spread in West Africa correlates with different historical MC rates. We suggest HIV-2 only formed early substantial foci in cities with substantial uncircumcised populations. Lack of MC in rural areas exposed to bushmeat may have had a role in successful HIV-2 emergence., Competing Interests: The authors have declared that no competing interests exist.

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., (Copyright © 2016 The Authors. Published by Elsevier B.V. All rights reserved.)

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.

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.

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., (Copyright © 2015 Elsevier B.V. All rights reserved.)

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., (© The Author 2014. Published by Oxford University Press on behalf of the Infectious Diseases Society of America. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Rilpivirine is a second-generation nonnucleoside reverse-transcriptase inhibitor (NNRTI) currently indicated for first-line therapy, but its clinical benefit for HIV-1 infected patients failing first-generation NNRTIs is largely undefined. This study quantified the extent of genotypic rilpivirine resistance in viral isolates from 1212 patients upon failure of efavirenz- or nevirapine-containing antiretroviral treatment, of whom more than respectively 80% and 90% showed high-level genotypic resistance to the failing NNRTI. Of all study patients, 47% showed a rilpivirine resistance-associated mutation (RPV-RAM), whereas preserved residual rilpivirine activity was predicted in half of the patients by three genotypic drug resistance interpretation algorithms. An NNRTI-dependent impact on rilpivirine resistance was detected. Compared with the use of nevirapine, the use of efavirenz was associated with a 32% lower risk of having a RPV-RAM and a 50% lower risk of predicted reduced rilpivirine susceptibility. Most prevalent RPV-RAMs after nevirapine experience were Y181C and H221Y, whereas L100I+K103N, Y188L and K101E occurred most in efavirenz-experienced patients. Predicted rilpivirine activity was not affected by HIV-1 subtype, although frequency of individual mutations differed across subtypes. In conclusion, this genotypic resistance analysis strongly suggests that the latest NNRTI, rilpivirine, may retain activity in a large proportion of HIV-1 patients in whom resistance failed while they were on an efavirenz- or nevirapine-containing regimen, and may present an attractive option for second-line treatment given its good safety profile and dosing convenience. However, prospective clinical studies assessing the effectiveness of rilpivirine for NNRTI-experienced patients are warranted to validate knowledge derived from genotypic and phenotypic drug resistance studies., (Copyright © 2015 The Authors. Published by Elsevier Ltd.. All rights reserved.)

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., (Copyright © 2014 Elsevier B.V. All rights reserved.)

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.

Objectives: The amino acid position 70 in HIV-1 reverse transcriptase (RT) plays an important role in nucleoside RT inhibitor (NRTI) resistance. K70R is part of the thymidine analog mutations, but also other amino acid changes have been associated with NRTI resistance, such as K70E and K70G. In this study, we investigated the in vivo selection of the HIV-1 RT mutations K70S and K70T and their in vitro effect on drug resistance and replication capacity., Methods: Recombinant viruses with RT mutations were generated to measure the in vitro drug susceptibility and replication capacity. Bayesian network analysis and three-dimensional modeling were performed to understand the selection and impact of the RT70 mutations., Results: K70S and K70T were found at a low frequency in RTI-experienced HIV-1 patients (0.10% and 0·20%). Baeyesian network learning identified no direct association with the in vivo exposure to any specific RTI. However, direct associations of K70S with mutations within the Q151M-complex and of K70T with K65R were observed. In vitro phenotypic testing revealed only minor effects of K70R/S/T as single mutations, associated with Q151M and within the context of the Q151M-complex., Discussion: These results suggest that the selection of K70S/T and their phenotypic impact are influenced by the presence of other mutations in RT. However, the low impact on in vitro phenotype here observed, alongside with the low in vivo prevalence, the exclusive direct association with known major RTI mutations and the unknown correlation with in vivo response, do not yet necessitate the inclusion of K70S/T in drug resistance interpretation systems.

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.

To characterize the HIV-2 integrase gene polymorphisms and the pathways to resistance of HIV-2 patients failing a raltegravir-containing regimen, we studied 63 integrase strand transfer inhibitors (INSTI)-naïve patients, and 10 heavily pretreated patients exhibiting virological failure while receiving a salvage raltegravir-containing regimen. All patients were infected by HIV-2 group A. 61.4% of the integrase residues were conserved, including the catalytic motif residues. No INSTI-major resistance mutations were detected in the virus population from naïve patients, but two amino acids that are secondary resistance mutations to INSTIs in HIV-1 were observed. The 10 raltegravir-experienced patients exhibited resistance mutations via three main genetic pathways: N155H, Q148R, and eventually E92Q - T97A. The 155 pathway was preferentially used (7/10 patients). Other mutations associated to raltegravir resistance in HIV-1 were also observed in our HIV-2 population (V151I and D232N), along with several novel mutations previously unreported. Data retrieved from this study should help build a more robust HIV-2-specific algorithm for the genotypic interpretation of raltegravir resistance, and contribute to improve the clinical monitoring of HIV-2-infected patients.

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., (Copyright © 2013 Elsevier B.V. All rights reserved.)

Quantitation of HIV-1 RNA levels in plasma has significant prognostic value since high viral load concentrations in plasma are associated with a faster disease progression. Viral load testing became one the most important tools for monitoring HIV patients. New real time methodologies to quantify HIV viral load had arisen in the last decade. HIV is a virus with a high genetic variability, with the potential to negatively affect the performance of the viral load assays. Consequently, any new assay should be challenged against, at least, the most prevalent HIV-1 genetic variants. In the present study, the new version of NucliSENS EasyQ(®) HIV-1 (Version 2.0) quantitative assay was compared with another ultra-sensitive test--Abbott RealTime HIV-1--using 175 plasma samples from patients infected with several HIV-1 subtypes and recombinant forms: subtype B (41, 23%), subtype C (19, 11%), subtype G (76, 44%), and CRF02&#95;AG (39, 22%). Overall, there was agreement between the assays in 95.43% of the samples. Both assays have a very good dynamic range [1.4-6.9] and [1.60-7.0] log10 copies/mL and excellent correlation in samples with various subtypes. Based on the fact that no clinically significant differences were observed in the viral load measurements by these two assays, HIV-1 subtypes are quantified equally by both assays. However due to HIV diversity, mainly in regions were non B subtypes are predominant more evaluations are needed, so we do not recommend to switch platform during longitudinal viral load monitoring., (Copyright © 2013 Elsevier B.V. All rights reserved.)

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&#95;AE was ~99% with statistical-based tools, with phylogenetic-based tools and with Stanford, one of the similarity based tools. CRF02&#95;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&#95;AE and CRF02&#95;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&#95;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., (Copyright © 2013 The Authors. Published by Elsevier B.V. All rights reserved.)

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.

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.

Objectives: Despite a decreasing mortality and morbidity in treated HIV-1 patients, highly active antiretroviral treatment (HAART) can still fail due to the development of drug resistance. Especially, multidrug-resistant viruses pose a threat to efficient therapy. We studied the changing prevalence of multidrug resistance (MDR) over time in a cohort of HIV-1-infected patients in Portugal., Patients and Methods: We used data of 8065 HIV-1-infected patients followed from July 2001 up to April 2012 in 22 hospitals located in Portugal. MDR at a specific date of sampling was defined as no more than one fully active drug (excluding integrase and entry inhibitors) at that time authorized by the Portuguese National Authority of Medicines and Health Products (INFARMED), as interpreted with the Rega algorithm version 8.0.2. A generalized linear mixed model was used to study the time trend of the prevalence of MDR., Results: We observed a statistically significant decrease in the prevalence of MDR over the last decade, from 6.9% (95% CI: 5.7-8.4) in 2001-03, 6.0% (95% CI: 4.9-7.2) in 2003-05, 3.7% (95% CI: 2.8-4.8) in 2005-07 and 1.6% (95% CI: 1.1-2.2) in 2007-09 down to 0.6% (95% CI: 0.3-0.9) in 2009-12 [OR=0.80 (95% CI: 0.75-0.86); P<0.001]. In July 2011 the last new case of MDR was seen., Conclusions: The prevalence of multidrug-resistant HIV-1 is decreasing over time in Portugal, reflecting the increasing efficiency of HAART and the availability of new drugs. Therefore, in designing a new drug, safety and practical aspects, e.g. less toxicity and ease of use, may need more attention than focusing mainly on efficacy against resistant strains.

Subtype-dependent selection of HIV-1 reverse transcriptase resistance mutation K65R was previously observed in cell culture and small clinical investigations. We compared K65R prevalence across subtypes A, B, C, F, G, and CRF02&#95;AG separately in a cohort of 3,076 patients on combination therapy including tenofovir. K65R selection was significantly higher in HIV-1 subtype C. This could not be explained by clinical and demographic factors in multivariate analysis, suggesting subtype sequence-specific K65R pathways.

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&#95;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.

Objectives: To investigate mutations selected in viruses from HIV-2-infected patients failing a highly active antiretroviral treatment (HAART) regimen including atazanavir/ritonavir., Methods: Twenty-eight HIV-2-infected patients previously exposed to atazanavir/ritonavir and failing therapy were studied. The protease (PR) gene was amplified and sequenced, and mutations emerging under atazanavir/ritonavir selective pressure were reported., Results: The I50L mutation emerged in 4 out of 28 HIV-2-infected patients failing a HAART regimen including atazanavir/ritonavir. Besides I50L, four PR mutations previously associated with protease inhibitor resistance (I54L, I64V, V71I and I82F) and six PR mutations of unknown impact (V10I, E37D, S43T, K45R, I75V and F85L) in HIV-2 were also identified in this small group of patients., Conclusions: Several mutations were associated with virological failure of a regimen including atazanavir/ritonavir in HIV-2-infected patients, including I50L for the first time. It should be included in HIV-2 algorithms for interpretation of genotypic resistance data, and taken into account when making therapeutic decisions for HIV-2-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.

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&#95;AG. These HS-linked polymorphisms were introduced in the BD6-15 CRF02&#95;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&#95;AG. In direct head-to-head virus competitions, CRF02&#95;AG harboring the 17E/64M polymorphisms also had higher replicative fitness than did the 17E or the 64M polymorphism in the CFR02&#95;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.

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.

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&#95;AG clade. By explicitly integrating prevalence counts and genetic population size estimates we date the epidemic emergence of CRF02&#95;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&#95;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&#95;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&#95;AG at a finer geographical detail. Phylogeographic dispersal in continuous space reveals that at least two distinct CRF02&#95;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&#95;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., (Copyright © 2011 Elsevier B.V. All rights reserved.)

HIV-2 is responsible for a limited epidemic in West Africa. Around 20% of all infected patients will progress to AIDS, and will need antiretroviral therapy. Unfortunately, antiretrovirals were developed to suppress HIV-1 replication; not all of them are active against HIV-2, e.g. all nonnucleoside reverse transcriptase inhibitors or fusion inhibitors. Moreover, only three protease inhibitors have the same activity in HIV-1 and HIV-2: lopinavir, saquinavir and darunavir. Even if all nucleoside and nucleotide reverse transcriptase inhibitors appear to be equally efficient against HIV-2, different resistance pathways and an increased facility of resistance selection make their use much more difficult than in HIV-1. Integrase inhibitors have a potent inhibitory effect on HIV-2 replication, but questions about the best timing for their use remain unanswered, as well as those regarding the use of entry inhibitors in this setting. The lack of reliable monitoring tools adds to the difficulty of treating HIV-2-infected patients, mostly because the viral load is not as useful as it is in HIV-1, and the incomplete knowledge about resistance pathways limits the clinical usefulness of resistance testing. With all these limitations, HIV-2 treatment remains a challenge. Further research is urgently needed, since antiretroviral therapy is now becoming available in countries where the HIV-2 prevalence is significant. The need for appropriate guidelines for HIV-2 treatment has become an emergency., (Copyright © 2012 S. Karger AG, Basel.)

The advent of new therapeutic approaches targeting env and the search for efficient anti-HIV-1 vaccines make it necessary to identify the number of recombinant forms using genomic regions that were previously not frequently sequenced. In this study, we have subtyped paired pol and env sequences from HIV-1 strains infecting 152 patients being clinically followed in Portugal. The percentage of strains in which we found discordant subtypes in pol and env was 25.7%. When the subtype in pol and env was concordant (65.1%), the most prevalent subtypes were subtype B (40.8%), followed by subtype C (17.8%) and subtype G (5.3%). The most prevalent recombinant form was CRF14&#95;BGpol/Genv (7.2%).

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.

Background: Failure on Highly Active Anti-Retroviral Treatment is often accompanied with development of antiviral resistance to one or more drugs included in the treatment. In general, the virus is more likely to develop resistance to drugs with a lower genetic barrier. Previously, we developed a method to reverse engineer, from clinical sequence data, a fitness landscape experienced by HIV-1 under nelfinavir (NFV) treatment. By simulation of evolution over this landscape, the individualized genetic barrier to NFV resistance may be estimated for an isolate., Results: We investigated the association of estimated genetic barrier with risk of development of NFV resistance at virological failure, in 201 patients that were predicted fully susceptible to NFV at baseline, and found that a higher estimated genetic barrier was indeed associated with lower odds for development of resistance at failure (OR 0.62 (0.45 - 0.94), per additional mutation needed, p = .02)., Conclusions: Thus, variation in individualized genetic barrier to NFV resistance may impact effective treatment options available after treatment failure. If similar results apply for other drugs, then estimated genetic barrier may be a new clinical tool for choice of treatment regimen, which allows consideration of available treatment options after virological failure.

A better understanding of human immunodeficiency virus type 1 drug-resistance evolution under the selective pressure of combination treatment is important for the design of long-term effective treatment strategies. We applied Bayesian network learning to sequences from patients treated with the reverse transcriptase inhibitor combination of zidovudine (AZT) and lamivudine (3TC) to identify the role of many treatment-selected mutations in the development of resistance. Based on the Bayesian network structure, an in vivo fitness landscape was built, reflecting the necessary selective pressure under treatment, to evolve naive sequences to sequences obtained from patients treated with the combination. This landscape, combined with an evolutionary model, was used to predict resistance evolution in longitudinal sequence pairs. In our analysis, mutations 41L, 70R, 184V and 215F/Y were identified as major resistance mutations to the combination of AZT and 3TC, as they were associated directly with treatment experience. The network also suggested a possible role in resistance development for a number of novel mutations. Estimated fitness, using the landscape, correlated significantly with in vitro resistance phenotype in genotype-phenotype pairs (R(2)=0.70). Variation in predicted evolution under selective pressure correlated significantly with observed in vivo evolution during AZT plus 3CT treatment. In conclusion, we confirmed current knowledge on resistance development to the combination of AZT and 3CT, but additional novel mutations were identified. Moreover, a model to predict resistance evolution during AZT and 3CT treatment has been built and validated.

The REGA HIV-1 subtyping tool is a phylogenetic-based method for subtyping HIV-1 genomic sequences that was published in 2005. The subtyping tool combines phylogenetic approaches with recombination detection methods. Recently, version 2 was released (http://www.bioafrica.net/rega-genotype/html/index.html) as an improvement of version 1. Version 2 implements a Decision-Tree-based algorithm that was not implemented in version 1. We wanted to compare the two versions on a large sequence dataset to assess the improvements of version 2 and to verify whether features lost during updating the tool needed to be recovered. We analysed the results of the two versions in the genotyping of 4676 HIV-1 pol sequences. We compared those results to a manual approach, used in previous studies. Our results show that version 2 has an overall better sensitivity but especially for the detection of subtypes A, B, D, F, G and CRF14&#95;BG and CRF06&#95;CPX. For the other subtypes, no significant differences were observed in the sensitivity of versions 1 and 2. The overall increase in sensitivity was however accompanied by a decrease in the specificity for the detection of subtype B. This is the main limitation of version 2. However, while the number of false negatives decreased by 53 samples, the number of false positives increased only by 5 samples from version 1 to 2. The performance of the REGA HIV-1 subtyping tool was considerably improved from one version to the other. Our results are very valuable and allow us to make suggestions for further improvement of the tool for a version 3 release., (2009 Elsevier B.V. All rights reserved.)