Switching between raltegravir resistance pathways analyzed by deep sequencing

Full-Length Envelope Analyzer (FLEA): A tool for longitudinal analysis of viral amplicons

Next generation sequencing of viral populations has advanced our understanding of viral population dynamics, the development of drug resistance, and escape from host immune responses. Many applications require complete gene sequences, which can be impossible to reconstruct from short reads. HIV env, the protein of interest for HIV vaccine studies, is exceptionally challenging for long-read sequencing and analysis due to its length, high substitution rate, and extensive indel variation. While long-read sequencing is attractive in this setting, the analysis of such data is not well handled by existing methods. To address this, we introduce FLEA (Full-Length Envelope Analyzer), which performs end-to-end analysis and visualization of long-read sequencing data. FLEA consists of both a pipeline (optionally run on a high-performance cluster), and a client-side web application that provides interactive results. The pipeline transforms FASTQ reads into high-quality consensus sequences (HQCSs) and uses them to build a codon-aware multiple sequence alignment. The resulting alignment is then used to infer phylogenies, selection pressure, and evolutionary dynamics. The web application provides publication-quality plots and interactive visualizations, including an annotated viral alignment browser, time series plots of evolutionary dynamics, visualizations of gene-wide selective pressures (such as dN/dS) across time and across protein structure, and a phylogenetic tree browser. We demonstrate how FLEA may be used to process Pacific Biosciences HIV env data and describe recent examples of its use. Simulations show how FLEA dramatically reduces the error rate of this sequencing platform, providing an accurate portrait of complex and variable HIV env populations. A public instance of FLEA is hosted at http://flea.datamonkey.org. The Python source code for the FLEA pipeline can be found at https://github.com/veg/flea-pipeline. The client-side application is available at https://github.com/veg/flea-web-app. A live demo of the P018 results can be found at http://flea.murrell.group/view/P018.

Viral populations constantly evolve and diversify. In this article we introduce a method, FLEA, for reconstructing and visualizing the details of evolutionary changes. FLEA specifically processes data from sequencing platforms that generate reads that are long, but error-prone. To study the evolutionary dynamics of entire genes during viral infection, data is collected via long-read sequencing at discrete time points, allowing us to understand how the virus changes over time. However, the experimental and sequencing process is imperfect, so the resulting data contain not only real evolutionary changes, but also mutations and other genetic artifacts caused by sequencing errors. Our method corrects most of these errors by combining thousands of erroneous sequences into a much smaller number of unique consensus sequences that represent biologically meaningful variation. The resulting high-quality sequences are used for further analysis, such as building an evolutionary tree that tracks and interprets the genetic changes in the viral population over time. FLEA is open source, and is freely available online.

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.

Drug Susceptibility and Viral Fitness of HIV-1 with Integrase Strand Transfer Inhibitor Resistance Substitution Q148R or N155H in Combination with Nucleoside/Nucleotide Reverse Transcriptase Inhibitor Resistance Substitutions

In clinical trials of coformulated elvitegravir (EVG), cobicistat (COBI), emtricitabine (FTC), and tenofovir disoproxil fumarate (TDF), emergent drug resistance predominantly involved the FTC resistance substitution M184V/I in reverse transcriptase (RT), with or without the tenofovir (TFV) resistance substitution K65R, accompanied by a primary EVG resistance substitution (E92Q, N155H, or Q148R) in integrase (IN). We previously reported that the RT-K65R, RT-M184V, and IN-E92Q substitutions lacked cross-class phenotypic resistance and replicative fitness compensation. As a follow-up, the in vitro characteristics of mutant HIV-1 containing RT-K65R and/or RT-M184V with IN-Q148R or IN-N155H were also evaluated, alone and in combination, for potential interactions. Single mutants displayed reduced susceptibility to their corresponding inhibitor classes, with no cross-class resistance. Viruses with IN-Q148R or IN-N155H exhibited reduced susceptibility to EVG (137- and 40-fold, respectively) that was not affected by the addition of RT-M184V or RT-K65R/M184V. All viruses containing RT-M184V were resistant to FTC (>1,000-fold). Mutants with RT-K65R had reduced susceptibility to TFV (3.3- to 3.6-fold). Without drugs present, the viral fitness of RT and/or IN mutants was diminished relative to that of the wild type in the following genotypic order: wild type > RT-M184V ≥ IN-N155H ≈ IN-Q148R ≥ RT-M184V + IN-N155H ≥ RT-M184V + IN-Q148R ≥ RT-K65R/M184V + IN-Q148R ≈ RT-K65R/M184V + IN-N155H. In the presence of drug concentrations approaching physiologic levels, drug resistance counteracted replication defects, allowing single mutants to outcompete the wild type with one drug present and double mutants to outcompete single mutants with two drugs present. These results suggest that during antiretroviral treatment with multiple drugs, the development of viruses with combinations of resistance substitutions may be favored despite diminished viral fitness.

