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

Contingency and Entrenchment of Drug-Resistance Mutations in HIV Viral Proteins

The ability of HIV-1 to rapidly mutate leads to antiretroviral therapy (ART) failure among infected patients. Drug-resistance mutations (DRMs), which cause a fitness penalty to intrinsic viral fitness, are compensated by accessory mutations with favorable epistatic interactions which cause an evolutionary trapping effect, but the kinetics of this overall process has not been well characterized. Here, using a Potts Hamiltonian model describing epistasis combined with kinetic Monte Carlo simulations of evolutionary trajectories, we explore how epistasis modulates the evolutionary dynamics of HIV DRMs. We show how the occurrence of a drug-resistance mutation is contingent on favorable epistatic interactions with many other residues of the sequence background and that subsequent mutations entrench DRMs. We measure the time-autocorrelation of fluctuations in the likelihood of DRMs due to epistatic coupling with the sequence background, which reveals the presence of two evolutionary processes controlling DRM kinetics with two distinct time scales. Further analysis of waiting times for the evolutionary trapping effect to reverse reveals that the sequences which entrench (trap) a DRM are responsible for the slower time scale. We also quantify the overall strength of epistatic effects on the evolutionary kinetics for different mutations and show these are much larger for DRM positions than polymorphic positions, and we also show that trapping of a DRM is often caused by the collective effect of many accessory mutations, rather than a few strongly coupled ones, suggesting the importance of multiresidue sequence variations in HIV evolution. The analysis presented here provides a framework to explore the kinetic pathways through which viral proteins like HIV evolve under drug-selection pressure.

The clarifying role of time series data in the population genetics of HIV

HIV can evolve remarkably quickly in response to antiretroviral therapies and the immune system. This evolution stymies treatment effectiveness and prevents the development of an HIV vaccine. Consequently, there has been a great interest in using population genetics to disentangle the forces that govern the HIV adaptive landscape (selection, drift, mutation, and recombination). Traditional population genetics approaches look at the current state of genetic variation and infer the processes that can generate it. However, because HIV evolves rapidly, we can also sample populations repeatedly over time and watch evolution in action. In this paper, we demonstrate how time series data can bound evolutionary parameters in a way that complements and informs traditional population genetic approaches. Specifically, we focus on our recent paper (Feder et al., 2016, eLife), in which we show that, as improved HIV drugs have led to fewer patients failing therapy due to resistance evolution, less genetic diversity has been maintained following the fixation of drug resistance mutations. Because soft sweeps of multiple drug resistance mutations spreading simultaneously have been previously documented in response to the less effective HIV therapies used early in the epidemic, we interpret the maintenance of post-sweep diversity in response to poor therapies as further evidence of soft sweeps and therefore a high population mutation rate (θ) in these intra-patient HIV populations. Because improved drugs resulted in rarer resistance evolution accompanied by lower post-sweep diversity, we suggest that both observations can be explained by decreased population mutation rates and a resultant transition to hard selective sweeps. A recent paper (Harris et al., 2018, PLOS Genetics) proposed an alternative interpretation: Diversity maintenance following drug resistance evolution in response to poor therapies may have been driven by recombination during slow, hard selective sweeps of single mutations. Then, if better drugs have led to faster hard selective sweeps of resistance, recombination will have less time to rescue diversity during the sweep, recapitulating the decrease in post-sweep diversity as drugs have improved. In this paper, we use time series data to show that drug resistance evolution during ineffective treatment is very fast, providing new evidence that soft sweeps drove early HIV treatment failure.

Drug Resistance Evolution in HIV in the Late 1990s: Hard Sweeps, Soft Sweeps, Clonal Interference and the Accumulation of Drug Resistance Mutations

