HIV-1 reverse transcriptase inhibitor resistance mutations and fitness: a view from the clinic and ex vivo

Adapting the rhizome concept to an extended definition of viral quasispecies and the implications for molecular evolution

The rhizome concept proposed by Gilles Deleuze and Félix Guattari offers a novel perspective on the organization and interdependence of complex constellations of heterogeneous entities, their mapping and their ruptures. The emphasis of the present study is placed on the dynamics of contacts and communication among such entities that arise from experimentation, without any favored hierarchy or origin. When applied to biological evolution, the rhizome concept integrates all types of heterogeneity resulting from "symbiotic" relationships among living beings (or their genomic material), horizontal genetic transfer, recombination and mutation, and breaks away from the approach that gives rise to the phylogenetic tree of life. It has already been applied to describe the dynamics and evolution of RNA viruses. Thus, here we introduce a novel framework for the interpretation the viral quasispecies concept, which explains the evolution of RNA virus populations as the result of dynamic interconnections and multifaceted interdependence between highly heterogeneous viral sequences and its inherently heterogeneous host cells. The rhizome network perspective underlines even further the medical implications of the broad mutant spectra of viruses that are in constant flow, given the multiple pathways they have available for fitness loss and gain.

Target-based drug design strategies to overcome resistance to antiviral agents: opportunities and challenges

Viral infections have a major impact in human health. Ongoing viral transmission and escalating selective pressure have the potential to favor the emergence of vaccine- and antiviral drug-resistant viruses. Target-based approaches for the design of antiviral drugs can play a pivotal role in combating drug-resistant challenges. Drug design computational tools facilitate the discovery of novel drugs. This review provides a comprehensive overview of current drug design strategies employed in the field of antiviral drug resistance, illustrated through the description of a series of successful applications. These strategies include technologies that enhance compound-target affinity while minimizing interactions with mutated binding pockets. Furthermore, emerging approaches such as virtual screening, targeted protein/RNA degradation, and resistance analysis during drug design have been harnessed to curtail the emergence of drug resistance. Additionally, host targeting antiviral drugs offer a promising avenue for circumventing viral mutation. The widespread adoption of these refined drug design strategies will effectively address the prevailing challenge posed by antiviral drug resistance.

Viral Fitness, Population Complexity, Host Interactions, and Resistance to Antiviral Agents

Fitness of viruses has become a standard parameter to quantify their adaptation to a biological environment. Fitness determinations for RNA viruses (and some highly variable DNA viruses) meet with several uncertainties. Of particular interest are those that arise from mutant spectrum complexity, absence of population equilibrium, and internal interactions among components of a mutant spectrum. Here, concepts, fitness measurements, limitations, and current views on experimental viral fitness landscapes are discussed. The effect of viral fitness on resistance to antiviral agents is covered in some detail since it constitutes a widespread problem in antiviral pharmacology, and a challenge for the design of effective antiviral treatments. Recent evidence with hepatitis C virus suggests the operation of mechanisms of antiviral resistance additional to the standard selection of drug-escape mutants. The possibility that high replicative fitness may be the driver of such alternative mechanisms is considered. New broad-spectrum antiviral designs that target viral fitness may curtail the impact of drug-escape mutants in treatment failures. We consider to what extent fitness-related concepts apply to coronaviruses and how they may affect strategies for COVID-19 prevention and treatment.

High Rate of HIV-1 Drug Resistance in Antiretroviral Therapy-Failure Patients in Liaoning Province, China

This study aimed to monitor the prevalence of HIV-1 drug resistance and risk factors associated with drug resistance in antiretroviral therapy (ART)-failure individuals in Liaoning Province, China. Plasma samples were collected from HIV-1-positive individuals who experienced ART failure in Liaoning Province between April 2018 and September 2019. Genotype resistance test was performed using an in-house assay on these collected samples. Factors associated with drug resistance were identified by logistic regression analysis. We collected a total of 468 ART-failure individuals, of which 256 were successfully included in the final study. Of these, the most predominant genotype was CRF01_AE, accounting for 77.73%. The resistance rate to any of the three classes of antiretroviral drugs (non-nucleoside reverse transcriptase inhibitors [NNRTIs], nucleoside reverse transcriptase inhibitors [NRTIs], and protease inhibitors [PIs]) was 64.84%. Among 256 ART-failure patients, 62.89% showed drug resistance to NNRTIs, 50.39% to NRTIs, and 3.13% to PIs. G190S (31.25%) and Y181C (25.78%) mutations were the most common NNRTIs resistance mutations. K65R (29.69%), M184V (28.52%) were the most common NRTIs resistance mutations. Factors associated with drug resistance included current ART regimen and viral load. The high drug resistance rate among ART-failure individuals in Liaoning Province needs more attention. Corresponding strategies for the risk factors associated with HIV drug resistance can better control and prevent the prevalence of resistance.

MrHAMER yields highly accurate single molecule viral sequences enabling analysis of intra-host evolution

Technical challenges remain in the sequencing of RNA viruses due to their high intra-host diversity. This bottleneck is particularly pronounced when interrogating long-range co-evolved genetic interactions given the read-length limitations of next-generation sequencing platforms. This has hampered the direct observation of these genetic interactions that code for protein-protein interfaces with relevance in both drug and vaccine development. Here we overcome these technical limitations by developing a nanopore-based long-range viral sequencing pipeline that yields accurate single molecule sequences of circulating virions from clinical samples. We demonstrate its utility in observing the evolution of individual HIV Gag-Pol genomes in response to antiviral pressure. Our pipeline, called Multi-read Hairpin Mediated Error-correction Reaction (MrHAMER), yields >1000s of viral genomes per sample at 99.9% accuracy, maintains the original proportion of sequenced virions present in a complex mixture, and allows the detection of rare viral genomes with their associated mutations present at <1% frequency. This method facilitates scalable investigation of genetic correlates of resistance to both antiviral therapy and immune pressure and enables the identification of novel host-viral and viral-viral interfaces that can be modulated for therapeutic benefit.

A Novel High Throughput, Parallel Infection Assay for Determining the Replication Capacities of 346 Primary HIV-1 Isolates of the Zurich Primary HIV-1 Infection Study in Primary Cells

HIV-1 replication capacity is an important characteristic to understand the replication competence of single variants or virus populations. It can further aid in the understanding of HIV-1 pathogenicity, disease progression, and drug resistance mutations. To effectively study RC, many assays have been established. However, there is still demand for a high throughput replication capacity assay using primary cells which is robust and reproducible. In this study, we established such an assay and validated it using 346 primary HIV-1 isolates from patients enrolled in the Zurich Primary HIV Infection study (ZPHI) and two control viruses, HIV-1 JR-CSFWT and HIV-1 JR-CSFK65R_M184V. Replication capacity was determined by measuring the viral growth on PBMCs over 10 days by longitudinally transferring cell culture supernatant to TZM-bl reporter cells. By utilizing the TZM-bl luciferase reporter assay, we determined replication capacity by measuring viral infectivity. The simplicity of the experimental setup allowed for all 346 primary HIV-1 isolates to be replicated at one time. Although the infectious input dose for each virus was normalized, a broad range of replication capacity values over 4 logs was observed. The approach was confirmed by two repeated experiments and we demonstrated that the reproducibility of the replication capacity values is statistically comparable between the two separate experiments. In summary, these results endorse our high throughput replication capacity assay as reproducible and robust and can be utilized for large scale HIV-1 replication capacity experiments in primary cells.

Evolutionary Genetics of Mycobacterium tuberculosis and HIV-1: "The Tortoise and the Hare"

