Comparative performance of high-density oligonucleotide sequencing and dideoxynucleotide sequencing of HIV type 1 pol from clinical samples

HIV-1 Drug Resistance Detected by Next-Generation Sequencing among ART-Naïve Individuals: A Systematic Review and Meta-Analysis

BACKGROUND

There are an increasing number of articles focused on the prevalence and clinical impact of pretreatment HIV drug resistance (PDR) detected by Sanger sequencing (SGS). PDR may contribute to the increased likelihood of virologic failure and the emergence of new resistance mutations. As SGS is gradually replaced by next-generation sequencing (NGS), it is necessary to assess the levels of PDR using NGS in ART-naïve patients systematically. NGS can detect the viral variants (low-abundance drug-resistant HIV-1 variants (LA-DRVs)) of virus quasi-species at levels below 20% that SGS may fail to detect. NGS has the potential to optimize current HIV drug resistance surveillance methods and inform future research directions. As the NGS technique has high sensitivity, it is highly likely that the level of pretreatment resistance would be underestimated using conventional techniques.

METHODS

For the systematic review and meta-analysis, we searched for original studies published in PubMed, Web of Science, Scopus, and Embase before 30 March 2023 that focused exclusively on the application of NGS in the detection of HIV drug resistance. Pooled prevalence estimates were calculated using a random effects model using the 'meta' package in R (version 4.2.3). We described drug resistance detected at five thresholds (>1%, 2%, 5%, 10%, and 20% of virus quasi-species). Chi-squared tests were used to analyze differences between the overall prevalence of PDR reported by SGS and NGS.

RESULTS

A total of 39 eligible studies were selected. The studies included a total of 15,242 ART-naïve individuals living with HIV. The prevalence of PDR was inversely correlated with the mutation detection threshold. The overall prevalence of PDR was 29.74% at the 1% threshold, 22.43% at the 2% threshold, 15.47% at the 5% threshold, 12.95% at the 10% threshold, and 11.08% at the 20% threshold. The prevalence of PDR to INSTIs was 1.22% (95%CI: 0.58-2.57), which is the lowest among the values for all antiretroviral drugs. The prevalence of LA-DRVs was 9.45%. At the 2% and 20% detection threshold, the prevalence of PDR was 22.43% and 11.08%, respectively. Resistance to PIs and INSTIs increased 5.52-fold and 7.08-fold, respectively, in those with a PDR threshold of 2% compared with those with PDR at 20%. However, resistance to NRTIs and NNRTIs increased 2.50-fold and 2.37-fold, respectively. There was a significant difference between the 2% and 5% threshold for detecting HIV drug resistance. There was no statistically significant difference between the results reported by SGS and NGS when using the 20% threshold for reporting resistance mutations.

CONCLUSION

In this study, we found that next-generation sequencing facilitates a more sensitive detection of HIV-1 drug resistance than SGS. The high prevalence of PDR emphasizes the importance of baseline resistance and assessing the threshold for optimal clinical detection using NGS.

HIV-1 Drug Resistance Assay Using Ion Torrent Next Generation Sequencing and On-Instrument End-to-End Analysis Software

HIV-1 antiretroviral therapy management requires sequencing the protease, reverse transcriptase, and integrase portions of the HIV-1 pol gene. Most resistance testing is performed with Sanger sequencing, which has limited ability to detect minor variants. Next generation sequencing (NGS) platforms enable variant detection at frequencies as low as 1% allowing for earlier detection of resistance and modification of therapy. Implementation of NGS assays in the clinical laboratory is hindered by complicated assay design, cumbersome wet bench procedures, and the complexity of data analysis and bioinformatics. We developed a complete NGS protocol and companion analysis and reporting pipeline using AmpliSeq multiplex PCR, Ion Torrent S5 XL sequencing, and Stanford's HIVdb resistance algorithm. Implemented as a Torrent Suite software plugin, the pipeline runs automatically after sequencing. An optimum variant frequency threshold of 10% was determined by comparing Sanger sequences of archived samples from ViroSeq testing, resulting in a sensitivity of 98.2% and specificity of 99.0%. The majority (91%) of drug resistance mutations were detected by both Sanger and NGS, with 1.7% only by Sanger and 7.3% only by NGS. Variant calls were highly reproducible and there was no cross-reactivity to VZV, HBV, CMV, EBV, and HCV. The limit of detection was 500 copies/mL. The NGS assay performance was comparable to ViroSeq Sanger sequencing and has several advantages, including a publicly available end-to-end analysis and reporting plugin. The assay provides a straightforward path for implementation of NGS for HIV drug resistance testing in the laboratory setting without additional investment in bioinformatics infrastructure and resources.

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.

Near point-of-care, point-mutation test to detect drug resistance in HIV-1: a validation study in a Mexican cohort

OBJECTIVE

Pretreatment HIV-drug resistance (PDR, HIVDR) to non-nucleoside reverse transcriptase inhibitors (NNRTIs) is increasing globally. NNRTIs continue to be used as first-line antiretroviral therapy (ART) in some communities due to the cost of dolutegravir-based ART or dolutegravir-associated adverse events. A simplified version of the oligonucleotide ligation assay (OLA) - 'OLA-Simple' - is a low-cost, near point-of-care assay that provides ready-to-use lyophilized reagents and reports HIVDR mutations as colored lines on lateral flow strips. Our objective was to design and validate OLA-Simple for a Mexican cohort.

DESIGN

OLA-Simple probes to detect K65R, K103N/S, Y181C, M184V, and G190A were optimized for HIV Mexican sequences. Sixty clinical plasma specimens were analyzed by OLA-Simple by technicians blinded to Illumina-MiSeq sequences, and HIVDR results were compared.

METHODS

Plasma RNA was tested using OLA-Simple kits. OLA-Simple lateral flow strips were read by in-house software, and were classified as mutant or wild-type at each codon. The comparison of results by OLA-Simple and Miseq was used to generate receiver-operating characteristic curves.

RESULTS

OLA-Simple PCR amplified 59 of 60 specimens and successfully genotyped 287 of 295 codons, with eight of 295 (2.7%) indeterminate results. Compared to MiSeq, OLA-Simple gave five of 295 (1.7%) false-positive and four of 295 (1.4%) false-negative results. Excluding indeterminate results, OLA-Simple classified mutant with an accuracy of 97.4 and 98.8% when using thresholds at 10 and 25% mutant within an individual's HIV quasispecies, respectively.

CONCLUSIONS

Compared to MiSeq, OLA-Simple detected HIVDR with high sensitivity and accuracy. OLA-Simple could expand access to affordable and rapid HIVDR testing to guide appropriate ART choices in populations using NNRTI-based ART.

Next-Generation Sequencing for HIV Drug Resistance Testing: Laboratory, Clinical, and Implementation Considerations

Higher accessibility and decreasing costs of next generation sequencing (NGS), availability of commercial kits, and development of dedicated analysis pipelines, have allowed an increasing number of laboratories to adopt this technology for HIV drug resistance (HIVDR) genotyping. Conventional HIVDR genotyping is traditionally carried out using population-based Sanger sequencing, which has a limited capacity for reliable detection of variants present at intra-host frequencies below a threshold of approximately 20%. NGS has the potential to improve sensitivity and quantitatively identify low-abundance variants, improving efficiency and lowering costs. However, some challenges exist for the standardization and quality assurance of NGS-based HIVDR genotyping. In this paper, we highlight considerations of these challenges as related to laboratory, clinical, and implementation of NGS for HIV drug resistance testing. Several sources of variation and bias occur in each step of the general NGS workflow, i.e., starting material, sample type, PCR amplification, library preparation method, instrument and sequencing chemistry-inherent errors, and data analysis options and limitations. Additionally, adoption of NGS-based HIVDR genotyping, especially for clinical care, poses pressing challenges, especially for resource-poor settings, including infrastructure and equipment requirements and cost, logistic and supply chains, instrument service availability, personnel training, validated laboratory protocols, and standardized analysis outputs. The establishment of external quality assessment programs may help to address some of these challenges and is needed to proceed with NGS-based HIVDR genotyping adoption.

Viral Diversity Based on Next-Generation Sequencing of HIV-1 Provides Precise Estimates of Infection Recency and Time Since Infection

BACKGROUND