The evolution of drug resistance in pathogens such as HIV is an important and widely known example in the field of evolutionary medicine. Here, we focus on a unique data set from the late 1990s with multiple viral sequences from multiple time points in 118 patients. We study patterns of evolutionary dynamics in the viral populations in these patients who were treated with Reverse Transcriptase Inhibitors and Protease Inhibitors in the late 1990s. Specifically, we aim to visualize and analyze examples of population genetic processes such as selective sweeps and clonal interference. The figures and descriptions in this paper can be used in evolution and population genetics classes. We show and analyze a wide variety of patterns, specifically: soft sweeps, hard sweeps, softening sweeps and hardening sweeps, simultaneous sweeps, accumulation of mutations and clonal interference.

HIV evolution and diversity in ART-treated patients

Characterizing HIV genetic diversity and evolution during antiretroviral therapy (ART) provides insights into the mechanisms that maintain the viral reservoir during ART. This review describes common methods used to obtain and analyze intra-patient HIV sequence data, the accumulation of diversity prior to ART and how it is affected by suppressive ART, the debate on viral replication and evolution in the presence of ART, HIV compartmentalization across various tissues, and mechanisms for the emergence of drug resistance. It also describes how CD4+ T cells that were likely infected with latent proviruses prior to initiating treatment can proliferate before and during ART, providing a renewable source of infected cells despite therapy. Some expanded cell clones carry intact and replication-competent proviruses with a small fraction of the clonal siblings being transcriptionally active and a source for residual viremia on ART. Such cells may also be the source for viral rebound after interrupting ART. The identical viral sequences observed for many years in both the plasma and infected cells of patients on long-term ART are likely due to the proliferation of infected cells both prior to and during treatment. Studies on HIV diversity may reveal targets that can be exploited in efforts to eradicate or control the infection without ART.

Characterization of HIV-1 Near Full-Length Proviral Genome Quasispecies from Patients with Undetectable Viral Load Undergoing First-Line HAART Therapy

Increased access to highly active antiretroviral therapy (HAART) by human immunodeficiency virus postive (HIV⁺) individuals has become a reality worldwide. In Brazil, HAART currently reaches over half of HIV-infected subjects. In the context of a remarkable HIV-1 genetic variability, highly related variants, called quasispecies, are generated. HIV quasispecies generated during infection can influence virus persistence and pathogenicity, representing a challenge to treatment. However, the clinical relevance of minority quasispecies is still uncertain. In this study, we have determined the archived proviral sequences, viral subtype and drug resistance mutations from a cohort of HIV⁺ patients with undetectable viral load undergoing HAART as first-line therapy using next-generation sequencing for near full-length virus genome (NFLG) assembly. HIV-1 consensus sequences representing NFLG were obtained for eleven patients, while for another twelve varying genome coverage rates were obtained. Phylogenetic analysis showed the predominance of subtype B (83%; 19/23). Considering the minority variants, 18 patients carried archived virus harboring at least one mutation conferring antiretroviral resistance; for six patients, the mutations correlated with the current ARVs used. These data highlight the importance of monitoring HIV minority drug resistant variants and their clinical impact, to guide future regimen switches and improve HIV treatment success.

Low Frequency of Drug-Resistant Variants Selected by Long-Acting Rilpivirine in Macaques Infected with Simian Immunodeficiency Virus Containing HIV-1 Reverse Transcriptase

Preexposure prophylaxis (PrEP) using antiretroviral drugs is effective in reducing the risk of human immunodeficiency virus type 1 (HIV-1) infection, but adherence to the PrEP regimen is needed. To improve adherence, a long-acting injectable formulation of the nonnucleoside reverse transcriptase (RT) inhibitor rilpivirine (RPV LA) has been developed. However, there are concerns that PrEP may select for drug-resistant mutations during preexisting or breakthrough infections, which could promote the spread of drug resistance and limit options for antiretroviral therapy. To address this concern, we administered RPV LA to macaques infected with simian immunodeficiency virus containing HIV-1 RT (RT-SHIV). Peak plasma RPV levels were equivalent to those reported in human trials and waned over time after dosing. RPV LA resulted in a 2-log decrease in plasma viremia, and the therapeutic effect was maintained for 15 weeks, until plasma drug concentrations dropped below 25 ng/ml. RT mutations E138G and E138Q were detected in single clones from plasma virus in separate animals only at one time point, and no resistance mutations were detected in viral RNA isolated from tissues. Wild-type and E138Q RT-SHIV displayed similar RPV susceptibilities in vitro, whereas E138G conferred 2-fold resistance to RPV. Overall, selection of RPV-resistant variants was rare in an RT-SHIV macaque model despite prolonged exposure to slowly decreasing RPV concentrations following injection of RPV LA.

