Prediction of Susceptibility to First-Line Tuberculosis Drugs by DNA Sequencing.

BACKGROUND: The World Health Organization recommends drug-susceptibility testing of Mycobacterium tuberculosis complex for all patients with tuberculosis to guide treatment decisions and improve outcomes. Whether DNA sequencing can be used to accurately predict profiles of susceptibility to first-line antituberculosis drugs has not been clear. METHODS: We obtained whole-genome sequences and associated phenotypes of resistance or susceptibility to the first-line antituberculosis drugs isoniazid, rifampin, ethambutol, and pyrazinamide for isolates from 16 countries across six continents. For each isolate, mutations associated with drug resistance and drug susceptibility were identified across nine genes, and individual phenotypes were predicted unless mutations of unknown association were also present. To identify how whole-genome sequencing might direct first-line drug therapy, complete susceptibility profiles were predicted. These profiles were predicted to be susceptible to all four drugs (i.e., pansusceptible) if they were predicted to be susceptible to isoniazid and to the other drugs or if they contained mutations of unknown association in genes that affect susceptibility to the other drugs. We simulated the way in which the negative predictive value changed with the prevalence of drug resistance. RESULTS: A total of 10,209 isolates were analyzed. The largest proportion of phenotypes was predicted for rifampin (9660 [95.4%] of 10,130) and the smallest was predicted for ethambutol (8794 [89.8%] of 9794). Resistance to isoniazid, rifampin, ethambutol, and pyrazinamide was correctly predicted with 97.1%, 97.5%, 94.6%, and 91.3% sensitivity, respectively, and susceptibility to these drugs was correctly predicted with 99.0%, 98.8%, 93.6%, and 96.8% specificity. Of the 7516 isolates with complete phenotypic drug-susceptibility profiles, 5865 (78.0%) had complete genotypic predictions, among which 5250 profiles (89.5%) were correctly predicted. Among the 4037 phenotypic profiles that were predicted to be pansusceptible, 3952 (97.9%) were correctly predicted. CONCLUSIONS: Genotypic predictions of the susceptibility of M. tuberculosis to first-line drugs were found to be correlated with phenotypic susceptibility to these drugs. (Funded by the Bill and Melinda Gates Foundation and others.).

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.

Rapid Microarray-Based Detection of Rifampin, Isoniazid, and Fluoroquinolone Resistance in Mycobacterium tuberculosis by Use of a Single Cartridge.

The rapid and robust identification of mutations in Mycobacterium tuberculosis complex (MTBC) strains mediating multidrug-resistant (MDR) and extensively drug-resistant (XDR) phenotypes is crucial to combating the MDR tuberculosis (TB) epidemic. Currently available molecular anti-TB drug susceptibility tests either are restricted to a single target or drug (i.e., the Xpert MTB/RIF test) or present a risk of cross-contamination due to the design limitations of the open platform (i.e., line probe assays). With a good understanding of the technical and commercial boundaries, we designed a test cartridge based on an oligonucleotide array into which dried reagents are introduced and which has the ability to identify MTBC strains resistant to isoniazid, rifampin, and the fluoroquinolones. The melting curve assay interrogates 43 different mutations in the rifampin resistance-determining region (RRDR) of rpoB, rpoB codon 572, katG codon 315, the inhA promoter region, and the quinolone resistance-determining region (QRDR) of gyrA in a closed cartridge system within 90 min. Assay performance was evaluated with 265 clinical MTBC isolates, including MDR/XDR, non-MDR, and fully susceptible isolates, from a drug resistance survey performed in Swaziland in 2009 and 2010. In 99.5% of the cases, the results were consistent with data previously acquired utilizing Sanger sequencing. The assay, which uses a closed cartridge system in combination with a battery-powered Alere q analyzer and which has the potential to extend the current gene target panel, could serve as a rapid and robust point-of-care test in settings lacking a comprehensive molecular laboratory infrastructure to differentiate TB patients infected with MDR and non-MDR strains and to assist clinicians with their early treatment decisions.

