Surveillance of HIV-1 transmitted integrase strand transfer inhibitor resistance in the UK.

BACKGROUND: HIV treatment guidelines have traditionally recommended that all HIV-positive individuals are tested for evidence of drug resistance prior to starting ART. Testing for resistance to reverse transcriptase inhibitors and PIs is well established in routine care. However, testing for integrase strand transfer inhibitor (InSTI) resistance is less consistent.
OBJECTIVES: To inform treatment guidelines by determining the prevalence of InSTI resistance in a national cohort of recently infected individuals. PATIENTS AND METHODS: Recent (within 4 months) HIV-1 infections were identified using a Recent Infection Testing Algorithm of new HIV-1 diagnoses in the UK. Resistance-associated mutations (RAMs) in integrase, protease and reverse transcriptase were detected by ultradeep sequencing, which allows for the sensitive estimation of the frequency of each resistant variant in a sample.
RESULTS: The analysis included 655 randomly selected individuals (median age = 33 years, 95% male, 83% MSM, 78% white) sampled in the period 2014 to 2016 and determined to have a recent infection. These comprised 320, 138 and 197 samples from 2014, 2015 and 2016, respectively. None of the samples had major InSTI RAMs occurring at high variant frequency (≥20%). A subset (25/640, 3.9%) had major InSTI RAMs occurring only as low-frequency variants (2%-20%). In contrast, 47/588 (8.0%) had major reverse transcriptase inhibitor and PI RAMs at high frequency.
CONCLUSIONS: Between 2014 and 2016, major InSTI RAMs were uncommon in adults with recent HIV-1 infection, only occurring as low-frequency variants of doubtful clinical significance. Continued surveillance of newly diagnosed patients for evidence of transmitted InSTI resistance is recommended to inform clinical practice.

Introduction

In 2007, raltegravir was the first integrase strand transfer inhibitor (InSTI) introduced into clinical practice, initially for treatment-experienced HIV-positive patients requiring rescue therapy and 2 years later for all patients, including treatment-naive patients. Raltegravir was followed by elvitegravir in 2012 as part of a fixed-dose, single-tablet combination (Stribild) that includes the booster cobicistat, emtricitabine and tenofovir disoproxil fumarate. Elvitegravir was subsequently reformulated in combination with tenofovir alafenamide instead of tenofovir disoproxil fumarate (Genvoya). Second-generation InSTIs comprise dolutegravir, which was approved in 2013, and more recently bictegravir coformulated with tenofovir alafenamide and emtricitabine (Biktarvy), which was approved in 2018. Large clinical trials have demonstrated that the second-generation InSTIs are potent suppressors of HIV replication and have good safety and high genetic barriers to the emergence of drug resistance. These features make them the preferred third agent for starting ART in combination with a backbone of two NRTIs., First-line regimens based on NNRTIs or boosted PIs (bPIs) are instead reserved for selected scenarios.

Mirroring European and American guidelines, the British HIV Association (BHIVA) guidelines for the treatment of HIV-1-positive adults recommend that resistance testing by viral partial genome sequencing be undertaken in all newly diagnosed patients prior to starting ART to allow the detection of transmitted drug resistance (TDR). Sequencing should be performed for reverse transcriptase and protease genes, based on studies showing that when the prevalence of TDR in the population exceeds a threshold of 1%–5% it is cost-effective to screen patients to guide treatment selection., In the UK, the prevalence of TDR affecting NRTIs, NNRTIs or PIs peaked at ∼14% in 2002 and has remained stable at 7%–9% since 2006., To date, there is no recommendation for baseline integrase sequencing as little evidence exists of the transmission of InSTI resistance-associated mutations (RAMs) in the UK and worldwide.

Most routine resistance testing is performed using conventional Sanger sequencing technology, which has a variant frequency detection threshold of ∼20% and hence fails to detect variants that are present below this threshold in a patient’s viral population. Next-generation sequencing (NGS) technologies allow the detection of variants present in a sample at a frequency as low as 1%. The clinical significance of low-frequency resistant variants remains under debate. It has been shown that low-frequency variants with mutations affecting the NNRTIs, and to a lesser extent the NRTIs, significantly reduce responses to first-line therapy with two NRTIs plus one NNRTI, while showing no appreciable effect on bPI-based regimens. However, transmission is unlikely to be a source of the large majority of these low-frequency variants in individuals who have recently acquired HIV and consequently would have minimal to no impact on treatment outcome as they would not have been selected under drug pressure.

