Origin of minority drug-resistant HIV-1 variants in primary HIV-1 infection

Cohort Profile: The Zurich Primary HIV Infection Study

The Zurich Primary HIV Infection (ZPHI) study is a longitudinal cohort study established in 2002, aiming to study the clinical, epidemiological, and biological characteristics of primary HIV infection. The ZPHI enrolls individuals with documented primary HIV-1 infection. At the baseline and thereafter, the socio-demographic, clinical, and laboratory data are systematically collected, and regular blood sampling is performed for biobanking. By the end of December 2022, 486 people were enrolled, of which 353 were still undergoing active follow-up. Of the 486 participants, 86% had an acute infection, and 14% a recent HIV-1 infection. Men who have sex with men accounted for 74% of the study population. The median time from the estimated date of infection to diagnosis was 32 days. The median time from diagnosis to the initiation of antiretroviral therapy was 11 days, and this has consistently decreased over the last two decades. During the seroconversion phase, 447 (92%) patients reported having symptoms, of which only 73% of the patients were classified as having typical acute retroviral syndrome. The ZPHI study is a well-characterized cohort belonging to the most extensively studied primary HIV infection cohort. Its findings contribute to advancing our understanding of the early stages of HIV infection and pathogenesis, and it is paving the way to further improve HIV translational research and HIV medicine.

From viral democratic genomes to viral wild bunch of quasispecies

The tremendous majority of RNA genomes from pathogenic viruses analyzed and deposited in databases are consensus or "democratic" genomes. They represent the genomes most frequently found in the clinical samples of patients but do not account for the huge genetic diversity of coexisting genomes, which is better described as quasispecies. A viral quasispecies is defined as the dynamic distribution of nonidentical but closely related mutants, variants, recombinant, or reassortant viral genomes. Viral quasispecies have collective behavior and dynamics and are the subject of internal interactions that comprise interference, complementation, or cooperation. In the setting of SARS-CoV-2 infection, intrahost SARS-CoV-2 genetic diversity was recently notably reported for immunocompromised, chronically infected patients, for patients treated with monoclonal antibodies targeting the viral spike protein, and for different body compartments of a single patient. A question that deserves attention is whether such diversity is generated postinfection from a clonal genome in response to selection pressure or is already present at the time of infection as a quasispecies. In the present review, we summarize the data supporting that hosts are infected by a "wild bunch" of viruses rather than by multiple virions sharing the same genome. Each virion in the "wild bunch" may have different virulence and tissue tropisms. As the number of viruses replicated during host infections is huge, a viral quasispecies at any time of infection is wide and is also influenced by host-specific selection pressure after infection, which accounts for the difficulty in deciphering and predicting the appearance of more fit variants and the evolution of epidemics of novel RNA viruses.

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.

Transmitted HIV-1 drug resistance in a large international cohort using next-generation sequencing: results from the Strategic Timing of Antiretroviral Treatment (START) study

OBJECTIVES

The aim of this analysis was to characterize transmitted drug resistance (TDR) in Strategic Timing of Antiretroviral Treatment (START) study participants by next-generation sequencing (NGS), a sensitive assay capable of detecting low-frequency variants.

METHODS

Stored plasma from participants with entry HIV RNA > 1000 copies/mL were analysed by NGS (Illumina MiSeq). TDR was based on the WHO 2009 surveillance definition with the addition of reverse transcriptase (RT) mutations T215N and E138K, and integrase strand transfer inhibitor (INSTI) surveillance mutations (Stanford HIVdb). Drug resistance mutations (DRMs) detected at three thresholds are reported: > 2%, 5% and 20% of the viral population.

RESULTS

Between 2009 and 2013, START enrolled 4684 antiretroviral therapy (ART)-naïve individuals in 35 countries. Baseline NGS data at study entry were available for 2902 participants. Overall prevalence rates of TDR using a detection threshold of 2%/5%/20% were 9.2%/5.6%/3.2% for nucleoside reverse transcriptase inhibitors (NRTIs), 9.2%/6.6%/4.9% for non-NRTIs, 11.4%/5.5%/2.4% for protease inhibitors (PIs) and 3.5%/1.6%/0.1% for INSTI DRMs and varied by geographic region. Using the 2% detection threshold, individual DRMs with the highest prevalence were: PI M46IL (5.5%), RT K103NS (3.5%), RT G190ASE (3.1%), T215ISCDVEN (2.5%), RT M41L (2.2%), RT K219QENR (1.7%) and PI D30N (1.6%). INSTI DRMs were detected almost exclusively below the 20% detection threshold, most commonly Y143H (0.4%), Q148R (0.4%) and T66I (0.4%).

CONCLUSIONS

Use of NGS in this study population resulted in the detection of a large proportion of low-level variants which would not have been detected by traditional Sanger sequencing. Global surveillance studies utilizing NGS should provide a more comprehensive assessment of TDR prevalence in different regions of the world.

Low-Abundance Drug-Resistant HIV-1 Variants in Antiretroviral Drug-Naive Individuals: A Systematic Review of Detection Methods, Prevalence, and Clinical Impact

BACKGROUND

The presence of high-abundance drug-resistant HIV-1 jeopardizes success of antiretroviral therapy (ART). Despite numerous investigations, the clinical impact of low-abundance drug-resistant HIV-1 variants (LA-DRVs) at levels <15%-25% of the virus population in antiretroviral (ARV) drug-naive individuals remains controversial.

METHODS

We systematically reviewed 103 studies assessing prevalence, detection methods, technical and clinical detection cutoffs, and clinical significance of LA-DRVs in antiretroviral drug-naive adults.

RESULTS

In total, 14 919 ARV drug-naive individuals were included. Prevalence of LA-DRVs (ie, proportion of individuals harboring LA-DRVs) was 0%-100%. Technical detection cutoffs showed a 4 log range (0.001%-10%); 42/103 (40.8%) studies investigating the impact of LA-DRVs on ART; 25 studies included only individuals on first-line nonnucleoside reverse transcriptase inhibitor-based ART regimens. Eleven of those 25 studies (44.0%) reported a significantly association between preexisting LA-DRVs and risk of virological failure whereas 14/25 (56.0%) did not.

CONCLUSIONS

Comparability of the 103 studies is hampered by high heterogeneity of the studies' designs and use of different methods to detect LA-DRVs. Thus, evaluating clinical impact of LA-DRVs on first-line ART remains challenging. We, the WHO HIVResNet working group, defined central areas of future investigations to guide further efforts to implement ultrasensitive resistance testing in routine settings.

