Factors limiting the transmission of HIV mutations conferring drug resistance: fitness costs and genetic bottlenecks

Trends of pre-treatment drug resistance in antiretroviral-naïve people with HIV-1 in the era of second-generation integrase strand-transfer inhibitors in Taiwan

BACKGROUND

Monitoring the trends of pre-treatment drug resistance (PDR) and resistance-associated mutations (RAMs) among antiretroviral-naïve people with HIV (PWH) is important for the implementation of HIV treatment and control programmes. We analysed the trends of HIV-1 PDR after the introduction of second-generation integrase strand-transfer inhibitors (INSTIs) in 2016 in Taiwan, when single-tablet regimens of non-nucleoside reverse-transcriptase inhibitor (NNRTI-) and INSTI-based antiretroviral therapy became the preferred treatments.

MATERIALS AND METHODS

In this multicentre study, we included newly diagnosed, antiretroviral-naïve PWH who underwent tests for RAMs between 2016 and 2022. Pre-treatment genotypic resistance testing was performed, along with HIV-1 subtyping and determinations of plasma HIV RNA load and CD4 lymphocyte counts. RAMs were analysed using the Stanford University HIV Drug Resistance Database and only RAMs conferring at least low-level resistance were included.

RESULTS

From 2016 to 2022, pre-treatment blood samples from 3001 newly diagnosed PWH, which constituted 24.3% of newly diagnosed PWH in Taiwan during the study period, were tested. Of the PWH with analysable gene sequences, the HIV-1 PDR prevalence to NNRTIs, nucleoside reverse-transcriptase inhibitors (NRTIs), first- and second-generation INSTIs and PIs was 10.0%, 2.1%, 2.5%, 0.6% and 0.4%, respectively. While the trends of PDR remained stable for NRTIs, INSTIs and PIs, there was a significantly increasing trend of PDR to NNRTIs from 6.0% in 2016% to 13.1% in 2022 (P = 0.001).

CONCLUSIONS

After the introduction of second-generation INSTIs in Taiwan, the trends of HIV-1 PDR to NRTIs and INSTIs remained low. Furthermore, there was no significant decrease of the prevalence of PDR toward NNRTIs between 2016 and 2022.

Nationwide Study of Drug Resistance Mutations in HIV-1 Infected Individuals under Antiretroviral Therapy in Brazil

The success of antiretroviral treatment (ART) is threatened by the emergence of drug resistance mutations (DRM). Since Brazil presents the largest number of people living with HIV (PLWH) in South America we aimed at understanding the dynamics of DRM in this country. We analyzed a total of 20,226 HIV-1 sequences collected from PLWH undergoing ART between 2008–2017. Results show a mild decline of DRM over the years but an increase of the K65R reverse transcriptase mutation from 2.23% to 12.11%. This increase gradually occurred following alterations in the ART regimens replacing zidovudine (AZT) with tenofovir (TDF). PLWH harboring the K65R had significantly higher viral loads than those without this mutation (p < 0.001). Among the two most prevalent HIV-1 subtypes (B and C) there was a significant (p < 0.001) association of K65R with subtype C (11.26%) when compared with subtype B (9.27%). Nonetheless, evidence for K65R transmission in Brazil was found both for C and B subtypes. Additionally, artificial neural network-based immunoinformatic predictions suggest that K65R could enhance viral recognition by HLA-B27 that has relatively low prevalence in the Brazilian population. Overall, the results suggest that tenofovir-based regimens need to be carefully monitored particularly in settings with subtype C and specific HLA profiles.

Characterisation of HIV-1 molecular transmission clusters among newly diagnosed individuals infected with non-B subtypes in Italy

OBJECTIVE

We evaluated the characteristics of HIV-1 molecular transmission clusters (MTCs) in 1890 newly diagnosed individuals infected with non-B subtypes between 2005 and 2017 in Italy.

METHODS

Phylogenetic analyses were performed on pol sequences to characterise subtypes/circulating recombinant forms and identify MTCs. MTCs were divided into small (SMTCs, 2-3 sequences), medium (MMTCs, 4-9 sequences) and large (LMTCs, ≥10 sequences). Factors associated with MTCs were evaluated using logistic regression analysis.

RESULTS