Several studies have reported no evidence of InSTI major RAMs in treatment-naive or recently infected HIV-1-positive populations using Sanger sequencing. The few studies reporting apparent transmission of InSTI RAMs included mutations that are polymorphic among ART-naive patients (e.g. L74IM, T97A and E157Q)., One notable exception is a study from Taiwan that observed InSTI major RAMs (e.g. Q148HKR and Y143R) in 1.2% of 1307 ART-naive individuals under a specific epidemiological circumstance where there was a large reservoir of InSTI resistance among the treated population.

The UK national reference laboratory receives blood samples from half of all newly diagnosed cases of HIV-1 infection for incidence testing using a Recent Infection Testing Algorithm (RITA). Recently infected individuals are the most relevant population as they are treatment naive and detection of resistance is most likely due to transmission and prior to natural decay. Using this resource, we determined the national prevalence of TDR to InSTIs, NRTIs, NNRTIs and bPIs by performing NGS on samples from newly diagnosed patients identified as infected in the previous 4 months. The findings will inform the clinical utility of baseline resistance testing for integrase in the UK.

Patients and methods

Ethics

PHE has Section 251 approval, which is reviewed annually and provides the legal basis for the collection of HIV patient-level data for public health monitoring purposes. In addition, the HIV surveillance dataset is reviewed annually by the PHE Caldicott Panel to ensure compliance with information governance policies.

Study population

The UK national reference laboratory applies RITA to new HIV-1 diagnoses as a sentinel national surveillance programme. Blood samples from ART-naive individuals collected at HIV or Genitourinary Medicine Clinics in the UK are identified as likely recent infections (within 4 months of sample collection) using a limiting-antigen (LAg) avidity assay with an OD index <1.5. The assay differentiates likely recent from long-standing infection by the strength of HIV-specific antibody–antigen binding. The assay has a misclassification rate of long-standing HIV infections as recent of <1% when RITA is applied and samples close to the OD index cut-off values are more likely to be misclassified. The RITA algorithm also includes matching the sample to individual HIV records of the HIV and AIDS Reporting System (HARS). Individuals with an OD index <1.5 must also have a CD4+ cell count of >200 cells/mm3 and viral load of >1000 copies/mL to be assigned as ‘recent infections’. In the period 2014–16, RITA was applied to 8379 (47.9%) of new diagnoses and 1765 (21.1%) were identified as recent infections (Figure 1). We randomly selected 655 (37.1%) of these plasma specimens with residual volume for NGS analysis, comprising, by year, 320, 138 and 197 samples collected in 2014, 2015 and 2016, respectively. Linked demographic and clinical information was extracted from HARS. Prescription data were used to determine the use of InSTIs in clinical practice in England between 2010 and 2015.

Figure 1.

Flow chart showing study participant selection. Patients were eligible for inclusion if they were newly diagnosed and identified to be recently infected (within 4 months of sampling) using RITA. PR, protease; RT, reverse transcriptase.

NGS

Samples collected in 2014 were extracted using a QIAsymphony Virus/Pathogen DSP Mini Kit (QIAGEN) using 200 μL of plasma, eluted in a final volume of 60 μL and processed using a previously described PCR amplicon-based NGS assay for protease/reverse transcriptase and integrase was amplified in a nested PCR reaction. Briefly, cDNA was generated using 20 μL of RNA, a QIAGEN OneStep RT–PCR Kit (QIAGEN) and primers H10F2 (5′-GCACAYAARGGRATTGGAGGAAATGA-3′) and H10R3 (5′-CCTAGTGGRATGTGTACTTCTGA-3′), both at 1 μM, under the following cycling conditions: 50°C for 40 min; 95°C for 15 min; 35 cycles of 95°C for 30 s, 53°C for 30 s and 72°C for 1 min; and then a final elongation step at 72°C for 4 min. Two microlitres of cDNA was then used in a semi-nested PCR using a Platinum Taq DNA Polymerase Kit (Invitrogen) and 0.4 μM primer H10F2 and 1.6 μM primer H10R2 (5′-CATATGRTGYTTTACTAAACTHTTCCA-3′) under the following cycling conditions: 95°C for 5 min; 35 cycles of 94°C for 30 s, 53°C for 30 s and 72°C for 1 min; and then a final elongation step at 72°C for 2 min. Amplicons for protease/reverse transcriptase and integrase were pooled in equimolar concentration and then sequenced as previously described.