Determining the Origins of Human Immunodeficiency Virus Type 1 Drug-resistant Minority Variants in People Who Are Recently Infected Using Phylogenetic Reconstruction

BACKGROUND

Drug-resistant minority variants (DRMinVs) detected in patients who recently acquired human immunodeficiency virus type 1 (HIV-1) can be transmitted, generated de novo through virus replication, or technical errors. The first form is likely to persist and result in treatment failure, while the latter two could be stochastic and transient.

METHODS

Ultradeep sequencing of plasma samples from 835 individuals with recent HIV-1 infection in the United Kingdom was performed to detect DRMinVs at a mutation frequency between 2% and 20%. Sequence alignments including >110 000 HIV-1 partial pol consensus sequences from the UK HIV Drug Resistance Database (UK-HDRD), linked to epidemiological and clinical data from the HIV and AIDS Reporting System, were used for transmission cluster analysis. Transmission clusters were identified using Cluster Picker with a clade support of >90% and maximum genetic distances of 4.5% or 1.5%, the latter to limit detection to likely direct transmission events.

RESULTS

Drug-resistant majority variants (DRMajVs) were detected in 66 (7.9%) and DRMinVs in 84 (10.1%) of the recently infected individuals. High levels of clustering to sequences in UK-HDRD were observed for both DRMajV (n = 48; 72.7%) and DRMinV (n = 63; 75.0%) sequences. Of these, 43 (65.2%) with DRMajVs were in a transmission cluster with sequences that harbored the same DR mutation compared to only 3 (3.6%) sequences with DRMinVs (P < .00001, Fisher exact test). Evidence of likely direct transmission of DRMajVs was observed for 25/66 (37.9%), whereas none were observed for the DRMinVs (P < .00001).

CONCLUSIONS

Using a densely sampled HIV-infected population, we show no evidence of DRMinV transmission among recently infected individuals.

Pretreatment HIV drug resistance spread within transmission clusters in Mexico City

BACKGROUND

Pretreatment HIV drug resistance (HIVDR) to NNRTIs has consistently increased in Mexico City during the last decade.

OBJECTIVES

To infer the HIV genetic transmission network in Mexico City to describe the dynamics of the local HIV epidemic and spread of HIVDR.

PATIENTS AND METHODS

HIV pol sequences were obtained by next-generation sequencing from 2447 individuals before initiation of ART at the largest HIV clinic in Mexico City (April 2016 to June 2018). Pretreatment HIVDR was estimated using the Stanford algorithm at a Sanger-like threshold (≥20%). Genetic networks were inferred with HIV-TRACE, establishing putative transmission links with genetic distances <1.5%. We examined demographic associations among linked individuals with shared drug resistance mutations (DRMs) using a ≥ 2% threshold to include low-frequency variants.

RESULTS

Pretreatment HIVDR reached 14.8% (95% CI 13.4%-16.2%) in the cohort overall and 9.6% (8.5%-10.8%) to NNRTIs. Putative links with at least one other sequence were found for 963/2447 (39%) sequences, forming 326 clusters (2-20 individuals). The inferred network was assortative by age and municipality (P < 0.001). Clustering individuals were younger [adjusted OR (aOR) per year = 0.96, 95% CI 0.95-0.97, P < 0.001] and less likely to include women (aOR = 0.46, 95% CI 0.28-0.75, P = 0.002). Among clustering individuals, 175/963 (18%) shared DRMs (involving 66 clusters), of which 66/175 (38%) shared K103N/S (24 clusters). Eight municipalities (out of 75) harboured 65% of persons sharing DRMs. Among all persons sharing DRMs, those sharing K103N were younger (aOR = 0.93, 95% CI 0.88-0.98, P = 0.003).

CONCLUSIONS

Our analyses suggest age- and geographically associated transmission of DRMs within the HIV genetic network in Mexico City, warranting continuous monitoring and focused interventions.

Caution is needed in interpreting HIV transmission chains by ultradeep sequencing

OBJECTIVES

Molecular epidemiology is applied to various aspects of HIV transmission analyses. With ultradeep sequencing (UDS), in-depth characterization of transmission episodes involving minority variants is permitted. We explored HIV-1 epidemiological linkage and evaluated characteristics of transmission dynamics and transmitted drug resistance (TDR) detection through the added value of UDS.

DESIGN

HIV pol gene fragments were sequenced by UDS and Sanger sequencing on samples of 70 HIV-1-infected, treatment-naive recently diagnosed MSM.

METHODS

Pairwise genetic distances and maximum likelihood phylogenies were computed. Transmission events were identified as clades with branch support at least 70% and intraclade genetic difference less than 4.5%. TDR mutations were recognized from the TDR consensus list. Transmission directionality, directness and inoculum size were inferred from tree topologies.

RESULTS

Both datasets concurred in the identification of seven transmission pairs and one cluster of three patients. With UDS, direction of transmission was inferred in four out of eight chains. Evidence for multiple founder viruses was found in two out of eight chains. No transmission of minority-resistant variants was evidenced. TDR mutations prevalence in protease and reverse transcriptase fragments was 4.3% with Sanger sequencing and 18.6% with UDS.

CONCLUSION

Although Sanger sequencing and UDS identified the same transmission chains, UDS provided additional information on founder viruses, direction of transmission and levels of TDR. Nevertheless, topology of clusters was not always consistent across gene fragments, calling for a cautious interpretation of the data. Moreover, unobserved intermediary links cannot be excluded. Phylogenetic analysis use as a forensic technique for HIV transmission investigations is risky.

Focus on recently developed assays for detection of resistance/sensitivity to reverse transcriptase inhibitors

The biology of HIV is rather complex due to high rate of replication, frequent recombination, and introduction of mutations. This gives rise to a number of distinct variants referred as quasispecies. In addition, the latency within reservoir allows the periodic reactivation of virus replication. The rapid replication of HIV allows immune response escape and establishment of resistance to therapy that can be acquired through drug selection and/or transmitted among individuals. This prompted, over the years, the development of a range of assays aimed to determine drug resistance and sensitivity, to be used both in clinical practice and in antiviral research. Reverse transcriptase (RT) inhibitors have an eminent place among the anti-HIV drugs, being constantly present from the beginning until today in the most commonly used antiviral regimens. This mini-review seeks to provide an up-to-date overview of recent efforts in developing even more reliable and simple methods, of both genotypic and phenotypic types, for specifically detecting drug resistance and sensitivity to RT inhibitors.

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.

Lower pre-ART intra-participant HIV-1 pol diversity may not be associated with virologic failure in adults

Background