Human immunodeficiency virus type 1 (HIV-1) genetic diversity increases over the course of infection and can be used to infer the time since infection and, consequently, infection recency, which are crucial for HIV-1 surveillance and the understanding of viral pathogenesis.

METHODS

We considered 313 HIV-infected individuals for whom reliable estimates of infection dates and next-generation sequencing (NGS)-derived nucleotide frequency data were available. Fractions of ambiguous nucleotides, obtained by population sequencing, were available for 207 samples. We assessed whether the average pairwise diversity calculated using NGS sequences provided a more exact prediction of the time since infection and classification of infection recency (<1 year after infection), compared with the fraction of ambiguous nucleotides.

RESULTS

NGS-derived average pairwise diversity classified an infection as recent with a sensitivity of 88% and a specificity of 85%. When considering only the 207 samples for which fractions of ambiguous nucleotides were available, the NGS-derived average pairwise diversity exhibited a higher sensitivity (90% vs 78%) and specificity (95% vs 67%) than the fraction of ambiguous nucleotides. Additionally, the average pairwise diversity could be used to estimate the time since infection with a mean absolute error of 0.84 years, compared with 1.03 years for the fraction of ambiguous nucleotides.

CONCLUSIONS

Viral diversity based on NGS data is more precise than that based on population sequencing in its ability to predict infection recency and provides an estimated time since infection that has a mean absolute error of <1 year.

Performance comparison of next generation sequencing analysis pipelines for HIV-1 drug resistance testing

Next generation sequencing (NGS) is a trending new standard for genotypic HIV-1 drug resistance (HIVDR) testing. Many NGS HIVDR data analysis pipelines have been independently developed, each with variable outputs and data management protocols. Standardization of such analytical methods and comparison of available pipelines are lacking, yet may impact subsequent HIVDR interpretation and other downstream applications. Here we compared the performance of five NGS HIVDR pipelines using proficiency panel samples from NIAID Virology Quality Assurance (VQA) program. Ten VQA panel specimens were genotyped by each of six international laboratories using their own in-house NGS assays. Raw NGS data were then processed using each of the five different pipelines including HyDRA, MiCall, PASeq, Hivmmer and DEEPGEN. All pipelines detected amino acid variants (AAVs) at full range of frequencies (1~100%) and demonstrated good linearity as compared to the reference frequency values. While the sensitivity in detecting low abundance AAVs, with frequencies between 1~20%, is less a concern for all pipelines, their specificity dramatically decreased at AAV frequencies <2%, suggesting that 2% threshold may be a more reliable reporting threshold for ensured specificity in AAV calling and reporting. More variations were observed among the pipelines when low abundance AAVs are concerned, likely due to differences in their NGS read quality control strategies. Findings from this study highlight the need for standardized strategies for NGS HIVDR data analysis, especially for the detection of minority HIVDR variants.

OLA-Simple: A software-guided HIV-1 drug resistance test for low-resource laboratories

Background

HIV drug resistance (HIVDR) testing can assist clinicians in selecting treatments. However, high complexity and cost of genotyping assays limit routine testing in settings where HIVDR prevalence has reached high levels.

Methods

The oligonucleotide ligation assay (OLA)-Simple kit was developed for detection of HIVDR against first-line non-nucleoside/nucleoside reverse transcriptase inhibitors and validated on 672 codons (168 specimens) from subtypes A, B, C, D, and AE. The kit uses dry reagents to facilitate assay setup, lateral flow devices for visual HIVDR detections, and in-house software with an interface for guiding users and analyzing results.

Findings

HIVDR analysis of specimens by OLA-Simple compared to Sanger sequencing revealed 99.6 ± 0.3% specificity and 98.2 ± 0.9% sensitivity, and compared to high-sensitivity assays, 99.6 ± 0.6% specificity and 86.2 ± 2.5% sensitivity, with 2.6 ± 0.9% indeterminate results. OLA-Simple was performed more rapidly compared to Sanger sequencing (<4 h vs. 35–72 h). Forty-one untrained volunteers blindly tested two specimens each with 96.8 ± 0.8% accuracy.

Interpretation

OLA-Simple compares favorably with HIVDR genotyping by Sanger and sensitive comparators. Instructional software enabled inexperienced, first-time users to perform the assay with high accuracy. The reduced complexity, cost, and training requirements of OLA-Simple could improve access to HIVDR testing in low-resource settings and potentially allow same-day selection of appropriate antiretroviral therapy.

Fund

USA National Institutes of Health R01; the Clinical and Retrovirology Research Core and the Molecular Profiling and Computational Biology Core of the UW CFAR; Seattle Children's Research Institute; UW Holloman Innovation Challenge Award; Pilcher Faculty Fellowship.

A MiSeq-HyDRA platform for enhanced HIV drug resistance genotyping and surveillance

Conventional HIV drug resistance (HIVDR) genotyping utilizes Sanger sequencing (SS) methods, which are limited by low data throughput and the inability of detecting low abundant drug resistant variants (LADRVs). Here we present a next generation sequencing (NGS)-based HIVDR typing platform that leverages the advantages of Illumina MiSeq and HyDRA Web. The platform consists of a fully validated sample processing protocol and HyDRA web, an open web portal that allows automated customizable NGS-based HIVDR data processing. This platform was characterized and validated using a panel of HIV-spiked plasma representing all major HIV-1 subtypes, pedigreed plasmids, HIVDR proficiency specimens and clinical specimens. All examined major HIV-1 subtypes were consistently amplified at viral loads of ≥1,000 copies/ml. The gross error rate of this platform was determined at 0.21%, and minor variations were reliably detected down to 0.50% in plasmid mixtures. All HIVDR mutations identifiable by SS were detected by the MiSeq-HyDRA protocol, while LADRVs at frequencies of 1~15% were detected by MiSeq-HyDRA only. As compared to SS approaches, the MiSeq-HyDRA platform has several notable advantages including reduced cost and labour, and increased sensitivity for LADRVs, making it suitable for routine HIVDR monitoring for both patient care and surveillance purposes.

Next-generation sequencing of HIV-1 single genome amplicons

The analysis of HIV-1 sequences has helped understand the viral molecular epidemiology, monitor the development of antiretroviral drug resistance, and design candidate vaccines. The introduction of single genome amplification (SGA) has been a major advancement in the field, allowing for the characterization of multiple sequences per patient while preserving linkage among polymorphisms in the same viral genome copy. Sequencing of SGA amplicons is performed by capillary Sanger sequencing, which presents low throughput, requires a high amount of template, and is highly sensitive to template/primer mismatching. In order to meet the increasing demand for HIV-1 SGA amplicon sequencing, we have developed a platform based on benchtop next-generation sequencing (NGS) (IonTorrent) accompanied by a bioinformatics pipeline capable of running on computer resources commonly available at research laboratories. During assay validation, the NGS-based sequencing of 10 HIV-1 env SGA amplicons was fully concordant with Sanger sequencing. The field test was conducted on plasma samples from 10 US Navy and Marine service members with recent HIV-1 infection (sampling interval: 2005–2010; plasma viral load: 5,884–194,984 copies/ml). The NGS analysis of 101 SGA amplicons (median: 10 amplicons/individual) showed within-individual viral sequence profiles expected in individuals at this disease stage, including individuals with highly homogeneous quasispecies, individuals with two highly homogeneous viral lineages, and individuals with heterogeneous viral populations. In a scalability assessment using the Ion Chef automated system, 41/43 tested env SGA amplicons (95%) multiplexed on a single Ion 318 chip showed consistent gene-wide coverage >50×. With lower sample requirements and higher throughput, this approach is suitable to support the increasing demand for high-quality and cost-effective HIV-1 sequences in fields such as molecular epidemiology, and development of preventive and therapeutic strategies.

Bioinformatic data processing pipelines in support of next-generation sequencing-based HIV drug resistance testing: the Winnipeg Consensus

INTRODUCTION

