Open Access

Stable HIV-1 integrase diversity during initial HIV-1 RNA Decay suggests complete blockade of plasma HIV-1 replication by effective raltegravir-containing salvage therapy

  • Marc Noguera-Julian1Email author,
  • Maria Casadellà1,
  • Christian Pou1,
  • Cristina Rodríguez1,
  • Susana Pérez-Álvarez1,
  • Jordi Puig2,
  • Bonaventura Clotet1, 2, 3 and
  • Roger Paredes1, 2, 3
Virology Journal201310:350

https://doi.org/10.1186/1743-422X-10-350

Received: 1 October 2013

Accepted: 27 November 2013

Published: 5 December 2013

Abstract

Background

There is legitimate concern that minority drug-resistant mutants may be selected during the initial HIV-1 RNA decay phase following antiretroviral therapy initiation, thus undermining efficacy of treatment. The goal of this study was to characterize viral resistance emergence and address viral population evolution during the first phase of viral decay after treatment containing initiation.

Findings

454 sequencing was used to characterize viral genetic diversity and polymorphism composition of the HIV-1 integrase gene during the first two weeks following initiation of raltegravir-containing HAART in four ART-experienced subjects. No low-prevalence Raltegravir (RAL) drug resistance mutations (DRM) were found at baseline. All patients undergoing treatment received a fully active ART according to GSS values (GSS ≥ 3.5). No emergence of DRM after treatment initiation was detected. Longitudinal analysis showed no evidence of any other polymorphic mutation emergence or variation in viral diversity indexes.

Conclusions

This suggests that fully active salvage antiretroviral therapy including raltegravir achieves a complete blockade of HIV-1 replication in plasma. It is unlikely that raltegravir-resistant HIV-1 may be selected in plasma during the early HIV-1 RNA decay after treatment initiation if the administered therapy is active enough.

Keywords

HIV-1RaltegravirDynamicsDeep sequencing

Findings

There is legitimate concern that minority drug-resistant mutants may be selected during the initial HIV-1 RNA decay phase following antiretroviral therapy initiation, when decreasing but still detectable HIV-1 RNA levels transiently coexist with suboptimal drug levels. Studies using allele-specific PCR suggested that drug-resistant HIV-1 might be selected and replicate in the first weeks of suppressive antiretroviral therapy (ART) [1, 2]. If they were confirmed, such observations would support using intensified ART regimens to accelerate viremia suppression and curtail drug resistance evolution, but not only in treatment-naïve subjects initiating first-line ART. Achieving fast blockade of HIV replication would even be more important in subjects with drug-resistant HIV-1 initiating salvage ART, given that they are more vulnerable to drug resistance evolution and are closer to exhausting their treatment options.

In the past, however, first-line treatment with 4, 5, or more antiretrovirals did not produce better virological outcomes than standard 3-drug ART [3]. Indeed, evolutionary studies in macaque models [4] suggested that viral replication was completely shut down by therapy and that development of antiretroviral drug resistance during the initial viremia decay was unlikely to occur, provided that treatments were active enough and resistant variants did not pre-exist. If HIV replication could be shut down with standard ART, addition of further drugs to ART regimens would not provide additional protection against drug resistance.

Mathematical models of HIV resistance kinetics predict faster selection of resistant viruses with more potent regimens. On the other hand, previous studies correlated higher resistance selection with longer time to achieving undetectable HIV-1 RNA levels [1, 2]. raltegravir (RAL) achieves fast viremia suppression when combined with an optimized backbone treatment (OBT) [5], but has a low-genetic barrier and HIV-1 variants with single raltegravir drug resistance mutations could pre-exist at low levels [6, 7].

We thus sought to explore in a proof-of-concept study whether raltegravir-resistant mutants were selected during the initial HIV-1 RNA decay following initiation of raltegravir-containing ART in 4 treatment-experienced subjects with HIV-1 resistance to at least 2 antiretroviral drug classes using ultrasensitive, quantitative 454 sequencing. We also sought to characterize the diversity and viral population shifts in HIV-1 integrase during this initial HIV-RNA decay phase.

Study design

The study included 4 ART-experienced HIV-1-infected subjects who initiated salvage ART including raltegravir (RAL) with more than 3 active drugs (Table 1). All subjects provided signed informed consent to participate into the study. The Ethics Committee of the “Germans Trias i Pujol” Hospital, approved the study on 21 December 2007, reference #: AC-07-107. Integrase sequencing was performed using 454 Genome Sequencer FLX (Life Sequencing/Roche) in consecutive plasma samples with detectable viral loads on days 0 to 14 following RAL initiation. 454 sequencing of integrase was performed using a four-amplicon design for robustness (see Additional file 1). A homology filter was applied to rule out pNL43 contamination. AVA v2.8 software was used to obtain mutation prevalence and unique consensus sequences, which were further processed. Short consensus sequences (less than 90% of the expected length) were filtered out. Aminoacid variant calling was realized using in-house perl code on pairwise alignments over HXB2R reference sequence, discarding those sequences with in-frame STOP codons. In order to discard strand-dependent sequencing errors, only variants which were present both in the forward and reverse strand and presented a forward/reverse prevalence frequency ratio within 1 log were accepted. GSS scores were calculated using Stanford HIVdb guidelines. MAFFT (v6.815b) software [8] was then used to perform multiple alignment on consensus sequences. Shannon entropy was calculated for each amplicon using the alignment obtained and the prevalence information kept for each position. Mean pairwise distances (both intra amplicon and relative to an HXB2 reference) were calculated using a Kimura-80 model and sequence prevalence information with R/ape package [9].
Table 1