The already enormous burden caused by Mycobacterium tuberculosis and Human Immunodeficiency Virus type 1 (HIV-1) alone is aggravated by co-infection. Despite obvious differences in the rate of evolution comparing these two human pathogens, genetic diversity plays an important role in the success of both. The extreme evolutionary dynamics of HIV-1 is in the basis of a robust capacity to evade immune responses, to generate drug-resistance and to diversify the population-level reservoir of M group viral subtypes. Compared to HIV-1 and other retroviruses, M. tuberculosis generates minute levels of genetic diversity within the host. However, emerging whole-genome sequencing data show that the M. tuberculosis complex contains at least nine human-adapted phylogenetic lineages. This level of genetic diversity results in differences in M. tuberculosis interactions with the host immune system, virulence and drug resistance propensity. In co-infected individuals, HIV-1 and M. tuberculosis are likely to co-colonize host cells. However, the evolutionary impact of the interaction between the host, the slowly evolving M. tuberculosis bacteria and the HIV-1 viral "mutant cloud" is poorly understood. These evolutionary dynamics, at the cellular niche of monocytes/macrophages, are also discussed and proposed as a relevant future research topic in the context of single-cell sequencing.

Purifying Selection in Human Immunodeficiency Virus-1 pol Gene in Perinatally Human Immunodeficiency Virus-1-Infected Children Harboring Discordant Immunological Response and Virological Nonresponse to Long-Term Antiretroviral Therapy

Background

Biological monitoring of antiretroviral treatment (ART) in human immunodeficiency virus (HIV)-infected pediatric population remains challenging. The aim of the present study was to assess the long-term HIV-1 genetic diversity in pol gene in HIV-1-infected children in virological failure under antiretroviral regimen adapted according to the successive World Health Organization (WHO) guidelines for resource-constrained settings.

Methods

HIV-1 diversity in pol gene was assessed in HIV-1-infected children and adolescents born from HIV-infected mothers (median age at follow-up: 13.8 years) in virological failure (VF⁺) despite long-term regimen recommended by the WHO. The numbers of nonsynonymous substitutions per potential nonsynonymous site (dN) and of synonymous substitutions at potential synonymous sites (dS) in HIV-1 pol gene and the dN/dS ratios were used to estimate the selective pressure on circulating HIV-1.

Results

The immunological responses to ART basically corresponded to: 1) Full therapeutic failure with immunological (I^-) and virological nonresponses in one-quarter (24.6%) of study children ((I^-, VF⁺) subgroup); 2) Discordant immunovirological responses with paradoxical high CD4 T cell counts (I⁺) and high HIV-1 RNA load in the remaining cohort patients (75.4%) ((I⁺, VF⁺) subgroup). The mean dS was 1.8-fold higher in (I⁺, VF⁺) than (I^-, VF⁺) subgroup (25.9 ± 18.4 vs. 14.3 ± 10.8). In the (I⁺, VF⁺) subgroup, the mean dS was 1.6-fold higher than the mean dN. Finally, the mean dN/dS ratio was 2.1-fold lower in (I⁺, VF⁺) than (I^-, VF⁺) subgroup (0.6 ± 0.3 vs. 1.3 ± 0.7), indicating purifying selection in the immunovirological discordant (I⁺, VF⁺) subgroup and positive selection in the immunovirological failure (I^-, VF⁺) subgroup.

Conclusions

Children and adolescents in immunovirological therapeutic failure harbor positive selection of HIV-1 strains favoring diversifying in pol-encoded amino acids. In contrast, children with persistent discordant immunovirological responses show accumulation of mutations and purifying selection in pol gene sequences, indicating limited genetic evolution and likely suggesting genetic adaptation of viruses to host functional constraints.

Escalating and sustained immunovirological dissociation among antiretroviral drug-experienced perinatally human immunodeficiency virus-1-infected children and adolescents living in the Central African Republic: A STROBE-compliant study

Sub-Saharan Africa has the vast majority (∼90%) of new pediatric acquired immunodeficiency syndrome cases worldwide. Biologically monitoring HIV-infected pediatric populations remains challenging. The differential interest of human immunodeficiency virus (HIV)-1 RNA loads and CD4 T-cell counts is debated for the treatment of pediatric acquired immunodeficiency syndrome patients.Long-term antiretroviral treatment (ART) outcomes regarding immunological and virological surrogate markers were longitudinally evaluated between 2009 and 2014 (over 57 months) in 245 perinatally HIV-1-infected children and adolescents born from HIV-infected mothers, treated at inclusion for at least 6 months by the World Health Organization-recommended ART in Bangui, Central African Republic.Patients were monitored over time biologically for CD4 T-cell counts, HIV-1 RNA loads, and drug resistance mutation genotyping.Children lost to follow-up totaled 6%. Four categories of immunovirological responses to ART were observed. At baseline, therapeutic success with sustained immunological and virological responses was observed in 80 (32.6%) children; immunological and virologic nonresponses occurred in 32 (13.0%) children; finally, the majority (133; 54.2%) of the remaining children showed discordant immunovirological responses. Among them, 33 (13.4%) children showed rapid virological responses to ART with an undetectable viral load, whereas immunological responses remained absent after 6 months of treatment and increased progressively over time in most of the cases, suggesting slow immunorestoration. Notably, nearly half of the children (40.8% at baseline and 48.2% at follow-up) harbored discordant immunovirological responses with a paradoxically high CD4 T-cell count and HIV-1 RNA load, which are always associated with high levels of drug resistance mutations. The latter category showed a significant increase over time, with a growth rate of 1.23% per year of follow-up.Our STROBE-compliant study demonstrates the high heterogeneity of biological responses under ART in children with frequent passage from 1 category to another over time. Close biological evaluation with access to routine plasma HIV-1 RNA load monitoring is crucial for adapting the complex outcomes of ART in HIV-infected children born from infected mothers.

Viral fitness: history and relevance for viral pathogenesis and antiviral interventions

The quasispecies dynamics of viral populations (continuous generation of variant genomes and competition among them) has as one of its frequent consequences variations in overall multiplication capacity, a major component of viral fitness. This parameter has multiple implications for viral pathogenesis and viral disease control, some of them unveiled thanks to deep sequencing of viral populations. Darwinian fitness is an old concept whose quantification dates back to the early developments of population genetics. It was later applied to viruses (mainly to RNA viruses) to quantify relative multiplication capacities of individual mutant clones or complex populations. The present article reviews the fitness concept and its relevance for the understanding of the adaptive dynamics of viruses in constant and changing environments. Many studies have addressed the fitness cost of escape mutations (to antibodies, cytotoxic T cells or inhibitors) as an influence on the efficacy of antiviral interventions. Here, we summarize the evidence that the basal fitness level can be a determinant of inhibitor resistance.

Quasispecies dynamics in disease prevention and control

Medical interventions to prevent and treat viral disease constitute evolutionary forces that may modify the genetic composition of viral populations that replicate in an infected host and influence the genomic composition of those viruses that are transmitted and progress at the epidemiological level. Given the adaptive potential of viruses in general and the RNA viruses in particular, the selection of viral mutants that display some degree of resistance to inhibitors or vaccines is a tangible challenge. Mutant selection may jeopardize control of the viral disease. Strategies intended to minimize vaccination and treatment failures are proposed and justified based on fundamental features of viral dynamics explained in the preceding chapters. The recommended use of complex, multiepitopic vaccines, and combination therapies as early as possible after initiation of infection falls under the general concept that complexity cannot be combated with simplicity. It also follows that sociopolitical action to interrupt virus replication and spread as soon as possible is as important as scientifically sound treatment designs to control viral disease on a global scale.

HIV-1 Reverse Transcriptase Promotes Tumor Growth and Metastasis Formation via ROS-Dependent Upregulation of Twist