New Mycobacterium tuberculosis Complex Sublineage, Brazzaville, Congo.

Tuberculosis is a leading cause of illness and death in Congo. No data are available about the population structure and transmission dynamics of the Mycobacterium tuberculosis complex strains prevalent in this central Africa country. On the basis of single-nucleotide polymorphisms detected by whole-genome sequencing, we phylogenetically characterized 74 MTBC isolates from Brazzaville, the capital of Congo. The diversity of the study population was high; most strains belonged to the Euro-American lineage, which split into Latin American Mediterranean, Uganda I, Uganda II, Haarlem, X type, and a new dominant sublineage named Congo type (n = 26). Thirty strains were grouped in 5 clusters (each within 12 single-nucleotide polymorphisms), from which 23 belonged to the Congo type. High cluster rates and low genomic diversity indicate recent emergence and transmission of the Congo type, a new Euro-American sublineage of MTBC.

Role of Alanine Racemase Mutations in Mycobacterium tuberculosis d-Cycloserine Resistance.

A screening of more than 1,500 drug-resistant strains of Mycobacterium tuberculosis revealed evolutionary patterns characteristic of positive selection for three alanine racemase (Alr) mutations. We investigated these mutations using molecular modeling, in vitro MIC testing, as well as direct measurements of enzymatic activity, which demonstrated that these mutations likely confer resistance to d-cycloserine.

A joint cross-border investigation of a cluster of multidrug-resistant tuberculosis in Austria, Romania and Germany in 2014 using classic, genotyping and whole genome sequencing methods: lessons learnt.

Molecular surveillance of multidrug-resistant tuberculosis (MDR-TB) using 24-loci MIRU-VNTR in the European Union suggests the occurrence of international transmission. In early 2014, Austria detected a molecular MDR-TB cluster of five isolates. Links to Romania and Germany prompted the three countries to investigate possible cross-border MDR-TB transmission jointly. We searched genotyping databases, genotyped additional isolates from Romania, used whole genome sequencing (WGS) to infer putative transmission links, and investigated pairwise epidemiological links and patient mobility. Ten isolates from 10 patients shared the same 24-loci MIRU-VNTR pattern. Within this cluster, WGS defined two subgroups of four patients each. The first comprised an MDR-TB patient from Romania who had sought medical care in Austria and two patients from Austria. The second comprised patients, two of them epidemiologically linked, who lived in three different countries but had the same city of provenance in Romania. Our findings strongly suggested that the two cases in Austrian citizens resulted from a newly introduced MDR-TB strain, followed by domestic transmission. For the other cases, transmission probably occurred in the same city of provenance. To prevent further MDR-TB transmission, we need to ensure universal access to early and adequate therapy and collaborate closely in tuberculosis care beyond administrative borders.

Drug resistance and population structure of M.tuberculosis isolates from prisons and communities in Ethiopia.

BACKGROUND: The population structure and drug resistance pattern of Mycobacterium tuberculosis complex (MTBC) isolates in Ethiopian prisons and some communities is still unknown. METHODS: A comparative cross sectional study was conducted on 126 MTBC strains isolated from prisons and communities in southwestern, southern and eastern Ethiopia. Phenotypic drug susceptibility testing was performed with the MGIT960 system. Combined 24-loci Mycobacterium interspersed repetitive unit-variable number tandem repeat and spacer oligonucleotide typing methods were used to study the MTBC population structure. The obtained data from prisons and communities were compared using statistical tests and regression analysis. RESULTS: A diverse population structure with 11 different lineages and sub-lineages was identified. The predominant strains were the recently described Ethiopia_H37Rv like (27.52%) and Ethiopia_3 (16.51%) with equal lineage distribution between prisons and communities. 28.57% of prison strains and 31.82% of community strains shared the identical genotype with at least one other strain. The multidrug-resistance (MDR) prevalence of the community was 2.27% whereas that of prisons was 9.52%. The highest mono resistance was seen against streptomycin (15.89%). CONCLUSION: Tuberculosis in communities and prisons is caused by a variety of MTBC lineages with predominance of local Ethiopian lineages. The increasing prevalence of MDR MTBC strains is alarming. These findings suggest the need for new approaches for control of MDR tuberculosis in Ethiopia.

