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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Large cluster outbreaks sustain the HIV epidemic among MSM in Quebec.

OBJECTIVE: HIV-1 epidemics among MSM remain unchecked despite advances in treatment and prevention paradigms. This study combined viral phylogenetic and behavioural risk data to better understand underlying factors governing the temporal growth of the HIV epidemic among MSM in Quebec (2002-2015). METHODS: Phylogenetic analysis of pol sequences was used to deduce HIV-1 transmission dynamics (cluster size, size distribution and growth rate) in first genotypes of treatment-naïve MSM (2002-2015, n = 3901). Low sequence diversity of first genotypes (0-0.44% mixed base calls) was used as an indication of early-stage infection. Behavioural risk data were obtained from the Montreal rapid testing site and primary HIV-1-infection cohorts. RESULTS: Phylogenetic analyses uncovered high proportion of clustering of new MSM infections. Overall, 27, 45, 48, 53 and 57% of first genotypes within one (singleton, n = 1359), 2-4 (n = 692), 5-9 (n = 367), 10-19 (n = 405) and 20+ (n = 1277) cluster size groups were early infections (<0.44% diversity). Thirty viruses within large 20+ clusters disproportionately fuelled the epidemic, representing 13, 25 and 42% of infections, first genotyped in 2004-2007 (n = 1314), 2008-2011 (n = 1356) and 2012-2015 (n = 1033), respectively. Of note, 35, 21 and 14% of MSM belonging to 20+, 2-19 and one (singleton) cluster groups were under 30 years of age, respectively. Half of persons seen at the rapid testing site (2009-2011, n = 1781) were untested in the prior year. Poor testing propensity was associated with fewer reported partnerships. CONCLUSION: Addressing the heterogeneity in transmission dynamics among HIV-1-infected MSM populations may help guide testing, treatment and prevention strategies.

Development of a G118R mutation in HIV-1 integrase following a switch to dolutegravir monotherapy leading to cross-resistance to integrase inhibitors.

OBJECTIVES: Dolutegravir shows a high barrier to resistance with no previously reported cases of acquired integrase mutations during first-line therapy. In this study, rapid development of the G118R mutation arose following a switch from first-line elvitegravir/cobicistat/tenofovir disoproxil fumarate/emtricitabine to dolutegravir monotherapy. The G118R mutation also arose in a treatment-experienced patient switched to dolutegravir monotherapy. The genetic basis for G118R selection and potential phenotypic outcome was ascertained. PATIENT AND METHODS: Genotypic analysis was performed on patients with virological failure (<1000 copies/mL) on dolutegravir-containing regimens. The Los Alamos database was queried for glycine codon 118 polymorphisms. Cell culture selections and phenotypic drug susceptibility assays assessed resistance via the G118R pathway. RESULTS: We report on two patients who developed viral failure while on dolutegravir monotherapy. Both patients had been on a current or previous regimen containing integrase inhibitors. Virological failure (<1000 copies/mL) emerged early within 2 months following the dolutegravir switch. The appearance of G118R in these two cases and subtype C and CRF02_AG in vitro selections were related to a rare GGA natural polymorphism at codon 118 (1.5% prevalence), facilitating a GGA to AGA transition. Cell culture selections were used to assess the in vitro progression of the G118R pathway leading to cross-resistance to all integrase inhibitors. CONCLUSIONS: Although resistance to dolutegravir is typically rare, genetic polymorphisms and monotherapy can facilitate the acquisition of G118R.

The dolutegravir R263K resistance mutation in HIV-1 integrase is incompatible with the emergence of resistance against raltegravir.

OBJECTIVE: Although the integrase inhibitor dolutegravir (DTG) has demonstrated greater resilience than other antiretroviral drugs at withstanding the emergence of HIV-1 resistance mutations, such substitutions can develop, albeit rarely, in treatment-experienced integrase inhibitor-naïve individuals. The most common substitution in integrase under those circumstances is R263K whereas another substitution that was selected against DTG in tissue culture was G118R. The objective of this study was to determine the effects of these DTG-specific resistance substitutions on the ability of HIV-1 to become resistant against either of two other integrase inhibitors, raltegravir (RAL) and elvitegravir (EVG). DESIGN AND METHODS: We performed tissue culture selection experiments using DTG-resistant viruses containing integrase substitutions at positions R263K, H51Y/R263K, E138K/R263K, G118R and H51Y/G118R in the presence of increasing concentrations of either RAL or EVG. Changes in integrase sequences were monitored by genotyping. RESULTS: The presence of the R263K substitution delayed the emergence of resistance against RAL whereas the simultaneous presence of either the H51Y or E138K secondary substitutions in combination with R263K somewhat mitigated this inhibitory effect. In contrast, resistance against EVG appeared earlier than in wild-type virus in viruses containing the R263K and E138K/R263K DTG-associated resistance substitutions. CONCLUSION: The DTG-resistant R263K substitution antagonized the development of HIV-1 resistance against RAL while partially facilitating the occurrence of resistance against EVG.