145 MTCs were identified and involved 666 individuals (35.2%); 319 of them (16.9%) were included in 13 LMTCs, 111 (5.9%) in 20 MMTCs and 236 (12.5%) in 112 SMTCs. Compared with individuals out of MTCs, individuals involved in MTCs were prevalently Italian (72.7% vs 30.9%, p<0.001), male (82.9% vs 62.3%, p<0.001) and men who have sex with men (MSM) (43.5% vs 14.5%, p<0.001). Individuals in MTCs were also younger (median (IQR) years: 41 (35-49) vs 43 (36-51), p<0.001) and had higher CD4 cell count in comparison with individuals out of MTCs (median (IQR): 10⁹/L: 0.4 (0.265-0.587) vs 0.246 (0.082-0.417), p<0.001). The viral load remained stable between the two groups (median (IQR) log₁₀ copies/mL: 4.8 (4.2-5.5) vs 5.0 (4.3-5.5), p=0.87). Logistic regression confirmed that certain factors such as being MSM, of Italian origin, younger age and higher CD4 cell count were significantly associated with MTCs.

CONCLUSIONS

Our findings show that HIV-1 newly diagnosed individuals infected with non-B subtypes are involved in several MTCs in Italy. These MTCs include mainly Italians and MSM and highlight the complex phenomenon characterising the HIV-1 spread. This is important especially in view of monitoring the HIV epidemic and guiding the public health response.

Transmission of HIV-1 drug resistance mutations within partner-pairs: A cross-sectional study of a primary HIV infection cohort

BACKGROUND

Transmission of human immunodeficiency virus type 1 (HIV-1) drug resistance mutations, particularly that of minority drug-resistant variants, remains poorly understood. Population-based studies suggest that drug-resistant HIV-1 is less transmissible than drug-susceptible viruses. We compared HIV-1 drug-resistant genotypes among partner-pairs in order to assess the likelihood of transmission of drug resistance mutations and investigate the role of minority variants in HIV transmission.

METHODS AND FINDINGS

From 1992-2010, 340 persons with primary HIV-1 infection and their partners were enrolled into observational research studies at the University of Washington Primary Infection Clinic (UWPIC). Out of 50 partner-pairs enrolled, 36 (72%) transmission relationships were confirmed by phylogenetic distance analysis of HIV-1 envelope (env) sequences, and 31 partner-pairs enrolled after 1995 met criteria for this study. Drug resistance mutations in the region of the HIV-1 polymerase gene (pol) that encodes protease and reverse transcriptase were assessed by 454-pyrosequencing. In 25 partner-pairs where the transmission direction could be determined, 12 (48%) transmitters had 1-4 drug resistance mutations (23 total) detected in their HIV-1 populations at a median frequency of 6.0% (IQR 1.5%-98.7%, range 1.0%-99.6%). Of 10 major mutations detected in five transmitters at a frequency >95%, 100% (95% CI 69.2%-100%) were detected in recipients. All of these transmitters were antiretroviral (ARV)-naïve at the time of specimen collection. Fourteen mutations (eight major mutations and six accessory mutations) were detected in nine transmitters at low frequencies (1.0%-11.8%); four of these transmitters had previously received ARV therapy. Two (14% [95% CI 1.8%-42.8%]) G73S accessory mutations were detected in both transmitter and recipient. This number is not significantly different from the number expected based on the observed frequencies of drug-resistant viruses in transmitting partners. Limitations of this study include the small sample size and uncertainties in determining the timing of virus transmission and mutation history.

CONCLUSIONS

Drug-resistant majority variants appeared to be commonly transmitted by ARV-naïve participants in our analysis and may contribute significantly to transmitted drug resistance on a population level. When present at low frequency, no major mutation was observed to be shared between partner-pairs; identification of accessory mutations shared within a pair could be due to transmission, laboratory artifact, or apolipoprotein B mRNA-editing enzyme, catalytic polypeptides (APOBECs), and warrants further study.

Quantifying the fitness cost of HIV-1 drug resistance mutations through phylodynamics

