Analysis of low-frequency mutations associated with drug resistance to raltegravir before antiretroviral treatment

Primary resistance to integrase strand transfer inhibitors in Spain using ultrasensitive HIV-1 genotyping

BACKGROUND

Transmission of resistance mutations to integrase strand transfer inhibitors (INSTIs) in HIV-infected patients may compromise the efficacy of first-line antiretroviral regimens currently recommended worldwide. Continued surveillance of transmitted drug resistance (TDR) is thus warranted.

OBJECTIVES

We evaluated the rates and effects on virological outcomes of TDR in a 96 week prospective multicentre cohort study of ART-naive HIV-1-infected subjects initiating INSTI-based ART in Spain between April 2015 and December 2016.

METHODS

Pre-ART plasma samples were genotyped for integrase, protease and reverse transcriptase resistance using Sanger population sequencing or MiSeq™ using a ≥ 20% mutant sensitivity cut-off. Those present at 1%-19% of the virus population were considered to be low-frequency variants.

RESULTS

From a total of 214 available samples, 173 (80.8%), 210 (98.1%) and 214 (100.0%) were successfully amplified for integrase, reverse transcriptase and protease genes, respectively. Using a Sanger-like cut-off, the overall prevalence of any TDR, INSTI-, NRTI-, NNRTI- and protease inhibitor (PI)-associated mutations was 13.1%, 1.7%, 3.8%, 7.1% and 0.9%, respectively. Only three (1.7%) subjects had INSTI TDR (R263K, E138K and G163R), while minority variants with integrase TDR were detected in 9.6% of subjects. There were no virological failures during 96 weeks of follow-up in subjects harbouring TDR as majority variants.

CONCLUSIONS

Transmitted INSTI resistance remains rare in Spain and, to date, is not associated with virological failure to first-line INSTI-based regimens.

Recent Advances in the Development of Integrase Inhibitors for HIV Treatment

PURPOSE OF THE REVIEW

The complex multistep life cycle of HIV allows it to proliferate within the host and integrate its genome in to the host chromosomal DNA. This provirus can remain dormant for an indefinite period. The process of integration, governed by integrase (IN), is highly conserved across the Retroviridae family. Hence, targeting integration is not only expected to block HIV replication but may also reveal new therapeutic strategies to treat HIV as well as other retrovirus infections.

RECENT FINDINGS

HIV integrase (IN) has gained attention as the most promising therapeutic target as there are no equivalent homologues of IN that has been discovered in humans. Although current nano-formulated long-acting IN inhibitors have demonstrated the phenomenal ability to block HIV integration and replication with extraordinary half-life, they also have certain limitations. In this review, we have summarized the current literature on clinically established IN inhibitors, their mechanism of action, the advantages and disadvantages associated with their therapeutic application, and finally current HIV cure strategies using these inhibitors.

Prevalence and clinical impact of minority resistant variants in patients failing an integrase inhibitor-based regimen by ultra-deep sequencing

Background

Integrase strand transfer inhibitors (INSTIs) are recommended by international guidelines as first-line therapy in antiretroviral-naive and -experienced HIV-1-infected patients.

Objectives

This study aimed at evaluating the prevalence at failure of INSTI-resistant variants and the impact of baseline minority resistant variants (MiRVs) on the virological response to an INSTI-based regimen.

Methods

Samples at failure of 134 patients failing a raltegravir-containing (n = 65), an elvitegravir-containing (n = 20) or a dolutegravir-containing (n = 49) regimen were sequenced by Sanger sequencing and ultra-deep sequencing (UDS). Baseline samples of patients with virological failure (VF) (n = 34) and of those with virological success (VS) (n = 31) under INSTI treatment were sequenced by UDS. Data were analysed using the SmartGene platform, and resistance was interpreted according to the ANRS algorithm version 27.

Results

At failure, the prevalence of at least one INSTI-resistant variant was 39.6% by Sanger sequencing and 57.5% by UDS, changing the interpretation of resistance in 17/134 (13%) patients. Among 53 patients harbouring at least one resistance mutation detected by both techniques, the most dominant INSTI resistance mutations were N155H (45%), Q148H/K/R (23%), T97A (19%) and Y143C (11%). There was no difference in prevalence of baseline MiRVs between patients with VF and those with VS. MiRVs found at baseline in patients with VF were not detected at failure either in majority or minority mutations.

Conclusions

UDS is more sensitive than Sanger sequencing at detecting INSTI MiRVs at treatment failure. The presence of MiRVs at failure could be important to the decision to switch to other INSTIs. However, there was no association between the presence of baseline MiRVs and the response to INSTI-based therapies in our study.

Analysis of transmitted HIV drug resistance from 2005 to 2015 in Victoria, Australia: a comparison of the old and the new

Background Baseline genotyping is part of standard-of-care treatment. It reveals that transmitted drug resistance (TDR) continues to be important for the management of HIV infection. Attention is typically focused on determining whether resistance to the protease inhibitors (PI) and reverse transcriptase inhibitors (RTI) occurs. However, the increasing use of integrase inhibitors (INIs) raises a concern that TDR to this class of antiretroviral drug may also occur.

METHODS

PI and RTI drug resistance genotyping was performed on blood samples collected between 2005 and 2015 from 772 treatment-naïve Victorian patients infected with HIV within the previous 12 months. Integrase genotyping was performed on 461 of the 485 patient samples collected between 2010 and 2015.

RESULTS

In the period 2005-10, 39 of 343 patients (11.4%) had at least one PI- or RTI-associated mutation, compared with 34 of 429 (7.9%) during the period 2011-15. Compared with 2005-10, during 2011-15 there was a significant decline in the prevalence of the non-nucleoside-associated mutation K103N and the nucleoside-associated mutations at codons M41 and T215. One patient was detected with a major INI resistance mutation, namely G118R. However, this mutation is rare and its effect on susceptibility is unclear. A small number of patients (n=12) was infected with HIV containing accessory resistance mutations in the integrase gene.

CONCLUSIONS

The lack of transmitted resistance to INIs is consistent with a low level of resistance to this class of drugs in the treated population. However, continued surveillance in the newly infected population is warranted as the use of INIs increases.

Primary resistance to integrase strand transfer inhibitors in patients infected with diverse HIV-1 subtypes in sub-Saharan Africa

Objectives

