Evolution of lamivudine resistance in human immunodeficiency virus type 1-infected individuals: the relative roles of drift and selection

Neurosymptomatic HIV-1 CSF escape is associated with replication in CNS T cells and inflammation

During antiretroviral therapy (ART), most people living with HIV-1 have undetectable HIV-1 RNA in their plasma. However, they occasionally present with new or progressive neurologic deficits and detectable HIV-1 RNA in the cerebrospinal fluid (CSF), a condition defined as neurosymptomatic HIV-1 CSF escape (NSE). We explored the source of neuropathogenesis and HIV-1 RNA in the CSF during NSE by characterizing HIV-1 populations and inflammatory biomarkers in CSF from 25 individuals with NSE. HIV-1 populations in the CSF were uniformly drug resistant and adapted to replication in CD4+ T cells, but differed greatly in genetic diversity, with some having low levels of diversity similar to those observed during untreated primary infection and others having high levels like those during untreated chronic infection. Higher diversity in the CSF during NSE was associated with greater CNS inflammation. Finally, optimization of ART regimen was associated with viral suppression and improvement of neurologic symptoms. These results are consistent with CNS inflammation and neurologic injury during NSE being driven by replication of partially drug-resistant virus in CNS CD4+ T cells. This is unlike nonsuppressible viremia in the plasma during ART, which typically lacks clinical consequences and is generated by virus expression without replication.

Oral HIV pre-exposure prophylaxis use and resistance-associated mutations among men who have sex with men and transgender persons newly diagnosed with HIV in the Netherlands: results from the ATHENA cohort, 2018 to 2022

BackgroundIn the Netherlands, HIV pre-exposure prophylaxis (PrEP) has been available since 2019. However, the extent of PrEP use prior to HIV diagnosis and development of PrEP-resistance-associated mutations (RAMs) is not known.AimWe assessed prior PrEP use and potential transmission of PrEP RAMs among men who have sex with men (MSM) and transgender persons (TGP) with a new HIV diagnosis in the Netherlands.MethodsData on prior PrEP use between 1 January 2018 and 31 December 2022 were available from the Dutch national ATHENA cohort. We assessed proportion of prior PrEP use, detected PrEP associated RAMs and assessed potential onward transmission of RAMs between 2010 and 2022 using a maximum likelihood tree.ResultsData on prior PrEP use were available for 583/1,552 (36.3%) individuals, with 16% (94/583) reporting prior PrEP use. In 489 individuals reporting no prior PrEP use, 51.5% did not use PrEP due to: low HIV-risk perception (29%), no access (19.1%), personal preference (13.1%), and being unaware of PrEP (19.1%). For PrEP users, 13/94 (13.8%) harboured a M184V/I mutation, of whom two also harboured a K65R mutation. In people with a recent HIV infection, detection of PrEP RAMs increased from 0.23% (2/862) before 2019 to 4.11% (9/219) from 2019. We found no evidence of onward transmission of PrEP RAMs.ConclusionThe prevalence of PrEP-associated RAMs has increased since PrEP became available in the Netherlands. More widespread access to PrEP and retaining people in PrEP programmes when still at substantial risk is crucial to preventing new HIV infections.

QuickFit: A High-Throughput RT-qPCR-Based Assay to Quantify Viral Growth and Fitness In Vitro

Quantifying viral growth rates is key to understanding evolutionary dynamics and the potential for mutants to escape antiviral drugs. Defining evolutionary escape paths and their impact on viral fitness allows for the development of drugs that are resistant to escape. In the case of HIV, combination antiretroviral therapy can successfully prevent or treat infection, but it relies on strict adherence to prevent escape. Here, we present a method termed QuickFit that enables the quantification of viral fitness by employing large numbers of parallel viral cultures to measure growth rates accurately. QuickFit consistently recapitulated HIV growth measurements obtained by traditional approaches, but with significantly higher throughput and lower rates of error. This method represents a promising tool for rapid and consistent evaluation of viral fitness.

Drug resistance and influencing factors in HIV-1-infected individuals under antiretroviral therapy in Guangxi, China

OBJECTIVES

To assess the profiles and determinants of drug resistance in HIV-1-infected individuals undergoing ART in Guangxi.

METHODS

Samples and data were collected from HIV-1-infected individuals experiencing virological failure post-ART from 14 cities in Guangxi. Sequencing of the HIV-1 pol gene was conducted, followed by analysis for drug resistance mutations using the Stanford University HIV Drug Resistance Database. Logistic regression was employed to identify potential risk factors associated with both HIV drug resistance and mortality.

RESULTS

A total of 8963 individuals with pol sequences were included in this study. The overall prevalence of HIV-1 drug resistance (HIVDR) was 42.43% (3808/8963), showing a decrease from 59.62% to 41.40% from 2016 to 2023. Factors such as being aged ≥50 years, male, Han nationality, lower education levels, occupations including workers, peasants and children, AIDS, pre-treatment CD4 T cell counts <200 cells/mm3, infection with CRF01_AE and CRF55_01B subtypes, and ART regimen lamivudine/zidovudine/nevirapine were associated with higher susceptibility to HIVDR. The common mutations were M184V (17.38%) and K103N (22.14%). Additionally, the prevalence of M184V, S68G, M41L and G190A were different between the Han and Zhuang populations. Factors including age, gender, ethnicity, education level, occupation, infectious route, clinical stage, viral load, subtype, ART regimen and HIVDR showed significant associations with mortality.

CONCLUSIONS

The factors contributing to drug resistance in the HIV-1 ART individuals in Guangxi appear to be notably intricate. Continuous reinforcement of drug resistance surveillance is imperative, accompanied by the optimization of ART regimens to mitigate virological failures effectively.

Bictegravir/Emtricitabine/Tenofovir Alafenamide for HIV-1: What is the Hidden Potential of This Emerging Treatment?

Bictegravir/emtricitabine/tenofovir alafenamide (B/F/TAF) is a single-tablet antiretroviral therapy regimen. B/F/TAF has become a popular treatment choice because of its small tablet size, high barrier to resistance, favorable tolerability, and limited drug-drug interaction profile. Continued research on B/F/TAF has revealed additional potential for this regimen. This review presents recent literature supporting the use of B/F/TAF as an option for consolidating therapy and maintaining virologic suppression in individuals despite M184V/I mutations. Additionally, children are a unique patient population with limited antiviral options. Standard dose B/F/TAF has demonstrated similar drug exposure in children and adolescents as adults, and low-dose B/F/TAF is approved for children living with HIV greater than two years of age and weighing at least 14 kg. Data supporting this recommendation is described in this review. Finally, despite a lack of prospective data, B/F/TAF may have a role in the future of pre- and post-exposure prophylaxis. This review discusses these discoveries and the continued exploration of the hidden potential of B/F/TAF.

High-level dolutegravir resistance can emerge rapidly from few variants and spread by recombination: implications for integrase strand transfer inhibitor salvage therapy

The integrase strand transfer inhibitor (INSTI) dolutegravir is commonly used in combination antiretroviral therapy regimens and retains strong potency even with primary resistance mutations to some other INSTIs. Acquisition of accessory mutations to primary mutations results in significant increases in dolutegravir resistance. Previously, we reported that addition of the secondary mutation T97A can result in rapid treatment failure in individuals with INSTI mutations at positions 140 and 148. Here, we conducted a detailed case study of one of these individuals and find that T97A-containing HIV emerged from a large replicating population from only a few (≤4) viral lineages. When combined with primary INSTI resistance mutations, T97A provides a strong selective advantage; the finding that T97A-containing variants spread by replication and recombination, and persisted for months after discontinuing dolutegravir, has important implications as dolutegravir is rolled out worldwide.

High efficacy of switching to bictegravir/emtricitabine/tenofovir alafenamide in people with suppressed HIV and preexisting M184V/I

OBJECTIVE

We investigated the prevalence of preexisting M184V/I and associated risk factors among clinical trial participants with suppressed HIV and evaluated the impact of M184V/I on virologic response after switching to bictegravir/emtricitabine/tenofovir alafenamide (B/F/TAF).

DESIGN

Participant data were pooled from six clinical trials investigating the safety and efficacy of switching to B/F/TAF in virologically suppressed people with HIV.

METHODS

