Rapid Identification of Drug-Resistant Tuberculosis Genes Using Direct PCR Amplification and Oxford Nanopore Technology Sequencing

Zhao, Kaishun; Tu, Chunlin; Chen, Wei; Liang, Haiying; Zhang, Wenjing; Wang, Yilei; Jin, Ye; Hu, Jianrong; Sun, Yameng; Xu, Jun; Yu, Yanfang

doi:https://doi.org/10.1155/2022/7588033

Canadian Journal of Infectious Diseases and Medical Microbiology

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 7588033 | https://doi.org/10.1155/2022/7588033

Rapid Identification of Drug-Resistant Tuberculosis Genes Using Direct PCR Amplification and Oxford Nanopore Technology Sequencing

Kaishun Zhao,¹Chunlin Tu,¹Wei Chen,²Haiying Liang,²Wenjing Zhang,¹Yilei Wang,¹Ye Jin,²Jianrong Hu,²Yameng Sun,³Jun Xu,³and Yanfang Yu¹

Academic Editor: Keke Zhang

Received29 Oct 2021

Revised27 Jan 2022

Accepted14 Mar 2022

Published28 Mar 2022

Abstract

Mycobacterium tuberculosis antimicrobial resistance has been continually reported and is a major public health issue worldwide. Rapid prediction of drug resistance is important for selecting appropriate antibiotic treatments, which significantly increases cure rates. Gene sequencing technology has proven to be a powerful strategy for identifying relevant drug resistance information. This study established a sequencing method and bioinformatics pipeline for resistance gene analysis using an Oxford Nanopore Technologies sequencer. The pipeline was validated by Sanger sequencing and exhibited 100% concordance with the identified variants. Turnaround time for the nanopore sequencing workflow was approximately 12 h, facilitating drug resistance prediction several weeks earlier than that of traditional phenotype drug susceptibility testing. This study produced a customized gene panel assay for rapid bacterial identification via nanopore sequencing, which improves the timeliness of tuberculosis diagnoses and provides a reliable method that may have clinical application.

1. Introduction

Tuberculosis (TB) is one of the top 10 causes of death worldwide and the most common cause of death from a single infectious agent, ranking above human immunodeficiency virus/acquired immunodeficiency syndrome (HIV/AIDS) [1]. There are more than 9 million new cases of TB annually and 2 million deaths [2]. The worldwide percentage of people who developed TB in China is 8.4%, which ranks third behind India (26%) and Indonesia (8.5%) [1]. Drug-resistant TB remains a major threat to public health. In China, approximately 5.7% of new cases of TB and 25.6% of previously treated cases are multidrug-resistant (MDR) [3].

Culture-based drug susceptibility testing (DST) methods using solid or liquid media are currently the gold standard for detecting drug resistance, but it often takes several days or weeks to obtain results [4, 5]. Moreover, culture-based DST depends on culture in laboratories, and thus require a rigorous laboratory environment. As parts of drug resistance loci are distributed in the bacterial genome, molecular tests have become increasingly popular for detection of drug-resistant TB. Engström et al. [6] and Campbell et al. [7] developed pyrosequencing and ABI 3130xl-based sequencing methods, respectively, to detect resistant mutations for several drugs. Whole genome sequencing (WGS) has also been used to genotype potentially resistant single nucleotide polymorphisms (SNPs) [8]. However, this method is dependent on laboratory conditions that require a high investment to purchase a sequencer and level of bioinformatic knowledge. Consistent with the above noted rapid molecular tests, target region enrichment is usually suitable for analysis of low abundance nucleic acids, and amplicon sequencing can rapidly provide comprehensive information on drug resistance across multiple gene regions or multiple sites in the same gene. With the critical importance of rapid drug resistance analysis, there is an urgent need to develop new molecular techniques. The nanopore sequencing platform from Oxford Nanopore Technologies (ONT) offers real-time sequencing and a rapid processing time [9, 10]. The relatively low equipment cost and lack of laboratory requirements and experience makes it an attractive approach.

To date, several studies have examined the capabilities of nanopore platforms with respect to antimicrobial resistance (AMR) detection. The portable MinION sequencing device developed by ONT was used to determine the gene organization of the AMR cassette in Salmonella typhi [11]. A recent study by Golparian et al. [12] used a nanopore sequencer to sequence Neisseria gonorrhoeae and predict susceptibility and resistance to assist in recommending appropriate antimicrobials. The sequencing capacity of nanopore sequencers has reached 450 bp/s, making this sequencing technique and its speed suitable for real-time Klebsiella pneumonia resistome detection and a means to accurately estimate gene expression levels [13].