Effect of suberoylanilide hydroxamic acid (SAHA) administration on the residual virus pool in a model of combination antiretroviral therapy-mediated suppression in SIVmac239-infected indian rhesus macaques

Nonhuman primate models are needed for evaluations of proposed strategies targeting residual virus that persists in HIV-1-infected individuals receiving suppressive combination antiretroviral therapy (cART). However, relevant nonhuman primate (NHP) models of cART-mediated suppression have proven challenging to develop. We used a novel three-class, six-drug cART regimen to achieve durable 4.0- to 5.5-log reductions in plasma viremia levels and declines in cell-associated viral RNA and DNA in blood and tissues of simian immunodeficiency virus SIVmac239-infected Indian-origin rhesus macaques, then evaluated the impact of treatment with the histone deacetylase inhibitor (HDACi) suberoylanilide hydroxamic acid (SAHA; Vorinostat) on the residual virus pool. Ex vivo SAHA treatment of CD4(+) T cells obtained from cART-suppressed animals increased histone acetylation and viral RNA levels in culture supernatants. cART-suppressed animals each received 84 total doses of oral SAHA. We observed SAHA dose-dependent increases in acetylated histones with evidence for sustained modulation as well as refractoriness following prolonged administration. In vivo virologic activity was demonstrated based on the ratio of viral RNA to viral DNA in peripheral blood mononuclear cells, a presumptive measure of viral transcription, which significantly increased in SAHA-treated animals. However, residual virus was readily detected at the end of treatment, suggesting that SAHA alone may be insufficient for viral eradication in the setting of suppressive cART. The effects observed were similar to emerging data for repeat-dose SAHA treatment of HIV-infected individuals on cART, demonstrating the feasibility, utility, and relevance of NHP models of cART-mediated suppression for in vivo assessments of AIDS virus functional cure/eradication approaches.

Development and customization of a color-coded microbeads-based assay for drug resistance in HIV-1 reverse transcriptase

Background

Drug resistance (DR) of HIV-1 can be examined genotypically or phenotypically. Although sequencing is the gold standard of the genotypic resistance testing (GRT), high-throughput GRT targeted to the codons responsible for DR may be more appropriate for epidemiological studies and public health research.

Methods

We used a Japanese database to design and synthesize sequence-specific oligonucleotide probes (SSOP) for the detection of wild-type sequences and 6 DR mutations in the clade B HIV-1 reverse transcriptase region. We coupled SSOP to microbeads of the Luminex 100 xMAP system and developed a GRT based on the polymerase chain reaction (PCR)-SSOP-Luminex method.

Results

Sixteen oligoprobes for discriminating DR mutations from wild-type sequences at 6 loci were designed and synthesized, and their sensitivity and specificity were confirmed using isogenic plasmids. The PCR-SSOP-Luminex DR assay was then compared to direct sequencing using 74 plasma specimens from treatment-naïve patients or those on failing treatment. In the majority of specimens, the results of the PCR-SSOP-Luminex DR assay were concordant with sequencing results: 62/74 (83.8%) for M41, 43/74 (58.1%) for K65, 70/74 (94.6%) for K70, 55/73 (75.3%) for K103, 63/73 (86.3%) for M184 and 68/73 (93.2%) for T215. There were a number of specimens without any positive signals, especially for K65. The nucleotide position of A2723G, A2747G and C2750T were frequent polymorphisms for the wild-type amino acids K65, K66 and D67, respectively, and 14 specimens had the D67N mutation encoded by G2748A. We synthesized 14 additional oligoprobes for K65, and the sensitivity for K65 loci improved from 43/74 (58.1%) to 68/74 (91.9%).

