Checklist for studies of HIV drug resistance prevalence or incidence: rationale and recommended use.

HIV drug resistance (HIVDR) is a major challenge to the effectiveness of antiretroviral therapy. Global efforts in addressing HIVDR require clear, transparent, and replicable reporting in HIVDR studies. We describe the rationale and recommended use of a checklist that should be included in reports of HIVDR incidence and prevalence. After preliminary consultations with experts on HIVDR and establishing the need for guidance on HIVDR reporting, we used a sequential, explanatory, mixed methods approach to create the checklist; together with the accompanying articles, the checklist was reviewed by the authors and validated externally. The checklist for studies on HIVDR prevalence or incidence (CEDRIC-HIV) includes 15 recommended items that would enhance transparency and facilitate interpretation, comparability, and replicability of HIVDR studies. CEDRIC-HIV will help authors of HIVDR studies prepare research reports and assist reviewers and editors in assessments of completeness of reporting. The checklist will also facilitate statistical pooling and interpretation of HIVDR data.

Doravirine responses to HIV-1 viruses bearing mutations to NRTIs and NNRTIs under in vitro selective drug pressure.

OBJECTIVES: The NNRTI doravirine has been recently approved for the first-line treatment of HIV-infected patients, eliciting favourable responses against viruses bearing the K103N, Y181C and G190A mutations. This study used in vitro drug selections to elaborate the breadth of doravirine responses against viruses bearing NNRTI and NRTI resistance-associated mutations (RAMs). METHODS: WT clinical isolates (n = 6) and viruses harbouring common NRTI and NNRTI RAMs (n = 6) were serially passaged in escalating concentrations of doravirine, doravirine/islatravir, doravirine/lamivudine and rilpivirine over 24 weeks. Genotypic analysis ascertained the appearance and accumulation of NNRTI RAMs. Phenotypic drug susceptibility assays assessed resistance conferred by acquired NNRTI RAMs. RESULTS: For WT viruses, doravirine pressure led to the appearance of V108I or V106A/I/M RAMs after 8 weeks, conferring low-level (∼2-fold) resistance. After 24 weeks, the accumulation of three to six secondary RAMs, including F227L, M230L, L234I and/or Y318, resulted in high-level (>100-fold) resistance to doravirine. Notably, viruses with these doravirine RAMs remained susceptible to rilpivirine and efavirenz. This contrasted with rilpivirine where acquisition of E138K, L100I and/or K101E resulted in >50-fold cross-resistance to all NNRTIs. Doravirine selection of viruses bearing common NRTI and NNRTI RAMs showed delayed acquisition of RAMs compared with WT virus. Pairing doravirine with islatravir or lamivudine attenuated the development of NNRTI RAMs. CONCLUSIONS: Doravirine showed favourable resistance profiles against viruses harbouring NRTI and NNRTI RAMs. The high barrier to resistance to doravirine coupled with the long intracellular half-life of islatravir may provide the opportunity for long-acting treatment options.

Long Dissociation of Bictegravir from HIV-1 Integrase-DNA Complexes.

The HIV integrase (IN) strand transfer inhibitor (INSTI) bictegravir (BIC) has a long dissociation half-life (t1/2) from wild-type IN-DNA complexes: BIC 163 hr > dolutegravir (DTG) 96 hr > raltegravir (RAL) 10 hr > elvitegravir (EVG) 3.3 hr. In cells, BIC had more durable antiviral activity against wild-type HIV after drug washout than RAL or EVG. BIC also had a longer t1/2 and maintained longer antiviral activity after drug washout than DTG with the clinically relevant resistance IN mutant G140S+Q148H. Structural analyses indicate that BIC makes more contacts with the IN-DNA complex than DTG mainly via its bicyclic ring system which may contribute to more prolonged residence time and resilience against many resistance mutations.

Phylogenetic Insights on Patterns of HIV-1 Spread and the Design of Epidemic Control Measures.