Samples collected in 2015 and 2016 were processed using a sequence-capture WGS assay. Briefly, 350 μL of plasma was extracted using the NucliSENS system on the easyMag platform (bioMérieux) and eluted into a volume of 25 μL, all of which was subjected to DNAse digestion with 0.25 U of TURBO DNase (Thermo Fisher Scientific) in 30 μL reactions incubated for 30 min at 37°C. Digestion products were cleaned up using 2× AMPure XP Beads (Beckman Coulter) following the manufacturer’s instructions, with a final elution volume of 10 μL nuclease-free water. The 10 μL volume of DNAse-digested RNA extracts was used to generate DNA libraries, using the KAPA RNA HyperPrep Kit (Roche). DNA libraries were pooled in a total of 500 ng and hybridized using 120 nt HIV-specific biotinylated oligonucleotide probes and NimbleGen SeqCap target enrichment reagents (Roche) following the manufacturer’s specifications. Following hybridization, the HIV DNA libraries bound to the biotinylated probes were partitioned using magnetic streptavidin-coated beads and subjected to a further 14 cycles of PCR amplification. The concentration of the final pool was quantified using the KAPA SYBR FAST Universal qPCR Kit for Illumina libraries (KAPA Biosystems) on a 7500 Real-Time PCR System (Applied Biosystems) and analysed for fragment size distribution using the High Sensitivity DNA Kit (Agilent) on a 2100 Bioanalyser Instrument, following the specifications of both manufacturers. Sequencing was performed on an Illumina MiSeq instrument using the MiSeq Reagent Kit V2 (300 cycles) (Illumina) according to the manufacturer’s guidelines, with the following minor modifications. The final pool was diluted to 2 nM and denatured with 0.2 N sodium hydroxide for 2 min, incubated for 4 min at 95°C and diluted in kit reagent HT1 to produce 1 mL of a 20 pM solution. This was further diluted to make 700 μL of a 9 pM solution, of which 10% was substituted with 12.5 pM PhiX (Illumina). A total of 600 μL of this final solution was loaded onto the MiSeq cartridge.

Bioinformatic analysis