Next-generation sequencing (NGS) has several advantages over conventional Sanger sequencing for HIV drug resistance (HIVDR) genotyping, including detection and quantitation of low-abundance variants bearing drug resistance mutations (DRMs). However, the high HIV genomic diversity, unprecedented large volume of data, complexity of analysis and potential for error pose significant challenges for data processing. Several NGS analysis pipelines have been developed and used in HIVDR research; however, the absence of uniformity in data processing strategies results in lack of consistency and comparability of outputs from different pipelines. To fill this gap, an international symposium on bioinformatic strategies for NGS-based HIVDR testing was held in February 2018 in Winnipeg, Canada, convening laboratory scientists, bioinformaticians and clinicians involved in four recently developed, publicly available NGS HIVDR pipelines. The goal of this symposium was to establish a consensus on effective bioinformatic strategies for NGS data management and its use for HIVDR reporting.

DISCUSSION

Essential functionalities of an NGS HIVDR pipeline were divided into five analytic blocks: (1) NGS read quality control (QC)/quality assurance (QA); (2) NGS read alignment and reference mapping; (3) HIV variant calling and variant QC; (4) NGS HIVDR reporting; and (5) extended data applications and additional considerations for data management. The consensuses reached among the participants on all major aspects of these blocks are summarized here. They encompass not only recommended data management and analysis strategies, but also detailed bioinformatic approaches that help ensure accuracy of the derived HIVDR analysis outputs for both research and potential clinical use.

CONCLUSIONS

While NGS is being adopted more broadly in HIVDR testing laboratories, data processing is often a bottleneck hindering its generalized application. The proposed standardization of NGS read QC/QA, read alignment and reference mapping, variant calling and QC, HIVDR reporting and relevant data management strategies in this "Winnipeg Consensus" may serve as a starting guideline for NGS HIVDR data processing that informs the refinement of existing pipelines and those yet to be developed. Moreover, the bioinformatic strategies presented here may apply more broadly to NGS data analysis of microbes harbouring significant genomic diversity.

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

Background

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

Objective

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

Methods

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

Results

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

Conclusion

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

Molecular evolution of HIV-1 integrase during the 20 years prior to the first approval of integrase inhibitors

BACKGROUND

Detailed knowledge of the evolutionary potential of polymorphic sites in a viral protein is important for understanding the development of drug resistance in the presence of an inhibitor. We therefore set out to analyse the molecular evolution of the HIV-1 subtype B integrase at the inter-patient level in Germany during a 20-year period prior to the first introduction of integrase strand inhibitors (INSTIs).

METHODS

We determined 337 HIV-1 integrase subtype B sequences (amino acids 1-278) from stored plasma samples of antiretroviral treatment-naïve individuals newly diagnosed with HIV-1 between 1986 and 2006. Shannon entropy was calculated to determine the variability at each amino acid position. Time trends in the frequency of amino acid variants were identified by linear regression. Direct coupling analysis was applied to detect covarying sites.

RESULTS

Twenty-two time trends in the frequency of amino acid variants demonstrated either single amino acid exchanges or variation in the degree of polymorphy. Covariation was observed for 17 amino acid variants with a temporal trend. Some minor INSTI resistance mutations (T124A, V151I, K156 N, T206S, S230 N) and some INSTI-selected mutations (M50I, L101I, T122I, T124 N, T125A, M154I, G193E, V201I) were identified at overall frequencies >5%. Among these, the frequencies of L101I, T122I, and V201I increased over time, whereas the frequency of M154I decreased. Moreover, L101I, T122I, T124A, T125A, M154I, and V201I covaried with non-resistance-associated variants.

CONCLUSIONS

Time-trending, covarying polymorphisms indicate that long-term evolutionary changes of the HIV-1 integrase involve defined clusters of possibly structurally or functionally associated sites independent of selective pressure through INSTIs at the inter-patient level. Linkage between polymorphic resistance- and non-resistance-associated sites can impact the selection of INSTI resistance mutations in complex ways. Identification of these sites can help in improving genotypic resistance assays, resistance prediction algorithms, and the development of new integrase inhibitors.

High prevalence of HIV-1 transmitted drug-resistance mutations from proviral DNA massively parallel sequencing data of therapy-naïve chronically infected Brazilian blood donors

BACKGROUND

An improved understanding of the prevalence of low-abundance transmitted drug-resistance mutations (TDRM) in therapy-naïve HIV-1-infected patients may help determine which patients are the best candidates for therapy. In this study, we aimed to obtain a comprehensive picture of the evolving HIV-1 TDRM across the massive parallel sequences (MPS) of the viral entire proviral genome in a well-characterized Brazilian blood donor naïve to antiretroviral drugs.

MATERIALS AND METHODS

The MPS data from 128 samples used in the analysis were sourced from Brazilian blood donors and were previously classified by less-sensitive (LS) or "detuned" enzyme immunoassay as non-recent or longstanding HIV-1 infections. The Stanford HIV Resistance Database (HIVDBv 6.2) and IAS-USA mutation lists were used to interpret the pattern of drug resistance. The minority variants with TDRM were identified using a threshold of ≥ 1.0% and ≤ 20% of the reads sequenced. The rate of TDRM in the MPS data of the proviral genome were compared with the corresponding published consensus sequences of their plasma viruses.

RESULTS

No TDRM were detected in the integrase or envelope regions. The overall prevalence of TDRM in the protease (PR) and reverse transcriptase (RT) regions of the HIV-1 pol gene was 44.5% (57/128), including any mutations to the nucleoside analogue reverse transcriptase inhibitors (NRTI) and non-nucleoside analogue reverse transcriptase inhibitors (NNRTI). Of the 57 subjects, 43 (75.4%) harbored a minority variant containing at least one clinically relevant TDRM. Among the 43 subjects, 33 (76.7%) had detectable minority resistant variants to NRTIs, 6 (13.9%) to NNRTIs, and 16 (37.2%) to PR inhibitors. The comparison of viral sequences in both sources, plasma and cells, would have detected 48 DNA provirus disclosed TDRM by MPS previously missed by plasma bulk analysis.

CONCLUSION

Our findings revealed a high prevalence of TDRM found in this group, as the use of MPS drastically increased the detection of these mutations. Sequencing proviral DNA provided additional information about TDRM, which may impact treatment decisions. The overall results emphasize the importance of continuous monitoring.

The HIV-1 transmission bottleneck

It is well established that most new systemic infections of HIV-1 can be traced back to one or a limited number of founder viruses. Usually, these founders are more closely related to minor HIV-1 populations in the blood of the presumed donor than to more abundant lineages. This has led to the widely accepted idea that transmission selects for viral characteristics that facilitate crossing the mucosal barrier of the recipient’s genital tract, although the specific selective forces or advantages are not completely defined. However, there are other steps along the way to becoming a founder virus at which selection may occur. These steps include the transition from the donor’s general circulation to the genital tract compartment, survival within the transmission fluid, and establishment of a nascent stable local infection in the recipient’s genital tract. Finally, there is the possibility that important narrowing events may also occur during establishment of systemic infection. This is suggested by the surprising observation that the number of founder viruses detected after transmission in intravenous drug users is also limited. Although some of these steps may be heavily selective, others may result mostly in a stochastic narrowing of the available founder pool. Collectively, they shape the initial infection in each recipient.

Simplified Paper Format for Detecting HIV Drug Resistance in Clinical Specimens by Oligonucleotide Ligation

Human immunodeficiency virus (HIV) is a chronic infection that can be managed by antiretroviral treatment (ART). However, periods of suboptimal viral suppression during lifelong ART can select for HIV drug resistant (DR) variants. Transmission of drug resistant virus can lessen or abrogate ART efficacy. Therefore, testing of individuals for drug resistance prior to initiation of treatment is recommended to ensure effective ART. Sensitive and inexpensive HIV genotyping methods are needed in low-resource settings where most HIV infections occur. The oligonucleotide ligation assay (OLA) is a sensitive point mutation assay for detection of drug resistance mutations in HIV pol. The current OLA involves four main steps from sample to analysis: (1) lysis and/or nucleic acid extraction, (2) amplification of HIV RNA or DNA, (3) ligation of oligonucleotide probes designed to detect single nucleotide mutations that confer HIV drug resistance, and (4) analysis via oligonucleotide surface capture, denaturation, and detection (CDD). The relative complexity of these steps has limited its adoption in resource-limited laboratories. Here we describe a simplification of the 2.5-hour plate-format CDD to a 45-minute paper-format CDD that eliminates the need for a plate reader. Analysis of mutations at four HIV-1 DR codons (K103N, Y181C, M184V, and G190A) in 26 blood specimens showed a strong correlation of the ratios of mutant signal to total signal between the paper CDD and the plate CDD. The assay described makes the OLA easier to perform in low resource laboratories.

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.

