Abstract

Objectives

To further characterize HIV-1 viruses of patients experiencing unexplained virological failure (VF) on PI-containing regimens, ultradeep sequencing was performed on protease, gag and gp41 genes in patients failing a first-line treatment.

Methods

All naive patients initiating an antiretroviral treatment based on boosted darunavir, atazanavir or lopinavir and experiencing VF without any transmitted drug resistance mutation detected by Sanger sequencing on protease and reverse transcriptase genes were selected. Ultradeep sequencing (IlluminaTM Nextera®) was performed on protease, gag and gp41 genes in plasma before initiation of treatment and at VF to identify emergent mutations.

Results

Among the 32 patients included in the study, emergent and previously undescribed mutations in the viral protease gene were identified in five patients at VF: 64M (1 CRF02_AG), 64M/70R with mutation 15V (2 CRF02_AG), 79A (1 CRF06_cpx) and 79A with mutation 15V (1 CRF02_AG). Two patients showed the emergence of R286K in the gag region, outside of cleavage sites (2 CRF02_AG). In the gp41 region, the V321I mutation emerged inside the cytoplasmic tail (1 subtype A and 1 subtype B). All these patients were treated with a darunavir/ritonavir-based regimen.

Conclusions

In some cases of VF to PIs, we observed the emergence of protease, Gag or Gp41 mutations that had not previously been associated with VF or PI resistance. These mutations should be further studied, in particular the 15V/64M/70R pattern in the protease gene identified among CRF02_AG viruses.

Introduction

Selection of drug-resistant HIV during viral replication under treatment pressure is a major reason for therapy failure. While the introduction of PIs has ensured robust combination therapy, virological failure (VF) still occurs in clinical practice. PIs are known to suppress proteolytic cleavage of HIV Gag and Pol polyproteins that include essential structural and enzymatic components of the virus.1 This prevents the conversion of HIV particles into their mature infectious form. Major PI resistance mutations in the protease gene lead to decreased PI binding.2–5 In parallel, compensatory mutations outside the substrate-binding pocket of protease can at least partially restore viral replication.6,7

The viruses of most patients on PIs who experience VF do not harbour any resistance associated-mutations in the protease gene. Studies of mutations inside or outside gag gene cleavage sites (CSs) have contributed to the understanding of new PI resistance mechanisms,8–12 but the impact of some identified mutations has not been clearly defined.13 Moreover, it has been suggested that Gp41 mutations could also have implications for resistance to PIs. Indeed, approximately half the inhibitory potential of PIs occurs at the entry step, probably reflecting the interactions between the uncleaved Gag polyprotein and the cytoplasmic tail (CT) of the envelope transmembrane glycoprotein Gp41.14

To further understanding of VF in patients failing a first-line PI-based regimen, we used ultradeep sequencing (UDS) to study the protease, gag and gp41 regions.

Patients and methods

Patients

All naive patients initiating an antiretroviral treatment based on a boosted darunavir-, atazanavir- or lopinavir-containing regimen and experiencing VF without any transmitted drug resistance (TDR) mutation detected by Sanger sequencing at baseline on protease and reverse transcriptase genes according to the 2009 consensus TDR list were selected.15 VF was defined as two consecutive plasma viral loads >50 copies/mL after a plasma viral load <50 copies/mL or as a viral load decrease not satisfying French monitoring guidelines after treatment initiation.16 Plasma concentrations of all ART drugs were determined using UPLC coupled with tandem mass spectrometry (Waters Corp., Milford, MA, USA).

Patients were monitored by the Department of Infectious Diseases at Pitié-Salpêtrière (from January 2011 to January 2016) and Bichat-Claude Bernard hospitals (from January 2012 to March 2015) (Paris, France). Data were obtained from the existing electronic database or medical records. Patients were informed that their demographic and clinical data would be recorded during follow-up and could be used for retrospective studies, and gave their consent. This study was approved by the Agence Nationale de Recherches sur le SIDA et les hépatites virales (ANRS) AC11 Ethics Committee.

