Minority variants associated with resistance to HIV-1 nonnucleoside reverse transcriptase inhibitors during primary infection

Full-spectrum HIV drug resistance mutation detection by high-resolution complete pol gene sequencing

BACKGROUND

Drug resistance mutation testing is a key element for HIV clinical management, informing effective treatment regimens. However, resistance screening in current clinical practice is limited in reporting linked cross-class resistance mutations and minority variants, both of which may increase the risk of virological failure.

METHODS

To address these limitations, we obtained 358 full-length pol gene sequences from 52 specimens of 20 HIV infected individuals by combining microdroplet amplification, unique molecular identifier (UMI) labeling, and long-read high-throughput sequencing.

RESULTS

We conducted a rigorous assessment of the accuracy of our pipeline for precision drug resistance mutation detection, verifying that a sequencing depth of 35 high-throughput reads achieved complete, error-free pol gene sequencing. We detected 26 distinct drug resistance mutations to Protease Inhibitors (PIs), Nucleoside Reverse Transcriptase Inhibitors (NRTIs), Non-Nucleoside Reverse Transcriptase Inhibitors (NNRTIs), and Integrase Strand Transfer Inhibitors (INSTIs). We detected linked cross-class drug resistance mutations (PI+NRTI, PI+NNRTI, and NRTI+NNRTI) that confer cross-resistance to multiple drugs in different classes. Fourteen different types of minority mutations were also detected with frequencies ranging from 3.2% to 19%, and the presence of these mutations was verified by Sanger reference sequencing. We detected a putative transmitted drug resistance mutation (TDRM) in one individual that persisted for over seven months from the first sample collected at the acute stage of infection prior to seroconversion.

CONCLUSIONS

Our comprehensive drug resistance mutation profiling can advance clinical practice by reporting mutation linkage and minority variants to better guide antiretroviral therapy options.

Establishment and application of a method of tagged-amplicon deep sequencing for low-abundance drug resistance in HIV-1

In the latest HIV-1 global drug resistance report released by WHO, countries are advised to strengthen pre-treatment monitoring of drug resistance in AIDS patients. In this study, we established an NGS-based segmented amplification HIV-1 drug resistance mutation detection method. The pol region of HIV-1 was divided into three short fragments for NGS. The entire amplification and sequencing panel were more cost-effective and batched by using the barcode sequence corresponding to the sample. Each parameter was evaluated using samples with known resistance variants frequencies. The nucleotide sequence error rate, amino acid error rate, and noise value of the NGS-based segmented amplification method were both less than 1%. When the threshold was 2%, the consensus sequences of the HIV-1 NL4-3 strain were completely consistent with the Sanger sequences. This method can detect the minimum viral load of the sample at 10² copies/ml, and the input frequency and detection frequency of HIV-1 resistance mutations within the range of 1%–100% had good conformity (R² = 0.9963; R² = 0.9955). This method had no non-specific amplification for Hepatitis B and C. Under the 2% threshold, the incidence of surveillance drug resistance mutations in ART-naive HIV-infected patients was 20.69%, among which NRTIs class resistance mutations were mainly.

Diversity of HIV-1 subtypes and transmitted drug-resistance mutations among minority HIV-1 variants in a Turkish cohort

BACKGROUND

The World Health Organization (WHO) recommends the surveillance of transmitted drug resistance mutations (TDRMs) to ensure the effectiveness and sustainability of HIV treatment programs.

OBJECTIVE

Our aim was to determine the TDRMs and evaluate the distribution of HIV-1 subtypes using and compared next-generation sequencing (NGS) and Sanger-based sequencing (SBS) in a cohort of 44 antiretroviral treatment-naïve patients.

METHODS

All samples that were referred to the microbiology laboratory for HIV drug resistance analysis between December 2016 and February 2018 were included in the study. After exclusions, 44 treatment-naive adult patients with a viral load of >1000 copies/mL were analyzed. DNA sequencing for reverse transcriptase and protease regions was performed using both DeepChek ABL single round kit and Sanger-based ViroSeq HIV-1 Genotyping System. The mutations and HIV-1 subtypes were analyzed using the Stanford HIVdb version 8.6.1 Genotypic Resistance software, and TDRMs were assessed using the WHO surveillance drug-resistance mutation database. HIV-1 subtypes were confirmed by constructing a maximum-likelihood phylogenetic tree using Los Alamos IQ-Tree software.

RESULTS

NGS identified nucleos(t)ide reverse transcriptase inhibitor (NRTI)-TDRMs in 9.1% of the patients, non-nucleos(t)ide reverse transcriptase inhibitor (NNRTI)-TDRMs in 6.8% of the patients, and protease inhibitor (PI)-TDRMs in 18.2% of the patients at a detection threshold of ≥1%. Using SBS, 2.3% and 6.8% of the patients were found to have NRTI- and NNRTI-TDRMs, respectively, but no major PI mutations were detected. M41L, L74I, K65R, M184V, and M184I related to NRTI, K103N to NNRTI, and N83D, M46I, I84V, V82A, L24I, L90M, I54V to the PI sites were identified using NGS. Most mutations were found in low-abundance (frequency range: 1.0% - 4.7%) HIV-1 variants, except M41L and K103N. The subtypes of the isolates were found as follows; 61.4% subtype B, 18.2% subtype B/CRF02_AG recombinant, 13.6% subtype A, 4.5% CRF43_02G, and 2.3% CRF02_AG. All TDRMs, except K65R, were detected in HIV-1 subtype B isolates.

CONCLUSION

The high diversity of protease site TDRMs in the minority HIV-1 variants and prevalence of CRFs were remarkable in this study. All minority HIV-1 variants were missed by conventional sequencing. TDRM prevalence among minority variants appears to be decreasing over time at our center.

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.

Impact of Human Immunodeficiency Virus Type 1 Minority Variants on the Virus Response to a Rilpivirine-Based First-line Regimen

Background

