Search

Background: Numerous studies have shown that baseline drug resistance patterns may influence the outcome of antiretroviral therapy. Therefore, guidelines recommend drug resistance testing to guide the choice of initial regimen. In addition to optimizing individual patient management, these baseline resistance data enable transmitted drug resistance (TDR) to be surveyed for public health purposes. The SPREAD program systematically collects data to gain insight into TDR occurring in Europe since 2001., Methods: Demographic, clinical, and virological data from 4140 antiretroviral-naive human immunodeficiency virus (HIV)-infected individuals from 26 countries who were newly diagnosed between 2008 and 2010 were analyzed. Evidence of TDR was defined using the WHO list for surveillance of drug resistance mutations. Prevalence of TDR was assessed over time by comparing the results to SPREAD data from 2002 to 2007. Baseline susceptibility to antiretroviral drugs was predicted using the Stanford HIVdb program version 7.0., Results: The overall prevalence of TDR did not change significantly over time and was 8.3% (95% confidence interval, 7.2%-9.5%) in 2008-2010. The most frequent indicators of TDR were nucleoside reverse transcriptase inhibitor (NRTI) mutations (4.5%), followed by nonnucleoside reverse transcriptase inhibitor (NNRTI) mutations (2.9%) and protease inhibitor mutations (2.0%). Baseline mutations were most predictive of reduced susceptibility to initial NNRTI-based regimens: 4.5% and 6.5% of patient isolates were predicted to have resistance to regimens containing efavirenz or rilpivirine, respectively, independent of current NRTI backbones., Conclusions: Although TDR was highest for NRTIs, the impact of baseline drug resistance patterns on susceptibility was largest for NNRTIs. The prevalence of TDR assessed by epidemiological surveys does not clearly indicate to what degree susceptibility to different drug classes is affected., (© The Author 2015. Published by Oxford University Press for the Infectious Diseases Society of America.)

Genetic analyses play a central role in infectious disease research. Massively parallelized "mechanical cloning" and sequencing technologies were quickly adopted by HIV researchers in order to broaden the understanding of the clinical importance of minor drug-resistant variants. These efforts have, however, remained largely limited to small genomic regions. The growing need to monitor multiple genome regions for drug resistance testing, as well as the obvious benefit for studying evolutionary and epidemic processes makes complete genome sequencing an important goal in viral research. In addition, a major drawback for NGS applications to RNA viruses is the need for large quantities of input DNA. Here, we use a generic overlapping amplicon-based near full-genome amplification protocol to compare low-input enzymatic fragmentation (Nextera™) with conventional mechanical shearing for Roche 454 sequencing. We find that the fragmentation method has only a modest impact on the characterization of the population composition and that for reliable results, the variation introduced at all steps of the procedure--from nucleic acid extraction to sequencing--should be taken into account, a finding that is also relevant for NGS technologies that are now more commonly used. Furthermore, by applying our protocol to deep sequence a number of pre-therapy plasma and PBMC samples, we illustrate the potential benefits of a near complete genome sequencing approach in routine genotyping.

For many years, population-based Sanger sequencing has been the golden standard for drug resistance testing within the routine follow-up of HIV-1 infected patients in resource-rich settings. Often, the data generated within this framework were subsequently used for research and surveillance purposes: to understand therapy response and to gain insights into epidemiological processes. Sanger sequencing was however ill suited for diagnostic and prognostic use in resource-limited settings (RLS) and therefore not broadly implemented. Next-generation sequencing (NGS) technologies provide high-throughput approaches by the rapid acquisition of thousands to millions of short nucleotide sequences. Depending on the experimental design, the roll-out of NGS drug resistance testing at a larger scale is feasible, providing better characterization and understanding of the evolving population of viral variants within a patient and potentially improving the prognostic value of drug resistance testing. Whether the same will become true for RLS will largely depend on affordability and sample logistics, and this may affect other mutation-specific approaches., (Copyright © 2015 Elsevier B.V. All rights reserved.)

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.

Treatment with pan-genotypic direct-acting antivirals, targeting different viral proteins, is the best option for clearing hepatitis C virus (HCV) infection in chronically infected patients. However, the diversity of the HCV genome is a major obstacle for the development of antiviral drugs, vaccines, and genotyping assays. In this large-scale analysis, genome-wide diversity and selective pressure was mapped, focusing on positions important for treatment, drug resistance, and resistance testing. A dataset of 1415 full-genome sequences, including genotypes 1-6 from the Los Alamos database, was analyzed. In 44% of all full-genome positions, the consensus amino acid was different for at least one genotype. Focusing on positions sharing the same consensus amino acid in all genotypes revealed that only 15% was defined as pan-genotypic highly conserved (≥99% amino acid identity) and an additional 24% as pan-genotypic conserved (≥95%). Despite its large genetic diversity, across all genotypes, codon positions were rarely identified to be positively selected (0.23%-0.46%) and predominantly found to be under negative selective pressure, suggesting mainly neutral evolution. For NS3, NS5A, and NS5B, respectively, 40% (6/15), 33% (3/9), and 14% (2/14) of the resistance-related positions harbored as consensus the amino acid variant related to resistance, potentially impeding treatment. For example, the NS3 variant 80K, conferring resistance to simeprevir used for treatment of HCV1 infected patients, was present in 39.3% of the HCV1a strains and 0.25% of HCV1b strains. Both NS5A variants 28M and 30S, known to be associated with resistance to the pan-genotypic drug daclatasvir, were found in a significant proportion of HCV4 strains (10.7%). NS5B variant 556G, known to confer resistance to non-nucleoside inhibitor dasabuvir, was observed in 8.4% of the HCV1b strains. Given the large HCV genetic diversity, sequencing efforts for resistance testing purposes may need to be genotype-specific or geographically tailored.

Objective: To determine how HIV-1 risk groups impact transmitted diversity and the tempo of viral evolution at a population scale., Methods: We investigated a set of previously described transmission chains (n = 70) using a population genetic approach, and tested whether the expected differences in proportions of multivariant transmissions are reflected by varying proportions of transmitted diversity between men having sex with men (MSM) and heterosexual (HET) subpopulations - the largest contributors to HIV spread. To assess evolutionary rate differences among the different risk groups, we compiled risk group datasets for subtypes A1, B and CRF01_AE, and directly compared the absolute substitution rate and its synonymous and non-synonymous components., Results: There was sufficient demographic signal to inform the transmission model in Bayesian evolutionary analysis by sampling trees using env data to compare the transmission bottleneck size between the MSM and HET risk groups. We found no indications for a different proportion of transmitted genetic diversity at the population level between these groups. In the direct rate comparisons between the risk groups, however, we consistently recovered a higher evolutionary rate in the male-dominated risk group compared to the HET datasets., Conclusion: We find that the risk group composition affects the viral evolutionary rate and therefore potentially also the adaptation rate. In particular, risk group-specific sex ratios, and the variation in within-host evolutionary rates between men and women, impose evolutionary rate differences at the epidemic level, but we cannot exclude a role of varying transmission rates.

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.)