Subject clinical details

Subject

DS1

DS2

DS3

DS4

Gender

Male

Male

Male

Female

Age

61

39

49

44

Date of diagnostic

13/05/96

01/01/90

01/07/87

01/01/91

Transmission type

Heterosexual

IVDU

IVDU

IVDU

Previous ART lines

9

6

12

12

#drugs before RAL

13

9

13

12

New drugs started alongside RAL

DRV

DRV

DRV,ETV, T-20

DRV,MVC

Previous exposure

    

NRTI

YES

YES

YES

YES

NNRTI

YES

YES

YES

YES

PI

YES

YES

YES

YES

T-20

NO

NO

NO

NO

ART

ddI + 3TC + DRV/r + RAL

TDF/FTC + DRV/r + RAL

TDF/FTC + DRV/r + RAL + T-20

DRV/r + RAL + MVC

GSS

4

3.5

3.5

3.5

Baseline CD4 (cells/mm 3)

614

300

332

399

Intravenous Drug Use (IVDU); Genotypic Sensitivity Score (GSS, Stanford HIVdb), NRTI: Nucleoside-analogues Reverse Transcriptase Inhibitor; NNRTI non- Nucleoside-analogues Reverse Transcriptase Inhibitor; PI: Protease Inhibitor; Darunavir (DRV); Didanosine (DDI); Lamivudine (3TC); Raltregravir (RAL); Ritonavir (RTV); Emtricitabine (FTC); Tenofovir (TDF); Etravirine (ETV); Enfuvirtide (T-20); Maraviroc (MVC).

Results

The genotypic susceptibility score (GSS) of the salvage treatment was ≥3.5 in all patients according to population sequencing. Sensitivity scores were further confirmed considering 454 deep sequencing data obtained for integrase and protease/retrotranscriptase, when available (see Additional file 2: Table S1). At baseline, median (interquartile range) for HIV-1 ARN and CD4+ T cell counts were, 116.500 (57.750; 215.000) copies/mL and 399 (300;614) cells/mm3. Viral load decline at day 7, when available, was more than 1log copies/mL.

The median number of sequences obtained from 454 sequencing per amplicon was 3404 (IQR: 1497–5304). Initiation of treatment with raltegravir was associated with a rapid decline in HIV-1 RNA levels but no changes in integrase diversity or shifts relative to an HXB2R external reference (see Additional file 3: Figures S1, S2 and Additional file 2: Table S2). The frequency of major integrase polymorphisms remained stable during the viral load decay phase (Figure 1). In addition, longitudinal 454 sequencing of integrase gene did not result in the emergence of resistant mutations or any of the detected polymorphisms and, so, sensitivity to RAL remained unchanged during the first two weeks and none of the minor viral populations carrying minor polymorphisms were selected.
Figure 1

Longitudinal evolution of integrase polymorphisms relative to Consensus B. Each plot corresponds to one subject. Integrase polymorphisms are in the vertical axis. The horizontal axis shows the different longitudinal timepoints assessed during the initial HIV-1 RNA load decay. The frequency of each polymorphism is represented using a color scale, from 100% (dark red) to a 1.0% threshold (light blue). White indicates non detection of polymorphisms. A remarkable stability in polymorphism type and frequency is observed. Of note, none of the IAS-list (2013) integrase resistant mutations were found.

Discussion

We observed stability of HIV-1 integrase diversity and heterogeneity within the first phase of HIV-1 RNA decay. The type and frequency of polymorphisms in sequential samples during early HIV-1 RNA decay suggests that salvage antiretroviral therapy including raltegravir is able to achieve an effective blockade of HIV-1 replication in plasma, even in heavily treatment-experienced individuals, in a context where a potent antiretroviral regimen, defined by a GSS ≥ 3, can be constructed. Of note, the rapid decline on HIV-1 RNA levels may also affect the capability of detecting minor polymorphisms using 454. When the number of obtained sequences is higher than the number of expected initial viral templates, re-sequencing of the same minor variants may lead to an artifactual appearance of low prevalence polymorphisms. Interestingly, minority polymorphisms were more frequent and fluctuated at later time points, at very low levels. Thus, obtained minor variant frequencies must be interpreted accounting for sampling limitations. Conversely, sampling limitations are not expected to significantly affect viral diversity due to the much larger impact of highly prevalent variants over minority ones in these measures.