Preexisting drug resistance was assessed by historical genotypes and/or baseline proviral DNA genotyping. Virologic outcomes were determined by last available on-treatment HIV-1 RNA. Stepwise selection identified potential risk factors for M184V/I in a multivariate logistic regression model.

RESULTS

Altogether, 2034 participants switched treatment regimens to B/F/TAF and had follow-up HIV-1 RNA data, and 1825 of these participants had baseline genotypic data available. Preexisting M184V/I was identified in 182 (10%), mostly by baseline proviral DNA genotype ( n  = 167). Most substitutions were M184V ( n  = 161) or M184V/I mixtures ( n  = 10). Other resistance substitutions were often detected in addition to M184V/I ( n  = 147). At last on-treatment visit, 98% (179/182) with preexisting M184V/I and 99% (2012/2034) of all B/F/TAF-treated participants had HIV-1 RNA less than 50 copies/ml, with no treatment-emergent resistance to B/F/TAF. Among adult participants, factors associated with preexisting M184V/I included other resistance, black race, Hispanic/Latinx ethnicity, lower baseline CD4 + cell count, advanced HIV disease, longer duration of antiretroviral therapy, and greater number of prior third agents.

CONCLUSION

M184V/I was detected in 10% of virologically suppressed clinical trial participants at study baseline. Switching to B/F/TAF demonstrated durable efficacy in maintaining viral suppression, including in those with preexisting M184V/I.

Understanding drug resistance patterns across different classes of antiretrovirals used in HIV-1-infected treatment-Naïve and experienced patients in Mumbai, India

Background

The aim of this study is to find out the proportion of treatment-naïve (Tn) and treatment-experienced (Te) patients experiencing HIV drug resistance (DR) to different classes of antiretrovirals (ARVs) being used for HIV treatment and their in class DR correlation.

Methods

A cross-sectional study was done on 109 HIV patients enrolled at a private hospital in Thane, India, from 2014 to 2019. All patients were tested for CD4 count, viral load, and resistance to ARVs.

Results

Sixty-six patients were Tn and 43 patients were Te. Among Tn and Te patients, the percentage of high-level resistance (HLR) for nonnucleoside reverse transcriptase inhibitors (NNRTI) was 4.55% and 37.8%, respectively, for nucleoside reverse transcriptase inhibitors (NRTI) was 0.43% and 36.4%, respectively. No HLR was observed for protease inhibitors (PIs) among Tn patients, while Te patients showed 2.62% HLR. Tn and Te patients showed high susceptibility for Darunavir (98.48% and 95.34%, respectively) followed by Atazanavir and Lopinavir (96.96%, each and 90.69%, each). Tn patients showed HLR for Lamivudine and Emtricitabine (1.52%, each). Integrase Strand Transfer Inhibitors were susceptible (100%) in both Tn and Te patients. A positive correlation was observed for within class across ARVs.

Conclusion

An increased incidence of HLR was observed for NNRTI as compared to NRTI while PIs and integrase strand transfer inhibitors (INSTIs) demonstrated no HLR in either group of patients. When selecting a regimen for Tn patients consisting of NRTIs + NNRTIs genotypic DR test is essential. While with PIs or INSTIs its optional. Among Te patients, DR testing is recommended for all classes of drugs.

Characterizing HIV-1 Genetic Subtypes and Drug Resistance Mutations among Children, Adolescents and Pregnant Women in Sierra Leone

Human immunodeficiency virus (HIV) drug resistance (HIVDR) is widespread in sub-Saharan Africa. Children and pregnant women are particularly vulnerable, and laboratory testing capacity remains limited. We, therefore, used a cross-sectional design and convenience sampling to characterize HIV subtypes and resistance-associated mutations (RAMs) in these groups in Sierra Leone. In total, 96 children (age 2-9 years, 100% ART-experienced), 47 adolescents (age 10-18 years, 100% ART-experienced), and 54 pregnant women (>18 years, 72% ART-experienced) were enrolled. Median treatment durations were 36, 84, and 3 months, respectively, while the sequencing success rates were 45%, 70%, and 59%, respectively, among children, adolescents, and pregnant women. Overall, the predominant HIV-1 subtype was CRF02_AG (87.9%, 95/108), with minority variants constituting 12%. Among children and adolescents, the most common RAMs were M184V (76.6%, n = 49/64), K103N (45.3%, n = 29/64), Y181C/V/I (28.1%, n = 18/64), T215F/Y (25.0%, n = 16/64), and V108I (18.8%, n = 12/64). Among pregnant women, the most frequent RAMs were K103N (20.6%, n = 7/34), M184V (11.8%, n = 4/34), Y181C/V/I (5.9%, n = 2/34), P225H (8.8%, n = 3/34), and K219N/E/Q/R (5.9%, n = 2/34). Protease and integrase inhibitor-RAMs were relatively few or absent. Based on the genotype susceptibility score distributions, 73%, 88%, and 14% of children, adolescents, and pregnant women, respectively, were not susceptible to all three drug components of the WHO preferred first-line regimens per 2018 guidelines. These findings suggest that routine HIVDR surveillance and access to better ART choices may improve treatment outcomes in Sierra Leone.

Characteristics of drug resistance in HIV-1 CRF55_01B from ART-experienced patients in Guangdong, China

BACKGROUND

HIV-1 acquired drug resistance (ADR) has become a critical clinical and public health issue. Recently, HIV-1 CRF55_01B has been found more frequently in the MSM population.

OBJECTIVE

To investigate the characteristics of HIV-1 drug resistance mutations (DRMs) and the extent of changes in drug susceptibility among ART-experienced CRF55_01B-infected adults of Guangdong.

METHODS

ADR was tested for immediately in CRF55_01B-infected patients with virological failure. Demographic and epidemiological information was collected. DRMs and antiretroviral susceptibility were interpreted using the Stanford University HIV Drug Resistance Database HIVdb program.

RESULTS

Overall, 162 (4.78%) CRF55_01B isolates were identified from 2013 to 2018. Among DRMs, M184V (43.83%) was the most frequent NRTI DRM, followed by K65R (23.46%), and V179E (98.77%) was the most frequent NNRTI DRM, followed by K103N (47.53%) and Y181C (14.81%). According to the HIVdb program, 79.01% of the CRF55_01B-infected patients carried mutations conferring low-level or higher drug resistance to any of the three classes of ART drugs. Among PI DRMs, only one mutation affording low-level resistance to nelfinavir was found (0.62%). Among NRTI DRMs, a high proportion of high-level resistance to lamivudine (58.64%) and emtricitabine (58.02%) was found. As regards NNRTIs, more than 75% of patients carried efavirenz and nevirapine DRMs. The percentages of high-level resistance were 70.99%, 63.58%, 22.22%, 17.90% and 4.32% for nevirapine, efavirenz, rilpivirine, doravirine and etravirine, respectively.

CONCLUSIONS

High frequencies of DRMs and resistance were observed among CRF55_01B-infected patients failing ART in Guangdong, and interventions may be considered to minimize ecological contributions to ART.

Acute HIV at the Time of Initiation of Pre-exposure or Post-exposure Prophylaxis: Impact on Drug Resistance and Clinical Outcomes

BACKGROUND

Initiating pre-exposure or post-exposure prophylaxis (PrEP/PEP) in the setting of undiagnosed acute HIV (AHI) could cause antiretroviral resistance. We sought to characterize clinical outcomes and drug resistance mutations among individuals prescribed PrEP/PEP with undiagnosed AHI at a San Francisco sexually transmitted disease clinic.

SETTING

In our PrEP/PEP program, patients are tested for HIV using a point-of-care antibody test. If negative, patients are started on prophylaxis and screened for AHI using pooled HIV RNA (5-10 days turn-around). We used 2-drug PEP until 05/2016.

METHODS

We identified patients who had as-yet-undiagnosed AHI on the day of PrEP/PEP start between 2011 and 2018, then used our clinical record and surveillance data to describe HIV resistance and clinical outcomes.

RESULTS

Of 1758 PrEP and 2242 PEP starts, there were 7 AHI cases among PrEP users (0.40%) and 6 among PEP users (0.30%). Median times for linkage to HIV care, initiation of HIV treatment, and viral suppression were 7, 12, and 43 days. On initiation of HIV care, 3 patients (23%) were found to have an M184 mutation 7-12 days after starting PrEP/PEP. All 3 had genotyping performed on stored serum available from the date of PrEP/PEP start, each of which demonstrated wild-type virus. All 3 patients achieved durable viral suppression.