HIV-induced immune suppression results in the high prevalence of HIV/AIDS-associated malignancies including Kaposi sarcoma, non-Hodgkin lymphoma, and cervical cancer. HIV-infected people are also at an increased risk of “non-AIDS-defining” malignancies not directly linked to immune suppression but associated with viral infections. Their incidence is increasing despite successful antiretroviral therapy. The mechanism behind this phenomenon remains unclear. Here, we obtained daughter clones of murine mammary gland adenocarcinoma 4T1luc2 cells expressing consensus reverse transcriptase of HIV-1 subtype A FSU_A strain (RT_A) with and without primary mutations of drug resistance. In in vitro tests, mutations of resistance to nucleoside inhibitors K65R/M184V reduced the polymerase, and to nonnucleoside inhibitors K103N/G190S, the RNase H activities of RT_A. Expression of these RT_A variants in 4T1luc2 cells led to increased production of the reactive oxygen species (ROS), lipid peroxidation, enhanced cell motility in the wound healing assay, and upregulation of expression of Vimentin and Twist. These properties, particularly, the expression of Twist, correlated with the levels of expression RT_A and/or the production of ROS. When implanted into syngeneic BALB/C mice, 4T1luc2 cells expressing nonmutated RT_A demonstrated enhanced rate of tumor growth and increased metastatic activity, dependent on the level of expression of RT_A and Twist. No enhancement was observed for the clones expressing mutated RT_A variants. Plausible mechanisms are discussed involving differential interactions of mutated and nonmutated RTs with its cellular partners involved in the regulation of ROS. This study establishes links between the expression of HIV-1 RT, production of ROS, induction of EMT, and enhanced propagation of RT-expressing tumor cells. Such scenario can be proposed as one of the mechanisms of HIV-induced/enhanced carcinogenesis not associated with immune suppression.

The increasing impact of lethal mutagenesis of viruses

Selection of viral mutants resistant to compounds used in therapy is a major determinant of treatment failure, a problem akin to antibiotic resistance in bacteria. In this scenario, mutagenic base and nucleoside analogs have entered the picture because they increase the mutation rate of viral populations to levels incompatible with their survival. This antiviral strategy is termed lethal mutagenesis. It has found a major impulse with the observation that some antiviral agents, which initially were considered only inhibitors of virus multiplication, may in effect exert part of their antiviral activity through mutagenesis. Here, we review the conceptual basis of lethal mutagenesis, the evidence of virus extinction through mutagenic nucleotide analogs and prospects for application in antiviral designs.

Natural control of HIV infection in young women in South Africa: HPTN 068

BACKGROUND

Some individuals control HIV replication without antiretroviral (ARV) therapy.

OBJECTIVE

To analyze viral suppression in young women in rural South Africa enrolled in a trial evaluating a behavioral intervention for HIV prevention.

METHODS

Plasma samples were obtained from women ages 13-24 (81 infected at enrollment, 164 seroconverters). ARV testing was performed using an assay that detects 20 ARV drugs. Women were classified as viremic controllers if they were virally suppressed for ≥12 months with no ARV drug use.

RESULTS

Samples from 216/245 (88.2%) women had no ARV drugs detected at their first HIV-positive visit. Thirty-four (15.7%) of the 216 women had a viral load <2,000 copies/mL. Fifteen of the 34 women were followed for ≥12 months; 12 were virally suppressed with no ARV drugs detected during follow-up. These women were classified as viremic controllers (overall: 12/216 = 5.6%). The median CD4 cell count at the first HIV-positive visit was higher among the 12 controllers than among the 204 women who were not using ARV drugs (759 vs. 549 cells/mm³, p = 0.02). Some women had a viral load <40 copies/mL at a single study visit, but none were classified as elite controllers (viral load <40 copies/mL for ≥12 months with no ARV drug use).

CONCLUSIONS

In this cohort, 5.6% of women who were not using ARV drugs had sustained viral suppression. This represents a minimum estimate of the frequency of viremic controllers in this cohort, since some women were not followed long enough to meet the criteria for classification.

Estimating the mutational fitness effects distribution during early HIV infection

The evolution of HIV during acute infection is often considered a neutral process. Recent analysis of sequencing data from this stage of infection, however, showed high levels of shared mutations between independent viral populations. This suggests that selection might play a role in the early stages of HIV infection. We adapted an existing model for random evolution during acute HIV-infection to include selection. Simulations of this model were used to fit a global mutational fitness effects distribution to previously published sequencing data of the env gene of individuals with acute HIV infection. Measures of sharing between viral populations were used as summary statistics to compare the data to the simulations. We confirm that evolution during acute infection is significantly different from neutral. The distribution of mutational fitness effects is best fit by a distribution with a low, but significant fraction of beneficial mutations and a high fraction of deleterious mutations. While most mutations are neutral or deleterious in this model, about 5% of mutations are beneficial. These beneficial mutations will, on average, result in a small but significant increase in fitness. When assuming no epistasis, this indicates that, at the moment of transmission, HIV is near, but not on the fitness peak for early infection.

Resistance of high fitness hepatitis C virus to lethal mutagenesis

Viral fitness quantifies the degree of virus adaptation to a given environment. How viral fitness can influence the mutant spectrum complexity of a viral quasispecies subjected to lethal mutagenesis has not been investigated. Here we document that two high fitness hepatitis C virus populations display higher resistance to the mutagenic nucleoside analogues favipiravir and ribavirin than their parental, low fitness HCV. All populations, however, exhibited a mutation transition bias indicative of active mutagenesis. Resistance to the analogues was associated with a limited expansion of mutant spectrum complexity, as evidenced by several diversity indices used to characterize mutant spectra. The results are consistent with a replicative site-drug competition mechanism that was previously proposed for HCV fitness-associated resistance to non-mutagenic inhibitors. Other alternative, non-mutually exclusive mechanisms are considered. The results introduce viral fitness as a relevant parameter to evaluate the response of viruses to lethal mutagenesis, with implications for antiviral designs.

Viral Fitness Correlates with the Magnitude and Direction of the Perturbation Induced in the Host's Transcriptome: The Tobacco Etch Potyvirus-Tobacco Case Study

Determining the fitness of viral genotypes has become a standard practice in virology as it is essential to evaluate their evolutionary potential. Darwinian fitness, defined as the advantage of a given genotype with respect to a reference one, is a complex property that captures, in a single figure, differences in performance at every stage of viral infection. To what extent does viral fitness result from specific molecular interactions with host factors and regulatory networks during infection? Can we identify host genes in functional classes whose expression depends on viral fitness? Here, we compared the transcriptomes of tobacco plants infected with seven genotypes of tobacco etch potyvirus that differ in fitness. We found that the larger the fitness differences among genotypes, the more dissimilar the transcriptomic profiles are. Consistently, two different mutations, one in the viral RNA polymerase and another in the viral suppressor of RNA silencing, resulted in significantly similar gene expression profiles. Moreover, we identified host genes whose expression showed a significant correlation, positive or negative, with the virus' fitness. Differentially expressed genes which were positively correlated with viral fitness activate hormone- and RNA silencing-mediated pathways of plant defense. In contrast, those that were negatively correlated with fitness affect metabolism, reducing growth, and development. Overall, these results reveal the high information content of viral fitness and suggest its potential use to predict differences in genomic profiles of infected hosts.

HIV Resistance Prediction to Reverse Transcriptase Inhibitors: Focus on Open Data

Research and development of new antiretroviral agents are in great demand due to issues with safety and efficacy of the antiretroviral drugs. HIV reverse transcriptase (RT) is an important target for HIV treatment. RT inhibitors targeting early stages of the virus-host interaction are of great interest for researchers. There are a lot of clinical and biochemical data on relationships between the occurring of the single point mutations and their combinations in the pol gene of HIV and resistance of the particular variants of HIV to nucleoside and non-nucleoside reverse transcriptase inhibitors. The experimental data stored in the databases of HIV sequences can be used for development of methods that are able to predict HIV resistance based on amino acid or nucleotide sequences. The data on HIV sequences resistance can be further used for (1) development of new antiretroviral agents with high potential for HIV inhibition and elimination and (2) optimization of antiretroviral therapy. In our communication, we focus on the data on the RT sequences and HIV resistance, which are available on the Internet. The experimental methods, which are applied to produce the data on HIV-1 resistance, the known data on their concordance, are also discussed.

Quantifying the fitness cost of HIV-1 drug resistance mutations through phylodynamics