UDS

UDS of protease, gag and gp41 genes was performed in plasma before initiation of treatment and at VF. After RNA extraction (Qiagen® QIAamp viral RNA mini, Venlo, The Netherlands), a reverse transcription into DNA was performed (Transcriptor®, Roche, Basel, Switzerland) and the DNA was then amplified by PCR. A nested PCR was carried out with a high-fidelity Taq polymerase (Q5 High-Fidelity DNA Polymerase, New England Biolabs, Ipswich, MA, USA for the protease gene; Superscript III® one-step RT-PCR with platinum Taq high fidelity, Invitrogen, Carlsbad, CA, USA for the gag and gp41 genes). UDS was in accordance with the manufacturer’s recommendations (Nextera® kit, Illumina®, San Diego, CA, USA) with a cut-off of 1% at baseline. At failure, only emergent mutations present in more than 20% of the total viral population, detected in at least one of two samples and whose proportion increased at least 2.5-fold between treatment initiation and VF were taken into account. Geneious® 9.0 (Brisbane, Australia) was used to analyse UDS results.17

Results

During the period under study, 409 patients began antiretroviral treatment with reverse transcriptase inhibitors in combination with boosted darunavir (318 patients), atazanavir (72 patients) or lopinavir (19 patients).

Of the total number of patients, 32 experienced VF without any TDR mutation detected by Sanger sequencing of protease and reverse transcriptase genes (7.8%) and were retrospectively included in the study. Among this group, 53.1% were infected with HIV-1 recombinant CRF02_AG. The median time between baseline and VF was 495 days and the median nadir CD4 cell count was 220 cells/mm3. Antiretroviral plasma concentrations in samples taken at VF were consistent with adequate compliance for all patients. Patient characteristics and treatments received are presented in Table 1.

Table 1.

Characteristics of patients and treatments received

CharacteristicsTotal population (n = 32)
Age, years, median (IQR)45 (41–50)
Male sex (%)57.6
Subtype B (%)19.4
CRF02_AG (%)53.1
CD4 cell count at VF, cells/mm3, median (IQR)243 (103–390)
Nadir CD4 cell count, cells/mm3, median (IQR)220 (99–319)
Plasma HIV viral load at D0, log10 copies/mL, median (IQR)5.3 (5–5.5)
Plasma HIV viral load at VF, log10 copies/mL, median (IQR)2.2 (2–2.7)
Time to treatment failure, days, median (IQR)495 (366–735)
Treatments, number of patients
 2 NRTIs+PI/r31
  tenofovir/emtricitabine+PI/r27
 1 INSTI+PI/r1
Protease inhibitors, number of patients
 darunavir25 (20 q24h; 5 q12h)
 atazanavir6
 lopinavir1
CharacteristicsTotal population (n = 32)
Age, years, median (IQR)45 (41–50)
Male sex (%)57.6
Subtype B (%)19.4
CRF02_AG (%)53.1
CD4 cell count at VF, cells/mm3, median (IQR)243 (103–390)
Nadir CD4 cell count, cells/mm3, median (IQR)220 (99–319)
Plasma HIV viral load at D0, log10 copies/mL, median (IQR)5.3 (5–5.5)
Plasma HIV viral load at VF, log10 copies/mL, median (IQR)2.2 (2–2.7)
Time to treatment failure, days, median (IQR)495 (366–735)
Treatments, number of patients
 2 NRTIs+PI/r31
  tenofovir/emtricitabine+PI/r27
 1 INSTI+PI/r1
Protease inhibitors, number of patients
 darunavir25 (20 q24h; 5 q12h)
 atazanavir6
 lopinavir1

D0, day 0 (baseline); INSTI, integrase inhibitor; VF, virological failure.

Table 1.

Characteristics of patients and treatments received