CONCLUSIONS

Although rare (occurring <0.5% of the time), AHI in the setting of PrEP/2-drug PEP can result in an M184 within days. Even with M184, persons with AHI achieve viral suppression when rapidly linked to care and initiated on antiretroviral therapy. Providers should consider AHI screening when starting PrEP/PEP.

Dolutegravir response in antiretroviral therapy naïve and experienced patients with M184V/I: Impact in low-and middle-income settings

BACKGROUND

Dolutegravir (DTG) is now recommended to all HIV infected adults, adolescents, and children of right age by WHO. The low cost of $75 per year for generic DTG-based combination, has allowed 3.9 million people living with HIV (PLWH) in low and middle-income countries (LMICs) access to DTG. Lamivudine and emtricitabine associated M184V/I mutation is highly prevalent in PLWH and the majority of HIV infected individuals receiving DTG regimens may already be carrying M184V/I mutation.

DISCUSSION

Despite high prevalence of M184V/I in antiretroviral therapy (ART) experienced patients, DTG treatment outcomes will likely not be adversely affected by this mutation. The use of DTG in ART naïve has been largely characterised by rare emergence of resistance and virological failure. DTG-based regimens have to great extent been effective at maintaining viral suppression in treatment experienced PLWH carrying M184V/I.

CONCLUSIONS

Initiating patients on DTG may help preserve more treatment options for HIV infected individuals living in LMICs. High genetic barrier to the development of resistance associated with DTG and progressive viral suppression in patients switched to DTG-based therapy with M184V/I, may encourage better DTG outcomes and help in curbing increasing levels of HIV drug resistance in LMICs.

The clarifying role of time series data in the population genetics of HIV

HIV can evolve remarkably quickly in response to antiretroviral therapies and the immune system. This evolution stymies treatment effectiveness and prevents the development of an HIV vaccine. Consequently, there has been a great interest in using population genetics to disentangle the forces that govern the HIV adaptive landscape (selection, drift, mutation, and recombination). Traditional population genetics approaches look at the current state of genetic variation and infer the processes that can generate it. However, because HIV evolves rapidly, we can also sample populations repeatedly over time and watch evolution in action. In this paper, we demonstrate how time series data can bound evolutionary parameters in a way that complements and informs traditional population genetic approaches. Specifically, we focus on our recent paper (Feder et al., 2016, eLife), in which we show that, as improved HIV drugs have led to fewer patients failing therapy due to resistance evolution, less genetic diversity has been maintained following the fixation of drug resistance mutations. Because soft sweeps of multiple drug resistance mutations spreading simultaneously have been previously documented in response to the less effective HIV therapies used early in the epidemic, we interpret the maintenance of post-sweep diversity in response to poor therapies as further evidence of soft sweeps and therefore a high population mutation rate (θ) in these intra-patient HIV populations. Because improved drugs resulted in rarer resistance evolution accompanied by lower post-sweep diversity, we suggest that both observations can be explained by decreased population mutation rates and a resultant transition to hard selective sweeps. A recent paper (Harris et al., 2018, PLOS Genetics) proposed an alternative interpretation: Diversity maintenance following drug resistance evolution in response to poor therapies may have been driven by recombination during slow, hard selective sweeps of single mutations. Then, if better drugs have led to faster hard selective sweeps of resistance, recombination will have less time to rescue diversity during the sweep, recapitulating the decrease in post-sweep diversity as drugs have improved. In this paper, we use time series data to show that drug resistance evolution during ineffective treatment is very fast, providing new evidence that soft sweeps drove early HIV treatment failure.

Predominance of positive epistasis among drug resistance-associated mutations in HIV-1 protease

Drug-resistant mutations often have deleterious impacts on replication fitness, posing a fitness cost that can only be overcome by compensatory mutations. However, the role of fitness cost in the evolution of drug resistance has often been overlooked in clinical studies or in vitro selection experiments, as these observations only capture the outcome of drug selection. In this study, we systematically profile the fitness landscape of resistance-associated sites in HIV-1 protease using deep mutational scanning. We construct a mutant library covering combinations of mutations at 11 sites in HIV-1 protease, all of which are associated with resistance to protease inhibitors in clinic. Using deep sequencing, we quantify the fitness of thousands of HIV-1 protease mutants after multiple cycles of replication in human T cells. Although the majority of resistance-associated mutations have deleterious effects on viral replication, we find that epistasis among resistance-associated mutations is predominantly positive. Furthermore, our fitness data are consistent with genetic interactions inferred directly from HIV sequence data of patients. Fitness valleys formed by strong positive epistasis reduce the likelihood of reversal of drug resistance mutations. Overall, our results support the view that strong compensatory effects are involved in the emergence of clinically observed resistance mutations and provide insights to understanding fitness barriers in the evolution and reversion of drug resistance.

Impact of Suboptimal APOBEC3G Neutralization on the Emergence of HIV Drug Resistance in Humanized Mice

HIV diversification facilitates immune escape and complicates antiretroviral therapy. In this study, we take advantage of a humanized-mouse model to probe the contribution of APOBEC3 mutagenesis to viral evolution. Humanized mice were infected with isogenic HIV molecular clones (HIV-WT, HIV-45G, and HIV-ΔSLQ) that differ in their abilities to counteract APOBEC3G (A3G). Infected mice remained naive or were treated with the reverse transcriptase (RT) inhibitor lamivudine (3TC). Viremia, emergence of drug-resistant variants, and quasispecies diversification in the plasma compartment were determined throughout infection. While both HIV-WT and HIV-45G achieved robust infection, over time, HIV-45G replication was significantly reduced compared to that of HIV-WT in the absence of 3TC treatment. In contrast, treatment responses differed significantly between HIV-45G- and HIV-WT-infected mice. Antiretroviral treatment failed in 91% of HIV-45G-infected mice, while only 36% of HIV-WT-infected mice displayed a similar negative outcome. Emergence of 3TC-resistant variants and nucleotide diversity were determined by analyzing 155,462 single HIV reverse transcriptase gene (RT) and 6,985 vif sequences from 33 mice. Prior to treatment, variants with genotypic 3TC resistance (RT-M184I/V) were detected at low levels in over a third of all the animals. Upon treatment, the composition of the plasma quasispecies rapidly changed, leading to a majority of circulating viral variants encoding RT-184I. Interestingly, increased viral diversity prior to treatment initiation correlated with higher plasma viremia in HIV-45G-infected animals, but not in HIV-WT-infected animals. Taken together, HIV variants with suboptimal anti-A3G activity were attenuated in the absence of selection but displayed a fitness advantage in the presence of antiretroviral treatment.IMPORTANCE Both viral (e.g., RT) and host (e.g., A3G) factors can contribute to HIV sequence diversity. This study shows that suboptimal anti-A3G activity shapes viral fitness and drives viral evolution in the plasma compartment in humanized mice.

Elucidating molecular interactions of L-nucleotides with HIV-1 reverse transcriptase and mechanism of M184V-caused drug resistance

Emtricitabine (FTC) and lamivudine (3TC), containing an oxathiolane ring with unnatural (-)-stereochemistry, are widely used nucleoside reverse transcriptase inhibitors (NRTIs) in anti-HIV therapy. Treatment with FTC or 3TC primarily selects for the HIV-1 RT M184V/I resistance mutations. Here we provide a comprehensive kinetic and structural basis for inhibiting HIV-1 RT by (-)-FTC-TP and (-)-3TC-TP and drug resistance by M184V. (-)-FTC-TP and (-)-3TC-TP have higher binding affinities (1/Kd) for wild-type RT but slower incorporation rates than dCTP. HIV-1 RT ternary crystal structures with (-)-FTC-TP and (-)-3TC-TP corroborate kinetic results demonstrating that their oxathiolane sulfur orients toward the DNA primer 3'-terminus and their triphosphate exists in two different binding conformations. M184V RT displays greater (>200-fold) Kd for the L-nucleotides and moderately higher (>9-fold) Kd for the D-isomers compared to dCTP. The M184V RT structure illustrates how the mutation repositions the oxathiolane of (-)-FTC-TP and shifts its triphosphate into a non-productive conformation.