To investigate the prevalence and patterns of major and accessory resistance mutations associated with integrase strand transfer inhibitors (INSTIs), across diverse HIV-1 subtypes in sub-Saharan Africa.

Methods

pol gene sequences were obtained using Illumina next-generation sequencing from 425 INSTI-naive HIV-infected adults from Kenya (21.2%), Nigeria (7.3%), South Africa (22.8%), Uganda (25.2%) and Zambia (23.5%). Drug resistance interpretation was based on the IAS 2017 mutation list and accessory mutations from Stanford HIVdb with resistance penalty scores of ≥10 to at least 1 INSTI. Resistance was further classified based on sensitivity thresholds of ≥20% (Sanger sequencing) and 1%-20% for low-frequency variants (next-generation sequencing).

Results

Of 425 genotypes, 48.7% were subtype C, 28.5% A, 10.1% D, 2.8% G and 9.9% were recombinants. Major INSTI resistance mutations were detected only at <20% threshold, at a prevalence of 2.4% (2.5% in subtype A, 2.4% C, 0% D, 8.3% G and 2.4% in recombinants) and included T66A/I (0.7%), E92G (0.5%), Y143C/S (0.7%), S147G (0.2%) and Q148R (0.5%). Accessory mutations occurred at a prevalence of 15.1% at the ≥20% threshold (23.1% in subtype A, 8.7% C, 11.6% D, 25% G and 23.8% in recombinants), and included L74I/M (10.4%), Q95K (0.5%), T97A (4%), E157Q (0.7%) and G163R/K (0.7%).

Conclusions

Major INSTI resistance mutations were rare and only occurred at low-level resistance detection thresholds. INSTI-based regimens are expected to be effective across the different major HIV-1 subtypes in the region.

Sanger and Next Generation Sequencing Approaches to Evaluate HIV-1 Virus in Blood Compartments

The implementation of antiretroviral treatment combined with the monitoring of drug resistance mutations improves the quality of life of HIV-1 positive patients. The drug resistance mutation patterns and viral genotypes are currently analyzed by DNA sequencing of the virus in the plasma of patients. However, the virus compartmentalizes, and different T cell subsets may harbor distinct viral subsets. In this study, we compared the patterns of HIV distribution in cell-free (blood plasma) and cell-associated viruses (peripheral blood mononuclear cells, PBMCs) derived from ART-treated patients by using Sanger sequencing- and Next-Generation sequencing-based HIV assay. CD4⁺CD45RA⁻RO⁺ memory T-cells were isolated from PBMCs using a BD FACSAria instrument. HIV pol (protease and reverse transcriptase) was RT-PCR or PCR amplified from the plasma and the T-cell subset, respectively. Sequences were obtained using Sanger sequencing and Next-Generation Sequencing (NGS). Sanger sequences were aligned and edited using RECall software (beta v3.03). The Stanford HIV database was used to evaluate drug resistance mutations. Illumina MiSeq platform and HyDRA Web were used to generate and analyze NGS data, respectively. Our results show a high correlation between Sanger sequencing and NGS results. However, some major and minor drug resistance mutations were only observed by NGS, albeit at different frequencies. Analysis of low-frequency drugs resistance mutations and virus distribution in the blood compartments may provide information to allow a more sustainable response to therapy and better disease management.

Pretreatment drug resistance in a large countrywide Ethiopian HIV-1C cohort: a comparison of Sanger and high-throughput sequencing

Baseline plasma samples of 490 randomly selected antiretroviral therapy (ART) naïve patients from seven hospitals participating in the first nationwide Ethiopian HIV-1 cohort were analysed for surveillance drug resistance mutations (sDRM) by population based Sanger sequencing (PBSS). Also next generation sequencing (NGS) was used in a subset of 109 baseline samples of patients. Treatment outcome after 6– and 12–months was assessed by on-treatment (OT) and intention-to-treat (ITT) analyses. Transmitted drug resistance (TDR) was detected in 3.9% (18/461) of successfully sequenced samples by PBSS. However, NGS detected sDRM more often (24%; 26/109) than PBSS (6%; 7/109) (p = 0.0001) and major integrase strand transfer inhibitors (INSTI) DRMs were also found in minor viral variants from five patients. Patients with sDRM had more frequent treatment failure in both OT and ITT analyses. The high rate of TDR by NGS and the identification of preexisting INSTI DRMs in minor wild-type HIV-1 subtype C viral variants infected Ethiopian patients underscores the importance of TDR surveillance in low– and middle–income countries and shows added value of high-throughput NGS in such studies.

Impact of the HIV-1 genetic background and HIV-1 population size on the evolution of raltegravir resistance

BACKGROUND

Emergence of resistance against integrase inhibitor raltegravir in human immunodeficiency virus type 1 (HIV-1) patients is generally associated with selection of one of three signature mutations: Y143C/R, Q148K/H/R or N155H, representing three distinct resistance pathways. The mechanisms that drive selection of a specific pathway are still poorly understood. We investigated the impact of the HIV-1 genetic background and population dynamics on the emergence of raltegravir resistance. Using deep sequencing we analyzed the integrase coding sequence (CDS) in longitudinal samples from five patients who initiated raltegravir plus optimized background therapy at viral loads > 5000 copies/ml. To investigate the role of the HIV-1 genetic background we created recombinant viruses containing the viral integrase coding region from pre-raltegravir samples from two patients in whom raltegravir resistance developed through different pathways. The in vitro selections performed with these recombinant viruses were designed to mimic natural population bottlenecks.

RESULTS

Deep sequencing analysis of the viral integrase CDS revealed that the virological response to raltegravir containing therapy inversely correlated with the relative amount of unique sequence variants that emerged suggesting diversifying selection during drug pressure. In 4/5 patients multiple signature mutations representing different resistance pathways were observed. Interestingly, the resistant population can consist of a single resistant variant that completely dominates the population but also of multiple variants from different resistance pathways that coexist in the viral population. We also found evidence for increased diversification after stronger bottlenecks. In vitro selections with low viral titers, mimicking population bottlenecks, revealed that both recombinant viruses and HXB2 reference virus were able to select mutations from different resistance pathways, although typically only one resistance pathway emerged in each individual culture.

CONCLUSIONS

The generation of a specific raltegravir resistant variant is not predisposed in the genetic background of the viral integrase CDS. Typically, in the early phases of therapy failure the sequence space is explored and multiple resistance pathways emerge and then compete for dominance which frequently results in a switch of the dominant population over time towards the fittest variant or even multiple variants of similar fitness that can coexist in the viral population.