CharacteristicsTotal population (n = 32)
Age, years, median (IQR)45 (41–50)
Male sex (%)57.6
Subtype B (%)19.4
CRF02_AG (%)53.1
CD4 cell count at VF, cells/mm3, median (IQR)243 (103–390)
Nadir CD4 cell count, cells/mm3, median (IQR)220 (99–319)
Plasma HIV viral load at D0, log10 copies/mL, median (IQR)5.3 (5–5.5)
Plasma HIV viral load at VF, log10 copies/mL, median (IQR)2.2 (2–2.7)
Time to treatment failure, days, median (IQR)495 (366–735)
Treatments, number of patients
 2 NRTIs+PI/r31
  tenofovir/emtricitabine+PI/r27
 1 INSTI+PI/r1
Protease inhibitors, number of patients
 darunavir25 (20 q24h; 5 q12h)
 atazanavir6
 lopinavir1
CharacteristicsTotal population (n = 32)
Age, years, median (IQR)45 (41–50)
Male sex (%)57.6
Subtype B (%)19.4
CRF02_AG (%)53.1
CD4 cell count at VF, cells/mm3, median (IQR)243 (103–390)
Nadir CD4 cell count, cells/mm3, median (IQR)220 (99–319)
Plasma HIV viral load at D0, log10 copies/mL, median (IQR)5.3 (5–5.5)
Plasma HIV viral load at VF, log10 copies/mL, median (IQR)2.2 (2–2.7)
Time to treatment failure, days, median (IQR)495 (366–735)
Treatments, number of patients
 2 NRTIs+PI/r31
  tenofovir/emtricitabine+PI/r27
 1 INSTI+PI/r1
Protease inhibitors, number of patients
 darunavir25 (20 q24h; 5 q12h)
 atazanavir6
 lopinavir1

D0, day 0 (baseline); INSTI, integrase inhibitor; VF, virological failure.

Before initiation of treatment, samples from 26, 27 and 27 patients were analysed for protease, gag and gp41 genes, respectively. At failure, samples from 21, 13 and 14 patients were analysed for protease, gag and gp41 genes, respectively. Before initiation of treatment, the median number of reads was 118443, 195405 and 197368 for protease, gag and gp41 sequences, respectively. At VF the median number of reads was 289786, 360697 and 316552 for protease, gag and gp41 sequences, respectively.

At VF, five of the 32 patients displayed viruses harbouring emergent mutations in the protease gene that have not been previously described at VF: 64M (one patient infected with a CRF02_AG subtype), 64M/70R in combination with mutation 15V (two patients infected with a CRF02_AG subtype), 79A (one patient infected with a CRF06_cpx subtype), and 79A in combination with mutation 15V (one patient infected with a CRF02_AG subtype). All five patients were treated with a combination of two NRTIs and boosted darunavir q24h (Table 2).

Table 2.

Resistance mutations detected at baseline and mutations selected at virological failure (ultradeep sequencing) in protease, gag and gp41 genes