Importantly, the blockade in replication observed in this study must be understood as a result of the global antiretroviral treatment pressure against the whole viral quasi-species. Indeed, additional 454 data on protease and reverse transcriptase at baseline confirms the full potency of the prescribed regimen, as GSS scores were not modified by the few minor resistant variants additionally detected. Our findings contrast with the rapid emergence of raltegravir resistance when non- fully active regimens (GSS < 3) are used [6, 7]. Also, antiretroviral drug resistance might evolve in other compartments or reservoirs with low drug penetration in the presence of residual HIV-1 replication. It is also worth noting that the sample size in this study is low (n = 4), thus reducing our ability to observe emergence of resistance, and that GSS in our study ranged from 3.5 to 4, so we cannot conclude whether a GSS = 3 would achieve similar results.

In conclusion, raltegravir-resistant HIV-1 may not be selected in plasma during the early HIV-1 RNA decay after salvage treatment initiation, provided that subjects are given a fully effective regimen that is able to shut down plasma replication. Further studies are required to assess whether similar findings might also be observed with other drugs or regimens and to clarify what is the minimum GSS (i.e. 3, 4 or higher) required to block HIV-1 evolution during HIV-1 RNA decay following ART initiation.

Declarations

Acknowledgements

This study was supported through an unrestricted research grant from Merck Sharpe & Dohme, ‘CHAIN, Collaborative HIV and Anti-HIV Drug Resistance Network’, Integrated Project no. 223131, funded by the European Commission Framework 7 Program, the Spanish AIDS network ‘Red Temática Cooperativa de Investigación en SIDA’ (RD06/0006), and the ‘Gala contra la sida – Barcelona 2011’. Work Contained in this paper was presented in the 19th Conference on Retroviruses and Opportunistic Infections, March, 2012, Seattle, WA, USA.

Authors’ Affiliations

(1)
IrsiCaixa AIDS Research Institute-HIVACAT, Hospital Universitari Germans Trias I Pujol
(2)
Fundació Lluita contra la SIDA Hospital Universitari Germans Trias i Pujol
(3)
Universitat Autònoma de Barcelona

References

  1. Metzner KJ, Allers K, Rauch P, Harrer T: Rapid selection of drug-resistant HIV-1 during the first months of suppressive ART in treatment-naive patients. AIDS 2007, 21: 703-711. 10.1097/QAD.0b013e3280121ac6PubMedView ArticleGoogle Scholar
  2. Metzner KJ, et al.: Minority quasispecies of drug-resistant HIV-1 that lead to early therapy failure in treatment-naive and -adherent patients. Clin Infect Dis 2009, 48: 239-247. 10.1086/595703PubMedView ArticleGoogle Scholar
  3. Hirsch MS, et al.: Antiretroviral drug resistance testing in adult HIV-1 infection: 2008 recommendations of an International AIDS Society-USA panel. Clin Infect Dis 2008, 47: 266-285. 10.1086/589297PubMedView ArticleGoogle Scholar
  4. Kearney M, et al.: Genetic diversity of simian immunodeficiency virus encoding HIV-1 reverse transcriptase persists in macaques despite antiretroviral therapy. J Virol 2011, 85: 1067-1076. 10.1128/JVI.01701-10PubMedPubMed CentralView ArticleGoogle Scholar
  5. Markowitz M, et al.: Rapid and durable antiretroviral effect of the HIV-1 Integrase inhibitor raltegravir as part of combination therapy in treatment-naive patients with HIV-1 infection: results of a 48-week controlled study. J Acquir Immune Defic Syndr 2007, 46: 125-133. 10.1097/QAI.0b013e318157131cPubMedView ArticleGoogle Scholar
  6. Armenia D, et al.: Study of genotypic and phenotypic HIV-1 dynamics of integrase mutations during raltegravir treatment: a refined analysis by ultra-deep 454 pyrosequencing. J Infect Dis 2012, 205: 557-567. 10.1093/infdis/jir821PubMedPubMed CentralView ArticleGoogle Scholar
  7. Codoñer FM, et al.: Dynamic escape of pre-existing raltegravir-resistant HIV-1 from raltegravir selection pressure. Antiviral Res 2010, 88: 281-286. 10.1016/j.antiviral.2010.09.016PubMedView ArticleGoogle Scholar
  8. Katoh K, Standley DM: MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 2013, 30: 772-780. 10.1093/molbev/mst010PubMedPubMed CentralView ArticleGoogle Scholar
  9. Paradis E, Claude J, Strimmer K: APE: analyses of phylogenetics and evolution in R language. Bioinformatics 2004, 20: 289-290. 10.1093/bioinformatics/btg412PubMedView ArticleGoogle Scholar

Copyright

© Noguera-Julian et al.; licensee BioMed Central Ltd. 2013

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Advertisement