Minority resistant variants of human immunodeficiency virus type 1 (HIV-1) could influence the virological response to treatment based on nonnucleoside reverse transcriptase inhibitors (NNRTIs). Data on minority rilpivirine-resistant variants are scarce. This study used next-generation sequencing (NGS) to identify patients harboring minority resistant variants to nucleos(t)ide reverse transcriptase inhibitors and NNRTIs and to assess their influence on the virological response (VR).

Methods

All the subjects, 541 HIV-1-infected patients started a first-line regimen containing rilpivirine. VR was defined as a HIV-1 RNA load <50 copies/mL at month 6 with continued suppression at month 12. NGS was performed at baseline (retrospectively) on the 454 GS-FLX platform (Roche).

Results

NGS revealed resistance-associated mutations accounting for 1% to <5% of variants in 17.2% of samples, for 5%-20% in 5.7% of samples, and for >20% in 29% of samples. We identified 43 (8.8%) and 36 (7.4%) patients who harbored rilpivirine-resistant variants with a 1% sensitivity threshold according to the French National Agency for Research on AIDS and Viral Hepatitis and Stanford algorithms, respectively. The VR was 96.9% at month 12. Detection of minority rilpivirine resistant variants was not associated with virological failure (VF). Multivariate analysis indicated that VF at month 12 was associated with a CD4 count <250 cells/µL at baseline, a slower decrease in viral load at month 3, and rilpivirine resistance at baseline using the Stanford algorithm with a 20% threshold.

Conclusions

Minority resistant variants had no impact on the VR of treatment-naive patients to a rilpivirine-based regimen.

Insights about minority HIV-1 strains in transmitted drug resistance mutation dynamics and disease progression

Objectives

The presence of minority transmitted drug resistance mutations was assessed using ultra-deep sequencing and correlated with disease progression among recently HIV-1-infected individuals from Brazil.

Methods

Samples at baseline during recent infection and 1 year after the establishment of the infection were analysed. Viral RNA and proviral DNA from 25 individuals were subjected to ultra-deep sequencing of the reverse transcriptase and protease regions of HIV-1.

Results

Viral strains carrying transmitted drug resistance mutations were detected in 9 out of the 25 patients, for all major antiretroviral classes, ranging from one to five mutations per patient. Ultra-deep sequencing detected strains with frequencies as low as 1.6% and only strains with frequencies >20% were detected by population plasma sequencing (three patients). Transmitted drug resistance strains with frequencies <14.8% did not persist upon established infection. The presence of transmitted drug resistance mutations was negatively correlated with the viral load and with CD4+ T cell count decay.

Conclusions

Transmitted drug resistance mutations representing small percentages of the viral population do not persist during infection because they are negatively selected in the first year after HIV-1 seroconversion.

HIV-1 genotypic resistance testing using the Vela automated next-generation sequencing platform

Objectives

To evaluate the diagnostic performance of the Vela next-generation sequencing (NGS) system in conjunction with the Sentosa SQ HIV Genotyping Assay for genotyping HIV-1.

Methods

Plasma RNA was extracted and templates prepared with the Sentosa SX instrument before sequencing the HIV-1 polymerase on the Sentosa SQ301 Sequencer (PGM IonTorrent). The Vela NGS System was compared with direct sequencing and the 454 GS-FLX (Roche) and MiSeq (Illumina) systems for genotypic resistance testing on clinical samples.

Results

The Vela NGS system detected majority resistance mutations in subtype B and CRF02-AG samples at 500 copies/mL and minority variants with a sensitivity of 5% at 100 000 copies/mL. The Vela NGS system and direct sequencing identified resistance mutations with 97% concordance in 46 clinical samples. Vela identified 1/20 of the 1%-5% mutations identified by 454, 5/12 of the 5%-20% mutations and 60/61 of the >20% mutations. Vela identified 3/14 of the 1%-5% mutations identified by MiSeq, 0/2 of the 5%-20% mutations and 47/47 of the >20% mutations. The resistance mutation quantifications by Vela and 454 were concordant (bias: 2.31%), as were those by Vela and MiSeq (bias: 1.06%).

Conclusions

The Vela NGS system provides automated nucleic acid extraction, PCR reagent distribution, library preparation and bioinformatics analysis. The analytical performance was very good when compared with direct sequencing, but was less sensitive than two other NGS platforms for detecting minority variants.

Assessment of Transmitted HIV-1 Drug Resistance Mutations Using Ultra- Deep Pyrosequencing in a Turkish Cohort

BACKGROUND

Antiretroviral treatment (ART) reduces morbidity and mortality caused by human immunodeficiency virus (HIV) infection; however, the emergence of drug-resistant strains poses an important obstacle to treatment success. Using conventional sequencing methods to determine antiretroviral resistance, mutations present in ≥20% of quasispecies can be identified, but drug-resistant minority variants can lead to virologic failure.

OBJECTIVE

We aimed to assess transmitted drug resistance mutations (TDRMs) within minority variants using ultra-deep pyrosequencing (UDPS).

METHOD

Treatment-naive adult patients were included in this observational study. Surveillance TDRMs were classified as ≥20% or at minority variant level (≥2% - <20%). Genotypic sensitivity score calculated by using all pre-treatment drug resistance mutations (PDRMs) was also evaluated.

RESULTS

Thirty-six patients were analyzed. Any TDRM at ≥20% level was detected in 8.3% of the patients (n=3). This prevalence increased to 30.6% (n=11) with the inclusion of minority variants. All non-nucleoside reverse transcriptase inhibitor and protease inhibitor-related TDRMs were within minority variants. The genotypic sensitivity score of rilpivirine-based regimens was considerably diminished when minority variants were included in the PDRM analysis.

CONCLUSION

UDPS was used for the first time to assess TDRM in a Turkish HIV cohort and uncovered several mutations hidden within minority variants. UDPS may be preferred to detect PDRMs for avoiding virologic failure with rilpivirine-based ART regimens.