Single Genome Analysis for the Detection of Linked Multiclass Drug Resistance Mutations in HIV-1-Infected Children After Failure of Protease Inhibitor-Based First-Line Therapy

The WHO recommends protease inhibitor (PI)-based antiretroviral therapy (ART) for vertically infected children after failed nevirapine (NVP) prophylaxis. Emergence of PI resistance on the backdrop of preexisting non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance could compromise long-term treatment options in such children. We characterized multiclass drug resistance using single genome sequencing (SGS) in children with viremia while receiving PI-based ART. We applied SGS of HIV-1 protease (PR) and reverse transcriptase to longitudinal samples from a cohort of the Children with HIV Early Antiretroviral Therapy trial with viral loads >1000 copies per milliliter after 40 weeks of early ART. Bulk sequencing revealed NVP-selected resistance in 50% of these children, whereas SGS revealed NVP-selected resistance in 70%. Two children had baseline NRTI and PI mutations, suggesting previous maternal ART. Linked multiclass drug resistance after PI-based ART was detected by SGS in 2 of 10 children. In one child, the majority species contained M184V in reverse transcriptase linked to L10F, M46I/L, I54V, and V82A in PR and a triple-class drug-resistant variant with these mutations linked to the NNRTI mutation V108I. In the second child, the majority species contained M184V and V82A linked within viral genomes. We conclude that when PI-based ART is initiated soon after birth after single dose-NVP prophylaxis, PI and NRTI resistance can occur in the majority species as expected and also be selected on the same genomes as preexisting NNRTI-resistant mutations. These observations highlight a future therapeutic challenge for vertically infected children where antiretroviral drug classes are limited.

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.

Enhanced detection of viral diversity using partial and near full-length genomes of human immunodeficiency virus Type 1 provirus deep sequencing data from recently infected donors at four blood centers in Brazil

BACKGROUND

Here, we report application of high-throughput near full-length genome (NFLG) and partial human immunodeficiency virus Type 1 (HIV-1) proviral genome deep sequencing to characterize HIV in recently infected blood donors at four major blood centers in Brazil.

STUDY DESIGN AND METHODS

From 2007 to 2011, a total of 341 HIV+ blood donors from four blood centers were recruited to participate in a case-control study to identify HIV risk factors and motivations to donate. Forty-seven (17 from São Paulo, eight from Minas Gerais, 11 from Pernambuco, and 11 from Rio de Janeiro) were classified as recently infected based on testing by less-sensitive enzyme immunoassays. Five overlapping amplicons spanning the HIV genome were polymerase chain reaction amplified from peripheral blood mononuclear cells. The amplicons were molecularly barcoded, pooled, and sequenced by a paired-end protocol (Illumina).

RESULTS

Of the 47 recently infected donor samples studied, 39 (82.9%) NFLGs and six (12.7%) partial fragments were de novo assembled into contiguous sequences and successfully subtyped. Subtype B was the only nonrecombinant virus identified in this study and accounted for 62.2% (28/45) of samples. The remaining 37.8% (17/45) of samples showed various patterns of subtype discordance in different regions of HIV-1 genomes, indicating two to four circulating recombinant subtypes derived from Clades B, F, and C. Fourteen samples (31.1%) from this study harbored drug resistance mutations, indicating higher rate of drug resistance among Brazilian blood donors.

CONCLUSION

Our findings revealed a high proportion of HIV-1 recombinants among recently infected blood donors in Brazil, which has implications for future blood screening, diagnosis, therapy, and vaccine development.

Deep sequencing of HIV-1 near full-length proviral genomes identifies high rates of BF1 recombinants including two novel circulating recombinant forms (CRF) 70_BF1 and a disseminating 71_BF1 among blood donors in Pernambuco, Brazil

Background

The findings of frequent circulation of HIV-1 subclade F1 viruses and the scarcity of BF1 recombinant viruses based on pol subgenomic fragment sequencing among blood donors in Pernambuco (PE), Northeast of Brazil, were reported recently. Here, we aimed to determine whether the classification of these strains (n = 26) extends to the whole genome sequences.

Methods

Five overlapping amplicons spanning the HIV near full-length genomes (NFLGs) were PCR amplified from peripheral blood mononuclear cells (PBMCs) of 26 blood donors. The amplicons were molecularly bar-coded, pooled, and sequenced by Illumina paired-end protocol. The prevalence of viral variants containing drug resistant mutations (DRMs) was compared between plasma and PBMCs.

Results

Of the 26 samples studied, 20 NFLGs and 4 partial fragments were de novo assembled into contiguous sequences and successfully subtyped. Two distinct BF1 recombinant profiles designated CRF70_BF1 and CRF71_BF1, with 4 samples in profile I and 11 in profile II were detected and thus constitute two novel recombinant forms circulating in PE. Evidence of dual infections was detected in four patients co-infected with distinct HIV-1 subtypes. According to our estimate, the new CRF71_BF1 accounts for 10% of the HIV-1 circulating strains among blood donors in PE. Discordant data between the plasma and PBMCs-virus were found in 15 of 24 donors. Six of these strains displayed major DRMs only in PBMCs and four of which had detectable DRMs changes at prevalence between 1-20% of the sequenced population.

Conclusions

The high percentage of the new RF71_BF1 and other BF1 recombinants found among blood donors in Pernambuco, coupled with high rates of transmitted DRMs and dual infections confirm the need for effective surveillance to monitor the prevalence and distribution of HIV variants in a variety of settings in Brazil.

Minority drug-resistant HIV-1 variants in treatment naïve East-African and Caucasian patients detected by allele-specific real-time PCR

Objective

To assess the presence of two major non-nucleoside reverse transcriptase inhibitors (NNRTI) drug resistance mutations (DRMs), Y181C and K103N, in minor viral quasispecies of treatment naïve HIV-1 infected East-African and Swedish patients by allele-specific polymerase chain reaction (AS-PCR).

Methods

Treatment naïve adults (n = 191) with three epidemiological backgrounds were included: 92 Ethiopians living in Ethiopia; 55 East-Africans who had migrated to Sweden; and 44 Caucasians living in Sweden. The pol gene was analysed by standard population sequencing and by AS-PCR for the detection of Y181C and K103N.

Results

The Y181C was detected in the minority quasispecies of six Ethiopians (6.5%), in two Caucasians (4.5%), and in one East-African (1.8%). The K103N was detected in one East- African (1.8%), by both methods. The proportion of mutants ranged from 0.25% to 17.5%. Additional DRMs were found in all three treatment naïve patient groups by population sequencing.

Conclusions

Major NNRTI mutations can be found by AS-PCR in minor quasispecies of treatment naïve HIV-1 infected Ethiopians living in Ethiopia, in East-African and Caucasian patients living in Sweden in whom population sequencing reveal wild-type virus only. Surveys with standard sequencing are likely to underestimate transmitted drug resistance and the presence of resistant minor quasispecies in treatment naïve patients should be topic for future large scale studies.

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.

Decreased phenotypic susceptibility to etravirine in patients with predicted genotypic sensitivity

Background

A sensitive, phenotypic reverse transcriptase (RT)-based drug susceptibility assay for the detection of etravirine (ETR) resistance in patient isolates was developed and compared with the results from direct sequencing and ultra-deep pyrosequencing (UDPS).

Methods

Samples were obtained from 15 patients with antiretroviral therapy (ART) failure and from five non-nucleoside reverse transcriptase inhibitor (NNRTI)-naïve patients of whom four were infected by an NNRTI-resistant strain (transmitted drug resistance, TDR). In five patients, two consecutive samples (a and b) were taken for follow up of the virological response. HIV-1 RT was purified and drug susceptibility (IC₅₀) to ETR was estimated. Direct sequencing was performed in all samples and UDPS in samples from nine patients.

Results

