Next Article in Journal
Approach to Strain Selection and the Propagation of Viral Stocks for Venezuelan Equine Encephalitis Virus Vaccine Efficacy Testing under the Animal Rule
Next Article in Special Issue
Evolution of BACON Domain Tandem Repeats in crAssphage and Novel Gut Bacteriophage Lineages
Previous Article in Journal
Immunopathologic Effects of Prednisolone and Cyclosporine A on Feline Immunodeficiency Virus Replication and Persistence
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 11
Issue 9
10.3390/v11090806
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MiDRMpol: A High-Throughput Multiplexed Amplicon Sequencing Workflow to Quantify HIV-1 Drug Resistance Mutations against Protease, Reverse Transcriptase, and Integrase Inhibitors

by
Shambhu G. Aralaguppe
1,*,†,
Anoop T. Ambikan
1,†,
Manickam Ashokkumar
1,2,
Milner M. Kumar
3,
Luke Elizabeth Hanna
2,
Wondwossen Amogne
4,
Anders Sönnerborg
1,5 and
Ujjwal Neogi
1
1
Division of Clinical Microbiology, Department of Laboratory Medicine, Karolinska Institutet, Huddinge, 14186 Stockholm, Sweden
2
Department of HIV/AIDS, National Institute for Research in Tuberculosis, Indian Council of Medical Research, Chennai 600031, India
3
Blackbuck Technologies, Porur, Chennai 600116, India
4
Department of Internal Medicine, School of Medicine, Addis Ababa University, 2380 Addis Ababa, Ethiopia
5
Division of Infectious Diseases, Department of Medicine Huddinge, Karolinska Institutet, 14186 Stockholm, Sweden
*
Author to whom correspondence should be addressed.
†
Equal contribution.
Viruses 2019, 11(9), 806; https://doi.org/10.3390/v11090806
Submission received: 15 August 2019 / Accepted: 24 August 2019 / Published: 30 August 2019
(This article belongs to the Special Issue Computational Biology of Viruses: From Molecules to Epidemics)
Browse Figures
Review Reports Versions Notes

Abstract

:
The detection of drug resistance mutations (DRMs) in minor viral populations is of potential clinical importance. However, sophisticated computational infrastructure and competence for analysis of high-throughput sequencing (HTS) data lack at most diagnostic laboratories. Thus, we have proposed a new pipeline, MiDRMpol, to quantify DRM from the HIV-1 pol region. The gag-vpu region of 87 plasma samples from HIV-infected individuals from three cohorts was amplified and sequenced by Illumina HiSeq2500. The sequence reads were adapter-trimmed, followed by analysis using in-house scripts. Samples from Swedish and Ethiopian cohorts were also sequenced by Sanger sequencing. The pipeline was validated against the online tool PASeq (Polymorphism Analysis by Sequencing). Based on an error rate of <1%, a value of >1% was set as reliable to consider a minor variant. Both pipelines detected the mutations in the dominant viral populations, while discrepancies were observed in minor viral populations. In five HIV-1 subtype C samples, minor mutations were detected at the <5% level by MiDRMpol but not by PASeq. MiDRMpol is a computationally as well as labor efficient bioinformatics pipeline for the detection of DRM from HTS data. It identifies minor viral populations (<20%) of DRMs. Our method can be incorporated into large-scale surveillance of HIV-1 DRM.

1. Introduction

Genotypic resistance testing (GRT) to monitor antiretroviral therapy (ART) of human immunodeficiency virus type 1 (HIV-1) infection is challenging due to drug resistance mutations (DRMs) in minor viral populations not detected by Sanger sequencing, although the impact of such DRMs on clinical outcome is still under debate [1]. A major challenge for the implementation of high-throughput sequencing (HTS) in clinical practice is the need of bioinformatics expertise, since there is still no consensus on the standardization of an efficient and easy pipeline for data analysis. However, advances in HTS technology have facilitated sufficient throughput data and sensitivity to detect very rare viral DRMs [2]. On the other hand, HTS-based GRT has a few limitations that need to be addressed, including short and error-prone reads [2], low-abundance interpretation of variant data, and a lack of standardization of data analysis protocols/pipelines [3].
In this study, we evaluated a simplified bioinformatics pipeline, named MiDRMpol, for analysis of HTS data which is cost, labor, and computationally efficient. The Fastq files generated by HTS can be used as input to get easily inferable diagnostic reports and can be applied in large-scale surveillance. It can be used without prior knowledge of HTS or bioinformatics. Here, we restricted the analysis of DRMs to the pol region of the HIV-1 genome, including the targets of the vast majority of presently approved HIV-1 drugs: the reverse transcriptase, protease, and integrase.

2. Materials and Methods

2.1. Clinical Specimens, Plasmids, and Control Virus