Background: Regional and subtype-specific mutational patterns of HIV-1 transmitted drug resistance (TDR) are essential for informing first-line antiretroviral (ARV) therapy guidelines and designing diagnostic assays for use in regions where standard genotypic resistance testing is not affordable. We sought to understand the molecular epidemiology of TDR and to identify the HIV-1 drug-resistance mutations responsible for TDR in different regions and virus subtypes., Methods and Findings: We reviewed all GenBank submissions of HIV-1 reverse transcriptase sequences with or without protease and identified 287 studies published between March 1, 2000, and December 31, 2013, with more than 25 recently or chronically infected ARV-naïve individuals. These studies comprised 50,870 individuals from 111 countries. Each set of study sequences was analyzed for phylogenetic clustering and the presence of 93 surveillance drug-resistance mutations (SDRMs). The median overall TDR prevalence in sub-Saharan Africa (SSA), south/southeast Asia (SSEA), upper-income Asian countries, Latin America/Caribbean, Europe, and North America was 2.8%, 2.9%, 5.6%, 7.6%, 9.4%, and 11.5%, respectively. In SSA, there was a yearly 1.09-fold (95% CI: 1.05-1.14) increase in odds of TDR since national ARV scale-up attributable to an increase in non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance. The odds of NNRTI-associated TDR also increased in Latin America/Caribbean (odds ratio [OR] = 1.16; 95% CI: 1.06-1.25), North America (OR = 1.19; 95% CI: 1.12-1.26), Europe (OR = 1.07; 95% CI: 1.01-1.13), and upper-income Asian countries (OR = 1.33; 95% CI: 1.12-1.55). In SSEA, there was no significant change in the odds of TDR since national ARV scale-up (OR = 0.97; 95% CI: 0.92-1.02). An analysis limited to sequences with mixtures at less than 0.5% of their nucleotide positions—a proxy for recent infection—yielded trends comparable to those obtained using the complete dataset. Four NNRTI SDRMs—K101E, K103N, Y181C, and G190A—accounted for >80% of NNRTI-associated TDR in all regions and subtypes. Sixteen nucleoside reverse transcriptase inhibitor (NRTI) SDRMs accounted for >69% of NRTI-associated TDR in all regions and subtypes. In SSA and SSEA, 89% of NNRTI SDRMs were associated with high-level resistance to nevirapine or efavirenz, whereas only 27% of NRTI SDRMs were associated with high-level resistance to zidovudine, lamivudine, tenofovir, or abacavir. Of 763 viruses with TDR in SSA and SSEA, 725 (95%) were genetically dissimilar; 38 (5%) formed 19 sequence pairs. Inherent limitations of this study are that some cohorts may not represent the broader regional population and that studies were heterogeneous with respect to duration of infection prior to sampling., Conclusions: Most TDR strains in SSA and SSEA arose independently, suggesting that ARV regimens with a high genetic barrier to resistance combined with improved patient adherence may mitigate TDR increases by reducing the generation of new ARV-resistant strains. A small number of NNRTI-resistance mutations were responsible for most cases of high-level resistance, suggesting that inexpensive point-mutation assays to detect these mutations may be useful for pre-therapy screening in regions with high levels of TDR. In the context of a public health approach to ARV therapy, a reliable point-of-care genotypic resistance test could identify which patients should receive standard first-line therapy and which should receive a protease-inhibitor-containing regimen.

Background: The HIV-1 epidemic in Slovenia, a small Central European country, has some characteristics that make it an ideal model to study HIV-1 transmission. The epidemic is predominantly affecting men who have sex with men infected with subtype B (89% of all patients), has a low prevalence (less than 1/1000) and is growing slowly. The aim of the present study was to analyze in detail the evolutionary history and the determinants of transmission., Methods: A total of 223 partial pol gene sequences from therapy naïve individuals were included, representing 52% of all patients newly diagnosed in 13 years (2000-2012) and analyzed together with genetically similar worldwide sequences, selected in a BLAST search., Results: Combined analysis (maximum likelihood and Bayesian) of HIV-1 transmission chains revealed 8 major clusters (n ≥ 10 patients), 1 group of 4 patients, 2 trios and 12 transmission pairs, thus leaving only 43 (19.3%) Slovenian patients infected with subtype B without a local epidemiological link, indicating a predominance of local transmission of HIV-1 infection. Bayesian analysis performed on a full set of sequences estimated several introductions of HIV-1 into Slovenia, with the most recent common ancestor (tMRCA) of the earliest Slovenian cluster dated to the late 1980s, although tMRCAs obtained from separate independent analysis of each cluster showed considerably more recent estimates. These findings indicate inconsistencies in molecular clock estimation, which we further explored. We hypothesize that these inconsistent dating estimates across the tree could be caused by an evolutionary rate acceleration of HIV-1 after entering the Slovenia epidemic that is not taken into account by the molecular clock model. It could be caused by a lower transmission rate in this setting, as demonstrated by the low epidemic growth rate estimated by Bayesian skyline plot analysis., Conclusions: HIV-1 subtype B was introduced into Slovenia at several time points from the late 80s onward. The majority of patients had a local transmission link, indicating a closed HIV community. The observed slower epidemic rate suggests that individuals with a long-lasting infection are the driving force of the epidemic in this region.