Drug resistant HIV is a major threat to the long-term efficacy of antiretroviral treatment. Around 10% of ART-naïve patients in Europe are infected with drug-resistant HIV type 1. Hence it is important to understand the dynamics of transmitted drug resistance evolution. Thanks to routinely performed drug resistance tests, HIV sequence data is increasingly available and can be used to reconstruct the phylogenetic relationship among viral lineages. In this study we employ a phylodynamic approach to quantify the fitness costs of major resistance mutations in the Swiss HIV cohort. The viral phylogeny reflects the transmission tree, which we model using stochastic birth–death-sampling processes with two types: hosts infected by a sensitive or resistant strain. This allows quantification of fitness cost as the ratio between transmission rates of hosts infected by drug resistant strains and transmission rates of hosts infected by drug sensitive strains. The resistance mutations 41L, 67N, 70R, 184V, 210W, 215D, 215S and 219Q (nRTI-related) and 103N, 108I, 138A, 181C, 190A (NNRTI-related) in the reverse trancriptase and the 90M mutation in the protease gene are included in this study. Among the considered resistance mutations, only the 90M mutation in the protease gene was found to have significantly higher fitness than the drug sensitive strains. The following mutations associated with resistance to reverse transcriptase inhibitors were found to be less fit than the sensitive strains: 67N, 70R, 184V, 219Q. The highest posterior density intervals of the transmission ratios for the remaining resistance mutations included in this study all included 1, suggesting that these mutations do not have a significant effect on viral transmissibility within the Swiss HIV cohort. These patterns are consistent with alternative measures of the fitness cost of resistance mutations. Overall, we have developed and validated a novel phylodynamic approach to estimate the transmission fitness cost of drug resistance mutations.

The introduction of antiretroviral therapy (ART) has decreased mortality and morbidity rates among HIV-infected people, and improved their quality of life. In fact, the WHO states that antiretroviral therapy programmes averted an estimated 7.8 million deaths worldwide between 2000 and 2014. However, the antiretroviral regimen prescribed to a patient may be unable to control HIV infection. Factors that can contribute to treatment failure include drug resistance, drug toxicity, or poor treatment adherence. In this study we aim to understand the dynamics of transmitted drug resistance by analysing the viral sequence data that was collected for resistance testing. We present a novel approach to quantify how drug resistance impacts virus lineage transmissibility, how fast resistance mutations evolve in sensitive strains and how fast they revert back to the sensitive type. We apply our approach to the Swiss HIV cohort study, and obtain patterns of viral transmission fitness that are consistent with alternative, harder to obtain measures of fitness.

Transmission fitness of drug-resistant HIV revealed in a surveillance system transmission network

Test-and-treat programs are central to the global control of HIV, but transmitted drug resistance threatens the effectiveness of these programs. HIV mutations conferring resistance to antiretroviral drugs reduce replicative fitness in vitro, but their effect on propagation in vivo is less understood. Here, we estimate transmission fitness of these mutations in antiretroviral-naïve populations in the U.S. National HIV Surveillance System by comparing their frequency of clustering in a genetic transmission network relative with wild-type viruses. The large dataset (66,221 persons), comprising 30,196 antiretroviral-naïve persons, permitted the evaluation of sixty-nine resistance mutations. Decreased transmission fitness was demonstrated for twenty-three mutations, including M184V. In contrast, many high prevalence mutations (e.g. K103N, Y181C, and L90M) had transmission fitness that was indistinguishable from or exceeded wild-type fitness, permitting the establishment of large, self-sustaining drug resistance reservoirs. We highlight implications of these findings on strategies to preserve global treatment effectiveness.

A phylotype-based analysis highlights the role of drug-naive HIV-positive individuals in the transmission of antiretroviral resistance in the UK

OBJECTIVE

Antiretroviral-naive HIV-positive individuals contribute to the transmission of drug-resistant viruses, compromising first-line therapy. Using phylogenetic inference, we quantified the proportion of transmitted drug-resistance originating from a treatment-naive source.

METHODS

Using a novel phylotype-based approach, 24 550 HIV-1 subtype B partial pol gene sequences from the UK HIV Drug Resistance database were analysed. Ongoing transmission of drug resistance amongst HIV-positive individuals was identified as phylotypes of at least three sequences with at least one shared drug resistance mutation, a maximum intra-clade genetic distance of 4.0% and a basal branch support at least 90%. The time of persistence of the transmission chains was estimated using a fast least-squares molecular clock inference approach.

RESULTS

Around 70% of transmitted drug-resistance had a treatment-naive source. The most commonly transmitted mutations were L90M in the protease gene and K103N, T215D and T215S in reverse transcriptase. Reversion to wild type occurred at a low frequency and drug-independent reservoirs of resistance have persisted for up to 13 years.

CONCLUSION

These results illustrate the impact of viral fitness on the establishment of resistance reservoirs and support the notion that earlier diagnoses and treatment of HIV infections are warranted for counteracting the spread of antiretroviral resistance. Phylotype-based phylogenetic inference is an attractive approach for the routine surveillance of transmitted drug resistance in HIV as well as in other pathogens for which genotypic resistance data are available.

Persistence of HIV-1 transmitted drug resistance mutations

