Skip to main content

Rapid and economical drug resistance profiling with Nanopore MinION for clinical specimens with low bacillary burden of Mycobacterium tuberculosis

Abstract

Objective

We designed and tested a Nanopore sequencing panel for direct tuberculosis drug resistance profiling. The panel targeted 10 resistance-associated loci. We assessed the feasibility of amplifying and sequencing these loci from 23 clinical specimens with low bacillary burden.

Results

At least 8 loci were successfully amplified from the majority for predicting first- and second-line drug resistance (14/23, 60.87%), and the 12 specimens yielding all 10 targets were sequenced with Nanopore MinION and Illumina MiSeq. MinION sequencing data was corrected by Nanopolish and recurrent variants were filtered. A total of 67,082 bases across all consensus sequences were analyzed, with 67,019 bases called by both MinION and MiSeq as wildtype. For the 41 single nucleotide variants (SNVs) called by MiSeq with 100% variant allelic frequency (VAF), 39 (95.1%) were called by MinION. For the 22 mixed bases called by MiSeq, a SNV with the highest VAF (70%) was called by MinION. With short assay time, reasonable reagent cost as well as continuously improving sequencing chemistry and signal correction pipelines, this Nanopore method can be a viable option for direct tuberculosis drug resistance profiling in the near future.

Introduction

The increasing threat of tuberculosis (TB) drug resistance highlights the importance of prompt drug susceptibility test (DST) results for better patient care and infection control [1, 2]. Nevertheless, culture-dependent methods cannot provide a quick answer due to the fastidious nature of Mycobacterium tuberculosis (MTB). From literature, 24–61% of pulmonary TB cases were acid-fast bacilli (AFB) smear-negative [3, 4], with smear-negative, culture-positive TB accounting for 13% of TB transmission [5]. Novel diagnostic tools are needed for rapid detection of drug resistance from the technically demanding smear-negative specimens.

Recent advent of next-generation sequencing (NGS) has facilitated comprehensive evaluation of MTB genome for drug resistance prediction [6]. Among various options in NGS market, Nanopore sequencers are ideal for infectious disease diagnosis, which requires short sample-to-answer time (Fig. 1). Despite its inferior sequencing accuracy [7], several groups utilized Nanopore MinION and successfully identified single nucleotide variants (SNVs) in Plasmodium falciparum [8], dengue virus [9] and chronic lymphocytic leukemia [10, 11].

Fig. 1
figure1

Assay time comparison between MiSeq, MinION, line probe assay and culture-dependent workflows

The goal of this study was to design and test a Nanopore targeted panel for direct TB drug resistance profiling. First, we attempted to amplify 10 resistance-associated loci from clinical specimens with low MTB burden. Second, we sequenced these amplicons with both Nanopore MinION and Illumina MiSeq, and the sequencing data was collated.

Main text

Methods

Clinical specimens

Twenty-three specimens were collected between August 2016 and May 2017 (Table 1). They were AFB smear-negative or classified as MTB detected low/very low by Xpert MTB/RIF assay (Cepheid, Sunnyvale, CA, USA). TB culture and DST were performed by 2 local reference laboratories, with ethambutol, isoniazid, rifampicin and streptomycin for first-line drug testing.

Table 1 Details of clinical specimens, routine test results and amplification of genomic regions associated with TB drug resistance

DNA extraction

Standard laboratory practice was applied to minimize the risk of infection and contamination. NucliSENS easyMAG automated system (bioMérieux, Marcy, I’Etoile, France) was used for DNA extraction. A maximum of 1-mL pretreated specimen was homogenized in lysis buffer and incubated at 80 °C for 20 min. DNA extraction was performed according to manufacturer’s recommendations, with an elution volume of 25 μL.

Amplification of resistance-associated loci

The 10 targets were located in embB, gyrA, katG, pncA, rpoB, rpsL, rrs, tlyA and promoter regions of eis and inhA genes. They were amplified in separate 25-μL polymerase chain reactions (PCRs), each comprising 1X PCR Buffer II, 1.5 mM MgCl2, 0.2 mM dNTPs, 1 μM forward and reverse primers (Additional file 1: Table S1), 2.5 U AmpliTaq Gold DNA polymerase (Applied Biosystems, Foster City, CA, USA), 1 M betaine (Sigma-Aldrich, St. Louis, MO, USA) and 2 μL of DNA. PCR conditions were as shown in Additional file 1: Table S2. Second PCR was performed if first PCR failed to yield the desired amplicon.

