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

From viral democratic genomes to viral wild bunch of quasispecies

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

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

BACKGROUND

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Epidemiological data analysis of viral quasispecies in the next-generation sequencing era

The unprecedented coverage offered by next-generation sequencing (NGS) technology has facilitated the assessment of the population complexity of intra-host RNA viral populations at an unprecedented level of detail. Consequently, analysis of NGS datasets could be used to extract and infer crucial epidemiological and biomedical information on the levels of both infected individuals and susceptible populations, thus enabling the development of more effective prevention strategies and antiviral therapeutics. Such information includes drug resistance, infection stage, transmission clusters and structures of transmission networks. However, NGS data require sophisticated analysis dealing with millions of error-prone short reads per patient. Prior to the NGS era, epidemiological and phylogenetic analyses were geared toward Sanger sequencing technology; now, they must be redesigned to handle the large-scale NGS datasets and properly model the evolution of heterogeneous rapidly mutating viral populations. Additionally, dedicated epidemiological surveillance systems require big data analytics to handle millions of reads obtained from thousands of patients for rapid outbreak investigation and management. We survey bioinformatics tools analyzing NGS data for (i) characterization of intra-host viral population complexity including single nucleotide variant and haplotype calling; (ii) downstream epidemiological analysis and inference of drug-resistant mutations, age of infection and linkage between patients; and (iii) data collection and analytics in surveillance systems for fast response and control of outbreaks.

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.

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.

Ultradeep sequencing of B and non-B HIV-1 subtypes: Viral diversity and drug resistance mutations before and after one month of antiretroviral therapy in naive patients

BACKGROUND

Ultradeep pyrosequencing technologies permit an assessment of the genetic diversity and the presence and frequency of minority variants in a viral population. The effect of these parameters on the outcome of highly active antiretroviral therapy (HAART) in HIV-infected patients is poorly understood.

OBJECTIVES

The present study used the pyrosequencing Roche 454 prototype assay to determine whether antiretroviral efficacy is correlated with viral diversity and minority drug resistance mutations in HIV-infected treatment-naive patients and to compare assay performance in B and non-B subtypes.

STUDY DESIGN

The study included 30 HIV-1 infected naive patients (20 with subtype non-B and 10 with subtype B). Ultradeep pyrosequencing of protease and reverse transcriptase genes was performed at baseline and 1 month after HAART initiation. Plasma HIV VL was measured at 0 and after 1, 3, and 6 months of HAART.

RESULTS

Pre-HAART minority drug resistance mutations were observed to NRTI in 4 patients, to NNRTI in 6 patients, and to PI in 1 patient; there was no difference in HAART-induced VL decay between patients. Pre-HAART diversity was significantly correlated with the time elapsed since HIV-1 infection diagnosis, but not with the subtype, VL, or CD4 count. Patients with an undetectable VL after 3 months of HAART had a higher pre-HAART diversity. Pre- and post-HAART diversities were not statistically different. There was no difference in assay performance between subtype B and non-B.

CONCLUSIONS

A high pre-HAART viral diversity might have a positive effect on the outcome of HAART. Pre-therapeutic minority drug resistance mutations are uncommon in naive patients.

Development of a Qualitative Quantitative Polymerase Chain Reaction Test to Identify Patients Failing First-Line Therapy to Non-Nucleotide Reverse Transcriptase Inhibitor

Antiretroviral therapy (ART) can be compromised by selection of drug resistance strains, which can be promoted by lack of adherence during therapy and drug tolerance, and some of these drug-resistant strains can persist for years as minority populations. The K103N drug resistance mutation is selected by the use of non-nucleotide reverse transcriptase inhibitors, including nevirapine or efavirenz (EFV), used in low-income countries. Here we describe the use of a less expensive qualitative point mutation polymerase chain reaction (PMqPCR_K103N) targeting K103N mutation. To validate the use of this methodology, we tested previously sequenced samples from patients treated with highly active ART with viral loads above 2,000 copies/ml and compared the results of our assay with Illumina deep sequencing. Due to its low cost and high specificity, this test is particularly suitable for low-income countries to screen for pretreatment resistance in patients either initiating ART or failing first-line regimens containing EFV.

Fecal bacterial microbiome diversity in chronic HIV-infected patients in China

