Low-Frequency Drug Resistance in HIV-Infected Ugandans on Antiretroviral Treatment Is Associated with Regimen Failure

Rate of response to initial antiretroviral therapy according to level of pre-existing HIV-1 drug resistance detected by next-generation sequencing in the strategic timing of antiretroviral treatment (START) study

OBJECTIVES

The main objective of this analysis was to evaluate the impact of pre-existing drug resistance by next-generation sequencing (NGS) on the risk of treatment failure (TF) of first-line regimens in participants enrolled in the START study.

METHODS

Stored plasma from participants with entry HIV RNA >1000 copies/mL were analysed using NGS (llumina MiSeq). Pre-existing drug resistance was defined using the mutations considered by the Stanford HIV Drug Resistance Database (HIVDB v8.6) to calculate the genotypic susceptibility score (GSS, estimating the number of active drugs) for the first-line regimen at the detection threshold windows of >20%, >5%, and >2% of the viral population. Survival analysis was conducted to evaluate the association between the GSS and risk of TF (viral load >200 copies/mL plus treatment change).

RESULTS

Baseline NGS data were available for 1380 antiretroviral therapy (ART)-naïve participants enrolled over 2009-2013. First-line ART included a non-nucleoside reverse transcriptase inhibitor (NNRTI) in 976 (71%), a boosted protease inhibitor in 297 (22%), or an integrase strand transfer inhibitor in 107 (8%). The proportions of participants with GSS <3 were 7% for >20%, 10% for >5%, and 17% for the >2% thresholds, respectively. The adjusted hazard ratio of TF associated with a GSS of 0-2.75 versus 3 in the subset of participants with mutations detected at the >2% threshold was 1.66 (95% confidence interval 1.01-2.74; p = 0.05) and 2.32 (95% confidence interval 1.32-4.09; p = 0.003) after restricting the analysis to participants who started an NNRTI-based regimen.

CONCLUSIONS

Up to 17% of participants initiated ART with a GSS <3 on the basis of NGS data. Minority variants were predictive of TF, especially for participants starting NNRTI-based regimens.

Long-Read Sequencing with Hierarchical Clustering for Antiretroviral Resistance Profiling of Mixed Human Immunodeficiency Virus Quasispecies

BACKGROUND

HIV infections often develop drug resistance mutations (DRMs), which can increase the risk of virological failure. However, it has been difficult to determine if minor mutations occur in the same genome or in different virions using Sanger sequencing and short-read sequencing methods. Oxford Nanopore Technologies (ONT) sequencing may improve antiretroviral resistance profiling by allowing for long-read clustering.

METHODS

A new ONT sequencing-based method for profiling DRMs in HIV quasispecies was developed and validated. The method used hierarchical clustering of long amplicons that cover regions associated with different types of antiretroviral drugs. A gradient series of an HIV plasmid and 2 plasma samples was prepared to validate the clustering performance. The ONT results were compared to those obtained with Sanger sequencing and Illumina sequencing in 77 HIV-positive plasma samples to evaluate the diagnostic performance.

RESULTS

In the validation study, the abundance of detected quasispecies was concordant with the predicted result with the R2 of > 0.99. During the diagnostic evaluation, 59/77 samples were successfully sequenced for DRMs. Among 18 failed samples, 17 were below the limit of detection of 303.9 copies/μL. Based on the receiver operating characteristic analysis, the ONT workflow achieved an F1 score of 0.96 with a cutoff of 0.4 variant allele frequency. Four cases were found to have quasispecies with DRMs, in which 2 harbored quasispecies with more than one class of DRMs. Treatment modifications were recommended for these cases.

CONCLUSIONS

Long-read sequencing coupled with hierarchical clustering could differentiate the quasispecies resistance profiles in HIV-infected samples, providing a clearer picture for medical care.

HIV-1 Low-Frequency Variants Identified in Antiretroviral-Naïve Subjects with Virologic Failure after 12 Months of Follow-Up in Panama

Low-frequency mutations associated with drug resistance have been related to virologic failure in subjects with no history of pre-treatment and recent HIV diagnosis. In total, 78 antiretroviral treatment (ART)-naïve subjects with a recent HIV diagnosis were selected and followed by CD4+ T lymphocytes and viral load tests to detect virologic failure. We sequenced the basal samples retrospectively using next-generation sequencing (NGS), looking for low-frequency mutations that had not been detected before using the Sanger sequencing method (SSM) and describing the response to ART. Twenty-two subjects developed virologic failure (VF), and thirteen of them had at least one drug-resistance mutation associated with Reverse Transcriptase Inhibitors (RTI) and Protease Inhibitors (PIs) at frequency levels ≤ 1%, not detected previously in their basal genotyping test. No resistance mutations were observed to Integrase Strand Transfer Inhibitors (INSTIs). We identified a possible cause of VF in ART-naïve subjects with low-frequency mutations detected. To our knowledge, this is the first evaluation of pre-existing drug resistance for HIV-1 minority variants carried out on ART-naïve people living with HIV/AIDS (PLWHA) by analyzing the HIV-1 pol gene using NGS in the country.

TaME-seq2: tagmentation-assisted multiplex PCR enrichment sequencing for viral genomic profiling

BACKGROUND

Previously developed TaME-seq method for deep sequencing of HPV, allowed simultaneous identification of the human papillomavirus (HPV) DNA consensus sequence, low-frequency variable sites, and chromosomal integration events. The method has been successfully validated and applied to the study of five carcinogenic high-risk (HR) HPV types (HPV16, 18, 31, 33, and 45). Here, we present TaME-seq2 with an updated laboratory workflow and bioinformatics pipeline. The HR-HPV type repertoire was expanded with HPV51, 52, and 59. As a proof-of-concept, TaME-seq2 was applied on SARS-CoV-2 positive samples showing the method's flexibility to a broader range of viruses, both DNA and RNA.

RESULTS

Compared to TaME-seq version 1, the bioinformatics pipeline of TaME-seq2 is approximately 40× faster. In total, 23 HPV-positive samples and seven SARS-CoV-2 clinical samples passed the threshold of 300× mean depth and were submitted to further analysis. The mean number of variable sites per 1 kb was ~ 1.5× higher in SARS-CoV-2 than in HPV-positive samples. Reproducibility and repeatability of the method were tested on a subset of samples. A viral integration breakpoint followed by a partial genomic deletion was found in within-run replicates of HPV59-positive sample. Identified viral consensus sequence in two separate runs was > 99.9% identical between replicates, differing by a couple of nucleotides identified in only one of the replicates. Conversely, the number of identical minor nucleotide variants (MNVs) differed greatly between replicates, probably caused by PCR-introduced bias. The total number of detected MNVs, calculated gene variability and mutational signature analysis, were unaffected by the sequencing run.

CONCLUSION

TaME-seq2 proved well suited for consensus sequence identification, and the detection of low-frequency viral genome variation and viral-chromosomal integrations. The repertoire of TaME-seq2 now encompasses seven HR-HPV types. Our goal is to further include all HR-HPV types in the TaME-seq2 repertoire. Moreover, with a minor modification of previously developed primers, the same method was successfully applied for the analysis of SARS-CoV-2 positive samples, implying the ease of adapting TaME-seq2 to other viruses.