HIV-1 integrase genotyping is reliable and reproducible for routine clinical detection of integrase resistance mutations even in patients with low-level viraemia

OBJECTIVES

Integrase drug resistance monitoring deserves attention because of the increasing number of patients being treated with integrase strand-transfer inhibitors. Therefore, we evaluated the integrase genotyping success rate at low-level viraemia (LLV, 51-1000 copies/mL) and resistance in raltegravir-failing patients.

METHODS

An integrase genotypic resistance test (GRT) was performed on 1734 HIV-1 samples collected during 2006-13. Genotyping success rate was determined according to the following viraemia levels: 51-500, 501-1000, 1001-10 000, 10 001-100 000 and >100 000 copies/mL. The reproducibility of integrase GRT was evaluated in 41 plasma samples processed in duplicate in two reference centres. The relationship between LLV and resistance prevalence was evaluated in a subset of 120 raltegravir-failing patients.

RESULTS

Overall, the integrase genotyping success rate was 95.7%. For viraemia levels 51-500 and 501-1000 copies/mL, the rate of success was 82.1% and 94.0%, respectively. GRT was reproducible, producing sequences with a high similarity and an equal resistance profile regardless of the sequencing centre or viraemia level. Resistance was detected both at LLV and at viraemia >1000 copies/mL (51-500 copies/mL = 18.2%; 501-1000 = 37.5%; 1001-10 000 = 53.7%; 10 001-100 000 = 30.0%; and >100 000 = 30.8%). At viraemia ≤500 copies/mL, Q148H/K/R and N155H had the same prevalence (9.1%), while the Y143C/H/R was completely absent. At early genotyping (within 3 months of raltegravir treatment), Q148H/K/R and N155H mutations were detected regardless of the viraemia level, while Y143C/H/R was observed only in samples with viraemia >1000 copies/mL.

CONCLUSIONS

Our findings prove the reliability of HIV-1 integrase genotyping and reinforce the concept that this assay may be useful in the management of failures even at LLV.

Adaptation on a genomic scale

Sequencing the genome of Candida albicans as it evolves in a patient reveals the genetic changes that allow the yeast to adapt to its environment.

Efficacy and safety of switching to raltegravir plus atazanavir dual therapy in pretreated HIV-1-infected patients over 144 weeks: a cohort study

BACKGROUND

To decrease drug burden among HIV-1-positive adults, we need a new gold standard for antiretroviral therapy maintenance strategies.

METHODS

This retrospective study aimed to assess efficacy in maintenance strategy of atazanavir (ATV) and raltegravir (RAL) dual therapy. The proportion of patients with HIV-1 RNA < 40 copies/ml at specific time points was recorded. Immunological response, safety, and pharmacokinetics were assessed.

RESULTS

Overall, 39 patients were switched to a RAL/ATV (n = 32) or RAL/ATV plus ritonavir (n = 7) regimen. Almost all patients (95%) received RAL twice daily. Most patients (70%) received a 400 mg ATV dosing per day, once (26%) or twice daily (44%). The percentages of virological success at weeks 24, 48, 96, and 144 were 92% (95% confidence interval (CI), 83-10), 86% (95% CI, 74-98), 70% (95% CI, 52-88), and 63% (95% CI, 42-84), respectively. Overall, 12 (31%) patients stopped dual therapy: 7 patients because of adverse events, mostly clinical myositis (n = 3). Confirmed virological failure occurred in three patients; two of them developed RAL resistance patterns. A significant increase in the CD4+/CD8 + T-cell ratio was observed at week 48 (p < 0.005). Only grade 1-2 adverse events were observed. Trough plasma levels presented a wide variability. Suggested trough concentrations were achieved in 79% and 94% of patients for ATV and RAL, respectively. An unboosted 400 mg per day ATV dosing seemed to be appropriate, regarding the targeted levels achieved and the lack of grade 3 or 4 hyperbilirubinemia.

CONCLUSIONS

We demonstrated, on a 3-year follow-up, the efficacy and safety of RAL plus ATV maintenance dual therapy.