Increased IC₅₀ to ETR was found in samples from 13 patients where direct sequencing predicted resistance in only four. UDPS identified additional (N = 11) NNRTI resistance associated mutations (RAMs) in six of nine tested patients. During early failure, IC₅₀ increases were observed in three of six patients without any ETR-RAMs detected by direct sequencing. In further two patients, who stopped NNRTI before sampling, increased IC₅₀ values were found shortly after, despite absence of ETR-RAMs. In two patients who had stopped NNRTI for >1 year, a concordance between phenotype and genotypes was found. Two patients with TDR had increased IC₅₀ despite no ETR-RAMs were detected by direct sequencing. UDPS revealed additional ETR-RAMs in four patients with a discrepancy between phenotype and direct sequencing.

Conclusions

The RT-based phenotypic assay showed decreased ETR susceptibility in patients where direct sequencing predicted ETR-sensitive virus. This increased phenotypic sensitivity was to a large extent supported by UDPS and treatment history. Our method could be valuable for further studies on the phenotypic kinetics of NNRTI resistance. The clinical relevance remains to be studied in larger patient-populations.

Theoretical and experimental assessment of degenerate primer tagging in ultra-deep applications of next-generation sequencing

Primer IDs (pIDs) are random oligonucleotide tags used in next-generation sequencing to identify sequences that originate from the same template. These tags are produced by degenerate primers during the reverse transcription of RNA molecules into cDNA. The use of pIDs helps to track the number of RNA molecules carried through amplification and sequencing, and allows resolution of inconsistencies between reads sharing a pID. Three potential issues complicate the above applications. First, multiple cDNAs may share a pID by chance; we found that while preventing any cDNAs from sharing a pID may be unfeasible, it is still practical to limit the number of these collisions. Secondly, a pID must be observed in at least three sequences to allow error correction; as such, pIDs observed only one or two times must be rejected. If the sequencing product contains copies from a high number of RT templates but produces few reads, our findings indicate that rejecting such pIDs will discard a great deal of data. Thirdly, the use of pIDs could influence amplification and sequencing. We examined the effects of several intrinsic and extrinsic factors on sequencing reads at both the individual and ensemble level.

Highly-sensitive allele-specific PCR testing identifies a greater prevalence of transmitted HIV drug resistance in Japan

BACKGROUND

The transmission of drug-resistant HIV in newly identified infected populations has become an underlying epidemic which can be better assessed with sensitive resistance testing. Since minority drug resistant variants cannot be detected by bulk sequencing, methods with improved sensitivity are required. Thus, the goal of this study was to evaluate if transmitted drug resistance mutations at minority levels in Japanese patients could be identified using highly sensitive allele-specific PCR (AS-PCR).

MATERIALS AND METHODS

Samples were taken from newly diagnosed HIV/AIDS cases at the National Nagoya Hospital from January 2008 to December 2009. All samples were bulk sequenced for HIV protease and reverse transcriptase. To detect minority populations with drug resistance, we used AS-PCR with mutation-specific primers designed for seven reverse transcriptase inhibitor resistance mutations, M41L, K65R, K70R, K103N, Y181C, M184V, and T215F/Y, and for three protease inhibitor resistance mutations, M46I/L and L90M.

RESULTS

We studied 149 newly identified HIV cases. Bulk sequencing detected 8 cases with NRTI resistance mutations (one with A62V, one D67E, one T215D, one T215E, two with T215L and two T215S) and 15 with PI resistance mutations (one with N88D and 14 with M46I). Results obtained by AS-PCR and bulk sequencing demonstrated good concordance but the AS-PCR enabled the detection of seven additional drug-resistant cases (one M41L, two with K65R, two with K70R, and one M184V) in the RT region. Additionally, AS-PCR assays identified 15 additional cases with M46I, five with M46L and four cases with L90M in the protease region.

CONCLUSIONS

Using AS-PCR substantially increased the detection of transmitted drug resistance in this population from 15.4% to 26.8%, further supporting the benefit of sensitive testing among drug-naïve populations. Since the clinical impact of minority drug-resistant populations is not fully comprehended for all mutations, follow-up studies are needed to understand their significance for treatment.

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.

Transmitted HIV drug resistance in treatment-naive Romanian patients

Transmitted HIV drug resistance (TDR) remains an important concern for individuals unexposed to antiretroviral treatment. Data on the prevalence of TDR, available mainly for HIV-1 subtype B, are now also emerging for other subtypes. In Romania, a steady predominance of subtype F was reported among both long-term survivor children and newly infected adults. The pol gene of 61 drug-naïve patients infected with HIV, diagnosed between 1997 and 2011 was sequenced in order to analyze the prevalence of primary resistance mutations and to correlate these with the infecting genotype. Only 5/61 specimens were classified as infected recently using the BED-Capture Enzyme Immunoassay. Subtype F1 was prevalent (80.3%), however, other HIV-1 clades are increasingly identified, especially in the group of subjects infected recently. An HIV transmission cluster, associated to injecting drug use was identified by phylogenetic analysis. The overall prevalence of TDR was 14.75%, mainly associated with NRTI resistance (13.11%), TAMs and M184V being the most common mutations. A declining trend of TDR was recorded from 26.08% in 1997-2004 to 7.89% in 2005-2011. No primary resistance was identified among recent seroconvertors. All HIV-1 strains had minor mutations in the protease and RT genes, often detected at polymorphic positions. The declining rates of TDR might be related to the high efficacy of HAART and to the increasing number of treated patients with virological success who have a low risk of transmission. The recent increase of HIV-1 infections which involve other subtypes impose a continuous surveillance of the genetic composition of the epidemic.

Characterizing the emergence and persistence of drug resistant mutations in HIV-1 subtype C infections using 454 ultra deep pyrosequencing

BACKGROUND

The role of HIV-1 RNA in the emergence of resistance to antiretroviral therapies (ARTs) is well documented while less is known about the role of historical viruses stored in the proviral DNA. The primary focus of this work was to characterize the genetic diversity and evolution of HIV drug resistant variants in an individual's provirus during antiretroviral therapy using next generation sequencing.

METHODS

Blood samples were collected prior to antiretroviral therapy exposure and during the course of treatment from five patients in whom drug resistance mutations had previously been identified using consensus sequencing. The spectrum of viral variants present in the provirus at each sampling time-point were characterized using 454 pyrosequencing from multiple combined PCR products. The prevalence of viral variants containing drug resistant mutations (DRMs) was characterized at each time-point.

RESULTS

Low abundance drug resistant viruses were identified in 14 of 15 sampling time-points from the five patients. In all individuals DRMs against current therapy were identified at one or more of the sampling time-points. In two of the five individuals studied these DRMs were present prior to treatment exposure and were present at high prevalence within the amplified and sequenced viral population. DRMs to drugs other than those being currently used were identified in four of the five individuals.

CONCLUSION

The presence of DRMs in the provirus, regardless of their observed prevalence did not appear to have an effect on clinical outcomes in the short term suggesting that the drug resistant viral variants present in the proviral DNA do not appear to play a role in the short term in facilitating the emergence of drug resistance.

A combined screening platform for HIV treatment failure and resistance

Background

To develop a low cost method to screen for virologic failure of antiretroviral therapy (ART) and HIV-1 drug resistance, we performed a retrospective evaluation of a screening assay using serial dilutions of HIV-1 RNA-spiked blood plasma and samples from patients receiving >6 months of first-line ART.

Methods

Serial dilution testing was used to assess sensitivity of a simple PCR-based assay (targeted at ≥1,000 HIV RNA copies/mL). We created blood plasma minipools of five samples, extracted HIV RNA from the pools, PCR amplified the reverse transcriptase (RT) coding region of the HIV-1 pol gene from extracted RNA, sequenced PCR product of positive pools, and used sequences to determine drug resistance. Sensitivity, specificity, and predictive values were determined for different levels of virologic failure based on maximum viral loads of individual samples within a pool.

Results

Of 295 samples analyzed, 43 (15%) had virologic failure at ≥50 copies/mL (range 50–10,500 copies/mL, four at ≥1,000 copies/mL). The assay demonstrated 100% sensitivity to detect virus from these four samples, requiring only one round of PCR, and 56% and 89% sensitivity to detect samples with ≥50 and ≥500 copies/mL using two rounds. Amplified PCR products of all positive pools were successfully sequenced and 30% harbored ≥1 major resistance mutation. This method would have cost 10% of the combined costs of individual viral load and resistance testing.

Conclusions

