HIV-1 strains belonging to large phylogenetic clusters show accelerated escape from integrase inhibitors in cell culture compared with viral isolates from singleton/small clusters.

Objectives: Viral phylogenetics revealed two patterns of HIV-1 spread among MSM in Quebec. While most HIV-1 strains ( n  =   2011) were associated with singleton/small clusters (cluster size 1-4), 30 viral lineages formed large networks (cluster size 20-140), contributing to 42% of diagnoses between 2011 and 2015. Herein, tissue culture selections ascertained if large cluster lineages possessed higher replicative fitness than singleton/small cluster isolates, allowing for viral escape from integrase inhibitors.
Methods: Primary HIV-1 isolates from large 20+ cluster ( n  =   11) or singleton/small cluster ( n  =   6) networks were passaged in vitro in escalating concentrations of dolutegravir, elvitegravir and lamivudine for 24-36 weeks. Sanger and deep sequencing assessed genotypic changes under selective drug pressure.
Results: Large cluster HIV-1 isolates selected for resistance to dolutegravir, elvitegravir and lamivudine faster than HIV-1 strains forming small clusters. With dolutegravir, large cluster HIV-1 variants acquired solitary R263K ( n  =   7), S153Y ( n  =   1) or H51Y ( n  =   1) mutations as the dominant quasi-species within 8-12 weeks as compared with small cluster lineages where R263K ( n  =   1/6), S153Y (1/6) or WT species (4/6) were observed after 24 weeks. Interestingly, dolutegravir-associated mutations compromised viral replicative fitness, precluding escalations in concentrations beyond 5-10 nM. With elvitegravir, large cluster variants more rapidly acquired first mutations (T66I, A92G, N155H or S147G) by week 8 followed by sequential accumulation of multiple mutations leading to viral escape (>10 μM) by week 24. Conclusions: Further studies are needed to understand virological features of large cluster viruses that may favour their transmissibility, replicative competence and potential to escape selective antiretroviral drug pressure.

Introduction

ART has transformed HIV/AIDS from a deadly pandemic to a treatable and potentially preventable disease., Expanded access to ART has resulted in 25%–50% declines in global epidemics in heterosexual populations, raising optimism that treatment-as-prevention may arrest the pandemic by 2030 (www.unaids.org). The World Health Organization has advanced global ‘90–90–90’ targets for 2020 where 90% living with HIV-1 are informed of their status, 90% of infected individuals receive antiviral treatment and 90% of those treated achieve viral suppression (www.unaids.org). HIV-1 guidelines have been revised to meet these objectives, recommending annual HIV testing, immediate treatment initiation and pre-exposure prophylaxis for high risk populations.,,,

As more asymptomatic patients are started on ART, treatment goals will require more robust and durable first-line regimens to prevent the development and transmission of drug resistance., The mainstay of treatment is two NRTIs with a third core agent from one of the NNRTI, boosted PI or integrase strand transfer inhibitor (INSTI) drug classes. In developing countries, rates of virological failure (40%) and transmitted drug resistance (5%–16%) are rising, related to the use of thymidine analogues, particularly stavudine and zidovudine, and NNRTIs.,