Basis for early and preferential selection of the E138K mutation in HIV-1 reverse transcriptase.

E138K, a G→A mutation in HIV-1 reverse transcriptase (RT), is preferentially selected by etravirine (ETR) and rilpivirine over other substitutions at position E138 that offer greater drug resistance. We hypothesized that there was a mutational bias for the E138K substitution and designed an allele-specific PCR to monitor the emergence of E138A/G/K/Q/R/V during ETR selection experiments. We also performed competition experiments using mutated viruses and quantified the prevalence of E138 minority species in drug-naive patients. E138K, as well as E138G, consistently emerged first during ETR selection experiments, followed by E138A and E138Q; E138R was never selected. Surprisingly, E138K was identified as a tiny minority in 23% of drug-naive subtype B patients, a result confirmed by ultradeep sequencing (UDS). This result could reflect a low fitness cost of E138K; however, E138K was one of the least fit substitutions at codon E138, even after taking into account the deoxynucleoside triphosphate pools of the cells used in competition experiments. Further UDS analysis revealed other minority species in a pattern consistent with the mutational bias of HIV RT. There was no evidence of APOBEC3-hypermutation in these selection experiments or in patients. Our results confirm the mutational bias of HIV-1 in patients and highlight the importance of G→A mutations in HIV-1 drug resistance evolution.

HIV gp120 H375 is unique to HIV-1 subtype CRF01_AE and confers strong resistance to the entry inhibitor BMS-599793, a candidate microbicide drug.

BMS-599793 is a small molecule entry inhibitor that binds to human immunodeficiency virus type 1 (HIV-1) gp120, resulting in the inhibition of CD4-dependent entry into cells. Since BMS-599793 is currently considered a candidate microbicide drug, we evaluated its efficacy against a number of primary patient HIV isolates from different subtypes and circulating recombinant forms (CRFs) and showed that activity varied between ∼3 ρM and 7 µM at 50% effective concentrations (EC(50)s). Interestingly, CRF01_AE HIV-1 isolates consistently demonstrated natural resistance against this compound. Genotypic analysis of >1,600 sequences (Los Alamos HIV sequence database) indicated that a single amino acid polymorphism in Env, H375, may account for the observed BMS-599793 resistance in CRF01_AE HIV-1. Results of site-directed mutagenesis experiments confirmed this hypothesis, and in silico drug docking simulations identified a drug resistance mechanism at the molecular level. In addition, CRF01_AE viruses were shown to be resistant to multiple broadly neutralizing monoclonal antibodies. Thus, our results not only provide insight into how Env polymorphisms may contribute to entry inhibitor resistance but also may help to elucidate how HIV can evade some broadly neutralizing antibodies. Furthermore, the high frequency of H375 in CRF01_AE HIV-1, and its apparent nonoccurrence in other subtypes, could serve as a means for rapid identification of CRF01_AE infections.

In vitro resistance profile of the candidate HIV-1 microbicide drug dapivirine.

Antiretroviral-based microbicides may offer a means to reduce the sexual transmission of HIV-1. Suboptimal use of a microbicide may, however, lead to the development of drug resistance in users that are already, or become, infected with HIV-1. In such cases, the efficacy of treatments may be compromised since the same (or similar) antiretrovirals used in treatments are being developed as microbicides. To help predict which drug resistance mutations may develop in the context of suboptimal use, HIV-1 primary isolates of different subtypes and different baseline resistance profiles were used to infect primary cells in vitro in the presence of increasing suboptimal concentrations of the two candidate microbicide antiretrovirals dapivirine (DAP) and tenofovir (TFV) alone or in combination. Infections were ongoing for 25 weeks, after which reverse transcriptase genotypes were determined and scrutinized for the presence of any clinically recognized reverse transcriptase drug resistance mutations. Results indicated that suboptimal concentrations of DAP alone facilitated the emergence of common nonnucleoside reverse transcriptase inhibitor resistance mutations, while suboptimal concentrations of DAP plus TFV gave rise to fewer mutations. Suboptimal concentrations of TFV alone did not frequently result in the development of resistance mutations. Sensitivity evaluations for stavudine (d4T), nevirapine (NVP), and lamivudine (3TC) revealed that the selection of resistance as a consequence of suboptimal concentrations of DAP may compromise the potential for NVP to be used in treatment, a finding of potential relevance in developing countries.