Competitive fitness of nevirapine-resistant human immunodeficiency virus type 1 mutants

U-CAN-seq: A Universal Competition Assay by Nanopore Sequencing

RNA viruses quickly evolve subtle genotypic changes that can have major impacts on viral fitness and host range, with potential consequences for human health. It is therefore important to understand the evolutionary fitness of novel viral variants relative to well-studied genotypes of epidemic viruses. Competition assays are an effective and rigorous system with which to assess the relative fitness of viral genotypes. However, it is challenging to quickly and cheaply distinguish and quantify fitness differences between very similar viral genotypes. Here, we describe a protocol for using reverse transcription PCR in combination with commercial nanopore sequencing services to perform competition assays on untagged RNA viruses. Our assay, called the Universal Competition Assay by Nanopore Sequencing (U-CAN-seq), is relatively cheap and highly sensitive. We used a well-studied N24A mutation in the chikungunya virus (CHIKV) nsp3 gene to confirm that we could detect a competitive disadvantage using U-CAN-seq. We also used this approach to show that mutations to the CHIKV 5' conserved sequence element that disrupt sequence but not structure did not affect the fitness of CHIKV. However, similar mutations to an adjacent CHIKV stem loop (SL3) did cause a fitness disadvantage compared to wild-type CHIKV, suggesting that structure-independent, primary sequence determinants in this loop play an important role in CHIKV biology. Our novel findings illustrate the utility of the U-CAN-seq competition assay.

An MHV-68 Mutator Phenotype Mutant Virus, Confirmed by CRISPR/Cas9-Mediated Gene Editing of the Viral DNA Polymerase Gene, Shows Reduced Viral Fitness

Drug resistance studies on human γ-herpesviruses are hampered by the absence of an in vitro system that allows efficient lytic viral replication. Therefore, we employed murine γ-herpesvirus-68 (MHV-68) that efficiently replicates in vitro as a model to study the antiviral resistance of γ-herpesviruses. In this study, we investigated the mechanism of resistance to nucleoside (ganciclovir (GCV)), nucleotide (cidofovir (CDV), HPMP-5azaC, HPMPO-DAPy) and pyrophosphate (foscarnet (PFA)) analogues and the impact of these drug resistance mutations on viral fitness. Viral fitness was determined by dual infection competition assays, where MHV-68 drug-resistant viral clones competed with the wild-type virus in the absence and presence of antivirals. Using next-generation sequencing, the composition of the viral populations was determined at the time of infection and after 5 days of growth. Antiviral drug resistance selection resulted in clones harboring mutations in the viral DNA polymerase (DP), denoted Y383S^GCV, Q827R^HPMP-5azaC, G302W^PFA, K442T^PFA, G302W+K442T^PFA, C297W^HPMPO-DAPy and C981Y^CDV. Without antiviral pressure, viral clones Q827R^HPMP-5azaC, G302W^PFA, K442T^PFA and G302W+K442T^PFA grew equal to the wild-type virus. However, in the presence of antivirals, these mutants had a growth advantage over the wild-type virus that was moderately to very strongly correlated with antiviral resistance. The Y383S^GCV mutant was more fit than the wild-type virus with and without antivirals, except in the presence of brivudin. The C297W and C981Y changes were associated with a mutator phenotype and had a severely impaired viral fitness in the absence and presence of antivirals. The mutator phenotype caused by C297W in MHV-68 DP was validated by using a CRISPR/Cas9 genome editing approach.

Viral fitness of MHV-68 viruses harboring drug resistance mutations in the protein kinase or thymidine kinase

Murine γ-herpesvirus-68 (MHV-68), genetically and biologically related to human γ-herpesviruses Epstein-Barr virus and Kaposi's sarcoma-associated herpesvirus, can be easily propagated in vitro allowing drug resistance studies. Previously, we described specific changes in MHV-68 protein kinase (PK) or thymidine kinase (TK) associated with resistance to various purine or pyrimidine nucleoside analogues, respectively. To investigate how specific TK and PK mutations affect viral replication capacity, we performed dual infection competition assays in which wild-type and drug-resistant virus compete in absence or presence of antivirals in Vero cells. The composition of the mixed viral population was analyzed using next-generation sequencing and relative fitness of seven MHV-68 PK or TK mutants was calculated based on the frequency of viral variants at the time of infection and after 5-days growth. A MHV-68 mutant losing the PK function due to a 2-nucleotide deletion was less fit than the wild-type virus in absence of antivirals, consistent with the essential role of viral PKs during lytic replication, but overgrew the wild-type virus under pressure of purine nucleosides. TK mutant viruses, with frameshift or missense mutations, grew equal to wild-type virus in absence of antivirals, in accordance with the viral TK function only being essential in non-replicating or in TK-deficient cells, but were more fit when treated with pyrimidine nucleosides. Moreover, TK missense mutant viruses also increased fitness under pressure of antivirals other than pyrimidine nucleosides, indicating that MHV-68 TK mutations might influence viral fitness by acting on cellular and/or viral functions that are unrelated to nucleoside activation.

Trend of HIV transmitted drug resistance before and after implementation of HAART regimen restriction in the treatment of HIV-1 infected patients in southern Taiwan

Background

The use of fixed combination antiretroviral therapy with a low genetic barrier for the treatment of patients infected with human immunodeficiency virus (HIV) may affect the local HIV transmitted drug resistance (TDR) pattern. The present study aimed to investigate changes in the prevalence of HIV TDR following the implementation of a fixed regimen of HIV treatment in Taiwan in 2012.

Methods

TDR was measured in antiretroviral treatment-naïve HIV-1-infected individuals who participated in voluntary counseling and testing between 2007 and 2015 in southern Taiwan. Antiretroviral resistance mutations were interpreted using the HIVdb program from the Stanford University HIV Drug Resistance Database.

Results

Sequences were obtained from 377 consecutive individuals between 2007 and 2015. The overall prevalence rates of TDR HIV among the study population from 2007 to 2011 and 2012–2015 were 10.6 and 7.9%, respectively. Among the detected mutations, K103 N and V179D + K103R were more frequently observed after 2012. Four HIV-infected patients with K103 N variants were detected after 2012, and 4 of the 5 patients with V179D + K103R variants were found after 2012. No significant differences were observed in the TDRs among nucleoside reverse transcriptase inhibitors (NRTIs), non-NRTIs (NNRTIs), protease inhibitors, multiple drug resistance, and any drug resistance between period 1 (2007–2011) and period 2 (2012–2015).

Conclusions

A fixed treatment regimen with zidovudine/lamivudine + efavirenz or nevirapine as first-line therapy for treatment-naïve patients infected with HIV did not significantly increase the TDR during the 4-year follow-up period. Due to the increase in NNRTI resistance associated with mutations after 2012, a longer follow-up period and larger sample size are needed in future studies.

DISSEQT-DIStribution-based modeling of SEQuence space Time dynamics

