Colombe Chappey

Evidence for Positive Epistasis in HIV-1

Reproductive strategies such as sexual reproduction and recombination that involve the shuffling of parental genomes for the production of offspring are ubiquitous in nature. However, their evolutionary benefit remains unclear. Many theories have identified potential benefits, but progress is hampered by the scarcity of relevant data. One class of theories is based on the assumption that mutations affecting fitness exhibit negative epistasis. Retroviruses recombine frequently and thus provide a unique opportunity to test these theories. Using amino acid sequence data and fitness values from 9466 human immunodeficiency virus 1 (HIV-1) isolates, we find in contrast to these theories strong statistical evidence for a predominance of positive epistasis in HIV-1.

An Evolutionary Model-Based Algorithm for Accurate Phylogenetic Breakpoint Mapping and Subtype Prediction in HIV-1

Genetically diverse pathogens (such as Human Immunodeficiency virus type 1, HIV-1) are frequently stratified into phylogenetically or immunologically defined subtypes for classification purposes. Computational identification of such subtypes is helpful in surveillance, epidemiological analysis and detection of novel variants, e.g., circulating recombinant forms in HIV-1. A number of conceptually and technically different techniques have been proposed for determining the subtype of a query sequence, but there is not a universally optimal approach. We present a model-based phylogenetic method for automatically subtyping an HIV-1 (or other viral or bacterial) sequence, mapping the location of breakpoints and assigning parental sequences in recombinant strains as well as computing confidence levels for the inferred quantities. Our Subtype Classification Using Evolutionary ALgorithms (SCUEAL) procedure is shown to perform very well in a variety of simulation scenarios, runs in parallel when multiple sequences are being screened, and matches or exceeds the performance of existing approaches on typical empirical cases. We applied SCUEAL to all available polymerase (pol) sequences from two large databases, the Stanford Drug Resistance database and the UK HIV Drug Resistance Database. Comparing with subtypes which had previously been assigned revealed that a minor but substantial (≈5%) fraction of pure subtype sequences may in fact be within- or inter-subtype recombinants. A free implementation of SCUEAL is provided as a module for the HyPhy package and the Datamonkey web server. Our method is especially useful when an accurate automatic classification of an unknown strain is desired, and is positioned to complement and extend faster but less accurate methods. Given the increasingly frequent use of HIV subtype information in studies focusing on the effect of subtype on treatment, clinical outcome, pathogenicity and vaccine design, the importance of accurate, robust and extensible subtyping procedures is clear.

A Novel Substrate-Based HIV-1 Protease Inhibitor Drug Resistance Mechanism

Background HIV protease inhibitor (PI) therapy results in the rapid selection of drug resistant viral variants harbouring one or two substitutions in the viral protease. To combat PI resistance development, two approaches have been developed. The first is to increase the level of PI in the plasma of the patient, and the second is to develop novel PI with high potency against the known PI-resistant HIV protease variants. Both approaches share the requirement for a considerable increase in the number of protease mutations to lead to clinical resistance, thereby increasing the genetic barrier. We investigated whether HIV could yet again find a way to become less susceptible to these novel inhibitors. Methods and Findings We have performed in vitro selection experiments using a novel PI with an increased genetic barrier (RO033-4649) and demonstrated selection of three viruses 4- to 8-fold resistant to all PI compared to wild type. These PI-resistant viruses did not have a single substitution in the viral protease. Full genomic sequencing revealed the presence of NC/p1 cleavage site substitutions in the viral Gag polyprotein (K436E and/or I437T/V) in all three resistant viruses. These changes, when introduced in a reference strain, conferred PI resistance. The mechanism leading to PI resistance is enhancement of the processing efficiency of the altered substrate by wild-type protease. Analysis of genotypic and phenotypic resistance profiles of 28,000 clinical isolates demonstrated the presence of these NC/p1 cleavage site mutations in some clinical samples (codon 431 substitutions in 13%, codon 436 substitutions in 8%, and codon 437 substitutions in 10%). Moreover, these cleavage site substitutions were highly significantly associated with reduced susceptibility to PI in clinical isolates lacking primary protease mutations. Furthermore, we used data from a clinical trial (NARVAL, ANRS 088) to demonstrate that these NC/p1 cleavage site changes are associated with virological failure during PI therapy. Conclusions HIV can use an alternative mechanism to become resistant to PI by changing the substrate instead of the protease. Further studies are required to determine to what extent cleavage site mutations may explain virological failure during PI therapy.