Plasma samples were obtained from three cohorts originating from Sweden (n = 60), India (n = 10), and Ethiopia (n = 17). All patients were treatment naïve except 10 Swedish subjects, of whom only 2 were given ART at the time of plasma sampling. The majority (30/50; 60%) of the Swedish treatment-naïve subjects had been infected outside the country. Eight of the 10 Indian samples were from children who were infected through mother-to-child transmission (MTCT). Besides, to deduce the reliability of the method, a standard molecular clone (pMJ4) [4] was used as described in [5].

2.2. High-Throughput Sequencing

The cDNA was prepared using gene-specific primers and the HTS was carried out as described by us [5]. The F1 fragment (Gag-to-vpu) (HXB2, 790-6231) was amplified and gel-purified using the QIAamp gel extraction kit (Qiagen, Venlo, The Netherlands). Briefly, the purified amplicon was fragmented on the Covaris S200 at 300 bp for 75 s with peak power-50 and cycle/burst-200, and the library was prepared using NEBNext® Ultra™ DNA Library Prep Kit for Illumina® (New England Biolabs, Ipswich, MA, USA) with multiplexed NEB next adaptors. The samples were then pooled together in sets of 48 samples along with other unrelated nonviral indexed libraries. Paired-end sequences of read length 250 bp were carried out on the Illumina HiSeq2500.

2.3. Bioinformatics Pipeline for Analysis

The raw reads were adapter-trimmed using the Cutadapt v1.8 program using the default setting (error rate of <0.1) [6], followed by removal of the low-quality bases (Phred value score of <Q30) by Sickle version 1.33 [7]. Duplicate reads were removed using FastUniq [8]. The processed reads were aligned against the gene nucleotide sequence of POL from HIV. This alignment was performed in local mode using Bowtie2 [9] to identify reads originating from the POL gene. Then, these reads were further aligned with the POL protein sequence using the blastx program from Blast package [10]. For each read, the best blast hit was chosen and all amino acid substitutions were counted. Though there were multiple amino acid substitutions at each position in the POL protein sequence, only the most-frequent substitution was considered, and the fraction of the most-frequent amino acid substitution was estimated. The complete script is available at GitHub at https://github.com/neogilab/MiDRMPol.

2.4. Standard Population Sequencing

The ViroSeq HIV-1 genotyping system had been used for the Swedish patients as a part of routine clinical practice. HIV-1 GRT by population sequencing (GRT-PS) was also performed for the Ethiopian patients using an in-house method [11] but not for Indian patients due to insufficient sample availability.

2.5. HIV-1 Subtyping and Phylogenetic Analysis

Consensus pol sequences were generated from the base count data with the occurrence of a particular nucleotide of >50% in a given position and subtyped by REGAv3 [12] and COMET-HIV [13], followed by phylogenetic analysis using MEGA7.

2.6. Drug Resistance Mutations

For treatment-naïve individuals, DRMs was described using definitions for surveillance of transmitted HIV-1 drug resistance, 2009 Update (sDRM2009) [14]. Mutation lists provided in the Stanford University HIV Drug Resistance Database [15], accessed on 6 January 2018 [15,16], was used for describing major integrase mutations, including Q148H, E138K, T66I, and Q148H. The DRMs of treatment-experienced patients were interpreted based on the Stanford University HIV Drug Resistance Database [15].

2.7. Availability of Data

The GenBank IDs for the sequences are available with the following accession numbers: for the Indian cohort: KX069219-228; for the Swedish cohort: MF373124-206; and for the Ethiopian cohort: KP411823-826, 828 and KP411830-845. Fastq files: Submission ID: SUB5359871, BioProject ID: PRJNA529776.

2.8. Ethics Approval and Consent of Participants

Ethical permissions were obtained from the respective sites. Swedish samples: Regional Ethics Committee Stockholm (Dnr: 2006/1367-31/4); Ethiopian samples: the Ethiopian Science and Technology Agency (Ref. No. RPHE/126-83/08), and the Drug Administration and Control Authority of Ethiopia (Ref. No. 02/6/22/17); Indian samples: National Institute for Research in Tuberculosis Institutional Ethics Committee (NIRT IEC No. 2009009). The patients’ information was anonymized and delinked before analysis.

3. Results

3.1. Coverage and Error Rate

The amino acid coverage at each postion of the HIV-pol region was calculated in terms of quartiles. The range of first quartile was 3836x–12,725x whereas for third quartile, it was 15,403×–36,451x. All the DRM positions had median coverage >5000x (Figure 1a). To estimate the intrinsic error introduced by the laboratory workflow, we sequenced pMJ4, where the error rate was found to be <1% at the nucleotide level [5]. In addition, we determined the median error rate at the amino acid level, which was 0.2% (0–1.8%) (Figure 1b).