Sequencing by MiSeq

For each specimen, 2-μL aliquots of each amplicon were pooled for Nextera XT DNA library preparation (Illumina, San Diego, CA, USA). Indexed libraries were sequenced using MiSeq v3 Reagent Kit (Illumina, San Diego, CA, USA) in a 250–200 paired-end run.

Sequencing by MinION

Sequencing libraries were prepared using Ligation Sequencing Kit 1D R9 Version (Oxford Nanopore Technologies, Oxford, England). Purified libraries were sequenced on FLO-MIN 106 R9 flow cells (Oxford Nanopore Technologies, Oxford, England).

MiSeq data analysis

Sequencing reads were quality-filtered using Trimmomatic (Galaxy version 0.36.0) [12] and mapped to MTB reference sequence (NC_000962.3) using Burrows-Wheeler Aligner MEM algorithm (UGENE version 1.29.0) [13]. Samtools mpileup (Galaxy version 2.1.4) was used to generate pileups for BAM files [14]. SNVs, insertions and deletions (indels) were detected using VarScan mpileup (Galaxy version 2.4.3.1) with minimum coverage setting of 20, or by manual inspection using Integrative Genomics Viewer (IGV) version 2.4.13 with coverage allele-fraction threshold of 0.1 [15,16,17]. Sequencing coverage was estimated using Tablet version 1.17.08.17 [18]. SNVs and indels were correlated with DST results if available.

MinION data analysis

MinKNOW version 2.0 was used for live basecalling. Adaptor sequences were removed using Porechop (Galaxy version 0.2.3). Trimmed reads were mapped to MTB reference sequence (NC_000962.3) using minimap2 (Galaxy version 2.12) [19]. Sequencing error rate was estimated using Qualimap version 2.2.1 [20]. IGV was used for manual inspection of BAM datasets. Nanopolish variants (Galaxy version 0.1.0) was used for signal-level variant calling [21].

Results and discussion

Routine test results

Nineteen specimens were AFB smear-negative, and the only smear-positive specimen was ‘MTB detected low’ (Patient 19) (Table 1). Twenty specimens were tested with Xpert MTB/RIF assay, with mean threshold cycle (Ct) values of 23.74–34.68, and rifampicin resistance was not detected in 19 but indeterminate for a lung biopsy with the lowest mean Ct value (Patient 2). TB culture was performed for 20 specimens, 16 were positive. Twelve of these isolates were susceptible to all first-line drugs tested, with isoniazid and streptomycin resistance detected in 2 and 3 isolates, respectively. Mean turnaround time of DST was 69.5 days.

Amplification of genomic targets

The targets were amplified in 64.35% of reactions (148/230) (Table 1). Success rate was the highest for embB, gyrA, rpsL and tlyA genes (16/23, 69.57%), while pncA gene was the lowest (12/23, 52.17%). At least 8 loci were successfully amplified from 14 specimens for predicting first- and second-line drug resistance (14/23, 60.87%), which was comparable to an 8-gene pyrosequencing assay for smear-negative specimens (54.9%) [22]. Twelve of these specimens (12/23, 52.17%) had all 10 loci amplified, with mean Ct values of 23.74–30.08. For the 8 specimens that did not yield all targets, the mean Ct values ranged from 28.16 to 34.68. It appeared that a Ct value of 28 or less was required for successful amplification of all targets.

MiSeq sequencing results

A total of 11,350,165 reads (4623 Mb) were generated, averaging 945,847 reads (385.3 Mb) per sample (Additional file 1: Table S3). Mean coverage breadth and depth per sample were 99.85% and 21,441.9, respectively. Minimum depth across all targets was 1488.

A total of 63 SNVs were called, 18 were synonymous, 25 were polymorphic, and 17 with unknown significance to drug resistance (Additional file 1: Table S4). The remaining 3 were missense mutations associated with quinolone, streptomycin and isoniazid resistance (Specimen 14, 16 and 19). No indels were called.