The Role of HIV-1 Drug-Resistant Minority Variants in Treatment Failure

Human immunodeficiency virus type 1 (HIV-1) drug resistance genotyping is recommended to help in the selection of antiretroviral therapy and to prevent virologic failure. There are several ultrasensitive assays able to detect HIV-1 drug-resistance minority variants (DRMVs) not detectable by standard population sequencing-based HIV genotyping assays. Presence of these DRMVs has been shown to be clinically relevant, but its impact does not appear to be uniform across drug classes. In this review, we summarize key evidence for the clinical impact of DRMVs across drug classes for both antiretroviral treatment-naive and antiretroviral treatment-experienced patients, and highlight areas where more supporting evidence is needed.

Lack of impact of pre-existing T97A HIV-1 integrase mutation on integrase strand transfer inhibitor resistance and treatment outcome

T97A is an HIV-1 integrase polymorphism associated with integrase strand transfer inhibitor (INSTI) resistance. Using pooled data from 16 clinical studies, we investigated the prevalence of T97A (pre-existing and emergent) and its impact on INSTI susceptibility and treatment response in INSTI-naive patients who enrolled on elvitegravir (EVG)- or raltegravir (RAL)-based regimens. Prior to INSTI-based therapy, primary INSTI resistance-associated mutations (RAMs) were absent and T97A pre-existed infrequently (1.4%; 47 of 3367 integrase sequences); most often among non-B (5.3%) than B (0.9%) HIV-1 subtypes. During INSTI-based therapy, few patients experienced virologic failure with emergent INSTI RAMs (3%; 122 of 3881 patients), among whom T97A emerged infrequently in the presence (n = 6) or absence (n = 8) of primary INSTI RAMs. A comparison between pre-existing and emergent T97A patient populations (i.e., in the absence of primary INSTI RAMs) showed no significant differences in EVG or RAL susceptibility in vitro. Furthermore, among all T97A-containing viruses tested, only 38-44% exhibited reduced susceptibility to EVG and/or RAL (all of low magnitude; <11-fold), while all maintained susceptibility to dolutegravir. Of the patients with pre-existing T97A, 17 had available clinical follow-up: 16 achieved virologic suppression and 1 maintained T97A and INSTI sensitivity without further resistance development. Overall, T97A is an infrequent integrase polymorphism that is enriched among non-B HIV-1 subtypes and can confer low-level reduced susceptibility to EVG and/or RAL. However, detection of T97A does not affect response to INSTI-based therapy with EVG or RAL. These results suggest a very low risk of initiating INSTI-based therapy in patients with pre-existing T97A.

Long Single-Molecule Reads Can Resolve the Complexity of the Influenza Virus Composed of Rare, Closely Related Mutant Variants

As a result of a high rate of mutations and recombination events, an RNA-virus exists as a heterogeneous "swarm" of mutant variants. The long read length offered by single-molecule sequencing technologies allows each mutant variant to be sequenced in a single pass. However, high error rate limits the ability to reconstruct heterogeneous viral population composed of rare, related mutant variants. In this article, we present two single-nucleotide variants (2SNV), a method able to tolerate the high error rate of the single-molecule protocol and reconstruct mutant variants. 2SNV uses linkage between single-nucleotide variations to efficiently distinguish them from read errors. To benchmark the sensitivity of 2SNV, we performed a single-molecule sequencing experiment on a sample containing a titrated level of known viral mutant variants. Our method is able to accurately reconstruct clone with frequency of 0.2% and distinguish clones that differed in only two nucleotides distantly located on the genome. 2SNV outperforms existing methods for full-length viral mutant reconstruction.

Efficacy and safety of antiretroviral regimens including raltegravir to treat HIV-infected patients with hemophilia

When treating HIV-infected patients with hemophilia, adverse drug reactions and interactions and the effect of treatment on bleeding disorders must be considered. Raltegravir is the first HIV integrase inhibitor, but its use in patients with hemophilia is rarely reported. Nine HIV-positive patients with hemophilia were retrospectively studied with a focus on the virological response, changes in the CD4 count, the tendency to bleed, and the response to replacement therapy before and after raltegravir-based antiretroviral therapy (ART). The nine patients were highly treatment-experienced patients and they received raltegravir-based ART for at least nine months. The patients had their own reasons for changing to raltegravir-based ART. During treatment, the CD4 count increased progressively in four patients, with a median absolute increase of 233 cells/mm(3), while the count stabilized in the remaining five patients. Two previous recipients of lopinavir/ritonavir (LPV/r) who failed to respond to lamivudine (3TC) + zidovudine (ZDV) + efavirenz (EFV) had a viral rebound. Genotyping indicated multidrug resistance to nucleoside reverse transcriptase inhibitors (NRTIs) and non-nucleoside reverse transcriptase inhibitors (NNRTIs). A pattern of resistance to raltegravir was evident, including the primary mutation N155H and the secondary mutation T97A. In the two patients, the tendency to bleed decreased markedly and monthly usage of clotting factor VIII decreased significantly decreased. In the remaining seven patients, the viral load remained < 40 copies/mL, there was no evidence of an increased tendency to bleed, and no evidence of changes in the response to replacement therapy. All of the patients had a stable condition with no signs of disease progression and no serious adverse reactions. Results indicated that Raltegravir-based therapy offered a safe and well-tolerated option for HIV-positive patients with hemophilia.

Comparison of Major and Minor Viral SNPs Identified through Single Template Sequencing and Pyrosequencing in Acute HIV-1 Infection

Massively parallel sequencing (MPS) technologies, such as 454-pyrosequencing, allow for the identification of variants in sequence populations at lower levels than consensus sequencing and most single-template Sanger sequencing experiments. We sought to determine if the greater depth of population sampling attainable using MPS technology would allow detection of minor variants in HIV founder virus populations very early in infection in instances where Sanger sequencing detects only a single variant. We compared single nucleotide polymorphisms (SNPs) during acute HIV-1 infection from 32 subjects using both single template Sanger and 454-pyrosequencing. Pyrosequences from a median of 2400 viral templates per subject and encompassing 40% of the HIV-1 genome, were compared to a median of five individually amplified near full-length viral genomes sequenced using Sanger technology. There was no difference in the consensus nucleotide sequences over the 3.6kb compared in 84% of the subjects infected with single founders and 33% of subjects infected with multiple founder variants: among the subjects with disagreements, mismatches were found in less than 1% of the sites evaluated (of a total of nearly 117,000 sites across all subjects). The majority of the SNPs observed only in pyrosequences were present at less than 2% of the subject’s viral sequence population. These results demonstrate the utility of the Sanger approach for study of early HIV infection and provide guidance regarding the design, utility and limitations of population sequencing from variable template sources, and emphasize parameters for improving the interpretation of massively parallel sequencing data to address important questions regarding target sequence evolution.

