INTRODUCTION
Exposure to first-generation HIV-1 integrase strand transfer inhibitors (INSTIs), such as raltegravir (RAL) and elvitegravir (EVG), is frequently associated with the selection of resistance-associated mutations (RAMs) at virological failure. Recently developed second-generation INSTIs, such as dolutegravir (DTG; approved by FDA and European Medicines Agency [EMA] in 2013 and 2014, respectively), bictegravir (BIC; approved by FDA and EMA in 2018 as part of a single tablet regimen, including tenofovir alafenamide [TAF] and emtricitabine [FTC]), and cabotegravir (CAB; currently under phase III clinical investigation) have demonstrated a superior genetic barrier to resistance and a variable activity against viruses harboring INSTI resistance mutations selected by RAL and EVG. However, predicting the activity of second-generation INSTIs on HIV-1 variants with different combination of INSTI RAMs is not straightforward. This study aimed at clarifying cross-resistance to DTG, BIC, and CAB in a panel of INSTI-resistant strains isolated from patients previously exposed to RAL, EVG, and/or DTG.
Plasma samples from 19 patients were collected during routine drug resistance testing under virological failure of an INSTI regimen, based on the detection of at least one major INSTI RAM (T66A, E92Q, E138K/T, G140A/C/S, Y143C/R, S147G, Q148H/R, N155H, and R263K). Virological failure was defined at the individual clinics as confirmed viral load of >50 HIV-1 RNA copies/ml or one viral load of >1,000 HIV-1 RNA copies/ml. Access to residual samples for research purposes was regulated by informed consent approved by the Southeast Tuscany Ethics Committee. Plasma RNA was extracted by the DSP virus kit using the EZ1 workstation (Qiagen, Hilden, Germany). Population sequencing of the integrase region was performed by homebrew technology using primers previously described (
1). For
in vitro susceptibility testing, the integrase coding region was amplified through a 2-step PCR protocol, and the resulting amplicon was used for the generation of chimeric viruses by homologous recombination with an integrase-deleted pNL4-3 vector in 293T cells (
2). Integrase sequences of recombinant viruses were confirmed by Sanger sequencing.
Susceptibility to DTG, BIC, and CAB was assessed through a TZM-bl cell line-based phenotypic assay, shown to correlate well with the
de facto reference PhenoSense assay by Monogram Biosciences and expressed as fold change (FC) with respect to the reference wild type NL4-3 virus (
2). Differences in FC values among groups were tested by Friedman’s ANOVA, and multiple comparisons were performed by the Dunn’s test when appropriate.
The patients had been exposed to RAL only (
n = 9), EVG only (
n = 4), RAL and DTG (
n = 5), or RAL and EVG (
n = 1). By querying the COMET HIV-1 subtyping tool (
3), sequences were assigned to subtype B (
n = 13), subtype F1 (
n = 2), subtype G (
n = 2), CRF03_AB (
n = 1), and CRF06_cpx (
n = 1). Sequences were grouped according to the presence of major INSTI RAMs at codon 143, 148, or 155, alone or in combination with other INSTI RAMs. Four samples with INSTI RAMs E92Q (
n = 2), T66A plus S147SG, and R263K were grouped together as “others” (
Table 1). Median (interquartile range [IQR]) FC values for DTG, BIC, and CAB were 3.5 (1.2 to 7.3), 2.4 (1.4 to 5.4), and 2.3 (1.3 to 21.7), respectively, without statistically significant differences between groups. Median (IQR) DTG, BIC, and CAB FC values from patients not exposed to DTG were 2.9 (1.0 to 4.4), 2.1 (1.3 to 3.3), and 2.3 (1.2 to 4.7), respectively, while exposure to RAL and then DTG resulted in a FC of >100 in three cases (Q148H/K/R plus 2 INSTI RAMs) and FC of <3.5 in two cases (samples with E92Q and with R263K) for all three drugs (
Fig. 1A). Median (IQR) DTG, BIC, and CAB FC values were consistently, but not significantly, higher following exposure to RAL only versus EVG only (3.5 [1.6 to 5.9], 2.4 [1.6 to 3.9], and 2.3 [1.4 to 12.3] versus 1.1 [0.8 to 3.6], 1.8 [0.3 to 3.3], and 1.6 [0.8 to 8.0], respectively) (
Fig. 1B). According to the three major INSTI resistance pathways, median (IQR) DTG, BIC, and CAB FC values were 2.3 (1.0 to 3.5), 2.5 (2.4 to 2.5), and 1.7 (1.2 to 2.2), respectively, with Y143R/C alone (
n = 2); 3.6 (2.6 to 4.4), 1.8 (1.7 to 2.7), and 2.3 (1.5 to 2.6) with N155H alone (
n = 6); 0.5, 0.8, and 1.2 with Q148R alone (
n = 1); 7.3 (1.0 to 8.0), 5.4 (3.2 to 8.0), and 21.7 (9.9 to 66.3) with Q148R plus one additional INSTI RAM (
n = 3); and >100 for all the three drugs with Q148H/K/R plus two additional INSTI RAMs (
n = 3) (
Fig. 1C). In the N155H group, the BIC and CAB FC values were significantly lower than the DTG FC values (
P = 0.005). Among samples with other INSTI RAMs (two with E92Q, one with R263K, and one with T66A plus S147S/G), FC values of ≤3.5 were measured for all the three drugs. In particular, one sample from a failing DTG-based regimen harboring R263K had DTG and CAB FCs of 3.5, while BIC FC was 1.1.
These data are in line with and yet add further information to those reported with specific laboratory clones and with a few clinical samples (
4–7). In the largest data set analyzed comprising 47 patient-derived HIV-1 isolates, BIC appeared to retain more activity against INSTI-resistant mutants than DTG (
4). The median FC value for BIC was also lower than that for DTG in our case file; however, the difference was statistically significant only for the N155H pathway and yet not for the data overall. The Q148H/K/R plus G140S pathway was confirmed to be the most challenging for all second-generation INSTIs, with the addition of a third RAM (particularly at codon 138) substantially increasing the FC value, as shown also in previous studies (
4–7). Contrary to the results generated with molecular clones (
7), we did not detect better activity of DTG than with CAB, likely due to the overrepresentation of N155H single mutants in our data set.
In summary, DTG, BIC, and CAB retain comparable activity against major INSTI RAMs; however, adding one and, particularly, two additional INSTI RAMs to the Q148H/K/R pathway results in a substantial loss of susceptibility for all three drugs. The full spectrum of additional mutations playing such a role in different HIV-1 subtypes remains to be elucidated. More extensive in vivo data are required to fully define the correlation between HIV-1 genotype, FC values, and response to treatment with second-generation INSTIs.
Data availability.
Integrase sequences of study samples were submitted to GenBank under accession numbers
MN389749 to
MN389767.