The MTB isolate from Specimen 14 was resistant to streptomycin, with wildtype rpsL and rrs genes. From literature, 25–52% of streptomycin-resistant MTB isolates harbored wildtype rpsL and rrs genes [23,24,25]. Other molecular mechanisms, such as reduced cell permeability to aminoglycosides or presence of drug-modifying enzymes, might contribute to its streptomycin resistance. On the other hand, 281A > G was called in gyrA gene with variant allelic frequency (VAF) of 11%. The resulting D94G substitution is frequently found in fluoroquinolone-resistant strains [26], yet we could not confirm the resistance phenotype as second-line drug susceptibility data was not available.

For Specimen 16, 128A > G was called in rpsL gene. We did not have any streptomycin susceptibility data for interpretation as TB culture was not performed.

For Specimen 19, − 15C > T was called in inhA promoter region, which is prevalent among isoniazid-resistant MTB isolates [27].

MinION sequencing results

General features

Twelve pools of amplicons were sequenced separately on flow cells with 84–624 active pores (Additional file 1: Table S3). First FASTQ files were generated with an elapsed time of 5–86 min, 19.8 min in average, yielding 5.2–5.8 Mb of data and 2866–4685 ‘pass’ reads. Mean coverage breadth and depth across all targets were 100% and 282.8, respectively.

Sequencing error analysis

General alignment error rate was 10.47–15.12%, and mean insertion and deletion error rate (per 100 aligned bases) were 2.09% and 2.83%, respectively. The raw data error rate was comparable to previous studies for chronic lymphocytic leukemia [10, 11].

Comparison with MiSeq sequencing data

Nanopolish was used to improve the accuracy of MinION sequencing data from signal level. For rpsL amplicons of Patient 3 and 21, second and third FASTQ files were also included for data analysis to meet the default minimum depth requirement.

A total of 67,082 ‘nanopolished’ bases across all MinION consensus sequences were compared with MiSeq data (Additional file 1: Table S4 and S5). For the 67,019 wildtype bases called by MiSeq, 66,935 (99.9%) were called as wildtype and 84 (0.1%) as SNVs by MinION. For the 41 SNVs with 100% VAF called by MiSeq, 39 (95.1%) were called by MinION. For the 22 mixed bases called by MiSeq, a SNV (VAF: 70%) was called by MinION, yet the remainder were not called (VAF: 10–37%). There were 32 and 44 insertions and deletions from ‘nanopolished’ data, respectively.

Characteristics of discordant SNVs

The details are summarized in Additional file 1: Table S6, with following observations:

  1. 1.

    All discordant SNVs were wildtype bases by MiSeq, with allelic frequency of at least 99%.

  2. 2.

    A-G and C-T substitutions were involved, with higher frequency for the latter (Additional file 1: Fig. S1).

  3. 3.

    From MinION raw data, the ‘dominant’ bases (VAF: 51–95%) were identical to that of MiSeq, but were not called by Nanopolish.

  4. 4

    From ‘nanopolished’ data, most discordant SNVs were present in eis promoter (n = 49) and absent in gyrA, inhA promoter, rrs and tlyA sequences (Additional file 1: Fig. S2).

  5. 5.

    From ‘nanopolished’ data, most discordant SNVs (80/84, 95.24%) were called recurrently in different specimens, and in 16 sequence patterns (Additional file 1: Fig. S3).

Characteristics of indels

The results are summarized in Additional file 1: Fig. S2 and S3, with following observations:

  1. 1.

    Most insertions were present in eis promoter (n = 9), whereas none was found in inhA promoter and rpoB sequences.

  2. 2.

    Most deletions were present in katG (n = 15) and absent in gyrA, inhA promoter, pncA and rpsL sequences.

  3. 3.

    Most indels (insertions, 75%; deletions, 90.91%) were called recurrently in different specimens (insertions, n = 7; deletions, n = 6).

Filtering recurrent variants

From ‘nanopolished’ data, recurrent variants were called across different samples, which was also observed in other studies using Nanopolish [10], nanocorrect and Amplicon Long-read Error Correction (ALEC) python script [11] for data processing. In both studies, the authors filtered recurrent variants to improve sequencing accuracy. By applying this strategy, the concordance for wildtype bases increased from 99.9% (66,935/67,019) to 100%, and insertions and deletions were reduced from 32 to 8 and 44 to 4, respectively (Additional file 1: Table S5).

Conclusions