Targeted deep sequencing of HIV-1 using the IonTorrentPGM platform

The characterization of mixed HIV-1 populations is a key question in clinical and basic research settings. This can be achieved through targeted deep sequencing (TDS), where next-generation sequencing is used to examine in depth a sub-genomic region of interest. This study explores the suitability of IonTorrent PGM(LifeTechnologies) for the TDS-based analysis of HIV-1 evolution. Using laboratory reagents and primary specimens sampled at pre-peak viremia the error rates from misincorporation and in vitro recombination were <0.5%. The sequencing error rate was 2- to 3-fold higher in/around homopolymeric tracts, and could be discerned from true polymorphism using bidirectional sequencing. The limit of detection of complex variants was further lowered by using haplotyping. The application of this system was illustrated on primary samples from an individual infected with HIV-1 followed from pre-peak viremia through six months post-acquisition. TDS provided an augmented view of the extent of genetic diversity, the covariation among polymorphisms, the evolutionary pathways, and the boundaries of the mutational space explored by the viral swarm. Based on its performance, the system can be applied for the characterization of minor viral variants in support of studies of viral evolution, which can inform the rational design of the next generation of vaccines and therapeutics.

High-Throughput Sequencing, a VersatileWeapon to Support Genome-Based Diagnosis in Infectious Diseases: Applications to Clinical Bacteriology

The recent progresses of high-throughput sequencing (HTS) technologies enable easy and cost-reduced access to whole genome sequencing (WGS) or re-sequencing. HTS associated with adapted, automatic and fast bioinformatics solutions for sequencing applications promises an accurate and timely identification and characterization of pathogenic agents. Many studies have demonstrated that data obtained from HTS analysis have allowed genome-based diagnosis, which has been consistent with phenotypic observations. These proofs of concept are probably the first steps toward the future of clinical microbiology. From concept to routine use, many parameters need to be considered to promote HTS as a powerful tool to help physicians and clinicians in microbiological investigations. This review highlights the milestones to be completed toward this purpose.

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.

Next-generation sequencing technology in clinical virology

Recent advances in nucleic acid sequencing technologies, referred to as 'next-generation' sequencing (NGS), have produced a true revolution and opened new perspectives for research and diagnostic applications, owing to the high speed and throughput of data generation. So far, NGS has been applied to metagenomics-based strategies for the discovery of novel viruses and the characterization of viral communities. Additional applications include whole viral genome sequencing, detection of viral genome variability, and the study of viral dynamics. These applications are particularly suitable for viruses such as human immunodeficiency virus, hepatitis B virus, and hepatitis C virus, whose error-prone replication machinery, combined with the high replication rate, results, in each infected individual, in the formation of many genetically related viral variants referred to as quasi-species. The viral quasi-species, in turn, represents the substrate for the selective pressure exerted by the immune system or by antiviral drugs. With traditional approaches, it is difficult to detect and quantify minority genomes present in viral quasi-species that, in fact, may have biological and clinical relevance. NGS provides, for each patient, a dataset of clonal sequences that is some order of magnitude higher than those obtained with conventional approaches. Hence, NGS is an extremely powerful tool with which to investigate previously inaccessible aspects of viral dynamics, such as the contribution of different viral reservoirs to replicating virus in the course of the natural history of the infection, co-receptor usage in minority viral populations harboured by different cell lineages, the dynamics of development of drug resistance, and the re-emergence of hidden genomes after treatment interruptions. The diagnostic application of NGS is just around the corner.

Multiple genetic pathways involving amino acid position 143 of HIV-1 integrase are preferentially associated with specific secondary amino acid substitutions and confer resistance to raltegravir and cross-resistance to elvitegravir