Limited Marginal Utility of Deep Sequencing for HIV Drug Resistance Testing in the Age of Integrase Inhibitors

HIV drug resistance genotyping is a critical tool in the clinical management of HIV infections. Although resistance genotyping has traditionally been conducted using Sanger sequencing, next-generation sequencing (NGS) is emerging as a powerful tool due to its ability to detect low-frequency alleles.

HIV drug resistance genotyping is a critical tool in the clinical management of HIV infections. Although resistance genotyping has traditionally been conducted using Sanger sequencing, next-generation sequencing (NGS) is emerging as a powerful tool due to its ability to detect low-frequency alleles. However, the clinical value added from NGS approaches to antiviral resistance testing remains to be demonstrated. We compared the variant detection capacity of NGS versus Sanger sequencing methods for resistance genotyping in 144 drug resistance tests (105 protease-reverse transcriptase tests and 39 integrase tests) submitted to our clinical virology laboratory over a four-month period in 2016 for Sanger-based HIV drug resistance testing. NGS detected all true high-frequency drug resistance mutations (>20% frequency) found by Sanger sequencing, with greater accuracy in one instance of a Sanger-detected false positive. Freely available online NGS variant callers HyDRA and PASeq were superior to Sanger methods for interpretations of allele linkage and automated variant calling. NGS additionally detected low-frequency mutations (1 to 20% frequency) associated with higher levels of drug resistance in 30/105 (29%) protease-reverse transcriptase tests and 4/39 (10%) integrase tests. In clinical follow-up of 69 individuals for a median of 674 days, we did not find a difference in rates of virological failure between individuals with and without low-frequency mutations, although rates of virological failure were higher for individuals with drug-relevant low-frequency mutations. However, all 27 individuals who experienced virological failure reported poor adherence to their drug regimen during the preceding follow-up time, and all 19 who subsequently improved their adherence achieved viral suppression at later time points, consistent with a lack of clinical resistance. In conclusion, in a population with low antiviral resistance emergence, NGS methods detected numerous instances of minor alleles that did not result in subsequent bona fide virological failure due to antiviral resistance.

RNA and DNA Sanger sequencing versus next-generation sequencing for HIV-1 drug resistance testing in treatment-naive patients

Background

Sanger sequencing of plasma RNA is the standard method for HIV-1 drug resistance testing in treatment-naive patients, but is limited by the non-detection of resistance-associated mutations (RAMs) with prevalence below approximately 20%.

Objectives

We compared RNA and DNA Sanger sequencing (RSS and DSS) with RNA next-generation sequencing (NGS) for RAM detection in HIV-1 reverse transcriptase (RT), protease (PR) and integrase (IN) genes.

Methods

Sanger sequencing was performed on RNA and DNA, following the recommendations of the French Agency for AIDS Research (ANRS). NGS was performed on RNA using the HIV-1 Drug Resistance Assay, v. 3.0 (Roche) on the 454 GS Junior sequencer. The IAS-USA list was used to identify RAMs. ANRS, Rega and Stanford algorithms were used for drug resistance interpretation.

Results

The study included 48 ART-naive patients. The number of patients with at least one major RAM was 3, 3, 4 and 8 when using RSS, DSS, NGS 20% and NGS 5%, respectively. Numerous minor mutations were detected in patients, especially in the protease gene. None of the methods detected any major mutation in the integrase gene. Overall, the mutation detection rate was similar between RSS and DSS, and higher with NGS 20%. Differences in drug resistance interpretation were found between algorithms. No impact of the minority RAMs detected by NGS was found on the short-term treatment outcome.

Conclusions

DSS does not clearly improve the detection of RAMs in ART-naive patients, as compared with RSS. NGS allows detection of additional minority RAMs; however, their clinical relevance requires further investigation.

Performance comparison of deep sequencing platforms for detecting HIV-1 variants in the pol gene

The present study compares the performances of an in-house sequencing protocol developed on MiSeq, the Sanger method, and the 454 GS-FLX for detecting and quantifying drug-resistant mutations (DRMs) in the human immunodeficiency virus polymerase gene (reverse transcriptase [RT] and protease [PR]). MiSeq sequencing identified all the resistance mutations detected by bulk sequencing (n = 84). Both the MiSeq and 454 GS-FLX platforms identified 67 DRMs in the RT and PR regions, but a further 25 DRMs were identified by only one or other of them. Pearson's analysis showed good concordance between the percentage of drug-resistant variants determined by MiSeq and 454 GS-FLX sequencing (ρ = .77, P < .0001). The MiSeq platform is as accurate as the 454 GS-FLX Roche system for determining RT and PR DRMs and could be used for monitoring human immunodeficiency virus type 1 drug resistance.

Antiretroviral resistance among HIV-1 patients on first-line therapy attending a comprehensive care clinic in Kenyatta National Hospital, Kenya: a retrospective analysis

INTRODUCTION

Antiretroviral therapy plays a major role in reducing the impact of Human Immunodeficiency Virus/Acquired Immune Disease Syndrome, especially in resource-limited settings. However, without proper infrastructure, it has resulted in emergence of drug resistance mutations in infected populations. To determine drug resistance mutations among patients attending a comprehensive care facility in Nairobi, 65 blood samples were successfully sequenced.

METHODS

Whole blood samples were also tested for CD4+T-cell count and plasma HIV-1 RNA Viral load. Drug-resistance testing targeting the HIV-1 RT gene was determined. Patients were on first line ART that consisted of two NRTIs, and one NNRTI.

RESULTS

Females were younger (mean 42) than males (mean 45) and lower median CD4+ counts (139 cells/μl) than males (152 cells/μl). The prevalence of drug resistance mutations (any major mutation) in this population was 23.1% (15/65). Major NRTI mutations were detected in 11 patient samples, which included M184V (n = 6), M41L (n=3), D67N (n=2), K219Q (n=3) and T215F (n=2). Major NNRTI mutations were detected in 14 patient samples. They included K103N (n = 10), G190A (n = 1), Y181C (n = 1) and Y188L (n = 1).