Background: The HIV pandemic is characterized by extensive genetic variability, which has challenged the development of HIV drugs and vaccines. Although HIV genomes have been classified into different types, groups, subtypes and recombinants, a comprehensive study that maps HIV genome-wide diversity at the population level is still lacking to date. This study aims to characterize HIV genomic diversity in large-scale sequence populations, and to identify driving factors that shape HIV genome diversity., Results: A total of 2996 full-length genomic sequences from 1705 patients infected with 16 major HIV groups, subtypes and circulating recombinant forms (CRFs) were analyzed along with structural, immunological and peptide inhibitor information. Average nucleotide diversity of HIV genomes was almost 50% between HIV-1 and HIV-2 types, 37.5% between HIV-1 groups, 14.7% between HIV-1 subtypes, 8.2% within individual HIV-1 subtypes and less than 1% within single patients. Along the HIV genome, diversity patterns and compositions of nucleotides and amino acids were highly similar across different groups, subtypes and CRFs. Current HIV-derived peptide inhibitors were predominantly derived from conserved, solvent accessible and intrinsically ordered structures in the HIV-1 subtype B genome. We identified these conserved regions in Capsid, Nucleocapsid, Protease, Integrase, Reverse transcriptase, Vpr and the GP41 N terminus as potential drug targets. In the analysis of factors that impact HIV-1 genomic diversity, we focused on protein multimerization, immunological constraints and HIV-human protein interactions. We found that amino acid diversity in monomeric proteins was higher than in multimeric proteins, and diversified positions were preferably located within human CD4 T cell and antibody epitopes. Moreover, intrinsic disorder regions in HIV-1 proteins coincided with high levels of amino acid diversity, facilitating a large number of interactions between HIV-1 and human proteins., Conclusions: This first large-scale analysis provided a detailed mapping of HIV genomic diversity and highlighted drug-target regions conserved across different groups, subtypes and CRFs. Our findings suggest that, in addition to the impact of protein multimerization and immune selective pressure on HIV-1 diversity, HIV-human protein interactions are facilitated by high variability within intrinsically disordered structures.

As commercial human immunodeficiency virus type 1 drug resistance assays are expensive, they are not commonly used in resource-limited settings. Hence, a more affordable in-house procedure was set up taking into account the specific epidemiological and economic circumstances of Cuba. The performance characteristics of the in-house assay were evaluated using clinical samples with various subtypes and resistance patterns. The lower limit of amplification was determined on dilutions series of 20 clinical isolates and ranged from 84 to 529 RNA copies/mL. For the assessment of trueness, 14 clinical samples were analyzed and the ViroSeq HIV-1 Genotyping System v2.0 was used as the reference standard. The mean nucleotide sequence identity between the two assays was 98.7% ± 1.0. Additionally, 99.0% of the amino acids at drug resistance positions were identical. The sensitivity and specificity in detecting drug resistance mutations was respectively 94.1% and 99.5%. Only few discordances in drug resistance interpretation patterns were observed. The repeatability and reproducibility were evaluated using 10 clinical samples with 3 replicates per sample. The in-house test was very precise as nucleotide sequence identity among paired nucleotide sequences ranged from 98.7% to 99.9%. The acceptance criteria were met by the in-house test for all performance characteristics, demonstrating a high degree of accuracy. Subsequently, the applicability in routine clinical practice was evaluated on 380 plasma samples. The amplification success rate was 91% and good quality consensus sequences encoding the entire protease and the first 335 codons in reverse transcriptase could be obtained for 99% of the successful amplicons. The reagent cost per sample using the in-house procedure was around € 80 per genotyping attempt. Overall, the in-house assay provided good results, was feasible with equipment and reagents available in Cuba and was half as expensive as commercial assays.

Background: Clinicians reported an increasing trend of rapid progression (RP) (AIDS within 3 years of infection) in Cuba., Methods: Recently infected patients were prospectively sampled, 52 RP at AIDS diagnosis (AIDS-RP) and 21 without AIDS in the same time frame (non-AIDS). 22 patients were sampled at AIDS diagnosis (chronic-AIDS) retrospectively assessed as > 3 years infected. Clinical, demographic, virological, epidemiological and immunological data were collected. Pol and env sequences were used for subtyping, transmission cluster analysis, and prediction of resistance, co-receptor use and evolutionary fitness. Host, immunological and viral predictors of RP were explored through data mining., Findings: Subtyping revealed 26 subtype B strains, 6 C, 6 CRF18_cpx, 9 CRF19_cpx, 29 BG-recombinants and other subtypes/URFs. All patients infected with CRF19 belonged to the AIDS-RP group. Data mining identified CRF19, oral candidiasis and RANTES levels as the strongest predictors of AIDS-RP. CRF19 was more frequently predicted to use the CXCR4 co-receptor, had higher fitness scores in the protease region, and patients had higher viral load at diagnosis., Interpretation: CRF19 is a recombinant of subtype D (C-part of Gag, PR, RT and nef), subtype A (N-part of Gag, Integrase, Env) and subtype G (Vif, Vpr, Vpu and C-part of Env). Since subtypes D and A have been associated with respectively faster and slower disease progression, our findings might indicate a fit PR driving high viral load, which in combination with co-infections may boost RANTES levels and thus CXCR4 use, potentially explaining the fast progression. We propose that CRF19 is evolutionary very fit and causing rapid progression to AIDS in many newly infected patients in Cuba.

Unlabelled: Human T-cell lymphotropic virus type 1 (HTLV-1) and HTLV-2 encode auxiliary proteins that play important roles in viral replication, viral latency, and immune escape. The presence of auxiliary protein-encoding open reading frames (ORFs) in HTLV-3, the latest HTLV to be discovered, is unknown. Simian T-cell lymphotropic virus type 3 (STLV-3) is almost identical to HTLV-3. Given the lack of HTLV-3-infected cell lines, we took advantage of STLV-3-infected cells and of an STLV-3 molecular clone to search for the presence of auxiliary transcripts. Using reverse transcriptase PCR (RT-PCR), we first uncovered the presence of three unknown viral mRNAs encoding putative proteins of 5, 8, and 9 kDa and confirmed the presence of the previously reported RorfII transcript. The existence of these viral mRNAs was confirmed by using splice site-specific RT-PCR with ex vivo samples. We showed that p5 is distributed throughout the cell and does not colocalize with a specific organelle. The p9 localization is similar to that of HTLV-1 p12 and induced a strong decrease in the calreticulin signal, similarly to HTLV-1 p12. Although p8, RorfII, and Rex-3 share an N-terminal sequence that is predicted to contain a nucleolar localization signal (NoLS), only p8 is found in the nucleolus. The p8 location in the nucleolus is linked to a bipartite NoLS. p8 and, to a lesser extent, p9 repressed viral expression but did not alter Rex-3-dependent mRNA export. Using a transformation assay, we finally showed that none of the STLV-3 auxiliary proteins had the ability to induce colony formation, while both Tax-3 and antisense protein of HTLV-3 (APH-3) promoted cellular transformation. Altogether, these results complete the characterization of the newly described primate T-lymphotropic virus type 3 (PTLV-3)., Importance: Together with their simian counterparts, HTLVs form the primate T-lymphotropic viruses. HTLVs arose from interspecies transmission between nonhuman primates and humans. HTLV-1 and HTLV-2 encode auxiliary proteins that play important roles in viral replication, viral latency, and immune escape. The presence of ORFs encoding auxiliary proteins in HTLV-3 or STLV-3 genomes was unknown. Using in silico analyses, ex vivo samples, or in vitro experiments, we have uncovered the presence of 3 previously unknown viral mRNAs encoding putative proteins and confirmed the presence of a previously reported viral transcript. We characterized the intracellular localization of the four proteins. We showed that two of these proteins repress viral expression but that none of them have the ability to induce colony formation. However, both Tax and the antisense protein APH-3 promote cell transformation. Our results allowed us to characterize 4 new retroviral proteins for the first time., (Copyright © 2015, American Society for Microbiology. All Rights Reserved.)