PatientSubtypeARTTimeProtease (%)Gag (%)Gp41 (%)
1CRF02_AGDRV q24hD0S373P (97.7%), L449P (99%)
VFL449P (99.6%)
2CRF02_AGDRV q24hD0K436R (98.5%), L449P (98.7%)
VFK436R (99.1%), L449P (99.6%)K172R (97%)
3CRF02_AGATVD0L449P (98.2%)E109D (5%)
VFL449P (99.5%)E109D (47%)
4BATVD0E109D (18%)
VFNANAE109D (100%)
5BDRV q24hD0K436R (97.2%), L449P (98.4%), S451N (98.4%)T165S (32%), V321I (27%)
VFK436R (98.7%), L449P (98.3%), S451N (98.3%)T165S (98%), V321I (82%)
6ADRV q24hD0Y441S (29%), L449P (98.1%)
VFY441S (99%), L449P (99.6%)
7ADRV q24hD0NAK436R (98.9%), L449P (99.1%)K172R (3%)
VFK436R (98.8%), L449P (99.5%)K172R (100%), R277Del (95%)
8CRF02_AGDRV q24hD0Y132F (40%), L449P (99.2%)K106R (25%), T165S (25%)
VFY132F (99.7%), L449P (99.4%)K106R (99%), T165S (100%)
9CRF02_AGDRV q24hD064M (2%)S373P (50.1%), S373T (47.4%), L449P (98.4%), P453L (50.9%)
VF64M (100%)S373P (99.5%), L449P (99.2%), P453L (97.9%)NA 
11Cpx_06DRV q24hD0S373T (98.9%), K436R (98.7%), L449P (98.6%)− 
VF79A (100%)S373T (99.6%), K436R (98.4%), L449P (58.8%), S451N (58.8%) -
12CRF02_AGDRV q24hD015V (32%)Y132F (98.2%), R286K (32%), L449P (77.7%)
VF15V (100%), 79A (100%)Y132F (99.6%), R286K (99%), L449P (97.1%)
14CRF02_AGDRV q24hD0R286K (20%), L449P (96.6%)
VFR286K (100%), L449P (99.2%)NA
15CRF02_AGDRV q24hD0
VF15V (100%), 64M (100%), 70R (100%)NANA
19CRF02_AGDRV q24hD0K436R (99.4%)NA
VFK436R (99.7%)NA
23CRF02_AGDRV q24hD0
VF15V (98%), 64M (95.9%), 70R (95.8%)NA NA
25CRF45_cpxLPVD0S373T (97.6%), L449P (99.3%), P453L (98.9%)NA
VFS373T (99.4%), L449P (99.1%), P453L (98.7%)
27ADRV q24hD0
VFNANAV321I (100%)
31CRF02_AGATVD0NA-
VFNANAK106R (100%)
PatientSubtypeARTTimeProtease (%)Gag (%)Gp41 (%)
1CRF02_AGDRV q24hD0S373P (97.7%), L449P (99%)
VFL449P (99.6%)
2CRF02_AGDRV q24hD0K436R (98.5%), L449P (98.7%)
VFK436R (99.1%), L449P (99.6%)K172R (97%)
3CRF02_AGATVD0L449P (98.2%)E109D (5%)
VFL449P (99.5%)E109D (47%)
4BATVD0E109D (18%)
VFNANAE109D (100%)
5BDRV q24hD0K436R (97.2%), L449P (98.4%), S451N (98.4%)T165S (32%), V321I (27%)
VFK436R (98.7%), L449P (98.3%), S451N (98.3%)T165S (98%), V321I (82%)
6ADRV q24hD0Y441S (29%), L449P (98.1%)
VFY441S (99%), L449P (99.6%)
7ADRV q24hD0NAK436R (98.9%), L449P (99.1%)K172R (3%)
VFK436R (98.8%), L449P (99.5%)K172R (100%), R277Del (95%)
8CRF02_AGDRV q24hD0Y132F (40%), L449P (99.2%)K106R (25%), T165S (25%)
VFY132F (99.7%), L449P (99.4%)K106R (99%), T165S (100%)
9CRF02_AGDRV q24hD064M (2%)S373P (50.1%), S373T (47.4%), L449P (98.4%), P453L (50.9%)
VF64M (100%)S373P (99.5%), L449P (99.2%), P453L (97.9%)NA 
11Cpx_06DRV q24hD0S373T (98.9%), K436R (98.7%), L449P (98.6%)− 
VF79A (100%)S373T (99.6%), K436R (98.4%), L449P (58.8%), S451N (58.8%) -
12CRF02_AGDRV q24hD015V (32%)Y132F (98.2%), R286K (32%), L449P (77.7%)
VF15V (100%), 79A (100%)Y132F (99.6%), R286K (99%), L449P (97.1%)
14CRF02_AGDRV q24hD0R286K (20%), L449P (96.6%)
VFR286K (100%), L449P (99.2%)NA
15CRF02_AGDRV q24hD0
VF15V (100%), 64M (100%), 70R (100%)NANA
19CRF02_AGDRV q24hD0K436R (99.4%)NA
VFK436R (99.7%)NA
23CRF02_AGDRV q24hD0
VF15V (98%), 64M (95.9%), 70R (95.8%)NA NA
25CRF45_cpxLPVD0S373T (97.6%), L449P (99.3%), P453L (98.9%)NA
VFS373T (99.4%), L449P (99.1%), P453L (98.7%)
27ADRV q24hD0
VFNANAV321I (100%)
31CRF02_AGATVD0NA-
VFNANAK106R (100%)