A 28-Year History of HIV-1 Drug Resistance and Transmission in Washington, DC

Washington, DC consistently has one of the highest annual rates of new HIV-1 diagnoses in the United States over the last 10 years. To guide intervention and prevention strategies to combat DC HIV infection, it is helpful to understand HIV transmission dynamics in a historical context. Toward this aim, we conducted a retrospective study (years 1987-2015) of 3,349 HIV pol sequences (1,026 bp) from 1,996 individuals living in the DC area belonging to three different cohorts. We coupled HIV sequence data with clinical information (sex, risk factor, race/ethnicity, viral load, subtype, anti-retroviral regimen) to identify circulating drug resistant mutations (DRM) and transmission clusters and assess their persistence over time. Of the transmission clusters identified in the DC area, 78.0 and 31.7% involved MSM and heterosexuals, respectively. The longest spread of time for a single cluster was 5 years (2007-2012) using a distance-based network inference approach and 27 years (1987-2014) using a maximum likelihood phylogenetic approach. We found eight subtypes and nine recombinants. Genetic diversity increased steadily over time with a slight peak in 2009 and remained constant thereafter until 2015. Nucleotide diversity also increased over time while relative genetic diversity (BEAST) remained relatively steady over the last 28 years with slight increases since 2000 in subtypes B and C. Sequences from individuals on drug therapy contained the highest total number of DRMs (1,104-1,600) and unique DRMs (63-97) and the highest proportion (>20%) of resistant individuals. Heterosexuals (43.94%), MSM (40.13%), and unknown (44.26%) risk factors showed similar prevalence of DRMs, while injection drug users had a lower prevalence (33.33%). Finally, there was a 60% spike in the number of codons with DRMs between 2007 and 2010. Past patterns of HIV transmission and DRM accumulation over time described here will help to predict future efficacy of ART drugs based on DRMs persisting over time and identify risk groups of interest for prevention and intervention efforts within the DC population. Our results show how longitudinal data can help to understand the temporal dynamics of HIV-1 at the local level.

New resistance mutations to nucleoside reverse transcriptase inhibitors at codon 184 of HIV-1 reverse transcriptase (M184L and M184T)

Mutations at HIV-1 reverse transcriptase (RT) codon 184 such as M184V confer resistance to two nucleos(t)ide RT inhibitors (NRTI), lamivudine (3TC) and emtricitabine (FTC). The prevalence of mutations at HIV-1 RT codon 184 was evaluated using three independent RT sequence databases from treatment-experienced (TE) and treatment-naïve (TN) individuals. Data were collected retrospectively from three centers: one in Italy and two in France between 1997 and 2016. In order to highlight the role of these mutations in conferring drug resistance, structural and thermodynamic analyses were conducted by means of computational approaches. Among 32,440 RT sequences isolated from TE and 12,365 isolated from TN patients, the prevalence of HIV-1 RT codon 184 substitutions in each group was 31.21% and 0.72%, respectively. The mutations M184L and M184T have been observed only in TE patients. In all cases but four, M184L and M184T mutations were present during NRTI treatment. Molecular recognition studies on M184L and M184T structures showed both FTC and 3TC thermodynamic profiles unfavorable in comparison with the wild-type sequence, corroborated by molecular dynamic simulations (MDS). In this study, we highlighted two new resistance mutations in vivo for NRTI resistance. The low frequency of this pathway can be related to high impairment of replicative capacity mediated by these mutations.

Sanger and Next Generation Sequencing Approaches to Evaluate HIV-1 Virus in Blood Compartments

The implementation of antiretroviral treatment combined with the monitoring of drug resistance mutations improves the quality of life of HIV-1 positive patients. The drug resistance mutation patterns and viral genotypes are currently analyzed by DNA sequencing of the virus in the plasma of patients. However, the virus compartmentalizes, and different T cell subsets may harbor distinct viral subsets. In this study, we compared the patterns of HIV distribution in cell-free (blood plasma) and cell-associated viruses (peripheral blood mononuclear cells, PBMCs) derived from ART-treated patients by using Sanger sequencing- and Next-Generation sequencing-based HIV assay. CD4⁺CD45RA⁻RO⁺ memory T-cells were isolated from PBMCs using a BD FACSAria instrument. HIV pol (protease and reverse transcriptase) was RT-PCR or PCR amplified from the plasma and the T-cell subset, respectively. Sequences were obtained using Sanger sequencing and Next-Generation Sequencing (NGS). Sanger sequences were aligned and edited using RECall software (beta v3.03). The Stanford HIV database was used to evaluate drug resistance mutations. Illumina MiSeq platform and HyDRA Web were used to generate and analyze NGS data, respectively. Our results show a high correlation between Sanger sequencing and NGS results. However, some major and minor drug resistance mutations were only observed by NGS, albeit at different frequencies. Analysis of low-frequency drugs resistance mutations and virus distribution in the blood compartments may provide information to allow a more sustainable response to therapy and better disease management.

Tackling HIV and AIDS: contributions by non-human primate models

During the past three decades, non-human primate (NHP) models have gained an increasing importance in HIV basic and translational research. In contrast to natural host models, infection of macaques with virulent simian or simian-human immunodeficiency viruses (SIV, SHIV) results in a disease that closely resembles HIV infection and AIDS. Although there is no perfect animal model, and each of the available models has its benefits and limitations, carefully designed NHP studies with selection of experimental variables have unraveled important questions of basic pathogenesis and have provided the tools to explore and screen intervention strategies. For example, NHP studies have advanced our understanding of the crucial events during early infection, and have provided proof-of-concept of antiretroviral drug treatment and prevention strategies such as pre-exposure prophylaxis (PrEP) regimes that are increasingly used worldwide, and upon overcoming further barriers of implementation, have the potential to make the next generation AIDS-free. Remaining goals include the pursuit of an effective HIV vaccine, and HIV cure strategies that would allow HIV-infected people to ultimately stop taking antiretroviral drugs. Through a reiterative process with feed-back from results of human studies, NHP models can be further validated and strengthened to advance our scientific knowledge and guide clinical trials.

Molecular Epidemiology of HIV-1 Virus in Puerto Rico: Novel Cases of HIV-1 Subtype C, D, and CRF-24BG

HIV-1 subtype B virus is the most prevalent subtype in Puerto Rico (PR), accounting for about 90% of infection in the island. Recently, other subtypes and circulating recombinant forms (CRFs), including F(12_BF), A (01_BF), and CRF-39 BF-like, have been identified. The purpose of this study is to assess the distribution of drug resistance mutations and subtypes in PR. A total of 846 nucleotide sequences from the period comprising 2013 through 2017 were obtained from our "HIV Genotyping" test file. Phylogenetic and molecular epidemiology analyses were performed to evaluate the evolutionary dynamics and prevalence of drug resistance mutations. According to our results, we detected a decrease in the prevalence of protease inhibitor, nucleoside reverse transcriptase inhibitor (NRTI), and non-NRTI (NNRTI) resistance mutations over time. In addition, we also detected recombinant forms and, for the first time, identified subtypes C, D, and CRF-24BG in PR. Recent studies suggest that non-subtypes B are associated with a high risk of treatment failure and disease progression. The constant monitoring of viral evolution and drug resistance mutation dynamics is important to establish appropriate efforts for controlling viral expansion.

Neutral Theory and Rapidly Evolving Viral Pathogens

The evolution of viral pathogens is shaped by strong selective forces that are exerted during jumps to new hosts, confrontations with host immune responses and antiviral drugs, and numerous other processes. However, while undeniably strong and frequent, adaptive evolution is largely confined to small parts of information-packed viral genomes, and the majority of observed variation is effectively neutral. The predictions and implications of the neutral theory have proven immensely useful in this context, with applications spanning understanding within-host population structure, tracing the origins and spread of viral pathogens, predicting evolutionary dynamics, and modeling the emergence of drug resistance. We highlight the multiple ways in which the neutral theory has had an impact, which has been accelerated in the age of high-throughput, high-resolution genomics.

Detection of minority drug resistant mutations in Malawian HIV-1 subtype C-positive patients initiating and on first-line antiretroviral therapy

Background

Minority drug resistance mutations (DRMs) that are often missed by Sanger sequencing are clinically significant, as they can cause virologic failure in individuals treated with antiretroviral therapy (ART) drugs.