Background: A key challenge in the field of HIV-1 protein evolution is the identification of coevolving amino acids at the molecular level. In the past decades, many sequence-based methods have been designed to detect position-specific coevolution within and between different proteins. However, an ensemble coevolution system that integrates different methods to improve the detection of HIV-1 protein coevolution has not been developed., Results: We integrated 27 sequence-based prediction methods published between 2004 and 2013 into an ensemble coevolution system. This system allowed combinations of different sequence-based methods for coevolution predictions. Using HIV-1 protein structures and experimental data, we evaluated the performance of individual and combined sequence-based methods in the prediction of HIV-1 intra- and inter-protein coevolution. We showed that sequence-based methods clustered according to their methodology, and a combination of four methods outperformed any of the 27 individual methods. This four-method combination estimated that HIV-1 intra-protein coevolving positions were mainly located in functional domains and physically contacted with each other in the protein tertiary structures. In the analysis of HIV-1 inter-protein coevolving positions between Gag and protease, protease drug resistance positions near the active site mostly coevolved with Gag cleavage positions (V128, S373-T375, A431, F448-P453) and Gag C-terminal positions (S489-Q500) under selective pressure of protease inhibitors., Conclusions: This study presents a new ensemble coevolution system which detects position-specific coevolution using combinations of 27 different sequence-based methods. Our findings highlight key coevolving residues within HIV-1 structural proteins and between Gag and protease, shedding light on HIV-1 intra- and inter-protein coevolution.

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.)

The HIV Rev protein mediates the transport of partially and unspliced HIV mRNA from the nucleus to the cytoplasm. Rev multimerizes on a secondary stem-loop structure present in the viral intron-containing mRNA species and recruits the cellular karyopherin CRM1 to export viral mRNAs from the nucleus to the cytoplasm. Previously we have identified a single-domain intrabody (Nb(190)), derived from a llama heavy-chain antibody, which efficiently inhibits Rev multimerization and suppresses the production of infectious virus. We recently mapped the epitope of this nanobody and demonstrated that Rev residues K20 and Y23 are crucial for interaction while residues V16, H53 and L60 are important to a lesser extent. Here, we generated cell lines stably expressing Nb(190) and assessed the capacity of these cell lines to suppress the replication of different HIV-1 subtypes. These cells stably expressing the single-domain antibody are protected from virus-induced cytopathogenic effect even in the context of high multiplicity of infection. In addition, the replication of different subtypes of group M and one strain of group O is significantly suppressed in these cell lines. Next, we analysed the natural variations of Rev amino acids in sequence samples from HIV-1 infected patients worldwide and assessed the effect of Nb(190) on the most prevalent polymorphisms occurring at the key epitope positions (K20 and Y23) in Rev. We found that Nb(190) was able to suppress the function of these Rev variants except for the K20N mutant, which was present in only 0.7% of HIV-1 sequence populations (n = 4632). Cells stably expressing the single-domain intrabody Nb(190) are protected against virus-induced cytopathogenic effect and display a selective survival advantage upon infection. In addition, Nb(190) suppresses the replication of a wide range of different HIV-1 subtypes. Large-scale sequence analysis reveals that the Nb(190) epitope positions in Rev are well conserved across major HIV-1 subtypes and groups. Altogether, our results indicate that Nb(190) may have broad potential as a gene therapeutic agent against HIV-1., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Introduction: Emergence of HIV-1 drug resistance may limit the sustained benefits of antiretroviral therapy (ART) in settings with limited laboratory monitoring and drug options. The objective is to implement the surveillance of drug resistance and subtypes in HIV-1 patients failing ART in Cuba., Methods: This study compiled clinical and genotypic drug resistance data 588 ART-experienced HIV-1 patients attending a clinical center in Havana in 2009-2013. Drug resistance testing was performed as part of routine clinical care. Drug resistance mutations and levels were determined using Rega version 8.0.2., Results: Eighty-three percent received solely ART containing at least three drugs. Patients from 2009 to 2010 were longer treated (median: 4.9 vs 2.7 years) and exposed to more ART regimens (median: 4 vs 2 regimens) compared to patients from 2011-2013. Nucleoside reverse transcriptase inhibitor (NRTI), non-nucleoside RTI (NNRTI) and PI mutations were present in 83.5, 77.4 and 52.0%. Full-class resistance (FCR) to NRTI, NNRTI, PI and multidrug resistance (MDR) were detected in 25.0, 33.7, 11.4 and 6.3%. FCR to NRTI, NNRTI, PI and MDR were present in 12.8, 28.7, 0 and 0% after first-line failure (164 patients) and in 23.1, 34.6, 3.8 and 3.1% after second-line failure (130 patients). Subtype B (32.5%), BG recombinants (19.6%) and CRF19_cpx (16.2%) were the most prevalent genetic forms. Subtype distribution did not change significantly between 2009-2010 and 2011-2013, except for BG recombinants that increased from 12.2 to 21.3% (p=0.002)., Conclusions: Our study found a high prevalence of drug resistance and supports the need for appropriate laboratory monitoring in clinical practice and access to drug options in case of virological failure.

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.