Integrase inhibitor (INI) genotypic resistance in treatment-naive and raltegravir-experienced patients infected with diverse HIV-1 clades

OBJECTIVES

The aim of this study was to characterize the prevalence and patterns of genotypic integrase inhibitor (INI) resistance in relation to HIV-1 clade.

METHODS

The cohort comprised 533 INI-naive subjects and 255 raltegravir recipients with viraemia who underwent integrase sequencing in routine care across Europe, including 134/533 (25.1%) and 46/255 (18.0%), respectively, with non-B clades (A, C, D, F, G, CRF01, CRF02, other CRFs, complex).

RESULTS

No major INI resistance-associated mutations (RAMs) occurred in INI-naive subjects. Among raltegravir recipients with viraemia (median 3523 HIV-1 RNA copies/mL), 113/255 (44.3%) had one or more major INI RAMs, most commonly N155H (45/255, 17.6%), Q148H/R/K + G140S/A (35/255, 13.7%) and Y143R/C/H (12/255, 4.7%). In addition, four (1.6%) raltegravir recipients showed novel mutations at recognized resistance sites (E92A, S147I, N155D, N155Q) and novel mutations at other integrase positions that were statistically associated with raltegravir exposure (K159Q/R, I161L/M/T/V, E170A/G). Comparing subtype B with non-B clades, Q148H/R/K occurred in 42/209 (20.1%) versus 2/46 (4.3%) subjects (P = 0.009) and G140S/A occurred in 36/209 (17.2%) versus 1/46 (2.2%) subjects (P = 0.005). Intermediate- to high-level cross-resistance to twice-daily dolutegravir was predicted in 40/255 (15.7%) subjects, more commonly in subtype B versus non-B clades (39/209, 18.7% versus 1/46, 2.2%; P = 0.003). A glycine (G) to serine (S) substitution at integrase position 140 required one nucleotide change in subtype B and two nucleotide changes in all non-B clades.

CONCLUSIONS

No major INI resistance mutations occurred in INI-naive subjects. Reduced occurrence of Q148H/R/K + G140S/A was seen in non-B clades versus subtype B, and was explained by the higher genetic barrier to the G140S mutation observed in all non-B clades analysed.

Highly frequent HIV-1 minority resistant variants at baseline of the ANRS 139 TRIO trial had a limited impact on virological response

OBJECTIVES

To assess the prevalence of minority resistant variants (MRVs) at baseline and their impact on the virological response. The ANRS 139 TRIO trial evaluated the combination of raltegravir, etravirine and darunavir, plus an optimized background therapy, in 87% of cases. Patients were highly experienced and harboured multiresistant viruses, but were naive to the three drugs, and showed a high level of virological suppression.

METHODS

Ultra-deep sequencing of reverse transcriptase, protease and integrase regions was performed at the trial baseline, and sequences were interpreted according to the ANRS algorithm. MRVs were assessed using MiSeq and 454 technologies (limit of detection 1%).

RESULTS

At baseline, minority variants with at least one NRTI, one NNRTI, one PI, one major PI or an integrase inhibitor resistance-associated mutation were present in 46%, 45%, 68%, 24% and 13% of patients, respectively. When minority variants are taken into account, the prevalence of resistance to etravirine, darunavir and raltegravir at baseline was 29%, 40% and 9%, respectively. No difference was observed in the prevalence of MRVs between patients with virological failure and those with virological success, except a trend for patients exhibiting baseline etravirine MRVs (50% versus 26%, P = 0.09).

CONCLUSIONS

We have shown a high level of MRVs at baseline in highly pre-treated patients harbouring multiresistant viruses. However, these MRVs were not associated with an increased risk of virological failure, except for a trend for etravirine MRVs.

Extensive recombination due to heteroduplexes generates large amounts of artificial gene fragments during PCR

Artificial recombinants can be generated during PCR when more than two genetically distinct templates coexist in a single PCR reaction. These recombinant amplicons can lead to the false interpretation of genetic diversity and incorrect identification of biological phenotypes that do not exist in vivo. We investigated how recombination between 2 or 35 genetically distinct HIV-1 genomes was affected by different PCR conditions using the parallel allele-specific sequencing (PASS) assay and the next generation sequencing method. In a standard PCR condition, about 40% of amplicons in a PCR reaction were recombinants. The high recombination frequency could be significantly reduced if the number of amplicons in a PCR reaction was below a threshold of 10¹³–10¹⁴ using low thermal cycles, fewer input templates, and longer extension time. Heteroduplexes (each DNA strand from a distinct template) were present at a large proportion in the PCR products when more thermal cycles, more templates, and shorter extension time were used. Importantly, the majority of recombinants were identified in heteroduplexes, indicating that the recombinants were mainly generated through heteroduplexes. Since prematurely terminated extension fragments can form heteroduplexes by annealing to different templates during PCR amplification, recombination has a better chance to occur with samples containing different genomes when the number of amplicons accumulate over the threshold. New technologies are warranted to accurately characterize complex quasispecies gene populations.

Accurate viral population assembly from ultra-deep sequencing data

MOTIVATION

Next-generation sequencing technologies sequence viruses with ultra-deep coverage, thus promising to revolutionize our understanding of the underlying diversity of viral populations. While the sequencing coverage is high enough that even rare viral variants are sequenced, the presence of sequencing errors makes it difficult to distinguish between rare variants and sequencing errors.

RESULTS

In this article, we present a method to overcome the limitations of sequencing technologies and assemble a diverse viral population that allows for the detection of previously undiscovered rare variants. The proposed method consists of a high-fidelity sequencing protocol and an accurate viral population assembly method, referred to as Viral Genome Assembler (VGA). The proposed protocol is able to eliminate sequencing errors by using individual barcodes attached to the sequencing fragments. Highly accurate data in combination with deep coverage allow VGA to assemble rare variants. VGA uses an expectation-maximization algorithm to estimate abundances of the assembled viral variants in the population. RESULTS on both synthetic and real datasets show that our method is able to accurately assemble an HIV viral population and detect rare variants previously undetectable due to sequencing errors. VGA outperforms state-of-the-art methods for genome-wide viral assembly. Furthermore, our method is the first viral assembly method that scales to millions of sequencing reads.