Phylogenetics provides a unique structural framework to track the spread of viral diseases, such as HIV-1 (human immunodeficiency virus type 1), the causative agent of AIDS (acquired immunodeficiency syndrome) [...].

The Role of Phylogenetics in Unravelling Patterns of HIV Transmission towards Epidemic Control: The Quebec Experience (2002-2020).

Phylogenetics has been advanced as a structural framework to infer evolving trends in the regional spread of HIV-1 and guide public health interventions. In Quebec, molecular network analyses tracked HIV transmission dynamics from 2002-2020 using MEGA10-Neighbour-joining, HIV-TRACE, and MicrobeTrace methodologies. Phylogenetics revealed three patterns of viral spread among Men having Sex with Men (MSM, n = 5024) and heterosexuals (HET, n = 1345) harbouring subtype B epidemics as well as B and non-B subtype epidemics (n = 1848) introduced through migration. Notably, half of new subtype B infections amongst MSM and HET segregating as solitary transmissions or small cluster networks (2-5 members) declined by 70% from 2006-2020, concomitant to advances in treatment-as-prevention. Nonetheless, subtype B epidemic control amongst MSM was thwarted by the ongoing genesis and expansion of super-spreader large cluster variants leading to micro-epidemics, averaging 49 members/cluster at the end of 2020. The growth of large clusters was related to forward transmission cascades of untreated early-stage infections, younger at-risk populations, more transmissible/replicative-competent strains, and changing demographics. Subtype B and non-B subtype infections introduced through recent migration now surpass the domestic epidemic amongst MSM. Phylodynamics can assist in predicting and responding to active, recurrent, and newly emergent large cluster networks, as well as the cryptic spread of HIV introduced through migration.

Cell culture selections reveal favourable drug resistance profiles for doravirine and islatravir.

BACKGROUND: The newer generation NNRTIs, including doravirine and rilpivirine, were designed to show high potency and overcome K103N, Y181C and G190A resistance. OBJECTIVES: To assess emergent resistance to doravirine and rilpivirine, alone and paired with lamivudine or islatravir through in vitro drug selections. METHODS: Subtype B (n = 3), non-B subtype (n = 3), and pNL4.3 viral isolates were passaged in cord blood mononuclear cells with progressively increasing concentrations of drug(s). Genotypic analysis compared the acquisition and accumulation of drug resistance mutations at weeks 8 and 24 following drug pressure. Cell-based phenotypic assays assessed cross-resistance patterns to NNRTIs by acquired resistance mutations. RESULTS: Doravirine pressure resulted in the acquisition of V108I (6/7) and V106A/I/M (5/7) mutations at weeks 8, followed by F227L (4/7), Y318F (4/7), M230L (2/7) or L234I (2/7) by weeks 24. In contrast, rilpivirine resulted in E138K (5/7) followed by L100I (3/7), K101E (1/7), or M230L (1/7). Doravirine resistance pathways retained susceptibility to rilpivirine, whereas rilpivirine resistance conferred intermediate resistance (12-152-fold) to doravirine. Dual selections with islatravir or lamivudine delayed and diminished emergent resistance to doravirine, resulting in V108I (9/15) with fewer or no other changes at weeks 24. There was a lesser delay in emergent resistance to rilpivirine when combined with islatravir or lamivudine. The M184V mutation did not arise in dual selections with islatravir or lamivudine. CONCLUSIONS: Doravirine showed a more robust resistance profile compared with other NNRTIs. The long intracellular half-life of islatravir and delayed acquisition of resistance in dual selections provide an opportunity for long-acting treatment options.

Phylogenetic Clustering among Asylum Seekers with New HIV-1 Diagnoses in Montreal, QC, Canada.