HIV-1 Gag amino acid substitutions associated with protease inhibitor (PI) treatment have mainly been reported in subtype B, while information on other subtypes is scarce. Using sequences from 11613 patients infected with different HIV-1 subtypes, we evaluated the prevalence of 93 Gag amino acid substitutions and their association with genotypic PI resistance. A significant association was found for 13 Gag substitutions, including A431V in both subtype B and CRF01&#95;AE. K415R in subtype C and S451G in subtype B were newly identified. Most PI-associated Gag substitutions are located in the flexible C-terminal domain, revealing the key role this region plays in PI resistance.

Introduction: In Tanzania, the follow-up on antiretroviral therapy (ART) response is based on clinical outcomes. We investigated virological response and ARV resistance mutations in relation to clinical response in ARV-treated patients., Methodology: A cross-sectional study of a cohort of 150 patients taking first-line ART in Dar-es-Salaam was conducted. Data were collected using standardized questionnaires and patients' blood samples. HIV viral load testing and genotyping was performed on all viremic samples. Statistical analyses compared clinical responders and non-responders., Results: The median time on ART was 20 months; 71 (47%) patients were ART clinical responders. Clinical non-responders were more likely to have started ART with advanced disease with significantly lower median percentage weight gain (6% versus 20%) with respect to pre-treatment levels. Sixty-one (86%) and 64 (81%) of clinical responders and non-responders, respectively, had undetectable viral loads. Genotyping was successful in 24 (96%) virologically failing patients, among whom 83% had resistance mutations; 67% had dual nucleoside reverse transcriptase inhibitor (NRTI)/non-NRTI (NNRTI) resistance mutations. Seventeen (71%) and 19 (79%) patients had NRTI and NNRTI resistance mutations, respectively, which were related to the ART in use, with no difference between clinical responders and non-responders. The most prevalent subtypes were A and C, found in 9 (38%) and 7 (29%) patients, respectively., Conclusions: The observed virological response was high and did not correlate with clinical response. The prevalence of ARV resistance mutations was high in viraemic patients and was related to the ARV prescribed. We recommend use of viral load monitoring during ART in Tanzania.

We aimed to study epidemic trends and predictors for transmitted drug resistance (TDR) in our region, its clinical impact and its association with transmission clusters. We included 778 patients from the AIDS Reference Center in Leuven (Belgium) diagnosed from 1998 to 2012. Resistance testing was performed using population-based sequencing and TDR was estimated using the WHO-2009 surveillance list. Phylogenetic analysis was performed using maximum likelihood and Bayesian techniques. The cohort was predominantly Belgian (58.4%), men who have sex with men (MSM) (42.8%), and chronically infected (86.5%). The overall TDR prevalence was 9.6% (95% confidence interval (CI): 7.7-11.9), 6.5% (CI: 5.0-8.5) for nucleoside reverse transcriptase inhibitors (NRTI), 2.2% (CI: 1.4-3.5) for non-NRTI (NNRTI), and 2.2% (CI: 1.4-3.5) for protease inhibitors. A significant parabolic trend of NNRTI-TDR was found (p = 0.019). Factors significantly associated with TDR in univariate analysis were male gender, Belgian origin, MSM, recent infection, transmission clusters and subtype B, while multivariate and Bayesian network analysis singled out subtype B as the most predictive factor of TDR. Subtype B was related with transmission clusters with TDR that included 42.6% of the TDR patients. Thanks to resistance testing, 83% of the patients with TDR who started therapy had undetectable viral load whereas half of the patients would likely have received a suboptimal therapy without this test. In conclusion, TDR remained stable and a NNRTI up-and-down trend was observed. While the presence of clusters with TDR is worrying, we could not identify an independent, non-sequence based predictor for TDR or transmission clusters with TDR that could help with guidelines or public health measures.

A new multicopper oxidase gene AaMco1 was identified in Acidomyces acidophilus, a pigmented extremophile ascomycete originally isolated from acidic water. Sequence analysis revealed that it encodes a 682 amino acid protein with an apparent molecular mass of 85 kDa as determined by denaturing SDS-PAGE. Interestingly, AaMco1 has a predicted N-terminal transmembrane helix and no signal peptide. To obtain an active and soluble protein, AaMco1 was truncated at its N-terminal to remove the transmembrane helix, but even in this form the protein was found in the insoluble fraction. AaMco1 and its truncated form were then denatured, purified and renatured before characterization. Structural analysis and protein characterization by enzymatic assays indicate that AaMco1 has ferroxidase activity. AaMco1 is also able to oxidize the DMPPDA compound and could be part of a new phylogenetic cluster, the ascomycete MCOs family, described for the first time here., (Copyright © 2014 Elsevier Masson SAS. All rights reserved.)

We investigated the origin and the effect of insertion D67D-THGERDLGPA within HIV-1 RT from a patient failing antiviral therapy. The insertion developed within the context of pre-existing NRTI and NNRTI mutations (M41L, L210W, T215Y and N348I). Concurrently, the NRTI mutations T69I and V118I and the NNRTI mutations K103N and Y181C were detected for the first time. High-level drug resistance (fold-changes≥50) and a good replication capacity (87% of wild-type) were observed, significantly higher than for the previous virus without insertion. The insertion was very similar to a region within human chromosome 17 (31/34 nucleotide identity), and had already been detected independently in a Japanese HIV-1 isolate. These results suggest that a particular sequence within human chromosome 17 is prone to horizontal gene transfer into the HIV-1 RT finger subdomain. This insertion confers selective advantage to HIV-1 by its contribution to multi-drug resistance and restoration of impaired replication capacity., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Background: In Europe, a continuous programme (SPREAD) has been in place for ten years to study transmission of drug resistant HIV. We analysed time trends of transmitted drug resistance mutations (TDRM) in relation to the risk behaviour reported., Methods: HIV-1 patients newly diagnosed in 27 countries from 2002 through 2007 were included. Inclusion was representative for risk group and geographical distribution in the participating countries in Europe. Trends over time were calculated by logistic regression., Results: From the 4317 patients included, the majority was men-having-sex-with-men -MSM (2084, 48%), followed by heterosexuals (1501, 35%) and injection drug users (IDU) (355, 8%). MSM were more often from Western Europe origin, infected with subtype B virus, and recently infected (<1 year) (p<0.001). The prevalence of TDRM was highest in MSM (prevalence of 11.1%), followed by heterosexuals (6.6%) and IDU (5.1%, p<0.001). TDRM was predominantly ascribed to nucleoside reverse transcriptase inhibitors (NRTI) with a prevalence of 6.6% in MSM, 3.3% in heterosexuals and 2.0% in IDU (p = 0.001). A significant increase in resistance to non- nucleoside reverse transcriptase inhibitors (NNRTIs) and a decrease in resistance to protease inhibitors was observed in MSM (p = 0.008 and p = 0.006, respectively), but not in heterosexual patients (p = 0.68 and p = 0.14, respectively)., Conclusions: MSM showed to have significantly higher TDRM prevalence compared to heterosexuals and IDU. The increasing NNRTI resistance in MSM is likely to negatively influence the therapy response of first-line therapy, as most include NNRTI drugs.