AVAILABILITY

Our tool VGA is freely available at http://genetics.cs.ucla.edu/vga/

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.

Comparison of illumina and 454 deep sequencing in participants failing raltegravir-based antiretroviral therapy

BACKGROUND

The impact of raltegravir-resistant HIV-1 minority variants (MVs) on raltegravir treatment failure is unknown. Illumina sequencing offers greater throughput than 454, but sequence analysis tools for viral sequencing are needed. We evaluated Illumina and 454 for the detection of HIV-1 raltegravir-resistant MVs.

METHODS

A5262 was a single-arm study of raltegravir and darunavir/ritonavir in treatment-naïve patients. Pre-treatment plasma was obtained from 5 participants with raltegravir resistance at the time of virologic failure. A control library was created by pooling integrase clones at predefined proportions. Multiplexed sequencing was performed with Illumina and 454 platforms at comparable costs. Illumina sequence analysis was performed with the novel snp-assess tool and 454 sequencing was analyzed with V-Phaser.

RESULTS

Illumina sequencing resulted in significantly higher sequence coverage and a 0.095% limit of detection. Illumina accurately detected all MVs in the control library at ≥0.5% and 7/10 MVs expected at 0.1%. 454 sequencing failed to detect any MVs at 0.1% with 5 false positive calls. For MVs detected in the patient samples by both 454 and Illumina, the correlation in the detected variant frequencies was high (R2 = 0.92, P<0.001). Illumina sequencing detected 2.4-fold greater nucleotide MVs and 2.9-fold greater amino acid MVs compared to 454. The only raltegravir-resistant MV detected was an E138K mutation in one participant by Illumina sequencing, but not by 454.

CONCLUSIONS

In participants of A5262 with raltegravir resistance at virologic failure, baseline raltegravir-resistant MVs were rarely detected. At comparable costs to 454 sequencing, Illumina demonstrated greater depth of coverage, increased sensitivity for detecting HIV MVs, and fewer false positive variant calls.

Sensitive deep-sequencing-based HIV-1 genotyping assay to simultaneously determine susceptibility to protease, reverse transcriptase, integrase, and maturation inhibitors, as well as HIV-1 coreceptor tropism

With 29 individual antiretroviral drugs available from six classes that are approved for the treatment of HIV-1 infection, a combination of different phenotypic and genotypic tests is currently needed to monitor HIV-infected individuals. In this study, we developed a novel HIV-1 genotypic assay based on deep sequencing (DeepGen HIV) to simultaneously assess HIV-1 susceptibilities to all drugs targeting the three viral enzymes and to predict HIV-1 coreceptor tropism. Patient-derived gag-p2/NCp7/p1/p6/pol-PR/RT/IN- and env-C2V3 PCR products were sequenced using the Ion Torrent Personal Genome Machine. Reads spanning the 3' end of the Gag, protease (PR), reverse transcriptase (RT), integrase (IN), and V3 regions were extracted, truncated, translated, and assembled for genotype and HIV-1 coreceptor tropism determination. DeepGen HIV consistently detected both minority drug-resistant viruses and non-R5 HIV-1 variants from clinical specimens with viral loads of ≥1,000 copies/ml and from B and non-B subtypes. Additional mutations associated with resistance to PR, RT, and IN inhibitors, previously undetected by standard (Sanger) population sequencing, were reliably identified at frequencies as low as 1%. DeepGen HIV results correlated with phenotypic (original Trofile, 92%; enhanced-sensitivity Trofile assay [ESTA], 80%; TROCAI, 81%; and VeriTrop, 80%) and genotypic (population sequencing/Geno2Pheno with a 10% false-positive rate [FPR], 84%) HIV-1 tropism test results. DeepGen HIV (83%) and Trofile (85%) showed similar concordances with the clinical response following an 8-day course of maraviroc monotherapy (MCT). In summary, this novel all-inclusive HIV-1 genotypic and coreceptor tropism assay, based on deep sequencing of the PR, RT, IN, and V3 regions, permits simultaneous multiplex detection of low-level drug-resistant and/or non-R5 viruses in up to 96 clinical samples. This comprehensive test, the first of its class, will be instrumental in the development of new antiretroviral drugs and, more importantly, will aid in the treatment and management of HIV-infected individuals.

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.

Association of preexisting drug-resistance mutations and treatment failure in hepatitis B patients

The role of preexisting minority drug-resistance mutations in treatment failure has not been fully understood in chronic hepatitis B patients. To understand mechanisms of drug resistance, we analyzed drug-resistance mutations in 46 treatment-failure patients and in 29 treatment-naïve patients and determined linkage patterns of the drug-resistance mutations in individual viral genomes using a highly sensitive parallel allele-specific sequencing (PASS) method. Lamivudine resistance (LAMr) mutations were predominant in treatment-failure patients, irrespective of the inclusion of LAM in the regimen. The primary LAMr mutations M204V and M204I were detected in 100% and 30% of the treatment-failure patients, respectively. Two secondary LAMr mutations (L180M and V173L) were also found in most treatment-failure patients (87% and 78%, respectively). The linkages containing these three mutations dominated the resistant viruses. Importantly, minority LAMr mutations present in <2% of the viral population were detected in 83% of the treatment-naïve patients. Moreover, the low-frequency same linked LAMr mutations (<0.15%) were detected in 24% of the treatment-naïve patients. Our results demonstrate that the selection of preexisting minority linked LAMr mutations may be an important mechanism for the rapid development of LAM resistance, caution the continuous use of LAM to treat drug-experienced and -naïve hepatitis B patients, and underline the importance of the detection of minority single and linked drug-resistance mutations before initiating antiviral therapy.

Emerging patterns and implications of HIV-1 integrase inhibitor resistance

PURPOSE OF REVIEW

This review highlights recent data on the pathways of resistance that impact the clinical activity of first-generation and second-generation integrase inhibitors.

RECENT FINDINGS