Conclusions

We developed a rapid high-throughput assay for clade B HIV-1 DR mutations, which could be customized by synthesizing oligoprobes suitable for the circulating viruses. The assay could be a useful tool especially for public health research in both resource-rich and resource-limited settings.

Robust suppression of env-SHIV viremia in Macaca nemestrina by 3-drug ART is independent of timing of initiation during chronic infection

BACKGROUND

Nonhuman primates (NHPs) are an important model organism for studies of HIV pathogenesis and preclinical evaluation of anti-HIV therapies. The successful translation of NHP-derived data to clinically relevant anti-HIV studies will require better understanding of the viral strains and NHP species used and their responses to existing antiretroviral therapies (ART).

METHODS

Five pigtailed macaques (Macaca nemestrina) were productively infected with the SIV/HIV chimeric virus SHIV-1157 ipd3N4 following intravenous challenge. After 8 or 27 weeks, ART (PMPA, FTC, raltegravir) was initiated. Viral load, T-cell counts, and production of SHIV-specific antibodies were monitored throughout the course of infection and ART.

RESULTS

ART led to a rapid and sustained decrease in plasma viral load. Suppression of plasma viremia by ART was independent of the timing of initiation during chronic infection.

CONCLUSIONS

We present a new NHP model of HIV infection on antiretroviral therapy, which should prove applicable to multiple clinically relevant anti-HIV approaches.

Exposure to MIV-150 from a high-dose intravaginal ring results in limited emergence of drug resistance mutations in SHIV-RT infected rhesus macaques

When microbicides used for HIV prevention contain antiretroviral drugs, there is concern for the potential emergence of drug-resistant HIV following use in infected individuals who are either unaware of their HIV infection status or who are aware but still choose to use the microbicide. Resistant virus could ultimately impact their responsiveness to treatment and/or result in subsequent transmission of drug-resistant virus. We tested whether drug resistance mutations (DRMs) would emerge in macaques infected with simian immunodeficiency virus expressing HIV reverse transcriptase (SHIV-RT) after sustained exposure to the potent non-nucleoside reverse transcriptase inhibitor (NNRTI) MIV-150 delivered via an intravaginal ring (IVR). We first treated 4 SHIV-RT-infected animals with daily intramuscular injections of MIV-150 over two 21 day (d) intervals separated by a 7 d drug hiatus. In all 4 animals, NNRTI DRMs (single and combinations) were detected within 14 d and expanded in proportion and diversity with time. Knowing that we could detect in vivo emergence of NNRTI DRMs in response to MIV-150, we then tested whether a high-dose MIV-150 IVR (loaded with >10 times the amount being used in a combination microbicide IVR in development) would select for resistance in 6 infected animals, modeling use of this prevention method by an HIV-infected woman. We previously demonstrated that this MIV-150 IVR provides significant protection against vaginal SHIV-RT challenge. Wearing the MIV-150 IVR for 56 d led to only 2 single DRMs in 2 of 6 animals (430 RT sequences analyzed total, 0.46%) from plasma and lymph nodes despite MIV-150 persisting in the plasma, vaginal fluids, and genital tissues. Only wild type virus sequences were detected in the genital tissues. These findings indicate a low probability for the emergence of DRMs after topical MIV-150 exposure and support the advancement of MIV-150-containing microbicides.

Loss and recovery of genetic diversity in adapting populations of HIV