Rapidly evolving microbes are a challenge to model because of the volatile, complex, and dynamic nature of their populations. We developed the DISSEQT pipeline (DIStribution-based SEQuence space Time dynamics) for analyzing, visualizing, and predicting the evolution of heterogeneous biological populations in multidimensional genetic space, suited for population-based modeling of deep sequencing and high-throughput data. The pipeline is openly available on GitHub (https://github.com/rasmushenningsson/DISSEQT.jl, accessed 23 June 2019) and Synapse (https://www.synapse.org/#!Synapse: syn11425758, accessed 23 June 2019), covering the entire workflow from read alignment to visualization of results. Our pipeline is centered around robust dimension and model reduction algorithms for analysis of genotypic data with additional capabilities for including phenotypic features to explore dynamic genotype-phenotype maps. We illustrate its utility and capacity with examples from evolving RNA virus populations, which present one of the highest degrees of genetic heterogeneity within a given population found in nature. Using our pipeline, we empirically reconstruct the evolutionary trajectories of evolving populations in sequence space and genotype-phenotype fitness landscapes. We show that while sequence space is vastly multidimensional, the relevant genetic space of evolving microbial populations is of intrinsically low dimension. In addition, evolutionary trajectories of these populations can be faithfully monitored to identify the key minority genotypes contributing most to evolution. Finally, we show that empirical fitness landscapes, when reconstructed to include minority variants, can predict phenotype from genotype with high accuracy.

Impact of antiretroviral resistance and virological failure on HIV-1 informational entropy

Objectives

The present study investigated the relationship between genomic variability and resistance of HIV-1 sequences in protease (PR) and reverse transcriptase (RT) regions of the pol gene. In addition, we analysed the resistance among 651 individuals presenting antiretroviral virological failure, from 2009 to 2011, in the state of São Paulo, Brazil.

Methods

The genomic variability was quantified by using informational entropy methods and the relationship between resistance and replicative fitness, as inferred by the residual viral load and CD4+ T cell count.

Results

The number of antiretroviral schemes is related to the number of resistance mutations in the HIV-1 PR (α = 0.2511, P = 0.0003, R2 = 0.8672) and the RT (α = 0.7892, P = 0.0001, R2 = 0.9141). Increased informational entropy rate is related to lower levels of HIV-1 viral loads (α = -0.0121, P = 0.0471, R2 = 0.7923), lower levels of CD4+ T cell counts (α = -0.0120, P = 0.0335, R2 = 0.8221) and a higher number of antiretroviral resistance-related mutations.

Conclusions

Less organized HIV genomes as inferred by higher levels of informational entropy relate to less competent host immune systems, lower levels of HIV replication and HIV genetic evolution as a consequence of antiretroviral resistance.

Increasing prevalence of K65K and K66K in HIV-1 subtype B reverse transcriptase

OBJECTIVE

Synonymous substitutions K65K/K66K in HIV-1 reverse transcriptase alleviate fitness and fidelity defects in HIV-1 molecular clones harboring thymidine analogue mutations (TAMs); however, their potential for transmission and persistence is unknown. Here, we investigated the temporal appearance of K65K/K66K relative to TAMs in a HIV-1 cohort, their prevalence over time, and their impact on viral fitness in the context of patient-derived reverse transcriptase sequences.

METHODS

Retrospective analyses of the temporal appearance and longitudinal prevalence of synonymous substitutions and drug resistance mutations were performed using the British Columbia Centre for Excellence in HIV/AIDS Drug Treatment Program (DTP) database. Plasma-derived HIV-1 from the DTP was used to generate infectious molecular clones. Growth competition assays were performed to determine viral fitness.

RESULTS

The prevalence of K65K/K66K in drug-naïve individuals tripled from 11% in 1997 to 37% in 2014 (P < 0.0001, n = 5221), with K66K mainly accounting for the increase. These mutations emerged in drug-treated individuals without TAMs in 14% of the cohort and conferred a fitness advantage in the context of patient-derived multidrug-resistant (MDR) virus in the absence of drug.

CONCLUSION

The appearance of K65K/K66K in drug-treated individuals was largely independent of TAMs, suggesting alternative factors are likely associated with their emergence. The increasing K65K/K66K prevalence to over a third of treatment-naïve individuals in the mostly subtype B DTP cohort and their ability to confer a fitness advantage to multidrug-resistant virus might explain the transmission and persistence of virus harbouring K65K/K66K in untreated individuals, and highlights their role in adaptive HIV-1 evolution.

Baseline antiretroviral resistance mutations and treatment-emergent resistance in HIV-1 RNA-suppressed patients switching to EVG/COBI/FTC/TDF or continuing on their PI-, NNRTI-, or RAL-based regimen

BACKGROUND

Stably suppressed HIV-1-infected patients that switched to elvitegravir/cobicistat/emtricitabine/tenofovir disoproxil fumarate (EVG/COBI/FTC/TDF) from regimens containing FTC/TDF plus a ritonavir-boosted protease inhibitor (PI + RTV), nonnucleoside reverse transcriptase inhibitor (NNRTI), or raltegravir in phase 3 studies STRATEGY-PI, STRATEGY-NNRTI, and GS-US-236-0123 maintained high rates of virologic suppression through 48 weeks. In this article, resistance analyses for these studies are described.

METHODS

HIV-1 historical genotypes obtained before therapy initiation were analyzed for preexisting/transmitted resistance (-R) in protease and reverse transcriptase (RT) and subtype. Patients with resistance to FTC/TDF were excluded. Viral isolates with HIV-1 RNA ≥400 copies per milliliter at confirmed virologic failure, discontinuation, or week 48 were analyzed for protease, RT, and integrase genotype and phenotype.

RESULTS

Historical genotypes from 626/628 subjects that switched to EVG/COBI/FTC/TDF indicated 25% had ≥1 primary resistance mutation in protease and/or RT. NNRTI-R was identified in 15%, NRTI-R in 8.3%, and PI-R in 3.7% of subjects. Week 48 virologic success rates (HIV-1 RNA <50 copies per milliliter) were 94% for all patients treated with EVG/COBI/FTC/TDF, 94% with preexisting resistance, 93% with subtype B, and 96% with non-B subtypes. Altogether, 2 subjects qualified for postbaseline resistance analyses. Neither had emergent resistance, and both resuppressed to HIV-1 RNA <50 copies per milliliter with no change in therapy.

CONCLUSIONS

Switching antiretroviral regimens to EVG/COBI/FTC/TDF in HIV-1 RNA-suppressed FTC/TDF-sensitive patients resulted in maintained virologic suppression through 48 weeks. Similar virologic success rates were achieved irrespective of the presence of preexisting resistance mutations or subtype. The lack of emergent resistance through 48 weeks supports utility of EVG/COBI/FTC/TDF for treatment-experienced patients seeking regimen modification or simplification.

A novel mutation, D404N, in the connection subdomain of reverse transcriptase of HIV-1 CRF08_BC subtype confers cross-resistance to NNRTIs

OBJECTIVES

Growing evidence suggests that mutations in the connection domain of the HIV-1 reverse transcriptase (RT) can contribute to viral resistance to RT inhibitors. This work was designed to determine the effects of a novel mutation, D404N, in the connection subdomain of RT of HIV-1 CRF08_BC subtype on drug resistance, viral replication capacity (RC) and RT activity.

METHODS

Mutation D404N, alone or together with the other reported mutations, was introduced into an HIV-1 CRF08_BC subtype infectious clone by site-directed mutagenesis. Viral susceptibility to nine RT inhibitors, viral RC and the DNA polymerase activity of viral RT of the constructed virus mutants were investigated. A modelling study using the server SWISS-MODEL was conducted to explore the possible structure-related drug resistance mechanism of the mutation D404N.

RESULTS

Single mutations D404N and H221Y conferred low-level resistance to nevirapine, efavirenz, rilpivirine and zidovudine. Double mutations Y181C/D404N and Y181C/H221Y significantly reduced susceptibility to NNRTIs. The most pronounced resistance to NNRTIs was observed with the triple mutation Y181C/D404N/H221Y. Virus containing D404N as the only mutation displayed ∼50% RC compared with the WT virus. The modelling study suggested that the D404N mutation might abolish the hydrogen bonds between residues 404 and K30 in p51 or K431 in p66, leading to impaired RT subunit structure and enhanced drug resistance.

CONCLUSIONS

These results indicate that D404N is a novel NNRTI-associated mutation in the HIV-1 subtype CRF08_BC and provides information valuable for the monitoring of clinical RTI resistance.

Diminished transmission of drug resistant HIV-1 variants with reduced replication capacity in a human transmission model

Background

Different patterns of drug resistance are observed in treated and therapy naïve HIV-1 infected populations. Especially the NRTI-related M184I/V variants, which are among the most frequently encountered mutations in treated patients, are underrepresented in the antiretroviral naïve population. M184I/V mutations are known to have a profound effect on viral replication and tend to revert over time in the new host. However it is debated whether a diminished transmission efficacy of HIV variants with a reduced replication capacity can also contribute to the observed discrepancy in genotypic patterns.

As dendritic cells (DCs) play a pivotal role in HIV-1 transmission, we used a model containing primary human Langerhans cells (LCs) and DCs to compare the transmission efficacy M184 variants (HIV-M184V/I/T) to HIV wild type (HIV-WT). As control, we used HIV harboring the NNRTI mutation K103N (HIV-K103N) which has a minor effect on replication and is found at a similar prevalence in treated and untreated individuals.

Results

In comparison to HIV-WT, the HIV-M184 variants were less efficiently transmitted to CCR5⁺ Jurkat T cells by both LCs and DCs. The transmission rate of HIV-K103N was slightly reduced to HIV-WT in LCs and even higher than HIV-WT in DCs. Replication experiments in CCR5⁺ Jurkat T cells revealed no apparent differences in replication capacity between the mutant viruses and HIV-WT. However, viral replication in LCs and DCs was in concordance with the transmission results; replication by the HIV-M184 variants was lower than replication by HIV-WT, and the level of replication of HIV-K103N was intermediate for LCs and higher than HIV-WT for DCs.

Conclusions

Our data demonstrate that drug resistant M184-variants display a reduced replication capacity in LCs and DCs which directly impairs their transmission efficacy. As such, diminished transmission efficacy may contribute to the lower prevalence of drug resistant variants in therapy naive individuals.

Persistence of frequently transmitted drug-resistant HIV-1 variants can be explained by high viral replication capacity

Background

In approximately 10% of newly diagnosed individuals in Europe, HIV-1 variants harboring transmitted drug resistance mutations (TDRM) are detected. For some TDRM it has been shown that they revert to wild type while other mutations persist in the absence of therapy. To understand the mechanisms explaining persistence we investigated the in vivo evolution of frequently transmitted HIV-1 variants and their impact on in vitro replicative capacity.

Results

We selected 31 individuals infected with HIV-1 harboring frequently observed TDRM such as M41L or K103N in reverse transcriptase (RT) or M46L in protease. In all these samples, polymorphisms at non-TDRM positions were present at baseline (median protease: 5, RT: 6). Extensive analysis of viral evolution of protease and RT demonstrated that the majority of TDRM (51/55) persisted for at least a year and even up to eight years in the plasma. During follow-up only limited selection of additional polymorphisms was observed (median: 1).

To investigate the impact of frequently observed TDRM on the replication capacity, mutant viruses were constructed with the most frequently encountered TDRM as site-directed mutants in the genetic background of the lab strain HXB2. In addition, viruses containing patient-derived protease or RT harboring similar TDRM were made. The replicative capacity of all viral variants was determined by infecting peripheral blood mononuclear cells and subsequently monitoring virus replication. The majority of site-directed mutations (M46I/M46L in protease and M41L, M41L + T215Y and K103N in RT) decreased viral replicative capacity; only protease mutation L90M did not hamper viral replication. Interestingly, most patient-derived viruses had a higher in vitro replicative capacity than the corresponding site-directed mutant viruses.

Conclusions

We demonstrate limited in vivo evolution of protease and RT harbouring frequently observed TDRM in the plasma. This is in line with the high in vitro replication capacity of patient-derived viruses harbouring TDRM compared to site-directed mutant viruses harbouring TDRM. As site-directed mutant viruses have a lower replication capacity than the patient-derived viruses with similar mutational patterns, we propose that (baseline) polymorphisms function as compensatory mutations improving viral replication capacity.

Baseline clinical HIV genotypes are a valid measure of transmitted drug resistance within the treatment-naive population

OBJECTIVE

To examine whether baseline clinical genotypes are equivalent to diagnostic serum genotypes for surveillance of HIV transmitted drug resistance (TDR).

DESIGN

Current HIV TDR surveillance in Canada is conducted through genotyping remnant diagnostic sera from new HIV diagnoses. As part of routine care, baseline genotyping is now conducted on all newly diagnosed HIV infections, with TDR data being generated a second time on the same patients.

METHODS

Surveillance genotyping, on HIV diagnostic serum, was performed on newly diagnosed HIV cases from 2007 to 2010 in Alberta, Canada. All subjects with a baseline clinical genotype result on file, and no evidence of antiretroviral therapy, were studied further. The HIV sequences from diagnosis and from the first clinical genotype were compared according to elapsed time between testing and by evaluating timing of infection based on BED capture enzyme immunoassay (BED-CEIA, abbreviated as BED in this article).

RESULTS

Eighty-seven genotype pairs were available for analysis, most of which were subtype B. The time between genotypes ranged from 0 to 755 days, with a median of 36 days and an interquartile range of 155.25 days. Genetic distance between genotypes varied between 0 and 0.03389 substitutions per site and did not correlate with sampling times. There was a tendency for the genotypes of infections classified as recent by BED to be more similar to their clinical genotypes but this effect was lost when adjusted for elapsed time between tests. There was no difference in the identified drug resistance.

CONCLUSIONS

Baseline clinical genotypes from treatment-naive patients may be used for HIV TDR surveillance.

Highly-sensitive allele-specific PCR testing identifies a greater prevalence of transmitted HIV drug resistance in Japan

BACKGROUND

The transmission of drug-resistant HIV in newly identified infected populations has become an underlying epidemic which can be better assessed with sensitive resistance testing. Since minority drug resistant variants cannot be detected by bulk sequencing, methods with improved sensitivity are required. Thus, the goal of this study was to evaluate if transmitted drug resistance mutations at minority levels in Japanese patients could be identified using highly sensitive allele-specific PCR (AS-PCR).

MATERIALS AND METHODS

Samples were taken from newly diagnosed HIV/AIDS cases at the National Nagoya Hospital from January 2008 to December 2009. All samples were bulk sequenced for HIV protease and reverse transcriptase. To detect minority populations with drug resistance, we used AS-PCR with mutation-specific primers designed for seven reverse transcriptase inhibitor resistance mutations, M41L, K65R, K70R, K103N, Y181C, M184V, and T215F/Y, and for three protease inhibitor resistance mutations, M46I/L and L90M.

RESULTS

We studied 149 newly identified HIV cases. Bulk sequencing detected 8 cases with NRTI resistance mutations (one with A62V, one D67E, one T215D, one T215E, two with T215L and two T215S) and 15 with PI resistance mutations (one with N88D and 14 with M46I). Results obtained by AS-PCR and bulk sequencing demonstrated good concordance but the AS-PCR enabled the detection of seven additional drug-resistant cases (one M41L, two with K65R, two with K70R, and one M184V) in the RT region. Additionally, AS-PCR assays identified 15 additional cases with M46I, five with M46L and four cases with L90M in the protease region.

CONCLUSIONS

Using AS-PCR substantially increased the detection of transmitted drug resistance in this population from 15.4% to 26.8%, further supporting the benefit of sensitive testing among drug-naïve populations. Since the clinical impact of minority drug-resistant populations is not fully comprehended for all mutations, follow-up studies are needed to understand their significance for treatment.

Ultrasensitive allele-specific PCR reveals rare preexisting drug-resistant variants and a large replicating virus population in macaques infected with a simian immunodeficiency virus containing human immunodeficiency virus reverse transcriptase

It has been proposed that most drug-resistant mutants, resulting from a single-nucleotide change, exist at low frequency in human immunodeficiency virus type 1 (HIV-1) and simian immunodeficiency virus (SIV) populations in vivo prior to the initiation of antiretroviral therapy (ART). To test this hypothesis and to investigate the emergence of resistant mutants with drug selection, we developed a new ultrasensitive allele-specific PCR (UsASP) assay, which can detect drug resistance mutations at a frequency of ≥0.001% of the virus population. We applied this assay to plasma samples obtained from macaques infected with an SIV variant containing HIV-1 reverse transcriptase (RT) (RT-simian-human immunodeficiency [SHIV](mne)), before and after they were exposed to a short course of efavirenz (EFV) monotherapy. We detected RT inhibitor (RTI) resistance mutations K65R and M184I but not K103N in 2 of 2 RT-SHIV-infected macaques prior to EFV exposure. After three doses over 4 days of EFV monotherapy, 103N mutations (AAC and AAT) rapidly emerged and increased in the population to levels of ∼20%, indicating that they were present prior to EFV exposure. The rapid increase of 103N mutations from <0.001% to 20% of the viral population indicates that the replicating virus population size in RT-SHIV-infected macaques must be 10(6) or more infected cells per replication cycle.

Clinical, virological and biochemical evidence supporting the association of HIV-1 reverse transcriptase polymorphism R284K and thymidine analogue resistance mutations M41L, L210W and T215Y in patients failing tenofovir/emtricitabine therapy

Background

Thymidine analogue resistance mutations (TAMs) selected under treatment with nucleoside analogues generate two distinct genotypic profiles in the HIV-1 reverse transcriptase (RT): (i) TAM1: M41L, L210W and T215Y, and (ii) TAM2: D67N, K70R and K219E/Q, and sometimes T215F. Secondary mutations, including thumb subdomain polymorphisms (e.g. R284K) have been identified in association with TAMs. We have identified mutational clusters associated with virological failure during salvage therapy with tenofovir/emtricitabine-based regimens. In this context, we have studied the role of R284K as a secondary mutation associated with mutations of the TAM1 complex.

Results

The cross-sectional study carried out with >200 HIV-1 genotypes showed that virological failure to tenofovir/emtricitabine was strongly associated with the presence of M184V (P < 10^-10) and TAMs (P < 10^-3), while K65R was relatively uncommon in previously-treated patients failing antiretroviral therapy. Clusters of mutations were identified, and among them, the TAM1 complex showed the highest correlation coefficients. Covariation of TAM1 mutations and V118I, V179I, M184V and R284K was observed. Virological studies showed that the combination of R284K with TAM1 mutations confers a fitness advantage in the presence of zidovudine or tenofovir. Studies with recombinant HIV-1 RTs showed that when associated with TAM1 mutations, R284K had a minimal impact on zidovudine or tenofovir inhibition, and in their ability to excise the inhibitors from blocked DNA primers. However, the mutant RT M41L/L210W/T215Y/R284K showed an increased catalytic rate for nucleotide incorporation and a higher RNase H activity in comparison with WT and mutant M41L/L210W/T215Y RTs. These effects were consistent with its enhanced chain-terminated primer rescue on DNA/DNA template-primers, but not on RNA/DNA complexes, and can explain the higher fitness of HIV-1 having TAM1/R284K mutations.

Conclusions

Our study shows the association of R284K and TAM1 mutations in individuals failing therapy with tenofovir/emtricitabine, and unveils a novel mechanism by which secondary mutations are selected in the context of drug-resistance mutations.

A hybrid stochastic-deterministic computational model accurately describes spatial dynamics and virus diffusion in HIV-1 growth competition assay

We present a new hybrid stochastic-deterministic, spatially distributed computational model to simulate growth competition assays on a relatively immobile monolayer of peripheral blood mononuclear cells (PBMCs), commonly used for determining ex vivo fitness of human immunodeficiency virus type-1 (HIV-1). The novel features of our approach include incorporation of viral diffusion through a deterministic diffusion model while simulating cellular dynamics via a stochastic Markov chain model. The model accounts for multiple infections of target cells, CD4-downregulation, and the delay between the infection of a cell and the production of new virus particles. The minimum threshold level of infection induced by a virus inoculum is determined via a series of dilution experiments, and is used to determine the probability of infection of a susceptible cell as a function of local virus density. We illustrate how this model can be used for estimating the distribution of cells infected by either a single virus type or two competing viruses. Our model captures experimentally observed variation in the fitness difference between two virus strains, and suggests a way to minimize variation and dual infection in experiments.

HIV-1 drug resistance in antiretroviral-naive individuals with HIV-1-associated tuberculous meningitis initiating antiretroviral therapy in Vietnam

BACKGROUND

Access to antiretroviral therapy (ART) for HIV-infected individuals in Vietnam is rapidly expanding, but there are limited data on HIV drug resistance (HIVDR) to guide ART strategies.

METHODS

We retrospectively conducted HIVDR testing in 220 ART-naive individuals recruited to a randomized controlled trial of immediate versus deferred ART in individuals with HIV-associated tuberculous meningitis in Ho Chi Minh City (HCMC) from 2005-2008. HIVDR mutations were identified by population sequencing of the HIV pol gene and were defined based on 2009 WHO surveillance drug resistance mutations (SDRMs).

RESULTS

We successfully sequenced 219/220 plasma samples of subjects prior to ART; 218 were subtype CRF01_AE and 1 was subtype B. SDRMs were identified in 14/219 (6.4%) subjects; 8/14 were resistant to nucleoside/nucleotide reverse transcriptase inhibitors (NRTIs; T69D, L74V, V75M, M184V/I and K219R), 5/14 to non-nucleoside reverse transcriptase inhibitors (NNRTIs; K103N, V106M, Y181C, Y188C and G190A), 1/14 to both NRTIs and NNRTIs (D67N and Y181C) and none to protease inhibitors. After 6 months of ART, eight subjects developed protocol-defined virological failure. HIVDR mutations were identified in 5/8 subjects. All five had mutations with high-level resistance to NNRTIs and three had mutations with high-level resistance to NRTIs. Due to a high early mortality rate (58%), the effect of pre-existing HIVDR mutations on treatment outcome could not be accurately assessed.

CONCLUSIONS

The prevalence of WHO SDRMs in ART-naive individuals with HIV-associated tuberculous meningitis in HCMC from 2005-2008 is 6.4%. The SDRMs identified conferred resistance to NRTIs and/or NNRTIs, reflecting the standard first-line ART regimens in Vietnam.

Exploring the complexity of the HIV-1 fitness landscape

Although fitness landscapes are central to evolutionary theory, so far no biologically realistic examples for large-scale fitness landscapes have been described. Most currently available biological examples are restricted to very few loci or alleles and therefore do not capture the high dimensionality characteristic of real fitness landscapes. Here we analyze large-scale fitness landscapes that are based on predictive models for in vitro replicative fitness of HIV-1. We find that these landscapes are characterized by large correlation lengths, considerable neutrality, and high ruggedness and that these properties depend only weakly on whether fitness is measured in the absence or presence of different antiretrovirals. Accordingly, adaptive processes on these landscapes depend sensitively on the initial conditions. While the relative extent to which mutations affect fitness on their own (main effects) or in combination with other mutations (epistasis) is a strong determinant of these properties, the fitness landscape of HIV-1 is considerably less rugged, less neutral, and more correlated than expected from the distribution of main effects and epistatic interactions alone. Overall this study confirms theoretical conjectures about the complexity of biological fitness landscapes and the importance of the high dimensionality of the genetic space in which adaptation takes place.

Evolutionary adaptation can be understood as populations moving uphill on landscapes, in which height corresponds to evolutionary fitness. Although such fitness landscapes are central to evolutionary theory, there is currently a lack of biologically realistic examples. Here we analyze large-scale fitness landscapes derived from in vitro fitness measurements of HIV-1. We find that these landscapes are very rugged and that, accordingly, adaptive processes on these landscapes depend sensitively on the initial conditions. Moreover, the landscapes contain large networks along which fitness changes only minimally. While the relative extent to which mutations affect fitness on their own or in combination with other mutations is a strong determinant of these properties, the fitness landscape of HIV-1 is considerably less rugged than expected from the individual and pair-wise effects of mutations. Overall this study confirms theoretical conjectures about the complexity of biological fitness landscapes and the importance of the high dimensionality of the genetic space in which adaptation takes place.

Transmitted drug resistance in the CFAR network of integrated clinical systems cohort: prevalence and effects on pre-therapy CD4 and viral load

Human immunodeficiency virus type 1 (HIV-1) genomes often carry one or more mutations associated with drug resistance upon transmission into a therapy-naïve individual. We assessed the prevalence and clinical significance of transmitted drug resistance (TDR) in chronically-infected therapy-naïve patients enrolled in a multi-center cohort in North America. Pre-therapy clinical significance was quantified by plasma viral load (pVL) and CD4+ cell count (CD4) at baseline. Naïve bulk sequences of HIV-1 protease and reverse transcriptase (RT) were screened for resistance mutations as defined by the World Health Organization surveillance list. The overall prevalence of TDR was 14.2%. We used a Bayesian network to identify co-transmission of TDR mutations in clusters associated with specific drugs or drug classes. Aggregate effects of mutations by drug class were estimated by fitting linear models of pVL and CD4 on weighted sums over TDR mutations according to the Stanford HIV Database algorithm. Transmitted resistance to both classes of reverse transcriptase inhibitors was significantly associated with lower CD4, but had opposing effects on pVL. In contrast, position-specific analyses of TDR mutations revealed substantial effects on CD4 and pVL at several residue positions that were being masked in the aggregate analyses, and significant interaction effects as well. Residue positions in RT with predominant effects on CD4 or pVL (D67 and M184) were re-evaluated in causal models using an inverse probability-weighting scheme to address the problem of confounding by other mutations and demographic or risk factors. We found that causal effect estimates of mutations M184V/I (-1.7 log₁₀pVL) and D67N/G (-2.1[³√CD4] and 0.4 log₁₀pVL) were compensated by K103N/S and K219Q/E/N/R. As TDR becomes an increasing dilemma in this modern era of highly-active antiretroviral therapy, these results have immediate significance for the clinical management of HIV-1 infections and our understanding of the ongoing adaptation of HIV-1 to human populations.

Evolutionary pathways of transmitted drug-resistant HIV-1

Several large studies in Europe and the USA revealed that approximately 10% of all newly diagnosed patients harbour HIV-1 variants with at least one major resistance-associated mutation. In this review we discuss the underlying mechanisms that drive the evolution of drug-resistant viruses after transmission to the new host. In a comprehensive literature search 12 papers describing the evolution of 58 cases of transmitted resistant HIV-1 variants were found. Based on observations in the literature we propose three pathways describing the evolution of resistant HIV-1 after transmission to a new host. Firstly, reversion of the resistance mutation towards wild-type may rapidly occur when drug resistance mutations severely impact replicative capacity. Alternatively, a second pathway involves replacement of transmitted drug resistance mutations by atypical amino acids that also improve viral replication capacity. In the third evolutionary pathway the resistance mutations persist either because they do not significantly affect viral replication capacity or evolution is constrained by fixation through compensatory mutations. In the near future ultra-sensitive resistance tests may provide more insight into the presence of archived and minority variants and their clinical relevance. Meanwhile, clinical guidelines advise population sequence analysis of the baseline plasma sample to identify transmission of resistance. Given the limited sensitivity of this technique for minority populations and the delay between the moment of infection and time of analysis, knowledge of the described evolutionary mechanisms of transmitted drug resistance patterns is essential for clinical management and public health strategies.

From resistance to stimulation: the evolution of a virus in the presence of a dominant lethal inhibitory scaffolding protein

By acquiring resistance to an inhibitor, viruses can become dependent on that inhibitor for optimal fitness. However, inhibitors rarely, if ever, stimulate resistant strain fitness to values that equal or exceed the uninhibited wild-type level. This would require an adaptive mechanism that converts the inhibitor into a beneficial replication factor. Using a plasmid-encoded inhibitory external scaffolding protein that blocks ϕX174 assembly, we previously demonstrated that such mechanisms are possible. The resistant strain, referred to as the evolved strain, contains four mutations contributing to the resistance phenotype. Three mutations confer substitutions in the coat protein, whereas the fourth mutation alters the virus-encoded external scaffolding protein. To determine whether stimulation by the inhibitory protein coevolved with resistance or whether it was acquired after resistance was firmly established, the strain temporally preceding the previously characterized mutant, referred to as the intermediary strain, was isolated and characterized. The results of the analysis indicated that the mutation in the virus-encoded external scaffolding protein was primarily responsible for stimulating strain fitness. When the mutation was placed in a wild-type background, it did not confer resistance. The mutation was also placed in cis with the plasmid-encoded dominant lethal mutation. In this configuration, the stimulating mutation exhibited no activity, regardless of the genotype (wild type, evolved, or intermediary) of the infecting virus. Thus, along with the coat protein mutations, stimulation required two external scaffolding protein genes: the once inhibitory gene and the mutant gene acquired during evolution.

Human immunodeficiency virus type 1 resistance or cross-resistance to nonnucleoside reverse transcriptase inhibitors currently under development as microbicides

Microbicides based on nonnucleoside reverse transcriptase inhibitors (NNRTIs) are currently being developed to protect women from HIV acquisition through sexual contact. However, the large-scale introduction of these products raises two major concerns. First, when these microbicides are used by undiagnosed HIV-positive women, they could potentially select for viral resistance, which may compromise subsequent therapeutic options. Second, NNRTI-based microbicides that are inactive against NNRTI-resistant strains might promote the selective transmission of these viruses. In order to address these concerns, drug resistance was selected in vitro by the serial passage of three viral isolates from subtypes B and C and CRF02_AG (a circulating recombinant form) in activated peripheral blood mononuclear cells (PBMCs) under conditions of increasing concentrations of three NNRTIs (i.e., TMC120, UC781, and MIV-160) that are currently being developed as candidate microbicides. TMC120 and MIV-160 displayed a high genetic barrier to resistance development, whereas resistance to UC781 emerged rapidly, similarly to efavirenz and nevirapine. Phenotypically, the selected viruses appeared to be highly cross-resistant to current first-line therapeutic NNRTIs (i.e., delavirdine, nevirapine, and efavirenz), although they retained some susceptibility to the more recently developed NNRTIs lersivirine and etravirine. The ability of UC781, TMC120, and MIV-160 to inhibit the in vitro-selected NNRTI-resistant viruses was also limited, although residual activity could be observed for the candidate microbicide NNRTI MIV-170. Interestingly, only four p2/p7/p1/p6/PR/RT/INT recombinant NNRTI-resistant viruses (i.e., TMC120-resistant VI829, EFV-resistant VI829, MIV-160-resistant VI829, and EFV-resistant MP568) showed impairments in replicative fitness. Overall, these in vitro analyses demonstrate that due to potential cross-resistance, the large-scale introduction of single-NNRTI-based microbicides should be considered with caution.

Replicative fitness costs of nonnucleoside reverse transcriptase inhibitor drug resistance mutations on HIV subtype C

Single-dose nevirapine (NVP) is quite effective in preventing transmission of the human immunodeficiency virus (HIV) from mother to child; however, many women develop resistance to NVP in this setting. Comparing outcomes of clinical studies reveals an increased amount of resistance in subtype C relative to that in other subtypes. This study investigates how nonnucleoside reverse transcriptase inhibitor (NNRTI) drug resistance mutations of subtype C affect replication capacity. The 103N, 106A, 106M, 181C, 188C, 188L, and 190A drug resistance mutations were placed in a reverse transcriptase (RT) that matches the consensus subtype C sequence as well as the HXB2 RT, as a subtype B reference. The replicative fitness of each mutant was compared with that of the wild type in a head-to-head competition assay. The 106A mutant of subtype C would not grow in the competition assay, making it the weakest virus tested. The effect of the 106M mutation was weaker than those of the 181C and 188C mutations in the consensus C RT, but in subtype B, this difference was not seen. To see if the 106A mutation in a different subtype C background would have a different replicative profile, the same NNRTI resistance mutations were added to the MJ4 RT, a reference subtype C molecular clone. In the context of MJ4 RT, the 106A mutant was not the only mutant that showed poor replicative fitness; the 106M, 188C, and 190A mutants also failed to replicate. These results suggest that NNRTIs may be a cost-effective alternative for salvage therapy if deleterious mutations are present in a subtype C setting.

Molecular and structural aspects of clinically relevant mutations related to the approved non-nucleoside inhibitors of HIV-1 reverse transcriptase

In recent years relevant progress has been made in the treatment of HIV-1 with a consequent decrease in mortality. The availability of potent antiretroviral drugs and the ability of viral load assays that accurately evaluate the true level of viral replication, have led to a better understanding of pathogenesis of the disease and how to obtain improved therapeutic profiles. The highly active antiretroviral therapy (HAART), based on a combination of three or more antiretroviral drugs, has radically changed the clinical outcome of HIV. In particular, reverse transcriptase non-nucleoside inhibitors (NNRTIs) play an essential role in most protocols and are often used in first line treatment. The high specificity of these inhibitors towards HIV-1 has increased the number of structural and molecular modeling studies of enzyme complexes and that have led to chemical syntheses of more selective second and third-generation NNRTIs. However, a considerable percentage of new HIV-1 infections are caused by the emergence of drug-resistant mutant viruses that complicate treatment strategies. In this review we discuss relevant clinical and structural aspects for the management of antiretroviral drug resistance, with detailed explanations of mechanisms and mutation patterns useful to better understand the relation between drug resistance and therapy failure.

Characterization of the E138K resistance mutation in HIV-1 reverse transcriptase conferring susceptibility to etravirine in B and non-B HIV-1 subtypes

We have selected for resistance to etravirine (ETR) and efavirenz (EFV) in tissue culture using three subtype B, three subtype C, and two CRF02_AG clinical isolates, grown in cord blood mononuclear cells. Genotypic analysis was performed at baseline and at various weeks of selection. Phenotypic resistance in regard to ETR, EFV, and nevirapine (NVP) was evaluated at weeks 25 to 30 for all ETR-selected viruses and in viral clones that contained specific resistance mutations that were inserted by site-directed mutagenesis into pNL-4.3 and AG plasmids. The results show that ETR selected mutations at positions V90I, K101Q, E138K, V179D/E/F, Y181C, V189I, G190E, H221H/Y, and M230L and that E138K was the first of these to emerge in most instances. The time to the emergence of resistance was longer in the case of ETR (18 weeks) compared to EFV (11 weeks), and no differences in the patterns of emergent mutations could be documented between the B and non-B subtypes. Viral clones containing E138K displayed low-level phenotypic resistance to ETR (3.8-fold) and modestly impaired replication capacity (2-fold) compared to wild-type virus. ETR-selected virus showed a high degree of cross-resistance to NVP but not to EFV. We identified K101Q, E138K, V179E, V189I, G190E, and H221Y as mutations not included among the 17 currently recognized resistance-associated mutations for ETR.

Determinants of individual variation in intracellular accumulation of anti-HIV nucleoside analog metabolites

Individual variation in response to antiretroviral therapy is well-known, but it is not clear if demographic characteristics such as gender, age, and ethnicity are responsible for the variation. To optimize anti-HIV therapy and guide antiretroviral drug discovery, determinants that cause variable responses to therapy need to be evaluated. We investigated the determinants of intracellular concentrations of nucleoside analogs using peripheral blood mononuclear cells from 40 healthy donors. We observed individual differences in the concentrations of the intracellular nucleoside analogs; the mean concentrations of the triphosphate metabolite of ethynylstavudine (4'-Ed4T), zidovudine (AZT), and lamivudine (3TC) were 0.71 pmol/10(6) cells (minimum and maximum, 0.10 and 3.00 pmol/10(6) cells, respectively), 0.88 pmol/10(6) cells (minimum and maximum, 0.10 and 15.18 pmol/10(6) cells, respectively), and 1.70 pmol/10(6) cells (minimum and maximum, 0.20 and 7.73 pmol/10(6) cells, respectively). Gender and ethnicity had no effect on the concentration of 4'-Ed4T and 3TC metabolites. There was a trend for moderation of the concentrations of AZT metabolites by gender (P = 0.17 for gender·metabolite concentration). We observed variability in the activity and expression of cellular kinases. There was no statistically significant correlation between thymidine kinase 1 (TK-1) activity or expression and thymidine analog metabolite concentrations. The correlation between the activity of deoxycytidine kinase (dCK) and the 3TC monophosphate metabolite concentration showed a trend toward significance (P = 0.1). We observed an inverse correlation between the multidrug-resistant protein 2 (MRP2) expression index and the concentrations of AZT monophosphate, AZT triphosphate, and total AZT metabolites. Our findings suggest that the observed variation in clinical response to nucleoside analogs may be due partly to the individual differences in the intracellular concentrations, which in turn may be affected by the cellular kinases involved in the phosphorylation pathway and ATP-binding cassette (ABC) transport proteins.

HIV-1 diversity after a class switch failure

The purpose of this study is to evaluate whether the choice of a PI- or an efavirenz (EFV)-based HAART initial regimen impacts on the viral diversity after failure from a second, class-switch salvage regimen. Sequential HAART failures after a class switch were identified for which the genotypes showed evidence of signature mutations at each failure. Each second failure was required to be from a viral burden <400 RNA c/ml. Thirteen cases of sequential failure from an initial EFV-containing to a PI-containing regimen (EP), and 19 sequential failures from an initial PI-containing to an EFV-containing regimen (PE) were identified. The persistence of signature mutations from the first failure were evaluated at second failure and compared between the EP and PE groups. Phylogenetic trees were constructed for a subgroup of cases from existing genetic sequence information and branch length analysis was used to determine evidence of viral diversity between groups. For EP sequential therapy, 10 of 12 cases carried forward a key non-nucleoside reverse transcriptase inhibitor (NNRTI) mutation in the second failure compared to 5 of 13 cases for PE sequential therapy (p = 0.041). Phylogenetic analysis demonstrated that there was more viral diversity in the PE group as compared to the EP group, consistent with the interpretation that mutations at the second failure added to an ancestral virus closer to baseline rather than to the dominant virus at first failure. The development of HIV viral diversity after multiple HAART failures is determined by the sequence in which the regimens are ordered.

Human immunodeficiency virus type 1 protease inhibitor drug-resistant mutants give discordant results when compared in single-cycle and multiple-cycle fitness assays

The replication fitness of HIV-1 drug-resistant mutants has been measured using either multiple-cycle or single-cycle assays (MCAs or SCAs); these assays have not been systematically compared. We developed an MCA and an SCA that utilized either intact or env-deleted recombinant viral vectors, respectively, in which virus-infected cells were detected by flow cytometry of a reporter gene product. Fitness was measured using each assay for 11 protease mutants, 9 reverse transcriptase mutants, and two mutants with mutations in gag p6, which is important for the release of virus particles from the cell membrane. In the SCA, fitness (replication capacity [RC]) was defined as the proportion of cells infected by the mutant compared to the wild type 40 h after infection. MCA fitness (1+s) was determined by comparing the changes in the relative proportions of cells infected by the mutant and the wild type between 3 and 5 days after infection. Five protease mutants showed statistically different fitness values by the MCA versus the SCA: the D30N, G48V, I50V, I54L, and I54M mutants. When all the mutants were ranked in order from most to least fit for both assays, 4 protease mutants moved more than 5 positions in rank: the D30N, I54L, I54M, and V82A mutants. There were no significant differences in fitness for the gag p6 or reverse transcriptase mutants. We propose that discordant results in the MCA and SCA are due to alterations in late events in the virus life cycle that are not captured in an SCA, such as burst size, cell-to-cell transmission, or infected-cell life span.

Computational models of HIV-1 resistance to gene therapy elucidate therapy design principles

Gene therapy is an emerging alternative to conventional anti-HIV-1 drugs, and can potentially control the virus while alleviating major limitations of current approaches. Yet, HIV-1's ability to rapidly acquire mutations and escape therapy presents a critical challenge to any novel treatment paradigm. Viral escape is thus a key consideration in the design of any gene-based technique. We develop a computational model of HIV's evolutionary dynamics in vivo in the presence of a genetic therapy to explore the impact of therapy parameters and strategies on the development of resistance. Our model is generic and captures the properties of a broad class of gene-based agents that inhibit early stages of the viral life cycle. We highlight the differences in viral resistance dynamics between gene and standard antiretroviral therapies, and identify key factors that impact long-term viral suppression. In particular, we underscore the importance of mutationally-induced viral fitness losses in cells that are not genetically modified, as these can severely constrain the replication of resistant virus. We also propose and investigate a novel treatment strategy that leverages upon gene therapy's unique capacity to deliver different genes to distinct cell populations, and we find that such a strategy can dramatically improve efficacy when used judiciously within a certain parametric regime. Finally, we revisit a previously-suggested idea of improving clinical outcomes by boosting the proliferation of the genetically-modified cells, but we find that such an approach has mixed effects on resistance dynamics. Our results provide insights into the short- and long-term effects of gene therapy and the role of its key properties in the evolution of resistance, which can serve as guidelines for the choice and optimization of effective therapeutic agents.

A primary obstacle to the success of any anti-HIV treatment is HIV's ability to rapidly resist it by generating new viral strains whose vulnerability to the treatment is reduced. Gene therapies represent a novel class of treatments for HIV infection that may supplement or replace present therapies, as they alleviate some of their major shortcomings. The design of gene therapeutic agents that effectively reduce viral resistance can be aided by a quantitative elucidation of the processes by which resistance is acquired following therapy initiation. We developed a computational model that describes a patient's response to therapy and used it to quantify the influence of therapy parameters and strategies on the development of viral resistance. We find that gene therapy induces different clinical conditions and a much slower viral response than present therapies. These dictate different design principles such as a greater significance to the virus' competence in the absence of therapy. We also show that one can effectively delay emergence of resistance by delivering distinct therapeutic genes into separate cell populations. Our results highlight the differences between traditional and gene therapies and provide a basic understanding of how key controllable parameters and strategies affect resistance development.

Estimating frequencies of minority nevirapine-resistant strains in chronically HIV-1-infected individuals naive to nevirapine by using stochastic simulations and a mathematical model

Nevirapine forms the mainstay of our efforts to curtail the pediatric AIDS epidemic through prevention of mother-to-child transmission of HIV-1. A key limitation, however, is the rapid selection of HIV-1 strains resistant to nevirapine following the administration of a single dose. This rapid selection of resistance suggests that nevirapine-resistant strains preexist in HIV-1 patients and may adversely affect outcomes of treatment. The frequencies of nevirapine-resistant strains in vivo, however, remain poorly estimated, possibly because they exist as a minority below current assay detection limits. Here, we employ stochastic simulations and a mathematical model to estimate the frequencies of strains carrying different combinations of the common nevirapine resistance mutations K103N, V106A, Y181C, Y188C, and G190A in chronically infected HIV-1 patients naïve to nevirapine. We estimate the relative fitness of mutant strains from an independent analysis of previous competitive growth assays. We predict that single mutants are likely to preexist in patients at frequencies ( approximately 0.01% to 0.001%) near or below current assay detection limits (>0.01%), emphasizing the need for more-sensitive assays. The existence of double mutants is subject to large stochastic variations. Triple and higher mutants are predicted not to exist. Our estimates are robust to variations in the recombination rate, cellular superinfection frequency, and the effective population size. Thus, with 10(7) to 10(8) infected cells in HIV-1 patients, even when undetected, nevirapine-resistant genomes may exist in substantial numbers and compromise efforts to prevent mother-to-child transmission of HIV-1, accelerate the failure of subsequent antiretroviral treatments, and facilitate the transmission of drug resistance.

Divergent evolution in reverse transcriptase (RT) of HIV-1 group O and M lineages: impact on structure, fitness, and sensitivity to nonnucleoside RT inhibitors

Natural evolution in primate lentiviral reverse transcriptase (RT) appears to have been constrained by the necessity to maintain function within an asymmetric protein composed of two identical primary amino acid sequences (66 kDa), of which one is cleaved (51 kDa). In this study, a detailed phylogenetic analysis now segregates groups O and M into clusters based on a cysteine or tyrosine residue located at position 181 of RT and linked to other signature residues. Divergent evolution of two group O (C181 or Y181) and the main (Y181 only) HIV-1 lineages did not appreciably impact RT activity or function. Group O RT structural models, based on group M subtype B RT crystal structures, revealed that most evolutionarily linked amino acids appear on a surface-exposed region of one subunit while in a noncatalytic RT pocket of the other subunit. This pocket binds nonnucleoside RT inhibitors (NNRTI); therefore, NNRTI sensitivity was used to probe enzyme differences in these group O and M lineages. In contrast to observations showing acquired drug resistance associated with fitness loss, the C181Y mutation in the C181 group O lineage resulted in a loss of intrinsic NNRTI resistance and was accompanied by fitness loss. Other mutations linked to the NNRTI-resistant C181 lineage also resulted in altered NNRTI sensitivity and a net fitness cost. Based on RT asymmetry and conservation of the intricate reverse transcription process, millions of years of divergent primate lentivirus evolution may be constrained to discrete mutations that appear primarily in the nonfunctional, solvent-accessible NNRTI binding pocket.

The non-nucleoside reverse transcriptase inhibitor efavirenz stimulates replication of human immunodeficiency virus type 1 harboring certain non-nucleoside resistance mutations

We measured the effects of non-nucleoside reverse transcriptase (RT) inhibitor-resistant mutations K101E+G190S, on replication fitness and EFV-resistance of HIV(NL4-3). K101E+G190S reduced fitness in the absence of EFV and increased EFV resistance, compared to either single mutant. Unexpectedly, K101E+G190S also replicated more efficiently in the presence of EFV than in its absence. Addition of the nucleoside resistance mutations L74V or M41L+T215Y to K101E+G190S improved fitness and abolished EFV-dependent stimulation of replication. D10, a clinical RT backbone containing M41L+T215Y and K101E+G190S, also demonstrated EFV-dependent stimulation that was dependent on the presence of K101E. These studies demonstrate that non-nucleoside reverse transcriptase inhibitors can stimulate replication of NNRTI-resistant HIV-1 and that nucleoside-resistant mutants can abolish this stimulation. The ability of EFV to stimulate NNRTI-resistant mutants may contribute to the selection of HIV-1 mutants in vivo. These studies have important implications regarding the treatment of HIV-1 with combination nucleoside and non-nucleoside therapies.

Ultrasensitive detection of minor drug-resistant variants for HIV after nevirapine exposure using allele-specific PCR: clinical significance

HIV-1 drug resistance mutations have been detected at low frequencies after single-dose nevirapine (sdNVP) for prevention of mother-to-child transmission (PMTCT). We investigated the relationship between these "minor variant" NVP-resistant viruses and clinical outcome with NVP-containing antiretroviral therapy (ART). An allele-specific quantitative PCR (ASPCR) assay was used to quantify the pre-ART frequency of K103N and Y181C in 26 women who had received sdNVP. The cohort was composed of 7 patients who experienced virologic failure and 19 control patients who maintained virologic suppression on NVP-containing ART; all were negative for resistance by standard genotyping. NVP resistance mutations were found in 17 of 26 (65%) patients using ASPCR. The frequency of NVP-resistant viruses ranged from 0.1% to 4.11%. Receiver operating characteristics (ROC) analysis identified a clinical threshold frequency of 0.19% for the ASPCR assay. Application of this threshold demonstrated minor variant resistance in 6 of 7 patients (86%) who failed treatment compared to 6 of 19 patients (32%) who were successful (OR = 13; 95% CI 1.27-133). ASPCR provides a means of detecting minor variant drug-resistant viruses that may impact subsequent treatment response. These data suggest a clinical role for highly sensitive assays to detect and quantify resistant viruses at low frequencies.

Low frequency nonnucleoside reverse-transcriptase inhibitor-resistant variants contribute to failure of efavirenz-containing regimens in treatment- experienced patients

BACKGROUND

The contribution of low frequency drug-resistant human immunodeficiency virus type 1 (HIV-1) variants to failure of antiretroviral therapy is not well defined in treatment-experienced patients. We sought to detect minor nonnucleoside reverse-transcriptase inhibitor (NNRTI)-resistant variants at the initiation of multidrug efavirenz-containing therapy in both NNRTI-naive and NNRTI-experienced patients and to determine their association with virologic response.

METHODS

Plasma samples at entry and at time of virologic failure from patients enrolled in the AIDS Clinical Trials Group study 398 were analyzed by standard genotype, single-genome sequencing and allele-specific polymerase chain reaction (K103N and Y181C) to detect and quantify minor NNRTI-resistant variants.

RESULTS

Minor populations of NNRTI-resistant variants that were missed by standard genotype were detected more often at study entry in NNRTI-experienced patients than NNRTI-naive patients by both single-genome sequencing (8 of 12 vs 3 of 15; P = .022) and allele-specific polymerase chain reaction (11% Y181C, 5 of 22 vs 3 of 72, respectively; P = .016). K103N variants at frequencies 11% were associated with inferior HIV-1 RNA response to efavirenz-containing therapy between entry and week 24 (change in HIV-1 RNA level, +0.5 vs -1.1 log(10) copies/mL; P < .001).

CONCLUSIONS

Minor NNRTI-resistant variants were more prevalent in NNRTI-experienced patients and were associated with reduced virologic response to efavirenz-containing multidrug regimens.

HIV Genetic Diversity and Drug Resistance

Most of the current knowledge on antiretroviral (ARV) drug development and resistance is based on the study of subtype B of HIV-1, which only accounts for 10% of the worldwide HIV infections. Cumulative evidence has emerged that different HIV types, groups and subtypes harbor distinct biological properties, including the response and susceptibility to ARV. Recent laboratory and clinical data highlighting such disparities are summarized in this review. Variations in drug susceptibility, in the emergence and selection of specific drug resistance mutations, in viral replicative capacity and in the dynamics of resistance acquisition under ARV selective pressure are discussed. Clinical responses to ARV therapy and associated confounding factors are also analyzed in the context of infections by distinct HIV genetic variants.

Infection with "escaped" virus variants impairs control of simian immunodeficiency virus SIVmac239 replication in Mamu-B*08-positive macaques

An understanding of the mechanism(s) by which some individuals spontaneously control human immunodeficiency virus (HIV)/simian immunodeficiency virus replication may aid vaccine design. Approximately 50% of Indian rhesus macaques that express the major histocompatibility complex (MHC) class I allele Mamu-B*08 become elite controllers after infection with simian immunodeficiency virus SIVmac239. Mamu-B*08 has a binding motif that is very similar to that of HLA-B27, a human MHC class I allele associated with the elite control of HIV, suggesting that SIVmac239-infected Mamu-B*08-positive (Mamu-B*08+) animals may be a good model for the elite control of HIV. The association with MHC class I alleles implicates CD8+ T cells and/or natural killer cells in the control of viral replication. We therefore introduced point mutations into eight Mamu-B*08-restricted CD8+ T-cell epitopes to investigate the contribution of epitope-specific CD8+ T-cell responses to the development of the control of viral replication. Ten Mamu-B*08+ macaques were infected with this mutant virus, 8X-SIVmac239. We compared immune responses and viral loads of these animals to those of wild-type SIVmac239-infected Mamu-B*08+ macaques. The five most immunodominant Mamu-B*08-restricted CD8+ T-cell responses were barely detectable in 8X-SIVmac239-infected animals. By 48 weeks postinfection, 2 of 10 8X-SIVmac239-infected Mamu-B*08+ animals controlled viral replication to <20,000 viral RNA (vRNA) copy equivalents (eq)/ml plasma, while 10 of 15 wild-type-infected Mamu-B*08+ animals had viral loads of <20,000 vRNA copy eq/ml (P = 0.04). Our results suggest that these epitope-specific CD8+ T-cell responses may play a role in establishing the control of viral replication in Mamu-B*08+ macaques.

Screening for hepatitis C virus non-nucleotide resistance mutations in treatment-naive women

OBJECTIVES

Hepatitis C virus (HCV) non-nucleoside inhibitors (NNIs) target the viral RNA-dependent RNA polymerase encoded by the NS5B gene. Several NNIs share a similar allosteric binding site, and their antiviral efficacy is attenuated by a cysteine-to-tyrosine mutation at amino acid 316 (C316Y). In the current study, we assessed NS5B resistance mutations in treatment-naive individuals from a prospective natural history study of viral infections in women.

METHODS

Partial NS5B sequences from HCV-positive women were amplified by RT-PCR. Additionally, subcloning was performed to evaluate intrapatient variability in selected samples.

RESULTS

HCV NS5B genotypes were 45 genotype 1a (57.0%), 11 genotype 1b (13.9%), 5 genotype 2a (6.3%), 3 genotype 2b (3.8%), 9 genotype 3a (11.4%) and 6 genotype 4a (7.6%). One HCV genotype 1a-infected patient was found to have the C316Y mutation (1.3%). Clonal analysis further revealed that all NS5B sequences from this individual--representing three serum samples collected 4 years apart--contained the C316Y mutation. In contrast, the S282T resistance mutation was not found in any samples.

CONCLUSIONS

The C316Y polymerase resistance mutation was found in 1.3% of samples from HCV-infected women. The presence of this mutation over time suggests significant replicative fitness of this variant and has implications for development of new specifically targeted antiviral therapies against HCV (STAT-C) targeting this region.

Select resistance-associated mutations in blood are associated with lower CSF viral loads and better neuropsychological performance

BACKGROUND

When antiretroviral therapy does not fully suppress HIV replication, suboptimal levels of antiretrovirals can select for antiretroviral resistant variants of HIV. These variants may exhibit reduced replication capacity and result in lower viral loads in blood. Our study evaluated whether antiretroviral resistance was associated with viral loads in the cerebrospinal fluid (CSF) and better neuropsychological (NP) performance.

METHODS

We enrolled 94 participants and each participant underwent a comprehensive neuromedical evaluation that used structured clinical assessments of medical history, ART and other medication use, comprehensive NP testing, and neurological and general physical signs of disease. Blood was collected by venipuncture, and all participants were offered lumbar puncture. Univariate and multivariate statistical methods were used to analyze the relationship between antiretroviral resistance, blood and CSF HIV RNA levels, substance use, and NP performance.

RESULTS

Antiretroviral resistance, detected in blood, was associated with lower CSF viral loads (p<0.01) and better NP performance (p=0.04) in multivariate analyses, independent of past and current ARV use and blood viral loads (model: p<0.01). However, HIV RNA levels in CSF did not independently correlate with NP performance. Low viral loads in the CSF limited our ability to investigate the relationship between antiretroviral resistance detected in CSF and NP performance.

CONCLUSIONS

Even in the absence of ART, antiretroviral resistance-associated mutations correlate with better NP performance possibly because these mutations reflect reduced neurovirulence compared with wild-type HIV.

Viral adaptation to an antiviral protein enhances the fitness level to above that of the uninhibited wild type

Viruses often evolve resistance to antiviral agents. While resistant strains are able to replicate in the presence of the agent, they generally exhibit lower fitness than the wild-type strain in the absence of the inhibitor. In some cases, resistant strains become dependent on the antiviral agent. However, the agent rarely, if ever, elevates dependent strain fitness above the uninhibited wild-type level. This would require an adaptive mechanism to convert the antiviral agent into a beneficial growth factor. Using an inhibitory scaffolding protein that specifically blocks phiX174 capsid assembly, we demonstrate that such mechanisms are possible. To obtain the quintuple-mutant resistant strain, the wild-type virus was propagated for approximately 150 viral life cycles in the presence of increasing concentrations of the inhibitory protein. The expression of the inhibitory protein elevated the strain's fitness significantly above the uninhibited wild-type level. Thus, selecting for resistance coselected for dependency, which was characterized and found to operate on the level of capsid nucleation. To the best of our knowledge, this is the first report of a virus evolving a mechanism to productively utilize an antiviral agent to stimulate its fitness above the uninhibited wild-type level. The results of this study may be predictive of the types of resistant phenotypes that could be selected by antiviral agents that specifically target capsid assembly.

Impact of novel human immunodeficiency virus type 1 reverse transcriptase mutations P119S and T165A on 4'-ethynylthymidine analog resistance profile

2',3'-Didehydro-3'-deoxy-4'-ethynylthymidine (4'-Ed4T), a derivative of stavudine (d4T), has potent activity against human immunodeficiency virus and is much less inhibitory to mitochondrial DNA synthesis and cell growth than its progenitor, d4T. 4'-Ed4T triphosphate was a better reverse transcriptase (RT) inhibitor than d4T triphosphate, due to the additional binding of the 4'-ethynyl group at a presumed hydrophobic pocket in the RT active site. Previous in vitro selection for 4'-Ed4T-resistant viral strains revealed M184V and P119S/T165A/M184V mutations on days 26 and 81, respectively; M184V and P119S/T165A/M184V conferred 3- and 130-fold resistance to 4'-Ed4T, respectively. We investigated the relative contributions of these mutations, engineered into the strain NL4-3 background, to drug resistance, RT activity, and viral growth. Viral variants with single RT mutations (P119S or T165A) did not show resistance to 4'-Ed4T; however, M184V and P119S/T165A/M184V conferred three- and fivefold resistance, respectively, compared with that of the wild-type virus. The P119S/M184V and T165A/M184V variants showed about fourfold resistance to 4'-Ed4T. The differences in the growth kinetics of the variants were not more than threefold. The purified RT of mutants with the P119S/M184V and T165A/M184V mutations were inhibited by 4'-Ed4TTP with 8- to 13-fold less efficiency than wild-type RT. M184V may be the primary resistance-associated mutation of 4'-Ed4T, and P119S and T165A are secondary mutations. On the basis of our findings and the results of structural modeling, a virus with a high degree of resistance to 4'-Ed4T (e.g., more than 50-fold resistance) will be difficult to develop. The previously observed 130-fold resistance of the virus with P119S/T165A/M184V to 4'-Ed4T may be partly due to mutations both in the RT sequence and outside the RT sequence.

Molecular basis of human immunodeficiency virus drug resistance: an update

Antiretroviral therapy has led to a significant decrease in human immunodeficiency virus (HIV)-related mortality. Approved antiretroviral drugs target different steps of the viral life cycle including viral entry (coreceptor antagonists and fusion inhibitors), reverse transcription (nucleoside and non-nucleoside inhibitors of the viral reverse transcriptase), integration (integrase inhibitors) and viral maturation (protease inhibitors). Despite the success of combination therapies, the emergence of drug resistance is still a major factor contributing to therapy failure. Viral resistance is caused by mutations in the HIV genome coding for structural changes in the target proteins that can affect the binding or activity of the antiretroviral drugs. This review provides an overview of the molecular mechanisms involved in the acquisition of resistance to currently used and promising investigational drugs, emphasizing the structural role of drug resistance mutations. The optimization of current antiretroviral drug regimens and the development of new drugs are still challenging issues in HIV chemotherapy. This article forms part of a special issue of Antiviral Research marking the 25th anniversary of antiretroviral drug discovery and development, Vol 85, issue 1, 2010.

Replicative capacity differences of thymidine analog resistance mutations in subtype B and C human immunodeficiency virus type 1

In order to understand the impact of zidovudine resistance and thymidine analog mutations (TAMs) on subtype C human immunodeficiency virus type 1, we created mutants in subtype C reverse transcriptase (RT). The subtype B RT was placed in a subtype C backbone to act as a control. Mutants and wild-type (WT) virus were competed in a head-to-head competition assay to determine how different clones grew in the same culture. Different viruses were distinguished by sequence tags in nef and a quantitative-PCR assay. The 67N and 70R accessory mutations gave an advantage over the WT in subtype C, but these mutations in subtype B had replication capacities similar to that of the WT. Of the triple mutants examined, the TAM-1 types, 41L210W215Y, were the most fit in both subtypes, but only in subtype C was the replication capacity the same as that of the WT. The TAM-2 mutants, 67N70R215F, had the slowest replication in both clones. The mixed TAM pathway mutant, 67N70R215Y, in subtype C had a significant advantage over the TAM-2 mutant, but this was not seen in subtype B. When the WT viruses were competed with each other, the subtype B RT had enhanced replication relative to subtype C. The increased capacities of the 67N and 70R mutations may indicate that there will be greater transmitted resistance and persistence in a subtype C setting than what is known for subtype B.

Variability in the Plasma Concentration of Efavirenz and Nevirapine is Associated with Genotypic Resistance after Treatment Interruption

Background

Selection and persistence of non-nucleoside reverse transcriptase inhibitor (NNRTI)-associated mutations during treatment interruptions (TIs) has been attributed to the long plasma half-life of these drugs. However, little is known about the contribution of variable NNRTI plasma levels before a TI. We evaluated the selection of NNRTI-related mutations and the coefficient of variation of NNRTI plasma concentrations during different TI periods.

Methods

The selection of NNRTI-related mutations was examined in 50 HIV type-1 (HIV-1)-infected patients on a virologically suppressive regimen who underwent TI guided by CD4+ T-cell counts and plasma viraemia. Population and clone-based sequencing of the reverse transcriptase coding region was performed using plasma HIV-1 RNA samples during TI and proviral DNA from peripheral blood mononuclear cells before TI. NNRTI plasma concentrations were determined by HPLC.

Results

In 7/50 treated patients, de novo and transient NNRTI-related mutations appeared when treatment was interrupted. Emergence of resistant variants (including K103N, Y181C or G190S) after interruption was associated with a higher coefficient of variation in NNRTI plasma concentrations during the treatment period. Moreover, minority HIV-1 variants containing different resistance patterns (V106I/A, K103R/E or Y188C/D/H) were detected regardless of NNRTI concentrations.

Conclusions

The emergence of NNRTI-associated mutations during TI appears to be associated with the variation of NNRTI plasma concentrations during the preceding treatment period. The selection of minority HIV-1 variants with different patterns of NNRTI resistance in the absence of drug pressure should be considered for the efficacy of future NNRTI-containing antiretroviral regimens.

Quantification of simian immunodeficiency virus cytotoxic T lymphocyte escape mutant viruses

Escape from cytotoxic T-lymphocyte (CTL) pressure is common in HIV-1 infection of humans and simian immunodeficiency virus (SIV) infections of macaques. CTL escape typically incurs a fitness cost as reversion back to wild-type can occur upon transmission. We utilized sequence-specific primers and DNA probes with real-time polymerase chain reaction (PCR) to sensitively and specifically track wild-type and escape mutant viremia at the Mane-A*17-restricted SIV Gag(371379) epitope AF9 in pigtail macaques. The generation of minor escape mutant populations is detected by the real-time PCR 2 weeks earlier than observed using standard sequencing techniques. We passaged the AF9 CTL escape mutant virus into two naïve Mane-A*17-negative pigtail macaques and showed that reversion to wild-type was rapid during acute infection and then slowed considerably at later stages of the infection. These data help refine our understanding of how CTL escape mutant viruses evolve.

Optimal timing and duration of induction therapy for HIV-1 infection

The tradeoff between the need to suppress drug-resistant viruses and the problem of treatment toxicity has led to the development of various drug-sparing HIV-1 treatment strategies. Here we use a stochastic simulation model for viral dynamics to investigate how the timing and duration of the induction phase of induction–maintenance therapies might be optimized. Our model suggests that under a variety of biologically plausible conditions, 6–10 mo of induction therapy are needed to achieve durable suppression and maximize the probability of eradicating viruses resistant to the maintenance regimen. For induction regimens of more limited duration, a delayed-induction or -intensification period initiated sometime after the start of maintenance therapy appears to be optimal. The optimal delay length depends on the fitness of resistant viruses and the rate at which target-cell populations recover after therapy is initiated. These observations have implications for both the timing and the kinds of drugs selected for induction–maintenance and therapy-intensification strategies.

Clinicians treating HIV infection must balance the need to suppress viral replication against the harmful side effects and significant cost of antiretroviral therapy. Inadequate therapy often results in the emergence of resistant viruses and treatment failure. These difficulties are especially acute in resource-poor settings, where antiretroviral agents are limited. This has prompted an interest in induction–maintenance (IM) treatment strategies, in which brief intensive therapy is used to reduce host viral levels. Induction is followed by a simplified and more easily tolerated maintenance regimen. IM approaches remain an unproven concept in HIV therapy. We have developed a mathematical model to simulate clinical responses to antiretroviral drug therapy. We account for latent infection, partial drug efficacy, cross-resistance, viral recombination, and other factors. This model accurately reflects expected outcomes under single, double, and standard three-drug antiretroviral therapy. When applied to IM therapy, we find that (1) IM is expected to be successful beyond 3 y under a variety of conditions; (2) short-term induction therapy is optimally started several days to weeks after the start of maintenance; and (3) IM therapy may eradicate some preexisting drug-resistant viral strains from the host. Our simulations may help develop new treatment strategies and optimize future clinical trials.

Modeling and estimation of kinetic parameters and replicative fitness of HIV-1 from flow-cytometry-based growth competition experiments

Growth competition assays have been developed to quantify the relative fitness of HIV-1 mutants. In this article, we develop mathematical models to describe viral/cellular dynamic interactions in the assay system from which the competitive fitness indices or parameters are defined. In our previous HIV-viral fitness experiments, the concentration of uninfected target cells was assumed to be constant (Wu et al. 2006). But this may not be true in some experiments. In addition, dual infection may frequently occur in viral fitness experiments and may not be ignorable. Here, we relax these two assumptions and extend our earlier viral fitness model (Wu et al. 2006). The resulting models then become nonlinear ODE systems for which closed-form solutions are not achievable. In the new model, the viral relative fitness is a function of time since it depends on the target cell concentration. First, we studied the structure identifiability of the nonlinear ODE models. The identifiability analysis showed that all parameters in the proposed models are identifiable from the flow-cytometry-based experimental data that we collected. We then employed a global optimization approach (the differential evolution algorithm) to directly estimate the kinetic parameters as well as the relative fitness index in the nonlinear ODE models using nonlinear least square regression based on the experimental data. Practical identifiability was investigated via Monte Carlo simulations.