We developed a Nanopore targeted panel for direct detection of TB drug resistance, with full-set sequence data retrieved from about half of the specimens with low MTB burden. The assay time was 6–9 h, which improved an average of 69.5 days by DST, and was comparable to commercial methods like line probe assay (Fig. 1). The method was merited by reasonable reagent cost (64 USD per sample for 24-plex workflow), which could be further lowered with the 96-barcode option and Flongle flow cells [28, 29]. With continuously improving sequencing chemistry, more sophisticated signal correction pipelines and the above-mentioned merits, we envision that Nanopore sequencing can be a viable option for ‘end TB’ in the near future.

Limitations

  1. 1.

    An increased sample size may better estimate the PCR success rate.

  2. 2.

    We could not correlate the data with second-line DST, which was not routinely performed for non-multidrug-resistant MTB isolates.

  3. 3.

    As the MinION flow cells possessed suboptimal number of active pores, the sequencing time might be overestimated.

  4. 4.

    Nanopore sequencing accuracy might be improved with ‘1D2′ chemistry, at the expense of lower sequencing depth [30].

  5. 5.

    Nanopolish could not improve the detection of minor variants, which might be caused by high error rate of Nanopore raw data. Sequencing error might originate from inhomogeneous translocation speed of DNA, low signal-to-noise ratio and simultaneous passage of multiple nucleotides through nanopores [31]. Nevertheless, substitution-type miscalls might also be present in MiSeq data, which might arise from similar emission spectra of fluorophores, phasing and pre-phasing phenomena, and index PCR of library preparation [32, 33]. In addition, polymerase error accumulated in target amplification could affect the accuracy of variant calling by both MiSeq and MinION.

  6. 6.

    Therefore, we should be vigilant that the variants called by MiSeq could be false whereas the variants called by MinION could be true in some occasions. It is advisable to confirm the presence of these variants by repeating the target amplification step followed by Sanger sequencing.

Availability of data and materials

The MinION and MiSeq sequencing reads were deposited in NCBI Sequence Read Archive (SRA) under the BioProject accession number PRJNA648373.

Abbreviations

AFB:

Acid-fast bacilli

ALEC:

Amplicon long-read error correction

Ct:

Threshold cycle

DST:

Drug susceptibility testing

IGV:

Integrative genomics viewer

Indels:

Insertions and deletions

MTB:

Mycobacterium tuberculosis

NGS:

Next-generation sequencing

PCR:

Polymerase chain reaction

SNV:

Single nucleotide variant

TB:

Tuberculosis

VAF:

Variant allelic frequency