ATV, atazanavir; D0, day 0 (baseline); DRV, darunavir; LPV, lopinavir; NA, not amplified, VF, virological failure.

Emerging mutations are shown in bold.

Table 2.

Resistance mutations detected at baseline and mutations selected at virological failure (ultradeep sequencing) in protease, gag and gp41 genes

PatientSubtypeARTTimeProtease (%)Gag (%)Gp41 (%)
1CRF02_AGDRV q24hD0S373P (97.7%), L449P (99%)
VFL449P (99.6%)
2CRF02_AGDRV q24hD0K436R (98.5%), L449P (98.7%)
VFK436R (99.1%), L449P (99.6%)K172R (97%)
3CRF02_AGATVD0L449P (98.2%)E109D (5%)
VFL449P (99.5%)E109D (47%)
4BATVD0E109D (18%)
VFNANAE109D (100%)
5BDRV q24hD0K436R (97.2%), L449P (98.4%), S451N (98.4%)T165S (32%), V321I (27%)
VFK436R (98.7%), L449P (98.3%), S451N (98.3%)T165S (98%), V321I (82%)
6ADRV q24hD0Y441S (29%), L449P (98.1%)
VFY441S (99%), L449P (99.6%)
7ADRV q24hD0NAK436R (98.9%), L449P (99.1%)K172R (3%)
VFK436R (98.8%), L449P (99.5%)K172R (100%), R277Del (95%)
8CRF02_AGDRV q24hD0Y132F (40%), L449P (99.2%)K106R (25%), T165S (25%)
VFY132F (99.7%), L449P (99.4%)K106R (99%), T165S (100%)
9CRF02_AGDRV q24hD064M (2%)S373P (50.1%), S373T (47.4%), L449P (98.4%), P453L (50.9%)
VF64M (100%)S373P (99.5%), L449P (99.2%), P453L (97.9%)NA 
11Cpx_06DRV q24hD0S373T (98.9%), K436R (98.7%), L449P (98.6%)− 
VF79A (100%)S373T (99.6%), K436R (98.4%), L449P (58.8%), S451N (58.8%) -
12CRF02_AGDRV q24hD015V (32%)Y132F (98.2%), R286K (32%), L449P (77.7%)
VF15V (100%), 79A (100%)Y132F (99.6%), R286K (99%), L449P (97.1%)
14CRF02_AGDRV q24hD0R286K (20%), L449P (96.6%)
VFR286K (100%), L449P (99.2%)NA
15CRF02_AGDRV q24hD0
VF15V (100%), 64M (100%), 70R (100%)NANA
19CRF02_AGDRV q24hD0K436R (99.4%)NA
VFK436R (99.7%)NA
23CRF02_AGDRV q24hD0
VF15V (98%), 64M (95.9%), 70R (95.8%)NA NA
25CRF45_cpxLPVD0S373T (97.6%), L449P (99.3%), P453L (98.9%)NA
VFS373T (99.4%), L449P (99.1%), P453L (98.7%)
27ADRV q24hD0
VFNANAV321I (100%)
31CRF02_AGATVD0NA-
VFNANAK106R (100%)
PatientSubtypeARTTimeProtease (%)Gag (%)Gp41 (%)
1CRF02_AGDRV q24hD0S373P (97.7%), L449P (99%)
VFL449P (99.6%)
2CRF02_AGDRV q24hD0K436R (98.5%), L449P (98.7%)
VFK436R (99.1%), L449P (99.6%)K172R (97%)
3CRF02_AGATVD0L449P (98.2%)E109D (5%)
VFL449P (99.5%)E109D (47%)
4BATVD0E109D (18%)
VFNANAE109D (100%)
5BDRV q24hD0K436R (97.2%), L449P (98.4%), S451N (98.4%)T165S (32%), V321I (27%)
VFK436R (98.7%), L449P (98.3%), S451N (98.3%)T165S (98%), V321I (82%)
6ADRV q24hD0Y441S (29%), L449P (98.1%)
VFY441S (99%), L449P (99.6%)
7ADRV q24hD0NAK436R (98.9%), L449P (99.1%)K172R (3%)
VFK436R (98.8%), L449P (99.5%)K172R (100%), R277Del (95%)
8CRF02_AGDRV q24hD0Y132F (40%), L449P (99.2%)K106R (25%), T165S (25%)
VFY132F (99.7%), L449P (99.4%)K106R (99%), T165S (100%)
9CRF02_AGDRV q24hD064M (2%)S373P (50.1%), S373T (47.4%), L449P (98.4%), P453L (50.9%)
VF64M (100%)S373P (99.5%), L449P (99.2%), P453L (97.9%)NA 
11Cpx_06DRV q24hD0S373T (98.9%), K436R (98.7%), L449P (98.6%)− 
VF79A (100%)S373T (99.6%), K436R (98.4%), L449P (58.8%), S451N (58.8%) -
12CRF02_AGDRV q24hD015V (32%)Y132F (98.2%), R286K (32%), L449P (77.7%)
VF15V (100%), 79A (100%)Y132F (99.6%), R286K (99%), L449P (97.1%)
14CRF02_AGDRV q24hD0R286K (20%), L449P (96.6%)
VFR286K (100%), L449P (99.2%)NA
15CRF02_AGDRV q24hD0
VF15V (100%), 64M (100%), 70R (100%)NANA
19CRF02_AGDRV q24hD0K436R (99.4%)NA
VFK436R (99.7%)NA
23CRF02_AGDRV q24hD0
VF15V (98%), 64M (95.9%), 70R (95.8%)NA NA
25CRF45_cpxLPVD0S373T (97.6%), L449P (99.3%), P453L (98.9%)NA
VFS373T (99.4%), L449P (99.1%), P453L (98.7%)
27ADRV q24hD0
VFNANAV321I (100%)
31CRF02_AGATVD0NA-
VFNANAK106R (100%)