CONCLUSION

Presence of major mutations in this study calls for proper laboratory infrastructure to monitor treatment as well as regular appraisals of available regimens.

Ultra-Deep Sequencing Analysis on HIV Drug-Resistance-Associated Mutations Among HIV-Infected Individuals: First Report from the Philippines

A sharp increase in the number of people living with HIV has been documented in the Philippines. In response, the government has instituted antiretroviral therapy (ART) nationwide through HIV treatment hubs. However, no data presently exist on the status of ART drug-resistance-associated mutations (DRMs). In this study, we aim at analyzing DRM profiles in the Philippines and at providing comprehensive data on DRMs to guide treatment decisions and prevent viral failures. We conducted a cross-sectional study in 119 volunteers who tested positive for HIV from more than 8,000 participants screened for HIV across the nation through the 2013 Integrated HIV Behavioral and Serologic Surveillance (IHBSS) program. Amplicons were generated from plasma RNA by using primers designed to analyze diverse HIV-1 isolates targeting the reverse transcriptase region and sequenced on a 454 ultra-deep sequencing (UDS) platform to assess DRMs. DRMs were defined by using the Stanford HIV drug resistance database, and we found only 2 from 110 evaluable individuals with major HIV variants (>20% prevalence) that were highly resistant to the non-nucleoside reverse transcriptase inhibitor (NNRTI: efavirenz and nevirapine). However, a larger fraction of individuals harbored minority drug-resistant HIV variants (0.5%-20% prevalence) and they were highly resistant to NNRTI nevirapine (89/110), rilpivirine (5/110), and efavirenz (49/110). This study is the first report on the presence of HIV drug resistance in the Philippines and demonstrates the utility of UDS in assisting the detection of HIV minor variants. Monitoring for ART-DRMs will assist in improving HIV management strategies in curtailing the evolving epidemic in the Philippines.

A Follow-Up of the Multicenter Collaborative Study on HIV-1 Drug Resistance and Tropism Testing Using 454 Ultra Deep Pyrosequencing

BACKGROUND

Ultra deep sequencing is of increasing use not only in research but also in diagnostics. For implementation of ultra deep sequencing assays in clinical laboratories for routine diagnostics, intra- and inter-laboratory testing are of the utmost importance.

METHODS

A multicenter study was conducted to validate an updated assay design for 454 Life Sciences' GS FLX Titanium system targeting protease/reverse transcriptase (RTP) and env (V3) regions to identify HIV-1 drug-resistance mutations and determine co-receptor use with high sensitivity. The study included 30 HIV-1 subtype B and 6 subtype non-B samples with viral titers (VT) of 3,940-447,400 copies/mL, two dilution series (52,129-1,340 and 25,130-734 copies/mL), and triplicate samples. Amplicons spanning PR codons 10-99, RT codons 1-251 and the entire V3 region were generated using barcoded primers. Analysis was performed using the GS Amplicon Variant Analyzer and geno2pheno for tropism. For comparison, population sequencing was performed using the ViroSeq HIV-1 genotyping system.

RESULTS

The median sequencing depth across the 11 sites was 1,829 reads per position for RTP (IQR 592-3,488) and 2,410 for V3 (IQR 786-3,695). 10 preselected drug resistant variants were measured across sites and showed high inter-laboratory correlation across all sites with data (P<0.001). The triplicate samples of a plasmid mixture confirmed the high inter-laboratory consistency (mean% ± stdev: 4.6 ±0.5, 4.8 ±0.4, 4.9 ±0.3) and revealed good intra-laboratory consistency (mean% range ± stdev range: 4.2-5.2 ± 0.04-0.65). In the two dilutions series, no variants >20% were missed, variants 2-10% were detected at most sites (even at low VT), and variants 1-2% were detected by some sites. All mutations detected by population sequencing were also detected by UDS.

CONCLUSIONS

This assay design results in an accurate and reproducible approach to analyze HIV-1 mutant spectra, even at variant frequencies well below those routinely detectable by population sequencing.

Evaluation of Automatic Analysis of Ultradeep Pyrosequencing Raw Data to Determine Percentages of HIV Resistance Mutations in Patients Followed-Up in Hospital

A major obstacle to using next generation sequencing (NGS) technology in clinical routine practice is reliable data analysis. Thousands of sequences need to be aligned and validated, to exclude sequencing artifacts and generate accurate results. We compared two analysis pipelines for Roche 454 ultradeep pyrosequencing (UDPS) raw data generated from HIV-1 clinical samples: a commercial and fully automated Web-based software NGS HIV-1 Module (SmartGene, Zug, Switzerland) vs. the Amplicon Variant Analyzer software (AVA, 454 Life Sciences; Roche). Results were also compared to those obtained with Sanger sequencing. HIV-1 reverse transcriptase and protease genes from 34 plasma samples were submitted to Sanger sequencing and GS Junior UDPS. Raw UDPS data (sff files) from all samples were analyzed with AVA 2.7 software plus manual review of the alignments and the fully automated SmartGene NGS HIV-1 Module prototype (SMG). Results obtained with both analysis pipelines showed good correlation (85.0%). Divergent results were mainly observed at homopolymer positions, such as K101, where the frame-aware alignment and error corrections of the automated approach were more efficient and more accurate, both in terms of detecting and quantifying drug resistance mutations. Our study shows that NGS data can easily be analyzed via a fully automated analysis pipeline, here the SmartGene NGS HIV-1 Module, thus minimizing the need for manual review of alignments by the user, otherwise essential to ensure accurate results. Such automated analysis pipelines may facilitate the adoption of NGS platforms in the routine clinical laboratory.

Drug resistance and tropism as markers of the dynamics of HIV-1 DNA quasispecies in blood cells of heavily pretreated patients who achieved sustained virological suppression

OBJECTIVES

