A pragmatic pipeline for drug resistance identification in Mycobacterium tuberculosis using whole genome sequencing

BackgroundDelays in accurate diagnosis of drug resistant tuberculosis (DR-TB) can hinder treatment. Whole genome sequencing (WGS) provides more information than standard molecular and phenotypic testing, but commonly used platforms are expensive to implement, and data interpretation requires significant expertise. AimsWe aimed to optimise a TB WGS diagnostic pipeline balancing user-friendliness, cost- effectiveness and time to results, whilst ensuring accuracy. Materials and methodsGrowth conditions, DNA extraction protocols and Oxford Nanopore Technologies (ONT) library preparation kits were compared. Software for basecalling and analysis were evaluated to find the most accurate resistance SNP and lineage predictor. ResultsOptimally, a spin-column CTAB DNA extraction method was combined with the RBK110.96 library preparation kit, high accuracy basecalling and data analysis using TB-Profiler. Compared with Illumina, the pipeline was concordant for 16/17 (94%) isolates (lineage) and for 17/17 (100%) isolates (resistance SNPs). Our pipeline was 71% (12/17) concordant with phenotypic drug susceptibility test (DST) results. Time-to-diagnosis was around four weeks. ConclusionsThis optimised TB sequencing pipeline requires less time expertise to run and analyse than Illumina, takes less time than phenotypic DSTs and the results are comparable with Illumina. The cost per sample is comparable with other methods. These features make it an important tool for incorporating into routine DR-TB diagnostic pipelines and larger scale drug resistance surveillance in all settings.