Drug resistant HIV is a major threat to the long-term efficacy of antiretroviral treatment. Around 10% of ART-naïve patients in Europe are infected with drug-resistant HIV type 1. Hence it is important to understand the dynamics of transmitted drug resistance evolution. Thanks to routinely performed drug resistance tests, HIV sequence data is increasingly available and can be used to reconstruct the phylogenetic relationship among viral lineages. In this study we employ a phylodynamic approach to quantify the fitness costs of major resistance mutations in the Swiss HIV cohort. The viral phylogeny reflects the transmission tree, which we model using stochastic birth–death-sampling processes with two types: hosts infected by a sensitive or resistant strain. This allows quantification of fitness cost as the ratio between transmission rates of hosts infected by drug resistant strains and transmission rates of hosts infected by drug sensitive strains. The resistance mutations 41L, 67N, 70R, 184V, 210W, 215D, 215S and 219Q (nRTI-related) and 103N, 108I, 138A, 181C, 190A (NNRTI-related) in the reverse trancriptase and the 90M mutation in the protease gene are included in this study. Among the considered resistance mutations, only the 90M mutation in the protease gene was found to have significantly higher fitness than the drug sensitive strains. The following mutations associated with resistance to reverse transcriptase inhibitors were found to be less fit than the sensitive strains: 67N, 70R, 184V, 219Q. The highest posterior density intervals of the transmission ratios for the remaining resistance mutations included in this study all included 1, suggesting that these mutations do not have a significant effect on viral transmissibility within the Swiss HIV cohort. These patterns are consistent with alternative measures of the fitness cost of resistance mutations. Overall, we have developed and validated a novel phylodynamic approach to estimate the transmission fitness cost of drug resistance mutations.

The introduction of antiretroviral therapy (ART) has decreased mortality and morbidity rates among HIV-infected people, and improved their quality of life. In fact, the WHO states that antiretroviral therapy programmes averted an estimated 7.8 million deaths worldwide between 2000 and 2014. However, the antiretroviral regimen prescribed to a patient may be unable to control HIV infection. Factors that can contribute to treatment failure include drug resistance, drug toxicity, or poor treatment adherence. In this study we aim to understand the dynamics of transmitted drug resistance by analysing the viral sequence data that was collected for resistance testing. We present a novel approach to quantify how drug resistance impacts virus lineage transmissibility, how fast resistance mutations evolve in sensitive strains and how fast they revert back to the sensitive type. We apply our approach to the Swiss HIV cohort study, and obtain patterns of viral transmission fitness that are consistent with alternative, harder to obtain measures of fitness.

Effects of Mutations on Replicative Fitness and Major Histocompatibility Complex Class I Binding Affinity Are Among the Determinants Underlying Cytotoxic-T-Lymphocyte Escape of HIV-1 Gag Epitopes

Certain “protective” major histocompatibility complex class I (MHC-I) alleles, such as B*57 and B*27, are associated with long-term control of HIV-1 in vivo mediated by the CD8⁺ cytotoxic-T-lymphocyte (CTL) response. However, the mechanism of such superior protection is not fully understood. Here we combined high-throughput fitness profiling of mutations in HIV-1 Gag, in silico prediction of MHC-peptide binding affinity, and analysis of intraperson virus evolution to systematically compare differences with respect to CTL escape mutations between epitopes targeted by protective MHC-I alleles and those targeted by nonprotective MHC-I alleles. We observed that the effects of mutations on both viral replication and MHC-I binding affinity are among the determinants of CTL escape. Mutations in Gag epitopes presented by protective MHC-I alleles are associated with significantly higher fitness cost and lower reductions in binding affinity with respect to MHC-I. A linear regression model accounting for the effect of mutations on both viral replicative capacity and MHC-I binding can explain the protective efficacy of MHC-I alleles. Finally, we found a consistent pattern in the evolution of Gag epitopes in long-term nonprogressors versus progressors. Overall, our results suggest that certain protective MHC-I alleles allow superior control of HIV-1 by targeting epitopes where mutations typically incur high fitness costs and small reductions in MHC-I binding affinity.

Understanding the mechanism of viral control achieved in long-term nonprogressors with protective HLA alleles provides insights for developing functional cure of HIV infection. Through the characterization of CTL escape mutations in infected persons, previous researchers hypothesized that protective alleles target epitopes where escape mutations significantly reduce viral replicative capacity. However, these studies were usually limited to a few mutations observed in vivo. Here we utilized our recently developed high-throughput fitness profiling method to quantitatively measure the fitness of mutations across the entirety of HIV-1 Gag. The data enabled us to integrate the results with in silico prediction of MHC-peptide binding affinity and analysis of intraperson virus evolution to systematically determine the differences in CTL escape mutations between epitopes targeted by protective HLA alleles and those targeted by nonprotective HLA alleles. We observed that the effects of Gag epitope mutations on HIV replicative fitness and MHC-I binding affinity are among the major determinants of CTL escape.

How to win the HIV-1 drug resistance hurdle race: running faster or jumping higher?

Infections by the human immunodeficiency virus type 1 (HIV-1), the causative agent of the acquired immunodeficiency syndrome (AIDS), are still totaling an appalling 36.7 millions worldwide, with 1.1 million AIDS deaths/year and a similar number of yearly new infections. All this, in spite of the discovery of HIV-1 as the AIDS etiological agent more than 30 years ago and the introduction of an effective combinatorial antiretroviral therapy (cART), able to control disease progression, more than 20 years ago. Although very effective, current cART is plagued by the emergence of drug-resistant viral variants and most of the efforts in the development of novel direct-acting antiviral agents (DAAs) against HIV-1 have been devoted toward the fighting of resistance. In this review, rather than providing a detailed listing of all the drugs and the corresponding resistance mutations, we aim, through relevant examples, at presenting to the general reader the conceptual shift in the approaches that are being taken to overcome the viral resistance hurdle. From the classic 'running faster' strategy, based on the development of novel DAAs active against the mutant viruses selected by the previous drugs and/or presenting to the virus a high genetic barrier toward the development of resilience, to a 'jumping higher' approach, which looks at the cell, rather than the virus, as a source of valuable drug targets, in order to make the cellular environment non-permissive toward the replication of both wild-type and mutated viruses.

Drug resistance-related mutations T369V/I in the connection subdomain of HIV-1 reverse transcriptase severely impair viral fitness

Fitness is a key parameter in the measurement of transmission capacity of individual drug-resistant HIV. Drug-resistance related mutations (DRMs) T369V/I and A371V in the connection subdomain (CN) of reverse transcriptase (RT) occur at higher frequencies in the individuals experiencing antiretroviral therapy failure. Here, we evaluated the effects of T369V/I and A371V on viral fitness, in the presence or in the absence of thymidine analogue resistance-associated mutations (TAMs) and assessed the effect of potential RT structure-related mechanism on change in viral fitness. Mutations T369V/I, A371V, alone or in combination with TAMs were introduced into a modified HIV-1 infectious clone AT1 by site-directed mutagenesis. Then, experiments on mutant and wild-type virus AT2 were performed separately using a growth-competition assay, and then the relative fitness was calculated. Structural analysis of RT was conducted using Pymol software. Results showed that T369V/I severely impaired the relative virus fitness, and A371V compensated for the viral fitness reduction caused by TAMs. Structural modeling of RT suggests that T369V/I substitutions disrupt powerful hydrogen bonds formed by T369 and V365 in p51 and p66. This study indicates that the secondary DRMs within CN might efficiently damage viral fitness, and provides valuable information for clinical surveillance and prevention of HIV-1 strains carrying these DRMs.

Prevalence of Transmitted Drug-Resistance Mutations and Polymorphisms in HIV-1 Reverse Transcriptase, Protease, and gp41 Sequences Among Recent Seroconverters in Southern Poland

Background

Monitoring of drug resistance-related mutations among patients with recent HIV-1 infection offers an opportunity to describe current patterns of transmitted drug resistance (TDR) mutations.

Material/Methods