We present a novel method that can screen for both virologic failure of first-line ART and drug resistance. The method is much less expensive than current methods, which may offer sustainability in resource-limited settings.

Highly sensitive and specific detection of rare variants in mixed viral populations from massively parallel sequence data

Viruses diversify over time within hosts, often undercutting the effectiveness of host defenses and therapeutic interventions. To design successful vaccines and therapeutics, it is critical to better understand viral diversification, including comprehensively characterizing the genetic variants in viral intra-host populations and modeling changes from transmission through the course of infection. Massively parallel sequencing technologies can overcome the cost constraints of older sequencing methods and obtain the high sequence coverage needed to detect rare genetic variants (<1%) within an infected host, and to assay variants without prior knowledge. Critical to interpreting deep sequence data sets is the ability to distinguish biological variants from process errors with high sensitivity and specificity. To address this challenge, we describe V-Phaser, an algorithm able to recognize rare biological variants in mixed populations. V-Phaser uses covariation (i.e. phasing) between observed variants to increase sensitivity and an expectation maximization algorithm that iteratively recalibrates base quality scores to increase specificity. Overall, V-Phaser achieved >97% sensitivity and >97% specificity on control read sets. On data derived from a patient after four years of HIV-1 infection, V-Phaser detected 2,015 variants across the ∼10 kb genome, including 603 rare variants (<1% frequency) detected only using phase information. V-Phaser identified variants at frequencies down to 0.2%, comparable to the detection threshold of allele-specific PCR, a method that requires prior knowledge of the variants. The high sensitivity and specificity of V-Phaser enables identifying and tracking changes in low frequency variants in mixed populations such as RNA viruses.

New sequencing technologies provide unprecedented resolution to study pathogen populations, such as the single stranded RNA viruses HIV, dengue (DENV), and West Nile (WNV), and how they evolve within infected individuals in response to immune, therapeutic, and vaccine pressures. While these new technologies provide high volumes of data, these data contain process errors. To detect biological variants, especially those occurring at low frequencies in the population, these technologies require a method to differentiate biological variants from process errors with high sensitivity and specificity. To address this challenge, we introduce the V-Phaser algorithm, which distinguished the covariation of biological variants from that of process errors. We validate the method by measuring how frequently it correctly identifies variants and errors on actual read sets with known variation. Further, using data derived from a patient following four years of HIV-1 infection, we show that V-Phaser can detect biological variants at frequencies comparable to approaches that require prior knowledge. V-Phaser is available for download at: http://www.broadinstitute.org/scientific-community/software.

Polymorphic mutations associated with the emergence of the multinucleoside/tide resistance mutations 69 insertion and Q151M

BACKGROUND

We hypothesized that polymorphic mutations exist that are associated with the emergence of the multinucleoside resistance mutations (MNR), 69 insertion and Q151M.

METHODS

The Swiss HIV Cohort Study was screened, and the frequencies of polymorphic mutations in HIV-1 (subtype B) were compared between patients detected with the 69 insertion (n = 17), Q151M (n = 29), ≥2 thymidine analogue mutations (TAM) 1 (n = 400) or ≥2 TAM 2 (n = 249). Logistic regressions adjusted for the antiretroviral treatment history were performed to analyze the association of the polymorphic mutations with MNR.

RESULTS

The 69 insertion and TAM 1 were strongly associated and occurred in 94.1% (16 of 17) together. The 69 insertion seemed to emerge as a consequence of the TAM 1 pathway (median years until detection: 6.8 compared with 4.4 for ≥2 TAM 1, P Wilcoxon = 0.009). Frequencies of 8 polymorphic mutations (K43E, V60I, S68G, S162C, T165I, I202V, R211K, F214L) were significantly different between groups. Logistic regression showed that F214L and V60I were associated with the emergence of Q151M/TAM 2 opposed to 69 insertion/TAM 1. S68G, T165I, and I202V were associated with Q151M instead of TAM 2.

CONCLUSIONS

Besides antiretroviral therapy, polymorphic mutations may contribute to the emergence of specific MNR mutations.

Ultra-deep sequencing of HIV-1 reverse transcriptase before start of an NNRTI-based regimen in treatment-naive patients

There are conflicting data on the impact of low frequency HIV-1 drug-resistant mutants on the response of first-line highly active antiretroviral therapy (HAART), more specifically containing a NNRTI. As population sequencing does not detect resistant viruses representing less than 15-25% of the viral population, more sensitive techniques have been developed but still need clinical validation. We evaluated ultra-deep sequencing (UDPS), recently more available and affordable, as a tool for the detection of HIV-1 minority species carrying drug resistant mutation (DRM) in a clinical setting. A retrospective analysis of the reverse transcriptase (RT) gene of plasma HIV-1 from 70 patients starting a NNRTI based regimen was performed. Minority populations were defined as representing > 1% and < 20% of the total viral population. Using UDPS, we could not confirm an association between the presence of low minority variants harbouring RT mutations at the start of therapy and primary or secondary therapeutic failure.

Accurate sampling and deep sequencing of the HIV-1 protease gene using a Primer ID

Viruses can create complex genetic populations within a host, and deep sequencing technologies allow extensive sampling of these populations. Limitations of these technologies, however, potentially bias this sampling, particularly when a PCR step precedes the sequencing protocol. Typically, an unknown number of templates are used in initiating the PCR amplification, and this can lead to unrecognized sequence resampling creating apparent homogeneity; also, PCR-mediated recombination can disrupt linkage, and differential amplification can skew allele frequency. Finally, misincorporation of nucleotides during PCR and errors during the sequencing protocol can inflate diversity. We have solved these problems by including a random sequence tag in the initial primer such that each template receives a unique Primer ID. After sequencing, repeated identification of a Primer ID reveals sequence resampling. These resampled sequences are then used to create an accurate consensus sequence for each template, correcting for recombination, allelic skewing, and misincorporation/sequencing errors. The resulting population of consensus sequences directly represents the initial sampled templates. We applied this approach to the HIV-1 protease (pro) gene to view the distribution of sequence variation of a complex viral population within a host. We identified major and minor polymorphisms at coding and noncoding positions. In addition, we observed dynamic genetic changes within the population during intermittent drug exposure, including the emergence of multiple resistant alleles. These results provide an unprecedented view of a complex viral population in the absence of PCR resampling.

HIV drug resistance (HIVDR) in antiretroviral therapy-naïve patients in Tanzania not eligible for WHO threshold HIVDR survey is dramatically high

Background

The World Health Organization (WHO) has recommended guidelines for a HIV drug resistance (HIVDR) survey for resource-limited countries. Eligibility criteria for patients include age below 25 years in order to focus on the prevalence of transmitted HIVDR (tHIVDR) in newly-infected individuals. Most of the participating sites across Africa have so far reported tHIVDR prevalences of below 5%. In this study we investigated whether the rate of HIVDR in patients <25 years is representative for HIVDR in the rest of the therapy-naïve population.

Methods and Findings

HIVDR was determined in 88 sequentially enrolled ART-naïve patients from Mwanza, Tanzania (mean age 35.4 years). Twenty patients were aged <25 years and 68 patients were aged 25–63 years. The frequency of HIVDR in the study population was 14.8% (95%; CI 0.072–0.223) and independent of NVP-resistance induced by prevention of mother-to-child transmission programs. Patients >25 years had a significantly higher HIVDR frequency than younger patients (19.1%; 95% CI 0.095–0.28) versus 0%, P = 0.0344). In 2 out of the 16 patients with HIVDR we found traces of antiretrovirals (ARVs) in plasma.

Conclusions

ART-naïve patients aged over 25 years exhibited significantly higher HIVDR than younger patients. Detection of traces of ARVs in individuals with HIVDR suggests that besides transmission, undisclosed misuse of ARVs may constitute a significant factor in the generation of the observed high HIVDR rate. The current WHO tHIVDR survey that is solely focused on the transmission of HIVDR and that excludes patients over 25 years of age may therefore result in substantial underestimation of the prevalence of HIVDR in the therapy-naïve population. Similar studies should be performed also in other areas to test whether the so far reported optimistic picture of low HIVDR prevalence in young individuals is really representative for the rest of the ART-naïve HIV-infected population.

Synthetic long oligonucleotides to generate artificial templates for use as positive controls in molecular assays: drug resistance mutations in influenza virus as an example