Objective

This study aimed to estimate the prevalence of minor DRMs among patients enrolled in a Malawi HIV drug resistance monitoring survey at baseline and at one year after initiation of ART.

Methods

Forty-one plasma specimens collected from HIV-1 subtype C-positive patients and seven clonal control samples were analysed using ultra-deep sequencing technology.

Results

Deep sequencing identified all 72 DRMs detected by Sanger sequencing at the level of ≥20% and 79 additional minority DRMs at the level of < 20% from the 41 Malawian clinical specimens. Overall, DRMs were detected in 85% of pre-ART and 90.5% of virologic failure patients by deep sequencing. Among pre-ART patients, deep sequencing identified a statistically significant higher prevalence of DRMs to nucleoside reverse transcriptase inhibitors (NRTIs) compared with Sanger sequencing. The difference was mainly due to the high prevalence of minority K65R and M184I mutations. Most virologic failure patients harboured DRMs against both NRTIs and non-nucleoside reverse transcriptase inhibitors (NNRTIs). These minority DRMs contributed to the increased or enhanced virologic failures in these patients.

Conclusion

The results revealed the presence of minority DRMs to NRTIs and NNRTIs in specimens collected at baseline and virologic failure time points. These minority DRMs not only increased resistance levels to NRTIs and NNRTIs for the prescribed ART, but also expanded resistance to additional major first-line ART drugs. This study suggested that drug resistance testing that uses more sensitive technologies, is needed in this setting.

Combinatorial CRISPR-Cas9 and RNA Interference Attack on HIV-1 DNA and RNA Can Lead to Cross-Resistance

Many potent antiviral drugs have been developed against HIV-1, and their combined action is usually successful in achieving durable virus suppression in infected individuals. This success is based on two effects: additive or even synergistic virus inhibition and an increase in the genetic threshold for development of drug resistance. More recently, several genetic approaches have been developed to attack the HIV-1 genome in a gene therapy setting. We set out to test the combinatorial possibilities for a therapy based on the CRISPR-Cas9 and RNA interference (RNAi) mechanisms that attack the viral DNA and RNA, respectively. When two different sites in the HIV-1 genome were targeted, either with dual CRISPR-Cas9 antivirals or with a combination of CRISPR-Cas9 and RNAi antivirals, we observed additive inhibition, much like what was reported for antiviral drugs. However, when the same or overlapping viral sequence was attacked by the antivirals, rapid escape from a CRISPR-Cas9 antiviral, assisted by the error-prone nonhomologous end joining (NHEJ) DNA repair machinery, accelerated the development of cross-resistance to the other CRISPR-Cas9 or RNAi antiviral. Thus, genetic antiviral approaches can be combined, but overlap should be avoided.

Mutation and recombination in pathogen evolution: Relevance, methods and controversies

Mutation and recombination drive the evolution of most pathogens by generating the genetic variants upon which selection operates. Those variants can, for example, confer resistance to host immune systems and drug therapies or lead to epidemic outbreaks. Given their importance, diverse evolutionary studies have investigated the abundance and consequences of mutation and recombination in pathogen populations. However, some controversies persist regarding the contribution of each evolutionary force to the development of particular phenotypic observations (e.g., drug resistance). In this study, we revise the importance of mutation and recombination in the evolution of pathogens at both intra-host and inter-host levels. We also describe state-of-the-art analytical methodologies to detect and quantify these two evolutionary forces, including biases that are often ignored in evolutionary studies. Finally, we present some of our former studies involving pathogenic taxa where mutation and recombination played crucial roles in the recovery of pathogenic fitness, the generation of interspecific genetic diversity, or the design of centralized vaccines. This review also illustrates several common controversies and pitfalls in the analysis and in the evaluation and interpretation of mutation and recombination outcomes.

HIV drug resistance in HIV positive individuals under antiretroviral treatment in Shandong Province, China

The efficacy of antiretroviral drugs is limited by the development of drug resistance. Therefore, it is important to examine HIV drug resistance following the nationwide implementation of drug resistance testing in China since 2009. We conducted drug resistance testing in patients who were already on or new to HIV antiretroviral therapy (ART) in Shandong Province, China, from 2011 to 2013, and grouped them based on the presence or absence of drug resistance to determine the effects of age, gender, ethnicity, marital status, educational level, route of transmission and treatment status on drug resistance. We then examined levels of drug resistance the following year. The drug resistance rates of HIV patients on ART in Shandong from 2011 to 2013 were 3.45% (21/608), 3.38% (31/916), and 4.29% (54/1259), per year, respectively. M184V was the most frequently found point mutation, conferring resistance to the nucleoside reverse transcriptase inhibitor, while Y181C, G190A, K103N and V179D/E/F were the most frequent point mutations conferring resistance to the non-nucleoside reverse transcriptase inhibitor. In addition, the protease inhibitor drug resistance mutations I54V and V82A were identified for the first time in Shandong Province. Primary resistance accounts for 20% of the impact factors for drug resistance. Furthermore, it was found that educational level and treatment regimen were high-risk factors for drug resistance in 2011 (P<0.05), while treatment regimen was a high risk factor for drug resistance in 2012 and 2013 (P<0.05). Among the 106 drug-resistant patients, 77 received immediate adjustment of treatment regimen following testing, and 69 (89.6%) showed a reduction in drug resistance the following year. HIV drug resistance has a low prevalence in Shandong Province. However, patients on second line ART regimens and those with low educational level need continuous monitoring. Active drug resistance testing can effectively prevent the development of drug resistance.

HIV-1 drug resistance and resistance testing

The global scale-up of antiretroviral (ARV) therapy (ART) has led to dramatic reductions in HIV-1 mortality and incidence. However, HIV drug resistance (HIVDR) poses a potential threat to the long-term success of ART and is emerging as a threat to the elimination of AIDS as a public health problem by 2030. In this review we describe the genetic mechanisms, epidemiology, and management of HIVDR at both individual and population levels across diverse economic and geographic settings. To describe the genetic mechanisms of HIVDR, we review the genetic barriers to resistance for the most commonly used ARVs and describe the extent of cross-resistance between them. To describe the epidemiology of HIVDR, we summarize the prevalence and patterns of transmitted drug resistance (TDR) and acquired drug resistance (ADR) in both high-income and low- and middle-income countries (LMICs). We also review to two categories of HIVDR with important public health relevance: (i) pre-treatment drug resistance (PDR), a World Health Organization-recommended HIVDR surveillance metric and (ii) and pre-exposure prophylaxis (PrEP)-related drug resistance, a type of ADR that can impact clinical outcomes if present at the time of treatment initiation. To summarize the implications of HIVDR for patient management, we review the role of genotypic resistance testing and treatment practices in both high-income and LMIC settings. In high-income countries where drug resistance testing is part of routine care, such an understanding can help clinicians prevent virological failure and accumulation of further HIVDR on an individual level by selecting the most efficacious regimens for their patients. Although there is reduced access to diagnostic testing and to many ARVs in LMIC, understanding the scientific basis and clinical implications of HIVDR is useful in all regions in order to shape appropriate surveillance, inform treatment algorithms, and manage difficult cases.

Population genomics of intrapatient HIV-1 evolution

Many microbial populations rapidly adapt to changing environments with multiple variants competing for survival. To quantify such complex evolutionary dynamics in vivo, time resolved and genome wide data including rare variants are essential. We performed whole-genome deep sequencing of HIV-1 populations in 9 untreated patients, with 6-12 longitudinal samples per patient spanning 5-8 years of infection. The data can be accessed and explored via an interactive web application. We show that patterns of minor diversity are reproducible between patients and mirror global HIV-1 diversity, suggesting a universal landscape of fitness costs that control diversity. Reversions towards the ancestral HIV-1 sequence are observed throughout infection and account for almost one third of all sequence changes. Reversion rates depend strongly on conservation. Frequent recombination limits linkage disequilibrium to about 100 bp in most of the genome, but strong hitch-hiking due to short range linkage limits diversity.

Recombination in viruses: mechanisms, methods of study, and evolutionary consequences