MiSeq paired-end FASTQs were trimmed for quality with Trimmomatic (v0.39, with LEADING and TRAILING set to 30 and MINLEN to 50) and human reads removed through BWA (v0.7.17) mapping to the human genome (GRCh37), retaining unmapped pairs. Contigs generated from dehumanized reads by de novo assembly using SPAdes (v3.13.1) were split into fragments of approximately 500 nt (depending on length) and BLASTed against 2427 reference genome sequences annotated and aligned by LANL (http://www.hiv.lanl.gov/). An in-house Python script was used to build a draft sequence from the contigs and their locations within the genome alignment, filling gaps from the reference sequence(s) informing the flanking termini (in a fashion somewhat analogous to LASTZ). Two rounds of BWA mapping and consensus derivation (using an in-house C++ script—QuasiBAM) were performed to obtain the final nucleotide frequency table, from which an in-house Perl script derived consensus genome sequences at 20% and 2% nucleotide frequencies at a minimum read depth of 30, the latter frequency previously established as the minimum threshold of the assay.

The positive percentage agreement (PPA) between amplicon and sequence-capture methods was 99.8% (95% CI = 99.7%–99.9%) and 99.0% (95% CI = 98.8%–99.2%) at 20% and 2% variant frequency thresholds, respectively, at the nucleotide level (Tables S1 and S2 and Figure S1, available as Supplementary data at JAC Online).

Drug resistance analysis

The generated consensus sequences were analysed for surveillance drug resistance mutations (SDRMs) using the Calibrated Population Resistance (CPR) tool that uses the WHO 2009 list of SDRMs for PIs and reverse transcriptase inhibitors and the recently proposed 2019 list of SDRMs for InSTIs: T66AIK, E92GQ, G118R, F121Y, E138AKT, G140ACS, Y143CHRS, S147G, Q148HRK, N155H, S230R and R263K. The integrase sequences were also analysed for the presence of the following InSTI accessory mutations: H51Y, Q95K, T97A, A128T, V151A, S153FY, E157Q and G163KR. These mutations have minimal, if any, effect on InSTI susceptibility when present alone, but may contribute to reduced susceptibility in combination with InSTI major RAMs. For determination of mutational load, viral load and CD4+ cell count data were only used if performed within 30 days of the sample used for RITA and sequencing. Additional mutational load data for PIs, NRTIs and NNRTIs were obtained from previously reported reverse transcriptase and protease sequencing of recently infected individuals from 2011 to 2013.

Statistical analyses

Descriptive statistics (medians and IQRs) are provided for continuous variables, whereas frequency distributions are provided for categorical variables. Mutational load datasets were compared using the Mann–Whitney U-test with the significance level set at P < 0.05.

Sequence data

Consensus HIV-1 pol sequences from this study have been submitted to GenBank and may be accessed by the following accession numbers: MT570368–MT571329.

Results

Study population characteristics

The proportion of adults using InSTIs increased significantly from a low of 2.4% in London in 2010 to between 17.4% and 25.7% in the Midlands & East of England and the North of England, respectively, in 2015 (Figure 2).

Figure 2.

Proportion of people on ART including InSTIs in England between 2010 and 2015. The data were estimated from ART prescribing data from NHS England and are stratified by the PHE regions of London, Midlands & East of England (M&E), North of England and South of England. Data for London in 2015 do not include figures for the month of December. Hospitals in the former North West Strategic Health Authority (SHA) preferred to use raltegravir over bPIs as the third agent for a significant portion of the survey period. Figures may include a source of overestimation as raltegravir was used as first-line post-exposure prophylaxis following sexual exposure (PEPSE). NA, data not available.

Sequencing was performed on plasma samples from 655 recently infected individuals collected between 2014 and 2016. The characteristics of the study population are summarized in Table 1. The majority were male (94.5%), of white ethnic background (77.9%) and with a risk factor for HIV infection being MSM (82.9%). Most of the recent infections were from the London region (57.1%) and the median age of the study population was 33 years (IQR = 26.5–41). Median viral load was 5.18 log10 copies/mL (IQR = 4.61–6.03) and median CD4+ cell count was 545 cells/mm3 (IQR = 408–723), in keeping with recency of infection. Subtyping using the pol gene showed most individuals were infected with subtype B (67.5%).

Table 1.

Characteristics of the study population

Gender, n (%)
male	619 (94.5)
female	36 (5.5)
Risk exposure, n (%)
MSM	543 (82.9)
heterosexual male	48 (7.3)
heterosexual female	33 (5.0)
IVDU	3 (0.5)
other/unknown	28 (4.3)
Ethnicity, n (%)
white	510 (77.9)
black (African/Caribbean/other)	41 (6.3)
other/unknown	104 (15.9)
Region, n (%)
London	374 (57.1)
North of England	113 (17.3)
Midlands & East of England	83 (12.7)
South of England	78 (11.9)
Northern Ireland	6 (0.9)
Wales	1 (0.2)
Subtype, n (%)
A	24 (3.7)
B	442 (67.5)
C	38 (5.8)
D	3 (0.5)
F	27 (4.1)
G	4 (0.6)
CRF01_AE	28 (4.3)
CFR02_AG	54 (8.2)
CRF06_cpx	18 (2.7)
CRF07_BC	1 (0.2)
CRF11_cpx	3 (0.5)
CRF12_BF	3 (0.5)
CRF24_BG	1 (0.2)
complex recombinants	9 (1.4)
Age (years), median (IQR)	33 (26.5–41)
CD4+ cell count (cells/mm3), median (IQR)	545 (408–723)a
Viral load (log10 copies/mL), median (IQR)	5.18 (4.61–6.03)b

n = 602, number of samples with CD4+ cell counts done within 30 days of date of collection of the sample used for RITA and sequencing.

n = 366, number of samples with viral load done within 30 days of date of collection of the sample used for RITA and sequencing.

Prevalence of integrase RAMs

Of the 655 samples sequenced, 640 (97.7%) generated good-quality integrase gene sequence data (complete gene coverage at minimum read depth of 100): 316, 132 and 192 in 2014, 2015 and 2016, respectively. No InSTI major RAMs were detected in the 640 sequences as high-frequency variants (≥20%). A total of 25 (3.9%) sequences contained major InSTI RAMs as low-frequency variants occurring at a frequency between 2% and 20% (Figure 3a). By year, 18 (5.7%), 3 (2.3%) and 4 (2.1%) sequences contained InSTI major RAMs as low-frequency variants in 2014, 2015 and 2016, respectively. In contrast, 39 (6.1%) sequences contained InSTI accessory mutations as high-frequency variants and 18 (2.8%) as low-frequency variants (Figure 3a). By year, the numbers of sequences containing InSTI accessory mutations were 19 (6.0%), 11 (8.3%) and 9 (4.7%) as high-frequency variants and 8 (2.5%), 7 (5.3%) and 3 (1.6%) as low-frequency variants in 2014, 2015 and 2016, respectively.

Figure 3.

Prevalence of InSTI RAMs in recently infected individuals in the UK between 2014 and 2016. Stacked column graphs show the prevalence of InSTI major and accessory RAMs (a) and types of InSTI RAMs (b) at high (>20%) and low (2%–20%) variant frequency. This was determined from 640 integrase sequences from recently infected individuals consisting of 316, 132 and 192 sequences in 2014, 2015 and 2016, respectively. Major and accessory InSTI RAMs were determined using the Stanford HIV Drug Resistance Database 2019 InSTI SDRM list and InSTI accessory mutations as of October 2019.

The InSTI major RAMs and accessory mutations detected are shown in Figure 3(b). The most common InSTI major RAMs detected as low-frequency variants were E138K (8/25; 32.0%) and Q148K (7/25; 28.0%). On the other hand, the most common InSTI accessory mutation detected as a high-frequency variant was E157Q (29/39; 74.4%), whereas V151A was the most common InSTI accessory mutation detected as a low-frequency variant (5/18; 27.8%). Most of the InSTI accessory mutations present as high-frequency variants were associated with non-B subtypes (22/39; 56.4%) with the majority associated with the circulating recombinant forms CRF02_AG (n = 9) and CRF06_cpx (n = 8). The presence of low-frequency RAMs was confirmed by read-based RAM analysis using the variant frequency files produced by the QuasiBam software (Table S3).

Prevalence of reverse transcriptase and protease RAMs

Of the 655 samples sequenced, 619, 593 and 588 generated good-quality sequence data for protease, reverse transcriptase and both gene regions, respectively. By year, 302, 129 and 188 protease and 295, 116 and 182 reverse transcriptase sequences were generated for 2014, 2015 and 2016, respectively. TDR prevalence for all drug classes for the period 2014–16 was 8.0% (47/588) for high-frequency variants and 10.9% (64/588) for low-frequency variants. TDR mutations were detected as high-frequency variants in 15 (2.4%), 20 (3.4%) and 17 (2.9%) sequences against PIs, NRTIs and NNRTIs, respectively (Figure 4a). The overall prevalence of TDR low-frequency variants against PIs, NRTIs and NNRTIs was 27 (4.4%), 29 (4.9%) and 12 (2.0%), respectively (Figure 4a and Table S3). The most common TDR mutations detected as high-frequency variants were L90M (7/15; 46.7%), T215rev (11/20; 55.0%) and K103N (12/17; 70.6%) against PIs, NRTIs and NNRTIs, respectively (Figure 4b). In contrast, the most common TDR mutations detected as low-frequency variants were M46IL (13/27; 48.1%), D67GNE (13/29; 44.8%) and G190E (6/12; 50.0%) against PIs, NRTIs and NNRTIs, respectively (Figure 4b).

Figure 4.

Prevalence of PI and reverse transcriptase inhibitor TDR mutations in recently infected individuals in the UK between 2014 and 2016. Stacked column graphs show the prevalence of PI and reverse transcriptase inhibitor TDR mutations (a) and types of PI, NRTI and NNRTI RAMs (b) at high (>20%) and low (2%–20%) variant frequency. This was determined from 619 protease and 593 reverse transcriptase sequences from recently infected individuals consisting of 302, 129 and 188 protease sequences and 295, 116 and 182 reverse transcriptase sequences in 2014, 2015 and 2016, respectively. PI, NRTI and NNRTI TDR mutations were determined using the WHO 2019 SDRM list.

Of the 655 sequences, 581 (88.7%) generated sufficient and good-quality data in all three polymerase gene regions. Two (0.3%) of the samples had a low-frequency InSTI major RAM and a high-frequency RAM in protease or reverse transcriptase: Q148K + G190A (NNRTI) and E92G + T215S (NRTI). Both were subtype B. Three samples (0.5%) had a low-frequency InSTI major RAM (E138K) and low-frequency RAM in protease and/or reverse transcriptase: D67N (NRTI) + M46I (PI), M46I (PI) and D67G (NRTI) (they were subtype F, B and CRF01_AE, respectively).

Mutational load of low-frequency RAMs

We determined the mutational load of low-frequency RAMs, as previously described. The median mutational load of low-frequency InSTI RAMs was 3833 copies/mL (IQR = 895–14 733) (Figure 5). The median mutational load of low-frequency InSTI RAMs was similar to that of low-frequency PI and NRTI RAMs at 5914 copies/mL (IQR = 1470–31 552) and 2706 copies/mL (IQR = 1387–11 726), respectively. In contrast, the median mutational load for low-frequency NNRTI RAMs was slightly higher and had a broad range at 10 188 copies/mL (IQR = 2170–95 313); however, the difference was not statistically significant (P > 0.05, Mann–Whitney U-test).

Figure 5.

Mutational load of low-frequency RAMs. Box-and-whisker plots of mutational load of low-frequency InSTI (n = 14), PI (n = 44), NRTI (n = 28) and NNRTI (n = 14) RAMs. Mutational load is defined as mutation frequency × viral load (copies/mL).

Discussion

Surveillance of transmitted InSTI resistance among 655 individuals with recent HIV infection in the UK who were sampled between 2014 and 2016 showed no evidence of major InSTI RAMs when considering mutants present at high frequency in the individuals’ samples. In contrast, at approximately 8%, the prevalence of TDR to PIs and reverse transcriptase inhibitors remains steady compared with last reported figures in 2014. This prevalence is still higher than the recommended threshold of 1%–5% where baseline resistance testing is considered of benefit at the population level.

On the other hand, accessory InSTI RAMs were detected as high-frequency variants in 6.1% of the study population. The most common accessory InSTI RAMs were T97A and E157Q; both are polymorphic and observed at a high prevalence (up to 7%) in InSTI-naive individuals infected with non-B subtypes, e.g. CRF02_AG. They are also selected in patients experiencing treatment failure with first-generation InSTIs, raltegravir and elvitegravir; however, they have little effect on InSTI susceptibility when present alone.

Ultradeep sequencing allowed the detection of mutations below the Sanger sequencing variant frequency threshold of ∼20%. Major InSTI RAMs were detected as low-frequency variants in 3.9% of the study population, at a variant frequency between 2% and 20%. It has been argued that the mutational load of low-frequency RAMs, especially for NNRTIs, could play a role in treatment failure. The mutational load of low-frequency InSTI RAMs was comparable to that of NRTI and PI RAMs, but was lower and had a very narrow range compared with that of NNRTI RAMs. Compared with the data for low-frequency NNRTI RAMs, there is less compelling evidence that low-frequency PI and NRTI RAMs contribute to treatment failure. It is likely that the mutational load of low-frequency RAMs is associated with the impact on virus replication fitness, therefore InSTI, NRTI and PI RAMs that have a high impact on virus replication fitness are unlikely to accumulate to high absolute levels compared with NNRTI RAMs. In addition, we recently showed evidence that the majority of low-frequency RAMs to PIs and reverse transcriptase inhibitors in recently infected individuals are not a result of a transmission event and thus would not have been selected under drug pressure. Furthermore, a recent study showed no association between the presence of low-frequency InSTI RAMs prior to initiation of treatment and treatment outcomes. Taken together, these data suggest that the low-frequency InSTI RAMs in recently infected individuals are less likely to affect treatment outcome, especially as current second-generation InSTIs, dolutegravir and bictegravir, are highly effective and have very high barriers to resistance. However, low-frequency InSTI RAMs may still have an impact in individuals with poor adherence or those who harbour resistance to other components of their ART regimen.

The proportion of individuals on an ART regimen that included an InSTI was approximately 20% during the period covered by the study, as estimated using prescription data from NHS England. The use of InSTIs as part of first-line regimens is anticipated to continue to rise in the UK, reflecting national and international treatment guidelines. The use of InSTIs in the UK increased from less than 10% in 2014 to over 20% in 2015; thus, the effect of this and further projected increases in InSTI use may not be captured in this surveillance study. Nonetheless, the virological suppression rate for people on InSTI-based therapy in the UK is very high (>95%) and the likelihood of the emergence of drug resistance for those failing dolutegravir or bictegravir plus two NRTIs in first-line ART is negligible., However, raltegravir and elvitegravir have been used for longer than dolutegravir and bictegravir in ART-naive and ART-experienced patients and these drugs are more likely to result in treatment failure with resistance selection. Thus, these groups may have generated a pool of potential transmitters, which may later contribute to transmitted InSTI resistance. All these factors necessitate continued surveillance of InSTI TDR in the coming years.

Reflecting the focus on recent infection, another limitation of the study is that most of the sampled population was male of white ethnic background from England whose probable route of HIV exposure was sex between men and who were infected with subtype B virus. This is because gay and bisexual men are more likely to have recently acquired infection at HIV diagnosis. Thus, these findings may not be generalizable to the whole of the UK population living with HIV, particularly women and those infected with non-B subtypes. In addition, the frequency of InSTI use may be different in Scotland, Wales and Northern Ireland, data that were not captured in this study.

Two different sequencing methods were used in this study, an amplicon-based approach using gene-specific nested PCR followed by DNA library prep and a sequence-capture approach that is dependent on RNA library prep followed by enrichment using HIV-1-specific probe baits. Overall, the consensus sequence generated by both methods was highly concordant (>99%) using the 20% and 2% variant frequency at nucleotide level (see Table S1). Discordances were at mixed-base positions where one method detected only one of the mixed bases, with a tendency for low-frequency variants detected by the amplicon method being infrequently detected by the sequence-capture method, whereas the opposite was true (see Table S2 and Figure S1). This could either be due to an overcall of low-frequency variants by the amplicon method or a decreased sensitivity for detection of low-frequency variants by the sequence-capture method. This requires further investigation using standardized reference or control material with well-characterized low-frequency variants at specific frequencies. The validation of NGS methods to accurately reflect in vivo low-frequency variants is essential to determine their effect on clinical outcomes.

In conclusion, this study shows no evidence of transmitted InSTI resistance in the recently infected population in the UK. However, performing baseline integrase resistance testing is still important, especially for national reference laboratories in order to provide surveillance data and possibly in selected patients in routine clinical practice. One consideration is that the use of InSTIs as part of first-line combination ART (cART) is anticipated to continue to increase worldwide following WHO recommendations. The large-scale use of dolutegravir in resource-limited settings is likely to take place with limited viral load monitoring and thus could result in significant increases in transmitted InSTI resistance from a more global perspective. In parallel, prospective cohort studies to assess treatment outcomes in recently infected individuals harbouring low-frequency RAMs would best inform their clinical significance and diagnostic utility. The use of WGS adopted from 2015 onwards will also be useful in analysing other regions of the HIV-1 genome that have been postulated to be involved in development of resistance to second-generation InSTIs, e.g. envelope and 3′ polypurine tract (PPT).,

Click here for additional data file.