In the current study, we designed an ONT nanopore sequencer-based system for amplicon sequencing to rapidly analyze seven drug resistance regions as a means of detecting TB drug resistance in different sputum isolates. A convenient clinical report was used to simply present raw sequence datasets and AMR information. The feasibility and accuracy of this method was verified using Sanger sequencing as the standard. The ability to obtain AMR results directly from clinical specimens provides greater open access and is faster than culture-based methods.

2. Materials and Methods

2.1. Clinical Specimens and DNA Extraction

Two hundred sputum specimens were collected from patients diagnosed with pulmonary TB at the Tuberculosis Laboratory of Shanghai Jiading Central Hospital between January 2019 and December 2020. We randomly selected 20 of these 200 specimens for use in this study. The selected specimen included 15 collected from men and 5 collected from women, aged 24–55 years (see Supplementary Material 1). All specimens were stored at −70°C until processed for DNA extraction.

Genomic DNA (gDNA) was extracted from the clinical specimens after pretreatment [14, 15]. Before starting library preparation, the gDNA samples were quantified using Qubit 3.0 fluorometer (Life Technologies, USA) and their purity estimated using a Nanodrop spectrophotometer (Thermo Fisher Scientific, USA). Quality requirements included OD 260/280 ratios of approximately 1.8 and OD 260/230 values of 2.0–2.2.

2.2. Selection of Multiplex Polymerase Chain Reaction (MPCR) Region

To decipher the genetic variation for resistance to the majority of first- and second-generation drugs, a selected panel of seven genes (rpoB, katG, inhA, eis, rrs, gyrA, and gyrB) involved in resistance to five drugs were amplified in each sample using a 20 μL PCR system. Rifampicin resistance was diagnosed based on rpoB sequence. Isoniazid resistance was diagnosed based on katG and inhA sequences. Fluoroquinolone (FQ) resistance was diagnosed based on gyrA and gyrB sequences. Levoxloxacin and moxifloxacin were classified as FQ drugs in the current study. Resistance related to amikacin was based on rrs sequence. Capreomycin resistance was diagnosed using rrs and eis sequences. For analysis of the 20 study samples, we usually prepared the PCR mix at a reaction volume for 21 samples according to the Platinum II Taq Hot-Start DNA Polymerase PCR kit (Thermo Fisher). The reaction mix included 8.4 μL each of 10 μM forward and reverse primers(see Table 1), 84 μL 5X Platinum II PCR Buffer, 8.4 μL 25 mM dNTP MIX, 84 µL Platinum GC Enhancer, 80.64 µL ddH2O, 3.36 μL Platinum II Taq Hot-Start DNA Polymerase, and 2 μL gDNA from each sample. The PCR cycling profile consisted of 94°C for 2 min, followed by 35 cycles of 94°C for 15 s, 60°C for 15 s, and 68°C for 15 s, and completing with a 4°C hold. PCR products were quantified at approximately 100–200 fmol for each sample using a gradient dilution. Equimolar amounts of the PCR products were mixed and a portion sent to Sangon Biotech (Shanghai) for Sanger sequencing. The Sanger sequencing analysis usually required several days due to transportation and processing times.

2.3. Nanopore Library Preparation and Sequencing

Multiplex PCR amplicons of the 20 study samples were prepared using a Ligation Sequencing Kit (SQK-LSK109; ONT, Oxford, England) and Native Barcoding Kit (EXP-NBD104 and EXP-NBD114; ONT). End-prep and native barcode ligation were performed for approximately 3 h using a 100–200 fmol sample diluted in 65 μL nuclease-free water according to the Native Barcoding Kit amplicon protocol. An adapter ligation and cleaning step was performed using NEB ligation and Agencourt AMPure XP beads (Beckman Coulter, USA) with the final adapter-ligated DNA library being 50–100 fmol. The library was loaded into a R9.4 flow cell (ONT) with 851 effective pores and then sequenced using a GridION instrument (ONT). After the sequencing run was completed, the flow cell was cleaned using a Flow Cell Wash Kit (EXP-WSH002; ONT) according to the manufacturer’s protocol and stored at 4°C until any subsequent use.

2.4. Nanopore Data Analysis