References

  1. 1.

    Drug-resistant TB: global situation. World Health Organization. 2019. https://www.who.int/tb/areas-of-work/drug-resistant-tb/global-situation/en/. Accessed 3 April 2019.

  2. 2.

    Dookie N, Rambaran S, Padayatchi N, Mahomed S, Naidoo K. Evolution of drug resistance in Mycobacterium tuberculosis: a review on the molecular determinants of resistance and implications for personalized care. J Antimicrob Chemother. 2018;73(5):1138–51.

    CAS  Article  Google Scholar 

  3. 3.

    Zhang ZX, Sng LH, Yong Y, Lin LM, Cheng TW, Seong NH, et al. Delays in diagnosis and treatment of pulmonary tuberculosis in AFB smear-negative patients with pneumonia. Int J Tuberc Lung Dis. 2017;21(5):544–9.

    CAS  Article  Google Scholar 

  4. 4.

    Getahun H, Harrington M, O’Brien R, Nunn P. Diagnosis of smear-negative pulmonary tuberculosis in people with HIV infection or AIDS in resource-constrained settings: informing urgent policy changes. Lancet. 2007;369(9578):2042–9.

    Article  Google Scholar 

  5. 5.

    Tostmann A, Kik SV, Kalisvaart NA, Sebek MM, Verver S, Boeree MJ, et al. Tuberculosis transmission by patients with smear-negative pulmonary tuberculosis in a large cohort in the Netherlands. Clin Infect Dis. 2008;47(9):1135–42.

    Article  Google Scholar 

  6. 6.

    Satta G, Lipman M, Smith GP, Arnold C, Kon OM, McHugh TD. Mycobacterium tuberculosis and whole-genome sequencing: how close are we to unleashing its full potential? Clin Microbiol Infect. 2018;24(6):604–9.

    CAS  Article  Google Scholar 

  7. 7.

    Mikheyev AS, Tin MM. A first look at the Oxford Nanopore MinION sequencer. Mol Ecol Resour. 2014;14(6):1097–102.

    CAS  Article  Google Scholar 

  8. 8.

    Imai K, Tarumoto N, Runtuwene LR, Sakai J, Hayashida K, Eshita Y, et al. An innovative diagnostic technology for the codon mutation C580Y in kelch13 of Plasmodium falciparum with MinION nanopore sequencer. Malar J. 2018;17(1):217.

    Article  Google Scholar 

  9. 9.

    Yamagishi J, Runtuwene LR, Hayashida K, Mongan AE, Thi LAN, Thuy LN, et al. Serotyping dengue virus with isothermal amplification and a portable sequencer. Sci Rep. 2017;7(1):3510.

    Article  Google Scholar 

  10. 10.

    Orsini P, Minervini CF, Cumbo C, Anelli L, Zagaria A, Minervini A, et al. Design and MinION testing of a nanopore targeted gene sequencing panel for chronic lymphocytic leukemia. Sci Rep. 2018;8(1):11798.

    Article  Google Scholar 

  11. 11.

    Minervini CF, Cumbo C, Orsini P, Brunetti C, Anelli L, Zagaria A, et al. TP53 gene mutation analysis in chronic lymphocytic leukemia by nanopore MinION sequencing. Diagn Pathol. 2016;11(1):96.

    Article  Google Scholar 

  12. 12.

    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    CAS  Article  Google Scholar 

  13. 13.

    Okonechnikov K, Golosova O, Fursov M, UGENE team. Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics. 2012;28(8):1166–7.

    CAS  Article  Google Scholar 

  14. 14.

    Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011;27(21):2987–93.

    CAS  Article  Google Scholar 

  15. 15.

    Koboldt DC, Larson DE, Wilson RK. Using VarScan 2 for germline variant calling and somatic mutation detection. Curr Protoc Bioinformati. 2013;44(1):15.

    Google Scholar 

  16. 16.

    Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29(1):24–6.

    CAS  Article  Google Scholar 

  17. 17.

    Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14(2):178–92.

    Article  Google Scholar 

  18. 18.

    Milne I, Stephen G, Bayer M, Cock PJ, Pritchard L, Cardle L, et al. Using Tablet for visual exploration of second-generation sequencing data. Brief Bioinform. 2013;14(2):193–202.

    CAS  Article  Google Scholar 

  19. 19.

    Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018;34(18):3094–100.

    CAS  Article  Google Scholar 

  20. 20.

    García-Alcalde F, Okonechnikov K, Carbonell J, Cruz LM, Götz S, Tarazona S, et al. Qualimap: evaluating next-generation sequencing alignment data. Bioinformatics. 2012;28(20):2678–9.

    Article  Google Scholar 

  21. 21.

    Simpson JT, Workman RE, Zuzarte PC, David M, Dursi LJ, Timp W. Detecting DNA cytosine methylation using nanopore sequencing. Nat Methods. 2017;14(4):407–10.

    CAS  Article  Google Scholar 

  22. 22.

    Georghiou SB, Seifert M, Lin SY, Catanzaro D, Garfein RS, Jackson RL, et al. Shedding light on the performance of a pyrosequencing assay for drug-resistant tuberculosis diagnosis. BMC Infect Dis. 2016;16(1):458.

    Article  Google Scholar 

  23. 23.

    Sreevatsan S, Pan X, Stockbauer KE, Williams DL, Kreiswirth BN, Musser JM. Characterization of rpsL and rrs mutations in streptomycin-resistant Mycobacterium tuberculosis isolates from diverse geographic localities. Antimicrob Agents Chemother. 1996;40(4):1024–6.

    CAS  Article  Google Scholar 

  24. 24.

    Cuevas-Córdoba B, Cuellar-Sánchez A, Pasissi-Crivelli A, Santana-Álvarez CA, Hernández-Illezcas J, Zenteno-Cuevas R. rrs and rpsL mutations in streptomycin-resistant isolates of Mycobacterium tuberculosis from Mexico. J Microbiol Immunol Infect. 2013;46(1):30–4.

    Article  Google Scholar 

  25. 25.

    Villellas C, Aristimuño L, Vitoria MA, Prat C, Blanco S, García de Viedma D, et al. Analysis of mutations in streptomycin-resistant strains reveals a simple and reliable genetic marker for identification of the Mycobacterium tuberculosis Beijing genotype. J Clin Microbiol. 2013;51(7):2124–30.

    CAS  Article  Google Scholar 

  26. 26.

    Maruri F, Sterling TR, Kaiga AW, Blackman A, van der Heijden YF, Mayer C, et al. A systematic review of gyrase mutations associated with fluoroquinolone-resistant Mycobacterium tuberculosis and a proposed gyrase numbering system. J Antimicrob Chemother. 2012;67(4):819–31.

    CAS  Article  Google Scholar 

  27. 27.

    Caws M, Duy PM, Tho DQ, Lan NT, Hoa DV, Farrar J. Mutations prevalent among rifampin- and isoniazid-resistant Mycobacterium tuberculosis isolates from a hospital in Vietnam. J Clin Microbiol. 2006;44(7):2333–7.

    CAS  Article  Google Scholar 

  28. 28.

    Native Barcoding Expansion 96. Oxford Nanopore Technologies. 2020. https://store.nanoporetech.com/native-barcoding-expansion-96.html. Accessed 27 July 2020.

  29. 29.

    Flongle Flow Cell (R9.4.1). Oxford Nanopore Technologies. 2020. https://store.nanoporetech.com/flowcells/flongle-flow-cell.html. Accessed 27 July 2020.

  30. 30.

    Cornelis S, Gansemans Y, Vander Plaetsen AS, Weymaere J, Willems S, Deforce D, et al. Forensic tri-allelic SNP genotyping using nanopore sequencing. Forensic science international. Genetics. 2019;38:204–10.

    CAS  Google Scholar 

  31. 31.

    Rang FJ, Kloosterman WP, de Ridder J. From squiggle to basepair: computational approaches for improving nanopore sequencing read accuracy. Genome Biol. 2018;19(1):90.

    Article  Google Scholar 

  32. 32.

    Schirmer M, Ijaz UZ, D’Amore R, Hall N, Sloan WT, Quince C. Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform. Nucleic Acids Res. 2015;43(6):e37.

    Article  Google Scholar 

  33. 33.

    Pfeiffer F, Gröber C, Blank M, Händler K, Beyer M, Schultze JL, et al. Systematic evaluation of error rates and causes in short samples in next-generation sequencing. Sci Rep. 2018;8(1):10950.

    Article  Google Scholar 