Identifying pre-ART factors associated with the emergence of HIV-1 drug resistance is critical for optimizing strategies to prevent virologic failure. A previous study reported that lower pre-ART HIV-1 pol diversity was associated with higher risk of virologic failure in HIV-1-infected children. To investigate this association in adults, we measured HIV-1 diversity with deep sequencing in pre-ART samples from adults with well-characterized virologic outcomes in a study (A5142) of initial ART conducted by the AIDS Clinical Trials Group (ACTG).

Methods

We identified 22 cases in ACTG A5142 who experienced virologic failure with drug resistance mutations in RT and 44 matched controls who did not experience virologic failure. cDNA was synthesized from plasma HIV-1 RNA. Each cDNA molecule was tagged with a unique primer ID and RT codons 41–103 were amplified and deep sequenced. Sequences with the same tag were aligned and a consensus was generated to reduce PCR and sequencing errors. Diversity was calculated by measuring average pairwise distance (APD) of the consensus sequences. An exact conditional logistic regression model with percent APD as the risk factor estimated the odds ratio for VF and the corresponding 95% confidence interval.

Results

Consensus single-genome sequences and diversity estimates of pol were obtained for pre-ART samples from 21 cases and 42 controls. The median (IQR) pre-ART percent APD was 0.71 (0.31–1.13) in cases and 0.58 (0.32–0.94) in controls. A possible trend was found for higher diversity being associated with greater risk of virologic failure in adults (OR = 2.2 per one percent APD increase, 95% CI = [0.8, 7.2]; p = 0.15).

Conclusions

This study in adults suggests there is a positive association between higher pre-ART pol diversity and the risk of virologic failure in adults rather than an inverse relationship reported in children.

No Substantial Evidence for Sexual Transmission of Minority HIV Drug Resistance Mutations in Men Who Have Sex with Men

During primary HIV infection, the presence of minority drug resistance mutations (DRM) may be a consequence of sexual transmission, de novo mutations, or technical errors in identification. Baseline blood samples were collected from 24 HIV-infected antiretroviral-naive, genetically and epidemiologically linked source and recipient partners shortly after the recipient's estimated date of infection. An additional 32 longitudinal samples were available from 11 recipients. Deep sequencing of HIV reverse transcriptase (RT) was performed (Roche/454), and the sequences were screened for nucleoside and nonnucleoside RT inhibitor DRM. The likelihood of sexual transmission and persistence of DRM was assessed using Bayesian-based statistical modeling. While the majority of DRM (>20%) were consistently transmitted from source to recipient, the probability of detecting a minority DRM in the recipient was not increased when the same minority DRM was detected in the source (Bayes factor [BF] = 6.37). Longitudinal analyses revealed an exponential decay of DRM (BF = 0.05) while genetic diversity increased. Our analysis revealed no substantial evidence for sexual transmission of minority DRM (BF = 0.02). The presence of minority DRM during early infection, followed by a rapid decay, is consistent with the "mutation-selection balance" hypothesis, in which deleterious mutations are more efficiently purged later during HIV infection when the larger effective population size allows more efficient selection. Future studies using more recent sequencing technologies that are less prone to single-base errors should confirm these results by applying a similar Bayesian framework in other clinical settings.IMPORTANCE The advent of sensitive sequencing platforms has led to an increased identification of minority drug resistance mutations (DRM), including among antiretroviral therapy-naive HIV-infected individuals. While transmission of DRM may impact future therapy options for newly infected individuals, the clinical significance of the detection of minority DRM remains controversial. In the present study, we applied deep-sequencing techniques within a Bayesian hierarchical framework to a cohort of 24 transmission pairs to investigate whether minority DRM detected shortly after transmission were the consequence of (i) sexual transmission from the source, (ii) de novo emergence shortly after infection followed by viral selection and evolution, or (iii) technical errors/limitations of deep-sequencing methods. We found no clear evidence to support the sexual transmission of minority resistant variants, and our results suggested that minor resistant variants may emerge de novo shortly after transmission, when the small effective population size limits efficient purge by natural selection.

Commonly Transmitted HIV-1 Drug Resistance Mutations in Reverse-Transcriptase and Protease in Antiretroviral Treatment-Naive Patients and Response to Regimens Containing Tenofovir Disoproxil Fumarate or Tenofovir Alafenamide

Background

The presence of transmitted drug resistance mutations (TDRMs) in antiretroviral treatment (ART)-naive patients can adversely affect the outcome of ART.

Methods

Resistance testing was conducted in 6704 ART-naive subjects predominantly from the United States and Europe in 9 clinical studies conducted by Gilead Sciences from 2000 to 2013.

Results

The presence of TDRMs increased during this period (from 5.2% to 11.4%), primarily driven by an increase in nonnucleoside reverse-transcriptase (RT) inhibitor (NNRTI) resistance mutations (from 0.3% to 7.1%), particularly K103N/S (increase from 0.3% to 5.3%). Nucleoside/nucleotide RT inhibitor mutations were found in 3.1% of patients. Only 1 patient had K65R (0.01%) and 7 had M184V/I (0.1%), despite high use of tenofovir disoproxil fumarate (TDF), emtricitabine, and lamivudine and potential transmission of resistance to these drugs. At least 1 thymidine-analogue mutations was present in 2.7% of patients with 0.07% harboring T215Y/F and 2.7% harboring T215 revertant mutations (T215rev). Patients with the combination of M41L + L210W + T215rev showed full human immunodeficiency virus RNA suppression while receiving a TDF- or tenofovir alafenamide-containing regimen.

Conclusions

There was an overall increase of TDRMs among patients enrolling in clinical trials from 2000 through 2013, driven primarily by an increase in NNRTI resistance. However, the presence of common TDRMs, including thymidine-analogue mutations/T215rev, showed no impact on response to TDF- or tenofovir alafenamide-containing regimens.

Defining HIV-1 transmission clusters based on sequence data

Understanding HIV-1 transmission dynamics is relevant to both screening and intervention strategies of HIV-1 infection. Commonly, HIV-1 transmission chains are determined based on sequence similarity assessed either directly from a sequence alignment or by inferring a phylogenetic tree. This review is aimed at both nonexperts interested in understanding and interpreting studies of HIV-1 transmission, and experts interested in finding the most appropriate cluster definition for a specific dataset and research question. We start by introducing the concepts and methodologies of how HIV-1 transmission clusters usually have been defined. We then present the results of a systematic review of 105 HIV-1 molecular epidemiology studies summarizing the most common methods and definitions in the literature. Finally, we offer our perspectives on how HIV-1 transmission clusters can be defined and provide some guidance based on examples from real life datasets.