Migrants are at an increased risk of HIV acquisition. We aimed to use phylogenetics to characterize transmission clusters among newly-diagnosed asylum seekers and to understand the role of networks in local HIV transmission. Retrospective chart reviews of asylum seekers linked to HIV care between 1 June 2017 and 31 December 2018 at the McGill University Health Centre and the Jewish General Hospital in Montreal were performed. HIV-1 partial pol sequences were analyzed among study participants and individuals in the provincial genotyping database. Trees were reconstructed using MEGA10 neighbor-joining analysis. Clustering of linked viral sequences was based on a strong bootstrap support (>97%) and a short genetic distance (<0.01). Overall, 10,645 provincial sequences and 105 asylum seekers were included. A total of 13/105 participant sequences (12%; n = 7 males) formed part of eight clusters. Four clusters (two to three people) included only study participants (n = 9) and four clusters (two to three people) included four study participants clustered with six individuals from the provincial genotyping database. Six (75%) clusters were HIV subtype B. We identified the presence of HIV-1 phylogenetic clusters among asylum seekers and at a population-level. Our findings highlight the complementary role of cohort data and population-level genotypic surveillance to better characterize transmission clusters in Quebec.

In-depth Sampling of High-risk Populations to Characterize HIV Transmission Epidemics Among Young MSM Using PrEP in France and Quebec.

BACKGROUND: A better understanding of HIV transmission dynamics among populations at high risk is important for development of prevention strategies. We determined HIV transmission networks from infected individuals enrolled in the pre-exposure prophylaxis (PrEP) IPERGAY trial in combination with the ANRS PRIMO and Montreal PHI cohorts to identify and characterize active clusters of transmission in this high-risk population. METHODS: Genotypic resistance tests were performed on plasma samples from 31 IPERGAY participants. Reverse transcriptase sequences were analyzed in combination with unique HIV pol sequences from 1351 individuals enrolled in the PRIMO ANRS cohort (1999-2014) and 511 individuals enrolled in the Montreal PHI cohort (1996-2016). Network analyses were performed to infer putative relationships between all participants. RESULTS: Overall, 1893 participants were included. Transmission network analyses revealed that 14 individuals (45.2%) from the IPERGAY trial were involved in 13 clusters sampled over a median period (interquartile range) of 2 (0.3-7.8) years, including 7 dyads and 6 larger clusters ranging from 4 to 28 individuals. When comparing characteristics between clustering individuals enrolled in the PRIMO cohort (n = 377) and in IPERGAY (n = 14), we found that IPERGAY participants had a higher viral load (5.93 vs 5.20 log10 copies/mL, P = .032) and reported a higher number of partners in the last 2 months (P < .01). CONCLUSIONS: These results demonstrate high rates of HIV transmission clustering among young high-risk MSM enrolled in the IPERGAY trial. In-depth sampling of high-risk populations may help to uncover unobserved transmission intermediaries and improve prevention efforts that could be targeted to the most active clusters.

Assessing the role of transmission chains in the spread of HIV-1 among men who have sex with men in Quebec, Canada.

BACKGROUND: Phylogenetics has been used to investigate HIV transmission among men who have sex with men. This study compares several methodologies to elucidate the role of transmission chains in the dynamics of HIV spread in Quebec, Canada. METHODS: The Quebec Human Immunodeficiency Virus (HIV) genotyping program database now includes viral sequences from close to 4,000 HIV-positive individuals classified as Men who have Sex with Men (MSMs), collected between 1996 and early 2016. Assessment of chain expansion may depend on the partitioning scheme used, and so, we produce estimates from several methods: the conventional Bayesian and maximum likelihood-bootstrap methods, in combination with a variety of schemes for applying a maximum distance criterion, and two other algorithms, DM-PhyClus, a Bayesian algorithm that produces a measure of uncertainty for proposed partitions, and the Gap Procedure, a fast non-phylogenetic approach. Sequences obtained from individuals in the Primary HIV Infection (PHI) stage serve to identify incident cases. We focus on the period ranging from January 1st 2012 to February 1st 2016. RESULTS AND CONCLUSION: The analyses reveal considerable overlap between chain estimates obtained from conventional methods, thus leading to similar estimates of recent temporal expansion. The Gap Procedure and DM-PhyClus suggest however moderately different chains. Nevertheless, all estimates stress that longer older chains are responsible for a sizeable proportion of the sampled incident cases among MSMs. Curbing the HIV epidemic will require strategies aimed specifically at preventing such growth.