A systems analysis of mutational effects in HIV-1 protease and reverse transcriptase

The development of a quantitative understanding of viral evolution and the fitness landscape in HIV-1 drug resistance is a formidable challenge given the large number of available drugs and drug resistance mutations. We analyzed a dataset measuring the in vitro fitness of 70,081 virus samples isolated from HIV-1 subtype B infected individuals undergoing routine drug resistance testing. We assayed virus samples for in vitro replicative capacity in the absence of drugs as well as in the presence of 15 individual drugs. We employed a generalized kernel ridge regression to estimate main fitness effects and epistatic interactions of 1,859 single amino acid variants found within the HIV-1 protease and reverse transcriptase sequences. Models including epistatic interactions predict an average of 54.8% of the variance in replicative capacity across the 16 different environments and substantially outperform models based on main fitness effects only. We find that the fitness landscape of HIV-1 protease and reverse transcriptase is characterized by strong epistasis.

Comparison of standard PCR/cloning to single genome sequencing for analysis of HIV-1 populations.

To compare standard PCR/cloning and single genome sequencing (SGS) in their ability to reflect actual intra-patient polymorphism of HIV-1 populations, a total of 530 HIV-1 pro-pol sequences obtained by both sequencing techniques from a set of 17 ART naïve patient specimens was analyzed. For each specimen, 12 and 15 sequences, on average, were characterized by the two techniques. Using phylogenetic analysis, tests for panmixia and entropy, and Bland-Altman plots, no difference in population structure or genetic diversity was shown in 14 of the 17 subjects. Evidence of sampling bias by the presence of subsets of identical sequences was found by either method. Overall, the study shows that neither method was more biased than the other, and providing that an adequate number of PCR templates is analyzed, and that the bulk sequencing captures the diversity of the viral population, either method is likely to provide a similar measure of population diversity.

Diversity of the V3 region of HIV in Paris, France.

Objective:To carry out, within France, a large-scale molecular epidemiological investigation on the principal neutralizing determinant of HIV-1, located in the third variable region (V3) of the envelope protein. Such investigations are of the utmost importance in the identification and monitoring of the distribution and spread of different viral strains internationally. Design:Using polymerase chain reaction (PCR), we examined the genetic variation of the V3 region sequences of 28 HIV-infected patients from Paris, France. Results:Comparison of the Parisian V3 loop sequences with other published data indicates that the range of diversity in France is included within that of a large group that contains sequences from North America, the rest of Europe, Japan, India and Africa. Variability appears to be lower in the V3 loop than in its flanking regions. Five out of the six putative N-linked glycosylation sites show preferential alterations to charged amino acids. We report two motifs at the tip of the loop that have not been described previously. Conclusions:The structural homogeneity and the wide geographic representation of the major V3 group suggests that a common strategy could be applied to a large proportion of isolates in the development of a broad-spectrum HIV vaccine.

Antiretroviral drug susceptibility among drug-naive adults with recent HIV infection in Rakai, Uganda

OBJECTIVEnTo analyze antiretroviral drug susceptibility in HIV from recently infected adults in Rakai, Uganda, prior to the availability of antiretroviral drug treatment.nnnMETHODSnSamples obtained at the time of HIV seroconversion (1998-2003) were analyzed using the GeneSeq HIV and PhenoSense HIV assays (Monogram Biosciences, Inc., South San Francisco, California, USA).nnnRESULTSnTest results were obtained for 104 samples (subtypes: 26A, 1C, 66D, 9A/D, 1C/D, 1 intersubtype recombinant). Mutations used for genotypic surveillance of transmitted antiretroviral drug resistance were identified in six samples: three had nucleoside reverse transcriptase inhibitor (NRTI) surveillance mutations (two had M41L, one had K219R), and three had protease inhibitor surveillance mutations (I47V, F53L, N88D); none had nonnucleoside reverse transcriptase inhibitor (NNRTI) surveillance mutations. Other resistance-associated mutations were identified in some samples. However, none of the samples had a sufficient number of mutations to predict reduced antiretroviral drug susceptibility. Ten (9.6%) of the samples had reduced phenotypic susceptibility to at least one drug (one had partial susceptibility to didanosine, one had nevirapine resistance, and eight had resistance or partial susceptibility to at least one protease inhibitor). Fifty-three (51%) of the samples had hypersusceptibility to at least one drug (seven had zidovudine hypersusceptibility, 28 had NNRTI hypersusceptibility, 34 had protease inhibitor hypersusceptibility). Delavirdine hypersusceptibility was more frequent in subtype A than D. In subtype D, efavirenz hypersusceptibility was associated with substitutions at codon 11 in HIV-reverse transcriptase.nnnCONCLUSIONnPhenotyping detected reduced antiretroviral drug susceptibility and hypersusceptibility in HIV from some antiretroviral-naive Ugandan adults that was not predicted by genotyping. Phenotyping may complement genotyping for analysis of antiretroviral drug susceptibility in populations with nonsubtype B HIV infection.