An Efficient Microarray-Based Genotyping Platform for the Identification of Drug-Resistance Mutations in Majority and Minority Subpopulations of HIV-1 Quasispecies

The response of human immunodeficiency virus type 1 (HIV-1) quasispecies to antiretroviral therapy is influenced by the ensemble of mutants that composes the evolving population. Low-abundance subpopulations within HIV-1 quasispecies may determine the viral response to the administered drug combinations. However, routine sequencing assays available to clinical laboratories do not recognize HIV-1 minority variants representing less than 25% of the population. Although several alternative and more sensitive genotyping techniques have been developed, including next-generation sequencing (NGS) methods, they are usually very time consuming, expensive and require highly trained personnel, thus becoming unrealistic approaches in daily clinical practice. Here we describe the development and testing of a HIV-1 genotyping DNA microarray that detects and quantifies, in majority and minority viral subpopulations, relevant mutations and amino acid insertions in 42 codons of the pol gene associated with drug- and multidrug-resistance to protease (PR) and reverse transcriptase (RT) inhibitors. A customized bioinformatics protocol has been implemented to analyze the microarray hybridization data by including a new normalization procedure and a stepwise filtering algorithm, which resulted in the highly accurate (96.33%) detection of positive/negative signals. This microarray has been tested with 57 subtype B HIV-1 clinical samples extracted from multi-treated patients, showing an overall identification of 95.53% and 89.24% of the queried PR and RT codons, respectively, and enough sensitivity to detect minority subpopulations representing as low as 5–10% of the total quasispecies. The developed genotyping platform represents an efficient diagnostic and prognostic tool useful to personalize antiviral treatments in clinical practice.

Deep sequencing for HIV-1 clinical management

The emerging HIV-1 resistance epidemic is threatening the impressive global advances in HIV-1 infection treatment and prevention achieved in the last decade. Next-generation sequencing is improving our ability to understand, diagnose and prevent HIV-1 resistance, being increasingly cost-effective and more accessible. However, NGS still faces a number of limitations that need to be addressed to enable its widespread use. Here, we will review the main NGS platforms available for HIV-1 diagnosis, the factors affecting the clinical utility of NGS testing and the evidence supporting -or not- ultrasensitive genotyping over Sanger sequencing for routine HIV-1 diagnosis. Now that global HIV-1 eradication might be within our reach, making NGS accessible also to LMICs has become a priority. Reductions in sequencing costs, particularly in library preparation, and accessibility to low-cost, robust but simplified automated bioinformatic analyses of NGS data will remain essential to end the HIV-1 pandemic.

Enhanced surveillance of HIV-1 drug resistance in recently infected MSM in the UK

OBJECTIVES

To determine the prevalence of inferred low-frequency HIV-1 transmitted drug resistance (TDR) in MSM in the UK and its predicted effect on first-line therapy.

METHODS

The HIV-1 pol gene was amplified from 442 newly diagnosed MSM identified as likely recently infected by serological avidity testing in 2011-13. The PCR products were sequenced by next-generation sequencing with a mutation frequency threshold of >2% and TDR mutations defined according to the 2009 WHO surveillance drug resistance mutations list.

RESULTS

The majority (75.6%) were infected with subtype B and 6.6% with rare complex or unique recombinant forms. At a mutation frequency threshold of >20%, 7.2% (95% CI 5.0%-10.1%) of the sequences had TDR and this doubled to 15.8% (95% CI 12.6%-19.6%) at >2% mutation frequency (P < 0.0001). The majority (26/42, 62%) of low-frequency variants were against PIs. The most common mutations detected at >20% and 2%-20% mutation frequency differed for each drug class, these respectively being: L90M (n = 7) and M46IL (n = 10) for PIs; T215rev (n = 9) and D67GN (n = 4) for NRTIs; and K103N (n = 5) and G190E (n = 2) for NNRTIs. Combined TDR was more frequent in subtype B than non-B (OR = 0.38; 95% CI = 0.17-0.88; P = 0.024) and had minimal predicted effect on recommended first-line therapies.

CONCLUSIONS

The data suggest differences in the types of low-frequency compared with majority TDR variants that require a better understanding of the origins and clinical significance of low-frequency variants. This will better inform diagnostic and treatment strategies.

Prevalence of integrase inhibitor resistance mutations in Austrian patients recently diagnosed with HIV from 2008 to 2013

PURPOSE

Treatment guidelines often do not advocate testing for integrase inhibitor resistance associated mutations (IRAM) before initiation of first line ART given the extremely low prevalence of mutations found in older surveillance studies. We aimed to describe the prevalence of IRAM in Austrian patients recently diagnosed with HIV in the 5 years following introduction of integrase inhibitors and to analyse trends and factors associated with their detection.

METHODS

Samples of antiretroviral treatment (ART) naïve patients recently diagnosed with HIV in Austria between 2008 and 2013 were analysed for the existence of IRAM and drug penalty scores were calculated to estimate response to drugs. Demographic and virological data were extracted from a database. Descriptive and comparative statistics were used.

RESULTS

A total of 303 samples were analysed. 78 % were male and mean age was 38 years. Overall prevalence of IRAM was 2.3 %. Six percent had at least potentially low-level resistance to raltegravir or elvitegravir, versus 1 % for dolutegravir. One primary mutation was observed (F121Y) in a patient sample from 2012 leading to 5-10-fold reduced susceptibility to raltegravir and elvitegravir. Two patients carried the accessory mutations E138K and G140A, respectively, where both lie on the Q148 pathway. No temporal trend of IRAM prevalence was observed (p = 0.16).

DISCUSSION

Primary IRAM are still rarely found despite the increasing use of INSTI in Austria, but there is a potential for reduced susceptibility to these drugs in selected patients. Routine resistance testing seems prudent to avoid the consequences including accumulation of further mutations and therapeutic failure.

Preexposure prophylaxis-selected drug resistance decays rapidly after drug cessation

OBJECTIVE

Resistance to emtricitabine plus tenofovir disoproxil fumarate (FTC/TDF) or TDF alone used as preexposure prophylaxis (PrEP) has been detected in individuals who initiated PrEP during unrecognized acute HIV infection and, rarely, in PrEP breakthrough infections. PrEP-selected resistance could alter future treatment options, and therefore we sought to determine how long resistance persisted after PrEP cessation.

METHODS