Transmission lies at the interface of human immunodeficiency virus type 1 (HIV-1) evolution within and among hosts and separates distinct selective pressures that impose differences in both the mode of diversification and the tempo of evolution. In the absence of comprehensive direct comparative analyses of the evolutionary processes at different biological scales, our understanding of how fast within-host HIV-1 evolutionary rates translate to lower rates at the between host level remains incomplete. Here, we address this by analyzing pol and env data from a large HIV-1 subtype C transmission chain for which both the timing and the direction is known for most transmission events. To this purpose, we develop a new transmission model in a Bayesian genealogical inference framework and demonstrate how to constrain the viral evolutionary history to be compatible with the transmission history while simultaneously inferring the within-host evolutionary and population dynamics. We show that accommodating a transmission bottleneck affords the best fit our data, but the sparse within-host HIV-1 sampling prevents accurate quantification of the concomitant loss in genetic diversity. We draw inference under the transmission model to estimate HIV-1 evolutionary rates among epidemiologically-related patients and demonstrate that they lie in between fast intra-host rates and lower rates among epidemiologically unrelated individuals infected with HIV subtype C. Using a new molecular clock approach, we quantify and find support for a lower evolutionary rate along branches that accommodate a transmission event or branches that represent the entire backbone of transmitted lineages in our transmission history. Finally, we recover the rate differences at the different biological scales for both synonymous and non-synonymous substitution rates, which is only compatible with the 'store and retrieve' hypothesis positing that viruses stored early in latently infected cells preferentially transmit or establish new infections upon reactivation.

Background: Human T-cell lymphotropic virus (HTLV-1) is the causative agent of the incapacitating, neuroinflammatory disease HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). Currently, there are no disease-modifying therapies with long-term clinical benefits or validated biomarkers for clinical follow-up in HAM/TSP. Although CD80 and CD86 costimulatory molecules play prominent roles in immune regulation and reflect disease status in multiple sclerosis (MS), data in HAM/TSP are lacking., Methods: Using flow cytometry, we quantified ex vivo and in vitro expression of CD80 and CD86 in PBMCs of healthy controls, HTLV-1-infected individuals with and without HAM/TSP, and MS patients. We hypothesized ex vivo CD80 and CD86 expressions and their in vitro regulation by interferon (IFN)-α/β mirror similarities between HAM/TSP and MS and hence might reveal clinically useful biomarkers in HAM/TSP., Results: Ex vivo expression of CD80 and CD86 in T and B cells increased in all HTLV-1 infected individuals, but with a selective defect for B cell CD86 upregulation in HAM/TSP. Despite decreased total B cells with increasing disease duration (p = 0.0003, r = -0.72), CD80+ B cells positively correlated with disease severity (p = 0.0017, r = 0.69) in HAM/TSP. B cell CD80 expression was higher in women with HAM/TSP, underscoring that immune markers can reflect the female predominance observed in most autoimmune diseases. In contrast to MS patients, CD80+ (p = 0.0001) and CD86+ (p = 0.0054) lymphocytes expanded upon in vitro culture in HAM/TSP patients. The expansion of CD80+ and CD86+ T cells but not B cells was associated with increased proliferation in HTLV-1 infection. In vitro treatment with IFN-β but not IFN-α resulted in a pronounced increase of B cell CD86 expression in healthy controls, as well as in patients with neuroinflammatory disease (HAM/TSP and MS), similar to in vivo treatment in MS., Conclusions: We propose two novel biomarkers, ex vivo CD80+ B cells positively correlating to disease severity and CD86+ B cells preferentially induced by IFN-β, which restores defective upregulation in HAM/TSP. This study suggests a role for B cells in HAM/TSP pathogenesis and opens avenues to B cell targeting (with proven clinical benefit in MS) in HAM/TSP but also CD80-directed immunotherapy, unprecedented in both HAM/TSP and MS.

Background: Superinfection with drug resistant HIV strains could potentially contribute to compromised therapy in patients initially infected with drug-sensitive virus and receiving antiretroviral therapy. To investigate the importance of this potential route to drug resistance, we developed a bioinformatics pipeline to detect superinfection from routinely collected genotyping data, and assessed whether superinfection contributed to increased drug resistance in a large European cohort of viremic, drug treated patients., Methods: We used sequence data from routine genotypic tests spanning the protease and partial reverse transcriptase regions in the Virolab and EuResist databases that collated data from five European countries. Superinfection was indicated when sequences of a patient failed to cluster together in phylogenetic trees constructed with selected sets of control sequences. A subset of the indicated cases was validated by re-sequencing pol and env regions from the original samples., Results: 4425 patients had at least two sequences in the database, with a total of 13816 distinct sequence entries (of which 86% belonged to subtype B). We identified 107 patients with phylogenetic evidence for superinfection. In 14 of these cases, we analyzed newly amplified sequences from the original samples for validation purposes: only 2 cases were verified as superinfections in the repeated analyses, the other 12 cases turned out to involve sample or sequence misidentification. Resistance to drugs used at the time of strain replacement did not change in these two patients. A third case could not be validated by re-sequencing, but was supported as superinfection by an intermediate sequence with high degenerate base pair count within the time frame of strain switching. Drug resistance increased in this single patient., Conclusions: Routine genotyping data are informative for the detection of HIV superinfection; however, most cases of non-monophyletic clustering in patient phylogenies arise from sample or sequence mix-up rather than from superinfection, which emphasizes the importance of validation. Non-transient superinfection was rare in our mainly treatment experienced cohort, and we found a single case of possible transmitted drug resistance by this route. We therefore conclude that in our large cohort, superinfection with drug resistant HIV did not compromise the efficiency of antiretroviral treatment.