Raltegravir (RAL) and elvitegravir (EVG) are highly efficacious in first-line antiretroviral therapy, with small numbers of virological failures observed in clinical trials. Durable activity in treatment-experienced patients requires a fully supportive background regimen. RAL and EVG show a low-to-moderate genetic barrier to resistance and extensive cross-resistance, which preclude their sequential use. Resistance to dolutegravir (DTG) is not selected as readily in vitro and has not emerged in studies of treatment-naïve patients to date. Both in vitro and in vivo, DTG retains activity against several RAL and EVG resistant strains, but susceptibility is variably impaired by multiple mutations within the G148 pathway, which are common after RAL or EVG failure. Cross-resistance can be partially overcome by doubling DTG dosing to twice daily, but durability of responses remains dependent on a supportive background regimen. There is variability in the integrase gene of circulating HIV strains, which does not appear to reduce drug activity, although it may influence the emergence and evolution of integrase resistance. Transmission of integrase resistance remains rare but surveillance is required.

SUMMARY

Integrase inhibitors provide a potent option for the treatment of HIV infection. Drug resistance remains a challenge, which may be partially overcome by the introduction of second-generation compounds. Prompt management of RAL and EVG failure is required to prevent the accumulation of multiple resistance mutations that reduce DTG susceptibility.

Clinical implications of HIV-1 minority variants

Technologic advances in human immunodeficiency virus type 1 (HIV-1) sequencing have revolutionized the study of antiretroviral drug resistance and are increasingly moving from the laboratory to clinical practice. These techniques are able to detect HIV-1 drug resistance mutations present at low frequencies not detectable by current HIV-1 genotyping assays. For a number of commonly used antiretroviral medications, such as nonnucleoside reverse transcriptase inhibitors, the detection of these drug-resistant minority variants significantly increases the risk of treatment failure. The level of evidence, however, is insufficient to determine the impact of HIV-1 minority variants for several other classes of antiretroviral medications. Clinicians should be aware of the novel technologies that are moving into routine clinical use and the clinical implications of HIV-1 minority variants. Additional studies are needed to determine the optimal platform for clinical application of these new technologies and to provide guidance to clinicians on the type and frequency of clinically important HIV-1 minority variants.

Factors associated with virological success with raltegravir-containing regimens and prevalence of raltegravir-resistance-associated mutations at failure in the ARCA database

Raltegravir (RAL) is the only licensed human immunodeficiency virus (HIV) integrase inhibitor. The factors associated with the virological response to RAL-containing regimens and the prevalence of integrase mutations associated with RAL failure deserve further investigation. From the Antiretroviral Resistance Cohort Analysis database, we selected triple-class-experienced subjects failing their current treatment with complete treatment history available. Selection criteria included HIV-RNA, CD4 count and HIV genotype within 3 months of RAL initiation. Factors associated with 24-week response were analysed; genotypic sensitivity scores (GSS) and weighted-GSS were evaluated. Virological response was achieved in 74.3% of 105 subjects. Mutations associated with RAL failure were detected in 12/24 subjects with an integrase genotype, with the prevalence of Q148H + G140S. Each extra unit of GSS (p 0.05, OR 2.62; 95% CI 1.00-6.87). was found to be a associated with response. Weighted-GSS had borderline statistical significance (p 0.063, OR 2.04; 95% CI 0.96-4.33) When stratifying for different cut-offs (<1 as reference, 1-1.49, ≥1.5), a borderline significant increase in the probability of response appeared for GSS ≥1.5 (p 0.053, OR 4.00; 95% CI 0.98-16.25). GSS ≥1 showed the highest sensitivity, 82.6%. Receiver operating characteristic curves depicted the widest area under the curve (0.663, p 0.054) of GSS ≥1. Unresponsiveness to RAL-containing regimens among triple-class-experienced subjects was low. The activity of the background regimen was strongly associated with response. Although few integrase genotypes were available at failure, half of these were without integrase resistance mutations. The substantial rate of RAL failure in the absence of known RAL-resistance mutations may be associated with adherence issues and this issue warrants further analysis in longer observations.

Impact of minority nonnucleoside reverse transcriptase inhibitor resistance mutations on resistance genotype after virologic failure

Drug-resistant human immunodeficiency virus type 1 (HIV-1) minority variants increase the risk of virologic failure for first-line nonnucleoside reverse transcriptase inhibitor (NNRTI)-based regimens. We performed a pooled analysis to evaluate the relationship between NNRTI-resistant minority variants and the likelihood and types of resistance mutations detected at virologic failure. In multivariable logistic regression analysis, higher NNRTI minority variant copy numbers, non-white race, and nevirapine use were associated with a higher risk of NNRTI resistance at virologic failure. Among participants on efavirenz, K103N was the most frequently observed resistance mutation at virologic failure regardless of the baseline minority variant. However, the presence of baseline Y181C minority variant was associated with a higher probability of Y181C detection after virologic failure. NNRTI regimen choice and preexisting NNRTI-resistant minority variants were both associated with the probability and type of resistance mutations detected after virologic failure.

Drug resistance in HIV-1

PURPOSE OF THE REVIEW

Changing antiretroviral regimens and the introduction of new antiretroviral drugs have altered drug resistance patterns in human immunodeficiency virus type 1 (HIV-1). This review summarizes recent information on antiretroviral drug resistance.

RECENT FINDINGS

As tenofovir and abacavir have replaced zidovudine and stavudine in antiretroviral regimens, thymidine analog resistance mutations have become less common in patients failing antiretroviral therapy in developed countries. Similarly, the near universal use of ritonavir-boosted protease inhibitors (PI) in place of unboosted PIs has made the selection of PI resistance mutations uncommon in patients failing a first-line or second-line PI regimen. The challenge of treating patients with multidrug-resistant HIV-1 has largely been addressed by the advent of newer PIs, second-generation non-nucleoside reverse transcriptase inhibitors and drugs in novel classes, including integrase inhibitors and CCR5 antagonists. Resistance to these newer agents can emerge, however, resulting in the appearance of novel drug resistance mutations in the HIV-1 polymerase, integrase and envelope genes.

SUMMARY

New drugs make possible the effective treatment of multidrug-resistant HIV-1, but the activity of these drugs may be limited by the appearance of novel drug resistance mutations.

Use of four next-generation sequencing platforms to determine HIV-1 coreceptor tropism