The evolution of drug resistance in HIV occurs by the fixation of specific, well-known, drug-resistance mutations, but the underlying population genetic processes are not well understood. By analyzing within-patient longitudinal sequence data, we make four observations that shed a light on the underlying processes and allow us to infer the short-term effective population size of the viral population in a patient. Our first observation is that the evolution of drug resistance usually occurs by the fixation of one drug-resistance mutation at a time, as opposed to several changes simultaneously. Second, we find that these fixation events are accompanied by a reduction in genetic diversity in the region surrounding the fixed drug-resistance mutation, due to the hitchhiking effect. Third, we observe that the fixation of drug-resistance mutations involves both hard and soft selective sweeps. In a hard sweep, a resistance mutation arises in a single viral particle and drives all linked mutations with it when it spreads in the viral population, which dramatically reduces genetic diversity. On the other hand, in a soft sweep, a resistance mutation occurs multiple times on different genetic backgrounds, and the reduction of diversity is weak. Using the frequency of occurrence of hard and soft sweeps we estimate the effective population size of HIV to be 1.5 x 10(5) (95% confidence interval [0.8 x 10(5),4.8 x 10(5)]). This number is much lower than the actual number of infected cells, but much larger than previous population size estimates based on synonymous diversity. We propose several explanations for the observed discrepancies. Finally, our fourth observation is that genetic diversity at non-synonymous sites recovers to its pre-fixation value within 18 months, whereas diversity at synonymous sites remains depressed after this time period. These results improve our understanding of HIV evolution and have potential implications for treatment strategies.

Low-frequency nevirapine (NVP)-resistant HIV-1 variants are not associated with failure of antiretroviral therapy in women without prior exposure to single-dose NVP

BACKGROUND

Low-frequency nevirapine (NVP)-resistant variants have been associated with virologic failure (VF) of initial NVP-based combination antiretroviral therapy (cART) in women with prior exposure to single-dose NVP (sdNVP). We investigated whether a similar association exists in women without prior sdNVP exposure.

METHODS

Pre-cART plasma was analyzed by allele-specific polymerase chain reaction to quantify NVP-resistant mutants in human immunodeficiency virus-infected African women without prior sdNVP who were starting first-line NVP-based cART in the OCTANE/A5208 trial 2. Associations between NVP-resistant mutants and VF or death were determined and compared with published results from women participating in the OCTANE/A5208 trial 1 who had taken sdNVP and initiated NVP-based cART.

RESULTS

Pre-cART NVP-resistant variants were detected in 18% (39/219) of women without prior sdNVP exposure, compared to 45% (51/114) with prior sdNVP exposure (P < .001). Among women without prior sdNVP exposure, 8 of 39 (21%) with NVP-resistant variants experienced VF or death vs 31 of 180 (17%) without such variants (P = .65); this compares with 21 of 51 (41%) vs 9 of 63 (14%) among women with prior exposure (P = .001).

CONCLUSIONS

The risk of VF on NVP-based cART from NVP-resistant variants differs between sdNVP-exposed and -unexposed women. This difference may be driven by drug-resistance mutations emerging after sdNVP exposure that are linked on the same viral genome.

CLINICAL TRIALS REGISTRATION

NCT00089505.

Inferring HIV Escape Rates from Multi-Locus Genotype Data

Cytotoxic T-lymphocytes (CTLs) recognize viral protein fragments displayed by major histocompatibility complex molecules on the surface of virally infected cells and generate an anti-viral response that can kill the infected cells. Virus variants whose protein fragments are not efficiently presented on infected cells or whose fragments are presented but not recognized by CTLs therefore have a competitive advantage and spread rapidly through the population. We present a method that allows a more robust estimation of these escape rates from serially sampled sequence data. The proposed method accounts for competition between multiple escapes by explicitly modeling the accumulation of escape mutations and the stochastic effects of rare multiple mutants. Applying our method to serially sampled HIV sequence data, we estimate rates of HIV escape that are substantially larger than those previously reported. The method can be extended to complex escapes that require compensatory mutations. We expect our method to be applicable in other contexts such as cancer evolution where time series data is also available.

Quantifying selection against synonymous mutations in HIV-1 env evolution