Endo-inulinase INU2 from Aspergillus ficuum belongs to glycosidase hydrolase family 32 (GH32) that degrades inulin into fructo oligosaccharides consisting mainly of inulotriose and inulotetraose. The 3D structure of INU2 was recently obtained (Pouyez et al., 2012, Biochimie, 94, 2423-2430). An enlarged cavity compared to exo-inulinase formed by the conserved motif W-M(I)-N-D(E)-P-N-G, the so-called loop 1 and the loop 4, was identified. In the present study we have characterized the importance of 12 residues situated around the enlarged cavity. These residues were mutated by site-directed mutagenesis. Comparative activity analysis was done by plate, spectrophotometric and thin-layer chromatography assay. Most of the mutants were less active than the wild-type enzyme. Most interestingly, mutant N42G differed in the size distribution of the FOS synthesized.

Background: HIV-1 replication can be successfully blocked by targeting gag gene products, offering a promising strategy for new drug classes that complement current HIV-1 treatment options. However, naturally occurring polymorphisms at drug binding sites can severely compromise HIV-1 susceptibility to gag inhibitors in clinical and experimental studies. Therefore, a comprehensive understanding of gag natural diversity is needed., Findings: We analyzed the degree of functional conservation in 10862 full-length gag sequences across 8 major HIV-1 subtypes and identified the impact of natural variation on known drug binding positions targeted by more than 20 gag inhibitors published to date. Complete conservation across all subtypes was detected in 147 (29%) out of 500 gag positions, with the highest level of conservation observed in capsid protein. Almost half (41%) of the 136 known drug binding positions were completely conserved, but all inhibitors were confronted with naturally occurring polymorphisms in their binding sites, some of which correlated with HIV-1 subtype. Integration of sequence and structural information revealed one drug binding pocket with minimal genetic variability, which is situated at the N-terminal domain of the capsid protein., Conclusions: This first large-scale analysis of full-length HIV-1 gag provided a detailed mapping of natural diversity across major subtypes and highlighted the considerable variation in current drug binding sites. Our results contribute to the optimization of gag inhibitors in rational drug design, given that drug binding sites should ideally be conserved across all HIV-1 subtypes.

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.)

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.)

Genotypic drug resistance testing is routine practice in HIV-1 clinical care. The visual interpretation of sequencing electropherograms is labour-intensive and subject to intra- and inter-assay variability because decisions are based on operators' judgments. In this study the performance of the automatic editing tool RECall was compared to the current standard of editing manually and editing using the tool ViroSeq. Using RECall a consensus sequence could be generated for 97% of the V3 loop and for 79% of the pol experiments. By comparison, using manual editing a consensus sequence could be reached for 87% of the V3 and 87% of the pol experiments. Using ViroSeq, a consensus sequence was generated for 68% of the pol experiments. On a predefined dataset, manual editing displayed the highest probability to accurately assign mixtures (0.91 vs. 0.88 by ViroSeq vs. 0.76 by RECall) and the lowest probability to inaccurately assign a mixture to a pure base call (0.002 vs. 0.019 by ViroSeq vs. 0.002 by RECall). As differences in base calling have little impact on drug resistance interpretation and hands-on-time could be substantially reduced, RECall could be a valuable tool for the standardization and acceleration of the editing process., (Copyright © 2013 Elsevier B.V. All rights reserved.)

Objectives: This study evaluates the impact of specific HIV-1 protease-compensatory mutations (wild-type amino acids in non-B subtypes) on virological response to a first-line lopinavir/ritonavir-containing regimen in an HIV-1 subtype B-infected population., Patients and Methods: The prevalence of protease-compensatory mutations from 1997 to 2011 was calculated in 3063 drug-naive HIV-1 B-infected patients. The role of these mutations on virological outcome is estimated in a subgroup of 201 patients starting their first lopinavir/ritonavir-containing regimen by covariation and docking analyses., Results: The number of HIV-1 B-infected patients with at least one protease-compensatory mutation increased over time (from 86.4% prior to 2001 to 92.6% after 2009, P = 0.02). Analysing 201 patients starting first-line lopinavir/ritonavir, the median time to virological failure was shorter in patients with at least one protease-compensatory mutation than in patients with no protease-compensatory mutations. By covariation and docking analyses, specific mutations were found to affect lopinavir affinity for HIV-1 protease and to impact virological failure. Specifically, the L10V + I13V + L63P + I93L cluster, related to fast virological failure, correlated with a decreased drug affinity for the enzyme in comparison with wild-type (ΔGmut = -30.0 kcal/mol versus ΔGwt = -42.3 kcal/mol)., Conclusions: Our study shows an increased prevalence of specific protease-compensatory mutations in an HIV-1 B-infected population and confirms that their copresence can affect the virological outcome in patients starting a lopinavir/ritonavir-containing regimen.

A large and complex outbreak of hepatitis C virus in Valencia, Spain that began 25 years ago led to the prosecution and conviction of an anesthetist who was accused of infecting hundreds of his patients. Evolutionary analyses of viral gene sequences were presented as evidence in the trial, and these are now described in detail by González-Candelas and colleagues in a paper published in BMC Biology. Their study illustrates the challenges and opportunities that arise from the use of phylogenetic inference in criminal trials concerning virus transmission.

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.

In Cuba, antiretroviral therapy rollout started in 2001 and antiretroviral therapy coverage has reached almost 40% since then. The objectives of this study were therefore to analyze subtype distribution, and level and patterns of drug resistance in therapy-naive HIV-1 patients. Four hundred and one plasma samples were collected from HIV-1 therapy-naive patients in 2003 and in 2007-2011. HIV-1 drug resistance genotyping was performed in the pol gene and drug resistance was interpreted according to the WHO surveillance drug-resistance mutations list, version 2009. Potential impact on first-line therapy response was estimated using genotypic drug resistance interpretation systems HIVdb version 6.2.0 and Rega version 8.0.2. Phylogenetic analysis was performed using Neighbor-Joining. The majority of patients were male (84.5%), men who have sex with men (78.1%) and from Havana City (73.6%). Subtype B was the most prevalent subtype (39.3%), followed by CRF20-23-24&#95;BG (19.5%), CRF19&#95;cpx (18.0%) and CRF18&#95;cpx (10.3%). Overall, 29 patients (7.2%) had evidence of drug resistance, with 4.0% (CI 1.6%-4.8%) in 2003 versus 12.5% (CI 7.2%-14.5%) in 2007-2011. A significant increase in drug resistance was observed in recently HIV-1 diagnosed patients, i.e. 14.8% (CI 8.0%-17.0%) in 2007-2011 versus 3.8% (CI 0.9%-4.7%) in 2003 (OR 3.9, CI 1.5-17.0, p=0.02). The majority of drug resistance was restricted to a single drug class (75.8%), with 55.2% patients displaying nucleoside reverse transcriptase inhibitor (NRTI), 10.3% non-NRTI (NNRTI) and 10.3% protease inhibitor (PI) resistance mutations. Respectively, 20.7% and 3.4% patients carried viruses containing drug resistance mutations against NRTI+NNRTI and NRTI+NNRTI+PI. The first cases of resistance towards other drug classes than NRTI were only detected from 2008 onwards. The most frequent resistance mutations were T215Y/rev (44.8%), M41L (31.0%), M184V (17.2%) and K103N (13.8%). The median genotypic susceptibility score for the commonly prescribed first-line therapies was 2.5. This analysis emphasizes the need to perform additional surveillance studies to accurately assess the level of transmitted drug resistance in Cuba, as the extent of drug resistance might jeopardize effectiveness of first-line regimens prescribed in Cuba and might necessitate the implementation of baseline drug resistance testing., (Copyright © 2013 Elsevier B.V. All rights reserved.)