HIV-1 coreceptor tropism assays are required to rule out the presence of CXCR4-tropic (non-R5) viruses prior treatment with CCR5 antagonists. Phenotypic (e.g., Trofile™, Monogram Biosciences) and genotypic (e.g., population sequencing linked to bioinformatic algorithms) assays are the most widely used. Although several next-generation sequencing (NGS) platforms are available, to date all published deep sequencing HIV-1 tropism studies have used the 454™ Life Sciences/Roche platform. In this study, HIV-1 co-receptor usage was predicted for twelve patients scheduled to start a maraviroc-based antiretroviral regimen. The V3 region of the HIV-1 env gene was sequenced using four NGS platforms: 454™, PacBio® RS (Pacific Biosciences), Illumina®, and Ion Torrent™ (Life Technologies). Cross-platform variation was evaluated, including number of reads, read length and error rates. HIV-1 tropism was inferred using Geno2Pheno, Web PSSM, and the 11/24/25 rule and compared with Trofile™ and virologic response to antiretroviral therapy. Error rates related to insertions/deletions (indels) and nucleotide substitutions introduced by the four NGS platforms were low compared to the actual HIV-1 sequence variation. Each platform detected all major virus variants within the HIV-1 population with similar frequencies. Identification of non-R5 viruses was comparable among the four platforms, with minor differences attributable to the algorithms used to infer HIV-1 tropism. All NGS platforms showed similar concordance with virologic response to the maraviroc-based regimen (75% to 80% range depending on the algorithm used), compared to Trofile (80%) and population sequencing (70%). In conclusion, all four NGS platforms were able to detect minority non-R5 variants at comparable levels suggesting that any NGS-based method can be used to predict HIV-1 coreceptor usage.

Impact of immune escape mutations on HIV-1 fitness in the context of the cognate transmitted/founder genome

BACKGROUND

A modest change in HIV-1 fitness can have a significant impact on viral quasispecies evolution and viral pathogenesis, transmission and disease progression. To determine the impact of immune escape mutations selected by cytotoxic T lymphocytes (CTL) on viral fitness in the context of the cognate transmitted/founder (T/F) genome, we developed a new competitive fitness assay using molecular clones of T/F genomes lacking exogenous genetic markers and a highly sensitive and precise parallel allele-specific sequencing (PASS) method.

RESULTS

The T/F and mutant viruses were competed in CD4+ T-cell enriched cultures, relative proportions of viruses were assayed after repeated cell-free passage, and fitness costs were estimated by mathematical modeling. Naturally occurring HLA B57-restricted mutations involving the TW10 epitope in Gag and two epitopes in Tat/Rev and Env were assessed independently and together. Compensatory mutations which restored viral replication fitness were also assessed. A principal TW10 escape mutation, T242N, led to a 42% reduction in replication fitness but V247I and G248A mutations in the same epitope restored fitness to wild-type levels. No fitness difference was observed between the T/F and a naturally selected variant carrying the early CTL escape mutation (R355K) in Env and a reversion mutation in the Tat/Rev overlapping region.

CONCLUSIONS

These findings reveal a broad spectrum of fitness costs to CTL escape mutations in T/F viral genomes, similar to recent findings reported for neutralizing antibody escape mutations, and highlight the extraordinary plasticity and adaptive potential of the HIV-1 genome. Analysis of T/F genomes and their evolved progeny is a powerful approach for assessing the impact of composite mutational events on viral fitness.

Development of elvitegravir resistance and linkage of integrase inhibitor mutations with protease and reverse transcriptase resistance mutations

Failure of antiretroviral regimens containing elvitegravir (EVG) and raltegravir (RAL) can result in the appearance of integrase inhibitor (INI) drug-resistance mutations (DRMs). While several INI DRMs have been identified, the evolution of EVG DRMs and the linkage of these DRMs with protease inhibitor (PI) and reverse transcriptase inhibitor (RTI) DRMs have not been studied at the clonal level. We examined the development of INI DRMs in 10 patients failing EVG-containing regimens over time, and the linkage of INI DRMs with PI and RTI DRMs in these patients plus 6 RAL-treated patients. A one-step RT-nested PCR protocol was used to generate a 2.7 kB amplicon that included the PR, RT, and IN coding region, and standard cloning and sequencing techniques were used to determine DRMs in 1,277 clones (mean 21 clones per time point). Results showed all patients had multiple PI, NRTI, and/or NNRTI DRMs at baseline, but no primary INI DRM. EVG-treated patients developed from 2 to 6 strains with different primary INI DRMs as early as 2 weeks after initiation of treatment, predominantly as single mutations. The prevalence of these strains fluctuated and new strains, and/or strains with new combinations of INI DRMs, developed over time. Final failure samples (weeks 14 to 48) typically showed a dominant strain with multiple mutations or N155H alone. Single N155H or multiple mutations were also observed in RAL-treated patients at virologic failure. All patient strains showed evidence of INI DRM co-located with single or multiple PI and/or RTI DRMs on the same viral strand. Our study shows that EVG treatment can select for a number of distinct INI-resistant strains whose prevalence fluctuates over time. Continued appearance of new INI DRMs after initial INI failure suggests a potent, highly dynamic selection of INI resistant strains that is unaffected by co-location with PI and RTI DRMs.

Ultra-deep sequencing of HIV-1 reverse transcriptase before start of an NNRTI-based regimen in treatment-naive patients

There are conflicting data on the impact of low frequency HIV-1 drug-resistant mutants on the response of first-line highly active antiretroviral therapy (HAART), more specifically containing a NNRTI. As population sequencing does not detect resistant viruses representing less than 15-25% of the viral population, more sensitive techniques have been developed but still need clinical validation. We evaluated ultra-deep sequencing (UDPS), recently more available and affordable, as a tool for the detection of HIV-1 minority species carrying drug resistant mutation (DRM) in a clinical setting. A retrospective analysis of the reverse transcriptase (RT) gene of plasma HIV-1 from 70 patients starting a NNRTI based regimen was performed. Minority populations were defined as representing > 1% and < 20% of the total viral population. Using UDPS, we could not confirm an association between the presence of low minority variants harbouring RT mutations at the start of therapy and primary or secondary therapeutic failure.

Study of genotypic and phenotypic HIV-1 dynamics of integrase mutations during raltegravir treatment: a refined analysis by ultra-deep 454 pyrosequencing

BACKGROUND

The dynamics of raltegravir-resistant variants and their impact on virologic response in 23 HIV-1-infected patients, who started a salvage raltegravir-containing regimen, were investigated.