Frequency matters: comparison of drug resistance mutation detection by Sanger and next-generation sequencing in HIV-1

BACKGROUND

Next-generation sequencing (NGS) is gradually replacing Sanger sequencing (SS) as the primary method for HIV genotypic resistance testing. However, there are limited systematic data on comparability of these methods in a clinical setting for the presence of low-abundance drug resistance mutations (DRMs) and their dependency on the variant-calling thresholds.

METHODS

To compare the HIV-DRMs detected by SS and NGS, we included participants enrolled in the Swiss HIV Cohort Study (SHCS) with SS and NGS sequences available with sample collection dates ≤7 days apart. We tested for the presence of HIV-DRMs and compared the agreement between SS and NGS at different variant-calling thresholds.

RESULTS

We included 594 pairs of SS and NGS from 527 SHCS participants. Males accounted for 80.5% of the participants, 76.3% were ART naive at sample collection and 78.1% of the sequences were subtype B. Overall, we observed a good agreement (Cohen's kappa >0.80) for HIV-DRMs for variant-calling thresholds ≥5%. We observed an increase in low-abundance HIV-DRMs detected at lower thresholds [28/417 (6.7%) at 10%-25% to 293/812 (36.1%) at 1%-2% threshold]. However, such low-abundance HIV-DRMs were overrepresented in ART-naive participants and were in most cases not detected in previously sampled sequences suggesting high sequencing error for thresholds <3%.

CONCLUSIONS

We found high concordance between SS and NGS but also a substantial number of low-abundance HIV-DRMs detected only by NGS at lower variant-calling thresholds. Our findings suggest that a substantial fraction of the low-abundance HIV-DRMs detected at thresholds <3% may represent sequencing errors and hence should not be overinterpreted in clinical practice.

HIV-1 Genotypic Resistance Testing Using Sanger and Next-Generation Sequencing in Adults with Low-Level Viremia in China

OBJECTIVE

In this study, we aimed to determine drug-resistance mutations (DRMs) in HIV-1 patients with low-level viremia (LLV) and explored the performance of next-generation sequencing (NGS) in detecting HIV DRMs by using LLV samples.

METHODS

Overall, 80 samples with LLV were amplified and sequenced using a commercial Sanger sequencing (SS) genotyping kit. Furthermore, 51 samples successfully sequenced using SS were simultaneously subjected to NGS. Genotyping success rates of various viremia categories by two sequencing methods were calculated. Stanford HIV-1 drug-resistance database (HIVdb version 8.9) was used to analyze the DRMs. In the NGS assay, a threshold of 5% was considered for reporting low-frequency variants, and the DRMs detected using SS and NGS were compared.

RESULTS

The overall success rate of PR/RT regions was 88.1% (67/80) using SS and 86.3% (44/51) using NGS. Furthermore, a significant linear trend was noted between viral load and the genotyping success rate. A total of 38.8% (26/67) participants harbored at least one mutation, as revealed through SS. Moreover, the prevalence of DRMs in persistent LLV was significantly higher than that in intermittent LLV (62.1% vs. 21.1%; P < 0.05). A total of 69 DRMs were detected using the two sequencing methods at the threshold of 5%. Moreover, 10 DRMs missed by SS were detected using NGS, whereas 8 DRMs missed by NGS were detected by SS.

CONCLUSION

Our data suggested that the genotyping resistance testing is necessary to guide antiretroviral therapy optimization in LLV patients.

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.

Low-frequency HIV-1 drug resistance mutations in antiretroviral naïve individuals in Botswana

BACKGROUND

Individuals living with human immunodeficiency virus (HIV) who experience virological failure (VF) after combination antiretroviral therapy (cART) initiation may have had low-frequency drug resistance mutations (DRMs) at cART initiation. There are no data on low-frequency DRMs among cART-naïve HIV-positive individuals in Botswana.

METHODS

We evaluated the prevalence of low-frequency DRMs among cART-naïve individuals previously sequenced using Sanger sequencing. The generated pol amplicons were sequenced by next-generation sequencing.

RESULTS

We observed low-frequency DRMs (detected at <20% in 33/103 (32%) of the successfully sequenced individuals, of whom four also had mutations detected at >20%. K65R was the most common low-frequency DRM detected in 8 individuals. Eighty-two of the 103 individuals had follow-up viral load data while on cART. Twenty-seven of the 82 individuals harbored low-frequency DRMs. Only 12 of 82 individuals experienced VF. The following low-frequency DRMs were observed in four individuals experiencing VF: K65R, K103N, V108I, and Y188C. No statistically significant difference was observed in the prevalence of low-frequency DRMs between individuals experiencing VF (4/12) and those not experiencing VF (23/70) (P = .97). However, individuals with non-nucleoside reverse transcriptase inhibitors-associated low-frequency DRMs were 2.68 times more likely to experience VF (odds ratio, 2.68; 95% confidential interval, 0.4-13.9) compared with those without (P = .22).

CONCLUSION

Next-generation sequencing was able to detect low-frequency DRMs in this cohort in Botswana, but these DRMs did not contribute significantly to VF.

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.

Reduced and highly diverse peripheral HIV-1 reservoir in virally suppressed patients infected with non-B HIV-1 strains in Uganda

BACKGROUND

Our understanding of the peripheral human immunodeficiency virus type 1 (HIV-1) reservoir is strongly biased towards subtype B HIV-1 strains, with only limited information available from patients infected with non-B HIV-1 subtypes, which are the predominant viruses seen in low- and middle-income countries (LMIC) in Africa and Asia.

RESULTS

In this study, blood samples were obtained from well-suppressed ART-experienced HIV-1 patients monitored in Uganda (n = 62) or the U.S. (n = 50), with plasma HIV-1 loads < 50 copies/ml and CD4+ T-cell counts > 300 cells/ml. The peripheral HIV-1 reservoir, i.e., cell-associated HIV-1 RNA and proviral DNA, was characterized using our novel deep sequencing-based EDITS assay. Ugandan patients were slightly younger (median age 43 vs 49 years) and had slightly lower CD4+ counts (508 vs 772 cells/ml) than U.S. individuals. All Ugandan patients were infected with non-B HIV-1 subtypes (31% A1, 64% D, or 5% C), while all U.S. individuals were infected with subtype B viruses. Unexpectedly, we observed a significantly larger peripheral inducible HIV-1 reservoir in U.S. patients compared to Ugandan individuals (48 vs. 11 cell equivalents/million cells, p < 0.0001). This divergence in reservoir size was verified measuring proviral DNA (206 vs. 88 cell equivalents/million cells, p < 0.0001). However, the peripheral HIV-1 reservoir was more diverse in Ugandan than in U.S. individuals (8.6 vs. 4.7 p-distance, p < 0.0001).

CONCLUSIONS

The smaller, but more diverse, peripheral HIV-1 reservoir in Ugandan patients might be associated with viral (e.g., non-B subtype with higher cytopathicity) and/or host (e.g., higher incidence of co-infections or co-morbidities leading to less clonal expansion) factors. This highlights the need to understand reservoir dynamics in diverse populations as part of ongoing efforts to find a functional cure for HIV-1 infection in LMICs.

Absence of Resistance Mutations in the Integrase Coding Region among ART-Experienced Patients in the Republic of the Congo

BACKGROUND

HIV infects around one hundred thousand patients in the Republic of the Congo. Approximately 25% of them receive an antiretroviral treatment; current first-line regimens include two NRTIs and one NNRTI, reverse transcriptase inhibitors. Recently, protease inhibitors (PIs) were also introduced as second-line therapy upon clinical signs of treatment failure. Due to the limited number of molecular characterizations and amount of drug resistance data available in the Republic of the Congo, this study aims to evaluate the prevalence of circulating resistance mutations within the pol region.

METHODS

HIV-positive, ART-experienced patients have been enrolled in four semi-urban localities in the Republic of the Congo. Plasma samples were collected, and viral RNA was extracted. The viral load for each patient was evaluated by RT-qPCR, following the general diagnostic procedures of the University Hospital of Bordeaux. Finally, drug resistance genotyping and phylogenetic analysis were conducted following Sanger sequencing of the pol region.

RESULTS

A high diversity of HIV-1 strains was observed with many recombinant forms. Drug resistance mutations in RT and PR genes were determined and correlated to HAART. Because integrase inhibitors are rarely included in treatments in the Republic of the Congo, the prevalence of integrase drug resistance mutations before treatment was also determined. Interestingly, very few mutations were observed.

CONCLUSIONS

We confirmed a high diversity of HIV-1 in the Republic of the Congo. Most patients presented an accumulation of mutations conferring resistance against NRTIs, NNRTIs and PIs. Nonetheless, the absence of integrase mutations associated with drug resistance suggests that the introduction of integrase inhibitors into therapy will be highly beneficial to patients in the Republic of the Congo.

Deep sequencing of HIV-1 reveals extensive subtype variation and drug resistance after failure of first-line antiretroviral regimens in Nigeria

BACKGROUND

Deep sequencing could improve understanding of HIV treatment failure and viral population dynamics. However, this tool is often inaccessible in low- and middle-income countries.

OBJECTIVES

To determine the genetic patterns of resistance emerging in West African HIV-1 subtypes during first-line virological failure, and the implications for future antiretroviral options.

PATIENTS AND METHODS

Participants were selected from a Nigerian cohort of people living with HIV who had failed first-line ART and subsequently switched to second-line therapy. Whole HIV-1 genome sequences were generated from first-line virological failure samples with Illumina MiSeq. Mutations detected at ≥2% frequency were analysed and compared by subtype.

RESULTS

HIV-1 sequences were obtained from 101 participants (65% female, median age 30 years, median 32.9 months of nevirapine- or efavirenz-based ART). Thymidine analogue mutations (TAMs) were detected in 61%, other core NRTI mutations in 92% and NNRTI mutations in 99%. Minority variants (<20% frequency) comprised 18% of all mutations. K65R was more prevalent in CRF02_AG than G subtypes (33% versus 7%; P = 0.002), and ≥3 TAMs were more common in G than CRF02_AG (52% versus 24%; P = 0.004). Subtype G viruses also contained more RT cleavage site mutations. Cross-resistance to at least one of the newer NNRTIs, doravirine, etravirine or rilpivirine, was predicted in 81% of participants.

CONCLUSIONS

Extensive drug resistance had accumulated in people with West African HIV-1 subtypes, prior to second-line ART. Deep sequencing significantly increased the detection of resistance-associated mutations. Caution should be used if considering newer-generation NNRTI agents in this setting.

CODEHOP-Mediated PCR Improves HIV-1 Genotyping and Detection of Variants by MinION Sequencing

HIV-1 is genetically heterogeneous, having different subtypes and circulating recombinant forms (CRFs). HIV-1 genotyping is used to determine drug resistance profiles and is based on the use of a mixture of consensus and degenerate primers targeting the pol gene. However, the use of this type of primers is associated with either PCR bias or PCR failure. Consensus-degenerate hybrid oligonucleotide primers (CODEHOPs) can detect and identify unknown and distantly related gene sequences by PCR. CODEHOPs designed using different HIV-1 subtypes and CRFs were evaluated for HIV-1 genotyping by Sanger and MinION sequencing. A total of 321 plasma samples were used for the validation of CODEHOP-mediated HIV-1 genotyping. CODEHOP-mediated PCR showed 100% sensitivity and specificity, with limits of detection and genotyping below 200 copies/ml. The head-to-head evaluation of CODEHOP-mediated PCR and standard PCR showed 97 to 98% and 82 to 84% PCR success rates, respectively. There was 100% agreement between the CODEHOP and the reference method in the drug resistance profiles determined by Sanger-based sequencing. Using MinION sequencing, the CODEHOP-mediated PCR scheme resulted in better depth of genome coverage and detection of more drug resistance variants in the protease and reverse transcriptase genes than the standard amplification scheme. The overall prevalences of drug resistance mutations were 17.1% in treatment-experienced patients and 1.2% in treatment-naive patients. They were mainly associated with resistance to reverse transcriptase inhibitors and were linked to virological failure and the patient’s treatment history. Findings from this study suggest that the performance of HIV-1 genotyping is improved by using CODEHOP-mediated PCR.

IMPORTANCE HIV-1 drug resistance is the main cause of treatment failure. Regular surveillance of resistance-associated mutations in HIV-1 genomes is essential for the optimal management of HIV-1 infections. Due to HIV-1’s genetic diversity, different HIV-1 genotypes are circulating worldwide. Standard primers used in the amplification of HIV-1 RNA have not been designed to cover all HIV-1 genotypes and are the main cause of amplification and drug resistance test failure. In this study, new sets of PCR primers targeting the protease, reverse transcriptase, and integrase genes were designed using the CODEHOP approach. They were compared to primers recommended in part by WHO for drug resistance testing using in-house PCR. Unsuccessful HIV-1 RNA amplification was less likely to occur with CODEHOP primers, leading to fewer test failures and lower cost. Furthermore, CODEHOP primers were more effective than standard primers for the detection of minority resistant variants by MinION sequencing.

Drug resistance mutations in protease gene of HIV-1 subtype C infected patient population

Failure of antiretroviral therapy (ART) in HIV-1 infection is a critical issue for the physicians treating HIV patients. The major cause of drug failure is the development of resistance mutations in reverse transcriptase (RT) and/or protease (PR) genes. Mutations associated with drug resistance decrease drug effectiveness. This study was conducted to assess drug resistance profile of the entire PR gene in 90 HIV-1 patients consisting of 23 ART non-responsive, 32 ART responsive and 35 drug naive patients. It was observed that the majority of the sequences (94.4%) belonged to subtype C and (5.5%) to subtype A1. The ART non-responsive and responsive patients were treated with either first line of ART regimen (two NRTI and one NNRTI) or second line of ART regimen that included additional one protease inhibitor (PI). All the patients in each group except one responsive patient had various minor resistance mutations. Thus, drug failures in ART non-responsive patients may not always be due to drug resistance mutations instead other factors may also be responsible for drug failures such as non-compliance, suboptimal dose or drug interaction. The presence of minor drug resistance mutations in drug naive patients is suggestive of transmitted resistance mutations.

Impact of pre-existing drug resistance on risk of virological failure in South Africa

OBJECTIVES

There is conflicting evidence on the impact of pre-existing HIV drug resistance mutations (DRMs) in patients infected with non-B subtype virus.

METHODS

We performed a case-cohort substudy of the AIDS Drug Resistance Surveillance Study, which enrolled South African patients initiating first-line efavirenz/emtricitabine/tenofovir. Pre-ART DRMs were detected by Illumina sequencing of HIV pol and DRMs present at <20% of the viral population were labelled as minority variants (MVs). Weighted Cox proportional hazards models estimated the association between pre-ART DRMs and risk of virological failure (VF), defined as confirmed HIV-1 RNA ≥1000 copies/mL after ≥5 months of ART.

RESULTS

The evaluable population included 178 participants from a randomly selected subcohort (16 with VF, 162 without VF) and 83 additional participants with VF. In the subcohort, 16% of participants harboured ≥1 majority DRM. The presence of any majority DRM was associated with a 3-fold greater risk of VF (P = 0.002), which increased to 9.2-fold (P < 0.001) in those with <2 active drugs. Thirteen percent of participants harboured MV DRMs in the absence of majority DRMs. Presence of MVs alone had no significant impact on the risk of VF. Inclusion of pre-ART MVs with majority DRMs improved the sensitivity but reduced the specificity of predicting VF.

CONCLUSIONS

In a South African cohort, the presence of majority DRMs increased the risk of VF, especially for participants receiving <2 active drugs. The detection of drug-resistant MVs alone did not predict an increased risk of VF, but their inclusion with majority DRMs affected the sensitivity/specificity of predicting VF.

Impact of archived M184V/I mutation on the effectiveness of switch to co-formulated elvitegravir, cobicistat, emtricitabine and tenofovir alafenamide among virally suppressed people living with HIV

OBJECTIVES

Real-world experience regarding the effectiveness of co-formulated elvitegravir, cobicistat, emtricitabine and tenofovir alafenamide (EVG/C/FTC/TAF) as a switch regimen is sparse among people living with HIV (PLWH) harbouring the M184V/I mutation with or without thymidine analogue-associated mutations (TAMs).

METHODS

In this retrospective multicentre study, PLWH who were switched to EVG/C/FTC/TAF after having achieved viral suppression (plasma HIV RNA <200 copies/mL) for 6 months or longer were included. Patients with archived M184V/I mutation (case patients) were matched to controls without M184V/I mutation at a 1:4 ratio. Patients with a history of virological failure or resistance to elvitegravir were excluded. The primary endpoint was virological non-success (plasma HIV RNA ≥50 copies/mL) at Week 48 of switch using a modified FDA snapshot analysis.

RESULTS

Overall, 100 case patients with the M184V/I mutation were identified, including 6 (6.0%) with K65R and 13 (13.0%) with at least one TAM, and were matched to 400 controls in terms of gender, age (mean = 40.3 versus 39.7 years) and cumulative exposure duration to tenofovir disoproxil fumarate (median = 146 versus 143 weeks). At Week 48, the rate of virological non-success for the case patients and controls was 5.0% (5/100) and 3.3% (13/400), respectively (difference = 1.7%; 95% CI = -2.9%-6.3%), while the rate of virological success was 88.0% and 89.5% for the case patients and controls, respectively. The presence of the K65R mutation or TAMs was not associated with virological non-response.

CONCLUSIONS

Among virally suppressed PLWH, EVG/C/FTC/TAF is effective in maintaining viral suppression at Week 48 despite archived M184V/I mutation with or without TAMs.

Investigation of drug resistance against protease, reverse transcriptase, and integrase inhibitors by next-generation sequencing in HIV-positive patients

Our aim was to investigate the mutations in protease (PR), reverse transcriptase (RT), and integrase (IN) gene regions in human immunodeficiency virus (HIV) using a single amplicon via next-generation sequencing (NGS). The study included plasma samples from 49 HIV-1-positive patients, which were referred for HIV-1 drug resistance testing during 2017. A nested polymerase chain reaction (PCR) was performed after the RNA extraction and one-step reverse transcription stages. The sequencing of the HIV genome in the PR, RT, and IN gene regions was carried out using MiSeq NGS technology. Sanger sequencing (SS) was used to analyze resistance mutations in the PR and RT gene regions using a ViroSeq HIV-1 Genotyping System. PCR products were analyzed with an ABI3500 GeneticAnalyzer (Applied Biosystems). Resistance mutations detected with NGS at frequencies above 20% were identical to the SS results. Resistance to at least one antiretroviral (ARV) drug was 22.4% (11 of 49) with NGS and 10.2% (5 of 49) with SS. At least one low-frequency resistance mutation was detected in 18.3% (9 of 49) of the samples. Low-frequency resistance mutations resulted in virological failure in only one patient. The cost of the analyses was reduced by sample pooling and multiplex analysis using the MiSeq system. This is the first study in Turkey to use NGS technologies for the detection of resistance mutations in all three gene (PR, RT, IN) regions using a single amplicon. Our findings suggest that NGS is more sensitive and cost-effective than the SS method.

Detection of pretreatment minority HIV-1 reverse transcriptase inhibitor-resistant variants by ultra-deep sequencing has a limited impact on virological outcomes

OBJECTIVES

Ultra-deep sequencing (UDS) is a powerful tool for exploring the impact on virological outcome of minority variants with low frequencies, some even <1% of the virus population. Here, we compared HIV-1 minority variants at baseline, through plasma RNA and PBMC DNA analyses, and the dominant variants at the virological failure (VF) point, to evaluate the impact of minority drug-resistant variants (MDRVs) on virological outcomes.

METHODS

Single-molecule real-time sequencing (SMRTS) was performed on baseline RNA and DNA. The Stanford HIV-1 drug resistance database was used for the identification and evaluation of drug resistance-associated mutations (DRAMs).

RESULTS

We classified 50 patients into virological success (VS) and VF groups. We found that the rates of reverse transcriptase inhibitor (RTI) DRAMs determined by SMRTS did not differ significantly within or between groups, whether based on RNA or DNA analyses. There was no significant difference in the level of resistance to specific drugs between groups, in either DNA or RNA analyses, except for the DNA-based analysis of lamivudine, for which there was a trend towards a higher prevalence of intermediate/high-level resistance in the VF group. The RNA MDRVs corresponded to DNA MDRVs, except for M100I and Y188H. Sequencing from DNA appeared to be more sensitive than from RNA to detect MDRVs.

CONCLUSIONS

Detection of pretreatment minority HIV-1 RTI-resistant variants by UDS showed that MDRVs at baseline were not significantly associated with virological outcome. However, HIV-1 DNA sequencing by UDS was useful for detecting pretreatment drug resistance mutations in patients, potentially affecting virological responses, suggesting a potential clinical relevance for ultra-deep DNA sequencing.

Applying next-generation sequencing to unravel the mutational landscape in viral quasispecies

Next-generation sequencing (NGS) has revolutionized the scale and depth of biomedical sciences. Because of its unique ability for the detection of sub-clonal variants within genetically diverse populations, NGS has been successfully applied to analyze and quantify the exceptionally-high diversity within viral quasispecies, and many low-frequency drug- or vaccine-resistant mutations of therapeutic importance have been discovered. Although many works have intensively discussed the latest NGS approaches and applications in general, none of them has focused on applying NGS in viral quasispecies studies, mostly due to the limited ability of current NGS technologies to accurately detect and quantify rare viral variants. Here, we summarize several error-correction strategies that have been developed to enhance the detection accuracy of minority variants. We also discuss critical considerations for preparing a sequencing library from viral RNAs and for analyzing NGS data to unravel the mutational landscape.

A single variant sequencing method for sensitive and quantitative detection of HIV-1 minority variants

HIV drug resistance is a major threat to achieving long-term viral suppression in HIV-positive individuals. Drug resistant HIV variants, including minority variants, can compromise response to antiretroviral therapy. Many studies have investigated the clinical relevance of drug resistant minority variants, but the level at which minority variants become clinically relevant remains unclear. A combination of Primer-ID and deep sequencing is a promising approach that may quantify minority variants more accurately compared to standard deep sequencing. However, most studies that used the Primer-ID method have analyzed clinical samples directly. Thus, its sensitivity and quantitative accuracy have not been adequately validated using known controls. Here, we constructed defined proportions of artificial RNA and virus quasispecies and measured their relative proportions using the Primer-ID based, quantitative single-variant sequencing (qSVS) assay. Our results showed that minority variants present at 1% of quasispecies were detected reproducibly with minimal variations between technical replicates. In addition, the measured frequencies were comparable to the expected frequencies. These data validate the accuracy and reproducibility of the qSVS assay in quantifying authentic HIV minority variants, and support the use of this approach to examine the impacts of minority HIV variants on virologic response and clinical outcome.

Adherence to antiretroviral therapy among HIV infected pregnant women in public health sectors: a pilot of Chilenje level one Hospital Lusaka, Zambia

Introduction

regular use of Antiretroviral Therapy (ART) in pregnancy and breastfeeding reduces the odds of Mother-to-Child HIV Transmission (MTCT). However, adherence to ART is critical for MTCT to be successful. The present study investigated factors that influence adherence to ART among HIV infected pregnant women in Zambia.

Methods

a cross-sectional study design was conducted involving 71 HIV infected pregnant women who were advised to join the Prevention of Mother-to-Child HIV Transmission (PMTCT) program during their routine Antenatal clinic (ANC) visit and were on ART for more than six months. We used the Medication Possession Ratio (MPR) to quantify adherence levels. We used logistic regression to establish factors that influence adherence to ART.

Results

a total of 71 HIV infected pregnant women with a median age of 27years (IQR, 25-30) were enrolled in the study. There was evidence of a difference in adherence levels between pregnant women above 30 years and ones between 15 years and 30 years (P<0.001). Median adherence levels in this group were found to be at 96%(IQR 89-97). The main predictor of adherence in this population was marital status (being on separation) and age. The women who were on separation were 0.14 times less likely to adhere to option B+ compared to married women.

Conclusion

adherence to option B+ among pregnant women is low. Adherence was significantly influenced by marital status (being on separation) and age. Efforts to improve adherence should be directed towards women on separation and young adults (< 30 years of age).

A cross-sectional study to characterize local HIV-1 dynamics in Washington, DC using next-generation sequencing

Washington, DC continues to experience a generalized HIV-1 epidemic. We characterized the local phylodynamics of HIV-1 in DC using next-generation sequencing (NGS) data. Viral samples from 68 participants from 2016 through 2017 were sequenced and paired with epidemiological data. Phylogenetic and network inferences, drug resistant mutations (DRMs), subtypes and HIV-1 diversity estimations were completed. Haplotypes were reconstructed to infer transmission clusters. Phylodynamic inferences based on the HIV-1 polymerase (pol) and envelope genes (env) were compared. Higher HIV-1 diversity (n.s.) was seen in men who have sex with men, heterosexual, and male participants in DC. 54.0% of the participants contained at least one DRM. The 40-49 year-olds showed the highest prevalence of DRMs (22.9%). Phylogenetic analysis of pol and env sequences grouped 31.9-33.8% of the participants into clusters. HIV-TRACE grouped 2.9-12.8% of participants when using consensus sequences and 9.0-64.2% when using haplotypes. NGS allowed us to characterize the local phylodynamics of HIV-1 in DC more broadly and accurately, given a better representation of its diversity and dynamics. Reconstructed haplotypes provided novel and deeper phylodynamic insights, which led to networks linking a higher number of participants. Our understanding of the HIV-1 epidemic was expanded with the powerful coupling of HIV-1 NGS data with epidemiological data.

Prospective Evaluation of the Vela Diagnostics Next-Generation Sequencing Platform for HIV-1 Genotypic Resistance Testing

Genotypic antiretroviral drug resistance testing is a critical component of the global efforts to control the HIV-1 epidemic. This study investigates the semiautomated, next-generation sequencing (NGS)-based Vela Diagnostics Sentosa SQ HIV-1 Genotyping Assay in a prospective cohort of HIV-1-infected patients. Two-hundred sixty-nine samples were successfully sequenced by both NGS and Sanger sequencing. Among the 261 protease/reverse transcriptase (PR/RT) sequences, a mean of 0.37 drug resistance mutations were identified by both Sanger and NGS, 0.08 by NGS alone, and 0.03 by Sanger alone. Among the 50 integrase sequences, a mean of 0.3 drug resistance mutations were detected by both Sanger and NGS, and 0.08 by NGS alone. NGS estimated higher levels of drug resistance to one or more antiretroviral drugs for 6.5% of PR/RT sequences and 4.0% of integrase sequences, whereas Sanger estimated higher levels of drug resistance for 3.8% of PR/RT sequences. Although the samples successfully sequenced by the Sentosa SQ HIV Genotyping Assay demonstrated similar predicted resistance compared with Sanger, 44% of Sentosa runs failed quality control requiring 17 additional runs. This semi-automated NGS-based assay may aid in HIV-1 genotypic drug resistance testing, though numerous quality control issues were observed when this platform was used in a clinical laboratory setting. With additional refinement, the Sentosa SQ HIV-1 Genotyping Assay may contribute to the global efforts to control HIV-1.

Drug resistance and optimizing dolutegravir regimens for adolescents and young adults failing antiretroviral therapy

OBJECTIVES

The integrase strand inhibitor dolutegravir (DTG) combined with tenofovir and lamivudine (TLD) is a single tablet regimen recommended for 1st, 2nd and 3rd-line public health antiretroviral therapy (ART). We determined drug resistance mutations (DRMs) and evaluated the predictive efficacy of a TLD containing regimen for viremic adolescents and young adults in Harare, Zimbabwe.

METHODS

We sequenced plasma viral RNA from HIV-1-infected adolescents and young adults on 1st and 2nd-line ART with confirmed virologic failure (viral load >1000 copies/ml) and calculated total genotypic susceptibility scores to current 2nd, 3rd line and DTG regimens.

RESULTS

A total of 160 participants were genotyped; 112 (70%) on 1st line and 48 (30%) on 2nd line, median (interquartile range) age 18 (15-19) and duration of ART (interquartile range) was 6 (4-8) years. Major DRMs were present in 94 and 67% of 1st and 2nd-line failures, respectively (P < 0.001). Dual class resistance to nucleotide reverse transcriptase inhibitors and nonnucleotide reverse transcriptase inhibitors was detected in 96 (60%) of 1st-line failures; protease inhibitor DRMs were detected in a minority (10%) of 2nd-line failures. A total genotypic susceptibility score of 2 or less may risk protease inhibitor or DTG monotherapy in 11 and 42% of 1st-line failures switching to 2nd-line protease inhibitor and TLD respectively.

CONCLUSION

Among adolescents and young adults, current protease inhibitor-based 2nd-line therapies are poorly tolerated, more expensive and adherence is poor. In 1st-line failure, implementation of TLD for many adolescents and young adults on long-term ART may require additional active drug(s). Drug resistance surveillance and susceptibility scores may inform strategies for the implementation of TLD.

The urgent need for more potent antiretroviral therapy in low-income countries to achieve UNAIDS 90-90-90 and complete eradication of AIDS by 2030

Background

Over 90% of Human Immunodeficiency Virus (HIV) infected individuals will be on treatment by 2020 under UNAIDS 90–90-90 global targets. Under World Health Organisation (WHO) “Treat All” approach, this number will be approximately 36.4 million people with over 98% in low-income countries (LICs).

Main body

Pretreatment drug resistance (PDR) largely driven by frequently use of non-nucleoside reverse transcriptase inhibitors (NNRTIs), efavirenz and nevirapine, has been increasing with roll-out of combined antiretroviral therapy (cART) with 29% annual increase in some LICs countries. PDR has exceeded 10% in most LICs which warrants change of first line regimen to more robust classes under WHO recommendations. If no change in regimens is enforced in LICs, it’s estimated that over 16% of total deaths, 9% of new infections, and 8% of total cART costs will be contributed by HIV drug resistance by 2030. Less than optimal adherence, and adverse side effects associated with currently available drug regimens, all pose a great threat to achievement of 90% viral suppression and elimination of AIDS as a public health threat by 2030. This calls for urgent introduction of policies that advocate for voluntary and compulsory drug licensing of new more potent drugs which should also emphasize universal access of these drugs to all individuals worldwide.

Conclusions

The achievement of United Nations Programme on HIV and AIDS 2020 and 2030 targets in LICs depends on access to active cART with higher genetic barrier to drug resistance, better safety, and tolerability profiles. It’s also imperative to strengthen quality service delivery in terms of retention of patients to treatment, support for adherence to cART, patient follow up and adequate drug stocks to help achieve a free AIDS generation.

Electronic supplementary material

The online version of this article (10.1186/s40249-019-0573-1) contains supplementary material, which is available to authorized users.

Validation of publicly-available software used in analyzing NGS data for HIV-1 drug resistance mutations and transmission networks in a Washington, DC, Cohort

The DC Cohort is an ongoing longitudinal observational study of persons living with HIV. To better understand HIV-1 drug resistance and potential transmission clusters among these participants, we performed targeted, paired-end next-generation sequencing (NGS) of protease, reverse transcriptase and integrase amplicons. We elected to use free, publicly-available software (HyDRA Web, Stanford HIVdb and HIV-TRACE) for data analyses so that laboratory personnel without extensive bioinformatics expertise could use it; making the approach accessible and affordable for labs worldwide. With more laboratories transitioning away from Sanger-based chemistries to NGS platforms, lower frequency drug resistance mutations (DRMs) can be detected, yet their clinical relevance is uncertain. We looked at the impact choice in cutoff percentage had on number of DRMs detected and found an inverse correlation between the two. Longitudinal studies will be needed to determine whether low frequency DRMs are an early indicator of emerging resistance. We successfully validated this pipeline against a commercial pipeline, and another free, publicly-available pipeline. RT DRM results from HyDRA Web were compared to both SmartGene and PASeq Web; using the Mantel test, R2 values were 0.9332 (p<0.0001) and 0.9097 (p<0.0001), respectively. PR and IN DRM results from HyDRA Web were then compared with PASeq Web only; using the Mantel test, R2 values were 0.9993 (p<0.0001) and 0.9765 (p<0.0001), respectively. Drug resistance was highest for the NRTI drug class and lowest for the PI drug class in this cohort. RT DRM interpretation reports from this pipeline were also highly correlative compared to SmartGene pipeline; using the Spearman's Correlation, rs value was 0.97757 (p<0.0001). HIV-TRACE was used to identify potential transmission clusters to better understand potential linkages among an urban cohort of persons living with HIV; more individuals were male, of black race, with an HIV risk factor of either MSM or High-risk Heterosexual. Common DRMs existed among individuals within a cluster. In summary, we validated a comprehensive, easy-to-use and affordable NGS approach for tracking HIV-1 drug resistance and identifying potential transmission clusters within the community.

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.

Characterization of minority HIV-1 drug resistant variants in the United Kingdom following the verification of a deep sequencing-based HIV-1 genotyping and tropism assay

BACKGROUND

The widespread global access to antiretroviral drugs has led to considerable reductions in morbidity and mortality but, unfortunately, the risk of virologic failure increases with the emergence, and potential transmission, of drug resistant viruses. Detecting and quantifying HIV-1 drug resistance has therefore become the standard of care when designing new antiretroviral regimens. The sensitivity of Sanger sequencing-based HIV-1 genotypic assays is limited by its inability to identify minority members of the quasispecies, i.e., it only detects variants present above ~ 20% of the viral population, thus, failing to detect minority variants below this threshold. It is clear that deep sequencing-based HIV-1 genotyping assays are an important step change towards accurately monitoring HIV-infected individuals.

METHODS

We implemented and verified a clinically validated HIV-1 genotyping assay based on deep sequencing (DEEPGEN™) in two clinical laboratories in the United Kingdom: St. George's University Hospitals Healthcare NHS Foundation Trust (London) and at NHS Lothian (Edinburgh), to characterize minority HIV-1 variants in 109 plasma samples from ART-naïve or -experienced individuals.

RESULTS

Although subtype B HIV-1 strains were highly prevalent (44%, 48/109), most individuals were infected with non-B subtype viruses (i.e., A1, A2, C, D, F1, G, CRF02_AG, and CRF01_AE). DEEPGEN™ was able to accurately detect drug resistance-associated mutations not identified using standard Sanger sequencing-based tests, which correlated significantly with patient's antiretroviral treatment histories. A higher proportion of minority PI-, NRTI-, and NNRTI-resistance mutations was detected in NHS Lothian patients compared to individuals from St. George's, mainly M46I/L and I50 V (associated with PIs), D67 N, K65R, L74I, M184 V/I, and K219Q (NRTIs), and L100I (NNRTIs). Interestingly, we observed an inverse correlation between intra-patient HIV-1 diversity and CD4⁺ T cell counts in the NHS Lothian patients.

CONCLUSIONS

This is the first study evaluating the transition, training, and implementation of DEEPGEN™ between three clinical laboratories in two different countries. More importantly, we were able to characterize the HIV-1 drug resistance profile (including minority variants), coreceptor tropism, subtyping, and intra-patient viral diversity in patients from the United Kingdom, providing a rigorous foundation for basing clinical decisions on highly sensitive and cost-effective deep sequencing-based HIV-1 genotyping assays in the country.

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

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

Comparison of an In Vitro Diagnostic Next-Generation Sequencing Assay with Sanger Sequencing for HIV-1 Genotypic Resistance Testing

The ability of next-generation sequencing (NGS) technologies to detect low frequency HIV-1 drug resistance mutations (DRMs) not detected by dideoxynucleotide Sanger sequencing has potential advantages for improved patient outcomes. We compared the performance of an in vitro diagnostic (IVD) NGS assay, the Sentosa SQ HIV genotyping assay for HIV-1 genotypic resistance testing, with Sanger sequencing on 138 protease/reverse transcriptase (RT) and 39 integrase sequences. The NGS assay used a 5% threshold for reporting low-frequency variants. The level of complete plus partial nucleotide sequence concordance between Sanger sequencing and NGS was 99.9%. Among the 138 protease/RT sequences, a mean of 6.4 DRMs was identified by both Sanger and NGS, a mean of 0.5 DRM was detected by NGS alone, and a mean of 0.1 DRM was detected by Sanger sequencing alone. Among the 39 integrase sequences, a mean of 1.6 DRMs was detected by both Sanger sequencing and NGS and a mean of 0.15 DRM was detected by NGS alone. Compared with Sanger sequencing, NGS estimated higher levels of resistance to one or more antiretroviral drugs for 18.2% of protease/RT sequences and 5.1% of integrase sequences. There was little evidence for technical artifacts in the NGS sequences, but the G-to-A hypermutation was detected in three samples. In conclusion, the IVD NGS assay evaluated in this study was highly concordant with Sanger sequencing. At the 5% threshold for reporting minority variants, NGS appeared to attain a modestly increased sensitivity for detecting low-frequency DRMs without compromising sequence accuracy.

Molecular evolution methods to study HIV-1 epidemics

Nucleotide sequences of HIV isolates are obtained routinely to evaluate the presence of resistance mutations to antiretroviral drugs. But, beyond their clinical use, these and other viral sequences include a wealth of information that can be used to better understand and characterize the epidemiology of HIV in relevant populations. In this review, we provide a brief overview of the main methods used to analyze HIV sequences, the data bases where reference sequences can be obtained, and some caveats about the possible applications for public health of these analyses, along with some considerations about their limitations and correct usage to derive robust and reliable conclusions.

Absence of HIV-1 Drug Resistance Mutations Supports the Use of Dolutegravir in Uganda

To screen for drug resistance and possible treatment with Dolutegravir (DTG) in treatment-naive patients and those experiencing virologic failure during first-, second-, and third-line combined antiretroviral therapy (cART) in Uganda. Samples from 417 patients in Uganda were analyzed for predicted drug resistance upon failing a first- (N = 158), second- (N = 121), or third-line [all 51 involving Raltegravir (RAL)] treatment regimen. HIV-1 pol gene was amplified and sequenced from plasma samples. Drug susceptibility was interpreted using the Stanford HIV database algorithm and SCUEAL was used for HIV-1 subtyping. Frequency of resistance to nucleoside reverse transcriptase inhibitors (NRTIs) (95%) and non-NRTI (NNRTI, 96%) was high in first-line treatment failures. Despite lack of NNRTI-based treatment for years, NNRTI resistance remained stable in 55% of patients failing second-line or third-line treatment, and was also at 10% in treatment-naive Ugandans. DTG resistance (n = 366) was not observed in treatment-naive individuals or individuals failing first- and second-line cART, and only found in two patients failing third-line cART, while 47% of the latter had RAL- and Elvitegravir-resistant HIV-1. Secondary mutations associated with DTG resistance were found in 2%-10% of patients failing third-line cART. Of 14 drugs currently available for cART in Uganda, resistance was readily observed to all antiretroviral drugs (except for DTG) in Ugandan patients failing first-, second-, or even third-line treatment regimens. The high NNRTI resistance in first-line treatment in Uganda even among treatment-naive patients calls for the use of DTG to reach the UNAIDS 90:90:90 goals.

Characterization of HIV-1 Near Full-Length Proviral Genome Quasispecies from Patients with Undetectable Viral Load Undergoing First-Line HAART Therapy

Increased access to highly active antiretroviral therapy (HAART) by human immunodeficiency virus postive (HIV⁺) individuals has become a reality worldwide. In Brazil, HAART currently reaches over half of HIV-infected subjects. In the context of a remarkable HIV-1 genetic variability, highly related variants, called quasispecies, are generated. HIV quasispecies generated during infection can influence virus persistence and pathogenicity, representing a challenge to treatment. However, the clinical relevance of minority quasispecies is still uncertain. In this study, we have determined the archived proviral sequences, viral subtype and drug resistance mutations from a cohort of HIV⁺ patients with undetectable viral load undergoing HAART as first-line therapy using next-generation sequencing for near full-length virus genome (NFLG) assembly. HIV-1 consensus sequences representing NFLG were obtained for eleven patients, while for another twelve varying genome coverage rates were obtained. Phylogenetic analysis showed the predominance of subtype B (83%; 19/23). Considering the minority variants, 18 patients carried archived virus harboring at least one mutation conferring antiretroviral resistance; for six patients, the mutations correlated with the current ARVs used. These data highlight the importance of monitoring HIV minority drug resistant variants and their clinical impact, to guide future regimen switches and improve HIV treatment success.

The Role of HIV-1 Drug-Resistant Minority Variants in Treatment Failure

Human immunodeficiency virus type 1 (HIV-1) drug resistance genotyping is recommended to help in the selection of antiretroviral therapy and to prevent virologic failure. There are several ultrasensitive assays able to detect HIV-1 drug-resistance minority variants (DRMVs) not detectable by standard population sequencing-based HIV genotyping assays. Presence of these DRMVs has been shown to be clinically relevant, but its impact does not appear to be uniform across drug classes. In this review, we summarize key evidence for the clinical impact of DRMVs across drug classes for both antiretroviral treatment-naive and antiretroviral treatment-experienced patients, and highlight areas where more supporting evidence is needed.

Sensitive detection of HIV-1 resistance to Zidovudine and impact on treatment outcomes in low- to middle-income countries

Background

Thymidine analogs, namely AZT (Zidovudine or Retrovir™) and d4T (Stavudine or Zerit™) are antiretroviral drugs still employed in over 75% of first line combination antiretroviral therapy (cART) in Kampala, Uganda despite aversion to prescribing these drugs for cART in high income countries due in part to adverse events. For this study, we explored how the continued use of these thymidine analogs in cART could impact emergence of drug resistance and impact on future treatment success in Uganda, a low-income country.

Methods

We examined the drug resistance genotypes by Sanger sequencing of 262 HIV-infected patients failing a first line combined antiretroviral treatment containing either AZT or d4T, which represents approximately 5% of the patients at the Joint Clinical Research Center receiving a AZT or d4T containing treatment. Next generation sequencing (DEEPGEN™HIV) and multiplex oligonucleotide ligation assays (AfriPOLA) were then performed on a subset of patient samples to detect low frequency drug resistant mutations. CD4 cell counts, viral RNA loads, and treatment changes were analyzed in a cohort of treatment success and failures.

Results

Over 80% of patients failing first line AZT/d4T-containing cART had predicted drug resistance to 3TC (Lamivudine) and non-nucleoside RT inhibitors (NNRTIs) in the treatment regimen but only 45% had resistance AZT/d4T associated resistance mutations (TAMs). TAMs were however detected at low frequency within the patients HIV quasispecies (1–20%) in 21 of 34 individuals who were failing first-line AZT-containing cART and lacked TAMs by Sanger. Due to lack of TAMs by Sanger, AZT was typically maintained in second-line therapies and these patients had a low frequency of subsequent virologic success.

Conclusions

Our findings suggest that continued use of AZT and d4T in first-line treatment in low-to-middle income countries may lead to misdiagnosis of HIV-1 drug resistance and possibly enhance a succession of second- and third-line treatment failures.

Electronic supplementary material

The online version of this article (doi: 10.1186/s40249-017-0377-0) contains supplementary material, which is available to authorized users.

High level of HIV-1 drug resistance mutations in patients with unsuppressed viral loads in rural northern South Africa

Background

Combination antiretroviral therapy (cART) has significantly reduced HIV morbidity and mortality in both developed and developing countries. However, the sustainability of cART may be compromised by the emergence of viral drug resistance mutations (DRM) and the cellular persistence of proviruses carrying these DRM. This is potentially a more serious problem in resource limited settings.

Methods

DRM were evaluated in individuals with unsuppressed viral loads after first or multiple lines of cART at two sites in rural Limpopo, South Africa. Seventy-two patients with viral loads of >1000 copies/ml were recruited between March 2014 and December 2015. Complete protease (PR) and partial Reverse Transcriptase (RT) sequences were amplified from both plasma RNA and paired proviral DNA from 35 of these subjects. Amplicons were directly sequenced to determine subtype and DRM using the Stanford HIV Drug Resistance Interpretation algorithm.

Results

Among the 72 samples, 69 could be PCR amplified from RNA and 35 from both RNA and DNA. Sixty-five (94.2%) viruses were subtype C, while one was subtype B (1.4%), one recombinant K/C, one recombinant C/B and one unclassified. Fifty-eight (84%) sequences carried at least one DRM, while 11 (15.9%) displayed no DRM. DRM prevalence according to drug class was: NRTI 60.8% NNRTI 65.2%, and PI 5.8%. The most common DRMs were; M184V (51.7%), K103N (50%), V106M (20.6%), D67N (13.3%), K65R (12%). The frequency of the DRM tracked well with the frequency of use of medications to which the mutations were predicted to confer resistance. Interestingly, a significant number of subjects showed predicted resistance to the newer NNRTIs, etravirine (33%) and rilpivirine (42%), both of which are not yet available in this setting. The proportion of DRM in RNA and DNA were mostly similar with the exception of the thymidine analogue mutations (TAMs) D67N, K70R, K219QE; and K103N which were slightly more prevalent in DNA than RNA. Subjects who had received cART for at least 5 years were more likely to harbour >2 DRM (p < 0.05) compared to those treated for a shorter period. DRM were more prevalent in this rural setting compared to a neighbouring urban setting.

Conclusion

We found a very high prevalence of NRTI and NNRTI DRM in patients from rural Limpopo settings with different durations of treatment. The prevalence was significantly higher than those reported in urban settings in South Africa. The dominance of NNRTI based mutations late in treatment supports the use of PI based regimens for second line treatment in this setting. The slight dominance of TAMs in DNA from infected PBMCs compared to plasma virus requires further studies that should include cART subjects with suppressed virus. Such studies will improve our understanding of the pattern of drug resistance and dynamics of viral persistence in these rural settings.

Decreased emergence of HIV-1 drug resistance mutations in a cohort of Ugandan women initiating option B+ for PMTCT

BACKGROUND

Since 2012, WHO guidelines for the prevention of mother-to-child transmission (PMTCT) of HIV-1 in resource-limited settings recommend the initiation of lifelong antiretroviral combination therapy (cART) for all pregnant HIV-1 positive women independent of CD4 count and WHO clinical stage (Option B+). However, long-term outcomes regarding development of drug resistance are lacking until now. Therefore, we analysed the emergence of drug resistance mutations (DRMs) in women initiating Option B+ in Fort Portal, Uganda, at 12 and 18 months postpartum (ppm).

METHODS AND FINDINGS

124 HIV-1 positive pregnant women were enrolled within antenatal care services in Fort Portal, Uganda. Blood samples were collected at the first visit prior starting Option B+ and postpartum at week six, month six, 12 and 18. Viral load was determined by real-time RT-PCR. An RT-PCR covering resistance associated positions in the protease and reverse transcriptase HIV-1 genomic region was performed. PCR-positive samples at 12/18 ppm and respective baseline samples were analysed by next generation sequencing regarding HIV-1 drug resistant variants including low-frequency variants. Furthermore, vertical transmission of HIV-1 was analysed. 49/124 (39.5%) women were included into the DRM analysis. Virological failure, defined as >1000 copies HIV-1 RNA/ml, was observed in three and seven women at 12 and 18 ppm, respectively. Sequences were obtained for three and six of these. In total, DRMs were detected in 3/49 (6.1%) women. Two women displayed dual-class resistance against all recommended first-line regimen drugs. Of 49 mother-infant-pairs no infant was HIV-1 positive at 12 or 18 ppm.

CONCLUSION

Our findings suggest that the WHO-recommended Option B+ for PMTCT is effective in a cohort of Ugandan HIV-1 positive pregnant women with regard to the low selection rate of DRMs and vertical transmission. Therefore, these results are encouraging for other countries considering the implementation of lifelong cART for all pregnant HIV-1 positive women.

Clinical and biological insights from viral genome sequencing

Whole-genome sequencing (WGS) of pathogens is becoming increasingly important not only for basic research but also for clinical science and practice. In virology, WGS is important for the development of novel treatments and vaccines, and for increasing the power of molecular epidemiology and evolutionary genomics. In this Opinion article, we suggest that WGS of viruses in a clinical setting will become increasingly important for patient care. We give an overview of different WGS methods that are used in virology and summarize their advantages and disadvantages. Although there are only partially addressed technical, financial and ethical issues in regard to the clinical application of viral WGS, this technique provides important insights into virus transmission, evolution and pathogenesis.