DM-PhyClus: a Bayesian phylogenetic algorithm for infectious disease transmission cluster inference.

BACKGROUND: Conventional phylogenetic clustering approaches rely on arbitrary cutpoints applied a posteriori to phylogenetic estimates. Although in practice, Bayesian and bootstrap-based clustering tend to lead to similar estimates, they often produce conflicting measures of confidence in clusters. The current study proposes a new Bayesian phylogenetic clustering algorithm, which we refer to as DM-PhyClus (Dirichlet-Multinomial Phylogenetic Clustering), that identifies sets of sequences resulting from quick transmission chains, thus yielding easily-interpretable clusters, without using any ad hoc distance or confidence requirement. RESULTS: Simulations reveal that DM-PhyClus can outperform conventional clustering methods, as well as the Gap procedure, a pure distance-based algorithm, in terms of mean cluster recovery. We apply DM-PhyClus to a sample of real HIV-1 sequences, producing a set of clusters whose inference is in line with the conclusions of a previous thorough analysis. CONCLUSIONS: DM-PhyClus, by eliminating the need for cutpoints and producing sensible inference for cluster configurations, can facilitate transmission cluster detection. Future efforts to reduce incidence of infectious diseases, like HIV-1, will need reliable estimates of transmission clusters. It follows that algorithms like DM-PhyClus could serve to better inform public health strategies.

Selective resistance profiles emerging in patient-derived clinical isolates with cabotegravir, bictegravir, dolutegravir, and elvitegravir.

BACKGROUND: Integrase strand transfer inhibitors (INSTIs) are recommended for first-line HIV therapy based on their relatively high genetic barrier to resistance. Although raltegravir (RAL) and elvitegravir (EVG) resistance profiles are well-characterized, resistance patterns for dolutegravir (DTG), bictegravir (BIC), and cabotegravir (CAB) remain largely unknown. Here, in vitro drug selections compared the development of resistance to DTG, BIC, CAB, EVG and RAL using clinical isolates from treatment-naïve primary HIV infection (PHI) cohort participants (n = 12), and pNL4.3 recombinant strains encoding patient-derived Integrase with (n = 5) and without (n = 5) the E157Q substitution. RESULTS: Patient-derived viral isolates were serially passaged in PHA-stimulated cord blood mononuclear cells in the presence of escalating concentrations of INSTIs over the course of 36-46 weeks. Drug resistance arose more rapidly in primary clinical isolates with EVG (12/12), followed by CAB (8/12), DTG (8/12) and BIC (6/12). For pNL4.3 recombinant strains encoding patient-derived integrase, the comparative genetic barrier to resistance was RAL > EVG > CAB > DTG and BIC. The E157Q substitution in integrase delayed the advent of resistance to INSTIs. With EVG, T66I/A, E92G/V/Q, T97A or R263K (n = 16, 3, 2 and 1, respectively) arose by weeks 8-16, followed by 1-4 accessory mutations, conferring high-level resistance (> 100-fold) by week 36. With DTG and BIC, solitary R263K (n = 27), S153F/Y (n = 7) H51Y (n = 2), Q146 R (n = 3) or S147G (n = 1) mutations conferred low-level (< 3-fold) resistance at weeks 36-46. Similarly, most CAB selections (n = 18) resulted in R263K, S153Y, S147G, H51Y, or Q146L solitary mutations. However, three CAB selections resulted in Q148R/K followed by secondary mutations conferring high-level cross-resistance to all INSTIs. EVG-resistant viruses (T66I/R263K, T66I/E157Q/R263K, and S153A/R263K) retained residual susceptibility when switched to DTG, BIC or CAB, losing T66I by week 27. Two EVG-resistant variants developed resistance to DTG, BIC and CAB through the additional acquisition of E138A/Q148R and S230N, respectively. One EVG-resistant variant (T66I) acquired L74M/G140S/S147G, L74M/E138K/S147G and H51Y with DTG CAB and BIC, respectively. CONCLUSIONS: Second generation INSTIs show a higher genetic barrier to resistance than EVG and RAL. The potency of CAB was lower than BIC and DTG. The development of Q148R/K with CAB can result in high-level cross-resistance to all INSTIs.