BACKGROUND

Positive controls are an integral component of any sensitive molecular diagnostic tool, but this can be affected, if several mutations are being screened in a scenario of a pandemic or newly emerging disease where it can be difficult to acquire all the necessary positive controls from the host. This work describes the development of a synthetic oligo-cassette for positive controls for accurate and highly sensitive diagnosis of several mutations relevant to influenza virus drug resistance.

RESULTS

Using influenza antiviral drug resistance mutations as an example by employing the utility of synthetic paired long oligonucleotides containing complementary sequences at their 3' ends and utilizing the formation of oligonucleotide dimers and DNA polymerization, we generated ~170bp dsDNA containing several known specific neuraminidase inhibitor (NAI) resistance mutations. These templates were further cloned and successfully applied as positive controls in downstream assays.

CONCLUSION

This approach significantly improved the development of diagnosis of resistance mutations in terms of time, accuracy, efficiency and sensitivity, which are paramount to monitoring the emergence and spread of antiviral drug resistant influenza strains. Thus, this may have a significantly broader application in molecular diagnostics along with its application in rapid molecular testing of all relevant mutations in an event of pandemic.

Reappearance of minority K103N HIV-1 variants after interruption of ART initiated during primary HIV-1 infection

BACKGROUND

In the Zurich Primary HIV infection study (ZPHI), minority drug-resistant HIV-1 variants were detected in some acutely HIV-1-infected patients prior to initiation of early antiretroviral therapy (ART). Here, we investigated the reappearance of minority K103N and M184V HIV-1 variants in these patients who interrupted efficient early ART after 8-27 months according to the study protocol. These mutations are key mutations conferring drug resistance to reverse transcriptase inhibitors and they belong to the most commonly transmitted drug resistance mutations.

METHODOLOGY/PRINCIPAL FINDINGS

Early ART was offered to acutely HIV-1-infected patients enrolled in the longitudinal prospective ZPHI study. Six patients harboring and eleven patients not harboring drug-resistant viruses at low frequencies prior to ART were included in this substudy. Minority K103N and M184V HIV-1 variants were quantified in longitudinal plasma samples after treatment interruption by allele-specific real-time PCR. All 17 patients were infected with HIV-1 subtype B between 04/2003 and 09/2005 and received LPV/r+AZT+3TC during primary HIV-1 infection (PHI). Minority K103N HIV-1 variants reappeared after cessation of ART in two of four patients harboring this variant during PHI and even persisted in one of those patients at frequencies similar to the frequency observed prior to ART (<1%). The K103N mutation did not appear during treatment interruption in any other patient. Minority M184V HIV-1 variants were detected in two patients after ART interruption, one harboring and one not harboring these variants prior to ART.

CONCLUSION

Minority K103N HIV-1 variants, present in acutely HIV-1 infected patients prior to early ART, can reappear and persist after interruption of suppressive ART containing two nucleoside/nucleotide analogue reverse transcriptase inhibitors and a ritonavir-boosted protease inhibitor.

TRIAL REGISTRATION

Clinicaltrials.gov NCT00537966.

Within-host co-evolution of Gag P453L and protease D30N/N88D demonstrates virological advantage in a highly protease inhibitor-exposed HIV-1 case

To better understand the mechanism of HIV group-specific antigen (Gag) and protease (PR) co-evolution in drug-resistance acquisition, we analyzed a drug-resistance case by both bioinformatics and virological methods. We especially considered the quality of sequence data and analytical accuracy by introducing single-genome sequencing (SGS) and Spidermonkey/Bayesian graphical models (BGM) analysis, respectively. We analyzed 129 HIV-1 Gag-PR linkage sequences obtained from 8 time points, and the resulting sequences were applied to the Spidermonkey co-evolution analysis program, which identified ten mutation pairs as significantly co-evolving. Among these, we focused on associations between Gag-P453L, the P5' position of the p1/p6 cleavage-site mutation, and PR-D30N/N88D nelfinavir-resistant mutations, and attempted to clarify their virological significance in vitro by constructing recombinant clones. The results showed that P453L(Gag) has the potential to improve replication capacity and the Gag processing efficiency of viruses with D30N(PR)/N88D(PR) but has little effect on nelfinavir susceptibility. Homology modeling analysis suggested that hydrogen bonds between the 30th PR residue and the R452Gag are disturbed by the D30N(PR) mutation, but the impaired interaction is compensated by P453L(Gag) generating new hydrophobic interactions. Furthermore, database analysis indicated that the P453L(Gag)/D30N(PR)/N88D(PR) association was not specific only to our clinical case, but was common among AIDS patients.

Replicative phenotyping adds value to genotypic resistance testing in heavily pre-treated HIV-infected individuals--the Swiss HIV Cohort Study

Background

Replicative phenotypic HIV resistance testing (rPRT) uses recombinant infectious virus to measure viral replication in the presence of antiretroviral drugs. Due to its high sensitivity of detection of viral minorities and its dissecting power for complex viral resistance patterns and mixed virus populations rPRT might help to improve HIV resistance diagnostics, particularly for patients with multiple drug failures. The aim was to investigate whether the addition of rPRT to genotypic resistance testing (GRT) compared to GRT alone is beneficial for obtaining a virological response in heavily pre-treated HIV-infected patients.

Methods

Patients with resistance tests between 2002 and 2006 were followed within the Swiss HIV Cohort Study (SHCS). We assessed patients' virological success after their antiretroviral therapy was switched following resistance testing. Multilevel logistic regression models with SHCS centre as a random effect were used to investigate the association between the type of resistance test and virological response (HIV-1 RNA <50 copies/mL or ≥1.5log reduction).

Results

Of 1158 individuals with resistance tests 221 with GRT+rPRT and 937 with GRT were eligible for analysis. Overall virological response rates were 85.1% for GRT+rPRT and 81.4% for GRT. In the subgroup of patients with >2 previous failures, the odds ratio (OR) for virological response of GRT+rPRT compared to GRT was 1.45 (95% CI 1.00-2.09). Multivariate analyses indicate a significant improvement with GRT+rPRT compared to GRT alone (OR 1.68, 95% CI 1.31-2.15).

Conclusions

In heavily pre-treated patients rPRT-based resistance information adds benefit, contributing to a higher rate of treatment success.

Indinavir resistance evolution in one human immunodeficiency virus type 1 infected patient revealed by single-genome amplification

Human Immunodeficiency Virus Type 1 exists in vivo as quasispecies, and one of the genome's characteristics is its diversity. During the antiretroviral therapy, drug resistance is the main obstacle to effective viral prevention. Understanding the molecular evolution process is fundamental to analyze the mechanism of drug resistance and develop a strategy to minimize resistance.

OBJECTIVE

The molecular evolution of drug resistance of one patient who had received reverse transcriptase inhibitors for a long time and had treatment which replaced Nevirapine with Indinavir was analyzed, with the aim of observing the drug resistance evolution pathway.

METHODS

The patient, XLF, was followed-up for six successive times. The viral populations were amplified and sequenced by single-genome amplification. All the sequences were submitted to the Stanford HIV Drug Resistance Database for the analysis of genotypic drug resistance.

RESULTS

149 entire protease and 171 entire reverse transcriptase sequences were obtained from these samples, and all sequences were identified as subtype B. Before the patient received Indinavir, the viral population only had some polymorphisms in the protease sequences. After the patient began Indinavir treatment, the variants carrying polymorphisms declined while variants carrying the secondary mutation G73S gained the advantage. As therapy was prolonged, G73S was combined with M46I/L90M to form a resistance pattern M46I/G73S/L90M, which then became the dominant population. 97.9% of variants had the M46I/G73S/L90M pattern at XLF6. During the emergence of protease inhibitors resistance, reverse transcriptase inhibitors resistance maintained high levels.

CONCLUSION

Indinavir-resistance evolution was observed by single-genome amplification. During the course of changing the regimen to incorporate Indinavir, the G73S mutation occurred and was combined with M46I/L90M.

Minority variants of drug-resistant HIV