Of 298 individuals newly diagnosed from March 2008 to February 2014 in southern Poland, 47 were deemed to have recent HIV-1 infection by the limiting antigen avidity immunoassay. Proviral DNA was amplified and sequenced in the reverse transcriptase, protease, and gp41 coding regions. Mutations were interpreted according to the Stanford Database algorithm and/or the International Antiviral Society USA guidelines. TDR mutations were defined according to the WHO surveillance list.

Results

Among 47 patients with recent HIV-1 infection only 1 (2%) had evidence of TDR mutation. No major resistance mutations were found, but the frequency of strains with ≥1 accessory resistance-associated mutations was high, at 98%. Accessory mutations were present in 11% of reverse transcriptase, 96% of protease, and 27% of gp41 sequences. Mean number of accessory resistance mutations in the reverse transcriptase and protease sequences was higher in viruses with no compensatory mutations in the gp41 HR2 domain than in strains with such mutations (p=0.031).

Conclusions

Despite the low prevalence of strains with TDR mutations, the frequency of accessory mutations was considerable, which may reflect the history of drug pressure among transmitters or natural viral genetic diversity, and may be relevant for future clinical outcomes. The accumulation of the accessory resistance mutations within the pol gene may restrict the occurrence of compensatory mutations related to enfuvirtide resistance or vice versa.

Lesion-Directed Therapies and Monitoring Tumor Evolution Using Liquid Biopsies

Precision oncology relies on targeted drugs, such as kinase inhibitors, that are presently administered based on molecular profiles obtained from surgical or bioptic tissue samples. The inherent ability of human tumors to molecularly evolve in response to drug pressures represents a daunting diagnostic challenge. Circulating free DNA (cfDNA) released from primary and metastatic lesions can be used to draw molecular maps that can be continuously updated to match each tumor's evolution. We will present evidence that liquid biopsies can effectively interrogate how targeted therapies drive lesion-specific drug-resistance mechanisms. The impact of drug-induced molecular heterogeneity on subsequent lines of treatment will also be discussed.

Characterization of the Drug Resistance Profiles of Patients Infected with CRF07_BC Using Phenotypic Assay and Ultra-Deep Pyrosequencing

The usefulness of ultra-deep pyrosequencing (UDPS) for the diagnosis of HIV-1 drug resistance (DR) remains to be determined. Previously, we reported an explosive outbreak of HIV-1 circulating recombinant form (CRF) 07_BC among injection drug users (IDUs) in Taiwan in 2004. The goal of this study was to characterize the DR of CRF07_BC strains using different assays including UDPS. Seven CRF07_BC isolates including 4 from early epidemic (collected in 2004–2005) and 3 from late epidemic (collected in 2008) were obtained from treatment-naïve patient’s peripheral blood mononuclear cells. Viral RNA was extracted directly from patient’s plasma or from cultural supernatant and the pol sequences were determined using RT-PCR sequencing or UDPS. For comparison, phenotypic drug susceptibility assay using MAGIC-5 cells (in-house phenotypic assay) and Antivirogram were performed. In-house phenotypic assay showed that all the early epidemic and none of the late epidemic CRF07_BC isolates were resistant to most protease inhibitors (PIs) (4.4–47.3 fold). Neither genotypic assay nor Antivirogram detected any DR mutations. UDPS showed that early epidemic isolates contained 0.01–0.08% of PI DR major mutations. Furthermore, the combinations of major and accessory PI DR mutations significantly correlated with the phenotypic DR. The in-house phenotypic assay is superior to other conventional phenotypic assays in the detection of DR variants with a frequency as low as 0.01%.

In vivo mutation rates and the landscape of fitness costs of HIV-1

Mutation rates and fitness costs of deleterious mutations are difficult to measure in vivo but essential for a quantitative understanding of evolution. Using whole genome deep sequencing data from longitudinal samples during untreated HIV-1 infection, we estimated mutation rates and fitness costs in HIV-1 from the dynamics of genetic variation. At approximately neutral sites, mutations accumulate with a rate of 1.2 × 10-5 per site per day, in agreement with the rate measured in cell cultures. We estimated the rate from G to A to be the largest, followed by the other transitions C to T, T to C, and A to G, while transversions are less frequent. At other sites, mutations tend to reduce virus replication. We estimated the fitness cost of mutations at every site in the HIV-1 genome using a model of mutation selection balance. About half of all non-synonymous mutations have large fitness costs (>10 percent), while most synonymous mutations have costs <1 percent. The cost of synonymous mutations is especially low in most of pol where we could not detect measurable costs for the majority of synonymous mutations. In contrast, we find high costs for synonymous mutations in important RNA structures and regulatory regions. The intra-patient fitness cost estimates are consistent across multiple patients, indicating that the deleterious part of the fitness landscape is universal and explains a large fraction of global HIV-1 group M diversity.

Strong Selection Significantly Increases Epistatic Interactions in the Long-Term Evolution of a Protein

Epistatic interactions between residues determine a protein’s adaptability and shape its evolutionary trajectory. When a protein experiences a changed environment, it is under strong selection to find a peak in the new fitness landscape. It has been shown that strong selection increases epistatic interactions as well as the ruggedness of the fitness landscape, but little is known about how the epistatic interactions change under selection in the long-term evolution of a protein. Here we analyze the evolution of epistasis in the protease of the human immunodeficiency virus type 1 (HIV-1) using protease sequences collected for almost a decade from both treated and untreated patients, to understand how epistasis changes and how those changes impact the long-term evolvability of a protein. We use an information-theoretic proxy for epistasis that quantifies the co-variation between sites, and show that positive information is a necessary (but not sufficient) condition that detects epistasis in most cases. We analyze the “fossils” of the evolutionary trajectories of the protein contained in the sequence data, and show that epistasis continues to enrich under strong selection, but not for proteins whose environment is unchanged. The increase in epistasis compensates for the information loss due to sequence variability brought about by treatment, and facilitates adaptation in the increasingly rugged fitness landscape of treatment. While epistasis is thought to enhance evolvability via valley-crossing early-on in adaptation, it can hinder adaptation later when the landscape has turned rugged. However, we find no evidence that the HIV-1 protease has reached its potential for evolution after 9 years of adapting to a drug environment that itself is constantly changing. We suggest that the mechanism of encoding new information into pairwise interactions is central to protein evolution not just in HIV-1 protease, but for any protein adapting to a changing environment.

Evolution is often viewed as a process that occurs “mutation by mutation”, suggesting that the effect of each mutation is independent of that of others. However, in reality the effect of a mutation often depends on the context of other mutations, a dependence known as “epistasis”. Even though epistasis can constrain protein evolution, it is actually very common. Such interactions are particularly pervasive in proteins that evolve resistance to a drug via mutations that create defects, and that must be repaired with compensatory mutations. We study how epistasis between protein residues evolves over time in a new and changing environment, and compare these findings to protein evolution in a constant environment. We analyze the sequences of the human immunodeficiency virus type 1 (HIV-1) protease enzyme collected over a period of 9 years from patients treated with anti-viral drugs (as well as from patients that went untreated), and find that epistasis between residues continues to increase as more potent anti-viral drugs enter the market, while epistasis is unchanging in the proteins exposed to a constant environment. Yet, the proteins adapting to the changing landscape do not appear to be constrained by the epistatic interactions and continue to manage to evade new drugs.

Potent Inhibitors Active against HIV Reverse Transcriptase with K101P, a Mutation Conferring Rilpivirine Resistance

Catechol diether compounds have nanomolar antiviral and enzymatic activity against HIV with reverse transcriptase (RT) variants containing K101P, a mutation that confers high-level resistance to FDA-approved non-nucleoside inhibitors efavirenz and rilpivirine. Kinetic data suggests that RT (K101P) variants are as catalytically fit as wild-type and thus can potentially increase in the viral population as more antiviral regimens include efavirenz or rilpivirine. Comparison of wild-type structures and a new crystal structure of RT (K101P) in complex with a leading compound confirms that the K101P mutation is not a liability for the catechol diethers while suggesting that key interactions are lost with efavirenz and rilpivirine.

Longitudinal Detection and Persistence of Minority Drug-Resistant Populations and Their Effect on Salvage Therapy