PhyResSE: a Web Tool Delineating Mycobacterium tuberculosis Antibiotic Resistance and Lineage from Whole-Genome Sequencing Data.

Antibiotic-resistant tuberculosis poses a global threat, causing the deaths of hundreds of thousands of people annually. While whole-genome sequencing (WGS), with its unprecedented level of detail, promises to play an increasingly important role in diagnosis, data analysis is a daunting challenge. Here, we present a simple-to-use web service (free for academic use at http://phyresse.org). Delineating both lineage and resistance, it provides state-of-the-art methodology to life scientists and physicians untrained in bioinformatics. It combines elaborate data processing and quality control, as befits human diagnostics, with a treasure trove of validated resistance data collected from well-characterized samples in-house and worldwide.

rplC T460C identified as a dominant mutation in linezolid-resistant Mycobacterium tuberculosis strains.

The ribosomal L3 protein was identified as a novel target in linezolid (LZD)-resistant Mycobacterium tuberculosis strains. Next-generation sequencing confirmed rplC T460C as the sole mutation in an LZD-resistant M. tuberculosis H37Rv strain selected in vitro. Sequencing analysis revealed the rplC T460C mutation in eight further LZD-resistant isolates (three in vitro-selected mutants and five patient isolates, including isolates from three different patients that developed LZD resistance during treatment) but in none of the susceptible control strains (n = 84).

Origin and Global Expansion of Mycobacterium tuberculosis Complex Lineage 3.

Mycobacterium tuberculosis complex (MTBC) Lineage 3 (L3) strains are abundant in world regions with the highest tuberculosis burden. To investigate the population structure and the global diversity of this major lineage, we analyzed a dataset comprising 2682 L3 strains from 38 countries over 5 continents, by employing 24-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats genotyping (MIRU-VNTR) and drug susceptibility testing. We further combined whole-genome sequencing (WGS) and phylogeographic analysis for 373 strains representing the global L3 genetic diversity. Ancestral state reconstruction confirmed that the origin of L3 strains is located in Southern Asia and further revealed multiple independent introduction events into North-East and East Africa. This study provides a systematic understanding of the global diversity of L3 strains and reports phylogenetic variations that could inform clinical trials which evaluate the effectivity of new drugs/regimens or vaccine candidates.

Phylogenomics of Mycobacterium africanum reveals a new lineage and a complex evolutionary history.

Human tuberculosis (TB) is caused by members of the Mycobacterium tuberculosis complex (MTBC). The MTBC comprises several human-adapted lineages known as M. tuberculosis sensu stricto, as well as two lineages (L5 and L6) traditionally referred to as Mycobacterium africanum. Strains of L5 and L6 are largely limited to West Africa for reasons unknown, and little is known of their genomic diversity, phylogeography and evolution. Here, we analysed the genomes of 350 L5 and 320 L6 strains, isolated from patients from 21 African countries, plus 5 related genomes that had not been classified into any of the known MTBC lineages. Our population genomic and phylogeographical analyses showed that the unclassified genomes belonged to a new group that we propose to name MTBC lineage 9 (L9). While the most likely ancestral distribution of L9 was predicted to be East Africa, the most likely ancestral distribution for both L5 and L6 was the Eastern part of West Africa. Moreover, we found important differences between L5 and L6 strains with respect to their phylogeographical substructure and genetic diversity. Finally, we could not confirm the previous association of drug-resistance markers with lineage and sublineages. Instead, our results indicate that the association of drug resistance with lineage is most likely driven by sample bias or geography. In conclusion, our study sheds new light onto the genomic diversity and evolutionary history of M. africanum, and highlights the need to consider the particularities of each MTBC lineage for understanding the ecology and epidemiology of TB in Africa and globally.

MDR M. tuberculosis outbreak clone in Eswatini missed by Xpert has elevated bedaquiline resistance dated to the pre-treatment era.

BACKGROUND: Multidrug-resistant (MDR) Mycobacterium tuberculosis complex strains not detected by commercial molecular drug susceptibility testing (mDST) assays due to the RpoB I491F resistance mutation are threatening the control of MDR tuberculosis (MDR-TB) in Eswatini. METHODS: We investigate the evolution and spread of MDR strains in Eswatini with a focus on bedaquiline (BDQ) and clofazimine (CFZ) resistance using whole-genome sequencing in two collections ((1) national drug resistance survey, 2009-2010; (2) MDR strains from the Nhlangano region, 2014-2017). RESULTS: MDR strains in collection 1 had a high cluster rate (95%, 117/123 MDR strains) with 55% grouped into the two largest clusters (gCL3, n = 28; gCL10, n = 40). All gCL10 isolates, which likely emerged around 1993 (95% highest posterior density 1987-1998), carried the mutation RpoB I491F that is missed by commercial mDST assays. In addition, 21 (53%) gCL10 isolates shared a Rv0678 M146T mutation that correlated with elevated minimum inhibitory concentrations (MICs) to BDQ and CFZ compared to wild type isolates. gCL10 isolates with the Rv0678 M146T mutation were also detected in collection 2. CONCLUSION: The high clustering rate suggests that transmission has been driving the MDR-TB epidemic in Eswatini for three decades. The presence of MDR strains in Eswatini that are not detected by commercial mDST assays and have elevated MICs to BDQ and CFZ potentially jeopardizes the successful implementation of new MDR-TB treatment guidelines. Measures to limit the spread of these outbreak isolates need to be implemented urgently.

Rapid microarray-based assay for detection of pyrazinamide resistant Mycobacterium tuberculosis.

Pyrazinamide (PZA) is a key antibiotic for the treatment of drug susceptible tuberculosis. PZA-resistance is mainly mediated by mutations in the pncA gene; however the current gold standard is a phenotypic drug susceptibility test requiring a well-adjusted pH-value for reliable results. Our melting curve assay detects a non-wild type genotype in selected pncA regions in at least 3750 gene copies/mL within 2.5 hours. The prototype assay was further evaluated by analyzing 271 Mycobacterium tuberculosis complex isolates from Swaziland originating from a previously published drug resistance survey and including 118 isolates with pncA mutations. Sensitivity was 83% (95% CI 75-89%) and specificity was 100% (95% CI 98-100%). Under consideration of further improvements with regard to the target range our melting curve assay has the potential as a rapid rule-in test for PZA susceptibility (wild type pncA), however false resistant results (mutant pncA, but PZA susceptible) cannot be ruled out completely.

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.

The relationship between transmission time and clustering methods in Mycobacterium tuberculosis epidemiology.

BACKGROUND: Tracking recent transmission is a vital part of controlling widespread pathogens such as Mycobacterium tuberculosis. Multiple methods with specific performance characteristics exist for detecting recent transmission chains, usually by clustering strains based on genotype similarities. With such a large variety of methods available, informed selection of an appropriate approach for determining transmissions within a given setting/time period is difficult. METHODS: This study combines whole genome sequence (WGS) data derived from 324 isolates collected 2005-2010 in Kinshasa, Democratic Republic of Congo (DRC), a high endemic setting, with phylodynamics to unveil the timing of transmission events posited by a variety of standard genotyping methods. Clustering data based on Spoligotyping, 24-loci MIRU-VNTR typing, WGS based SNP (Single Nucleotide Polymorphism) and core genome multi locus sequence typing (cgMLST) typing were evaluated. FINDINGS: Our results suggest that clusters based on Spoligotyping could encompass transmission events that occurred almost 200 years prior to sampling while 24-loci-MIRU-VNTR often represented three decades of transmission. Instead, WGS based genotyping applying low SNP or cgMLST allele thresholds allows for determination of recent transmission events, e.g. in timespans of up to 10 years for a 5 SNP/allele cut-off. INTERPRETATION: With the rapid uptake of WGS methods in surveillance and outbreak tracking, the findings obtained in this study can guide the selection of appropriate clustering methods for uncovering relevant transmission chains within a given time-period. For high resolution cluster analyses, WGS-SNP and cgMLST based analyses have similar clustering/timing characteristics even for data obtained from a high incidence setting.

Outbreak of multidrug-resistant tuberculosis in South Africa undetected by WHO-endorsed commercial tests: an observational study.

BACKGROUND: Global roll-out of rapid molecular assays is revolutionising the diagnosis of rifampicin resistance, predictive of multidrug-resistance, in tuberculosis. However, 30% of the multidrug-resistant (MDR) strains in an eSwatini study harboured the Ile491Phe mutation in the rpoB gene, which is associated with poor rifampicin-based treatment outcomes but is missed by commercial molecular assays or scored as susceptible by phenotypic drug-susceptibility testing deployed in South Africa. We evaluated the presence of Ile491Phe among South African tuberculosis isolates reported as isoniazid-monoresistant according to current national testing algorithms. METHODS: We screened records of 37 644 Mycobacterium tuberculosis positive cultures from four South African provinces, diagnosed at the National Health Laboratory Service-Dr George Mukhari Tertiary Laboratory, to identify isolates with rifampicin sensitivity and isoniazid resistance according to Xpert MTB/RIF, GenoType MTBDRplus, and BACTEC MGIT 960. Of 1823 isolates that met these criteria, 277 were randomly selected and screened for Ile491Phe with multiplex allele-specific PCR and Sanger sequencing of rpoB. Ile491Phe-positive strains (as well as 17 Ile491Phe-bearing isolates from the eSwatini study) were then tested by Deeplex-MycTB deep sequencing and whole-genome sequencing to evaluate their patterns of extensive resistance, transmission, and evolution. FINDINGS: Ile491Phe was identified in 37 (15%) of 249 samples with valid multiplex allele-specific PCR and sequencing results, thus reclassifying them as MDR. All 37 isolates were additionally identified as genotypically resistant to all first-line drugs by Deeplex-MycTB. Six of the South African isolates harboured four distinct mutations potentially associated with decreased bedaquiline sensitivity. Consistent with Deeplex-MycTB genotypic profiles, whole-genome sequencing revealed concurrent silent spread in South Africa of a MDR tuberculosis strain lineage extending from the eSwatini outbreak and at least another independently emerged Ile491Phe-bearing lineage. Whole-genome sequencing further suggested acquisition of mechanisms compensating for the Ile491Phe fitness cost, and of additional bedaquiline resistance following the introduction of this drug in South Africa. INTERPRETATION: A substantial number of MDR tuberculosis cases harbouring the Ile491Phe mutation in the rpoB gene in South Africa are missed by current diagnostic strategies, resulting in ineffective first-line treatment, continued amplification of drug resistance, and concurrent silent spread in the community. FUNDING: VLIR-UOS, National Research Foundation (South Africa), and INNOVIRIS.

Comparative genomics of Mycobacterium africanum Lineage 5 and Lineage 6 from Ghana suggests distinct ecological niches.

Mycobacterium africanum (Maf) causes a substantial proportion of human tuberculosis in some countries of West Africa, but little is known on this pathogen. We compared the genomes of 253 Maf clinical isolates from Ghana, including N = 175 Lineage 5 (L5) and N = 78 Lineage 6 (L6). We found that the genomic diversity of L6 was higher than in L5 despite the smaller sample size. Regulatory proteins appeared to evolve neutrally in L5 but under purifying selection in L6. Even though over 90% of the human T cell epitopes were conserved in both lineages, L6 showed a higher ratio of non-synonymous to synonymous single nucleotide variation in these epitopes overall compared to L5. Of the 10% human T cell epitopes that were variable, most carried mutations that were lineage-specific. Our findings indicate that Maf L5 and L6 differ in some of their population genomic characteristics, possibly reflecting different selection pressures linked to distinct ecological niches.