The Partners PrEP Study was a randomized placebo-controlled trial of FTC/TDF or TDF as PrEP for HIV prevention. We previously reported that PrEP-related mutations (K65R, K70E or M184IV) were detected by 454 sequencing following seroconversion in nine individuals who acquired HIV during the Partners PrEP Study. In the current study, we used 454 sequencing to detect and quantify PrEP-related mutations in HIV RNA-positive plasma samples prior to seroconversion, as well as in plasma from 6, 12, and 24 months after PrEP cessation from these nine individuals.

RESULTS

HIV RNA-positive, antibody-negative samples were available prior to seroconversion for four of nine individuals with resistance detected at seroconversion. In all four cases, K65R, K70E and M184IV were not detected prior to seroconversion, suggesting PrEP-related resistance was selected and not transmitted. All PrEP-selected mutations were no longer detectable by 6 months after PrEP cessation and remained undetectable at 12 and 24 months in the absence of antiretroviral therapy.

CONCLUSION

Using highly sensitive assays, PrEP-selected resistance in plasma decays below detection by 6 months following drug cessation and remains undetectable for at least 24 months. Even high levels of resistance mutations during acute infection decay rapidly in the absence of ongoing PrEP exposure.

Contribution of APOBEC3G/F activity to the development of low-abundance drug-resistant human immunodeficiency virus type 1 variants