Y143C,R substitutions in HIV-1 integrase define one of three primary raltegravir (RAL) resistance pathways. Here we describe clinical isolates with alternative substitutions at position 143 (Y143A, Y143G, Y143H, and Y143S [Y143A,G,H,S]) that emerge less frequently, and we compare the genotypic and phenotypic profiles of these viruses to Y143C,R viruses to reconcile the preferential selection of Y143C,R variants during RAL treatment. Integrase amino acid sequences and RAL susceptibility were characterized in 117 patient isolates submitted for drug resistance testing and contained Y143 amino acid changes. The influence of specific Y143 substitutions on RAL susceptibility and their preferential association with particular secondary substitutions were further defined by evaluating the composition of patient virus populations along with a large panel of site-directed mutants. Our observations demonstrate that the RAL resistance profiles of Y143A,G,H,S viruses and their association with specific secondary substitutions are similar to the well-established Y143C profile but distinct from the Y143R profile. Y143R viruses differ from Y143A,C,G,H,S viruses in that Y143R confers a greater reduction in RAL susceptibility as a single substitution, consistent with a lower resistance barrier. Among Y143A,C,G,H,S viruses, the higher prevalence of Y143C viruses is the result of a lower genetic barrier than that of the Y143A,G,S viruses and a lower resistance barrier than that of the Y143H viruses. In addition, Y143A,C,G,H,S viruses require multiple secondary substitutions to develop large reductions in RAL susceptibility. Patient-derived viruses containing Y143 substitutions exhibit cross-resistance to elvitegravir.

Plasmodium falciparum-like parasites infecting wild apes in southern Cameroon do not represent a recurrent source of human malaria

Wild-living chimpanzees and gorillas harbor a multitude of Plasmodium species, including six of the subgenus Laverania, one of which served as the progenitor of Plasmodium falciparum. Despite the magnitude of this reservoir, it is unknown whether apes represent a source of human infections. Here, we used Plasmodium species-specific PCR, single-genome amplification, and 454 sequencing to screen humans from remote areas of southern Cameroon for ape Laverania infections. Among 1,402 blood samples, we found 1,000 to be Plasmodium mitochondrial DNA (mtDNA) positive, all of which contained human parasites as determined by sequencing and/or restriction enzyme digestion. To exclude low-abundance infections, we subjected 514 of these samples to 454 sequencing, targeting a region of the mtDNA genome that distinguishes ape from human Laverania species. Using algorithms specifically developed to differentiate rare Plasmodium variants from 454-sequencing error, we identified single and mixed-species infections with P. falciparum, Plasmodium malariae, and/or Plasmodium ovale. However, none of the human samples contained ape Laverania parasites, including the gorilla precursor of P. falciparum. To characterize further the diversity of P. falciparum in Cameroon, we used single-genome amplification to amplify 3.4-kb mtDNA fragments from 229 infected humans. Phylogenetic analysis identified 62 new variants, all of which clustered with extant P. falciparum, providing further evidence that P. falciparum emerged following a single gorilla-to-human transmission. Thus, unlike Plasmodium knowlesi-infected macaques in southeast Asia, African apes harboring Laverania parasites do not seem to serve as a recurrent source of human malaria, a finding of import to ongoing control and eradication measures.

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.

Impact of primary elvitegravir resistance-associated mutations in HIV-1 integrase on drug susceptibility and viral replication fitness

Elvitegravir (EVG) is an effective HIV-1 integrase (IN) strand transfer inhibitor (INSTI) in advanced clinical development. Primary INSTI resistance-associated mutations (RAMs) at six IN positions have been identified in HIV-1-infected patients failing EVG-containing regimens in clinical studies: T66I/A/K, E92Q/G, T97A, S147G, Q148R/H/K, and N155H. In this study, the effect of these primary IN mutations, alone and in combination, on susceptibility to the INSTIs EVG, raltegravir (RAL), and dolutegravir (DTG); IN enzyme activities; and viral replication fitness was characterized. Recombinant viruses containing the six most common mutations exhibited a range of reduced EVG susceptibility: 92-fold for Q148R, 30-fold for N155H, 26-fold for E92Q, 10-fold for T66I, 4-fold for S147G, and 2-fold for T97A. Less commonly observed primary IN mutations also showed a range of reduced EVG susceptibilities: 40- to 94-fold for T66K and Q148K and 5- to 10-fold for T66A, E92G, and Q148H. Some primary IN mutations exhibited broad cross-resistance between EVG and RAL (T66K, E92Q, Q148R/H/K, and N155H), while others retained susceptibility to RAL (T66I/A, E92G, T97A, and S147G). Dual combinations of primary IN mutations further reduced INSTI susceptibility, replication capacity, and viral fitness relative to either mutation alone. Susceptibility to DTG was retained by single primary IN mutations but reduced by dual mutation combinations with Q148R. Primary EVG RAMs also diminished IN enzymatic activities, concordant with their structural proximity to the active site. Greater reductions in viral fitness of dual mutation combinations may explain why some primary INSTI RAMs do not readily coexist on the same HIV-1 genome but rather establish independent pathways of resistance to EVG.

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.

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.

