MTBseq: a comprehensive pipeline for whole genome sequence analysis of Mycobacterium tuberculosis complex isolates

Analyzing whole-genome sequencing data of Mycobacterium tuberculosis complex (MTBC) isolates in a standardized workflow enables both comprehensive antibiotic resistance profiling and outbreak surveillance with highest resolution up to the identification of recent transmission chains. Here, we present MTBseq, a bioinformatics pipeline for next-generation genome sequence data analysis of MTBC isolates. Employing a reference mapping based workflow, MTBseq reports detected variant positions annotated with known association to antibiotic resistance and performs a lineage classification based on phylogenetic single nucleotide polymorphisms (SNPs). When comparing multiple datasets, MTBseq provides a joint list of variants and a FASTA alignment of SNP positions for use in phylogenomic analysis, and identifies groups of related isolates. The pipeline is customizable, expandable and can be used on a desktop computer or laptop without any internet connection, ensuring mobile usage and data security. MTBseq and accompanying documentation is available from https://github.com/ngs-fzb/MTBseq_source

What Is Resistance? Impact of Phenotypic versus Molecular Drug Resistance Testing on Therapy for Multi- and Extensively Drug-Resistant Tuberculosis.

Rapid and accurate drug susceptibility testing (DST) is essential for the treatment of multi- and extensively drug-resistant tuberculosis (M/XDR-TB). We compared the utility of genotypic DST assays with phenotypic DST (pDST) using Bactec 960 MGIT or Löwenstein-Jensen to construct M/XDR-TB treatment regimens for a cohort of 25 consecutive M/XDR-TB patients and 15 possible anti-TB drugs. Genotypic DST results from Cepheid GeneXpert MTB/RIF (Xpert) and line probe assays (LPAs; Hain GenoType MTBDRplus 2.0 and MTBDRsl 2.0) and whole-genome sequencing (WGS) were translated into individual algorithm-derived treatment regimens for each patient. We further analyzed if discrepancies between the various methods were due to flaws in the genotypic or phenotypic test using MIC results. Compared with pDST, the average agreement in the number of drugs prescribed in genotypic regimens ranged from just 49% (95% confidence interval [CI], 39 to 59%) for Xpert and 63% (95% CI, 56 to 70%) for LPAs to 93% (95% CI, 88 to 98%) for WGS. Only the WGS regimens did not contain any drugs to which pDST showed resistance. Importantly, MIC testing revealed that pDST likely underestimated the true rate of resistance for key drugs (rifampin, levofloxacin, moxifloxacin, and kanamycin) because critical concentrations (CCs) were too high. WGS can be used to rule in resistance even in M/XDR strains with complex resistance patterns, but pDST for some drugs is still needed to confirm susceptibility and construct the final regimens. Some CCs for pDST need to be reexamined to avoid systematic false-susceptible results in low-level resistant isolates

What Is Resistance? Impact of Phenotypic versus Molecular Drug Resistance Testing on Therapy for Multi- and Extensively Drug-Resistant Tuberculosis.

Rapid and accurate drug susceptibility testing (DST) is essential for the treatment of multi- and extensively drug-resistant tuberculosis (M/XDR-TB). We compared the utility of genotypic DST assays with phenotypic DST (pDST) using Bactec 960 MGIT or Löwenstein-Jensen to construct M/XDR-TB treatment regimens for a cohort of 25 consecutive M/XDR-TB patients and 15 possible anti-TB drugs. Genotypic DST results from Cepheid GeneXpert MTB/RIF (Xpert) and line probe assays (LPAs; Hain GenoType MTBDRplus 2.0 and MTBDRsl 2.0) and whole-genome sequencing (WGS) were translated into individual algorithm-derived treatment regimens for each patient. We further analyzed if discrepancies between the various methods were due to flaws in the genotypic or phenotypic test using MIC results. Compared with pDST, the average agreement in the number of drugs prescribed in genotypic regimens ranged from just 49% (95% confidence interval [CI], 39 to 59%) for Xpert and 63% (95% CI, 56 to 70%) for LPAs to 93% (95% CI, 88 to 98%) for WGS. Only the WGS regimens did not contain any drugs to which pDST showed resistance. Importantly, MIC testing revealed that pDST likely underestimated the true rate of resistance for key drugs (rifampin, levofloxacin, moxifloxacin, and kanamycin) because critical concentrations (CCs) were too high. WGS can be used to rule in resistance even in M/XDR strains with complex resistance patterns, but pDST for some drugs is still needed to confirm susceptibility and construct the final regimens. Some CCs for pDST need to be reexamined to avoid systematic false-susceptible results in low-level resistant isolates

Reply to Dookie et al., "Whole-Genome Sequencing To Guide the Selection of Treatment for Drug-Resistant Tuberculosis".

We would like to thank Dr. Navisha Dookie and colleagues for their comment and for sharing their view on genotypic drug-susceptibility testing (DST) to predict phenotypic drug-susceptibility for the management of patients with drug-resistant tuberculosis (1, 2). We agree that DST should ideally be conducted for all drugs in a regimen, particularly given the severe adverse effects of second-line anti-tuberculosis drugs. Yet, even whole-genome sequencing in not a one-size-fits-all solution. Instead, the goal has to be to optimally combine genotypic and phenotypic assays in a diagnostic algorithm that capitalises on the strengths of both approaches whilst including appropriate reflex and confirmatory testing to minimize the impact of their respective limitations, including systematic errors (3)

Reply to Dookie et al., "Whole-Genome Sequencing To Guide the Selection of Treatment for Drug-Resistant Tuberculosis".

We would like to thank Dr. Navisha Dookie and colleagues for their comment and for sharing their view on genotypic drug-susceptibility testing (DST) to predict phenotypic drug-susceptibility for the management of patients with drug-resistant tuberculosis (1, 2). We agree that DST should ideally be conducted for all drugs in a regimen, particularly given the severe adverse effects of second-line anti-tuberculosis drugs. Yet, even whole-genome sequencing in not a one-size-fits-all solution. Instead, the goal has to be to optimally combine genotypic and phenotypic assays in a diagnostic algorithm that capitalises on the strengths of both approaches whilst including appropriate reflex and confirmatory testing to minimize the impact of their respective limitations, including systematic errors (3)