Durable suppression of HIV-1 with resistance mutations to integrase inhibitors by dolutegravir following drug washout.

OBJECTIVES: Dolutegravir (DTG) has achieved better long-term suppression of HIV-1 replication than other integrase strand transfer inhibitors (INSTIs), such as raltegravir (RAL) and elvitegravir (EVG). In in-vitro drug washout experiments, we previously showed that removal of DTG from pretreated MT-2 cells infected with wild-type HIV-1 showed slower rebound in viral replication as compared to removal of RAL. Now, we performed DTG, EVG and RAL washout experiments to compare the recovery of viral integration and production of 2-long terminal repeat (LTR) circles using wild-type HIV-1 clones, R263K viruses with low-level resistance to DTG and viruses with G140S/Q148H mutations showing cross-resistance against all currently approved INSTIs. DESIGN AND METHODS: MT-2 cells infected with wild-type, R263K or G140S/Q148H HIV-1 clones were treated with DTG, RAL or EVG for 3 days. Viral rebound following drug washout was assessed, monitoring viral integration and 2-LTR circle production by qPCR. RESULTS: Viral integration did not resume for up to 8 days after DTG washout from the wild-type or R263K infections but increased soon after washout of either RAL or EVG. With the G140S/Q148H virus, levels of integration were not significantly affected by the presence of either RAL or EVG. With DTG, integration was much lower at 3 days after infection than for the no-drug control. At 8 days after DTG washout, viral integration resumed but remained relatively low. CONCLUSION: DTG antiretroviral activity in tissue culture is more durable than that of either RAL or EVG after drug washout and this is true for both wild-type and drug-resistant viruses.

The S230R Integrase Substitution Associated With Virus Load Rebound During Dolutegravir Monotherapy Confers Low-Level Resistance to Integrase Strand-Transfer Inhibitors.

Background: Dolutegravir (DTG) is an integrase strand-transfer inhibitor (INSTI) used for treatment of human immunodeficiency virus (HIV)-infected individuals. Owing to its high genetic barrier to resistance, DTG has been clinically investigated as maintenance monotherapy to maintain viral suppression and to reduce complication and healthcare costs. Our study aims to explain the underlying mechanism related to the emergence of a S230R substitution in patients who experienced virologic failure while using DTG monotherapy. Methods: We evaluated the effect of the S230R substitution in regard to integrase enzyme activity, viral infectivity, replicative capacity, and susceptibility to different INSTIs by biochemical and cell-based assays. Results: The S230R substitution conferred a 63% reduction in enzyme efficiency. S230R virus was 1.29-fold less infectious than wild-type virus but could replicate in PM1 cells without significant delay. Resistance levels against DTG, cabotegravir, raltegravir, and elvitegravir in tissue culture were 3.85-, 3.72-, 1.52-, and 1.21-fold, respectively, in virus with the S230R substitution. Conclusions: Our data indicate that the S230R substitution is comparable to the previously reported R263K substitution in some respects. Virologic failure during DTG monotherapy can occur through the development of the S230R or R263K mutation, without the need for high-level DTG resistance.

Genotypic and Phylogenetic Insights on Prevention of the Spread of HIV-1 and Drug Resistance in "Real-World" Settings.