The objective of this study was to address the dynamics of archived resistant quasispecies in cell-associated HIV-1 DNA over time in heavily ART-experienced patients with currently suppressed plasma HIV-1 RNA.

METHODS

Longitudinal ultra-deep sequencing (UDS) analysis of reverse transcriptase, protease and V3 Env regions was performed on blood-cell-associated HIV-1 DNA samples. Drug-resistance-associated mutations (DRAMs) and tropism were interpreted using the ANRS and Geno2Pheno algorithms. We analysed frozen blood cells from patients enrolled in the INNOVE and ANRS 123 ETOILE studies who achieved sustained viral suppression after salvage optimized ART (SOT).

RESULTS

Samples were available at baseline and 6 and ≥12 months after SOT initiation in 10 patients. V3 loop sequences displayed wide intra-individual dynamics over time. Viral variants harbouring DRAMs exhibited three non-exclusive scenarios. First, when SOT exerted the same selective pressure as previous failing regimens, some viral quasispecies still harboured the same DRAMs at the same level as at the time of virological failure. Thus, as DRAMs were mostly associated with the same viral variant, variants with a complete resistance pattern remained archived. Second, some viral variants harbouring DRAMs were no longer detected over time when SOT consisted of new antiretroviral classes or had resistance profiles distinct from those of previous failing regimens. Third, variants with new DRAMs associated with SOT emerged in blood cells during follow-up despite sustained virological control.

CONCLUSIONS

Using longitudinal UDS analysis and focusing on DRAMs and tropism as markers, we demonstrated that, despite sustained virological control, archived HIV-1 DNA quasispecies continued to evolve.

Emergent HIV-1 Drug Resistance Mutations Were Not Present at Low-Frequency at Baseline in Non-Nucleoside Reverse Transcriptase Inhibitor-Treated Subjects in the STaR Study

At Week 96 of the Single-Tablet Regimen (STaR) study, more treatment-naïve subjects that received rilpivirine/emtricitabine/tenofovir DF (RPV/FTC/TDF) developed resistance mutations compared to those treated with efavirenz (EFV)/FTC/TDF by population sequencing. Furthermore, more RPV/FTC/TDF-treated subjects with baseline HIV-1 RNA >100,000 copies/mL developed resistance compared to subjects with baseline HIV-1 RNA ≤100,000 copies/mL. Here, deep sequencing was utilized to assess the presence of pre-existing low-frequency variants in subjects with and without resistance development in the STaR study. Deep sequencing (Illumina MiSeq) was performed on baseline and virologic failure samples for all subjects analyzed for resistance by population sequencing during the clinical study (n = 33), as well as baseline samples from control subjects with virologic response (n = 118). Primary NRTI or NNRTI drug resistance mutations present at low frequency (≥2% to 20%) were detected in 6.6% of baseline samples by deep sequencing, all of which occurred in control subjects. Deep sequencing results were generally consistent with population sequencing but detected additional primary NNRTI and NRTI resistance mutations at virologic failure in seven samples. HIV-1 drug resistance mutations emerging while on RPV/FTC/TDF or EFV/FTC/TDF treatment were not present at low frequency at baseline in the STaR study.

Comparison of 454 Ultra-Deep Sequencing and Allele-Specific Real-Time PCR with Regard to the Detection of Emerging Drug-Resistant Minor HIV-1 Variants after Antiretroviral Prophylaxis for Vertical Transmission

BACKGROUND

Pregnant HIV-infected women were screened for the development of HIV-1 drug resistance after implementation of a triple-antiretroviral transmission prophylaxis as recommended by the WHO in 2006. The study offered the opportunity to compare amplicon-based 454 ultra-deep sequencing (UDS) and allele-specific real-time PCR (ASPCR) for the detection of drug-resistant minor variants in the HIV-1 reverse transcriptase (RT).

METHODS

Plasma samples from 34 Tanzanian women were previously analysed by ASPCR for key resistance mutations in the viral RT selected by AZT, 3TC, and NVP (K70R, K103N, Y181C, M184V, T215Y/F). In this study, the RT region of the same samples was investigated by amplicon-based UDS for resistance mutations using the 454 GS FLX System.

RESULTS

Drug-resistant HIV-variants were identified in 69% (20/29) of women by UDS and in 45% (13/29) by ASPCR. The absolute number of resistance mutations identified by UDS was twice that identified by ASPCR (45 vs 24). By UDS 14 of 24 ASPCR-detected resistance mutations were identified at the same position. The overall concordance between UDS and ASPCR was 61.0% (25/41). The proportions of variants quantified by UDS were approximately 2-3 times lower than by ASPCR. Amplicon generation from samples with viral loads below 20,000 copies/ml failed more frequently by UDS compared to ASPCR (limit of detection = 650 copies/ml), resulting in missing or insufficient sequence coverage.

CONCLUSIONS

Both methods can provide useful information about drug-resistant minor HIV-1 variants. ASPCR has a higher sensitivity than UDS, but is restricted to single resistance mutations. In contrast, UDS is limited by its requirement for high viral loads to achieve sufficient sequence coverage, but the sequence information reveals the complete resistance patterns within the genomic region analysed. Improvements to the UDS limit of detection are in progress, and UDS could then facilitate monitoring of drug-resistant minor variants in the HIV-1 quasispecies.

Minority resistant HIV-1 variants and the response to first-line NNRTI therapy

BACKGROUND

The presence of low-frequency HIV-1 variants with mutations making them resistant to non-nucleoside reverse-transcriptase inhibitors (NNRTI) could influence the virological response to first-line NNRTI therapy.

OBJECTIVES

This study was designed to describe the proportions and quantities of NRTI and NNRTI-resistant variants in patients with successful first-line NNRTI therapy.

STUDY DESIGN

We evaluated the presence of drug-resistance mutations (DRMs) prior to treatment initiation in 131 naive chronically HIV-1-infected patients initiating NNRTI-based first-line therapy. DRMs were detected by ultradeep pyrosequencing (UDPS) on a GS Junior instrument (Roche).