Recombination is a pervasive process generating diversity in most viruses. It joins variants that arise independently within the same molecule, creating new opportunities for viruses to overcome selective pressures and to adapt to new environments and hosts. Consequently, the analysis of viral recombination attracts the interest of clinicians, epidemiologists, molecular biologists and evolutionary biologists. In this review we present an overview of three major areas related to viral recombination: (i) the molecular mechanisms that underlie recombination in model viruses, including DNA-viruses (Herpesvirus) and RNA-viruses (Human Influenza Virus and Human Immunodeficiency Virus), (ii) the analytical procedures to detect recombination in viral sequences and to determine the recombination breakpoints, along with the conceptual and methodological tools currently used and a brief overview of the impact of new sequencing technologies on the detection of recombination, and (iii) the major areas in the evolutionary analysis of viral populations on which recombination has an impact. These include the evaluation of selective pressures acting on viral populations, the application of evolutionary reconstructions in the characterization of centralized genes for vaccine design, and the evaluation of linkage disequilibrium and population structure.

Drug resistance and plasma viral RNA level after ineffective use of oral pre-exposure prophylaxis in women

BACKGROUND

Pre-exposure prophylaxis (PrEP) with daily oral emtricitabine (FTC)/tenofovir disoproxil fumarate may select for drug resistance if there is low adherence.

METHODS

Plasma viral HIV-1 RNA level, CD4+ T-cell counts, and drug resistance were evaluated among seroconverting women in the FEM-PrEP trial (clinicaltrials.gov NCT00625404) using standard clinical tests, allele-specific PCR (ASPCR), and by deep sequencing. Tenofovir, FTC, and their intracellular metabolites were measured in plasma and cells.

RESULTS

There was no difference in plasma HIV-1 RNA level or CD4+ cell count among seroconverters in the active arm versus those receiving placebo. Tenofovir resistance was not observed. FTC resistance was detected using clinical assays in five seroconverters (four in the active arm and one in the placebo arm); two in the active arm occurred among women having moderate concentrations of PrEP drugs in the blood. The first evidence of infection occurred at the first postenrollment visit in three of the four with FTC resistance, although none had detectable viral nucleic acids at enrollment. FTC-resistant minor variants were detected in an additional four seroconverters (one in the active arm and three in the placebo arm).

CONCLUSIONS

Drug resistance detected during ineffective PrEP use had characteristics suggesting transmitted infection or incubating infection prior to starting PrEP.

Estimating HIV-1 fitness characteristics from cross-sectional genotype data

Despite the success of highly active antiretroviral therapy (HAART) in the management of human immunodeficiency virus (HIV)-1 infection, virological failure due to drug resistance development remains a major challenge. Resistant mutants display reduced drug susceptibilities, but in the absence of drug, they generally have a lower fitness than the wild type, owing to a mutation-incurred cost. The interaction between these fitness costs and drug resistance dictates the appearance of mutants and influences viral suppression and therapeutic success. Assessing in vivo viral fitness is a challenging task and yet one that has significant clinical relevance. Here, we present a new computational modelling approach for estimating viral fitness that relies on common sparse cross-sectional clinical data by combining statistical approaches to learn drug-specific mutational pathways and resistance factors with viral dynamics models to represent the host-virus interaction and actions of drug mechanistically. We estimate in vivo fitness characteristics of mutant genotypes for two antiretroviral drugs, the reverse transcriptase inhibitor zidovudine (ZDV) and the protease inhibitor indinavir (IDV). Well-known features of HIV-1 fitness landscapes are recovered, both in the absence and presence of drugs. We quantify the complex interplay between fitness costs and resistance by computing selective advantages for different mutants. Our approach extends naturally to multiple drugs and we illustrate this by simulating a dual therapy with ZDV and IDV to assess therapy failure. The combined statistical and dynamical modelling approach may help in dissecting the effects of fitness costs and resistance with the ultimate aim of assisting the choice of salvage therapies after treatment failure.

Mutations conferring drug resistance represent major threats to the therapeutic success of highly active antiretroviral therapy (HAART) against human immunodeficiency virus (HIV)-1 infection. Viral mutants differ in their fitness and assessing viral fitness is a challenging task. In this article, we estimate drug-specific mutational pathways by learning from clinical data using statistical techniques and incorporate these into mathematical models of in vivo viral infection dynamics. This approach enables us to estimate mutant fitness characteristics. We illustrate our method by predicting fitness characteristics of mutant genotypes for two different antiretroviral therapies with the drugs zidovudine and indinavir. We recover several established features of mutant fitnesses and quantify fitness characteristics both in the absence and presence of drugs. Our model extends naturally to multiple drugs and we illustrate this by simulating a dual therapy with ZDV and IDV to assess therapy failure. Additionally, our modelling approach relies only on cross-sectional clinical data. We believe that such an approach is a highly valuable tool in assisting the choice of salvage therapies after treatment failure.

Exposure to MIV-150 from a high-dose intravaginal ring results in limited emergence of drug resistance mutations in SHIV-RT infected rhesus macaques

When microbicides used for HIV prevention contain antiretroviral drugs, there is concern for the potential emergence of drug-resistant HIV following use in infected individuals who are either unaware of their HIV infection status or who are aware but still choose to use the microbicide. Resistant virus could ultimately impact their responsiveness to treatment and/or result in subsequent transmission of drug-resistant virus. We tested whether drug resistance mutations (DRMs) would emerge in macaques infected with simian immunodeficiency virus expressing HIV reverse transcriptase (SHIV-RT) after sustained exposure to the potent non-nucleoside reverse transcriptase inhibitor (NNRTI) MIV-150 delivered via an intravaginal ring (IVR). We first treated 4 SHIV-RT-infected animals with daily intramuscular injections of MIV-150 over two 21 day (d) intervals separated by a 7 d drug hiatus. In all 4 animals, NNRTI DRMs (single and combinations) were detected within 14 d and expanded in proportion and diversity with time. Knowing that we could detect in vivo emergence of NNRTI DRMs in response to MIV-150, we then tested whether a high-dose MIV-150 IVR (loaded with >10 times the amount being used in a combination microbicide IVR in development) would select for resistance in 6 infected animals, modeling use of this prevention method by an HIV-infected woman. We previously demonstrated that this MIV-150 IVR provides significant protection against vaginal SHIV-RT challenge. Wearing the MIV-150 IVR for 56 d led to only 2 single DRMs in 2 of 6 animals (430 RT sequences analyzed total, 0.46%) from plasma and lymph nodes despite MIV-150 persisting in the plasma, vaginal fluids, and genital tissues. Only wild type virus sequences were detected in the genital tissues. These findings indicate a low probability for the emergence of DRMs after topical MIV-150 exposure and support the advancement of MIV-150-containing microbicides.

rtM204Q may serve as a novel lamivudine-resistance-associated mutation of hepatitis B virus

Background and Aims

Lamivudine (LAM) is still widely used for anti-HBV therapy in China. The study aimed to clarify whether a newly-found rtM204Q mutation from patients was associated with the drug resistance.

Methods

HBV complete reverse-transcriptase region was screened by direct sequencing and verified by clonal sequencing. Replication-competent plasmids containing patient-derived 1.1mer mutant or wild-type viral genome were constructed and transfected into HepG2 cells. After cultured with or without serially-diluted antiviral drugs, intracellular HBV replicative intermediates were quantitated for calculating the 50% effective concentration of drug (EC₅₀).

Results

A total of 12,000 serum samples of 9,830 patients with chronic HBV infection were screened. rtM204Q mutation was detected in seven LAM-refractory patients. By contrast, rtM204I/rtM204V mutations were detected in 2,502 patients' samples. The rtM204Q emerged either alone or in concomitance with rtM204I/rtM204V, and all were accompanied with virologic breakthrough in clinical course. Clonal sequencing verified that rtM204Q mutant was predominant in viral quasispecies of these samples. Phenotypic analysis showed that rtM204Q mutant had 89.9% of replication capacity and 76-fold increased LAM EC₅₀ of the concomitant wild-type strain. By contrast, rtM204I mutant in the sample had lower replication capacity and higher LAM resistance (46.3% and 1396-fold increased LAM EC₅₀ of the wild-type strain) compared to rtM204Q mutant. rtM204Q mutant was susceptible to adefovir dipivoxil (ADV) in vitro and ADV/ADV+LAM rescue therapy in clinic.

Conclusion