HIV continues to spread among vulnerable heterosexual (HET), Men-having-Sex with Men (MSM) and intravenous drug user (IDU) populations, influenced by a complex array of biological, behavioral and societal factors. Phylogenetics analyses of large sequence datasets from national drug resistance testing programs reveal the evolutionary interrelationships of viral strains implicated in the dynamic spread of HIV in different regional settings. Viral phylogenetics can be combined with demographic and behavioral information to gain insights on epidemiological processes shaping transmission networks at the population-level. Drug resistance testing programs also reveal emergent mutational pathways leading to resistance to the 23 antiretroviral drugs used in HIV-1 management in low-, middle- and high-income settings. This article describes how genotypic and phylogenetic information from Quebec and elsewhere provide critical information on HIV transmission and resistance, Cumulative findings can be used to optimize public health strategies to tackle the challenges of HIV in "real-world" settings.

Antiviral Activity of Bictegravir and Cabotegravir against Integrase Inhibitor-Resistant SIVmac239 and HIV-1.

Animal models are essential to study novel antiretroviral drugs, resistance-associated mutations (RAMs), and treatment strategies. Bictegravir (BIC) is a novel potent integrase strand transfer inhibitor (INSTI) that has shown promising results against HIV-1 infection in vitro and in vivo and against clinical isolates with resistance against INSTIs. BIC has a higher genetic barrier to the development of resistance than two clinically approved INSTIs, termed raltegravir and elvitegravir. Another clinically approved INSTI, dolutegravir (DTG) also possesses a high genetic barrier to resistance, while a fourth compound, termed cabotegravir (CAB), is currently in late phases of clinical development. Here we report the susceptibilities of simian immunodeficiency virus (SIV) and HIV-1 integrase (IN) mutants containing various RAMs to BIC, CAB, and DTG. BIC potently inhibited SIV and HIV-1 in single cycle infection with 50% effective concentrations (EC50s) in the low nM range. In single cycle SIV infections, none of the E92Q, T97A, Y143R, or N155H substitutions had a significant effect on susceptibility to BIC (≤4-fold increase in EC50), whereas G118R and R263K conferred ∼14-fold and ∼6-fold increases in EC50, respectively. In both single and multiple rounds of HIV-1 infections, BIC remained active against the Y143R, N155H, R263K, R263K/M50I, and R263K/E138K mutants (≤4-fold increase in EC50). In multiple rounds of infection, the G140S/Q148H combination of substitutions decreased HIV-1 susceptibility to BIC 4.8-fold compared to 16.8- and 7.4-fold for CAB and DTG, respectively. BIC possesses an excellent resistance profile in regard to HIV and SIV and could be useful in nonhuman primate models of HIV infection.

M184I/V substitutions and E138K/M184I/V double substitutions in HIV reverse transcriptase do not significantly affect the antiviral activity of EFdA.

BACKGROUND: 4'-Ethynyl-2-fluoro-2'-deoxyadenosine (EFdA) is a potent nucleoside analogue inhibitor of HIV that has an unusually long half-life. Cell culture selections with either EFdA or lamivudine using both subtype B and non-B clinical isolates selected the M184I/V substitutions in reverse transcriptase (RT). Unlike lamivudine, however, EFdA retained significant activity against viruses containing the M184I/V substitutions. In other clinical trials that evaluated rilpivirine together with emtricitabine in first-line therapy, the emergence of both the M184I/V and E138K mutations in RT was demonstrated. Moreover, the M184I/V and E138K substitutions were shown to be compensatory for each other with regard to both efficiency of RT activity and viral replicative capacity. This creates concern that E138K might emerge as a compensatory mutation for M184I/V in the aftermath of the use of EFdA. OBJECTIVES: We wished to determine whether the E138K mutation in HIV RT together with M184I/V would compromise the activity of EFdA. METHODS: Recombinant viruses containing the M184I/V and/or E138K substitutions were generated by site-directed mutagenesis and evaluated in tissue culture for susceptibility to various nucleoside compounds, including EFdA. RESULTS: Susceptibility to EFdA was retained in M184I/V viruses that also contained the E138K substitution. Moreover, the E138K substitution was not generated in these studies under selection pressure with EFdA. CONCLUSIONS: These findings help to alleviate concern that EFdA may not be active against viruses that contain both the M184I/V and E138K substitutions in RT.