RESULTS

The mean HIV RNA concentration was 4.78 ± 0.74 log copies/mL and the mean CD4 cell count was 368 ± 184 CD4 cells/mm(3). Patients were mainly infected with subtype B (68%) and 96% were treated with efavirenz. The sensitivity threshold for each mutation was 0.13-1.05% for 2000 reads. Major NRTI-resistant or NNRTI-resistant mutations were detected in 40 patients (33.6%). The median frequency of major NRTI-resistant mutations was 1.37% [IQR: 0.39-84.1], i.e.: a median of 556 copies/mL [IQR: 123-37,553]. The median frequency of major NNRTI-resistant DRMs was 0.78% [IQR: 0.67-7.06], i.e.: a median of 715 copies/mL [IQR: 391-3452]. The genotypic susceptibility score (GSS) of 9 (7.3%) patients with mutations to given treatment detected by UDPS was 1.5 or 2.

CONCLUSIONS

First-line NNRTI-based treatment can produce virological success in naïve HIV-1-infected patients harboring low-frequency DRMs representing <1% of the viral quasispecies. Further studies are needed to determine the clinical cut-off of low-frequency resistant variants associated to virological failure.

Burden of nonnucleoside reverse transcriptase inhibitor resistance in HIV-1-infected patients: a systematic review and meta-analysis

The prevalence of HIV drug resistance varies with geographic location, year, and treatment exposure. This study generated yearly estimates of nonnucleoside reverse transcriptase inhibitor (NNRTI) resistance in treatment-naive (TN) and treatment-experienced (TE) patients in the United States (US), Europe (EU), and Canada. Studies reporting NNRTI resistance identified in electronic databases and 11 conferences were analyzed in three groups: (1) TN patients in one of four geographic regions [US, Canada, EU countries with larger surveillance networks ("EU1"), and EU countries with fewer data ("EU2")]; (2) TE patients from any region; and (3) TN patients failing NNRTI-based treatments in clinical trials. Analysis data included 158 unique studies from 22 countries representing 84 cohorts of TN patients, 21 cohorts of TE patients, and 8 trials reporting resistance at failure. From 1995 to 2000, resistance prevalence in TN patients increased in US and EU1 from 3.1% to 7.5% and 0.8% to 3.6%, respectively. Resistance in both regions stabilized in 2006 onward. Little resistance was identified in EU2 before 2000, and increased from 2006 (5.0%) to 2010 (13.7%). One TN Canadian study was identified and reported resistance of 8.1% in 2006. Half of TN clinical trial patients had resistance after treatment failure at weeks 48-144. Resistance in TE patients increased from 1998 (10.1%) to 2001 (44.0%), then decreased after 2004. Trends in NNRTI resistance among TN patients show an increased burden in the US and some EU countries compared to others. These findings signify a need for alternate first-line treatments in some regions.

Real-life rilpivirine resistance and potential emergence of an E138A-positive HIV strain in north-eastern France

OBJECTIVES

To assess the prevalence of resistance to rilpivirine and mutations at position 138 in reverse transcriptase and to identify associated epidemiological and biological characteristics.

METHODS

This retrospective study included 238 patients with available HIV-1 nucleotide sequences analysed at the Laboratory of Virology at the University Hospital of Nancy between January 2011 and June 2013. Resistance to non-nucleoside reverse transcriptase inhibitors (NNRTIs) was evaluated according to the ANRS algorithm (version 23) and correlated with clinico-epidemiological and therapeutic data. The virus strains were analysed by evaluating the distance and distribution of the phylogenetic tree (MEGAv5).

RESULTS

Among previously treated patients (111/238, 46.6%), 68/111 (61.3%) had received NNRTIs; all were rilpivirine-naive. The prevalence of rilpivirine resistance in the whole cohort was 12.6% (30/238), and was 10.2% (13/127) and 15.3% (17/111) in naive and pre-treated patients, respectively. The E138A mutation was the most frequent mutation associated with resistance to rilpivirine (P < 0.0001). The prevalence of the E138A mutation tended to increase over time, from 3.6% (2/55) during the first half of 2011 to 9.3% (4/43) during the first half of 2013 (P = 0.0614). Seven viral strains from seven naive male patients positive for the E138A mutation appeared in the same cluster.

CONCLUSIONS

In our cohort of patients, we observed significantly increased resistance to rilpivirine, mostly because of the E138A mutation, probably due to an E138A strain circulating in newly diagnosed men who have sex with men. Taken together, our results emphasize the need to investigate the prevalence of rilpivirine resistance-associated mutations in the coming years both in France and abroad.

Prevalence of WHO transmitted drug resistance mutations by deep sequencing in antiretroviral-naïve subjects in Hunan Province, China

BACKGROUND

There are few data on the prevalence of WHO transmitted drug resistance mutations (TDRs) that could affect treatment responses to first line antiretroviral therapy (ART) in Hunan Province, China.

OBJECTIVE

Determine the prevalence of WHO NRTI/NNRTI/PI TDRs in ART-naïve subjects in Hunan Province by deep sequencing.

METHODS

ART-naïve subjects diagnosed in Hunan between 2010-2011 were evaluated by deep sequencing for low-frequency HIV variants possessing WHO TDRs to 1% levels. Mutations were scored using the HIVdb.stanford.edu algorithm to infer drug susceptibility.

RESULTS

Deep sequencing was performed on samples from 90 ART-naïve subjects; 83.3% were AE subtype. All subjects had advanced disease (average CD4 count 134 cells/mm3). Overall 25.6%(23/90) of subjects had HIV with major WHO NRTI/NNRTI TDRs by deep sequencing at a variant frequency level ≥ 1%; 16.7%(15/90) had NRTI TDR and 12.2%(11/90) had a major NNRTI TDR. The majority of NRTI/NNRTI mutations were identified at variant levels <5%. Mutations were analyzed by HIVdb.stanford.edu and 7.8% of subjects had variants with high-level nevirapine resistance; 4.4% had high-level NRTI resistance. Deep sequencing identified 24(27.6%) subjects with variants possessing either a PI TDR or hivdb.stanford.edu PI mutation (algorithm value ≥ 15). 17(19.5%) had PI TDRs at levels >1%.

