Kristof Theys

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

BackgroundUnderstanding 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.ResultsWe investigated the subtype distribution of HIV-1 in Europe and Israel in a representative sample of patients diagnosed between 2002 and 2005 and related it to the demographic data available. 2793 PRO-RT sequences were subtyped either with the REGA Subtyping tool or by a manual procedure that included phylogenetic tree and recombination analysis. The most prevalent subtypes/CRFs in our dataset were subtype B (66.1%), followed by sub-subtype A1 (6.9%), subtype C (6.8%) and CRF02_AG (4.7%). Substantial differences in the proportion of new diagnoses with distinct subtypes were found between European countries: the lowest proportion of subtype B was found in Israel (27.9%) and Portugal (39.2%), while the highest was observed in Poland (96.2%) and Slovenia (93.6%). Other subtypes were significantly more diagnosed in immigrant populations. Subtype B was significantly more diagnosed in men than in women and in MSM > IDUs > heterosexuals. Furthermore, the subtype distribution according to continent of origin of the patients suggests they acquired their infection there or in Europe from compatriots.ConclusionsThe association of subtype with demographic parameters suggests highly compartmentalized epidemics, determined by social and behavioural characteristics of the patients.

Global epidemiology of drug resistance after failure of WHO recommended first-line regimens for adult HIV-1 infection: A multicentre retrospective cohort study

Summary Background Antiretroviral therapy (ART) is crucial for controlling HIV-1 infection through wide-scale treatment as prevention and pre-exposure prophylaxis (PrEP). Potent tenofovir disoproxil fumarate-containing regimens are increasingly used to treat and prevent HIV, although few data exist for frequency and risk factors of acquired drug resistance in regions hardest hit by the HIV pandemic. We aimed to do a global assessment of drug resistance after virological failure with first-line tenofovir-containing ART. Methods The TenoRes collaboration comprises adult HIV treatment cohorts and clinical trials of HIV drug resistance testing in Europe, Latin and North America, sub-Saharan Africa, and Asia. We extracted and harmonised data for patients undergoing genotypic resistance testing after virological failure with a first-line regimen containing tenofovir plus a cytosine analogue (lamivudine or emtricitabine) plus a non-nucleotide reverse-transcriptase inhibitor (NNRTI; efavirenz or nevirapine). We used an individual participant-level meta-analysis and multiple logistic regression to identify covariates associated with drug resistance. Our primary outcome was tenofovir resistance, defined as presence of K65R/N or K70E/G/Q mutations in the reverse transcriptase (RT) gene. Findings We included 1926 patients from 36 countries with treatment failure between 1998 and 2015. Prevalence of tenofovir resistance was highest in sub-Saharan Africa (370/654 [57%]). Pre-ART CD4 cell count was the covariate most strongly associated with the development of tenofovir resistance (odds ratio [OR] 1·50, 95% CI 1·27–1·77 for CD4 cell count <100 cells per μL). Use of lamivudine versus emtricitabine increased the risk of tenofovir resistance across regions (OR 1·48, 95% CI 1·20–1·82). Of 700 individuals with tenofovir resistance, 578 (83%) had cytosine analogue resistance (M184V/I mutation), 543 (78%) had major NNRTI resistance, and 457 (65%) had both. The mean plasma viral load at virological failure was similar in individuals with and without tenofovir resistance (145 700 copies per mL [SE 12 480] versus 133 900 copies per mL [SE 16 650; p=0·626]). Interpretation We recorded drug resistance in a high proportion of patients after virological failure on a tenofovir-containing first-line regimen across low-income and middle-income regions. Effective surveillance for transmission of drug resistance is crucial. Funding The Wellcome Trust.Summary Background Antiretroviral therapy (ART) is crucial for controlling HIV-1 infection through wide-scale treatment as prevention and pre-exposure prophylaxis (PrEP). Potent tenofovir disoproxil fumarate-containing regimens are increasingly used to treat and prevent HIV, although few data exist for frequency and risk factors of acquired drug resistance in regions hardest hit by the HIV pandemic. We aimed to do a global assessment of drug resistance after virological failure with first-line tenofovir-containing ART. Methods The TenoRes collaboration comprises adult HIV treatment cohorts and clinical trials of HIV drug resistance testing in Europe, Latin and North America, sub-Saharan Africa, and Asia. We extracted and harmonised data for patients undergoing genotypic resistance testing after virological failure with a first-line regimen containing tenofovir plus a cytosine analogue (lamivudine or emtricitabine) plus a non-nucleotide reverse-transcriptase inhibitor (NNRTI; efavirenz or nevirapine). We used an individual participant-level meta-analysis and multiple logistic regression to identify covariates associated with drug resistance. Our primary outcome was tenofovir resistance, defined as presence of K65R/N or K70E/G/Q mutations in the reverse transcriptase ( RT ) gene. Findings We included 1926 patients from 36 countries with treatment failure between 1998 and 2015. Prevalence of tenofovir resistance was highest in sub-Saharan Africa (370/654 [57%]). Pre-ART CD4 cell count was the covariate most strongly associated with the development of tenofovir resistance (odds ratio [OR] 1·50, 95% CI 1·27–1·77 for CD4 cell count Interpretation We recorded drug resistance in a high proportion of patients after virological failure on a tenofovir-containing first-line regimen across low-income and middle-income regions. Effective surveillance for transmission of drug resistance is crucial. Funding The Wellcome Trust.

Genetic Diversity and Selective Pressure in Hepatitis C Virus Genotypes 1–6: Significance for Direct-Acting Antiviral Treatment and Drug Resistance

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.

Functional conservation of HIV-1 Gag: implications for rational drug design.

BackgroundHIV-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.FindingsWe 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.ConclusionsThis 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.

HIV-1 Subtype Is an Independent Predictor of Reverse Transcriptase Mutation K65R in HIV-1 Patients Treated with Combination Antiretroviral Therapy Including Tenofovir

ABSTRACT 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_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.

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

BackgroundThe 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.ResultsGenotypic 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.ConclusionsThese 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.

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

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

HIV-1 genotypic drug resistance testing: digging deep, reaching wide?

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.

Resistance pathways of human immunodeficiency virus type 1 against the combination of zidovudine and lamivudine.

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.

HIV-1 gp120 N-linked glycosylation differs between plasma and leukocyte compartments

BackgroundN-linked glycosylation is a major mechanism for minimizing virus neutralizing antibody response and is present on the Human Immunodeficiency Virus (HIV) envelope glycoprotein. Although it is known that glycosylation changes can dramatically influence virus recognition by the host antibody, the actual contribution of compartmental differences in N-linked glycosylation patterns remains unclear.Methodology and Principal FindingsWe amplified the env gp120 C2-V5 region and analyzed 305 clones derived from plasma and other compartments from 15 HIV-1 patients. Bioinformatics and Bayesian network analyses were used to examine N-linked glycosylation differences between compartments. We found evidence for cellspecific single amino acid changes particular to monocytes, and significant variation was found in the total number of N-linked glycosylation sites between patients. Further, significant differences in the number of glycosylation sites were observed between plasma and cellular compartments. Bayesian network analyses showed an interdependency between N-linked glycosylation sites found in our study, which may have immense functional relevance.ConclusionOur analyses have identified single cell/compartment-specific amino acid changes and differences in N-linked glycosylation patterns between plasma and diverse blood leukocytes. Bayesian network analyses showed associations inferring alternative glycosylation pathways. We believe that these studies will provide crucial insights into the host immune response and its ability in controlling HIV replication in vivo. These findings could also have relevance in shielding and evasion of HIV-1 from neutralizing antibodies.