Intrapatient evolution of human immunodeficiency virus type 1 (HIV-1) is driven by the adaptive immune system resulting in rapid change of HIV-1 proteins. When cytotoxic CD8(+) T cells or neutralizing antibodies target a new epitope, the virus often escapes via nonsynonymous mutations that impair recognition. Synonymous mutations do not affect this interplay and are often assumed to be neutral. We test this assumption by tracking synonymous mutations in longitudinal intrapatient data from the C2-V5 part of the env gene. We find that most synonymous variants are lost even though they often reach high frequencies in the viral population, suggesting a cost to the virus. Using published data from SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) assays, we find that synonymous mutations that disrupt base pairs in RNA stems flanking the variable loops of gp120 are more likely to be lost than other synonymous changes: these RNA hairpins might be important for HIV-1. Computational modeling indicates that, to be consistent with the data, a large fraction of synonymous mutations in this genomic region need to be deleterious with a cost on the order of 0.002 per day. This weak selection against synonymous substitutions does not result in a strong pattern of conservation in cross-sectional data but slows down the rate of evolution considerably. Our findings are consistent with the notion that large-scale patterns of RNA structure are functionally relevant, whereas the precise base pairing pattern is not.

Simultaneous detection of major drug resistance mutations in the protease and reverse transcriptase genes for HIV-1 subtype C by use of a multiplex allele-specific assay

High-throughput, sensitive, and cost-effective HIV drug resistance (HIVDR) detection assays are needed for large-scale monitoring of the emergence and transmission of HIVDR in resource-limited settings. Using suspension array technology, we have developed a multiplex allele-specific (MAS) assay that can simultaneously detect major HIVDR mutations at 20 loci. Forty-five allele-specific primers tagged with unique 24-base oligonucleotides at the 5' end were designed to detect wild-type and mutant alleles at the 20 loci of HIV-1 subtype C. The MAS assay was first established and optimized with three plasmid templates (C-wt, C-mut1, and C-mut2) and then evaluated using 148 plasma specimens from HIV-1 subtype C-infected individuals. All the wild-type and mutant alleles were unequivocally distinguished with plasmid templates, and the limits of detection were 1.56% for K219Q and K219E, 3.13% for L76V, 6.25% for K65R, K70R, L74V, L100I, K103N, K103R, Q151M, Y181C, and I47V, and 12.5% for M41L, K101P, K101E, V106A, V106M, Y115F, M184V, Y188L, G190A, V32I, I47A, I84V, and L90M. Analyses of 148 plasma specimens revealed that the MAS assay gave 100% concordance with conventional sequencing at eight loci and >95% (range, 95.21% to 99.32%) concordance at the remaining 12 loci. The differences observed were caused mainly by 24 additional low-abundance alleles detected by the MAS assay. Ultradeep sequencing analysis confirmed 15 of the 16 low-abundance alleles. This multiplex, sensitive, and straightforward result-reporting assay represents a new efficient genotyping tool for HIVDR surveillance and monitoring.

Considerations in the development of nonhuman primate models of combination antiretroviral therapy for studies of AIDS virus suppression, residual virus, and curative strategies

PURPOSE OF REVIEW

Animal models will be critical for preclinical evaluations of novel HIV eradication and/or functional cure strategies in the setting of suppressive combination antiretroviral therapy (cART). Here, the strengths, limitations, and challenges of recent efforts to develop nonhuman primate (NHP) models of cART-mediated suppression for use in studies of persistent virus and curative approaches are discussed.

RECENT FINDINGS

Several combinations of NHP species and viruses that recapitulate key aspects of human HIV infection have been adapted for cART-mediated suppression studies. Different cART regimens incorporating drugs targeting multiple different steps of the viral replication cycle have provided varying levels of virologic suppression, dependent in part upon the host species, virus, drug regimen and timing, and virologic monitoring assay sensitivity. New, increasingly sensitive virologic monitoring approaches for measurements of plasma viral RNA, cell-associated and tissue-associated viral RNA and DNA, and the replication-competent residual viral pool in the setting of cART in NHP models are being developed to allow for the assessment of persistent virus on cART and to evaluate the impact of viral induction/eradication strategies in vivo.

SUMMARY

Given the vagaries of each specific virus and host species, and cART regimen, each model will require further development and analysis to determine their appropriate application for addressing specific experimental questions.