The purpose of this study was to identify fecal bacterial microbiome changes in patients with chronic human immunodeficiency virus (HIV) infection in China. Bacterial 16S rRNA genes were amplified, sequenced (454 pyrosequencing), and clustered into operational taxonomic units using the QIIME software. Relative abundance at the phylum and genus levels were calculated. Alpha diversity was determined by Chao 1 and observed-species indices, and beta diversity was determined by double principal component analysis using the estimated phylogeny-based unweighted Unifrac distance matrices. Fecal samples of the patients with chronic HIV-infection tended to be enriched with bacteria of the phyla Firmicutes (47.20%±0.43 relative abundance) and Proteobacteria (37.21%±0.36) compared with those of the non-HIV infected controls (17.95%±0.06 and 3.81%±0.02, respectively). Members of the genus Bilophila were exclusively detected in samples of the non-HIV infected controls. Bacteroides and arabacteroides were more abundant in the chronic HIV-infected patients. Our study indicated that chronic HIV-infected patients in China have a fecal bacterial microbiome composition that is largely different from that found in non-HIV infected controls, and further study is needed to evaluate whether microbiome changes play a role in disease complications in the distal gut, including opportunistic infections.

Quantifying Next Generation Sequencing Sample Pre-Processing Bias in HIV-1 Complete Genome Sequencing

Genetic analyses play a central role in infectious disease research. Massively parallelized “mechanical cloning” and sequencing technologies were quickly adopted by HIV researchers in order to broaden the understanding of the clinical importance of minor drug-resistant variants. These efforts have, however, remained largely limited to small genomic regions. The growing need to monitor multiple genome regions for drug resistance testing, as well as the obvious benefit for studying evolutionary and epidemic processes makes complete genome sequencing an important goal in viral research. In addition, a major drawback for NGS applications to RNA viruses is the need for large quantities of input DNA. Here, we use a generic overlapping amplicon-based near full-genome amplification protocol to compare low-input enzymatic fragmentation (Nextera™) with conventional mechanical shearing for Roche 454 sequencing. We find that the fragmentation method has only a modest impact on the characterization of the population composition and that for reliable results, the variation introduced at all steps of the procedure—from nucleic acid extraction to sequencing—should be taken into account, a finding that is also relevant for NGS technologies that are now more commonly used. Furthermore, by applying our protocol to deep sequence a number of pre-therapy plasma and PBMC samples, we illustrate the potential benefits of a near complete genome sequencing approach in routine genotyping.

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

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

Multiplexed highly-accurate DNA sequencing of closely-related HIV-1 variants using continuous long reads from single molecule, real-time sequencing

Single Molecule, Real-Time (SMRT) Sequencing (Pacific Biosciences, Menlo Park, CA, USA) provides the longest continuous DNA sequencing reads currently available. However, the relatively high error rate in the raw read data requires novel analysis methods to deconvolute sequences derived from complex samples. Here, we present a workflow of novel computer algorithms able to reconstruct viral variant genomes present in mixtures with an accuracy of >QV50. This approach relies exclusively on Continuous Long Reads (CLR), which are the raw reads generated during SMRT Sequencing. We successfully implement this workflow for simultaneous sequencing of mixtures containing up to forty different >9 kb HIV-1 full genomes. This was achieved using a single SMRT Cell for each mixture and desktop computing power. This novel approach opens the possibility of solving complex sequencing tasks that currently lack a solution.

Risk of drug resistance among persons acquiring HIV within a randomized clinical trial of single- or dual-agent preexposure prophylaxis

BACKGROUND

Preexposure prophylaxis (PrEP) with emtricitabine plus tenofovir disoproxil fumarate (FTC/TDF) or TDF alone reduces the risk of human immunodeficiency virus (HIV) acquisition. Understanding the risk of antiretroviral resistance selected by PrEP during breakthrough infections is important because of the risk of treatment failure during subsequent antiretroviral use.

METHODS

Within the largest randomized trial of FTC/TDF versus TDF as PrEP, plasma samples were tested for HIV with resistance mutations associated with FTC (K65R and M184IV) and TDF (K65R and K70E), using 454 sequencing.

RESULTS

Of 121 HIV seroconverters, 25 received FTC/TDF, 38 received TDF, and 58 received placebo. Plasma drug levels in 26 individuals indicated PrEP use during or after HIV acquisition, of which 5 had virus with resistance mutations associated with their PrEP regimen. Among those with PrEP drug detected during infection, resistance was more frequent in the FTC/TDF arm (4 of 7 [57%]), compared with the TDF arm (1 of 19 [5.3%]; P = .01), owing to the FTC-associated mutation M184IV. Of these cases, 3 had unrecognized acute infection at PrEP randomization, and 2 were HIV negative at enrollment.