Plasma drug-resistant minority human immunodeficiency virus type 1 variants (DRMVs) increase the risk of virological failure to first-line non-nucleoside reverse transcriptase inhibitor antiretroviral therapy (ART). The origin of DRMVs in ART-naive patients, however, remains unclear. In a large pan-European case-control study investigating the clinical relevance of pre-existing DRMVs using 454 pyrosequencing, the six most prevalent plasma DRMVs detected corresponded to G-to-A nucleotide mutations (V90I, V106I, V108I, E138K, M184I and M230I). Here, we evaluated if such DRMVs could have emerged from apolipoprotein B mRNA editing enzyme, catalytic polypeptide 3G/F (APOBEC3G/F) activity. Out of 236 ART-naive subjects evaluated, APOBEC3G/F hypermutation signatures were detected in plasma viruses of 14 (5.9%) individuals. Samples with minority E138K, M184I, and M230I mutations, but not those with V90I, V106I or V108I, were significantly associated with APOBEC3G/F activity (Fisher's P < 0.005), defined as the presence of > 0.5% of sample sequences with an APOBEC3G/F signature. Mutations E138K, M184I and M230I co-occurred in the same sequence as APOBEC3G/F signatures in 3/9 (33%), 5/11 (45%) and 4/8 (50%) of samples, respectively; such linkage was not found for V90I, V106I or V108I. In-frame STOP codons were observed in 1.5% of all clonal sequences; 14.8% of them co-occurred with APOBEC3G/F signatures. APOBEC3G/F-associated E138K, M184I and M230I appeared within clonal sequences containing in-frame STOP codons in 2/3 (66%), 5/5 (100%) and 4/4 (100%) of the samples. In a re-analysis of the parent case control study, the presence of APOBEC3G/F signatures was not associated with virological failure. In conclusion, the contribution of APOBEC3G/F editing to the development of DRMVs is very limited and does not affect the efficacy of non-nucleoside reverse transcriptase inhibitor ART.

Longitudinal Detection and Persistence of Minority Drug-Resistant Populations and Their Effect on Salvage Therapy

BACKGROUND

Drug-resistant HIV are more prevalent and persist longer than previously demonstrated by bulk sequencing due to the ability to detect low-frequency variants. To clarify a clinical benefit to monitoring minority-level drug resistance populations as a guide to select active drugs for salvage therapy, we retrospectively analyzed the dynamics of low-frequency drug-resistant population in antiretroviral (ARV)-exposed drug resistant individuals.

MATERIALS AND METHODS

Six HIV-infected individuals treated with ARV for more than five years were analyzed. These individuals had difficulty in controlling viremia, and treatment regimens were switched multiple times guided by standard drug resistance testing using bulk sequencing. To detect minority variant populations with drug resistance, we used a highly sensitive allele-specific PCR (AS-PCR) with detection thresholds of 0.3-2%. According to ARV used in these individuals, we focused on the following seven reverse transcriptase inhibitor-resistant mutations: M41L, K65R, K70R, K103N, Y181C, M184V, and T215F/Y. Results of AS-PCR were compared with bulk sequencing data for concordance and presence of additional mutations. To clarify the genetic relationship between low-frequency and high-frequency populations, AS-PCR amplicon sequences were compared with bulk sequences in phylogenetic analysis.

RESULTS

The use of AS-PCR enabled detection of the drug-resistant mutations, M41L, K103N, Y181C, M184V and T215Y, present as low-frequency populations in five of the six individuals. These drug resistant variants persisted for several years without ARV pressure. Phylogenetic analysis indicated that pre-existing K103N and T215I variants had close genetic relationships with high-frequency K103N and T215I observed during treatment.

DISCUSSION AND CONCLUSION

Our results demonstrate the long-term persistence of drug-resistant viruses in the absence of drug pressure. The rapid virologic failures with pre-existing mutant viruses detectable by AS-PCR highlight the clinical importance of low-frequency drug-resistant viruses. Thus, our results highlight the usefulness of AS-PCR and support its expanded evaluation in ART clinical management.

The impact of vaccination on the breadth and magnitude of the antibody response to influenza A viruses in HIV-infected individuals

OBJECTIVE

HIV-positive individuals have lower antibody titers to influenza viruses than HIV-negative individuals, and the benefits of the annual vaccinations are controversially discussed. Also, there is no information about the breadth of the antibody response in HIV-infected individuals.

DESIGN

The binding and neutralizing antibody titers to various human and nonhuman influenza A virus strain were determined in sera from 146 HIV-infected volunteers: They were compared with those found in 305 randomly selected HIV-negative donors, and put in relation to HIV-specific parameters. Univariable and multivariable regression was used to identify HIV-specific parameters associated with the measured binding and neutralizing activity.

METHODS

Enzyme-linked immunosorbent assays and in-vitro neutralization assays were used to determine the binding and neutralizing antibodiy titers to homo and heterosubtypic influenza A subtypes.

RESULTS

We found that both homo and heterosubtypic antibody titers are lower in HIV-positive individuals. Vaccination promoted higher binding and neutralizing antibody titers to human but not to nonhuman isolates. HIV-induced immune damage (high viral load, low CD4 T-cell counts, and long untreated disease progression) is associated with impaired homosubtypic responses, but can have beneficial effects on the development of heterosubtypic antibodies, and an improved ratio of binding to neutralizing antibody titers to homosubtypic isolates.

CONCLUSIONS

Our results indicate that repetitive vaccinations in HIV-positive individuals enhance antibody titers to human isolates. Interestingly, development of antibody titers to conserved heterosubtypic epitopes paradoxically appeared to profit from HIV-induced immune damage, as did the ratio of binding to neutralizing antibodies.

Trends in Transmission of Drug Resistance and Prevalence of Non-B Subtypes in Patients with Acute or Recent HIV-1 Infection in Barcelona in the Last 16 Years (1997-2012)

Objectives

To evaluate the prevalence of transmitted drug resistance (TDR) and non-B subtypes in patients with acute/recent HIV-1 infection in Barcelona during the period 1997-2012.

Methods

Patients from the “Hospital Clínic Primary HIV-1 Infection Cohort” with a genotyping test performed within 180 days of infection were included. The 2009 WHO List of Mutations for Surveillance of Transmitted HIV-1 Drug Resistance was used for estimating the prevalence of TDR and phylogenetic analysis for subtype determination.

Results

189 patients with acute/recent HIV-1 infection were analyzed in 4 time periods (1997-2000, n=28; 2001-4, n=42; 2005-8, n=55 and 2009-12, n=64). The proportion of patients with acute/recent HIV-1 infection with respect to the total of newly HIV-diagnosed patients in our center increased over the time and was 2.18%, 3.82%, 4.15% and 4.55% for the 4 periods, respectively (p=0.005). The global prevalence of TDR was 9%, or 17.9%, 9.5%, 3.6% and 9.4% by study period (p=0.2). The increase in the last period was driven by protease-inhibitor and nucleoside-reverse-transcriptase-inhibitor resistance mutations while non-nucleoside-reverse-transcriptase inhibitor TDR and TDR of more than one family decreased. The overall prevalence of non-B subtypes was 11.1%, or 0%, 4.8%, 9.1% and 20.3 by study period (p=0.01). B/F recombinants, B/G recombinants and subtype F emerged in the last period. We also noticed an increase in the number of immigrant patients (p=0.052). The proportion of men-who-have-sex-with-men (MSM) among patients with acute/recent HIV-1 infection increased over the time (p=0.04).

Conclusions

The overall prevalence of TDR in patients with acute/recent HIV-1 infection in Barcelona was 9%, and it has stayed relatively stable in recent years. Non-B subtypes and immigrants proportions progressively increased.

Epitope-specific CD8+ T cell kinetics rather than viral variability determine the timing of immune escape in simian immunodeficiency virus infection

CD8(+) T cells are important for the control of chronic HIV infection. However, the virus rapidly acquires "escape mutations" that reduce CD8(+) T cell recognition and viral control. The timing of when immune escape occurs at a given epitope varies widely among patients and also among different epitopes within a patient. The strength of the CD8(+) T cell response, as well as mutation rates, patterns of particular amino acids undergoing escape, and growth rates of escape mutants, may affect when escape occurs. In this study, we analyze the epitope-specific CD8(+) T cells in 25 SIV-infected pigtail macaques responding to three SIV epitopes. Two epitopes showed a variable escape pattern and one had a highly monomorphic escape pattern. Despite very different patterns, immune escape occurs with a similar delay of on average 18 d after the epitope-specific CD8(+) T cells reach 0.5% of total CD8(+) T cells. We find that the most delayed escape occurs in one of the highly variable epitopes, and that this is associated with a delay in the epitope-specific CD8(+) T cells responding to this epitope. When we analyzed the kinetics of immune escape, we found that multiple escape mutants emerge simultaneously during the escape, implying that a diverse population of potential escape mutants is present during immune selection. Our results suggest that the conservation or variability of an epitope does not appear to affect the timing of immune escape in SIV. Instead, timing of escape is largely determined by the kinetics of epitope-specific CD8(+) T cells.

Persistence of transmitted HIV-1 drug resistance mutations associated with fitness costs and viral genetic backgrounds

Transmission of drug-resistant pathogens presents an almost-universal challenge for fighting infectious diseases. Transmitted drug resistance mutations (TDRM) can persist in the absence of drugs for considerable time. It is generally believed that differential TDRM-persistence is caused, at least partially, by variations in TDRM-fitness-costs. However, in vivo epidemiological evidence for the impact of fitness costs on TDRM-persistence is rare.

Here, we studied the persistence of TDRM in HIV-1 using longitudinally-sampled nucleotide sequences from the Swiss-HIV-Cohort-Study (SHCS). All treatment-naïve individuals with TDRM at baseline were included. Persistence of TDRM was quantified via reversion rates (RR) determined with interval-censored survival models. Fitness costs of TDRM were estimated in the genetic background in which they occurred using a previously published and validated machine-learning algorithm (based on in vitro replicative capacities) and were included in the survival models as explanatory variables.

In 857 sequential samples from 168 treatment-naïve patients, 17 TDRM were analyzed. RR varied substantially and ranged from 174.0/100-person-years;CI=[51.4, 588.8] (for 184V) to 2.7/100-person-years;[0.7, 10.9] (for 215D). RR increased significantly with fitness cost (increase by 1.6[1.3,2.0] per standard deviation of fitness costs). When subdividing fitness costs into the average fitness cost of a given mutation and the deviation from the average fitness cost of a mutation in a given genetic background, we found that both components were significantly associated with reversion-rates.

Our results show that the substantial variations of TDRM persistence in the absence of drugs are associated with fitness-cost differences both among mutations and among different genetic backgrounds for the same mutation.

The evolution of resistance is a universal challenge in antimicrobial chemotherapy. A key driver of resistance is that drug resistance mutations often persist even in the absence of drugs and despite the fact that resistance mutations are often associated with reduced pathogen replication (“fitness costs”). Such persistence may occur because fitness costs are low, especially if they are compensated by additional mutations in their “genetic background”. Here we assessed the role of fitness-cost and the genetic background for resistance in a real-world epidemiological setting by studying the persistence behavior of transmitted antiretroviral resistance mutations of HIV. This persistence behavior was associated with the predicted fitness cost of a given resistance mutation in the particular genetic background in which it occurred. We found that persistence behavior varied strongly across both mutation types and genetic backgrounds and that persistence was significantly associated with predicted fitness costs. In particular we found that even mutations of the same type tended to persist longer if they occurred in a genetic background where they caused weak fitness costs. Overall our results underline the variability of persistence behavior as well as the important role of fitness costs and the genetic background in the evolution of antimicrobial resistance.

Assessing the Paradox Between Transmitted and Acquired HIV Type 1 Drug Resistance Mutations in the Swiss HIV Cohort Study From 1998 to 2012

BACKGROUND

Transmitted human immunodeficiency virus type 1 (HIV) drug resistance (TDR) mutations are transmitted from nonresponding patients (defined as patients with no initial response to treatment and those with an initial response for whom treatment later failed) or from patients who are naive to treatment. Although the prevalence of drug resistance in patients who are not responding to treatment has declined in developed countries, the prevalence of TDR mutations has not. Mechanisms causing this paradox are poorly explored.

METHODS

We included recently infected, treatment-naive patients with genotypic resistance tests performed ≤ 1 year after infection and before 2013. Potential risk factors for TDR mutations were analyzed using logistic regression. The association between the prevalence of TDR mutations and population viral load (PVL) among treated patients during 1997-2011 was estimated with Poisson regression for all TDR mutations and individually for the most frequent resistance mutations against each drug class (ie, M184V/L90M/K103N).

RESULTS

We included 2421 recently infected, treatment-naive patients and 5399 patients with no response to treatment. The prevalence of TDR mutations fluctuated considerably over time. Two opposing developments could explain these fluctuations: generally continuous increases in the prevalence of TDR mutations (odds ratio, 1.13; P = .010), punctuated by sharp decreases in the prevalence when new drug classes were introduced. Overall, the prevalence of TDR mutations increased with decreasing PVL (rate ratio [RR], 0.91 per 1000 decrease in PVL; P = .033). Additionally, we observed that the transmitted high-fitness-cost mutation M184V was positively associated with the PVL of nonresponding patients carrying M184V (RR, 1.50 per 100 increase in PVL; P < .001). Such association was absent for K103N (RR, 1.00 per 100 increase in PVL; P = .99) and negative for L90M (RR, 0.75 per 100 increase in PVL; P = .022).

CONCLUSIONS

Transmission of antiretroviral drug resistance is temporarily reduced by the introduction of new drug classes and driven by nonresponding and treatment-naive patients. These findings suggest a continuous need for new drugs, early detection/treatment of HIV-1 infection.

Diminished transmission of drug resistant HIV-1 variants with reduced replication capacity in a human transmission model

Background

Different patterns of drug resistance are observed in treated and therapy naïve HIV-1 infected populations. Especially the NRTI-related M184I/V variants, which are among the most frequently encountered mutations in treated patients, are underrepresented in the antiretroviral naïve population. M184I/V mutations are known to have a profound effect on viral replication and tend to revert over time in the new host. However it is debated whether a diminished transmission efficacy of HIV variants with a reduced replication capacity can also contribute to the observed discrepancy in genotypic patterns.

As dendritic cells (DCs) play a pivotal role in HIV-1 transmission, we used a model containing primary human Langerhans cells (LCs) and DCs to compare the transmission efficacy M184 variants (HIV-M184V/I/T) to HIV wild type (HIV-WT). As control, we used HIV harboring the NNRTI mutation K103N (HIV-K103N) which has a minor effect on replication and is found at a similar prevalence in treated and untreated individuals.

Results

In comparison to HIV-WT, the HIV-M184 variants were less efficiently transmitted to CCR5⁺ Jurkat T cells by both LCs and DCs. The transmission rate of HIV-K103N was slightly reduced to HIV-WT in LCs and even higher than HIV-WT in DCs. Replication experiments in CCR5⁺ Jurkat T cells revealed no apparent differences in replication capacity between the mutant viruses and HIV-WT. However, viral replication in LCs and DCs was in concordance with the transmission results; replication by the HIV-M184 variants was lower than replication by HIV-WT, and the level of replication of HIV-K103N was intermediate for LCs and higher than HIV-WT for DCs.

Conclusions

Our data demonstrate that drug resistant M184-variants display a reduced replication capacity in LCs and DCs which directly impairs their transmission efficacy. As such, diminished transmission efficacy may contribute to the lower prevalence of drug resistant variants in therapy naive individuals.

Low-frequency drug-resistant HIV-1 and risk of virological failure to first-line NNRTI-based ART: a multicohort European case-control study using centralized ultrasensitive 454 pyrosequencing

Objectives

It is still debated if pre-existing minority drug-resistant HIV-1 variants (MVs) affect the virological outcomes of first-line NNRTI-containing ART.

Methods

This Europe-wide case–control study included ART-naive subjects infected with drug-susceptible HIV-1 as revealed by population sequencing, who achieved virological suppression on first-line ART including one NNRTI. Cases experienced virological failure and controls were subjects from the same cohort whose viraemia remained suppressed at a matched time since initiation of ART. Blinded, centralized 454 pyrosequencing with parallel bioinformatic analysis in two laboratories was used to identify MVs in the 1%–25% frequency range. ORs of virological failure according to MV detection were estimated by logistic regression.

Results

Two hundred and sixty samples (76 cases and 184 controls), mostly subtype B (73.5%), were used for the analysis. Identical MVs were detected in the two laboratories. 31.6% of cases and 16.8% of controls harboured pre-existing MVs. Detection of at least one MV versus no MVs was associated with an increased risk of virological failure (OR = 2.75, 95% CI = 1.35–5.60, P = 0.005); similar associations were observed for at least one MV versus no NRTI MVs (OR = 2.27, 95% CI = 0.76–6.77, P = 0.140) and at least one MV versus no NNRTI MVs (OR = 2.41, 95% CI = 1.12–5.18, P = 0.024). A dose–effect relationship between virological failure and mutational load was found.

Conclusions

Pre-existing MVs more than double the risk of virological failure to first-line NNRTI-based ART.

Limited clinical benefit of minority K103N and Y181C-variant detection in addition to routine genotypic resistance testing in antiretroviral therapy-naive patients

OBJECTIVE

The presence of minority nonnucleoside reverse transcriptase inhibitor (NNRTI)-resistant HIV-1 variants prior to antiretroviral therapy (ART) has been linked to virologic failure in treatment-naive patients.

DESIGN

We performed a large retrospective study to determine the number of treatment failures that could have been prevented by implementing minority drug-resistant HIV-1 variant analyses in ART-naïve patients in whom no NNRTI resistance mutations were detected by routine resistance testing.

METHODS

Of 1608 patients in the Swiss HIV Cohort Study, who have initiated first-line ART with two nucleoside reverse transcriptase inhibitors (NRTIs) and one NNRTI before July 2008, 519 patients were eligible by means of HIV-1 subtype, viral load and sample availability. Key NNRTI drug resistance mutations K103N and Y181C were measured by allele-specific PCR in 208 of 519 randomly chosen patients.

RESULTS

Minority K103N and Y181C drug resistance mutations were detected in five out of 190 (2.6%) and 10 out of 201 (5%) patients, respectively. Focusing on 183 patients for whom virologic success or failure could be examined, virologic failure occurred in seven out of 183 (3.8%) patients; minority K103N and/or Y181C variants were present prior to ART initiation in only two of those patients. The NNRTI-containing, first-line ART was effective in 10 patients with preexisting minority NNRTI-resistant HIV-1 variant.

CONCLUSION

As revealed in settings of case-control studies, minority NNRTI-resistant HIV-1 variants can have an impact on ART. However, the implementation of minority NNRTI-resistant HIV-1 variant analysis in addition to genotypic resistance testing (GRT) cannot be recommended in routine clinical settings. Additional associated risk factors need to be discovered.

Full-length haplotype reconstruction to infer the structure of heterogeneous virus populations

Next-generation sequencing (NGS) technologies enable new insights into the diversity of virus populations within their hosts. Diversity estimation is currently restricted to single-nucleotide variants or to local fragments of no more than a few hundred nucleotides defined by the length of sequence reads. To study complex heterogeneous virus populations comprehensively, novel methods are required that allow for complete reconstruction of the individual viral haplotypes. Here, we show that assembly of whole viral genomes of ∼8600 nucleotides length is feasible from mixtures of heterogeneous HIV-1 strains derived from defined combinations of cloned virus strains and from clinical samples of an HIV-1 superinfected individual. Haplotype reconstruction was achieved using optimized experimental protocols and computational methods for amplification, sequencing and assembly. We comparatively assessed the performance of the three NGS platforms 454 Life Sciences/Roche, Illumina and Pacific Biosciences for this task. Our results prove and delineate the feasibility of NGS-based full-length viral haplotype reconstruction and provide new tools for studying evolution and pathogenesis of viruses.

HIV reverse-transcriptase drug resistance mutations during early infection reveal greater transmission diversity than in envelope sequences

BACKGROUND

Drug resistance mutations (DRMs) can serve as distinct, nonpolymorphic markers for evaluating diversity of expressed HIV-1. We screened for DRMs during early-acute viremia and examined the diversity in reverse transcriptase (RT) relative to envelope (env) in cases of transmitted drug resistance.

METHODS

We evaluated 111 longitudinal plasma samples collected every 2-7 days from 15 individuals who seroconverted for HIV-1 infection in 1994-2000. The samples were screened with sensitive polymerase chain reaction assays for the commonly transmitted M41L and K70R mutations and for K65R, which was undetected by bulk sequencing. Mutation-positive samples were further characterized by clonal sequencing of RT and env V1-V3.

RESULTS

Drug resistance mutations were detected in 4 of 15 seroconverters at 5-50 days of viral nucleic acid expression; most mutations disappeared about the time of seroconversion. Clonal sequencing verified low-level K65R at frequencies of 0.4%-4.9%. In each case, K65R coexisted unlinked with variants carrying 2-5 thymidine analog mutations at frequencies of 1.6%-23.0%. In one seroconverter, variants with M184V and nonnucleoside RT inhibitor mutations were also identified at first RNA expression. Each seroconverter displayed a homogeneous V1-V3 env population.

CONCLUSIONS

Reverse-transcriptase DRMs demonstrate that the breadth of variants in transmission may be greater than what is reflected in envelope sequences.

Lack of resistance to integrase inhibitors among antiretroviral-naive subjects with primary HIV-1 infection, 2007-2013

BACKGROUND

US guidelines recommend genotyping for persons newly diagnosed with HIV infection to identify transmitted drug resistance mutations associated with decreased susceptibility to nucleoside reverse transcriptase inhibitors, non-nucleoside reverse transcriptase inhibitors and protease inhibitors. To date, testing for integrase strand transfer inhibitor (INSTI) mutations has not been routinely recommended. We aimed to evaluate the prevalence of transmitted INSTI mutations among persons with primary HIV-1 infection in Seattle, WA, USA.

METHODS

Persons with primary HIV-1 infection have enrolled in an observational cohort at the University of Washington Primary Infection Clinic since 1992. We performed a retrospective analysis of plasma specimens collected prospectively from the 82 antiretroviral-naive subjects who were enrolled from 2007-2013, after FDA-approval of the first INSTI. Resistance testing was performed by consensus sequencing.

RESULTS

Specimens for analysis had been obtained a median of 24 (IQR 18-41, range 8-108) days after the estimated date of HIV-1 infection. All subjects were infected with HIV-1 subtype B except for one subject infected with subtype C. Consensus sequencing identified no subjects with major INSTI mutations (T66I, E92Q, G140S, Y143C/H/R, S147G, Q148H/K/R, N155H). Using exact binomial CIs, the upper bound of the 95% CI was 4.4%.

CONCLUSIONS

Although our sample size was small, this study does not support the need at this time to evaluate integrase mutations as part of routine consensus sequencing among persons newly diagnosed with HIV-1 infection. However, it is likely that the prevalence of transmitted INSTI mutations may increase with the recent commercial introduction of additional INSTIs and presumably greater INSTI use among persons living with HIV-1.

Next-generation sequencing of HIV-1 RNA genomes: determination of error rates and minimizing artificial recombination

Next-generation sequencing (NGS) is a valuable tool for the detection and quantification of HIV-1 variants in vivo. However, these technologies require detailed characterization and control of artificially induced errors to be applicable for accurate haplotype reconstruction. To investigate the occurrence of substitutions, insertions, and deletions at the individual steps of RT-PCR and NGS, 454 pyrosequencing was performed on amplified and non-amplified HIV-1 genomes. Artificial recombination was explored by mixing five different HIV-1 clonal strains (5-virus-mix) and applying different RT-PCR conditions followed by 454 pyrosequencing. Error rates ranged from 0.04-0.66% and were similar in amplified and non-amplified samples. Discrepancies were observed between forward and reverse reads, indicating that most errors were introduced during the pyrosequencing step. Using the 5-virus-mix, non-optimized, standard RT-PCR conditions introduced artificial recombinants in a fraction of at least 30% of the reads that subsequently led to an underestimation of true haplotype frequencies. We minimized the fraction of recombinants down to 0.9-2.6% by optimized, artifact-reducing RT-PCR conditions. This approach enabled correct haplotype reconstruction and frequency estimations consistent with reference data obtained by single genome amplification. RT-PCR conditions are crucial for correct frequency estimation and analysis of haplotypes in heterogeneous virus populations. We developed an RT-PCR procedure to generate NGS data useful for reliable haplotype reconstruction and quantification.