ATV, atazanavir; D0, day 0 (baseline); DRV, darunavir; LPV, lopinavir; NA, not amplified, VF, virological failure.

Emerging mutations are shown in bold.

Two patients were found to have viruses containing the mutation R286K in Gag, outside the CS. Both patients were infected with a CRF02_AG virus and treated with a combination of two NRTIs and boosted darunavir q24h. At baseline, these two patients’ viruses also harboured CS Gag mutations in >90% of the total viral population: one patient’s virus harboured the mutation L449P and the other patient’s virus harboured the mutations Y132F and L449P. In addition, 13 other patients’ viruses harboured mutations of CS before treatment began (Table 2). Most of these mutations were present in more than 90% of total viral population: Y132F (p17/p24), S373P and S373T (p17/p24), K436R (p7/p1), L449P, S451N and P453L (p1/p6).

In the gp41 region, two patients’ viruses had selected a mutation inside the CT, V321I (one patient with subtype A and one with subtype B viruses; both patients were treated with a combination of two NRTIs and boosted darunavir q24h), and eight had selected mutations outside the CT: K106R (two patients with CRF02_AG viruses; one patient was treated with boosted darunavir q24h and the other with boosted atazanavir), E109D (one patient with CRF02_AG and one with subtype B viruses; both patients were treated with two NRTIs and boosted atazanavir), T165S (one patient with CRF02_AG and one with subtype B viruses; both patients were treated with two NRTIs and boosted darunavir q24h); and K172R (one patient with CRF02_AG and one with subtype A viruses; both patients were treated with two NRTIs and boosted darunavir q24h) (Table 2).