Prolonged and substantial discordance in prevalence of raltegravir-resistant HIV-1 in plasma versus PBMC samples revealed by 454 "deep" sequencing

The evolution of drug resistance mutations in plasma samples is relatively well-characterized. However, the viral population and diversity in other body compartments such as peripheral blood mononuclear cells (PBMC) remains poorly understood. Previous studies have mostly focused on protease and reverse transcriptase drug resistance mutations (DRMs). In this study, we used 454 "deep" sequencing technology to observe and quantify longitudinally the prevalence of resistance mutations associated with the integrase inhibitor, raltegravir, in plasma versus PBMC samples from a San Francisco-based cohort. Four heavily treatment-experienced subjects were monitored in this study over a median of 1.2 years since the initiation of raltegravir-containing regimens. We observed a consistent discordance in the prevalence of DRMs, but not resistance pathway(s), in the plasma versus PBMC viral populations. In the final paired samples that were tested while the subjects were on a raltegravir-containing regimen, DRM prevalence reached 100% in plasma but remained 1% in PBMC on day 177 post-therapy in Subject 3180 (Q148H/G140S), 100% in plasma and 36% in PBMC on day 224 in Subject 3242 (N155H), 78% in plasma and 11-12% in PBMC on day 338 in Subject 3501 (Q148H/G140S), and 100% in plasma and 0% in PBMC on day 197 in Subject 3508 (Y143R). Furthermore, absolute sequence homology comparison between the two compartments revealed that 21% - 99% of PBMC sequences had no match in plasma, whereas 14% - 100% of plasma sequences had no match in PBMC. Overall, our observations suggested that plasma and PBMC hosted drastically different HIV-1 populations even after a prolonged exposure to raltegravir selection pressure.

Longitudinal analysis of integrase N155H variants in heavily treated patients failing raltegravir-based regimens

OBJECTIVES

The mechanism of raltegravir (RAL)-resistant evolutions has not already been elucidated. Because the emergence of RAL resistance is usually initiated by the N155H mutant, we assessed the role of minor N155H-mutated variants in circulating RNA and archived DNA in five heavily treated patients experiencing long-term RAL therapy failure and harbouring three different resistance profiles determined by standard genotyping.

METHODS

Allele-specific polymerase chain reaction (AS-PCR) was used to detect N155H mutants in longitudinal stored plasma and whole-blood samples before, during and after RAL-based regimens in five patients infected with the HIV-1 B subtype.

RESULTS

No minor N155H-mutated variant was found by AS-PCR in either plasma or whole-blood samples collected at baseline and after RAL withdrawal in any of the five patients. During RAL failure, the mutation N155H was detected at different levels in three patients displaying the N155H pathway and gradually declined when the double mutant Q148H+G140S was selected in one patient. In two patients with the Q148H resistance pathway, no N155H variant was identified by AS-PCR in either viral RNA or DNA.

CONCLUSIONS

The N155H mutation present at various levels from minority to majority showed no relationship with the three RAL-associated resistance profiles, suggesting that this mutant may not play a role in determining different resistance profiles. Moreover, pre-existing N155H is very infrequent and, if selected during RAL failure, the N155H mutant disappears quickly after RAL withdrawal.

Challenges and opportunities in estimating viral genetic diversity from next-generation sequencing data

Many viruses, including the clinically relevant RNA viruses HIV (human immunodeficiency virus) and HCV (hepatitis C virus), exist in large populations and display high genetic heterogeneity within and between infected hosts. Assessing intra-patient viral genetic diversity is essential for understanding the evolutionary dynamics of viruses, for designing effective vaccines, and for the success of antiviral therapy. Next-generation sequencing (NGS) technologies allow the rapid and cost-effective acquisition of thousands to millions of short DNA sequences from a single sample. However, this approach entails several challenges in experimental design and computational data analysis. Here, we review the entire process of inferring viral diversity from sample collection to computing measures of genetic diversity. We discuss sample preparation, including reverse transcription and amplification, and the effect of experimental conditions on diversity estimates due to in vitro base substitutions, insertions, deletions, and recombination. The use of different NGS platforms and their sequencing error profiles are compared in the context of various applications of diversity estimation, ranging from the detection of single nucleotide variants (SNVs) to the reconstruction of whole-genome haplotypes. We describe the statistical and computational challenges arising from these technical artifacts, and we review existing approaches, including available software, for their solution. Finally, we discuss open problems, and highlight successful biomedical applications and potential future clinical use of NGS to estimate viral diversity.

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.