BACKGROUND

Drug-resistant HIV are more prevalent and persist longer than previously demonstrated by bulk sequencing due to the ability to detect low-frequency variants. To clarify a clinical benefit to monitoring minority-level drug resistance populations as a guide to select active drugs for salvage therapy, we retrospectively analyzed the dynamics of low-frequency drug-resistant population in antiretroviral (ARV)-exposed drug resistant individuals.

MATERIALS AND METHODS

Six HIV-infected individuals treated with ARV for more than five years were analyzed. These individuals had difficulty in controlling viremia, and treatment regimens were switched multiple times guided by standard drug resistance testing using bulk sequencing. To detect minority variant populations with drug resistance, we used a highly sensitive allele-specific PCR (AS-PCR) with detection thresholds of 0.3-2%. According to ARV used in these individuals, we focused on the following seven reverse transcriptase inhibitor-resistant mutations: M41L, K65R, K70R, K103N, Y181C, M184V, and T215F/Y. Results of AS-PCR were compared with bulk sequencing data for concordance and presence of additional mutations. To clarify the genetic relationship between low-frequency and high-frequency populations, AS-PCR amplicon sequences were compared with bulk sequences in phylogenetic analysis.

RESULTS

The use of AS-PCR enabled detection of the drug-resistant mutations, M41L, K103N, Y181C, M184V and T215Y, present as low-frequency populations in five of the six individuals. These drug resistant variants persisted for several years without ARV pressure. Phylogenetic analysis indicated that pre-existing K103N and T215I variants had close genetic relationships with high-frequency K103N and T215I observed during treatment.

DISCUSSION AND CONCLUSION

Our results demonstrate the long-term persistence of drug-resistant viruses in the absence of drug pressure. The rapid virologic failures with pre-existing mutant viruses detectable by AS-PCR highlight the clinical importance of low-frequency drug-resistant viruses. Thus, our results highlight the usefulness of AS-PCR and support its expanded evaluation in ART clinical management.

Genetic Diversity and Drug Resistance Among Antiretroviral Treatment-Failed Individuals from 2010 to 2012 in Honghe, China

The most common antiretroviral treatment (ART) received by individuals infected with HIV-1 in China is the combination therapy, comprised of nucleoside reverse transcriptase inhibitors (NRTIs) and nonnucleoside reverse transcriptase inhibitors (NNRTIs). To assess the prevalence of HIV-1 drug resistance and subtypes in Honghe of Yunnan, China, patient plasmas from ART-failed individuals were collected from January 2010 to December 2012. Genotyping was conducted using an in-house assay on patient plasmas. A total of 254 pol sequences were obtained. The prevalence of drug resistance was 47.2% in ART-failed individuals. Of these drug-resistant individuals, 51.7% harbored HIV strains dually resistant to NRTIs and NNRTIs or protease inhibitors (PIs) (34.2% for NNRTIs and 14.2% for NRTIs). Mutations such as M184V, A62V, T69Ins, K103N, Y181C, and G190A were common among the ART-failed individuals. The frequencies of M184V, A62V, and K103N were 20.5%, 11.0%, and 23.6%, respectively. The most common subtypes in Honghe were CRF08_BC (68.50%) and CRF07_BC (12.20%). The subtypes were almost consistent in different time points for one individual. When receiving ART for 6-12 months, the frequency of HIV-1 drug-resistant variants ranked first. This study shows that the high prevalence of HIV drug resistance observed among the ART-failed individuals should be of increasing concern (monitoring of resistance mutations) in ART regions and facilitate developing novel strategies for prevention and control of HIV infection in China.

The prevalence and determinants of drug-resistance-associated mutations in the HIV-1-infected MSM population of Henan Province in China

To estimate the prevalence of human immunodeficiency virus (HIV) drug resistance (DR) in a population of men who have sex with men (MSM) from Henan Province of China and to identify the DR-associated HIV-1 mutations in these MSM. The HIV-positive status of the MSM subjects in this study was confirmed using ELISA and Western blotting. The MSM subjects were classified into non-treatment group (n = 106) and treatment group (n = 313). CD4(+) T-lymphocyte counts were obtained by flow cytometry, and viral load was measured by branched DNA (bDNA) signal amplification assay. HIV-1 genotypic resistance tests were performed by sequence analysis of the HIV-1 protease and reverse transcriptase genes. In the non-treatment group, 15 patients (14.2 %) displayed DR to non-nucleoside reverse transcriptase inhibitor (NNRTI). In the treatment group, the failure rate of viral suppression was 38.33 % and the DR rate was 33.2 %, which was higher than the rate observed in the non-treatment group (P < 0.05). The incidence of mutations corresponding to NNRTI resistance was significantly higher than the incidence of mutations corresponding to nucleoside reverse transcriptase inhibitor (NRTI) resistance (32.9 % vs. 26.5 %) in the cohort. After antiretroviral therapy (ART), the frequencies of K103N, G190A, Y181C, and V106A mutations were highly elevated. Logistic regression analysis results showed that duration of treatment, poor treatment compliance, drug abuse and homosexual orientation are the major risk factors for DR in this MSM population (all P < 0.05). Our results showed that DR-associated mutations in the HIV-1-infected MSM population increased significantly after ART. Furthermore, duration of treatment, poor treatment compliance, drug abuse and homosexual orientation were identified as the risk factors for DR in the MSM population from Henan Province in China.

Different Effects of Nonnucleoside Reverse Transcriptase Inhibitor Resistance Mutations on Cytotoxic T Lymphocyte Recognition between HIV-1 Subtype B and Subtype A/E Infections

The effect of antiretroviral drug resistance mutations on cytotoxic T lymphocyte (CTL) recognition has been analyzed in HIV-1 subtype B infections, but it remains unclear in infections by other HIV-1 subtypes that are epidemic in countries where antiretroviral drugs are not effectively used. We investigated the effect of nonnucleoside reverse transcriptase (RT) inhibitor (NNRTI)-resistance mutations (Y181C, Y181I, and Y181V) on epitope recognition by CTLs specific for 3 different HIV-1 epitopes (HLA-A*02:01-restricted IV10, HLA-B*35:01-restricted NY9, and HLA-C*12:02-restricted KY9) in subtype B and subtype A/E infections and the accumulation of these mutations in treatment-naive Japanese and Vietnamese. These NNRTI-resistance mutations critically affected NY9-specific and KY9-specific T cell responses in the subtype B infections, whereas they showed a different effect on IV10-specific T cell responses among the subtype B-infected individuals. These mutations affected IV10-specific T cell responses but weakly affected NY9-specific T cell responses in the subtype A/E infections. The substitution at position 3 of NY9 epitope which was found in the subtype A/E virus differently influenced the peptide binding to HLA-B*35:01, suggesting that the differences in peptide binding may result in the differences in T cell recognition between the subtype B virus and A/E virus infections. The Y181C mutation was found to be accumulating in treatment-naive Vietnamese infected with the subtype A/E virus. The present study demonstrated different effects of NNRTI-resistance RT181 mutations on CTL responses between the 2 subtype infections. The Y181C mutation may influence HIV-1 control by the CTLs in Vietnam, since this mutation has been accumulating in treatment-naive Vietnamese.

IMPORTANCE

Antiretroviral therapy leads to the emergence of drug-resistant HIV-1, resulting in virological and clinical failures. Though HIV-1-specific CTLs play a critical role in HIV-1 infection, some of drug resistance mutations located in CTL epitopes are known to affect HIV-1-specific CTL responses. Nonnucleoside reverse transcriptase inhibitor (NNRTI)-resistance RT181 mutations are frequently observed in patients treated with NNRTIs. Such drug resistance mutations may have an influence on immune control by HIV-1-specific CTLs, especially in countries where antiretroviral drugs are not effectively used. We here investigated the effect of three NNRTI-resistance RT181 mutations on immune responses by HIV-1-specific CTLs and the recent accumulation of these mutations in treatment-naive Vietnamese infected with HIV-1 subtype A/E virus. RT181 mutations affected CTL recognition in both subtype A/E and B infections, while the RT Y181C mutation has been accumulating in treatment-naive Vietnamese. The results suggest that the Y181C mutation may influence HIV-1 control by CTLs in Vietnam.