Sanger sequencing of integrase and reverse transcriptase genes was performed at VF for 17/32 patients whose samples were still available. No mutation was observed in the integrase gene. For the reverse transcriptase gene, one mutation (M184I) was detected at failure in patient 12.

Discussion

In 14 patients failing a boosted PI first-line regimen we identified the emergence of protease, Gag and gp41 mutations that have not previously been described or clearly associated with PI resistance. Eleven were treated with boosted darunavir once daily.

These results highlight the emergence of a pattern of protease mutations (15V, 64M, 70R) among CRF02_AG viruses under darunavir pressure. These mutations may affect the virological response of HIV-1 non-B subtype to PIs, and when they are detected a more rigorous follow-up might be considered for these patients. The mutation 15V has not previously been associated with darunavir resistance but has already been associated with resistance to saquinavir. The mutations 64M and 70R have already been identified in CRF02_AG HIV as polymorphisms, but not as emergent mutations at failure.

We did not identify the emergence of the previously described CS Gag mutation at VF. However, one mutation emerged outside CS (R286K). This mutation has already been associated with a significantly higher rate of VF under a PI-based regimen and contributes to PI resistance.18 Before treatment commenced, mutations identified in CS were already present for many patients in >90% of the total viral population.9,10 Even if these mutations were not associated with VF of darunavir, they might have influenced the selection of treatment-associated substitutions.

To the best of our knowledge, data on emerging mutations in the gp41 gene in cases of VF of a PI regimen are reported here for the first time. It would be interesting to study further the mutation V321I located on CT which has been identified in two viruses’ gp41 gene at VF.14 This mutation could be linked to some interactions between the uncleaved Gag polyprotein and the CT of the Env protein, and may have a potential impact on virological response to PIs.14

Our study is limited by its sample size, mainly as a result of the strict criteria applied to patient selection. Moreover, Sanger sequencing of the reverse transcriptase gene could not be performed for every case. Nevertheless, the sequencing was performed on protease, gag and gp41 genes by a very sensitive technology, at baseline and at VF, in treatment-compliant patients.

In conclusion, this study shows selection of some protease, Gag and Gp41 mutations that have not previously been described and could be related to VF. Site-directed mutagenesis and phenotypic studies should be conducted to improve understanding of these mutations and better characterize their impact on susceptibility to PIs.

Acknowledgements

We thank Isabelle Malet for technical assistance.

Funding

This work was supported by the Agence Nationale de Recherche sur le SIDA et les hépatites virales (ANRS). This study was also funded by Janssen-Cilag and BMS. The study sponsors agreed with the study design, but otherwise had no role in data collection and analysis, decision to publish or preparation of the manuscript.

Transparency declarations

None to declare.

References

1

Pettit
SC
,
Moody
MD
,
Wehbie
RS
et al.
The p2 domain of human immunodeficiency virus type 1 Gag regulates sequential proteolytic processing and is required to produce fully infectious virions
.
J Virol
1994
;
68
:
8017
27
.

2

Chen
Z
,
Li
Y
,
Schock
HB
et al.
Three-dimensional structure of a mutant HIV-1 protease displaying cross-resistance to all protease inhibitors in clinical trials
.
J Biol Chem
1995
;
270
:
21433
6
.

3

Gulnik
SV
,
Suvorov
LI
,
Liu
B
et al.
Kinetic characterization and cross-resistance patterns of HIV-1 protease mutants selected under drug pressure
.
Biochemistry
1995
;
34
:
9282
7
.