There are few data on the persistence of individual human immunodeficiency virus type 1 (HIV-1) transmitted drug resistance (TDR) mutations in the absence of selective drug pressure. We studied 313 patients in whom TDR mutations were detected at their first resistance test and who had a subsequent test performed while ART-naive. The rate at which mutations became undetectable was estimated using exponential regression accounting for interval censoring. Most thymidine analogue mutations (TAMs) and T215 revertants (but not T215F/Y) were found to be highly stable, with NNRTI and PI mutations being relatively less persistent. Our estimates are important for informing HIV transmission models.

Detectable signals of episodic risk effects on acute HIV transmission: strategies for analyzing transmission systems using genetic data

Episodic high-risk sexual behavior is common and can have a profound effect on HIV transmission. In a model of HIV transmission among men who have sex with men (MSM), changing the frequency, duration and contact rates of high-risk episodes can take endemic prevalence from zero to 50% and more than double transmissions during acute HIV infection (AHI). Undirected test and treat could be inefficient in the presence of strong episodic risk effects. Partner services approaches that use a variety of control options will be likely to have better effects under these conditions, but the question remains: What data will reveal if a population is experiencing episodic risk effects? HIV sequence data from Montreal reveals genetic clusters whose size distribution stabilizes over time and reflects the size distribution of acute infection outbreaks (AIOs). Surveillance provides complementary behavioral data. In order to use both types of data efficiently, it is essential to examine aspects of models that affect both the episodic risk effects and the shape of transmission trees. As a demonstration, we use a deterministic compartmental model of episodic risk to explore the determinants of the fraction of transmissions during acute HIV infection (AHI) at the endemic equilibrium. We use a corresponding individual-based model to observe AIO size distributions and patterns of transmission within AIO. Episodic risk parameters determining whether AHI transmission trees had longer chains, more clustered transmissions from single individuals, or different mixes of these were explored. Encouragingly for parameter estimation, AIO size distributions reflected the frequency of transmissions from acute infection across divergent parameter sets. Our results show that episodic risk dynamics influence both the size and duration of acute infection outbreaks, thus providing a possible link between genetic cluster size distributions and episodic risk dynamics.

Inefficient vaginal transmission of tenofovir-resistant HIV-1

Transmission of drug-resistant HIV has been postulated to be a threat to current first-line antiretroviral therapy (ART) regimens and the efficacy of several antiretroviral-based preexposure prophylaxis (PrEP) strategies being tested. Here we evaluated the effect of the common tenofovir (TFV) resistance mutation K65R on vaginal HIV transmission. Our results demonstrate that despite no overt loss of overall replication competence in vivo, this mutation results in significantly reduced mucosal transmission. When transmitted, the mutant virus eventually reverted to the wild type in 2 of 3 animals examined.

Increase in transmitted HIV drug resistance among persons undergoing genotypic resistance testing in Ontario, Canada, 2002-09

OBJECTIVES

To characterize persons undergoing HIV genotypic resistance testing (GRT) while treatment naive and to estimate the prevalence of transmitted HIV drug resistance (TDR) among HIV-positive outpatients in Ontario, Canada.

METHODS

We analysed data from a multi-site cohort of persons receiving HIV care. Data were obtained from medical chart abstractions, interviews and record linkage with the Public Health Laboratories, Public Health Ontario. The analysis was restricted to 626 treatment-naive persons diagnosed in 2002-09. TDR mutations were identified using the calibrated population resistance tool. We used descriptive statistics and regression methods to characterize treatment-naive GRT test uptake and patterns of TDR.

RESULTS

Overall, 53.2% (333/626) of participants had baseline GRT. The proportion increased with year of HIV diagnosis from 30.0% in 2002 to 82.6% in 2009 (P < 0.0001). Among those tested, 13.6% (CI 9.9-17.3%) had one or more drug resistance mutations, and 8.8% (CI 5.7-11.8%), 4.8% (CI 2.5-7.2%) and 2.7% (CI 1.0-4.5%) had mutations conferring resistance to nucleoside/tide reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs) and protease inhibitors (PIs), respectively. TDR prevalence increased from 2002-07 to 2008-09 (adjusted OR 3.7, 95% CI 1.7-8.2), driven by a higher proportion with NRTI (18.2% versus 5.9%, P = 0.0009) and NNRTI mutations (11.7% versus 2.8%, P = 0.004) in the later time period. PI TDR remained unchanged.

CONCLUSIONS

Baseline GRT increased dramatically since 2002, but remains below 100%. The prevalence of overall TDR tripled due to increases in NRTI and NNRTI mutations. These findings highlight the value of routine baseline GRT for TDR surveillance and patient care.