Introduction: Clinically evaluating genotypic interpretation systems is essential to provide optimal guidance in designing potent individualized HIV-regimens. This study aimed at investigating the ability of the latest Rega algorithm to predict virological response on a short and longer period., Materials Methods: 9231 treatment changes episodes were extracted from an integrated patient database. The virological response after 8, 24 and 48 weeks was dichotomized to success and failure. Success was defined as a viral load below 50 copies/ml or alternatively, a 2 log decrease from the baseline viral load at 8 weeks. The predictive ability of Rega version 8 was analysed in comparison with that of previous evaluated version Rega 5 and two other algorithms (ANRS v2011.05 and Stanford HIVdb v6.0.11). A logistic model based on the genotypic susceptibility score was used to predict virological response, and additionally, confounding factors were added to the model. Performance of the models was compared using the area under the ROC curve (AUC) and a Wilcoxon signed-rank test., Results: Per unit increase of the GSS reported by Rega 8, the odds on having a successful therapy response on week 8 increased significantly by 81% (OR = 1.81, CI = [1.76-1.86]), on week 24 by 73% (OR = 1.73, CI = [1.69-1.78]) and on week 48 by 85% (OR = 1.85, CI = [1.80-1.91]). No significant differences in AUC were found between the performance of Rega 8 and Rega 5, ANRS v2011.05 and Stanford HIVdb v6.0.11, however Rega 8 had the highest sensitivity: 76.9%, 76.5% and 77.2% on 8, 24 and 48 weeks respectively. Inclusion of additional factors increased the performance significantly., Conclusion: Rega 8 is a significant predictor for virological response with a better sensitivity than previously, and with rules for recently approved drugs. Additional variables should be taken into account to ensure an effective regimen.

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: We investigated gender differences in treatment outcome during first line antiretroviral treatment (ART) in a hospital setting in Tanzania, assessing clinical, social demographic, virological and immunological factors., Methods: We conducted a cohort study involving HIV infected patients scheduled to start ART and followed up to 1 year on ART. Structured questionnaires and patients file review were used to collect information and blood was collected for CD4 and viral load testing. Gender differences were assessed using Kruskal-Wallis test and chi-square test for continuous and categorical data respectively. Survival distributions for male and female patients were estimated using the Kaplan-Meier method and compared using Cox proportional hazards models., Results: Of 234 patients recruited in this study, 70% were females. At baseline, women had significantly lower education level; lower monthly income, lower knowledge on ARV, less advanced HIV disease (33% women; 47% men started ART at WHO stage IV, p = 0.04), higher CD4 cell count (median 149 for women, 102 for men, p = 0.02) and higher BMI (p = 0.002). After 1 year of standard ART, a higher proportion of females survived although this was not significant, a significantly higher proportion of females had undetectable plasma viral load (69% women, 45% men, p = 0.003), however females ended at a comparable CD4 cell count (median CD4, 312 women; 321 men) signifying a worse CD4 cell increase (p = 0.05), even though they still had a higher BMI (p = 0.02). The unadjusted relative hazard for death for men compared to women was 1.94. After correcting for confounding factors, the Cox proportional hazards showed no significant difference in the survival rate (relative hazard 1.02)., Conclusion: We observed women were starting treatment at a less advanced disease stage, but they had a lower socioeconomical status. After one year, both men and women had similar clinical and immunological conditions. It is not clear why women lose their immunological advantage over men despite a better virological treatment response. We recommend continuous follow up of this and more cohorts of patients to better understand the underlying causes for these differences and whether this will translate also in longer term differences.

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: Emergence of HIV-1 drug resistance may limit the sustained benefits of antiretroviral therapy (ART) in settings with limited laboratory monitoring and drug options., Objectives: Surveillance of drug resistance and subtypes in HIV-1 patients failing ART in Cuba., Study Design: This study compiled data of ART-experienced HIV-1 patients attending a clinical center in Havana in 2003 and 2009-2011. The first period included results of a cross-sectional study, whereas in the second period genotyping was performed as part of routine care. Drug resistance mutations and levels were determined using HIVdb version 6.0.9., Results: Seventy-six percent received solely ART containing at least 3 drugs, of which 79.1% ever receiving unboosted protease inhibitors (PI). Patients from 2009 to 2011 were longer treated and exposed to more ART regimens. Subtype B (39%) and CRF19&#95;cpx (18%) were the most prevalent genetic forms. Subtype distribution did not change significantly between both periods, except for BG recombinants that increased from 6% to 14%. Nucleoside reverse transcriptase inhibitor (NRTI), non-nucleoside RTI (NNRTI) and PI mutations were present in 69.5%, 54.8% and 44.4%. Full-class resistance (FCR) to NRTI, NNRTI, PI and multidrug resistance (MDR) were detected in 31.8%, 37.9%, 18.5% and 15.4%. FCR to NRTI, NNRTI, PI and MDR were present in 9.8%, 14.1%, 0%, 0% after first-line failure and in 19.8%, 20.8%, 2.9% and 2.9% after second-line failure., Conclusions: Our study found a high prevalence of drug resistance and supports the need for appropriate laboratory monitoring in clinical practice and access to drug options in case of virological failure., (Copyright © 2012 Elsevier B.V. All rights reserved.)