CONCLUSIONS

ART-naïve subjects from Hunan Province China infected predominantly with subtype AE frequently possessed HIV variants with WHO NRTI/NNRTI TDRs by deep sequencing that would affect the first line ART used in the region. Specific mutations conferring nevirapine high-level resistance were identified in 7.8% of subjects. The majority of TDRs detected were at variant levels <5% likely due to subjects having advanced chronic disease at the time of testing. PI TDRs were identified frequently, but were found in isolation and at low variant frequency. As PI/r use is infrequent in Hunan, the existence of PI mutations likely represent AE subtype natural polymorphism at low variant level frequency.

Comparison of ultra-deep versus Sanger sequencing detection of minority mutations on the HIV-1 drug resistance interpretations after virological failure

OBJECTIVE

Drug-resistance mutations are routinely detected using standard Sanger sequencing, which does not detect minor variants with a frequency below 20%. The impact of detecting minor variants generated by ultra-deep sequencing (UDS) on HIV drug-resistance interpretations has not yet been studied.

DESIGN

Fifty HIV-1 patients who experienced virological failure were included in this retrospective study.

METHODS

The HIV-1 UDS protocol allowed the detection and quantification of HIV-1 protease and reverse transcriptase variants related to genotypes A, B, C, F and G. DeepChek-HIV simplified drug-resistance interpretation software was used to compare Sanger sequencing and UDS.

RESULTS

The total time required for the UDS protocol was found to be approximately three times longer than Sanger sequencing with equivalent reagent costs. UDS detected all of the mutations found by population sequencing and identified additional resistance variants in all patients. An analysis of drug resistance revealed a total of 643 and 224 clinically relevant mutations by UDS and Sanger sequencing, respectively. Three resistance mutations with more than 20% prevalence were detected solely by UDS: A98S (23%), E138A (21%) and V179I (25%). A significant difference in the drug-resistance interpretations for 19 antiretroviral drugs was observed between the UDS and Sanger sequencing methods. Y181C and T215Y were the most frequent mutations associated with interpretation differences.

CONCLUSION

A combination of UDS and DeepChek software for the interpretation of drug resistance results would help clinicians provide suitable treatments. A cut-off of 1% allowed a better characterization of the viral population by identifying additional resistance mutations and improving the drug-resistance interpretation.

Deep sequencing of HIV: clinical and research applications

Human immunodeficiency virus (HIV) exhibits remarkable diversity in its genomic makeup and exists in any given individual as a complex distribution of closely related but nonidentical genomes called a viral quasispecies, which is subject to genetic variation, competition, and selection. This viral diversity clinically manifests as a selection of mutant variants based on viral fitness in treatment-naive individuals and based on drug-selective pressure in those on antiretroviral therapy (ART). The current standard-of-care ART consists of a combination of antiretroviral agents, which ensures maximal viral suppression while preventing the emergence of drug-resistant HIV variants. Unfortunately, transmission of drug-resistant HIV does occur, affecting 5% to >20% of newly infected individuals. To optimize therapy, clinicians rely on viral genotypic information obtained from conventional population sequencing-based assays, which cannot reliably detect viral variants that constitute <20% of the circulating viral quasispecies. These low-frequency variants can be detected by highly sensitive genotyping methods collectively grouped under the moniker of deep sequencing. Low-frequency variants have been correlated to treatment failures and HIV transmission, and detection of these variants is helping to inform strategies for vaccine development. Here, we discuss the molecular virology of HIV, viral heterogeneity, drug-resistance mutations, and the application of deep sequencing technologies in research and the clinical care of HIV-infected individuals.

An international multicenter study on HIV-1 drug resistance testing by 454 ultra-deep pyrosequencing

The detection of mutant spectra within the viral quasispecies is critical for therapeutic management of HIV-1 infections. Routine clinical application of ultrasensitive genotyping requires reproducibility and concordance within and between laboratories. The goal of the study was to evaluate a new protocol on HIV-1 drug resistance testing by 454 ultra-deep pyrosequencing (454-UDS) in an international multicenter study. Sixteen blinded HIV-1 subtype B samples were provided for 454-UDS as both RNA and cDNA with viral titers of 88,600-573,000 HIV-1 RNA copies/ml. Eight overlapping amplicons spanning protease (PR) codons 10-99 and reverse transcriptase (RT) codons 1-251 were generated using molecular barcoded primers. 454-UDS was performed using the 454 Life Sciences/Roche GS FLX platform. PR and RT sequences were analyzed using 454 Life Sciences Amplicon Variant Analyzer (AVA) software. Quantified variation data were analyzed for intra-laboratory reproducibility and inter-laboratory concordance. Routine population sequencing was performed using the ViroSeq HIV-1 genotyping system. Eleven laboratories and the reference laboratory 454 Life Sciences sequenced the HIV-1 sample set. Data presented are derived from seven laboratories and the reference laboratory since severe study protocol execution errors occurred in four laboratories leading to exclusion. The median sequencing depth across all sites was 1364 reads per position (IQR=809-2065). 100% of the ViroSeq-reported mutations were also detected by 454-UDS. Minority HIV-1 drug resistance mutations, defined as HIV-1 drug resistance mutations identified at frequencies of 1-25%, were only detected by 454-UDS. Analysis of 10 preselected majority and minority mutations were consistently found across sites. The analysis of drug-resistance mutations detected between 1 and 10% demonstrated high intra- and inter-laboratory consistency in frequency estimates for both RNA and prepared cDNA samples, indicating robustness of the method. HIV-1 drug resistance testing using 454 ultra-deep pyrosequencing results in an accurate and highly reproducible, albeit complex, approach to the analysis of HIV-1 mutant spectra, even at frequencies well below those detected by routine population sequencing.