3.2. Subtyping and Phylogenetic Analysis

Based on the phylogenetic analysis, we found that the majority of the viral isolates (n = 78, 90%) were subtype C, of which 10 isolates were from treatment-experienced individuals. Seven (8%) isolates were subtype B. In addition, we identified one HIV-1A1 and one 01_AE. No cross contamination was found (Figure 2). Two Swedish samples (SE70 and SE73) clustered together and were epidemiologically linked.

3.3. Identification of DRM by GRT-HTS

DRMs were identified and quantified in 23/87 patients (Table 1). In the Swedish cohort, 9/50 (18%) therapy-naïve patients had DRMs (M46I: n = 3) (M184I, T215S, or K219R: n = 3); major--integrase strand transfer inhibitors (INSTI) mutation Q148H or E138K: n = 2)). Among Indian patients, 5/10 (50%) had reverse transcriptase inhibitors (RTI) DRM, and all 5 were MTCT patients. No DRMs were detected in horizontally infected patients. Among Ethiopian patients, 6/17 (35%) had DRMs, of whom 5 had a single-class DRM and one had both nucleoside reverse transcriptase inhibitors (NRTI) (M184I) and INSTI (Q148H) DRM.
Among the 10 treatment-experienced patients from the Swedish cohort, only two patients were on ART at the time of plasma sampling, although they had high viral RNA levels, indicating that they did not take the drugs and no DRMs were found. In three patients, DRMs were found. One patient (SE03; HIV-1C) who had terminated nevirapine, tenofovir, and raltegravir therapy one month earlier had D67G and L74I. This patient had earlier been given abacavir-containing therapy. Patients SE70 and SE73 (HIV-1C) had stopped cART ABC + DDI + EFV and d4T + 3TC + NVP, respectively, a few months before plasma sampling. These subjects were epidemiologically linked and were identified with the same set of RTI mutations: D67N, L210W, T215Y, K103N, and G190A.

3.4. Comparison of DRM Detected by GRT-PS and GRT-HTS

While comparing GRT-PS and GRT-HTS among the Swedish and Ethiopian patients, additional DRMs were detected by GRT-HTS (Table 1). Among the treatment-naïve individuals (n = 77), GRT-HTS detected DRMs in 26% (20/77) of the samples (e.g., M184I (n = 4), T125S (n = 1), and M230L (n = 1)), whereas GRT-PS detected them in only one naïve patient (M230L).
In treatment-experienced patients (n = 10), DRMs were detected in two and three samples by GRT-PS and GRT-HTS, respectively. In one sample (SE73), Y188YC was detected by GRT-PS and not by GRT-HTS. One of the Indian MTCT samples (IN07) had several major NRTI and NNRTI DRMs at frequencies of <20%. GRT-HTS also identified major INSTI mutations at low frequencies (e.g., Q148H, E138K, and T66I), which were undetectable by GRT-PS. In two samples (SE73 and SE31), GRT-PS did not detect DRMs which were >20%.
When comparing MiDRMpol and PASeq outputs for 23 samples, the same DRMs with comparable frequencies were reported for 13 sequences (Table 1). In six samples (SE03, SE07, ET115, ET122, ET159, and ET171), only MiDRMpol reported minor DRMs (frequency: 1.0%–4.8%). In four samples, PASeq did not deliver results due to the large file size (the file upload limit of PASeq is 1 GB).
In detail, the discrepancy was due to M184I (1.0%–1.5%), which is considered to be a DRM in MiDRMpol but not in PASeq, in three samples (ET115, ET159, and ET171). In sample SE03, D67G and L74I (2.0% and 1.9%, respectively); in samples SE07 and ET171, Q148H (1.1% and 1.2%, respectively); and in sample ET122, M230L (4.8%) were reported by MiDRMpol but not by PASeq. This is possibly because of the FastUniq tool used in MiDRMpol, which removes duplicate reads, due to which the minor populations of mutations are enriched.
In sample 24C, T215A (99.7%) and T69N (41.7%) were not detected in our pipeline because we do not consider T215A and T69N as resistance mutations against any of the RTIs. One of the mutations (D67E) was detected at 77.1% in PASeq but not detected in our pipeline.

4. Discussion