Minor drug-resistant variants exist in every patient infected with human immunodeficiency virus (HIV). Because these minority variants are usually present at very low levels, they cannot be detected and quantified using conventional genotypic and phenotypic tests. Recently, several assays have been developed to characterize these low-abundance drug-resistant variants in the large, genetically complex population that is present in every HIV-infected individual. The most important issue is what results generated by these assays can predict clinical or treatment outcomes and might guide the management of patients in clinical practice. Cutoff values for the detection of these low-abundance viral variants that predict an increased risk of treatment failure should be determined. These thresholds may be specific for each mutation and treatment regimen. In this review, we summarize the attributes and limitations of the currently available detection assays and review the existing information about both acquired and transmitted drug-resistant minority variants.

Secondary integrase resistance mutations found in HIV-1 minority quasispecies in integrase therapy-naive patients have little or no effect on susceptibility to integrase inhibitors

The goal of this study was to explore the presence of integrase strand transfer inhibitor (InSTI) resistance mutations in HIV-1 quasispecies present in InSTI-naïve patients and to evaluate their in vitro effects on phenotypic susceptibility to InSTIs and their replication capacities. The RT-RNase H-IN region was PCR amplified from plasma viral RNA obtained from 49 HIV-1 subtype B-infected patients (21 drug naïve and 28 failing highly active antiretroviral therapy [HAART] not containing InSTIs) and recombined with an HXB2-based backbone with RT and IN deleted. Recombinant viruses were tested against raltegravir and elvitegravir and for replication capacity. Three-hundred forty-four recombinant viruses from 49 patients were successfully analyzed both phenotypically and genotypically. The majority of clones were not phenotypically resistant to InSTIs: 0/344 clones showed raltegravir resistance, and only 3 (0.87%) showed low-level elvitegravir resistance. No primary resistance mutations for raltegravir and elvitegravir were found as major or minor species. The majority of secondary mutations were also absent or rarely present. Secondary mutations, such as T97A and G140S, found rarely and only as minority quasispecies, were present in the elvitegravir-resistant clones. A novel mutation, E92G, although rarely found in minority quasispecies, showed elvitegravir resistance. Preexisting genotypic and phenotypic raltegravir resistance was extremely rare in InSTI-naïve patients and confined to only a restricted minority of secondary variants. Overall, these results, together with others based on population and ultradeep sequencing, suggest that at this point IN genotyping in all patients before raltegravir treatment may not be cost-effective and should not be recommended until evidence of transmitted drug resistance to InSTIs or the clinical relevance of IN minor variants/polymorphisms is determined.

High frequency of integrase Q148R minority variants in HIV-infected patients naive of integrase inhibitors

BACKGROUND

Integrase positions 148 and 155 represent main determinants of resistance to integrase inhibitors. We assessed the prevalence of minority variants harboring such mutations in integrase-naive HIV-infected patients.

METHODS

Two groups of patients were studied: 40 heavily antiretroviral-experienced patients, initiating a raltegravir-based therapy and 51 antiretroviral-naive patients. Allele-specific real-time PCR (AS-PCR) systems, developed for Q148H, Q148R and N155H mutations, were performed at baseline for antiretroviral-experienced patients. Samples from antiretroviral-naive patients were tested with the Q148R AS-PCR assay.

RESULTS

The limits of detection of AS-PCR systems were 0.10, 0.10 and 0.05% for Q148H, Q148R and N155H mutations, respectively. AS-PCR systems were successful in 79 of 91 samples. In antiretroviral-experienced patients, Q148R minority variants were frequently detected (26/32 patients, 81%) at low-level frequency (median = 0.40%), whereas no minority variants exhibiting Q148H or N155H mutation were found. Twenty-four of 26 patients exhibiting Q148R variants were virological responders but four of them displayed a delayed virological response occurring between W18 and W36. Two patients exhibited virological failure under raltegravir, both harboring Q148R minority variants at baseline. However, we did not find any association between the presence of Q148R minority variants and an increased risk of virological failure. Q148R minority variants were also found in 86% of antiretroviral-naive patients, a prevalence significantly higher than that of K103N minority variants (26%).

CONCLUSION

Q148R variants were frequently detected, always at low-level, in antiretroviral-experienced and naive patients. Although their presence was not consistently associated with virological failure, their impact on long-term viral suppression needs to be further investigated.

Low-frequency K103N strengthens the impact of transmitted drug resistance on virologic responses to first-line efavirenz or nevirapine-based highly active antiretroviral therapy

BACKGROUND

There are conflicting data on the impact of low-frequency transmitted drug-resistant mutants on responses to first-line highly active antiretroviral therapy (HAART).

METHODS

Patients started nevirapine or efavirenz with two or more nucleoside/nucleotide reverse transcriptase inhibitors in 1998-2007 without a prior resistance test at a median 1.0 (interquartile range, 0.0-3.4) year after diagnosis and with a median 218 (interquartile range, 131-296) CD4 cells/mm3, and had at least 24 weeks of follow up. Pre-HAART plasma samples were tested retrospectively by bulk genotyping and sensitive real-time polymerase chain reaction targeting reverse transcriptase K65R, K103N, Y181C, M184V, and G190A (interpretative cutoff 0.3%-0.9%).

RESULTS

Among 93 patients, seven of 18 who experienced virologic failure and zero of 75 who maintained virologic suppression showed pre-HAART resistance, including three with high-frequency mutations detectable by bulk genotyping (two K103N, one G190A) and four with low-frequency K103N detectable only by polymerase chain reaction. Detection of either bulk (P = 0.006) or low-frequency (P = 0.001) resistance was significantly associated with the odds of virologic failure; combining the two markedly increased the strength of the association (P < 0.0001). At failure, the pre-HAART mutations were detected by bulk genotyping in five of seven patients alongside additional reverse transcriptase mutations.

CONCLUSIONS

Low-frequency K103N mutants were as prevalent as bulk-detectable variants before starting HAART. Both high- and low-frequency mutants were significantly associated with virologic failure.

The efficiency of single genome amplification and sequencing is improved by quantitation and use of a bioinformatics tool

Typically, population-based sequencing of HIV does not detect minority variants present at levels below 20-30%. Single genome amplification (SGA) and sequencing improves detection, but it requires many PCRs to find the optimal terminal dilution to use. A novel method for guiding the selection of a terminal dilution was developed and compared to standard methods. A quantitative real-time PCR (qRT-PCR) protocol was developed. HIV RNA was extracted, reverse transcribed, and quantitated. A bioinformatics web-based application was created for calculating the optimal concentration of cDNA to use based on results of a trial PCR using the dilution suggested by the qRT-PCR results. This method was compared to the standard. Using the standard protocol, the mean number of PCRs giving an average of 30 (26-34, SD=3) SGA per sample was 245 (218-266, SD=20) after an average of 8 trial dilutions. Using this method, 135 PCRs (135-135, SD=0) produced 30 (27-30, SD=1) SGA using exactly two dilutions. This new method reduced turnaround time from 8 to 2 days. Standard methods of SGA sequencing can be costly and both time- and labor-intensive. By choosing a terminal dilution concentration with the proposed method, the number of PCRs required is decreased and efficiency improved.

Sensitive detection of the K103N non-nucleoside reverse transcriptase inhibitor resistance mutation in treatment-naïve HIV-1 infected individuals by rolling circle amplification

Primary or transmitted antiretroviral drug resistance mutations pose a significant obstacle for optimizing antiviral treatment. When present at low-levels, resistance mutations are less likely to be detected by standard genotyping assays. This study utilizes a novel rolling circle amplification (RCA) method using padlock probes to achieve the sensitive, specific and low-level detection of the NNRTI resistance K103N from 59 HIV+ treatment-naïve patients from Beijing, China. Using standard genotyping methods, primary drug resistance mutations to either protease or RT inhibitors were found in 25% (15/59) of patients attending hospital clinics in Beijing. Among these 15 patients with antiretroviral (ARV) resistance mutations, standard sequence-based genotyping revealed that most (10/15) had the 103N. Using a highly sensitive RCA assay, 5 more patients among the 59 treatment-naïve cohort were found to have the 103N, but at low-levels, leading to an overall rate of 103N at 25.4% (15/59) in this population. The high prevalence of the 103N suggests that baseline resistance testing should be performed before treatment in this population. Importantly, the new RCA technology allows large-scale, sensitive detection of drug resistance mutations, including detection of minority populations with minimal equipment requirement.