Detection of low-frequency HIV type 1 reverse transcriptase drug resistance mutations by ultradeep sequencing in naive HIV type 1-infected individuals

Genotypic resistance testing is recommended to evaluate the susceptibility of HIV to antiretroviral drugs. These tests are based on bulk population sequencing and thus consider only variants representing more than 20% of the viral population, whereas next generation sequencing methods allow detection below this threshold. We aimed to evaluate the potential use of ultradeep pyrosequencing (UDPS) for genotypic resistance testing in clinical routine at the University Hospital of Bordeaux, France. We performed UDPS on reverse transcriptase (RT) from 47 HIV-1 individuals, naive of antiretroviral treatment and for whom genotypic resistance testing was requested for clinical management in 2011-2012. In 8.5% of the patients, only low-frequency variants harboring RT drug resistance mutations were detected raising the question of their clinical significance. Rilpivirine-associated resistance mutations were detected in 19.1% of our population study. To conclude, UDPS could become a routine tool for the evaluation of HIV-infected patients in hospital laboratories.

Role of the K101E substitution in HIV-1 reverse transcriptase in resistance to rilpivirine and other nonnucleoside reverse transcriptase inhibitors

Resistance to the recently approved nonnucleoside reverse transcriptase inhibitor (NNRTI) rilpivirine (RPV) commonly involves substitutions at positions E138K and K101E in HIV-1 reverse transcriptase (RT), together with an M184I substitution that is associated with resistance to coutilized emtricitabine (FTC). Previous biochemical and virological studies have shown that compensatory interactions between substitutions E138K and M184I can restore enzyme processivity and the viral replication capacity. Structural modeling studies have also shown that disruption of the salt bridge between K101 and E138 can affect RPV binding. The current study was designed to investigate the impact of K101E, alone or in combination with E138K and/or M184I, on drug susceptibility, viral replication capacity, and enzyme function. We show here that K101E can be selected in cell culture by the NNRTIs etravirine (ETR), efavirenz (EFV), and dapivirine (DPV) as well as by RPV. Recombinant RT enzymes and viruses containing K101E, but not E138K, were highly resistant to nevirapine (NVP) and delavirdine (DLV) as well as ETR and RPV, but not EFV. The addition of K101E to E138K slightly enhanced ETR and RPV resistance compared to that obtained with E138K alone but restored susceptibility to NVP and DLV. The K101E substitution can compensate for deficits in viral replication capacity and enzyme processivity associated with M184I, while M184I can compensate for the diminished efficiency of DNA polymerization associated with K101E. The coexistence of K101E and E138K does not impair either viral replication or enzyme fitness. We conclude that K101E can play a significant role in resistance to RPV.

Simultaneous detection of major drug resistance mutations in the protease and reverse transcriptase genes for HIV-1 subtype C by use of a multiplex allele-specific assay

High-throughput, sensitive, and cost-effective HIV drug resistance (HIVDR) detection assays are needed for large-scale monitoring of the emergence and transmission of HIVDR in resource-limited settings. Using suspension array technology, we have developed a multiplex allele-specific (MAS) assay that can simultaneously detect major HIVDR mutations at 20 loci. Forty-five allele-specific primers tagged with unique 24-base oligonucleotides at the 5' end were designed to detect wild-type and mutant alleles at the 20 loci of HIV-1 subtype C. The MAS assay was first established and optimized with three plasmid templates (C-wt, C-mut1, and C-mut2) and then evaluated using 148 plasma specimens from HIV-1 subtype C-infected individuals. All the wild-type and mutant alleles were unequivocally distinguished with plasmid templates, and the limits of detection were 1.56% for K219Q and K219E, 3.13% for L76V, 6.25% for K65R, K70R, L74V, L100I, K103N, K103R, Q151M, Y181C, and I47V, and 12.5% for M41L, K101P, K101E, V106A, V106M, Y115F, M184V, Y188L, G190A, V32I, I47A, I84V, and L90M. Analyses of 148 plasma specimens revealed that the MAS assay gave 100% concordance with conventional sequencing at eight loci and >95% (range, 95.21% to 99.32%) concordance at the remaining 12 loci. The differences observed were caused mainly by 24 additional low-abundance alleles detected by the MAS assay. Ultradeep sequencing analysis confirmed 15 of the 16 low-abundance alleles. This multiplex, sensitive, and straightforward result-reporting assay represents a new efficient genotyping tool for HIVDR surveillance and monitoring.

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

BACKGROUND

Drug-resistant human immunodeficiency virus type 1 (HIV-1) minority variants (MVs) are present in some antiretroviral therapy (ART)-naive patients. They may result from de novo mutagenesis or transmission. To date, the latter has not been proven.

METHODS

MVs were quantified by allele-specific polymerase chain reaction in 204 acute or recent seroconverters from the Zurich Primary HIV Infection study and 382 ART-naive, chronically infected patients. Phylogenetic analyses identified transmission clusters.

RESULTS

Three lines of evidence were observed in support of transmission of MVs. First, potential transmitters were identified for 12 of 16 acute or recent seroconverters harboring M184V MVs. These variants were also detected in plasma and/or peripheral blood mononuclear cells at the estimated time of transmission in 3 of 4 potential transmitters who experienced virological failure accompanied by the selection of the M184V mutation before transmission. Second, prevalence between MVs harboring the frequent mutation M184V and the particularly uncommon integrase mutation N155H differed highly significantly in acute or recent seroconverters (8.2% vs 0.5%; P < .001). Third, the prevalence of less-fit M184V MVs is significantly higher in acutely or recently than in chronically HIV-1-infected patients (8.2% vs 2.5%; P = .004).

CONCLUSIONS

Drug-resistant HIV-1 MVs can be transmitted. To what extent the origin-transmission vs sporadic appearance-of these variants determines their impact on ART needs to be further explored.