Phenotypic Susceptibility to Didanosine Is Associated with Antiviral Activity in Treatment-Experienced Patients with HIV-1 Infection

OBJECTIVEnWe investigated the relationship between human immunodeficiency virus (HIV) phenotypic susceptibility to didanosine and the antiviral activity of didanosine (ddI) in the JAGUAR study.nnnMETHODSnBaseline plasma HIV phenotypic susceptibility to ddI was assessed using a phenotype assay of patients randomized to receive ddI or placebo for 4 weeks in addition to their current regimen. Phenotypic susceptibility scores (PSSs) were then calculated for each sample. Associations between PSS and week 4 reductions in plasma HIV-1 RNA load or virologic response were assessed using linear regression and Jonckherres test and the Wilcoxon and Cochran-Armitage tests, respectively.nnnRESULTSnIn the ddI arm, a significant association between reduction in viral load and continuous PSS was observed (P<.0001). Using distinct categories, an increasing fold change (FC) in susceptibility to ddI was strongly associated with smaller reductions in plasma HIV-1 RNA load (P<.0001). The proportion of virologic responders was 83% (15/18) for patients with a ddI FC < or =1.3, 50% (33/66) for patients with an FC of 1.3-2.2, and 29% (4/14) for patients with an FC > or =2.2 (P=.0008). After we determined these findings, 3 ddI FC categories were defined using 1.3 and 2.2 as thresholds.nnnCONCLUSIONSnThe relationship between phenotypic susceptibility to ddI and reduction in plasma HIV-1 RNA load describes a continuum. The establishment of a lower clinical cutoff at 1.3 and an upper clinical cutoff at 2.2 are clinically relevant.

Differences in Reversion of Resistance Mutations to Wild-Type under Structured Treatment Interruption and Related Increase in Replication Capacity

Background The CPCRA 064 study examined the effect of structured treatment interruption (STI) of up to 4 months followed by salvage treatment in patients failing therapy with multi-drug resistant HIV. We examined the relationship between the reversion rate of major reverse transcriptase (RT) resistance-associated mutations and change in viral replication capacity (RC). The dataset included 90 patients with RC and genotypic data from virus samples collected at 0 (baseline), 2 and 4 months of STI. Principal Findings Rapid shift towards wild-type RC was observed during the first 2 months of STI. Median RC increased from 47.5% at baseline to 86.0% at 2 months and to 97.5% at 4 months. Between baseline and 2 months of STI, T215F had the fastest rate of reversion (41%) and the reversion of E44D and T69D was associated with the largest changes in RC. Among the most prevalent RT mutations, M184V had the fastest rate of reversion from baseline to 2 months (40%), and its reversion was associated with the largest increase in RC. Most rates of reversion increased between 2 months and 4 months, but the change in RC was more limited as it was already close to 100%. The highest frequency of concurrent reversion was found for L100I and K103N. Mutagenesis tree models showed that M184V, when present, was overall the first mutation to revert among all the RT mutations reported in the study. Conclusion Longitudinal analysis of combined phenotypic and genotypic data during STI showed a large amount of variability in prevalence and reversion rates to wild-type codons among the RT resistance-associated mutations. The rate of reversion of these mutations may depend on the extent of RC increase as well as the co-occurring reversion of other mutations belonging to the same mutational pathway.

Immune biomarkers associated with clinical benefit from atezolizumab (MPDL3280a; anti-PD-L1) in advanced urothelial bladder cancer (UBC)

Meeting abstractsnnAtezolizumab (anti-PD-L1) has demonstrated robust clinical activity in UBC [[1][1]]. Elevated PD-L1 expression on tumor-infiltrating immune cells (IC) is associated with increased clinical efficacy; however, the contribution of other immune biomarkers is unknown. In this study, we