METHODS

Integrase population sequencing and Ultra-Deep-454 Pyrosequencing (UDPS) were performed on plasma samples at baseline and at raltegravir failure. All integrase mutations detected at a frequency ≥1% were considered to be reliable for the UDPS analyses. Phylogenetic and phenotypic resistance analyses were also performed.

RESULTS

At baseline, primary resistance mutations were not detected by both population and UDPS genotypic assays; few secondary mutations (T97A-V151I-G163R) were rarely detected and did not show any statistically association either with virologic response at 24-weeks or with the development of resistant variants at failure. At UDPS, not all resistant variants appearing early during treatment evolved as major populations during failure; only specific resistance pathways (Y143R-Q148H/R-N155H) associated with an increased rate of fitness and phenotypic resistance were selected.

CONCLUSIONS

Resistance to raltegravir in integrase strand transfer inhibitor-naive patients remains today a rare event, which might be changed by future extensive use of such drugs. In our study, pathways of resistance at failure were not predicted by baseline mutations, suggesting that evolution plus stochastic selection plays a major role in the appearance of integrase-resistance mutations, whereas fitness and resistance are dominant factors acting for the late selection of resistant quasispecies.

Efficacy of a nucleoside-sparing regimen of darunavir/ritonavir plus raltegravir in treatment-naive HIV-1-infected patients (ACTG A5262)

OBJECTIVE

To explore darunavir/ritonavir (DRV/r) plus raltegravir (RAL) combination therapy in antiretroviral-naive patients.

DESIGN

Phase IIb, single-arm, open-label, multicenter study.

METHODS

One hundred and twelve antiretroviral-naive, HIV-1-infected patients received DRV/r 800/100 mg once daily and RAL 400 mg twice daily. Primary endpoint was virologic failure by week 24. Virologic failure was defined as confirmed viral load of 1000 copies/ml or more at week 12, or an increase of more than 0.5 log(10) copies/ml in viral load from week 4 to 12, or a confirmed viral load of more than 50 copies/ml at or after week 24. Protease and integrase genes were sequenced in patients experiencing virologic failure.

RESULTS

Virologic failure rate was 16% [95% confidence interval (CI) 10-24] by week 24 and 26% (95% CI 19-36) by week 48 in an intent-to-treat analysis. Viral load at virologic failure was 51-200 copies/ml in 17/28 failures. Adjusting for age and sex, virologic failure was associated with baseline viral load of more than 100,000 copies/ml [hazard ratio 3.76, 95% CI (1.52-9.31), P = 0.004] and lower CD4 cell count [0.77 per 100 cells/μl increase (95% CI 0.61-0.98), P = 0.037]. When trough RAL concentrations were included as a time-varying covariate in the analysis, virologic failure remained associated with baseline viral load more than 100,000 copies/ml [hazard ratio = 4.67 (95% CI 1.93-11.25), P < 0.001], whereas RAL level below detection limit in plasma at one or more previous visits was associated with increased hazard [hazard ratio = 3.42 (95% CI 1.41-8.26), P = 0.006]. All five participants with integrase mutations during virologic failure had baseline viral load more than 100,000 copies/ml.

CONCLUSION

DRV/r plus RAL was effective and well tolerated in most patients, but virologic failure and integrase resistance were common, particularly in patients with baseline viral load more than 100,000 copies/ml.

Primary infection by a human immunodeficiency virus with atypical coreceptor tropism

The great majority of human immunodeficiency virus type 1 (HIV-1) strains enter CD4+ target cells by interacting with one of two coreceptors, CCR5 or CXCR4. Here we describe a transmitted/founder (T/F) virus (ZP6248) that was profoundly impaired in its ability to utilize CCR5 and CXCR4 coreceptors on multiple CD4+ cell lines as well as primary human CD4+ T cells and macrophages in vitro yet replicated to very high titers (>80 million RNA copies/ml) in an acutely infected individual. Interestingly, the envelope (Env) glycoprotein of this clade B virus had a rare GPEK sequence in the crown of its third variable loop (V3) rather than the consensus GPGR sequence. Extensive sequencing of sequential plasma samples showed that the GPEK sequence was present in virtually all Envs, including those from the earliest time points after infection. The molecularly cloned (single) T/F virus was able to replicate, albeit poorly, in cells obtained from ccr5Δ32 homozygous donors. The ZP6248 T/F virus could also infect cell lines overexpressing the alternative coreceptors GPR15, APJ, and FPRL-1. A single mutation in the V3 crown sequence (GPEK->GPGK) of ZP6248 restored its infectivity in CCR5+ cells but reduced its ability to replicate in GPR15+ cells, indicating that the V3 crown motif played an important role in usage of this alternative coreceptor. These results suggest that the ZP6248 T/F virus established an acute in vivo infection by using coreceptor(s) other than CCR5 or CXCR4 or that the CCR5 coreceptor existed in an unusual conformation in this individual.

Evolution of drug-resistant viral populations during interruption of antiretroviral therapy

Analysis of a large number of HIV-1 genomes at multiple time points after antiretroviral treatment (ART) interruption allows determination of the evolution of drug-resistant viruses and viral fitness in vivo in the absence of drug selection pressure. Using a parallel allele-specific sequencing (PASS) assay, potential primary drug-resistant mutations in five individual patients were studied by analyzing over 18,000 viral genomes. A three-phase evolution of drug-resistant viruses was observed after termination of ART. In the first phase, viruses carrying various combinations of multiple-drug-resistant (MDR) mutations predominated with each mutation persisting in relatively stable proportions while the overall number of resistant viruses gradually increased. In the second phase, viruses with linked MDR mutations rapidly became undetectable and single-drug-resistant (SDR) viruses emerged as minority populations while wild-type viruses quickly predominated. In the third phase, low-frequency SDR viruses remained detectable as long as 59 weeks after treatment interruption. Mathematical modeling showed that the loss in relative fitness increased with the number of mutations in each viral genome and that viruses with MDR mutations had lower fitness than viruses with SDR mutations. No single viral genome had seven or more drug resistance mutations, suggesting that such severely mutated viruses were too unfit to be detected or that the resistance gain offered by the seventh mutation did not outweigh its contribution to the overall fitness loss of the virus. These data provide a more comprehensive understanding of evolution and fitness of drug-resistant viruses in vivo and may lead to improved treatment strategies for ART-experienced patients.