In this study, our in-house pipeline, MiDRMpol, was evaluated for the identification of DRMs in sequences obtained by HTS. Patients from Swedish, Indian, and Ethiopian cohorts, representing diverse subtypes, were included. The MiDRMpol pipeline was successful in detecting all major and minor DRMs. Identification of mutations such as M184I, which appears at the very early stage of the evolution of M184V, and T215S, a revertant mutation, by GRT-HTS shows that it can detect the early evolution of DRM and earlier exposure of cART, which can provide further information about, for example, pretreatment drug resistance (PDR) in global surveillances of the spread of resistant HIV-1. Although GRT-HTS detects HIV DRMs more sensitively than GRT-PS, the generalized application of such methods for clinical DR monitoring is not currently feasible due to the high cost per sample. However, the cost per megabase of raw data of DNA sequences has decreased 370-fold, from US$5.200 in 2001, to US$0.014 in 2015 [17]. Equipment and maintenance costs, however, remain unaffordable for many low-income countries. The most important limitation, however, is the lack of automated, validated, and robust but simplified bioinformatic analyses coupled with HIV-1 resistance interpretations to enable NGS use and interpretation by laboratory technicians, but even this is improving rapidly [18]. Currently, for our approach, the cost per sample was US$130–US$150 when 24 samples were pooled in a single run [5]. This can be further reduced by multiplexing a higher number of samples. MiDRMpol would cut the costs down by excluding the necessity of bioinformaticians to analyze HTS data.
DRMs were identified 8 times more often by GRT-HTS than by GRT-PS, which is in agreement with another study, which observed that 30%–50% of transmitted drug resistance (TDR) mutations were not identified by GRT-PS [19,20]. An investigation of the established pipelines PASeq and HyDRA showed minimal advantage when using HTS data over Sanger sequencing in identifying low-frequency DRMs (≥5%) in virological failure individuals treated with integrase inhibitors [21]. In our earlier study of an Ethiopian cohort, we identified an increased risk of treatment failure in patients who had pretherapy minor RTI mutations [22]. To the best of our knowledge, no bioinformatics pipeline interprets the success of ART. It is important that the outcome of the result generated through the pipeline is used to predict the treatment outcomes and guide patient management in clinical practice [23]. The impact of drug-resistant variants does not depend only on the frequency; rather, the mutational load (Equation (1)) determines more precisely the absolute burden of the resistant variants [24,25]. Thus, at a high viral load, even <1% of minor NNRTI DRM variants have been claimed to have the potential to overcome the drug pressure [26,27].
M u t a t i o n i a l   l o a d   =   H I V   R N A   c o p i e s   p e r   m L   × M u t a n t   v a r i a n t   f r e q u e n c y .
However, the pipeline MiDRMpol was designed to identify >1% of minor DRM populations, independent of the subtype. Although 10 of the patients from Sweden were treatment experienced, only 2 of them claimed to be undergoing treatment but had most likely not taken the drugs. Only 3 of the 10 patients had detectable DRMs despite earlier treatment failures. It should, however, be noted that ART had been terminated just a few months before sampling, suggesting that the DRMs had not yet had time to disappear. In other patients, the treatment had been terminated several months or years before. Thus, although no precise analysis can be done due to lack of samples during earlier treatment failure, our data indicate that even the use of a very sensitive HTS may fail to detect DRMs which had developed at earlier treatment failures.
In the Indian cohort, we identified minor DRMs in 2 out of 10 treatment-naïve individuals who were vertically infected. Eight minor RTI DRMs, including K65R (3.0%), M184V (7.5%), and V106M (3.9%), were detected in one infant (IN07), as a result of the transmission of the resistant virus from the mother. Thus, our data show that GRT-HTS improves detection of vertically transmitted DRMs in infants [28].
Similar algorithms for determining the frequency of DRM variants are available, for example, PASeq [29], a web service which uses various tools, and CoVaMa (Co-Variation Mapper), which is scripted in Python and takes HTS alignment data and populates large matrices of contingency tables [30]. Other available automated pipelines for analysis of sequences obtained through HTS have been summarized in a comparative study [31]. We compared our pipeline with PASeq and observed concordance at the DRM level, although our pipeline also detected M184I mutations in treatment-naïve individuals at minor populations.

5. Conclusions

In this study, we developed an efficient and simplified high-throughput amplicon sequence analysis pipeline named MiDRMpol, which does not require onsite bioinformatics expertise. This pipeline can detect minor DRMs in viral populations at a >1% cutoff and was successfully evaluated in three cohorts consisting of various subtypes. The reliability was analyzed by comparing the sequences generated through both GRT-HTS and GRT-PS. In all of our cohorts, MiDRMpol identified DRMs against the three drug classes in minor viral quasi-species at analysis of sequences obtained through HTS. The complete script is available at GitHub at https://github.com/neogilab/MiDRMPol. Our simplified online pipeline for HTS data of the HIV-1 pol gene would be a valuable tool in the use of diagnostics of drug resistance variants before initiation of ART and after ART failure in optimizing subsequent regimens.

Author Contributions