Download references

Acknowledgements

We thank the colleagues of the Department of Pathology, Hong Kong Sanatorium & Hospital, for their dedicated and professional work on routine laboratory diagnostics.

Funding

Not applicable.

Author information

Affiliations

Authors

Contributions

BSFT and WSC conceived and designed the study. BSFT, ESKM, TLC, WSC, CHA, YC, HCML, DNH, EYLW and TWL were involved in data collection and analysis. WSC wrote the first draft. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Bone Siu Fai Tang.

Ethics declarations

Ethics approval and consent to participate

This study was approved and consent to participate was waived by Research Ethics Committee (REC) of Hong Kong Sanatorium & Hospital under the reference number RC-2019-21. No patient-identifying data was collected.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Fig. S1

Types of discordant SNVs from ‘nanopolished’ consensus sequences. Fig. S2 Number of discordant SNVs and indels per gene in ‘nanopolished’ consensus sequences. Fig. S3 Recurrent variants in ‘nanopolished’ consensus sequences. The bases in square brackets denote the positions of SNVs and indels. Table S1 Primers used in this study. Table S2 PCR conditions. Table S3 General features of MiSeq and MinION sequencing runs. Table S4 Single nucleotide variants (SNVs) called by MiSeq and corresponding MinION results. Table S5 Comparison of MiSeq and MinION data. Table S6 Details of discordant SNVs by MinION.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chan, W.S., Au, C.H., Chung, Y. et al. Rapid and economical drug resistance profiling with Nanopore MinION for clinical specimens with low bacillary burden of Mycobacterium tuberculosis. BMC Res Notes 13, 444 (2020). https://doi.org/10.1186/s13104-020-05287-9

Download citation

Keywords

  • Antibiotic resistance
  • Illumina MiSeq
  • MDR-TB
  • Nanopore MinION
  • NGS
  • Tuberculosis
  • XDR-TB
  • Xpert MTB/RIF