Rapid and economical drug resistance profiling with Nanopore MinION for clinical specimens with low bacillary burden of Mycobacterium tuberculosis
BMC Research Notes volume 13, Article number: 444 (2020)
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.
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.
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 . 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 . 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 , several groups utilized Nanopore MinION and successfully identified single nucleotide variants (SNVs) in Plasmodium falciparum , dengue virus  and chronic lymphocytic leukemia [10, 11].
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.
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.
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)  and mapped to MTB reference sequence (NC_000962.3) using Burrows-Wheeler Aligner MEM algorithm (UGENE version 1.29.0) . Samtools mpileup (Galaxy version 2.1.4) was used to generate pileups for BAM files . SNVs, insertions and deletions (indels) were detected using VarScan mpileup (Galaxy version 22.214.171.124) 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 . 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) . Sequencing error rate was estimated using Qualimap version 2.2.1 . IGV was used for manual inspection of BAM datasets. Nanopolish variants (Galaxy version 0.1.0) was used for signal-level variant calling .
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%) . 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 , 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 .
MinION sequencing results
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:
All discordant SNVs were wildtype bases by MiSeq, with allelic frequency of at least 99%.
A-G and C-T substitutions were involved, with higher frequency for the latter (Additional file 1: Fig. S1).
From MinION raw data, the ‘dominant’ bases (VAF: 51–95%) were identical to that of MiSeq, but were not called by Nanopolish.
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).
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:
Most insertions were present in eis promoter (n = 9), whereas none was found in inhA promoter and rpoB sequences.
Most deletions were present in katG (n = 15) and absent in gyrA, inhA promoter, pncA and rpsL sequences.
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 , nanocorrect and Amplicon Long-read Error Correction (ALEC) python script  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).
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.
An increased sample size may better estimate the PCR success rate.
We could not correlate the data with second-line DST, which was not routinely performed for non-multidrug-resistant MTB isolates.
As the MinION flow cells possessed suboptimal number of active pores, the sequencing time might be overestimated.
Nanopore sequencing accuracy might be improved with ‘1D2′ chemistry, at the expense of lower sequencing depth .
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 . 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.
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.
Amplicon long-read error correction
Drug susceptibility testing
Integrative genomics viewer
Insertions and deletions
Polymerase chain reaction
Single nucleotide variant
Variant allelic frequency
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.
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.
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.
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.
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.
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.
Mikheyev AS, Tin MM. A first look at the Oxford Nanopore MinION sequencer. Mol Ecol Resour. 2014;14(6):1097–102.
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.
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.
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.
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.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Okonechnikov K, Golosova O, Fursov M, UGENE team. Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics. 2012;28(8):1166–7.
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.
Koboldt DC, Larson DE, Wilson RK. Using VarScan 2 for germline variant calling and somatic mutation detection. Curr Protoc Bioinformati. 2013;44(1):15.
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.
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.
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.
Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018;34(18):3094–100.
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.
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.
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.
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.
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.
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.
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.
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.
Native Barcoding Expansion 96. Oxford Nanopore Technologies. 2020. https://store.nanoporetech.com/native-barcoding-expansion-96.html. Accessed 27 July 2020.
Flongle Flow Cell (R9.4.1). Oxford Nanopore Technologies. 2020. https://store.nanoporetech.com/flowcells/flongle-flow-cell.html. Accessed 27 July 2020.
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.
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.
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.
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.
We thank the colleagues of the Department of Pathology, Hong Kong Sanatorium & Hospital, for their dedicated and professional work on routine laboratory diagnostics.
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
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
About this article
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
- Antibiotic resistance
- Illumina MiSeq
- Nanopore MinION
- Xpert MTB/RIF