CONCLUSIONS

These results suggest that resistance selected by PrEP is rare but can occur both with PrEP initiation during acute seronegative HIV infection and in PrEP breakthrough infections and that FTC is associated with a greater frequency of resistance mutations than TDF.

A near-full length genotypic assay for HCV1b

A near-full genome genotypic assay for HCV1b was developed, which may prove useful to investigate antiviral drug resistance, given new combination therapies for HCV1 infection. The assay consists of three partially overlapping PCRs followed by Sanger population or Illumina next-generation sequencing. Seventy-seven therapy-naïve samples, spanning the entire diversity range of currently known HCV1b, were used for optimization of PCRs, of which ten were sequenced using Sanger and of these ten, four using Illumina. The median detection limits for the three regions, 5'UTR-NS2, E2-NS5A and NS4B-NS5B, were 570, 5670 and 56,670 IU/ml respectively. The number of Illumina reads mapped varied according to the software used, Segminator II being the best performing (81%). Consensus Illumina and Sanger sequencing results accord largely (0.013% major discordances). Differences were due almost exclusively to a larger number of ambiguities (presumably minority variants) scored by Illumina (1.50% minor discordances). The assay is easy to perform in an equipped laboratory; nevertheless, it was difficult to reach high sensitivity and reproducibility, due to the high genetic viral variability. This assay proved to be suitable for detecting drug resistance mutations and can also be used for epidemiological research, even though only a limited set of samples was used for validation.

Contribution of human immunodeficiency virus type 1 minority variants to reduced drug susceptibility in patients on an integrase strand transfer inhibitor-based therapy

The role of HIV-1 minority variants on transmission, pathogenesis, and virologic failure to antiretroviral regimens has been explored; however, most studies of low-level HIV-1 drug-resistant variants have focused in single target regions. Here we used a novel HIV-1 genotypic assay based on deep sequencing, DEEPGEN (Gibson et al 2014 Antimicrob Agents Chemother 58∶2167) to simultaneously analyze the presence of minority variants carrying mutations associated with reduced susceptibility to protease (PR), reverse transcriptase (RT), and integrase strand transfer integrase inhibitors (INSTIs), as well as HIV-1 coreceptor tropism. gag-p2/NCp7/p1/p6/pol-PR/RT/INT and env/C2V3 PCR products were obtained from twelve heavily treatment-experienced patients experiencing virologic failure while participating in a 48-week dose-ranging study of elvitegravir (GS-US-183-0105). Deep sequencing results were compared with (i) virological response to treatment, (ii) genotyping based on population sequencing, (iii) phenotyping data using PhenoSense and VIRALARTS, and (iv) HIV-1 coreceptor tropism based on the phenotypic test VERITROP. Most patients failed the antiretroviral regimen with numerous pre-existing mutations in the PR and RT, and additionally newly acquired INSTI-resistance mutations as determined by population sequencing (mean 9.4, 5.3, and 1.4 PI- RTI-, and INSTI-resistance mutations, respectively). Interestingly, since DEEPGEN allows the accurate detection of amino acid substitutions at frequencies as low as 1% of the population, a series of additional drug resistance mutations were detected by deep sequencing (mean 2.5, 1.5, and 0.9, respectively). The presence of these low-abundance HIV-1 variants was associated with drug susceptibility, replicative fitness, and coreceptor tropism determined using sensitive phenotypic assays, enhancing the overall burden of resistance to all four antiretroviral drug classes. Further longitudinal studies based on deep sequencing tests will help to clarify (i) the potential impact of minority HIV-1 drug resistant variants in response to antiretroviral therapy and (ii) the importance of the detection of HIV minority variants in the clinical practice.

Deep sequencing: becoming a critical tool in clinical virology

Population (Sanger) sequencing has been the standard method in basic and clinical DNA sequencing for almost 40 years; however, next-generation (deep) sequencing methodologies are now revolutionizing the field of genomics, and clinical virology is no exception. Deep sequencing is highly efficient, producing an enormous amount of information at low cost in a relatively short period of time. High-throughput sequencing techniques have enabled significant contributions to multiples areas in virology, including virus discovery and metagenomics (viromes), molecular epidemiology, pathogenesis, and studies of how viruses to escape the host immune system and antiviral pressures. In addition, new and more affordable deep sequencing-based assays are now being implemented in clinical laboratories. Here, we review the use of the current deep sequencing platforms in virology, focusing on three of the most studied viruses: human immunodeficiency virus (HIV), hepatitis C virus (HCV), and influenza virus.