rtM204Q is suggested to be a novel LAM-resistance-associated mutation. It conferred a moderate resistance with higher competent natural replication capacity compared to rtM204I mutation.

Loss and recovery of genetic diversity in adapting populations of HIV

The evolution of drug resistance in HIV occurs by the fixation of specific, well-known, drug-resistance mutations, but the underlying population genetic processes are not well understood. By analyzing within-patient longitudinal sequence data, we make four observations that shed a light on the underlying processes and allow us to infer the short-term effective population size of the viral population in a patient. Our first observation is that the evolution of drug resistance usually occurs by the fixation of one drug-resistance mutation at a time, as opposed to several changes simultaneously. Second, we find that these fixation events are accompanied by a reduction in genetic diversity in the region surrounding the fixed drug-resistance mutation, due to the hitchhiking effect. Third, we observe that the fixation of drug-resistance mutations involves both hard and soft selective sweeps. In a hard sweep, a resistance mutation arises in a single viral particle and drives all linked mutations with it when it spreads in the viral population, which dramatically reduces genetic diversity. On the other hand, in a soft sweep, a resistance mutation occurs multiple times on different genetic backgrounds, and the reduction of diversity is weak. Using the frequency of occurrence of hard and soft sweeps we estimate the effective population size of HIV to be 1.5 x 10(5) (95% confidence interval [0.8 x 10(5),4.8 x 10(5)]). This number is much lower than the actual number of infected cells, but much larger than previous population size estimates based on synonymous diversity. We propose several explanations for the observed discrepancies. Finally, our fourth observation is that genetic diversity at non-synonymous sites recovers to its pre-fixation value within 18 months, whereas diversity at synonymous sites remains depressed after this time period. These results improve our understanding of HIV evolution and have potential implications for treatment strategies.

The role of backward mutations on the within-host dynamics of HIV-1

The quality of life for patients infected with human immunodeficiency virus (HIV-1) has been positively impacted by the use of antiretroviral therapy (ART). However, the benefits of ART are usually halted by the emergence of drug resistance. Drug-resistant strains arise from virus mutations, as HIV-1 reverse transcription is prone to errors, with mutations normally carrying fitness costs to the virus. When ART is interrupted, the wild-type drug-sensitive strain rapidly out-competes the resistant strain, as the former strain is fitter than the latter in the absence of ART. One mechanism for sustaining the sensitive strain during ART is given by the virus mutating from resistant to sensitive strains, which is referred to as backward mutation. This is important during periods of treatment interruptions as prior existence of the sensitive strain would lead to replacement of the resistant strain. In order to assess the role of backward mutations in the dynamics of HIV-1 within an infected host, we analyze a mathematical model of two interacting virus strains in either absence or presence of ART. We study the effect of backward mutations on the definition of the basic reproductive number, and the value and stability of equilibrium points. The analysis of the model shows that, thanks to both forward and backward mutations, sensitive and resistant strains co-exist. In addition, conditions for the dominance of a viral strain with or without ART are provided. For this model, backward mutations are shown to be necessary for the persistence of the sensitive strain during ART.

HIV populations are large and accumulate high genetic diversity in a nonlinear fashion

HIV infection is characterized by rapid and error-prone viral replication resulting in genetically diverse virus populations. The rate of accumulation of diversity and the mechanisms involved are under intense study to provide useful information to understand immune evasion and the development of drug resistance. To characterize the development of viral diversity after infection, we carried out an in-depth analysis of single genome sequences of HIV pro-pol to assess diversity and divergence and to estimate replicating population sizes in a group of treatment-naive HIV-infected individuals sampled at single (n = 22) or multiple, longitudinal (n = 11) time points. Analysis of single genome sequences revealed nonlinear accumulation of sequence diversity during the course of infection. Diversity accumulated in recently infected individuals at rates 30-fold higher than in patients with chronic infection. Accumulation of synonymous changes accounted for most of the diversity during chronic infection. Accumulation of diversity resulted in population shifts, but the rates of change were low relative to estimated replication cycle times, consistent with relatively large population sizes. Analysis of changes in allele frequencies revealed effective population sizes that are substantially higher than previous estimates of approximately 1,000 infectious particles/infected individual. Taken together, these observations indicate that HIV populations are large, diverse, and slow to change in chronic infection and that the emergence of new mutations, including drug resistance mutations, is governed by both selection forces and drift.

Computational inference methods for selective sweeps arising in acute HIV infection

During the first weeks of human immunodeficiency virus-1 (HIV-1) infection, cytotoxic T-lymphocytes (CTLs) select for multiple escape mutations in the infecting HIV population. In recent years, methods that use escape mutation data to estimate rates of HIV escape have been developed, thereby providing a quantitative framework for exploring HIV escape from CTL response. Current methods for escape-rate inference focus on a specific HIV mutant selected by a single CTL response. However, recent studies have shown that during the first weeks of infection, CTL responses occur at one to three epitopes and HIV escape occurs through complex mutation pathways. Consequently, HIV escape from CTL response forms a complex, selective sweep that is difficult to analyze. In this work, we develop a model of initial infection, based on the well-known standard model, that allows for a description of multi-epitope response and the complex mutation pathways of HIV escape. Under this model, we develop Bayesian and hypothesis-test inference methods that allow us to analyze and estimate HIV escape rates. The methods are applied to two HIV patient data sets, concretely demonstrating the utility of our approach.

Design requirements for interfering particles to maintain coadaptive stability with HIV-1

Defective interfering particles (DIPs) are viral deletion mutants lacking essential transacting or packaging elements and must be complemented by wild-type virus to propagate. DIPs transmit through human populations, replicating at the expense of the wild-type virus and acting as molecular parasites of viruses. Consequently, engineered DIPs have been proposed as therapies for a number of diseases, including human immunodeficiency virus (HIV). However, it is not clear if DIP-based therapies would face evolutionary blocks given the high mutation rates and high within-host diversity of lentiviruses. Divergent evolution of HIV and DIPs appears likely since natural DIPs have not been detected for lentiviruses, despite extensive sequencing of HIVs and simian immunodeficiency viruses (SIVs). Here, we tested if the apparent lack of lentiviral DIPs is due to natural selection and analyzed which molecular characteristics a DIP or DIP-based therapy would need to maintain coadaptive stability with HIV-1. Using a well-established mathematical model of HIV-1 in a host extended to include its replication in a single cell and interference from DIP, we calculated evolutionary selection coefficients. The analysis predicts that interference by codimerization between DIPs and HIV-1 genomes is evolutionarily unstable, indicating that recombination between DIPs and HIV-1 would be selected against. In contrast, DIPs that interfere via competition for capsids have the potential to be evolutionarily stable if the capsid-to-genome production ratio of HIV-1 is >1. Thus, HIV-1 variants that attempt to "starve" DIPs to escape interference would be selected against. In summary, the analysis suggests specific experimental measurements that could address the apparent lack of naturally occurring lentiviral DIPs and specifies how therapeutic approaches based on engineered DIPs could be evolutionarily robust and avoid recombination.

Significant reductions in Gag-protease-mediated HIV-1 replication capacity during the course of the epidemic in Japan

Human immunodeficiency virus type 1 (HIV-1) evolves rapidly in response to host immune selection pressures. As a result, the functional properties of HIV-1 isolates from earlier in the epidemic may differ from those of isolates from later stages. However, few studies have investigated alterations in viral replication capacity (RC) over the epidemic. In the present study, we compare Gag-Protease-associated RC between early and late isolates in Japan (1994 to 2009). HIV-1 subtype B sequences from 156 antiretroviral-naïve Japanese with chronic asymptomatic infection were used to construct a chimeric NL4-3 strain encoding plasma-derived gag-protease. Viral replication capacity was examined by infecting a long terminal repeat-driven green fluorescent protein-reporter T cell line. We observed a reduction in the RC of chimeric NL4-3 over the epidemic, which remained significant after adjusting for the CD4(+) T cell count and plasma virus load. The same outcome was seen when limiting the analysis to a single large cluster of related sequences, indicating that our results are not due to shifts in the molecular epidemiology of the epidemic in Japan. Moreover, the change in RC was independent of genetic distance between patient-derived sequences and wild-type NL4-3, thus ruling out potential temporal bias due to genetic similarity between patient and historic viral backbone sequences. Collectively, these data indicate that Gag-Protease-associated HIV-1 replication capacity has decreased over the epidemic in Japan. Larger studies from multiple geographical regions will be required to confirm this phenomenon.