Nanopore raw data (fast5) were analyzed using Guppy Version 4.5.2 software (ONT) with a q-score threshold of 9. The data were re-basecalled using the parameter “--config dna_r9.4.1_450 bps_hac.cfg--num_callers 4 --cpu_threads_per_caller 4”. The barcode was recognized using the parameter “--barcode_kits “EXP-NBD104 EXP-NBD114” and trimmed using the parameter “--config configuration.cfg--trim_barcodes”. The sequencing data were counted using NanoPlot v1.28.1 [16] and variant calls found using medaka v1.3.2 [17] (-m r941_min_high_g360). The raw reads were mapped to seven gene-region combinations (3,675 bp) and the trimmed reads were assembled into reference genomes using Genomics software (version 3.0; Hangzhou Baiyi Technology Co., Ltd.). Sequence depth of the ONT sequencing reads versus the seven gene-region combinations were then assessed using SAMtools [18], Minimap2 [19], and bamdst v1.0.9 (https://github.com/shiquan/bamdst).

2.5. Consensus Generation for Gene Variation Identification

Accuracy of the nanopore sequencing variants was determined by aligning the assembled nanopore sequences of the seven gene regions (rpoB, katG, inhA, eis, rrs, gyrA, and gyrB) to that of the Sanger reference sequences using ClustalW (https://www.ebi.ac.uk/Tools/msa/clustalo/). Percent identities were determined for each alignment to ascertain the accuracy of the nanopore sequencing. Significantly, because of the low quality, the Sanger sequences processed by trimming 40 bp from the head and 25 bp bases from the tail, including the primer binding sites [20–22].

Antimicrobial susceptibility testing and minimum inhibitory concentrations (MICs).

Drug susceptibility testing was performed using the assay from the Chinese Antituberculosis Association [23] and resistance ratios were determined [24]. To evaluate the accuracy of the biomolecular technology, we considered the consistency for each drug by comparing the sequencing and MIC results.

3. Results

3.1. Nanopore Sequencing Results

To overcome the time-consuming and tedious process of sample preparation for DNA sequencing, we attempted to amplify key genes directly from routine clinical specimens without a DNA purification step. The time-to-result analysis of the 20 study specimens took approximately 12 h, including gDNA extraction (3 h), MPCR amplification (1 h), library preparation (4 h), nanopore sequencing (3 h), and data analysis(1 h).

The flow cell had an enormous excess capacity for PCR amplicons in our size range, and the 20 sample amplicons were sequenced on a flow cell with approximately 100 activated pores. A total of 1.13 M reads (606.43 Mb) were generated in 3 h, averaging 18.3 Mb per sample. The quality of the trim barcode sequencing data from multiplex ONT sequencing experiments was analyzed using NanoPlot (see Supplementary Material 2). No major differences were noted when evaluating each sample sequencing output, except for samples Y183 and Y83 (see Supplementary Material 2), which may have been lost during library preparation. We found that the depth of coverage showed the same trend among the different samples. The depth of coverage findings for katG (mean depth 10,097) and gyrA (mean depth 10,331) suggested a large amount of data for these genes. Some biases were observed in the reads. For instance, gyrB (mean depth 350) and eis (mean depth 350) were detected at lower abundances compared to that of the other genes (see Figure 1). The mean read length varied much less (479.5–497.2 bp) than that of the depth of coverage. Meanwhile, sequencing quality was shown to be consistent among the 20 samples, with high mean base quality scores ranging from 13.0 to 13.1. We focused on variants in rpoB, katG, inhA, gyrA, gyrB, rrs, and eis, and calculated the depth at specific positions for these genes. Among all variants, the minimum depth was 162 × in gyrB (G1510 A) of sample Y105, and the maximum depth was 13,901 × in gyrA (GAC-94-GGC) of sample Y88 (see Supplementary Material 3). We confirmed 100% nucleotide identity for the ∼500 bp amplicons of the seven genes compared to that of the Sanger sequencing results (see Table 2).

3.2. Resistance Gene Identification

A total of 17 single nucleotide variants contained six unknown mutations, (gyrA : G61 C; eis: C257 T; gyrB : G1510 A, G1255 A; rpoB : A1291 G, and A1379 C; see Supplementary Material 3). The effectiveness of using the samples and nanopore sequencing for the identification of resistance genes was evaluated by comparison with the Sanger sequencing results. The concordance of variant calls between Sanger sequencing (see Supplementary Material 4) and nanopore sequencing was 100%. For example, the results revealed a consistent base call of “A” for the majority of the reads at position 1510 in gyrA of Y105 (see Figure 2). Seven variant types in rpoB were detected among the 20 samples: TCG-531-TTG (6/13), GAC-516-GGC (2/13), CAC-526-TGC (1/13), ATC-572-CTC (1/13), TCG-531-TTT (1/13), A1291 G (1/13), and A1379 C (1/13). A coding missense mutation at position 315 was the only change in katG, which was present in samples Y12, Y80, Y105, Y252, and Y256. Three gyrB variants were found in samples Y76, Y105, and Y252. Substitutions G61 C and G284 C in gyrA were detected in all the samples.

3.3. Consistency between MIC and Molecular Sequencing

Comparing the nanopore and Sanger sequencing results with that of the MIC results (see Supplementary Material 3 and Supplementary Material 5) revealed 100% agreement for FQ and 80% agreement for rifampicin. Unfortunately, there was only a 30%, 25%, and 20% agreement for amikacin, isoniazid, and capreomycin, respectively.

4. Discussion

This was a valuable study aimed at elucidating a method for the rapid identification of drug-resistant TB. In addition to the function capability of nanopore sequencing to easily and quickly generate drug resistance regions, sequencing can be performed in the laboratory using its removable ability. Sample preparation was simple and involved end repair, barcoding, and adapter ligation, all of which could be performed in a tube.

A comparison between the Sanger and ONT sequencing results for the study samples revealed 100% identity between the methods, indicating a potential for using ONT sequencing for drug susceptibility prediction. Consensus sequences obtained via ONT sequencing were longer than those obtained via Sanger sequencing (see Table 2), as the latter typically lacked a number of bases at the 5′ and 3′ ends due to lower base quality. However, this had minimal impact in that we were only concerned about the target regions, which were not affected by the truncated sequence data.

Currently, identification of clinically relevant drug resistance relies mainly on laboratory culture techniques. However, culture-dependent methods are time-consuming (average 11.5 d) [25] and tend to produce poorly reproducible test results because of the MIC of some drugs [26]. In our current study, samples (n = 20) were simultaneously multiplexed in one run required approximately 12 h. Moreover, flexibility of the study method design makes this strategy a variable sequencing purpose for genotyping. This speed of the ONT sequencing for identifying drug-resistant TB will decrease the hold time by weeks, which is critical for situations in an emergency phase.

Next-generation sequencing (NGS) is an alternative approach for detecting drug resistance-associated variants in clinical specimens [27–29]. While WGS is widely used to explore comprehensive genomic information [30, 31], in many cases, there are inadequate amounts of clinical specimens for isolating DNA or sequencing. Compared to WGS, amplicon sequencing, which employs PCR products to detect gene information, is an effective and accurate approach for determining genes of known pathogens and has been used for detecting Zika virus [32], polioviruses [33], and enteroviruses [34]. Meanwhile, although high-throughput sequencing offers a faster turnaround time than other detection techniques, purification of bacterial DNA from samples, which typically takes several hours, is a rate-limiting step in the workflow. To further improve the sample preparation process, we attempted to amplify genes directly from clinical specimens without the use of DNA purification. Moreover, instead of sequencing the entire genome (approximately 4.4 Mb), we only focused on recognized target areas to be sequenced (approximately 4 kb).

We analyzed 20 sputum specimens using multiple approaches and demonstrated consistent SNP results. Surprisingly, during this study, we identified a new variant in gyrA, G61 C, which may have been the reason for the observed resistance in the MIC test (Supplementary Material 5). Further evaluation of this variant is our next priority. We have also learned that G284 C (Ser95Thr) in gyrA does not lead to FQ resistance in M. tuberculosis [35]. However, variants in gyrB that were only detected in Y76 (GGG-551-AGG), Y105 (G1510 A), and Y252 (G1255 A) may be the novel causes for drug resistance, which differs from previous studies [36, 37]. The common variant TCG-531-TTG in rpoB and AGC-315-ACC in katG were revealed as the primary reasons for MTB resistance to rifampicin and isoniazid, which is consistent with the results of Solo et al. [38], Feizabadi et al. [39], and Sun et al. [40]. In terms of the 15 amikacin- and capreomycin-resistant MIC results, Y80 was the only sample that may have been had its drug resistance caused by the A1401 G variant, which conflicts with results from previous studies [41, 42]. Most of the inconformity with the MIC results may have been caused by the limited area of amplification and/or the presence of other unknown resistance genes. For example, variations in fabG1 or other katG mutations contribute to isoniazid resistance discordance [43]. Nanopore sequencing demonstrated lower consistency for the detection of amikacin and capreomycin resistance. One reason may have been it missed variants in additional genes, such as tlyA and gidB [44]. This serves as an important reminder that the mechanisms of resistance caused by variants in other genes we may be missed. In the future, it may be necessary to redesign the PCR to amplify gene region. Moreover, we suspect that insufficient outputs using this diagnostic could lead to false-negatives [45]. In view of the above-mentioned fact we just increased the output in another half an hour or one hour.

Although culture identification of M. tuberculosis is the gold standard for TB diagnosis, we need further argument like Illumina MiSeq sequencing because of its ability in false-positive results [46]. In addition, laboratory cross-contamination [47] should be suspected when MIC findings are inconsistent with molecular detection technology results. Moreover, sample gene expression [48] and drug efflux mechanisms [49] should not be ignored when formulating strategies for combatting drug resistance. Of these variants, the (C257 T) caused (Ala to Val) in eis refer to CPM detected by these two sequencing platforms seems to be a meaningless change in AMR.

5. Conclusion

In this study, we developed a nanopore-based panel for the molecular diagnosis of TB through the direct amplification of resistance genes in 20 clinical specimens in a single reaction. This approach provides an attractive option for the detection of antimicrobial susceptibility. Further improvements and the establishment of a relatively simple workflow for predicting drug resistance via nanopore sequencing would reduce the turnaround time of sample analysis and provide a viable method applicable to clinical settings. However, this technology has limited ability to detect variants in other unknown gene regions. That noted, this is a flexible diagnostic platform, and new panels can be added according to the specific demands by designing the appropriate primers.

Data Availability

The data in our study have been deposited in The National Center for Biotechnology Information (NCBI) under the accession number PRJNA766801 (https://dataview.ncbi.nlm.nih.gov/object/PRJNA766801?reviewer=tlpsmq4fkj689otofb67sho4lq).

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

Kaishun Zhao and Chunlin Tu are equal contributions.

Acknowledgments

This work was supported by grants from the Scientific Research Project of Shanghai Municipal Health Commission (201940315), Shanghai Jiading District Health System Scientific Research Project (2017ZD05), and Scientific Research Project of Jiading District Science and Technology Commission, Shanghai, China (JDKW-2019-W03).

Supplementary Materials

Supplementary Material 1: Details of the 20 Mycobacterium tuberculosis specimens. Supplementary Material 2: Summary of quality statistics of multiplexed trim sequencing data. Supplementary Material 3: Targeted mutations identified by nanopore sequencing of 20 Mycobacterium tuberculosis specimens. Supplementary Material 4: Sanger sequencing data for 20 Mycobacterium tuberculosis specimens. Supplementary Material 5: MIC diagnostic performance of 20 Mycobacterium tuberculosis specimens. (Supplementary Materials)

References

World Health Organization, Global Tuberculosis Report 2020 [R], World Health Organization, Geneva, Switzerland, 2016.
S. Montalvo-Quirós, S. Gómez-Graa, M. Vallet-Regí, R. C. Prados-Rosales, B. González, and J. L. Luque-Garcia, “Mesoporous silica nanoparticles containing silver as novel antimycobacterial agents against Mycobacterium tuberculosis,” Colloids and Surfaces B: Biointerfaces, vol. 197, Article ID 111405, 2021.
View at: Google Scholar
Y. Zhao, S. Xu, L. Wang et al., “National survey of drug-resistant tuberculosis in China,” New England Journal of Medicine, vol. 366, no. 23, pp. 2161–2170, 2012.
View at: Publisher Site | Google Scholar
J. Cervantes, N. Yokobori, and B.-Y. Hong, “Genetic identification and drug-resistance characterization of Mycobacterium tuberculosis using a portable sequencing device. a pilot study,” Antibiotics, vol. 9, no. 9, p. 548, 2020.
View at: Publisher Site | Google Scholar
A. Koch, H. Cox, and V. Mizrahi, “Drug-resistant tuberculosis: challenges and opportunities for diagnosis and treatment,” Current Opinion in Pharmacology, vol. 42, pp. 7–15, 2018.
View at: Publisher Site | Google Scholar
A. Engström, N. Morcillo, B. Imperiale, S. E. Hoffner, and P. Juréen, “Detection of first- and second-line drug resistance in Mycobacterium tuberculosis clinical isolates by pyrosequencing,” Journal of Clinical Microbiology, vol. 50, no. 6, pp. 2026–2033, 2012.
View at: Publisher Site | Google Scholar
P. J. Campbell, G. P. Morlock, R. D. Sikes et al., “Molecular detection of mutations associated with first-and second-line drug resistance compared with conventional drug susceptibility testing of Mycobacterium tuberculosis,” Antimicrobial Agents and Chemotherapy, vol. 55, no. 5, pp. 2032–2041, 2011.
View at: Publisher Site | Google Scholar
F. W. Gregory, “Clinical application of whole-genome sequencing,” JAMA: Journal of the American Medical Association, vol. 311, no. 10, pp. 1017–1019, 2014.
View at: Google Scholar
T. Magdy, H. H. Kuo, and P. W. Burridge, “Precise and cost-effective nanopore sequencing for post-GWAS fine-mapping and causal variant identification,” iScience, vol. 23, no. 4, Article ID 100971, 2020.
View at: Publisher Site | Google Scholar
L. R. Runtuwene, J. S. B. Tuda, A. E. Mongan, and Y. Suzuki, “On-site MinION sequencing,” Single Molecule and Single Cell Sequencing, vol. 1129, pp. 143–150, 2019.
View at: Publisher Site | Google Scholar
P. M. Ashton, S. Nair, T. Dallman et al., “MinION nanopore sequencing identifies the position and structure of a bacterial antibiotic resistance island,” Nature Biotechnology, vol. 33, no. 3, pp. 296–300, 2015.
View at: Publisher Site | Google Scholar
D. Golparian, V. Donà, L. Sánchez-Busó et al., “Antimicrobial resistance prediction and phylogenetic analysis of Neisseria gonorrhoeae isolates using the Oxford Nanopore MinION sequencer,” Scientific Reports, vol. 8, Article ID 17596, 2018.
View at: Publisher Site | Google Scholar
M. E. Pitt, S. H. Nguyen, T. P. S. Duarte et al., “Evaluating the genome and resistome of extensively drug-resistant Klebsiella pneumoniae using native DNA and RNA Nanopore sequencing,” GigaScience, vol. 9, no. 2, 2020.
View at: Publisher Site | Google Scholar
S. George, Y. Xu, G. Rodger et al., “DNA thermo-protection facilitates whole genome sequencing of mycobacteria direct from clinical samples by the nanopore platform,” Journal of Clinical Microbiology, vol. 58, no. 10, 2020.
View at: Publisher Site | Google Scholar
A. A. Votintseva, P. Bradley, L. Pankhurst et al., “Same-day diagnostic and surveillance data for tuberculosis via whole-genome sequencing of direct respiratory samples,” Journal of Clinical Microbiology, vol. 55, no. 5, pp. 1285–1298, 2017.
View at: Publisher Site | Google Scholar
D. C. Wouter, D. Svenn, D. T. Schultz, M. Cruts, and C. van Broeckhoven, “NanoPack: visualizing and processing long read sequencing data,” Bioinformatics, vol. 34, no. 15, pp. 2666–2669, 2108.
View at: Google Scholar
Geneva Switzerland Who Global Tuberculosis Programme, “Sequence correction provided by ONT Research,” 2021, https://github.com/nanoporetech/medaka.
View at: Google Scholar
H. Li, B. Handsaker, A. Wysoker et al., “The sequence alignment/map format and SAMtools,” Bioinformatics, vol. 25, no. 16, pp. 2078-2079, 2009.
View at: Publisher Site | Google Scholar
H. Li, “Minimap2: fast pairwise alignment for long nucleotide sequences,” Bioinformatics, vol. 34, no. 18, pp. 3094–3100, 2017.
View at: Google Scholar
J. T. T. Hill, B. L. Demarest, B. W. Bisgrove, Y. C. Su, M. Smith, and H. J. Yost, “Poly peak parser: method and software for identification of unknown indels using sanger sequencing of polymerase chain reaction products,” Developmental Dynamics, vol. 243, no. 12, pp. 1632–1636, 2014.
View at: Publisher Site | Google Scholar
L. Chen, Y. Cai, G. Zhou et al., “Rapid sanger sequencing of the 16S rRNA gene for identification of some common pathogens,” PLoS One, vol. 9, no. 2, Article ID e88886, 2014.
View at: Publisher Site | Google Scholar
D. Tewari, S. Cieply, and J. Livengood, “Identification of bacteria recovered from animals using the 16S ribosomal RNA gene with pyrosequencing and sanger sequencing,” Journal of Veterinary Diagnostic Investigation, vol. 23, no. 6, pp. 1104–1108, 2011.
View at: Publisher Site | Google Scholar
China Antituberculosis Association, “The laboratory science procedure of diagnostic bacteriology in tuberculosis,” Bulletin of Chinese Antituberculosis Association, vol. 18, 1996.
View at: Google Scholar
World Health Organization (WHO), “Global tuberculosis programme; international union against tuberculosis and lung disease,” Guidelines for Surveillance of Drug Resistance in Tuberculosis, World Health Organization, Geneva, Switzerland, 1997.
View at: Google Scholar
J. P. Libonati, C. E. Stager, J. R. Davis, and S. H. Siddiqi, “Direct antimicrobial drug susceptibility testing of Mycobacterium tuberculosis by the radiometric method,” Diagnostic Microbiology and Infectious Disease, vol. 10, no. 1, pp. 41–48, 1988.
View at: Publisher Site | Google Scholar
S. J. Kim, “Drug-susceptibility testing in tuberculosis: methods and reliability of results,” European Respiratory Journal, vol. 25, no. 3, pp. 564–569, 2005.
View at: Publisher Site | Google Scholar
H. Zhang, D. Li, L. Zhao et al., “Genome sequencing of 161 Mycobacterium tuberculosis isolates from China identifies genes and intergenic regions associated with drug resistance,” Nature Genetics, vol. 45, no. 10, pp. 1255–1260, 2013.
View at: Publisher Site | Google Scholar
I. Basu, J. E. Bower, S. A. Roberts et al., “Utility of whole genome sequencing for multidrug resistant Mycobacterium tuberculosis isolates in a reference TB laboratory in New Zealand,” New Zealand Medical Journal, vol. 131, no. 1487, pp. 15–22, 2018.
View at: Google Scholar
S. M. Regmi, A. Chaiprasert, S. Kulawonganunchai et al., “Whole genome sequence analysis of multidrug-resistant Mycobacterium tuberculosis Beijing isolates from an outbreak in Thailand,” Molecular Genetics and Genomics, vol. 290, no. 5, pp. 1933–1941, 2015.
View at: Publisher Site | Google Scholar
J. E. Phelan, D. M. O’Sullivan, D. Machado et al., “Integrating informatics tools and portable sequencing technology for rapid detection of resistance to anti-tuberculous drugs,” Genome Medicine, vol. 11, no. 1, p. 41, 2019.
View at: Publisher Site | Google Scholar
M. Ezewudo, A. Borens, Á. Chiner-Oms et al., “Integrating standardized whole genome sequence analysis with a global Mycobacterium tuberculosis antibiotic resistance knowledgebase,” Scientific Reports, vol. 8, no. 1, 2018.
View at: Publisher Site | Google Scholar
J. Quick, N. D. Grubaugh, S. T. Pullan et al., “Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples,” Nature Protocols, vol. 12, no. 6, pp. 1261–1276, 2017.
View at: Publisher Site | Google Scholar
A. G. Shaw, M. Majumdar, C. Troman et al., “Rapid and sensitive direct detection and identification of poliovirus from stool and environmental surveillance samples using nanopore sequencing,” Journal of Clinical Microbiology, vol. 58, no. 9, 2020.
View at: Publisher Site | Google Scholar
C. Grdel, M. Miani, C. Baumann et al., “Whole-genome sequencing of human enteroviruses from clinical samples by nanopore direct RNA sequencing,” Viruses, vol. 12, no. 8, 2020.
View at: Google Scholar
E. Y. Nosova, A. A. Bukatina, Y. D. Isaeva, M. V. Makarova, K. Y. Galkina, and A. M. Moroz, “Analysis of mutations in the gyrA and gyrB genes and their association with the resistance of Mycobacterium tuberculosis to levofloxacin, moxifloxacin and gatifloxacin,” Journal of Medical Microbiology, vol. 62, no. 1, pp. 108–113, 2013.
View at: Publisher Site | Google Scholar
S. Feuerriegel, H. S. Cox, N. Zarkua et al., “Sequence analyses of just four genes to detect extensively drug-resistant Mycobacterium tuberculosis strains in multidrug-resistant tuberculosis patients undergoing treatment,” Antimicrobial Agents and Chemotherapy, vol. 53, no. 8, pp. 3353–3356, 2009.
View at: Publisher Site | Google Scholar
A. Aubry, N. Veziris, E. Cambau, C. Truffot-Pernot, V. Jarlier, and L. M. Fisher, “Novel GyrAse mutations in quinolone-resistant and-hypersusceptible clinical isolates of Mycobacterium tuberculosis: functional analysis of mutant enzymes,” Antimicrobial Agents & Chemotherapy, vol. 50, no. 1, pp. 104–112, 2005.
View at: Google Scholar
E. S. Solo, C. Nakajima, T. Kaile et al., “Mutations in rpoB and katG genes and the inhA operon in multidrug-resistant Mycobacterium tuberculosis isolates from Zambia,” Journal of Global Antimicrobial Resistance, vol. 22, pp. 302–307, 2020.
View at: Publisher Site | Google Scholar
M. Feizabadi, M. Haeili, A. Fooladi et al., “Rapid screening of rpoB and katG mutations in Mycobacterium tuberculosis isolates by high-resolution melting curve analysis,” Indian Journal of Medical Microbiology, vol. 32, no. 4, pp. 398–403, 2014.
View at: Google Scholar
Y. Sun, S. Li, L. Zhou et al., “A rapid fluorescence polarization-based method for genotypic detection of drug resistance in Mycobacterium tuberculosis,” Applied Microbiology and Biotechnology, vol. 98, no. 9, pp. 4095–4105, 2014.
View at: Publisher Site | Google Scholar
Q. Du, G. Dai, Q. Long et al., “Mycobacterium tuberculosis rrs A1401G mutation correlates with high-level resistance to kanamycin, amikacin, and capreomycin in clinical isolates from mainland China,” Diagnostic Microbiology and Infectious Disease, vol. 77, no. 2, pp. 138–142, 2013.
View at: Publisher Site | Google Scholar
F. A. Sirgel, M. Tait, R. M. Warren et al., “Mutations in therrsA1401G gene and phenotypic resistance to amikacin and capreomycin inMycobacterium tuberculosis,” Microbial Drug Resistance, vol. 18, no. 2, pp. 193–197, 2012.
View at: Publisher Site | Google Scholar
J. N. Torres, L. V. Paul, T. C. Rodwell et al., “Novel katG mutations causing isoniazid resistance in clinical M. tuberculosis isolates,” Emerging Microbes & Infections, vol. 4, no. 7, p. e42, 2015.
View at: Publisher Site | Google Scholar
S. B. Georghiou, M. Magana, R. S. Garfein, D. G. Catanzaro, A. Catanzaro, and T. C. Rodwell, “Evaluation of genetic mutations associated with Mycobacterium tuberculosis resistance to amikacin, kanamycin and capreomycin: a systematic review,” PLoS One, vol. 7, no. 3, Article ID e33275, 2012.
View at: Publisher Site | Google Scholar
C. Smith, T. Halse, J. Shea et al., “Assessing Nanopore sequencing for clinical diagnostics: a comparison of NGS methods for Mycobacterium tuberculosis,” Journal of Clinical Microbiology, vol. 59, no. 1, 2020.
View at: Publisher Site | Google Scholar
J. Min, K. Kim, H. Choi et al., “Investigation of false-positive Mycobacterium tuberculosis culture tests using whole genome sequencing,” Annals of Thoracic Medicine, vol. 14, no. 1, pp. 90–93, 2019.
View at: Publisher Site | Google Scholar
L. Fitzpatrick, C. Braden, W. Cronin et al., “Investigation of laboratory cross-contamination of Mycobacterium tuberculosis cultures,” Clinical Infectious Diseases : An Official Publication of the Infectious Diseases Society of America, vol. 38, no. 6, pp. e52–e54, 2004.
View at: Publisher Site | Google Scholar
P. Freihofer, R. Akbergenov, Y. Teo, R. Juskeviciene, D. I. Andersson, and E. C. Böttger, “Nonmutational compensation of the fitness cost of antibiotic resistance in mycobacteria by overexpression of tlyA rRNA methylase,” RNA, vol. 22, no. 12, pp. 1836–1843, 2016.
View at: Publisher Site | Google Scholar
K. Poole, “Efflux-mediated antimicrobial resistance,” Journal of Antimicrobial Chemotherapy, vol. 56, no. 1, pp. 20–51, 2005.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Kaishun Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

2147

Downloads

1176

Citations