Integrase inhibitors, including raltegravir, elvitegravir and dolutegravir, are the preferred choice for first-line therapy, based on their improved tolerability, dosing (once daily) and higher barrier to resistance than other antiretroviral drugs (http://aidsinfo.nih.gov/guidelines)., Resistance to elvitegravir and raltegravir occurs via several mutational pathways, including: (i) N155H or G140A/G148RHQ pathways conferring raltegravir and elvitegravir cross-resistance; (ii) T66I or E92Q/G elvitegravir-specific pathways; or (iii) the Y143R/H/C raltegravir-specific resistance pathway. Dolutegravir retains antiviral activity against HIV-1 variants resistant to raltegravir and elvitegravir. Dolutegravir is far less prone to the development of resistance in first-line therapy with isolated reported cases of the emergence of R263K, N155H or G118R. These substitutions have been shown to arise duringin vitro cell culture selection studies., Transmitted resistance to integrase inhibitors remains anecdotal.,

Disturbingly, HIV-1 epidemics among MSM have not declined in resource-rich nations having access to the best antiretroviral treatments. Viral phylogenetic linkage has revealed that epidemics among MSM have been sustained by high rates of transmission clustering related in part to episodic risk effects during highly infectious primary HIV infection (PHI). A shift towards large cluster lineages has contributed to the growth of HIV epidemics among MSM in Quebec and the Netherlands., Whereas half of the HIV-1 epidemic in MSM in Quebec can be ascribed to 2011 viral lineages leading to singleton/small cluster transmissions (cluster size 1–4), 30 viral variants strains have led to large cluster outbreaks (cluster size 20–140), rising from 13% of new diagnoses in 2004 to 42% of new infections in 2015., The rise in large 20+ clusters has offset steady declines in other cluster groups over the last decade., Transmission clustering has been implicated in the spread of resistance to thymidine analogues and NNRTIs.,,

HIV-1 transmission constraints lead to a single monophyletic HIV-1 variant, termed the transmitted/founder virus, establishing most new infections. There is continued debate as to whether transmission bottlenecks are stochastic, randomly restricting all viruses, or selective, favouring specific transmission/founder viruses Observational data from the Montreal PHI cohort showed that cluster size and the skewed role of 30 founder viruses (1.7% of viral lineages) in 1200 forward transmissions were not directly related to patient epidemiological factors, including numbers of reported sexual partnerships and viral load. We postulated that HIV-1 strains associated with large clusters have unique phenotypic features in reverse transcriptase (RT) and integrase replicative processes that may have contributed to heightened infectiousness and cluster burst size.

To test this hypothesis, tissue culture selections with dolutegravir, elvitegravir and lamivudine were used to compare the barriers to resistance for viruses derived from 11 patients belonging to 10 large clusters (cluster size 20–140) and 6 persons associated with singleton/small clusters (cluster size 1–4). Sanger (population) and next-generation (deep) sequencing was performed to monitor genotypic changes under selective drug pressure over 36 weeks. HIV-1 isolates from large cluster lineages demonstrated a lower genetic barrier to resistance to dolutegravir, elvitegravir and lamivudine as compared with viruses from singleton/small cluster networks. However, the rapid acquisition of R263K or S153Y mutations with dolutegravir compromised viral replicative competence and hindered viral breakthrough. Taken together, our findings show a selection bias for large cluster viral variants showing higher replicative fitness under selective drug pressure.

Methods

Viral phylogenetic reconstruction of the HIV-1 epidemic in MSM

Viral phylogenetics was used to reconstruct patterns of HIV-1 spread among newly diagnosed treatment-naive MSM (n = 4039, 2002–15).,, HIV-1 genotyping was performed as previously described, generating partialpol sequences that span the viral protease and RT regions.,, Phylogenetic trees were built using MEGA7 integrated software (www.megasoftware.net).,,, Clustering of viral strains was defined by high bootstrap support (>95%) and short genetic distances (<1.5%). HIV-1 strains from large clusters were resequenced across the entire integrase region as previously described to compare clustering patterns across the protease, RT and integrase regions.

Isolation of viruses from MSM within large and small cluster groups

The Fonds de recherche Santé (FRQS) Réseau SIDA supports a cohort of newly infected persons with clinical indication of primary infection. In this study, HIV-1 strains were isolated from 17 MSM subjects, 11 of whom belonged to 10 large clusters (cluster size 20–140) and 6 from singleton/small cluster transmissions (cluster size 1–4). HIV-1 isolates were amplified as previously described through co-culture of patient CD8-cell-depleted peripheral blood mononuclear cells with primary human cord blood mononuclear cells (CBMCs)., HIV-1 strains, integrase natural polymorphisms, clinical features and GenBank accession numbers are summarized in Table1.

Table 1

Baseline natural polymorphisms in the integrase of subtype B HIV-1 isolates used for in vitro selections with elvitegravir, dolutegravir and lamivudine

Lab IDa	Cluster ID (cluster size)b	Visit 1 date	Months infected (partners)c	CD4 cells/μLc	Integrase polymorphisms	GenBank number
Large
14947	C185.1 (n = 45)	01/2014	2 (1–4)	60	E11D, R20K, V31I, S39N, M50I, I72V, S119T, T124N, G193E, V201I, T218ST, D286N	KT988126
14637	C185.09 (n = 45)	05/2013	3 (1–4)	550	E11D, R20K, V31I, S39N, M50I, S119T, T124N, G193E, V201I, D286N	KT988125
14997	C027.28 (n = 40)	02/2014	3 (1–4)	70	E11D, V31I, L101I, S119G, T122I T124A, T125A, F181L, I203M, K211Q, T218I, D232E, D256E, D286N	KT988128
14969	C053.16 (n = 24)	01/2014	2 (1–4)	220	S17N, S24D, D25E, S39C, V201I, S255Q	KT988127
10387	C099.10 (n = 35)	04/2008	2 (5–9)	580	V31I, L45Q, I72V, K111E, I113V, T125A, M154I, V201I, D207N, L234I	KX714013
10249	C099.08 103N (n = 34)	11/2007	5 (>20)	1010	V31I, L45Q, I72V, L101I, K111E, I113V, D207N, L234I	KX714014
12083	C050.37 G190A (n = 140)	02/2010	6 (>20)	1470	S17N, I72V, I73V, L101I, T124A, K136Q	KX714016
14624	C050.65 G190A (n = 140)	07/2013	9 (1–4)	187	S17N, A23V, L28I, S39C, L101I, T124A	KX714018
12608	C163.05 (n = 90)	05/2011	4 (1–4)	300	E11D, A21T, S39C, L45R, I72V, L101I, G163E, K244R	KX714017
10679	C075.19 (n = 21)	09/2008	3 (1–4)	1250	E11D, A21T, L101I, I113V, S119P, T125A, G193E, D278N	KX714015
5331	C089.01 (n = 39)	02/1997	5 (0)	191	E11D, E13D, I72V, L101I, I113V, S119P, T122I, D256E	KX714020
Small
5326	C196.02 (n = 4)	11/1998	—	—	E11D, A21T, L101I, I113V, S119P, T125A, G193E, D278N	KX714021
15366	C506.01 (n = 2)	12/2014	1 (1–4)	23	S17N, I72V, T124A, T206S, N254H, A265V, S283G	KT988129
14514	S630 (n = 1)	02/2014	5 (1–4)	430	K7R, S17N, M50I, I72V, K111Q, T112V, T124A, T125V, I220L, Y227F, D256E	KT988124
14380	C309 (n = 1)	11/2012	2 (1–4)	730	R20K, V31I, I84V, T112M, T124N, M154I, V201I, D278E	KT988123
5833	S347 (n = 1)	02/2001	—	—	T112V, G193E, V201I, S230N	KT988122
15600	C545.01 (n = 2)	07/2015	2 (1–4)	320	E11D, V31I, A124T, R127K, F181L, V201I, T218I, S230N, N232D, D256E	KX714019

Viruses were isolated from patients participating in the PHI cohort in Montreal.

Cluster size was based on phylogenetic analysis of the entire provincial genotyping programme.

Cohort data, including estimated time post-infection (months) based on biological assays, the reported number of partners 3 months prior to infection and CD4 counts, were determined at the first visit.

In vitro selection of HIV-1 strains with reduced susceptibility to antiretroviral drugs

Dolutegravir was provided by ViiV Healthcare Ltd (Research Triangle Park, NC, USA). Elvitegravir was obtained from Gilead Sciences Inc. (Foster City, CA, USA). Lamivudine was obtained from the NIH AIDS Research and Reference Reagent Program (http://www.aidsreagent.org). HIV-1 strains from 6 and 11 patients from singleton/small and large clusters, respectively, were used to infect CBMCs at a similar moi of 0.1, i.e. 1 × 106 TCID50 per 10 × 106 CBMCs, based on p24 and RT assays. Viruses were serial passaged in CBMCs in the presence of escalating concentrations of dolutegravir, elvitegravir and/or lamivudine over 24–36 weeks, as previously described., RT activity in culture supernatants was quantified weekly to monitor viral replication of individual viral variants in the presence and absence of drug pressure, adjusting incremental increases in drug concentrations if RT activity assays exceeded 100 000 cpm., At each passage, aliquots of cell-free supernatant were stored at −70°C for further genotypic analysis.

HIV-1 genotyping based on Sanger and deep sequencing

Genotypic analyses at selected passages evaluated the development of amino acid substitutions under selective drug pressure, conferring reduced susceptibility to antiretroviral drugs (i.e. drug resistance mutations). Sanger (population) sequencing of viral RNA or complementary DNA extracted from culture supernatants was performed as previously described., Deep sequencing was performed using a CLIA/CAP-validated HIV-1 genotyping assay (DEEPGEN™HIV) with a >1% cut-off for the detection of minority HIV-1 variants based on the intrinsic error rate of the entire system (0.3%). Ultradeep sequencing confirmed the homogeneity of founder variants, as well as the earliest time to emergence of resistance to integrase inhibitors.

Ethics approval

Ethics approvals for phylogenetic surveillance of the HIV-1 epidemic were obtained from the Ministry of Health (Quebec) and L’Institut national de santé publique du Québec. Blood samples from the Montreal primary infection cohort were obtained with informed consent and ethics approval from recruited individuals from participating hospitals or private clinics.

Results

Isolation of HIV-1 strains belonging to representative large and singleton/small clusters

Viral phylogenetics revealed two predominant patterns of HIV-1 spread among MSM in Quebec. Half of new infections were associated with singleton/small cluster networks (cluster size 1–4), while 30 viruses (1.7% of viral lineages) resulted in micro-epidemics, averaging 40 transmissions/cluster (Figure1). Viruses within both cluster groups were isolated to compare viral escape from antiretroviral pressure.

Figure 1

Patterns of spread of the HIV-1 epidemic among the MSM population in Quebec (2002–15). Phylogenetic tree analysis ascertained clustering of new infections genotyped between 2002 and 2015. The frequency of viral lineages (grey) and HIV-1 infections (black) associated with cluster groups, having 1 to 4, 5 to 9, 10 to 19, or 20 to 140 linked transmissions.

Eleven subtype B viruses belonging to 10 large 20+ clusters and 6 singleton/small clusters were derived from PHI cohort participants with estimated dates of infection of 3.7 ± 0.6 months based on Detuned/BED biological assays (Table1). The 11 large cluster strains were collectively associated with 600 transmissions in Quebec. The six singleton/small cluster variants included two reference clinical isolates (5326 and 5833) used in cell culture selections by our laboratory for over a decade. Table 1 shows clinical and genotypic features of the selected viruses, illustrating that cluster size was not related to the number of reported sexual partnerships by cohort participants 3 months prior to their infection, as previously reported for the entire MSM cohort.

Viruses associated with the selected 10 20+ large clusters were resequenced across the integrase-coding regions. Phylogenetic trees showed similar co-clustering patterns across the integrase region as previously observed across the protease/RT domains. It was, however, noteworthy that clusters 45, 118 and 185 (cluster size = 38, 43 and 45, respectively) had a shared common branch in the integrase tree although the protease/RT regions within these three clusters were distinct (Figure2).

Figure 2

Phylogenetic tree analysis shows 10 large clusters sequenced across the protease/RT- and integrase-coding regions. Note that viruses belonging to cluster 185 (n = 44), cluster 45 (n = 38) and cluster 118 (n = 43) share sequence homology across the integrase domain.

In vitro selection of HIV-1 variants resistant to dolutegravir, elvitegravir and lamivudine

Clinical isolates from newly infected persons typically use the CCR5 receptor, show T-cell tropism and do not grow in MT-2 cells. Tissue culture selections were performed using CBMCs to characterize possible differences in the development of drug resistance to lamivudine and integrase inhibitors., CBMCs infected with relevant HIV-1 strains (from large or singleton/small clusters) were serially passaged over the course of 24–36 weeks in the presence of increasing concentrations of dolutegravir, elvitegravir, lamivudine or dolutegravir + lamivudine dual pressure. These drugs are recommended components of choice in first-line regimens.

Cell culture selections using large cluster and singleton/small cluster HIV-1 strains were conducted in parallel under the same conditions. Stepwise drug dose escalations were based on weekly RT assays, performed for each isolate in the presence and absence of drug. As illustrated in Figure3, viruses in both cluster groups replicated equally well in the absence of drug pressure. The progress of viral outgrowth in the presence of stepwise increasing concentrations of inhibitors over the course of 36 weeks, comparing viruses belonging to cluster 185 (14947 and 14637), cluster 27 (14997) and cluster 53 (14969) and four singleton small cluster lineages (HIV-1 isolates 143809, 14515, 15336 and 5833), is depicted in Figure4.

Figure 3

Levels of viral RT activity following initial infection of CBMCs with large cluster (n = 4, filled symbols) and small cluster (n = 4, open symbols) viral lineages. Levels of RT activity remained similar in subsequent serial passage in the absence of drug pressure.

Figure 4

Drug resistance selections with lamivudine, elvitegravir, dolutegravir and dolutegravir + lamivudine. Stepwise increases in drug concentrations under selective drug pressure were based on levels of RT enzymatic activity at each weekly passage relative to a no-drug control. The mean weekly drug concentrations were averaged for the four large cluster viral lineages (black lines) belonging to clusters 27 (14997), 53 (14969) and 185 (14947 and 14637) as compared with four singleton/small cluster variants (isolates 143809, 14515, 15336 and 5833, grey lines) grown under selective pressure with lamivudine, elvitegravir, dolutegravir or dolutegravir + lamivudine over 36 weeks. The arrows represent viral breakthrough based on the first acquisition of new resistance mutations presented in Figures4–6. Regarding lamivudine, one large cluster virus passaged with lamivudine acquired M184V leading to earlier viral escape by week 16; the remaining three large clusters developed M184I. DTG, dolutegravir; EVG, elvitegravir; 3TC, lamivudine.

As illustrated in Figure4(a), drug dose escalations and viral breakthrough, under lamivudine selective pressure, progressed at faster rates in the case of large cluster lineages as compared with singleton/small cluster strains. Drug dose escalations with lamivudine in large cluster variants were associated with the appearance at week 8 of M184V (isolate 14947) or M184I (isolates 14637, 14997 and 14637). The latter isolates shifted from M184I to M184V by weeks 15–23, leading to higher fold resistance (>100-fold) and viral escape. Small cluster isolates developed resistance more slowly, acquiring M184I infections by week 23 (Figure3a). The rapid selection of M184V by isolate 14947 parallels that observed for CXCR4-tropic pNL4.3 reference strain grown in MT-2 cells or CBMCs. The M184I substitution conferred lower fold resistance, decreased processivity and replicative fitness than the M184V mutation., Large cluster variants showed a more rapid shift from M184I to M184V when compared with small cluster lineages where selection of M184I was attenuated.

Similarly, elvitegravir drug selections performed on four large cluster variants progressed more rapidly, leading to the first acquisition of the INSTI resistance mutations T66I (n = 2), N155H (n = 1) or S147G (n = 1) at week 8, followed by the further accumulation of resistance mutations leading to viral escape (>100-fold resistance) at weeks 23 and 36 (Figure3b). In contrast, elvitegravir selections performed on singleton/small cluster viral lineages were delayed with the first appearance of T66I as a singlet mutation occurring at week 23 in only two of four small cluster strains (Figure4b).

Large cluster strains also showed a lower barrier to resistance to dolutegravir with R263K (n = 3) or S153Y (n = 1) mutations appearing by week 8 in large cluster lineages; these mutations were absent in the four small cluster viruses at week 36 (Figure 4c). It was, however, noteworthy that the acquisition of R263K and S153Y had a negative impact on RT enzymatic activity, hindering escalations in dolutegravir concentrations beyond 5–10 nM (Figure4c). Dual selections with dolutegravir and lamivudine showed the appearance of resistance in only one large cluster isolate (14947) with the first appearance of R263K at week 8 preceding the acquisition of M184I/V and M184V at weeks 15 and 23, respectively (Figure4d). Drug concentrations of both drugs could not be increased following the acquisition of R263K.

Genotypic analysis of viral outgrowth under dolutegravir and elvitegravir pressure

HIV-1 genotypic analysis based on deep sequencing was performed at selected passages to monitor temporal dynamics of the emergence and accumulation of mutations conferring resistance to dolutegravir in eight large cluster strains, as compared with six small cluster variants (Figures5 and6).

Figure 5

Viral outgrowth of HIV-1 large cluster strains cultured in the presence of increasing concentrations of dolutegravir. HIV-1 isolates from eight large clusters (cluster size in parentheses) were cultured in the presence of increasing concentrations of dolutegravir. Cell-free supernatant was collected at selected time periods and viral RNA purified and genotyped across the integrase or RT domain. The frequency (%) of individual amino acid substitutions acquired in the virus population at the designated weekly passage was determined based on deep sequencing. Cluster 14947 (a and i) was selected with both dolutegravir and dolutegravir + lamivudine to ascertain the sequential appearance of R263K and M184V. DTG, dolutegravir; 3TC, lamivudine. An asterisk indicates Sanger sequencing only.

Figure 6

Viral outgrowth of HIV-1 singleton/small cluster strains cultured in the presence of increasing concentrations of dolutegravir. HIV-1 isolates from six singleton/small clusters (cluster size in parentheses) were cultured in the presence of increasing concentrations of dolutegravir. Cell-free supernatant was collected at selected time periods and viral RNA purified and genotyped across the integrase domain. The frequency (%) of individual amino acid substitutions acquired in the virus population at the designated weekly passage was determined based on deep sequencing. DTG, dolutegravir.

As depicted in Figure5, dolutegravir selections, performed on large cluster viruses, led to rapid appearance of R263K (n = 6/8) or S153Y (n = 1/8) as dominant members of quasi-species within 6–8 weeks, while H51Y (n = 1/8) appeared by week 12 (Figure5). In general, viral outgrowth under dolutegravir pressure resulted in R263K, S153Y or H51Y as singlet mutations, representing the dominant quasi-species (>96%) at weeks 24–36 based on ultradeep sequencing. The infrequent appearance of S230R, E157Q, G163R or G59E minority species (<10%) occurred in some large cluster selections at selected weeks of passage but failed to accumulate over time. Although isolate 14637, belonging to cluster 185, selected for S153Y (Figure 5b), three subsequent dolutegravir selections using this isolate resulted in R263K as the sole mutation (B. G. Brenner and M. Oliveira, unpublished results). Isolate 10387 showed the dual appearance of R263K and S153F as 52% and 10.9% of the quasi-species at week 8, although R263K (99%) became and remained the dominant quasi-species at week 26 (Figure5d).

There was a delay in the selection of dolutegravir resistance in small cluster as compared with large cluster lineages (Figure5). Four of six small cluster lineages (cluster sizes 1, 1, 1 and 2) failed to select any resistance to dolutegravir over 36 weeks. INSTI-associated mutations were present at really low frequency as determined by deep sequencing (R263K at 1.1% and Y143H at 1.6% in 14514 and 14380, respectively) (Figure5c and d). The remaining two small clusters (cluster sizes 2 and 4) selected R263K or S153Y as singleton mutations (Figure5a and b).

Incremental increases in dolutegravir concentrations (beyond 10 nM) following acquisition of R263K, S153Y and H51Y mutations could not be attempted because the pressure on viral survival was too great based on weekly RT assays (Figures4c and5). This decline in viral replicative fitness under drug pressure is consistent with cell culture selection studies and integrase enzymatic assays previously performed by our group using site-directed mutagenesis to insert resistance substitutions into the laboratory HIV-1 pNL4.3 strain. Notably, isolate 5326 acquiring S153Y showed no further accumulation of mutations after 128 weeks of culture in the presence of dolutegravir (Figure5b).

Cluster 185 isolates 14947 and 14637 harbouring R263K and S153Y at weeks 30 and 34, respectively, were grown in the absence of dolutegravir pressure for a subsequent 19 and 17 weeks, respectively (Figure4a and b). The release of drug pressure failed to lead to a reversion of R263K or S153Y, substantiating their presence as dominant quasi-species.

Interestingly, culturing HIV-1 isolate 14947 (cluster 185) with dolutegravir led to the rapid selection of variants carrying the R263K mutation, as rapidly as M184V selection with lamivudine (Figure4a and d). Again, the acquisition of M184V/I was faster for large cluster as compared with small cluster lineages. It is also significant that R263K arose more rapidly than M184V in dolutegravir + lamivudine dual selections performed with isolate 14947 and that both mutations can coexist (Figures4d and5i).

The development of resistance under elvitegravir pressure proceeded rapidly for large cluster lineages with a sequential accumulation of multiple mutations along different resistance pathways (Figure7a–d). Isolate 14947 acquired T66I and R263K within 8 weeks, representing 98.7% and 44.9% of the quasi-species, respectively (Figure6a). By week 24, the dominant virus coexpressed T66I (99.5%), R263K (98.7%) and E157Q (94.1%) (Figure6a). The second isolate, belonging to cluster 185 (14637), developed resistance with viruses coexpressing N155H (99.5%), Q95K (97.5%) and S230R (99.6%) as the dominant quasi-species (Figure6b). Isolate 14997 acquired a mixture of different viruses with T66I, E92G, E157Q, H51Y and/or S147G representing 85.9%, 19.8%, 99.7%, 15.0% and 9.9% of variants at week 36 (Figure6c). Isolate 14969 acquired S147G (92.1%) at week 8 followed by Q148R (99.7%), known to confer high levels of cross-resistance to elvitegravir and raltegravir (Figure7d). The accumulation of mutations along different pathways contributed to viral breakthrough and >100-fold resistance. Similarly, the acquisition of T66I followed by secondary mutations, including L74M, Q146K, and G163R, by small cluster viruses contributed to viral escape by week 24 (Figure7e and f).

Figure 7

Genotyping analysis based on Sanger and deep sequencing of HIV-1 small cluster strains cultured in the presence of increasing concentrations of elvitegravir. HIV-1 isolates from four large clusters and two singleton/small clusters were cultured in the presence of increasing concentrations of elvitegravir. Cell-free supernatant was collected at selected time periods and viral RNA purified and the RT- and/or integrase-coding region sequenced using Sanger or deep sequencing as previously described. The frequency (%) of individual amino acid substitutions acquired in the virus population at the designated weekly passage was determined based on deep sequencing. EVG, elvitegravir.

It was rather surprising that R263K and H51Y did not adversely impact on viral replicative competence with elvitegravir as observed for dolutegravir (Figure6a). Deep sequencing-based genotyping showed T66I selection occurred prior to R263K followed by E157Q, which may have compensated for the fitness compromise of R263K (Figure6a). The presence of secondary (accessory) mutations, including T66I, E157Q, G163R and S230R, can increase levels of resistance and/or viral fitness (Figures4b and6).

Discussion

Genetic diversity is the hallmark feature of HIV-1 contributing to the divergent spread of distinct HIV-1 subtypes and recombinant forms in different regional settings. In Quebec, the subtype B HIV-1 epidemic among MSM has not declined in the era of potent and simplified ART and improved treatment-as-prevention strategies. Viral phylogenetic cluster analysis has shown that half of the subtype B HIV-1 epidemic among MSM can be ascribed to the episodic genesis and expansion of 30 viral lineages, contributing to 13%, 25% and 42% of new diagnoses over the periods 2004–07, 2008–11 and 2012–15, respectively., Here we proposed that the disproportionate contribution of these 30 strains over more than a decade may be related to distinct virological features beyond the effects of risk behaviour and the high infectiousness of acute/recent stage infection. Novel findings observed in this study usingin vitro selections revealed that HIV-1 large cluster lineages show a lower barrier to resistance when compared with viruses from singleton/small cluster networks.

The rapid acquisition of R263K, S153Y or H51Y with dolutegravir was unexpected given the isolated number of reported cases of resistance in the clinic. Indeed, the appearance of these dolutegravir-related mutations within 6–12 weeks using HIV-1 strains from large clusters was far more rapid than previous studies by our group using laboratory strains where resistance arises after 20 weeks of culture with dolutegravir., Furthermore, viruses obtained by site-directed mutagenesis of the laboratory strain pNL4.3 with R263K or H51Y are severely compromised in their ability to acquire or coexist with other mutations, such as M184V and T66I. Findings reported herein revealed the emergence of viral variants coexpressing R263K/M184V and R263K/T66I as dominant (>98%) quasi-species under selective pressure with dolutegravir + lamivudine and elvitegravir, respectively. The acquisition of E157Q (94% and 99%) by two viruses (and) coexpressing T66I/R263K (>98%, isolate 14947) or a mixed T66I (85.9%) and H51Y (15.9%) quasi-species (isolate 14997) under elvitegravir pressure suggests that E157Q may serve as an accessory mutation that restores viral replicative fitness and increases resistance.

Given the high potency of dolutegravir in the clinic, it has proven difficult to document and characterize dolutegravir resistance profiles. Cell culture selections performed using replicatively resilient large cluster variants highlighted novel aspects regarding emergent resistance to dolutegravir that could not be otherwise ascertained. We observed that the rate of appearance of R263K with dolutegravir was as rapid as that observed for the appearance of M184V with lamivudine. The emergence of select mutations under dolutegravir pressure, including R263K (AGA→AAA or AGG→AAG), G163R (GGA→AGA or GGG→AGG) and E138K (GAA→AAA, GAG→AAG), may be related to APOBEC-mediated G→A hypermutation., Our previous studies also showed the development of G118R (GGA→AGA) in two patients failing dolutegravir monotherapy. These patients harboured a rare GGA natural polymorphism (1.5% prevalence) compared with the majority of clinical viral isolates harbouring GGC (79%) or GGT (18%) glycine motifs where G118R acquisition would require two point mutations (GGG→AGA or GGU→AGA) or single unfavourable transversions (GGG→CGG orGGU→CGT), respectively. The generation of M184I (ATG→ATA) can also arise from G to A hypermutation although M184V leads to higher fold drug resistance. Although APOBEC3G-mediated G to A hypermutation may contribute to the development of drug resistance, it can also inhibit viral replication by driving viruses into error catastrophe mode.

The development of resistance to antiretroviral drug regimensin vivo is predicated on the degree and fitness costs conferred by individual mutations, alone and in combination with other resistance mutations. Dolutegravir is a drug for which thein vitro acquisition of singlet substitutions, including R263K, S153Y and H51Y, interfered with further acquisition of secondary resistance mutations., These mutations confer minimal resistance (<3-fold) but impose enormous constraints on viral replicative fitness, decreasing viral infectiousness, replication capacity, DNA integration and integrase strand-transfer activity.,, The high potency for dolutegravir resulted in EC50s of 1 nM (0.4 ng/ml) and drug concentrations could not be increased beyond 5 to 10 nM with R263K or S153Y to sustainin vitro viral replication in CBMCs.,, Thesein vitro drug concentrations are consistent with therapeutic concentrations of unbound (17 ng/mL) and bound (3400 ng/mL) dolutegravir in plasma., This contrasts with elvitegravir selections, such as the cluster 185 isolate, 14947, where acquisition of T66I preceded R263K, resulting in higher levels of resistance, the acquisition of the compensatory E157Q and viral escape. The combination of R263K with viruses containing M184I or M184V mutations in dolutegravir + lamivudine selections further reduced viral replication competence compared with when only single mutations were present. Taken together, ourin vitro findings incorporating magnitude of resistance and fitness costs of individual mutations provide an explanation for the remarkable clinical success of dolutegravir and bictegravir in recent clinical trials. The high potency of second-generation integrase inhibitors may prevent or attenuate the development of replicatively fit large cluster lineages.

Large cluster outbreaks can contribute to the onward spread of transmitted drug resistance, principally involving NNRTIs., Indeed, the frequency of these resistance mutations in treatment-naive persons has surpassed that observed in persons failing therapy since 2011., Notably, cluster 50 (cluster size 140) and cluster 99-K103N (cluster size 34) remain persistent sub-epidemics harbouring G190A and K103N, respectively. Genetic differences between viral strains may contribute to their relative transmissibility, replicative competence and propensity towards the development of resistance. Our findings show that R263K and S153Y, once acquired, may not revert. The negative viral fitness costs of these mutations may, however, limit their development and potential transmissibility. In this regard the prevalence of resistance in patients failing integrase-based regimens was 8% in Quebec (January 2014–April 2016), largely associated with raltegravir- and elvitegravir-associated resistance mutations (I. Hardy and B. G. Brenner, unpublished results). Resistance to dolutegravir has, however, been limited to the G118R case and a novel resistance pathway from a raltegravir-experienced patient harbouring the N155H mutation.,

Bottlenecks constrain HIV-1 transmission, such that a single monophyletic founder variant is responsible for establishing individual infections. Studies have compared the ‘signature’ genetic motifs and phenotypic properties of newly transmitted viruses to those of chronic-stage viruses., Transmitted viruses are most often CCR5-tropic, possess signature leader sequences (codon 11), show less envelope glycosylation and have shorter variable envelope loops than chronic stage viruses., Here, we also show an accelerated development of resistance to integrase and RT inhibitors on the part of viruses belonging to 10 different large 20+ clusters, perhaps suggesting a potential replicative and/or enzymatic fitness advantage in the subset of large cluster viruses. Studies are ongoing to ascertain RT and integrase activity for WT and resistant (R263K or S153Y) cluster 185 isolates, compared with pNL4.3. Direct study of the infection and replication properties of clinically relevant large cluster viruses in cell culture may help define distinct features of transmitted viruses and early events following infection.

In summary, viral phylogenetic cluster analysis of the HIV-1 epidemic among MSM in Quebec shows an evolutionary shift towards large cluster lineages. These viruses appear to be able to overcome severe transmission bottlenecks and circumvent genetic barriers to resistance to dolutegravir.,, These findings have implications for novel treatment paradigms under investigation, including simplified dolutegravir + lamivudine dual therapy and dolutegravir monotherapy regimens., Large cluster viruses may assist in assessing as yet unknown pathways implicated in emergent resistance to new second-generation INSTIs in clinical trial, including bictegravir and cabotegravir.,,