Viral quasispecies evolution

Evolution of RNA viruses occurs through disequilibria of collections of closely related mutant spectra or mutant clouds termed viral quasispecies. Here we review the origin of the quasispecies concept and some biological implications of quasispecies dynamics. Two main aspects are addressed: (i) mutant clouds as reservoirs of phenotypic variants for virus adaptability and (ii) the internal interactions that are established within mutant spectra that render a virus ensemble the unit of selection. The understanding of viruses as quasispecies has led to new antiviral designs, such as lethal mutagenesis, whose aim is to drive viruses toward low fitness values with limited chances of fitness recovery. The impact of quasispecies for three salient human pathogens, human immunodeficiency virus and the hepatitis B and C viruses, is reviewed, with emphasis on antiviral treatment strategies. Finally, extensions of quasispecies to nonviral systems are briefly mentioned to emphasize the broad applicability of quasispecies theory.

Substitutions at amino acid positions 143, 148, and 155 of HIV-1 integrase define distinct genetic barriers to raltegravir resistance in vivo

Mutations at amino acids 143, 148, and 155 in HIV-1 integrase (IN) define primary resistance pathways in subjects failing raltegravir (RAL)-containing treatments. Although each pathway appears to be genetically distinct, shifts in the predominant resistant virus population have been reported under continued drug pressure. To better understand this dynamic, we characterized the RAL susceptibility of 200 resistant viruses, and we performed sequential clonal analysis for selected cases. Patient viruses containing Y143R, Q148R, or Q148H mutations consistently exhibited larger reductions in RAL susceptibility than patient viruses containing N155H mutations. Sequential analyses of virus populations from three subjects revealed temporal shifts in subpopulations representing N155H, Y143R, or Q148H escape pathways. Evaluation of molecular clones isolated from different time points demonstrated that Y143R and Q148H variants exhibited larger reductions in RAL susceptibility and higher IN-mediated replication capacity (RC) than N155H variants within the same subject. Furthermore, shifts from the N155H pathway to either the Q148R or H pathway or the Y143R pathway were dependent on the amino acid substitution at position 148 and the secondary mutations in Y143R- or Q148R- or H-containing variants and correlated with reductions in RAL susceptibility and restorations in RC. Our observations in patient viruses were confirmed by analyzing site-directed mutations. In summary, viruses that acquire mutations defining the 143 or 148 escape pathways are less susceptible to RAL and exhibit greater RC than viruses containing 155 pathway mutations. These selective pressures result in the displacement of N155H variants by 143 or 148 variants under continued drug exposure.

Raltegravir in HIV-1 infection: Safety and Efficacy in Treatment-naïve Patients

The hunt for a compound which inhibits the HIV-1 integrase had been painstakingly difficult. Integrase is essential for viral replication as it mediates the integration of the viral DNA genome into the host DNA resulting in the establishment of the permanent provirus. Persistent efforts have resulted in the discovery of Raltegravir (Isentress, MK-0518), the first integrase inhibitor approved by US Food and Drug Administration for the treatment in HIV-1 infected patients. Numerous clinical studies with raltegravir have found it to be safe and effective in treatment naïve as well as treatment experienced patients. Adverse events associated with raltegravir based therapy are milder compared to previously available regimens. Raltegravir is metabolized primarily via glucuronidation mediated by uridine diphosphate glucuronosyltransferase and has a favorable pharmacokinetics independent of age, gender, race, food, and drug-drug interactions. Within a short period of time of its introduction, raltegravir has been included as one of DHHS recommended preferred regimen for the treatment of HIV-1 infection in treatment naïve patients.

Applications of next-generation sequencing technologies to diagnostic virology

Novel DNA sequencing techniques, referred to as “next-generation” sequencing (NGS), provide high speed and throughput that can produce an enormous volume of sequences with many possible applications in research and diagnostic settings. In this article, we provide an overview of the many applications of NGS in diagnostic virology. NGS techniques have been used for high-throughput whole viral genome sequencing, such as sequencing of new influenza viruses, for detection of viral genome variability and evolution within the host, such as investigation of human immunodeficiency virus and human hepatitis C virus quasispecies, and monitoring of low-abundance antiviral drug-resistance mutations. NGS techniques have been applied to metagenomics-based strategies for the detection of unexpected disease-associated viruses and for the discovery of novel human viruses, including cancer-related viruses. Finally, the human virome in healthy and disease conditions has been described by NGS-based metagenomics.