4

Mahalingam
B
,
Louis
JM
,
Reed
CC
et al.
Structural and kinetic analysis of drug resistant mutants of HIV-1 protease
.
Eur J Biochem
1999
;
263
:
238
45
.

5

Nijhuis
M
,
Deeks
S
,
Boucher
C.
Implications of antiretroviral resistance on viral fitness
.
Curr Opin Infect Dis
2001
;
14
:
23
8
.

6

Borman
AM
,
Paulous
S
,
Clavel
F.
Resistance of human immunodeficiency virus type 1 to protease inhibitors: selection of resistance mutations in the presence and absence of the drug
.
J Gen Virol
1996
;
77
:
419
26
.

7

Nijhuis
M
,
Schuurman
R
,
de Jong
D
et al.
Increased fitness of drug resistant HIV-1 protease as a result of acquisition of compensatory mutations during suboptimal therapy
.
AIDS
1999
;
13
:
2349
59
.

8

Larrouy
L
,
Chazallon
C
,
Landman
R
et al.
Gag mutations can impact virological response to dual-boosted protease inhibitor combinations in antiretroviral-naïve HIV-infected patients
.
Antimicrob Agents Chemother
2010
;
54
:
2910
9
.

9

Nijhuis
M
,
van Maarseveen
NM
,
Lastere
S
et al.
A novel substrate-based HIV-1 protease inhibitor drug resistance mechanism
.
PLoS Med
2007
;
4
:
e36.

10

Gatanaga
H
,
Suzuki
Y
,
Tsang
H
et al.
Amino acid substitutions in Gag protein at non-cleavage sites are indispensable for the development of a high multitude of HIV-1 resistance against protease inhibitors
.
J Biol Chem
2002
;
277
:
5952
61
.

11

Kletenkov
K
,
Hoffmann
D
,
Böni
J
et al.
Role of Gag mutations in PI resistance in the Swiss HIV cohort study: bystanders or contributors?
J Antimicrob Chemother
2017
;
72
:
866
75
.

12

Malet
I
,
Roquebert
B
,
Dalban
C
et al.
Association of Gag cleavage sites to protease mutations and to virological response in HIV-1 treated patients
.
J Infect
2007
;
54
:
367
74
.

13

Maguire
MF
,
Guinea
R
,
Griffin
P
et al.
Changes in human immunodeficiency virus type 1 Gag at positions L449 and P453 are linked to I50V protease mutants in vivo and cause reduction of sensitivity to amprenavir and improved viral fitness in vitro
.
J Virol
2002
;
76
:
7398
406
.

14

Rabi
SA
,
Laird
GM
,
Durand
CM
et al.
Multi-step inhibition explains HIV-1 protease inhibitor pharmacodynamics and resistance
.
J Clin Invest
2013
;
123
:
3848
60
.

15

Bennett
DE
,
Camacho
RJ
,
Otelea
D
et al.
Drug resistance mutations for surveillance of transmitted HIV-1 drug resistance: 2009 update
.
PLos One
2009
;
4
:
e4724
.

16

Conseil national du Sida. Prise en charge médicale des patients vivant avec le VIH. https://cns.sante.fr/wp-content/uploads/2018/05/experts-vih_suivi.pdf.

17

Kearse
M
,
Moir
R
,
Wilson
A
et al.
Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data
.
Bioinformatics
2012
;
28
:
1647
9
.

18

Blanch-Lombarte
O
,
Santos
J
,
Peña
R
et al. Virological failure to Protease inhibitors in Monotherapy is linked to the presence of signature mutations in Gag without changes in HIV-1 replication. In: Abstracts of the 15th European Meeting on HIV & Hepatitis, Rome, Italy,
2017
. Abstract 03, p. 5.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://dbpia.nl.go.kr/journals/pages/open_access/funder_policies/chorus/standard_publication_model)