Pairwise growth competition assay for determining the replication fitness of human immunodeficiency viruses

In vitro fitness assays are essential tools for determining viral replication fitness for viruses such as HIV-1. Various measurements have been used to extrapolate viral replication fitness, ranging from the number of viral particles per infectious unit, growth rate in cell culture, and relative fitness derived from multiple-cycle growth competition assays. Growth competition assays provide a particularly sensitive measurement of fitness since the viruses are competing for cellular targets under identical growth conditions. There are several experimental factors to consider when conducting growth competition assays, including the multiplicity of infection (MOI), sampling times, and viral detection and fitness calculation methods. Each factor can affect the end result and hence must be considered carefully during the experimental design. The protocol presented here includes steps from constructing a new recombinant HIV-1 clone to performing growth competition assays and analyzing the experimental results. This protocol utilizes experimental parameter values previously shown to yield consistent and robust results. Alternatives are discussed, as some parameters need to be adjusted according to the cell type and viruses being studied. The protocol contains two alternative viral detection methods to provide flexibility as the availability of instruments, reagents and expertise varies between laboratories.

Persistence of transmitted HIV-1 drug resistance mutations associated with fitness costs and viral genetic backgrounds

Transmission of drug-resistant pathogens presents an almost-universal challenge for fighting infectious diseases. Transmitted drug resistance mutations (TDRM) can persist in the absence of drugs for considerable time. It is generally believed that differential TDRM-persistence is caused, at least partially, by variations in TDRM-fitness-costs. However, in vivo epidemiological evidence for the impact of fitness costs on TDRM-persistence is rare.

Here, we studied the persistence of TDRM in HIV-1 using longitudinally-sampled nucleotide sequences from the Swiss-HIV-Cohort-Study (SHCS). All treatment-naïve individuals with TDRM at baseline were included. Persistence of TDRM was quantified via reversion rates (RR) determined with interval-censored survival models. Fitness costs of TDRM were estimated in the genetic background in which they occurred using a previously published and validated machine-learning algorithm (based on in vitro replicative capacities) and were included in the survival models as explanatory variables.

In 857 sequential samples from 168 treatment-naïve patients, 17 TDRM were analyzed. RR varied substantially and ranged from 174.0/100-person-years;CI=[51.4, 588.8] (for 184V) to 2.7/100-person-years;[0.7, 10.9] (for 215D). RR increased significantly with fitness cost (increase by 1.6[1.3,2.0] per standard deviation of fitness costs). When subdividing fitness costs into the average fitness cost of a given mutation and the deviation from the average fitness cost of a mutation in a given genetic background, we found that both components were significantly associated with reversion-rates.

Our results show that the substantial variations of TDRM persistence in the absence of drugs are associated with fitness-cost differences both among mutations and among different genetic backgrounds for the same mutation.

The evolution of resistance is a universal challenge in antimicrobial chemotherapy. A key driver of resistance is that drug resistance mutations often persist even in the absence of drugs and despite the fact that resistance mutations are often associated with reduced pathogen replication (“fitness costs”). Such persistence may occur because fitness costs are low, especially if they are compensated by additional mutations in their “genetic background”. Here we assessed the role of fitness-cost and the genetic background for resistance in a real-world epidemiological setting by studying the persistence behavior of transmitted antiretroviral resistance mutations of HIV. This persistence behavior was associated with the predicted fitness cost of a given resistance mutation in the particular genetic background in which it occurred. We found that persistence behavior varied strongly across both mutation types and genetic backgrounds and that persistence was significantly associated with predicted fitness costs. In particular we found that even mutations of the same type tended to persist longer if they occurred in a genetic background where they caused weak fitness costs. Overall our results underline the variability of persistence behavior as well as the important role of fitness costs and the genetic background in the evolution of antimicrobial resistance.

Evolution of oncolytic viruses

Owing to their replicative capacity, oncolytic viruses (OVs) can evolve under the action of natural selection. Reversion to virulence and recombination with wild-type strains may compromise OV safety, therefore requiring evolutionary risk assessment studies. On the other hand, evolution can be directed in the laboratory to create more potent and safer OVs. Previous work in the experimental evolution field provides a background for OV directed evolution, and has identified interesting exploitable features. While genetic engineering has greatly advanced the field of oncolytic virotherapy, this approach is sometimes curtailed by the complexity and diversity of virus-host interactions. Directed evolution provides an alternative approach that may help to obtain new OVs without prejudice toward the underlying molecular mechanisms involved.

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.

Low-level Viremia early in HIV infection

HIV RNA levels are usually high early in HIV infection. In the HPTN 061 study, men were tested for HIV infection every 6 months; 6 (21.4%) of 28 men who acquired HIV infection during the study had low or undetectable HIV RNA at the time of HIV diagnosis. Antiretroviral drugs were not detected at the time of HIV diagnosis. False-negative HIV test results were obtained for 2 men using multiple assays. Antiretroviral drug resistance mutations were detected in HIV from 1 man. Additional studies are needed to identify factors associated with low HIV RNA levels during early HIV infection.

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.

In vitro resistance selection with doravirine (MK-1439), a novel nonnucleoside reverse transcriptase inhibitor with distinct mutation development pathways

Doravirine (DOR, formerly known as MK-1439) is a human immunodeficiency type 1 virus (HIV-1) nonnucleoside reverse transcriptase inhibitor (NNRTI) that is currently in phase 2b clinical trials. In vitro resistance selection of subtype B virus (MT4-green fluorescent protein [GFP] cells), as well as subtype A and C viruses (MT4-GFP/CCR5 cells) was conducted with DOR, rilpivirine (RPV), and efavirine (EFV) under low-multiplicity-of-infection conditions in a 96-well format. Resistance selection was performed with escalating concentrations of the NNRTIs ranging from the 95% effective concentration (1 × EC(95)) to 1,000 × EC(95) in the presence of 10% fetal bovine serum. In the resistance selection of subtype B virus with DOR, a V106A mutant virus led to two mutation pathways, followed by the emergence separately of either F227L or L234I. In the resistance selection of subtype A and C viruses, similar mutation development pathways were detected, in which a V106A or V106M mutant was also the starting virus in the pathways. Mutations that are commonly associated with RPV and EFV in clinical settings were also identified in subtype B viruses such as the E138K and K103N mutants, respectively, in this in vitro resistance selection study. The susceptibility of subtype B mutant viruses selected by DOR, RPV, and EFV to NNRTIs was evaluated. Results suggest that mutant viruses selected by DOR are susceptible to RPV and EFV and mutants selected by RPV and EFV are susceptible to DOR. When the replication capacity of the V106A mutant was compared with that of the wild-type (WT) virus, the mutant virus was 4-fold less fit than the WT virus.

Impact of human immunodeficiency virus type-1 sequence diversity on antiretroviral therapy outcomes

Worldwide circulating HIV-1 genomes show extensive variation represented by different subtypes, polymorphisms and drug-resistant strains. Reports on the impact of sequence variation on antiretroviral therapy (ART) outcomes are mixed. In this review, we summarize relevant published data from both resource-rich and resource-limited countries in the last 10 years on the impact of HIV-1 sequence diversity on treatment outcomes. The prevalence of transmission of drug resistant mutations (DRMs) varies considerably, ranging from 0% to 27% worldwide. Factors such as geographic location, access and availability to ART, duration since inception of treatment programs, quality of care, risk-taking behaviors, mode of transmission, and viral subtype all dictate the prevalence in a particular geographical region. Although HIV-1 subtype may not be a good predictor of treatment outcome, review of emerging evidence supports the fact that HIV-1 genome sequence-resulting from natural polymorphisms or drug-associated mutations-matters when it comes to treatment outcomes. Therefore, continued surveillance of drug resistant variants in both treatment-naïve and treatment-experienced populations is needed to reduce the transmission of DRMs and to optimize the efficacy of the current ART armamentarium.