Stochastic effects are important in intrahost HIV evolution even when viral loads are high

Blood plasma viral loads and the time to progress to AIDS differ widely among untreated HIV-infected humans. Although people with certain HLA (HLA-I) alleles are more likely to control HIV infections without therapy, the majority of such untreated individuals exhibit high viral loads and progress to AIDS. Stochastic effects are considered unimportant for evolutionary dynamics in HIV-infected people when viral load is high or when selective forces strongly drive mutation. We describe a computational study of host-pathogen interaction demonstrating that stochastic effects can have a profound influence on disease dynamics, even in cases of high viral load and strong selective pressure. These stochastic effects are pronounced when the virus must traverse a fitness "barrier" in sequence space to escape the host's cytotoxic T-lymphocyte (CTL) response, as often occurs when a fitness defect imposed by a CTL-driven mutation must be compensated for by other mutations. These "barrier-crossing" events are infrequent and stochastic, resulting in divergent disease outcomes in genetically identical individuals infected by the same viral strain. Our results reveal how genetic determinants of the CTL response control the probability with which an individual is able to control HIV infection indefinitely, and thus provide clues for vaccine design.

The origin of genetic diversity in HIV-1

One of the hallmarks of HIV infection is the rapid development of a genetically complex population (quasispecies) from an initially limited number of infectious particles. Genetic diversity remains one of the major obstacles to eradication of HIV. The viral quasispecies can respond rapidly to selective pressures, such as that imposed by the immune system and antiretroviral therapy, and frustrates vaccine design efforts. Two unique features of retroviral replication are responsible for the unprecedented variation generated during infection. First, mutations are frequently introduced into the viral genome by the error prone viral reverse transcriptase and through the actions of host cellular factors, such as the APOBEC family of nucleic acid editing enzymes. Second, the HIV reverse transcriptase can utilize both copies of the co-packaged viral genome in a process termed retroviral recombination. When the co-packaged viral genomes are genetically different, retroviral recombination can lead to the shuffling of mutations between viral genomes in the quasispecies. This review outlines the stages of the retroviral life cycle where genetic variation is introduced, focusing on the principal mechanisms of mutation and recombination. Understanding the mechanistic origin of genetic diversity is essential to combating HIV.

The HIV Epidemic: High-Income Countries

The HIV epidemic in higher-income nations is driven by receptive anal intercourse, injection drug use through needle/syringe sharing, and, less efficiently, vaginal intercourse. Alcohol and noninjecting drug use increase sexual HIV vulnerability. Appropriate diagnostic screening has nearly eliminated blood/blood product-related transmissions and, with antiretroviral therapy, has reduced mother-to-child transmission radically. Affected subgroups have changed over time (e.g., increasing numbers of Black and minority ethnic men who have sex with men). Molecular phylogenetic approaches have established historical links between HIV strains from central Africa to those in the United States and thence to Europe. However, Europe did not just receive virus from the United States, as it was also imported from Africa directly. Initial introductions led to epidemics in different risk groups in Western Europe distinguished by viral clades/sequences, and likewise, more recent explosive epidemics linked to injection drug use in Eastern Europe are associated with specific strains. Recent developments in phylodynamic approaches have made it possible to obtain estimates of sequence evolution rates and network parameters for epidemics.

Emerging complexities of APOBEC3G action on immunity and viral fitness during HIV infection and treatment

The enzyme APOBEC3G (A3G) mutates the human immunodeficiency virus (HIV) genome by converting deoxycytidine (dC) to deoxyuridine (dU) on minus strand viral DNA during reverse transcription. A3G restricts viral propagation by degrading or incapacitating the coding ability of the HIV genome. Thus, this enzyme has been perceived as an innate immune barrier to viral replication whilst adaptive immunity responses escalate to effective levels. The discovery of A3G less than a decade ago led to the promise of new anti-viral therapies based on manipulation of its cellular expression and/or activity. The rationale for therapeutic approaches has been solidified by demonstration of the effectiveness of A3G in diminishing viral replication in cell culture systems of HIV infection, reports of its mutational footprint in virions from patients, and recognition of its unusually robust enzymatic potential in biochemical studies in vitro. Despite its effectiveness in various experimental systems, numerous recent studies have shown that the ability of A3G to combat HIV in the physiological setting is severely limited. In fact, it has become apparent that its mutational activity may actually enhance viral fitness by accelerating HIV evolution towards the evasion of both anti-viral drugs and the immune system. This body of work suggests that the role of A3G in HIV infection is more complex than heretofore appreciated and supports the hypothesis that HIV has evolved to exploit the action of this host factor. Here we present an overview of recent data that bring to light historical overestimation of A3G's standing as a strictly anti-viral agent. We discuss the limitations of experimental systems used to assess its activities as well as caveats in data interpretation.

The dynamics of HIV-1 adaptation in early infection

Human immunodeficiency virus type 1 (HIV-1) undergoes a severe population bottleneck during sexual transmission and yet adapts extremely rapidly to the earliest immune responses. The bottleneck has been inferred to typically consist of a single genome, and typically eight amino acid mutations in viral proteins spread to fixation by the end of the early chronic phase of infection in response to selection by CD8(+) T cells. Stochastic simulation was used to examine the effects of the transmission bottleneck and of potential interference among spreading immune-escape mutations on the adaptive dynamics of the virus in early infection. If major viral population genetic parameters are assigned realistic values that permit rapid adaptive evolution, then a bottleneck of a single genome is not inconsistent with the observed pattern of adaptive fixations. One requirement is strong selection by CD8(+) T cells that decreases over time. Such selection may reduce effective population sizes at linked loci through genetic hitchhiking. However, this effect is predicted to be minor in early infection because the transmission bottleneck reduces the effective population size to such an extent that the resulting strong selection and weak mutation cause beneficial mutations to fix sequentially and thus avoid interference.

The impact of molecular manipulation in residue 114 of human immunodeficiency virus type-1 reverse transcriptase on dNTP substrate binding and viral replication

Human immunodeficiency virus type-1 (HIV-1) reverse transcriptase (RT) has a unique tight binding to dNTP substrates. Structural modeling of Ala-114 of HIV-1 RT suggests that longer side chains at this residue can reduce the space normally occupied by the sugar moiety of an incoming dNTP. Indeed, mutations at Ala-114 decrease the ability of RT to synthesize DNA at low dNTP concentrations and reduce the dNTP-binding affinity (K(d)) of RT. However, the K(d) values of WT and A114C RT remained equivalent with an acyclic dNTP substrate. Finally, mutant A114 RT HIV-1 vectors displayed a greatly reduced transduction in nondividing human lung fibroblasts (HLFs), while WT HIV-1 vector efficiently transduced both dividing and nondividing HLFs. Together these data support that the A114 residue of HIV-1 RT plays a key mechanistic role in the dNTP binding of HIV-1 RT and the unique viral infectivity of target cell types with low dNTP pools.

The use of nonhuman primate models of HIV infection for the evaluation of antiviral strategies

Several nonhuman primate models are used in HIV/AIDS research. In contrast to natural host models, infection of macaques with virulent simian immunodeficiency virus (SIV) isolates results in a disease (simian AIDS) that closely resembles HIV infection and AIDS. Although there is no perfect animal model, and each of the available models has its limitations, a carefully designed study allows experimental approaches that are not feasible in humans, but that can provide better insights in disease pathogenesis and proof-of-concept of novel intervention strategies. In the early years of the HIV pandemic, nonhuman primate models played a minor role in the development of antiviral strategies. Since then, a better understanding of the disease and the development of better compounds and assays to monitor antiviral effects have increased the usefulness and relevance of these animal models in the preclinical development of HIV vaccines, microbicides, and antiretroviral drugs. Several strategies that were first discovered to have efficacy in nonhuman primate models are now increasingly used in humans. Recent trends include the use of nonhuman primate models to explore strategies that could reduce viral reservoirs and, ultimately, attempt to cure infection. Ongoing comparison of results obtained in nonhuman primate models with those observed in human studies will lead to further validation and improvement of these animal models so they can continue to advance our scientific knowledge and guide clinical trials.