Conceptualization, U.N.; methodology, U.N., S.G.A., and M.A.; software, A.T.A. and M.M.K.; validation, S.G.A., M.A., and A.T.A.; formal analysis, U.N., S.G.A., and A.T.A.; investigation, S.G.A. and M.A.; resources, L.E.H., W.A., and A.S.; data curation, A.T.A.; writing—original draft preparation, S.G.A.; writing—review and editing, S.G.A., A.T.A., M.A., M.M.K., L.E.H., W.A., A.S., and U.N.; visualization, S.G.A., A.T.A., and U.N.; supervision, U.N.; project administration, U.N.; funding acquisition, U.N.

Funding

The study was funded by the Swedish Research Council (2017-01330, UN; 2016-01675, AS), Karolinska Institutet Doctoral Student Funding (KID2015-154), the European Union´s Horizon 2020 research and innovation program under grant agreement no. 825673, and the Stockholm County Council (ALF 20160074, AS).

Acknowledgments

The following reagent was obtained through the NIH AIDS Reagent Program, Division of AIDS, NIAID, NIH: HIV-1 MJ4 Infectious Molecular Clone (pMJ4) from Drs. Thumbi Ndung’u, Boris Renjifo, and Max Essex. M.A. would like to acknowledge ICMR, India for providing the senior research fellowship.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Inzaule, S.C.; Hamers, R.L.; Noguera-Julian, M.; Casadella, M.; Parera, M.; Kityo, C.; Steegen, K.; Naniche, D.; Clotet, B.; Rinke de Wit, T.F.; et al. Clinically relevant thresholds for ultrasensitive hiv drug resistance testing: A multi-country nested case-control study. Lancet HIV 2018, 5, e638–e646. [Google Scholar] [CrossRef]
  2. Posada-Cespedes, S.; Seifert, D.; Beerenwinkel, N. Recent advances in inferring viral diversity from high-throughput sequencing data. Virus Res. 2017, 239, 17–32. [Google Scholar] [CrossRef] [PubMed]
  3. Ji, H.; Enns, E.; Brumme, C.J.; Parkin, N.; Howison, M.; Lee, E.R.; Capina, R.; Marinier, E.; Avila-Rios, S.; Sandstrom, P.; et al. Bioinformatic data processing pipelines in support of next-generation sequencing-based hiv drug resistance testing: The winnipeg consensus. J. Int. AIDS Soc. 2018, 21, e25193. [Google Scholar] [CrossRef] [PubMed]
  4. Ndung’u, T.; Renjifo, B.; Essex, M. Construction and analysis of an infectious human immunodeficiency virus type 1 subtype c molecular clone. J. Virol. 2001, 75, 4964–4972. [Google Scholar] [CrossRef] [PubMed]
  5. Aralaguppe, S.G.; Siddik, A.B.; Manickam, A.; Ambikan, A.T.; Kumar, M.M.; Fernandes, S.J.; Amogne, W.; Bangaruswamy, D.K.; Hanna, L.E.; Sonnerborg, A.; et al. Multiplexed next-generation sequencing and de novo assembly to obtain near full-length hiv-1 genome from plasma virus. J. Virol. Methods 2016, 236, 98–104. [Google Scholar] [CrossRef] [PubMed]
  6. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  7. Joshi, N.; Fass, J. Sickle: A sliding-window, adaptive, quality-based trimming tool for fastq files (version 1.33) [software]. 2011. Available online: https://github.com/najoshi/sickle (accessed on 2 January 2019).
  8. Xu, H.; Luo, X.; Qian, J.; Pang, X.; Song, J.; Qian, G.; Chen, J.; Chen, S. Fastuniq: A fast de novo duplicates removal tool for paired short reads. PLoS ONE 2012, 7, e52249. [Google Scholar] [CrossRef] [PubMed]
  9. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with bowtie 2. Nat Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed]
  10. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  11. Amogne, W.; Bontell, I.; Grossmann, S.; Aderaye, G.; Lindquist, L.; Sonnerborg, A.; Neogi, U. Phylogenetic analysis of ethiopian hiv-1 subtype c near full-length genomes reveals high intrasubtype diversity and a strong geographical cluster. AIDS Res. Hum. Retroviruses 2016, 32, 471–474. [Google Scholar] [CrossRef]
  12. Regav3. Available online: http://dbpartners.stanford.edu:8080/RegaSubtyping/stanford-hiv/typingtool/ (accessed on 30 March 2016).
  13. Comet-HIV. Available online: https://comet.lih.lu/ (accessed on 30 March 2016).
  14. Bennett, D.E.; Camacho, R.J.; Otelea, D.; Kuritzkes, D.R.; Fleury, H.; Kiuchi, M.; Heneine, W.; Kantor, R.; Jordan, M.R.; Schapiro, J.M.; et al. Drug resistance mutations for surveillance of transmitted hiv-1 drug-resistance: 2009 update. PLoS ONE 2009, 4, e4724. [Google Scholar] [CrossRef] [PubMed]
  15. Stanford hiv db. Available online: https://hivdb.stanford.edu/assets/media/resistance-mutation-handout-Feb2017.516aee6f.pdf (accessed on 2 January 2018).
  16. Rhee, S.Y.; Gonzales, M.J.; Kantor, R.; Betts, B.J.; Ravela, J.; Shafer, R.W. Human immunodeficiency virus reverse transcriptase and protease sequence database. Nucleic Acids Res. 2003, 31, 298–303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. DNA sequencing costs: Data from the nhgri genome sequencing program (gsp). Available online: https://www.genome.gov/about-genomics/fact-sheets/DNA-Sequencing-Costs-Data (accessed on 13 August 2019).
  18. Casadella, M.; Paredes, R. Deep sequencing for hiv-1 clinical management. Virus Res. 2017, 239, 69–81. [Google Scholar] [CrossRef] [PubMed]
  19. Ekici, H.; Rao, S.D.; Sonnerborg, A.; Ramprasad, V.L.; Gupta, R.; Neogi, U. Cost-efficient hiv-1 drug resistance surveillance using multiplexed high-throughput amplicon sequencing: Implications for use in low- and middle-income countries. J. Antimicrob. Chemother. 2014, 69, 3349–3355. [Google Scholar] [CrossRef] [PubMed]
  20. Arias, A.; López, P.; Sánchez, R.; Yamamura, Y.; Rivera-Amill, V. Sanger and next generation sequencing approaches to evaluate hiv-1 virus in blood compartments. Int. J. Environ. Res. Public Health 2018, 15, 1697. [Google Scholar] [CrossRef]
  21. Dalmat, R.R.; Makhsous, N.; Pepper, G.G.; Magaret, A.; Jerome, K.R.; Wald, A.; Greninger, A.L. Limited marginal utility of deep sequencing for hiv drug resistance testing in the age of integrase inhibitors. J. Clin. Microbiol. 2018, 56. [Google Scholar] [CrossRef]
  22. Telele, N.F.; Kalu, A.W.; Gebre-Selassie, S.; Fekade, D.; Abdurahman, S.; Marrone, G.; Neogi, U.; Tegbaru, B.; Sonnerborg, A. Pretreatment drug resistance in a large countrywide ethiopian hiv-1c cohort: A comparison of sanger and high-throughput sequencing. Sci. Rep. 2018, 8, 7556. [Google Scholar] [CrossRef]
  23. Gianella, S.; Richman, D.D. Minority variants of drug-resistant hiv. J. Infect. Dis. 2010, 202, 657–666. [Google Scholar] [CrossRef]
  24. Kozal, M.J.; Chiarella, J.; St John, E.P.; Moreno, E.A.; Simen, B.B.; Arnold, T.E.; Lataillade, M. Prevalence of low-level hiv-1 variants with reverse transcriptase mutation k65r and the effect of antiretroviral drug exposure on variant levels. Antivir. Ther. 2011, 16, 925–929. [Google Scholar] [CrossRef]
  25. Lataillade, M.; Chiarella, J.; Yang, R.; DeGrosky, M.; Uy, J.; Seekins, D.; Simen, B.; St John, E.; Moreno, E.; Kozal, M. Virologic failures on initial boosted-pi regimen infrequently possess low-level variants with major pi resistance mutations by ultra-deep sequencing. PLoS ONE 2012, 7, e30118. [Google Scholar] [CrossRef]
  26. Li, J.Z.; Paredes, R.; Ribaudo, H.J.; Svarovskaia, E.S.; Metzner, K.J.; Kozal, M.J.; Hullsiek, K.H.; Balduin, M.; Jakobsen, M.R.; Geretti, A.M.; et al. Low-frequency hiv-1 drug resistance mutations and risk of nnrti-based antiretroviral treatment failure: A systematic review and pooled analysis. JAMA J. Am. Med Assoc. 2011, 305, 1327–1335. [Google Scholar] [CrossRef] [PubMed]
  27. Li, J.Z.; Paredes, R.; Ribaudo, H.J.; Kozal, M.J.; Svarovskaia, E.S.; Johnson, J.A.; Geretti, A.M.; Metzner, K.J.; Jakobsen, M.R.; Hullsiek, K.H.; et al. Impact of minority nonnucleoside reverse transcriptase inhibitor resistance mutations on resistance genotype after virologic failure. J. Infect. Dis. 2013, 207, 893–897. [Google Scholar] [CrossRef] [PubMed]
  28. Fisher, R.G.; Smith, D.M.; Murrell, B.; Slabbert, R.; Kirby, B.M.; Edson, C.; Cotton, M.F.; Haubrich, R.H.; Kosakovsky Pond, S.L.; Van Zyl, G.U. Next generation sequencing improves detection of drug resistance mutations in infants after pmtct failure. J. Clin. Virol. 2015, 62, 48–53. [Google Scholar] [CrossRef] [PubMed]
  29. Available online: https://www.paseq.org/ (accessed on 2 January 2019).
  30. Routh, A.; Chang, M.W.; Okulicz, J.F.; Johnson, J.E.; Torbett, B.E. Covama: Co-variation mapper for disequilibrium analysis of mutant loci in viral populations using next-generation sequence data. Methods 2015, 91, 40–47. [Google Scholar] [CrossRef] [PubMed]
  31. Noguera-Julian, M.; Edgil, D.; Harrigan, P.R.; Sandstrom, P.; Godfrey, C.; Paredes, R. Next-generation human immunodeficiency virus sequencing for patient management and drug resistance surveillance. J. Infect. Dis. 2017, 216, S829–S833. [Google Scholar] [CrossRef] [PubMed]