Synthesis and antiviral evaluation of 2-amino-6-carbamoylpurine dioxolane nucleoside derivatives and their phosphoramidates prodrugs

The synthesis of 9-(β-d-1,3-dioxolan-4-yl)2,6-diaminopurine nucleoside phosphoramidate prodrugs as well as various 2-amino-6-carbamoylpurine dioxolane derivatives and their phosphoramidates prodrugs is reported. Their ability to block HIV and HBV replication along with their cytotoxicity toward HepG2, human lymphocyte, CEM and Vero cells was also assessed.

Characterization of amino acids Arg, Ser and Thr at position 70 within HIV-1 reverse transcriptase

OBJECTIVES

The amino acid position 70 in HIV-1 reverse transcriptase (RT) plays an important role in nucleoside RT inhibitor (NRTI) resistance. K70R is part of the thymidine analog mutations, but also other amino acid changes have been associated with NRTI resistance, such as K70E and K70G. In this study, we investigated the in vivo selection of the HIV-1 RT mutations K70S and K70T and their in vitro effect on drug resistance and replication capacity.

METHODS

Recombinant viruses with RT mutations were generated to measure the in vitro drug susceptibility and replication capacity. Bayesian network analysis and three-dimensional modeling were performed to understand the selection and impact of the RT70 mutations.

RESULTS

K70S and K70T were found at a low frequency in RTI-experienced HIV-1 patients (0.10% and 0·20%). Baeyesian network learning identified no direct association with the in vivo exposure to any specific RTI. However, direct associations of K70S with mutations within the Q151M-complex and of K70T with K65R were observed. In vitro phenotypic testing revealed only minor effects of K70R/S/T as single mutations, associated with Q151M and within the context of the Q151M-complex.

DISCUSSION

These results suggest that the selection of K70S/T and their phenotypic impact are influenced by the presence of other mutations in RT. However, the low impact on in vitro phenotype here observed, alongside with the low in vivo prevalence, the exclusive direct association with known major RTI mutations and the unknown correlation with in vivo response, do not yet necessitate the inclusion of K70S/T in drug resistance interpretation systems.

Increased replicative fitness can lead to decreased drug sensitivity of hepatitis C virus

Passage of hepatitis C virus (HCV) in human hepatoma cells resulted in populations that displayed partial resistance to alpha interferon (IFN-α), telaprevir, daclatasvir, cyclosporine, and ribavirin, despite no prior exposure to these drugs. Mutant spectrum analyses and kinetics of virus production in the absence and presence of drugs indicate that resistance is not due to the presence of drug resistance mutations in the mutant spectrum of the initial or passaged populations but to increased replicative fitness acquired during passage. Fitness increases did not alter host factors that lead to shutoff of general host cell protein synthesis and preferential translation of HCV RNA. The results imply that viral replicative fitness is a mechanism of multidrug resistance in HCV. Importance: Viral drug resistance is usually attributed to the presence of amino acid substitutions in the protein targeted by the drug. In the present study with HCV, we show that high viral replicative fitness can confer a general drug resistance phenotype to the virus. The results exclude the possibility that genomes with drug resistance mutations are responsible for the observed phenotype. The fact that replicative fitness can be a determinant of multidrug resistance may explain why the virus is less sensitive to drug treatments in prolonged chronic HCV infections that favor increases in replicative fitness.

Experimental evolution of an oncolytic vesicular stomatitis virus with increased selectivity for p53-deficient cells

Experimental evolution has been used for various biotechnological applications including protein and microbial cell engineering, but less commonly in the field of oncolytic virotherapy. Here, we sought to adapt a rapidly evolving RNA virus to cells deficient for the tumor suppressor gene p53, a hallmark of cancer cells. To achieve this goal, we established four independent evolution lines of the vesicular stomatitis virus (VSV) in p53-knockout mouse embryonic fibroblasts (p53-/- MEFs) under conditions favoring the action of natural selection. We found that some evolved viruses showed increased fitness and cytotoxicity in p53-/- cells but not in isogenic p53+/+ cells, indicating gene-specific adaptation. However, full-length sequencing revealed no obvious or previously described genetic changes associated with oncolytic activity. Half-maximal effective dose (EC50) assays in mouse p53-positive colon cancer (CT26) and p53-deficient breast cancer (4T1) cells indicated that the evolved viruses were more effective against 4T1 cells than the parental virus or a reference oncolytic VSV (MΔ51), but showed no increased efficacy against CT26 cells. In vivo assays using 4T1 syngeneic tumor models showed that one of the evolved lines significantly delayed tumor growth compared to mice treated with the parental virus or untreated controls, and was able to induce transient tumor suppression. Our results show that RNA viruses can be specifically adapted typical cancer features such as p53 inactivation, and illustrate the usefulness of experimental evolution for oncolytic virotherapy.

HIV-1 quasispecies delineation by tag linkage deep sequencing

Trade-offs between throughput, read length, and error rates in high-throughput sequencing limit certain applications such as monitoring viral quasispecies. Here, we describe a molecular-based tag linkage method that allows assemblage of short sequence reads into long DNA fragments. It enables haplotype phasing with high accuracy and sensitivity to interrogate individual viral sequences in a quasispecies. This approach is demonstrated to deduce ∼2000 unique 1.3 kb viral sequences from HIV-1 quasispecies in vivo and after passaging ex vivo with a detection limit of ∼0.005% to ∼0.001%. Reproducibility of the method is validated quantitatively and qualitatively by a technical replicate. This approach can improve monitoring of the genetic architecture and evolution dynamics in any quasispecies population.

The contribution of viral genotype to plasma viral set-point in HIV infection

Disease progression in HIV-infected individuals varies greatly, and while the environmental and host factors influencing this variation have been widely investigated, the viral contribution to variation in set-point viral load, a predictor of disease progression, is less clear. Previous studies, using transmission-pairs and analysis of phylogenetic signal in small numbers of individuals, have produced a wide range of viral genetic effect estimates. Here we present a novel application of a population-scale method based in quantitative genetics to estimate the viral genetic effect on set-point viral load in the UK subtype B HIV-1 epidemic, based on a very large data set. Analyzing the initial viral load and associated pol sequence, both taken before anti-retroviral therapy, of 8,483 patients, we estimate the proportion of variance in viral load explained by viral genetic effects to be 5.7% (CI 2.8-8.6%). We also estimated the change in viral load over time due to selection on the virus and environmental effects to be a decline of 0.05 log10 copies/mL/year, in contrast to recent studies which suggested a reported small increase in viral load over the last 20 years might be due to evolutionary changes in the virus. Our results suggest that in the UK epidemic, subtype B has a small but significant viral genetic effect on viral load. By allowing the analysis of large sample sizes, we expect our approach to be applicable to the estimation of the genetic contribution to traits in many organisms.

Viral factors in non-progression

Research has undergone considerable development in understanding a small subset of human immunodeficiency virus type 1 (HIV-1)-infected, therapy-naive individuals who maintain a favorable course of infection surviving for longer periods of time. Although, viral, host genetic, and immunological factors have been analyzed in many previous studies in order to delineate mechanisms that contribute to non-progressive HIV disease, there appears to be a no clear cut winner and the non-progressive HIV disease in <1% of HIV-infected individuals appears to be a complex interplay between viral and host factors. Therefore, it is important to review them separately to signify their potential contribution to non-progressive HIV disease. With respect to virological features, genomic sequencing of HIV-1 strains derived from long-term non-progressors has shown that some individuals are infected with attenuated strains of HIV-1 and harbor mutations from single nucleotide polymorphisms to large deletions in HIV-1 structure, regulatory, and accessory genes. The elucidation of functional attributes of defective/attenuated HIV strains may provide better understanding of viral pathogenesis and the discovery of new therapeutic strategies against HIV. This review mainly focuses on the defects in viral genes that possibly guide non-progressive HIV disease.

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.