Viruses 11 00806 g001 550
Figure 1. Coverage plot and error calculation: (a) The median coverage in each of the positions is indicated by the a dark green line. The light-green shade is the range of the coverage. (b) The error was calculated from the amino acid sequences of pMJ4 (NCBI# AF321523). Red dots indicate the drug resistance mutation (DRM) positions.
Figure 1. Coverage plot and error calculation: (a) The median coverage in each of the positions is indicated by the a dark green line. The light-green shade is the range of the coverage. (b) The error was calculated from the amino acid sequences of pMJ4 (NCBI# AF321523). Red dots indicate the drug resistance mutation (DRM) positions.
Viruses 11 00806 g001
Viruses 11 00806 g002 550
Figure 2. The phylogenetic analysis of the gag-pol sequences. A total of 100 sequences were used in the analysis (87 samples and 13 consensus sequences of gag-vpu from diverse HIV-1 subtypes downloaded from HIVdb). Green circles represent the Swedish cohort (n = 60), blue triangles represent the Indian cohort (n = 10), and red squares represent the Ethiopian cohort (n = 17). The reference sequences of gag-pol from various genetic subtypes—A1, A2, B, C, D, F1, F2, G, H, AE, and AG—were downloaded from the Los Alamos HIV database. After gap-stripping, there were a total of 3096 nucleic acid positions in the final dataset. The bootstrap consensus tree inferred from 1000 replicates is taken to represent the evolutionary history of the taxa analyzed. The tree was drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the maximum composite likelihood method and are in the units of the number of base substitutions per site.
Figure 2. The phylogenetic analysis of the gag-pol sequences. A total of 100 sequences were used in the analysis (87 samples and 13 consensus sequences of gag-vpu from diverse HIV-1 subtypes downloaded from HIVdb). Green circles represent the Swedish cohort (n = 60), blue triangles represent the Indian cohort (n = 10), and red squares represent the Ethiopian cohort (n = 17). The reference sequences of gag-pol from various genetic subtypes—A1, A2, B, C, D, F1, F2, G, H, AE, and AG—were downloaded from the Los Alamos HIV database. After gap-stripping, there were a total of 3096 nucleic acid positions in the final dataset. The bootstrap consensus tree inferred from 1000 replicates is taken to represent the evolutionary history of the taxa analyzed. The tree was drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the maximum composite likelihood method and are in the units of the number of base substitutions per site.
Viruses 11 00806 g002
Table
Table 1. Comparison of DRMs identified by high-throughput sequencing or population sequencing as well as in PASeq and MiDRMpol outputs.
Table 1. Comparison of DRMs identified by high-throughput sequencing or population sequencing as well as in PASeq and MiDRMpol outputs.
CohortPIDTreatmentMiDRMpolPASeqGRT-PS
PINRTINNRTIINIPINRTINNRTIININRTINNRTI
SwedishSE03ExperiencedNoneD67G (2.0%)NoneNoneNoneNoneK103R (16.9%)NoneNoneNone
L74I (1.9%)
SE70ExperiencedNoneD67N (99.7%)K103N (99.8%)NoneNoneD67N (99.9%)K103N (99.9%)NoneD67NK103N
L210W (99.5%)G190A (99.6%) L210W (99.8%)G190A (99.7%) L210WG190A
T215Y (99.4%) T215Y (99.8%) T215Y
SE73ExperiencedNoneD67N (44.2%)K103N (95.9%)NoneNoneD67N (49.2%)K103N (96.1%)NoneNoneK103N
L210W (23.4%)G190A (35.2%) L210W (46.1%)G190A (49.6%) Y188C
T215Y (19.7%) T215Y (45.9%)
19CNaïveM46I (4.6%)NoneNoneNoneFailedNoneNone
20CNaïveNoneT215S (2.1%)NoneNoneFailedNoneNone
26CNaïveNoneM184I (1.1%)NoneNoneFailedNoneNone
05CNaïveNoneK219R (2.3%)NoneNoneFailedNoneNone
SE07NaïveNoneNoneNoneQ148H (1.1%)NoneNoneNoneNoneNoneNone
SE14NaïveNoneNoneNoneE138K (5.5%)NoneNoneNoneE138K (1.6%)NoneNone
SE31NaïveM46I (30.2%)NoneNoneNoneM46I (28.9%)NoneNoneNoneNoneNone
SE97NaïveD30N (13.7%)NoneNoneNoneD30N (12.2%)NoneNoneNoneNoneNone
24CNaïveNoneT69D (99.7%)M230L (99.7%)NoneNoneD67E (77.2%)M230L (99.8%)NoneT69DM230L
T69D (98.1%)
T215A (99.4%)
IndianIN04NaïveNoneNoneK101E (99.7%)NoneNoneNoneK101E (99.9%)NoneNANA
IN05NaïveNoneV75M (7.8%)NoneNoneNoneV75M (8.6%)NoneNoneNANA
IN06NaïveNoneNoneV106M (98.1%)NoneNoneNoneV106M (100.0%)NoneNANA
Y181C (22.8%) Y181C (23.6%)
IN07NaïveNoneK65R (3.0%)V106M (3.9%)NoneNoneK65R (2.8%)V106M (3.6%)NoneNANA
K70E (3.4%)G190S (95.6%) K70E (3.2%)G190S (96.1%)
L74V (66.1%) L74V (67.6%)
M184V (7.5%) M184V (8.4%)
M184I (3.7%) M184I (3.2%)
T215I (3.4%) T215I (3.3%)
IN08NaïveNoneNoneY181C (99.5%)NoneNoneNoneY181C (99.9%)NoneNANA
EthiopianET160NaïveL76V (4.5%)NoneNoneNoneL76V (4.4%)NoneNoneNoneNoneNone
ET115NaïveNoneM184I (1.1%)NoneNoneNoneNoneNoneNoneNoneNone
ET122NaïveNoneNoneM230L (4.8%)NoneNoneNoneNoneNoneNoneNone
ET155NaïveNoneNoneNoneT66I (3.3%)NoneNoneNoneT66I (5.8%)NoneNone
ET159NaïveNoneM184I (1.0%)NoneNoneNoneNoneNoneNoneNoneNone
ET171NaïveNoneM184I (1.5%)NoneQ148H (1.2%)NoneNoneNoneNoneNoneNone
None → no DRM detected; Failed → failure due to upload limit of HTS file; NA → not available.

Share and Cite

MDPI and ACS Style

Aralaguppe, S.G.; Ambikan, A.T.; Ashokkumar, M.; Kumar, M.M.; Hanna, L.E.; Amogne, W.; Sönnerborg, A.; Neogi, U. MiDRMpol: A High-Throughput Multiplexed Amplicon Sequencing Workflow to Quantify HIV-1 Drug Resistance Mutations against Protease, Reverse Transcriptase, and Integrase Inhibitors. Viruses 2019, 11, 806. https://doi.org/10.3390/v11090806

AMA Style

Aralaguppe SG, Ambikan AT, Ashokkumar M, Kumar MM, Hanna LE, Amogne W, Sönnerborg A, Neogi U. MiDRMpol: A High-Throughput Multiplexed Amplicon Sequencing Workflow to Quantify HIV-1 Drug Resistance Mutations against Protease, Reverse Transcriptase, and Integrase Inhibitors. Viruses. 2019; 11(9):806. https://doi.org/10.3390/v11090806

Chicago/Turabian Style

Aralaguppe, Shambhu G., Anoop T. Ambikan, Manickam Ashokkumar, Milner M. Kumar, Luke Elizabeth Hanna, Wondwossen Amogne, Anders Sönnerborg, and Ujjwal Neogi. 2019. "MiDRMpol: A High-Throughput Multiplexed Amplicon Sequencing Workflow to Quantify HIV-1 Drug Resistance Mutations against Protease, Reverse Transcriptase, and Integrase Inhibitors" Viruses 11, no. 9: 806. https://doi.org/10.3390/v11090806

APA Style

Aralaguppe, S. G., Ambikan, A. T., Ashokkumar, M., Kumar, M. M., Hanna, L. E., Amogne, W., Sönnerborg, A., & Neogi, U. (2019). MiDRMpol: A High-Throughput Multiplexed Amplicon Sequencing Workflow to Quantify HIV-1 Drug Resistance Mutations against Protease, Reverse Transcriptase, and Integrase Inhibitors. Viruses, 11(9), 806. https://doi.org/10.3390/v11090806

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop