Refined understanding of the impact of the Mycobacterium tuberculosis complex diversity on the intrinsic susceptibility to pretomanid

Previous work reported unprecedented differences in the intrinsic in vitro susceptibility of the Mycobacterium tuberculosis complex (MTBC) to pretomanid (Pa) using the Mycobacteria Growth Indicator Tube (MGIT) system. We tested 125 phylogenetically diverse strains from all known MTBC lineages (1-9) without known Pa resistance mutations and four strains with known resistance mutations as controls. This confirmed that MTBC, unlike most bacteria-antimicrobial combinations, displayed substantial differences in the intrinsic susceptibility relative to the technical variation of Pa MIC testing. This was also the case for the Middlebrook 7H11 (7H11) medium, demonstrating that these differences were not specific to MGIT. Notably, lineage 1 was confirmed to have intrinsically elevated MICs compared with lineages 2, 3, 4, and 7 (L2-4/7), underlining the urgent need for WHO to publish its decision of whether lineage 1 should be deemed treatable by BPaL(M), the now preferred all-oral regimen for treating rifampin-resistant tuberculosis. Lineages 5 and 6, which are most frequent in West Africa, responded differently to Pa, with lineage 5 being more similar to L2-4/7 and lineage 6 being more susceptible. More data are needed to determine whether 7H11 MICs are systematically lower than those in MGIT.

IMPORTANCE

This study confirmed that the Mycobacterium tuberculosis complex lineage 1, responsible for 28% of global tuberculosis cases, is less susceptible to pretomanid (Pa). It also refined the understanding of the intrinsic susceptibilities of lineages 5 and 6, most frequent in West Africa, and lineages 8 and 9. Regulators must review whether these in vitro differences affect the clinical efficacy of the WHO-recommended BPaL(M) regimen and set breakpoints for antimicrobial susceptibility testing accordingly. Notably, regulators should provide detailed justifications for their decisions to facilitate public scrutiny.

Quantitative measurement of antibiotic resistance in Mycobacterium tuberculosis reveals genetic determinants of resistance and susceptibility in a target gene approach

The World Health Organization has a goal of universal drug susceptibility testing for patients with tuberculosis. However, molecular diagnostics to date have focused largely on first-line drugs and predicting susceptibilities in a binary manner (classifying strains as either susceptible or resistant). Here, we used a multivariable linear mixed model alongside whole genome sequencing and a quantitative microtiter plate assay to relate genomic mutations to minimum inhibitory concentration (MIC) in 15,211 Mycobacterium tuberculosis clinical isolates from 23 countries across five continents. We identified 492 unique MIC-elevating variants across 13 drugs, as well as 91 mutations likely linked to hypersensitivity. Our results advance genetics-based diagnostics for tuberculosis and serve as a curated training/testing dataset for development of drug resistance prediction algorithms.

Two-stage tuberculosis diagnostics: combining centrifugal microfluidics to detect TB infection and Inh and Rif resistance at the point of care with subsequent antibiotic resistance profiling by targeted NGS

Globally, tuberculosis (TB) remains the deadliest bacterial infectious disease, and spreading antibiotic resistances is the biggest challenge for combatting the disease. Rapid and comprehensive diagnostics including drug susceptibility testing (DST) would assure early treatment, reduction of morbidity and the interruption of transmission chains. To date, rapid genetic resistance testing addresses only one to four drug groups while complete DST is done phenotypically and takes several weeks. To overcome these limitations, we developed a two-stage workflow for rapid TB diagnostics including DST from a single sputum sample that can be completed within three days. The first stage is qPCR detection of M. tuberculosis complex (MTBC) including antibiotic resistance testing against the first-line antibiotics, isoniazid (Inh) and rifampicin (Rif). The test is automated by centrifugal microfluidics and designed for point of care (PoC). Furthermore, enriched MTBC DNA is provided in a detachable sample tube to enable the second stage: if the PCR detects MTBC and resistance to either Inh or Rif, the MTBC DNA is shipped to specialized facilities and analyzed by targeted next generation sequencing (tNGS) to assess the complete resistance profile. Proof-of-concept testing of the PoC test revealed an analytical sensitivity of 44.2 CFU ml-1, a diagnostic sensitivity of 96%, and a diagnostic specificity of 100% for MTBC detection. Coupled tNGS successfully provided resistance profiles, demonstrated for samples from 17 patients. To the best of our knowledge, the presented combination of PoC qPCR with tNGS allows for the fastest comprehensive TB diagnostics comprising decentralized pathogen detection with subsequent resistance profiling in a facility specialized in tNGS.

In vitro activity of new combinations of β-lactam and β-lactamase inhibitors against the Mycobacterium tuberculosis complex

As meropenem-clavulanic acid is recommended for the treatment of drug-resistant tuberculosis, the repurposing of new carbapenem combinations may provide new treatment options, including oral alternatives. Therefore, we studied the in vitro activities of meropenem-vaborbactam, meropenem-clavulanic acid, and tebipenem-clavulanic acid. One hundred nine Mycobacterium tuberculosis complex (MTBC) clinical isolates were tested, of which 69 were pan-susceptible and the remaining pyrazinamide- or multidrug-resistant. Broth microdilution MICs were determined using the EUCAST reference method. Meropenem and tebipenem were tested individually and in combination with vaborbactam 8 mg/L and clavulanic-acid 2 and 4 mg/L, respectively. Whole-genome sequencing was performed to explore resistance mechanisms. Clavulanic acid lowered the modal tebipenem MIC approximately 16-fold (from 16 to 1 mg/L). The modal meropenem MIC was reduced twofold by vaborbactam compared with an approximately eightfold decrease by clavulanic acid. The only previously described high-confidence carbapenem resistance mutation, crfA T62A, was shared by a subgroup of lineage 4.3.4.1 isolates and did not correlate with elevated MICs. The presence of a β-lactamase inhibitor reduced the MTBC MICs of tebipenem and meropenem. The resulting MIC distribution was lowest for the orally available drugs tebipenem-clavulanic acid. Whether this in vitro activity translates to similar or greater clinical efficacy of tebipenem-clavulanic acid compared with the currently WHO-endorsed meropenem-clavulanic acid requires clinical studies. IMPORTANCE Repurposing of already approved antibiotics, such as β-lactams in combination with β-lactamase inhibitors, may provide new treatment alternatives for drug-resistant tuberculosis. Meropenem-clavulanic acid was more active in vitro compared to meropenem-vaborbactam. Notably, tebipenem-clavulanic acid showed even better activity, raising the potential of an all-oral treatment option. Clinical data are needed to investigate whether the better in vitro activity of tebipenem-clavulanic acid correlates with greater clinical efficacy compared with the currently WHO-endorsed meropenem-clavulanic acid.

Transmission Dynamics of a Mycobacterium tuberculosis Complex Outbreak in an Indigenous Population in the Colombian Amazon Region

Whole genome sequencing (WGS) has become the main tool for studying the transmission of Mycobacterium tuberculosis complex (MTBC) strains; however, the clonal expansion of one strain often limits its application in local MTBC outbreaks. The use of an alternative reference genome and the inclusion of repetitive regions in the analysis could potentially increase the resolution, but the added value has not yet been defined. Here, we leveraged short and long WGS read data of a previously reported MTBC outbreak in the Colombian Amazon Region to analyze possible transmission chains among 74 patients in the indigenous setting of Puerto Nariño (March to October 2016). In total, 90.5% (67/74) of the patients were infected with one distinct MTBC strain belonging to lineage 4.3.3. Employing a reference genome from an outbreak strain and highly confident single nucleotide polymorphisms (SNPs) in repetitive genomic regions, e.g., the proline-glutamic acid/proline-proline-glutamic-acid (PE/PPE) gene family, increased the phylogenetic resolution compared to a classical H37Rv reference mapping approach. Specifically, the number of differentiating SNPs increased from 890 to 1,094, which resulted in a more granular transmission network as judged by an increasing number of individual nodes in a maximum parsimony tree, i.e., 5 versus 9 nodes. We also found in 29.9% (20/67) of the outbreak isolates, heterogenous alleles at phylogenetically informative sites, suggesting that these patients are infected with more than one clone. In conclusion, customized SNP calling thresholds and employment of a local reference genome for a mapping approach can improve the phylogenetic resolution in highly clonal MTBC populations and help elucidate within-host MTBC diversity. IMPORTANCE The Colombian Amazon around Puerto Nariño has a high tuberculosis burden with a prevalence of 1,267/100,000 people in 2016. Recently, an outbreak of Mycobacterium tuberculosis complex (MTBC) bacteria among the indigenous populations was identified with classical MTBC genotyping methods. Here, we employed a whole-genome sequencing-based outbreak investigation in order to improve the phylogenetic resolution and gain new insights into the transmission dynamics in this remote Colombian Amazon Region. The inclusion of well-supported single nucleotide polymorphisms in repetitive regions and a de novo-assembled local reference genome provided a more granular picture of the circulating outbreak strain and revealed new transmission chains. Multiple patients from different settlements were possibly infected with at least two different clones in this high-incidence setting. Thus, our results have the potential to improve molecular surveillance studies in other high-burden settings, especially regions with few clonal multidrug-resistant (MDR) MTBC lineages/clades.

Limited Nosocomial Transmission of Drug-Resistant Tuberculosis, Moldova

Applying whole-genome-sequencing, we aimed to detect transmission events of multidrug-resistant/rifampin-resistant strains of Mycobacterium tuberculosis complex at a tuberculosis hospital in Chisinau, Moldova. We recorded ward, room, and bed information for each patient and monitored in-hospital transfers over 1 year. Detailed molecular and patient surveillance revealed only 2 nosocomial transmission events.

Bedaquiline and clofazimine resistance in Mycobacterium tuberculosis: an in-vitro and in-silico data analysis

BACKGROUND

Bedaquiline is a core drug for the treatment of multidrug-resistant tuberculosis; however, the understanding of resistance mechanisms is poor, which is hampering rapid molecular diagnostics. Some bedaquiline-resistant mutants are also cross-resistant to clofazimine. To decipher bedaquiline and clofazimine resistance determinants, we combined experimental evolution, protein modelling, genome sequencing, and phenotypic data.

METHODS

For this in-vitro and in-silico data analysis, we used a novel in-vitro evolutionary model using subinhibitory drug concentrations to select bedaquiline-resistant and clofazimine-resistant mutants. We determined bedaquiline and clofazimine minimum inhibitory concentrations and did Illumina and PacBio sequencing to characterise selected mutants and establish a mutation catalogue. This catalogue also includes phenotypic and genotypic data of a global collection of more than 14 000 clinical Mycobacterium tuberculosis complex isolates, and publicly available data. We investigated variants implicated in bedaquiline resistance by protein modelling and dynamic simulations.

FINDINGS

We discerned 265 genomic variants implicated in bedaquiline resistance, with 250 (94%) variants affecting the transcriptional repressor (Rv0678) of the MmpS5-MmpL5 efflux system. We identified 40 new variants in vitro, and a new bedaquiline resistance mechanism caused by a large-scale genomic rearrangement. Additionally, we identified in vitro 15 (7%) of 208 mutations found in clinical bedaquiline-resistant isolates. From our in-vitro work, we detected 14 (16%) of 88 mutations so far identified as being associated with clofazimine resistance and also seen in clinically resistant strains, and catalogued 35 new mutations. Structural modelling of Rv0678 showed four major mechanisms of bedaquiline resistance: impaired DNA binding, reduction in protein stability, disruption of protein dimerisation, and alteration in affinity for its fatty acid ligand.

INTERPRETATION

Our findings advance the understanding of drug resistance mechanisms in M tuberculosis complex strains. We have established an extended mutation catalogue, comprising variants implicated in resistance and susceptibility to bedaquiline and clofazimine. Our data emphasise that genotypic testing can delineate clinical isolates with borderline phenotypes, which is essential for the design of effective treatments.

FUNDING

Leibniz ScienceCampus Evolutionary Medicine of the Lung, Deutsche Forschungsgemeinschaft, Research Training Group 2501 TransEvo, Rhodes Trust, Stanford University Medical Scientist Training Program, National Institute for Health and Care Research Oxford Biomedical Research Centre, Oxford University Hospitals NHS Foundation Trust, Bill & Melinda Gates Foundation, Wellcome Trust, and Marie Skłodowska-Curie Actions.

Clinical implications of molecular drug resistance testing for Mycobacterium tuberculosis: a 2023 TBnet/RESIST-TB consensus statement

Drug-resistant tuberculosis is a substantial health-care concern worldwide. Despite culture-based methods being considered the gold standard for drug susceptibility testing, molecular methods provide rapid information about the Mycobacterium tuberculosis mutations associated with resistance to anti-tuberculosis drugs. This consensus document was developed on the basis of a comprehensive literature search, by the TBnet and RESIST-TB networks, about reporting standards for the clinical use of molecular drug susceptibility testing. Review and the search for evidence included hand-searching journals and searching electronic databases. The panel identified studies that linked mutations in genomic regions of M tuberculosis with treatment outcome data. Implementation of molecular testing for the prediction of drug resistance in M tuberculosis is key. Detection of mutations in clinical isolates has implications for the clinical management of patients with multidrug-resistant or rifampicin-resistant tuberculosis, especially in situations when phenotypic drug susceptibility testing is not available. A multidisciplinary team including clinicians, microbiologists, and laboratory scientists reached a consensus on key questions relevant to molecular prediction of drug susceptibility or resistance to M tuberculosis, and their implications for clinical practice. This consensus document should help clinicians in the management of patients with tuberculosis, providing guidance for the design of treatment regimens and optimising outcomes.

Investigating resistance in clinical Mycobacterium tuberculosis complex isolates with genomic and phenotypic antimicrobial susceptibility testing: a multicentre observational study

BACKGROUND

Whole-genome sequencing (WGS) of Mycobacterium tuberculosis complex has become an important tool in diagnosis and management of drug-resistant tuberculosis. However, data correlating resistance genotype with quantitative phenotypic antimicrobial susceptibility testing (AST) are scarce.

METHODS

In a prospective multicentre observational study, 900 clinical M tuberculosis complex isolates were collected from adults with drug-resistant tuberculosis in five high-endemic tuberculosis settings around the world (Georgia, Moldova, Peru, South Africa, and Viet Nam) between Dec 5, 2014, and Dec 12, 2017. Minimum inhibitory concentrations (MICs) and resulting binary phenotypic AST results for up to nine antituberculosis drugs were determined and correlated with resistance-conferring mutations identified by WGS.

FINDINGS

Considering WHO-endorsed critical concentrations as reference, WGS had high accuracy for prediction of resistance to isoniazid (sensitivity 98·8% [95% CI 98·5-99·0]; specificity 96·6% [95% CI 95·2-97·9]), levofloxacin (sensitivity 94·8% [93·3-97·6]; specificity 97·1% [96·7-97·6]), kanamycin (sensitivity 96·1% [95·4-96·8]; specificity 95·0% [94·4-95·7]), amikacin (sensitivity 97·2% [96·4-98·1]; specificity 98·6% [98·3-98·9]), and capreomycin (sensitivity 93·1% [90·0-96·3]; specificity 98·3% [98·0-98·7]). For rifampicin, pyrazinamide, and ethambutol, the specificity of resistance prediction was suboptimal (64·0% [61·0-67·1], 83·8% [81·0-86·5], and 40·1% [37·4-42·9], respectively). Specificity for rifampicin increased to 83·9% when borderline mutations with MICs overlapping with the critical concentration were excluded. Consequently, we highlighted mutations in M tuberculosis complex isolates that are often falsely identified as susceptible by phenotypic AST, and we identified potential novel resistance-conferring mutations.

INTERPRETATION

The combined analysis of mutations and quantitative phenotypes shows the potential of WGS to produce a refined interpretation of resistance, which is needed for individualised therapy, and eventually could allow differential drug dosing. However, variability of MIC data for some M tuberculosis complex isolates carrying identical mutations also reveals limitations of our understanding of the genotype and phenotype relationships (eg, including epistasis and strain genetic background).

FUNDING

Bill & Melinda Gates Foundation, German Centre for Infection Research, German Research Foundation, Excellence Cluster Precision Medicine of Inflammation (EXC 2167), and Leibniz ScienceCampus EvoLUNG.

Transcontinental spread and evolution of Mycobacterium tuberculosis W148 European/Russian clade toward extensively drug resistant tuberculosis

Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a “perfect storm” that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.

An outbreak of the multidrug-resistant Mycobacterium tuberculosis lineage W148 has spread widely across Russia, Central Asia and Europe. Here, the authors use whole genome sequences of ~700 isolates of this lineage collected over ~20 years to analyze its spread, evolution of drug resistance, and impact of compensatory mutations.

High fluoroquinolone resistance proportions among multidrug-resistant tuberculosis driven by dominant L2 Mycobacterium tuberculosis clones in the Mumbai Metropolitan Region

Background

Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains are a serious health problem in India, also contributing to one-fourth of the global MDR tuberculosis (TB) burden. About 36% of the MDR MTBC strains are reported fluoroquinolone (FQ) resistant leading to high pre-extensively drug-resistant (pre-XDR) and XDR-TB (further resistance against bedaquiline and/or linezolid) rates. Still, factors driving the MDR/pre-XDR epidemic in India are not well defined.

Methods

In a retrospective study, we analyzed 1852 consecutive MTBC strains obtained from patients from a tertiary care hospital laboratory in Mumbai by whole genome sequencing (WGS). Univariate and multivariate statistics was used to investigate factors associated with pre-XDR. Core genome multi locus sequence typing, time scaled haplotypic density (THD) method and homoplasy analysis were used to analyze epidemiological success, and positive selection in different strain groups, respectively.

Results

In total, 1016 MTBC strains were MDR, out of which 703 (69.2%) were pre-XDR and 45 (4.4%) were XDR. Cluster rates were high among MDR (57.8%) and pre-XDR/XDR (79%) strains with three dominant L2 (Beijing) strain clusters (Cl 1–3) representing half of the pre-XDR and 40% of the XDR-TB cases. L2 strains were associated with pre-XDR/XDR-TB (P < 0.001) and, particularly Cl 1–3 strains, had high first-line and FQ resistance rates (81.6–90.6%). Epidemic success analysis using THD showed that L2 strains outperformed L1, L3, and L4 strains in short- and long-term time scales. More importantly, L2 MDR and MDR + strains had higher THD success indices than their not-MDR counterparts. Overall, compensatory mutation rates were highest in L2 strains and positive selection was detected in genes of L2 strains associated with drug tolerance (prpB and ppsA) and virulence (Rv2828c). Compensatory mutations in L2 strains were associated with a threefold increase of THD indices, suggesting improved transmissibility.

Conclusions

Our data indicate a drastic increase of FQ resistance, as well as emerging bedaquiline resistance which endangers the success of newly endorsed MDR-TB treatment regimens. Rapid changes in treatment and control strategies are required to contain transmission of highly successful pre-XDR L2 strains in the Mumbai Metropolitan region but presumably also India-wide.

Prediction of drug resistance by Sanger sequencing of Mycobacterium tuberculosis complex strains isolated from multidrug resistant tuberculosis suspect patients in Ethiopia

Background

Ethiopia is one of the high multidrug-resistant tuberculosis (MDR-TB) burden countries. However, phenotypic drug susceptibility testing can take several weeks due to the slow growth of Mycobacterium tuberculosis complex (MTBC) strains. In this study, we assessed the performance of a Sanger sequencing approach to predict resistance against five anti-tuberculosis drugs and the pattern of resistance mediating mutations.

Methods

We enrolled 226 MTBC culture-positive MDR-TB suspects and collected sputum specimens and socio-demographic and TB related data from each suspect between June 2015 and December 2016 in Addis Ababa, Ethiopia. Phenotypic drug susceptibility testing (pDST) for rifampicin, isoniazid, pyrazinamide, ethambutol, and streptomycin using BACTEC MGIT 960 was compared with the results of a Sanger sequencing analysis of seven resistance determining regions in the genes rpoB, katG, fabG-inhA, pncA, embB, rpsL, and rrs.

Result

DNA isolation for Sanger sequencing was successfully extracted from 92.5% (209/226) of the MTBC positive cultures, and the remaining 7.5% (17/226) strains were excluded from the final analysis. Based on pDST results, drug resistance proportions were as follows: isoniazid: 109/209 (52.2%), streptomycin: 93/209 (44.5%), rifampicin: 88/209 (42.1%), ethambutol: 74/209 (35.4%), and pyrazinamide: 69/209 (33.0%). Resistance against isoniazid was mainly mediated by the mutation katG S315T (97/209, 46.4%) and resistance against rifampicin by rpoB S531L (58/209, 27.8%). The dominating resistance-conferring mutations for ethambutol, streptomycin, and pyrazinamide affected codon 306 in embB (48/209, 21.1%), codon 88 in rpsL (43/209, 20.6%), and codon 65 in pncA (19/209, 9.1%), respectively. We observed a high agreement between phenotypic and genotypic DST, such as 89.9% (at 95% confidence interval [CI], 84.2%–95.8%) for isoniazid, 95.5% (95% CI, 91.2%–99.8%) for rifampicin, 98.6% (95% CI, 95.9–100%) for ethambutol, 91.3% (95% CI, 84.6–98.1%) for pyrazinamide and 57.0% (95% CI, 46.9%–67.1%) for streptomycin.

Conclusion

We detected canonical mutations implicated in resistance to rifampicin, isoniazid, pyrazinamide, ethambutol, and streptomycin. High agreement with phenotypic DST results for all drugs renders Sanger sequencing promising to be performed as a complementary measure to routine phenotypic DST in Ethiopia. Sanger sequencing directly from sputum may accelerate accurate clinical decision-making in the future.

Updating the approaches to define susceptibility and resistance to anti-tuberculosis agents: implications for diagnosis and treatment

Approximately 85 000 deaths globally in 2019 were due to drug-resistant tuberculosis (TB), which corresponds to 7% of global deaths attributable to bacterial antimicrobial resistance [1]. Yet concerns have been mounting that drug-resistant TB was being underestimated because the approaches to define susceptibility and resistance to anti-TB agents had not kept up with those used for other major bacterial pathogens [2–9]. Here, we outline the recent, evidence-based initiatives spearheaded by the World Health Organization (WHO) and others to update breakpoints (traditionally referred to as critical concentrations (CCs)) that are used for phenotypic antimicrobial susceptibility testing (AST), also called drug susceptibility testing in the TB literature.

Inappropriately high breakpoints have resulted in systematic false-susceptible AST results to anti-TB drugs. MIC, PK/PD and clinical outcome data should be combined when setting breakpoints to minimise the emergence and spread of antimicrobial resistance. https://bit.ly/3i43wb6

Emergence of bedaquiline resistance in a high tuberculosis burden country

RATIONALE

Bedaquiline has been classified as a group A drug for the treatment of multidrug-resistant tuberculosis (MDR-TB) by the World Health Organization; however, globally emerging resistance threatens the effectivity of novel MDR-TB treatment regimens.

OBJECTIVES

We analysed pre-existing and emerging bedaquiline resistance in bedaquiline-based MDR-TB therapies, and risk factors associated with treatment failure and death.

METHODS

In a cross-sectional cohort study, we employed patient data, whole-genome sequencing (WGS) and phenotyping of Mycobacterium tuberculosis complex (MTBC) isolates. We could retrieve baseline isolates from 30.5% (62 out of 203) of all MDR-TB patients who received bedaquiline between 2016 and 2018 in the Republic of Moldova. This includes 26 patients for whom we could also retrieve a follow-up isolate.

MEASUREMENTS AND MAIN RESULTS

At baseline, all MTBC isolates were susceptible to bedaquiline. Among 26 patients with available baseline and follow-up isolates, four (15.3%) patients harboured strains which acquired bedaquiline resistance under therapy, while one (3.8%) patient was re-infected with a second bedaquiline-resistant strain. Treatment failure and death were associated with cavitary disease (p=0.011), and any additional drug prescribed in the bedaquiline-containing regimen with WGS-predicted resistance at baseline (OR 1.92 per unit increase, 95% CI 1.15-3.21; p=0.012).

CONCLUSIONS

MDR-TB treatments based on bedaquiline require a functional background regimen to achieve high cure rates and to prevent the evolution of bedaquiline resistance. Novel MDR-TB therapies with bedaquiline require timely and comprehensive drug resistance monitoring.

Whole genome sequencing-based classification of human-related Haemophilus species and detection of antimicrobial resistance genes

Background

Bacteria belonging to the genus Haemophilus cause a wide range of diseases in humans. Recently, H. influenzae was classified by the WHO as priority pathogen due to the wide spread of ampicillin resistant strains. However, other Haemophilus spp. are often misclassified as H. influenzae. Therefore, we established an accurate and rapid whole genome sequencing (WGS) based classification and serotyping algorithm and combined it with the detection of resistance genes.

Methods

A gene presence/absence-based classification algorithm was developed, which employs the open-source gene-detection tool SRST2 and a new classification database comprising 36 genes, including capsule loci for serotyping. These genes were identified using a comparative genome analysis of 215 strains belonging to ten human-related Haemophilus (sub)species (training dataset). The algorithm was evaluated on 1329 public short read datasets (evaluation dataset) and used to reclassify 262 clinical Haemophilus spp. isolates from 250 patients (German cohort). In addition, the presence of antibiotic resistance genes within the German dataset was evaluated with SRST2 and correlated with results of traditional phenotyping assays.

Results

The newly developed algorithm can differentiate between clinically relevant Haemophilus species including, but not limited to, H. influenzae, H. haemolyticus, and H. parainfluenzae. It can also identify putative haemin-independent H. haemolyticus strains and determine the serotype of typeable Haemophilus strains. The algorithm performed excellently in the evaluation dataset (99.6% concordance with reported species classification and 99.5% with reported serotype) and revealed several misclassifications. Additionally, 83 out of 262 (31.7%) suspected H. influenzae strains from the German cohort were in fact H. haemolyticus strains, some of which associated with mouth abscesses and lower respiratory tract infections.

Resistance genes were detected in 16 out of 262 datasets from the German cohort. Prediction of ampicillin resistance, associated with bla_TEM-1D, and tetracycline resistance, associated with tetB, correlated well with available phenotypic data.

Conclusions

Our new classification database and algorithm have the potential to improve diagnosis and surveillance of Haemophilus spp. and can easily be coupled with other public genotyping and antimicrobial resistance databases. Our data also point towards a possible pathogenic role of H. haemolyticus strains, which needs to be further investigated.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01017-x.

OUP accepted manuscript

Background

Rifampicin- or multidrug-resistant (RR/MDR) Mycobacterium tuberculosis complex (MTBC) strains account for considerable morbidity and mortality globally. WGS-based prediction of drug resistance may guide clinical decisions, especially for the design of RR/MDR-TB therapies.

Methods

We compared WGS-based drug resistance-predictive mutations for 42 MTBC isolates from MDR-TB patients in Tanzania with the MICs of 14 antibiotics measured in the Sensititre™ MycoTB assay. An isolate was phenotypically categorized as resistant if it had an MIC above the epidemiological-cut-off (ECOFF) value, or as susceptible if it had an MIC below or equal to the ECOFF.

Results

Overall, genotypically non-wild-type MTBC isolates with high-level resistance mutations (gNWT-R) correlated with isolates with MIC values above the ECOFF. For instance, the median MIC value (mg/L) for rifampicin-gNWT-R strains was >4.0 (IQR 4.0–4.0) compared with 0.5 (IQR 0.38–0.50) in genotypically wild-type (gWT-S, P < 0.001); isoniazid-gNWT-R >4.0 (IQR 2.0–4.0) compared with 0.25 (IQR 0.12–1.00) among gWT-S (P = 0.001); ethionamide-gNWT-R 15.0 (IQR 10.0–20.0) compared with 2.50 (IQR; 2.50–5.00) among gWT-S (P < 0.001). WGS correctly predicted resistance in 95% (36/38) and 100% (38/38) of the rifampicin-resistant isolates with ECOFFs >0.5 and >0.125 mg/L, respectively. No known resistance-conferring mutations were present in genes associated with resistance to fluoroquinolones, aminoglycosides, capreomycin, bedaquiline, delamanid, linezolid, clofazimine, cycloserine, or p-amino salicylic acid.

Conclusions

WGS-based drug resistance prediction worked well to rule-in phenotypic drug resistance and the absence of second-line drug resistance-mediating mutations has the potential to guide the design of RR/MDR-TB regimens in the future.

OUP accepted manuscript

Objectives

To develop a robust phenotypic antimicrobial susceptibility testing (AST) method with a correctly set breakpoint for pretomanid (Pa), the most recently approved anti-tuberculosis drug.

Methods

The Becton Dickinson Mycobacterial Growth Indicator Tube™ (MGIT) system was used at six laboratories to determine the MICs of a phylogenetically diverse collection of 356 Mycobacterium tuberculosis complex (MTBC) strains to establish the epidemiological cut-off value for pretomanid. MICs were correlated with WGS data to study the genetic basis of differences in the susceptibility to pretomanid.

Results

We observed ancient differences in the susceptibility to pretomanid among various members of MTBC. Most notably, lineage 1 of M. tuberculosis, which is estimated to account for 28% of tuberculosis cases globally, was less susceptible than lineages 2, 3, 4 and 7 of M. tuberculosis, resulting in a 99th percentile of 2 mg/L for lineage 1 compared with 0.5 mg/L for the remaining M. tuberculosis lineages. Moreover, we observed that higher MICs (≥8 mg/L), which probably confer resistance, had recently evolved independently in six different M. tuberculosis strains. Unlike the aforementioned ancient differences in susceptibility, these recent differences were likely caused by mutations in the known pretomanid resistance genes.

Conclusions

In light of these findings, the provisional critical concentration of 1 mg/L for MGIT set by EMA must be re-evaluated. More broadly, these findings underline the importance of considering the global diversity of MTBC during clinical development of drugs and when defining breakpoints for AST.

Evaluation of whole-genome sequence data analysis approaches for short- and long-read sequencing of Mycobacterium tuberculosis

Whole-genome sequencing (WGS) of Mycobacterium tuberculosis (MTB) isolates can be used to get an accurate diagnosis, to guide clinical decision making, to control tuberculosis (TB) and for outbreak investigations. We evaluated the performance of long-read (LR) and/or short-read (SR) sequencing for anti-TB drug-resistance prediction using the TBProfiler and Mykrobe tools, the fraction of genome recovery, assembly accuracies and the robustness of two typing approaches based on core-genome SNP (cgSNP) typing and core-genome multi-locus sequence typing (cgMLST). Most of the discrepancies between phenotypic drug-susceptibility testing (DST) and drug-resistance prediction were observed for the first-line drugs rifampicin, isoniazid, pyrazinamide and ethambutol, mainly with LR sequence data. Resistance prediction to second-line drugs made by both TBProfiler and Mykrobe tools with SR- and LR-sequence data were in complete agreement with phenotypic DST except for one isolate. The SR assemblies were more accurate than the LR assemblies, having significantly (P <0.05) fewer indels and mismatches per 100 kbp. However, the hybrid and LR assemblies had slightly higher genome fractions. For LR assemblies, Canu followed by Racon, and Medaka polishing was the most accurate approach. The cgSNP approach, based on either reads or assemblies, was more robust than the cgMLST approach, especially for LR sequence data. In conclusion, anti-TB drug-resistance prediction, particularly with only LR sequence data, remains challenging, especially for first-line drugs. In addition, SR assemblies appear more accurate than LR ones, and reproducible phylogeny can be achieved using cgSNP approaches.

Transmission patterns of rifampicin resistant Mycobacterium tuberculosis complex strains in Cameroon: a genomic epidemiological study

Background

Determining factors affecting the transmission of rifampicin (RR) and multidrug-resistant (MDR) Mycobacterium tuberculosis complex strains under standardized tuberculosis (TB) treatment is key to control TB and prevent the evolution of drug resistance.

Methods

We combined bacterial whole genome sequencing (WGS) and epidemiological investigations for 37% (n = 195) of all RR/MDR-TB patients in Cameroon (2012–2015) to identify factors associated with recent transmission.

Results

Patients infected with a strain resistant to high-dose isoniazid, and ethambutol had 7.4 (95% CI 2.6–21.4), and 2.4 (95% CI 1.2–4.8) times increased odds of being in a WGS-cluster, a surrogate for recent transmission. Furthermore, age between 30 and 50 was positively correlated with recent transmission (adjusted OR 3.8, 95% CI 1.3–11.4). We found high drug-resistance proportions against three drugs used in the short standardized MDR-TB regimen in Cameroon, i.e. high-dose isoniazid (77.4%), ethambutol (56.9%), and pyrazinamide (43.1%). Virtually all strains were susceptible to fluoroquinolones, kanamycin, and clofazimine, and treatment outcomes were mostly favourable (87.5%).

Conclusion

Pre-existing resistance to high-dose isoniazid, and ethambutol is associated with recent transmission of RR/MDR strains in our study. A possible contributing factor for this observation is the absence of universal drug susceptibility testing in Cameroon, likely resulting in prolonged exposure of new RR/MDR-TB patients to sub-optimal or failing first-line drug regimens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-06593-8.

Microevolution of Mycobacterium tuberculosis Subpopulations and Heteroresistance in a Patient Receiving 27 Years of Tuberculosis Treatment in Germany

Preexisting and newly emerging resistant pathogen subpopulations (heteroresistance) are potential risk factors for treatment failure of multi/extensively drug resistant (MDR/XDR) tuberculosis (TB). Intrapatient evolutionary dynamics of Mycobacterium tuberculosis complex (Mtbc) strains and their implications on treatment outcomes are still not completely understood. To elucidate how Mtbc strains escape therapy, we analyzed 13 serial isolates from a German patient by whole-genome sequencing. Sequencing data were compared with phenotypic drug susceptibility profiles and the patient's collective 27-year treatment history to further elucidate factors fostering intrapatient resistance evolution. The patient endured five distinct TB episodes, ending in resistance to 16 drugs and a nearly untreatable XDR-TB infection. The first isolate obtained, during the patient's 5th TB episode, presented fixed resistance mutations to 7 anti-TB drugs, including isoniazid, rifampin, streptomycin, pyrazinamide, prothionamide, para-aminosalicylic acid, and cycloserine-terizidone. Over the next 13 years, a dynamic evolution with coexisting, heterogeneous subpopulations was observed in 6 out of 13 sequential bacterial isolates. The emergence of drug-resistant subpopulations coincided with frequent changes in treatment regimens, which often included two or fewer active compounds. This evolutionary arms race between competing subpopulations ultimately resulted in the fixation of a single XDR variant. Our data demonstrate the complex intrapatient microevolution of Mtbc subpopulations during failing MDR/XDR-TB treatment. Designing effective treatment regimens based on rapid detection of (hetero) resistance is key to avoid resistance development and treatment failure.

Perspectives for systems biology in the management of tuberculosis

Standardised management of tuberculosis may soon be replaced by individualised, precision medicine-guided therapies informed with knowledge provided by the field of systems biology. Systems biology is a rapidly expanding field of computational and mathematical analysis and modelling of complex biological systems that can provide insights into mechanisms underlying tuberculosis, identify novel biomarkers, and help to optimise prevention, diagnosis and treatment of disease. These advances are critically important in the context of the evolving epidemic of drug-resistant tuberculosis. Here, we review the available evidence on the role of systems biology approaches - human and mycobacterial genomics and transcriptomics, proteomics, lipidomics/metabolomics, immunophenotyping, systems pharmacology and gut microbiomes - in the management of tuberculosis including prediction of risk for disease progression, severity of mycobacterial virulence and drug resistance, adverse events, comorbidities, response to therapy and treatment outcomes. Application of the Grading of Recommendations, Assessment, Development and Evaluation (GRADE) approach demonstrated that at present most of the studies provide "very low" certainty of evidence for answering clinically relevant questions. Further studies in large prospective cohorts of patients, including randomised clinical trials, are necessary to assess the applicability of the findings in tuberculosis prevention and more efficient clinical management of patients.

Design of multidrug-resistant tuberculosis treatment regimens based on DNA sequencing

BACKGROUND

Comprehensive and reliable drug susceptibility testing (DST) is urgently needed to provide adequate treatment regimens for patients with multidrug-resistant/rifampicin-resistant tuberculosis (MDR/RR-TB). We investigated if next generation sequencing (NGS) analysis of Mycobacterium tuberculosis complex isolates and genes implicated in drug resistance can guide the design of effective MDR/RR-TB treatment regimens.

METHODS

NGS-based genomic DST predictions of M. tuberculosis complex isolates from MDR/RR-TB patients admitted to a TB reference center in Germany between 01/01/2015 and 04/30/2019 were compared with phenotypic DST results of Mycobacteria growth indicator tubes (MGIT). Standardized treatment algorithms were applied to design individualized therapies based on either genomic or phenotypic DST results, and discrepancies were further evaluated by determination of minimum inhibitory drug concentrations (MIC) using Sensititre MYCOTBI and UKMYC microtiter plates.

RESULTS

In 70 patients with MDR/RR-TB, agreement among 1048 pairwise comparisons of genomic and phenotypic DST was 86.3%; 76 (7.2%) results were discordant, and 68 (6.5%) could not be evaluated due to presence of polymorphisms with yet unknown implications for drug resistance. Importantly, 549/561 (97.9%) predictions of drug susceptibility were phenotypically confirmed in MGIT, and 27/64 (42.2%) false positive results were linked to previously described mutations mediating a low or moderate MIC increase. Virtually all drugs (99.0%) used in combination therapies that were inferred from genomic DST, were confirmed to be susceptible by pDST.

CONCLUSIONS

NGS-based genomic DST can reliably guide the design of effective MDR/RR-TB treatment regimens.

Evaluation of the Roche cobas MTB and MTB-RIF/INH Assays in Samples from Germany and Sierra Leone

The Roche cobas MTB and MTB-RIF/INH assays allow for detection of Mycobacterium tuberculosis complex (MTBC) nucleic acid and rifampicin (RIF) and isoniazid (INH) resistance-associated mutations in an automated, high-throughput workflow. In this study, we evaluated the performance of these assays, employing samples from settings of low and high tuberculosis (TB) burdens. A total of 325 frozen, leftover respiratory samples collected from treatment-naive patients with presumptive TB in Germany (n = 280) and presumptive RIF-resistant TB in Sierra Leone (n = 45) were used in this study. cobas MTB results for detection of MTBC DNA from N-acetyl-l-cysteine-sodium hydroxide (NALC-NaOH)-treated samples were compared to culture results. Predictions of RIF and INH resistance by the cobas MTB-RIF/INH assay were compared to a composite reference standard (phenotypic drug susceptibility testing and line probe assay). Whole-genome sequencing was used to resolve discordances. The overall sensitivity of cobas MTB for detection of MTBC DNA in culture-positive samples (n = 102) was 89.2% (95% confidence interval [CI], 81.7 to 93.9%). The specificity of cobas MTB was 98.6% (95% CI, 96.1 to 99.5%). Sensitivity and specificity for detection of RIF and INH resistance were 88.4% (95% CI, 75.5 to 94.9%) and 97.6% (95% CI, 87.4 to 99.6%) and 76.6% (95% CI, 62.8 to 86.4%) and 100.0% (95% CI, 90.8 to 100.0%), respectively. Discordant results for RIF and INH resistance were mainly due to uncommon mutations in samples from Sierra Leone that were not covered by the cobas MTB-RIF/INH assay. In conclusion, cobas MTB and MTB-RIF/INH assays provide accurate detection of MTBC DNA and resistance-associated mutations in respiratory samples. The influence of regional variations in the prevalence of resistance-conferring mutations requires further investigation.

Mycobacterium tuberculosis Complex Lineage 3 as Causative Agent of Pulmonary Tuberculosis, Eastern Sudan1

Pathogen-based factors associated with tuberculosis (TB) in eastern Sudan are not well defined. We investigated genetic diversity, drug resistance, and possible transmission clusters of Mycobacterium tuberculosis complex (MTBC) strains by using a genomic epidemiology approach. We collected 383 sputum specimens at 3 hospitals in 2014 and 2016 from patients with symptoms suggestive of TB; of these, 171 grew MTBC strains. Whole-genome sequencing could be performed on 166 MTBC strains; phylogenetic classification revealed that most (73.4%; n = 122) belonged to lineage 3 (L3). Genome-based cluster analysis showed that 76 strains (45.9%) were grouped into 29 molecular clusters, comprising 2–8 strains/patients. Of the strains investigated, 9.0% (15/166) were multidrug resistant (MDR); 10 MDR MTBC strains were linked to 1 large MDR transmission network. Our findings indicate that L3 strains are the main causative agent of TB in eastern Sudan; MDR TB is caused mainly by transmission of MDR L3 strains.

Multidrug- and Extensively Drug-Resistant Mycobacterium tuberculosis Beijing Clades, Ukraine, 2015

Multidrug-resistant (MDR) and extensively drug-resistant (XDR) tuberculosis (TB) is an emerging threat to TB control in Ukraine, a country with the third highest XDR TB burden globally. We used whole-genome sequencing of a convenience sample to identify bacterial genetic and patient-related factors associated with MDR/XDR TB in this country. MDR/XDR TB was associated with 3 distinct Mycobacterium tuberculosis complex lineage 2 (Beijing) clades, Europe/Russia W148 outbreak, Central Asia outbreak, and Ukraine outbreak, which comprised 68.9% of all MDR/XDR TB strains from southern Ukraine. MDR/XDR TB was also associated with previous treatment for TB and urban residence. The circulation of Beijing outbreak strains harboring broad drug resistance, coupled with constraints in drug supply and limited availability of phenotypic drug susceptibility testing, needs to be considered when new TB management strategies are implemented in Ukraine.

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.

Perspective for Precision Medicine for Tuberculosis

Tuberculosis is a bacterial infectious disease that is mainly transmitted from human to human via infectious aerosols. Currently, tuberculosis is the leading cause of death by an infectious disease world-wide. In the past decade, the number of patients affected by tuberculosis has increased by ~20 percent and the emergence of drug-resistant strains of Mycobacterium tuberculosis challenges the goal of elimination of tuberculosis in the near future. For the last 50 years, management of patients with tuberculosis has followed a standardized management approach. This standardization neglects the variation in human susceptibility to infection, immune response, the pharmacokinetics of drugs, and the individual duration of treatment needed to achieve relapse-free cure. Here we propose a package of precision medicine-guided therapies that has the prospect to drive clinical management decisions, based on both host immunity and M. tuberculosis strains genetics. Recently, important scientific discoveries and technological advances have been achieved that provide a perspective for individualized rather than standardized management of patients with tuberculosis. For the individual selection of best medicines and host-directed therapies, personalized drug dosing, and treatment durations, physicians treating patients with tuberculosis will be able to rely on these advances in systems biology and to apply them at the bedside.

Insertion and deletion evolution reflects antibiotics selection pressure in a Mycobacterium tuberculosis outbreak

In genome evolution, genetic variants are the source of diversity, which natural selection acts upon. Treatment of human tuberculosis (TB) induces a strong selection pressure for the emergence of antibiotic resistance-conferring variants in the infecting Mycobacterium tuberculosis (MTB) strains. MTB evolution in response to treatment has been intensively studied and mainly attributed to point substitutions. However, the frequency and contribution of insertions and deletions (indels) to MTB genome evolution remains poorly understood. Here, we analyzed a multi-drug resistant MTB outbreak for the presence of high-quality indels and substitutions. We find that indels are significantly enriched in genes conferring antibiotic resistance. Furthermore, we show that indels are inherited during the outbreak and follow a molecular clock with an evolutionary rate of 5.37e-9 indels/site/year, which is 23 times lower than the substitution rate. Inherited indels may co-occur with substitutions in genes along related biological pathways; examples are iron storage and resistance to second-line antibiotics. This suggests that epistatic interactions between indels and substitutions affect antibiotic resistance and compensatory evolution in MTB.

Evolutionary Approaches to Combat Antibiotic Resistance: Opportunities and Challenges for Precision Medicine

The rise of antimicrobial resistance (AMR) in bacterial pathogens is acknowledged by the WHO as a major global health crisis. It is estimated that in 2050 annually up to 10 million people will die from infections with drug resistant pathogens if no efficient countermeasures are implemented. Evolution of pathogens lies at the core of this crisis, which enables rapid adaptation to the selective pressures imposed by antimicrobial usage in both medical treatment and agriculture, consequently promoting the spread of resistance genes or alleles in bacterial populations. Approaches developed in the field of Evolutionary Medicine attempt to exploit evolutionary insight into these adaptive processes, with the aim to improve diagnostics and the sustainability of antimicrobial therapy. Here, we review the concept of evolutionary trade-offs in the development of AMR as well as new therapeutic approaches and their impact on host-microbiome-pathogen interactions. We further discuss the possible translation of evolution-informed treatments into clinical practice, considering both the rapid cure of the individual patients and the prevention of AMR.

Phylogenetically informative mutations in genes implicated in antibiotic resistance in Mycobacterium tuberculosis complex

BACKGROUND

A comprehensive understanding of the pre-existing genetic variation in genes associated with antibiotic resistance in the Mycobacterium tuberculosis complex (MTBC) is needed to accurately interpret whole-genome sequencing data for genotypic drug susceptibility testing (DST).

METHODS

We investigated mutations in 92 genes implicated in resistance to 21 anti-tuberculosis drugs using the genomes of 405 phylogenetically diverse MTBC strains. The role of phylogenetically informative mutations was assessed by routine phenotypic DST data for the first-line drugs isoniazid, rifampicin, ethambutol, and pyrazinamide from a separate collection of over 7000 clinical strains. Selected mutations/strains were further investigated by minimum inhibitory concentration (MIC) testing.

RESULTS

Out of 547 phylogenetically informative mutations identified, 138 were classified as not correlating with resistance to first-line drugs. MIC testing did not reveal a discernible impact of a Rv1979c deletion shared by M. africanum lineage 5 strains on resistance to clofazimine. Finally, we found molecular evidence that some MTBC subgroups may be hyper-susceptible to bedaquiline and clofazimine by different loss-of-function mutations affecting a drug efflux pump subunit (MmpL5).

CONCLUSIONS

Our findings underline that the genetic diversity in MTBC has to be studied more systematically to inform the design of clinical trials and to define sound epidemiologic cut-off values (ECOFFs) for new and repurposed anti-tuberculosis drugs. In that regard, our comprehensive variant catalogue provides a solid basis for the interpretation of mutations in genotypic as well as in phenotypic DST assays.

A landscape of genomic alterations at the root of a near-untreatable tuberculosis epidemic

BACKGROUND

Atypical Beijing genotype Mycobacterium tuberculosis strains are widespread in South Africa and have acquired resistance to up to 13 drugs on multiple occasions. It is puzzling that these strains have retained fitness and transmissibility despite the potential fitness cost associated with drug resistance mutations.

METHODS

We conducted Illumina sequencing of 211 Beijing genotype M. tuberculosis isolates to facilitate the detection of genomic features that may promote acquisition of drug resistance and restore fitness in highly resistant atypical Beijing forms. Phylogenetic and comparative genomic analysis was done to determine changes that are unique to the resistant strains that also transmit well. Minimum inhibitory concentration (MIC) determination for streptomycin and bedaquiline was done for a limited number of isolates to demonstrate a difference in MIC between isolates with and without certain variants.

RESULTS

Phylogenetic analysis confirmed that two clades of atypical Beijing strains have independently developed resistance to virtually all the potent drugs included in standard (pre-bedaquiline) drug-resistant TB treatment regimens. We show that undetected drug resistance in a progenitor strain was likely instrumental in this resistance acquisition. In this cohort, ethionamide (ethA A381P) resistance would be missed in first-line drug-susceptible isolates, and streptomycin (gidB L79S) resistance may be missed due to an MIC close to the critical concentration. Subsequent inadequate treatment historically led to amplification of resistance and facilitated spread of the strains. Bedaquiline resistance was found in a small number of isolates, despite lack of exposure to the drug. The highly resistant clades also carry inhA promoter mutations, which arose after ethA and katG mutations. In these isolates, inhA promoter mutations do not alter drug resistance, suggesting a possible alternative role.

CONCLUSION

The presence of the ethA mutation in otherwise susceptible isolates from ethionamide-naïve patients demonstrates that known exposure is not an adequate indicator of drug susceptibility. Similarly, it is demonstrated that bedaquiline resistance can occur without exposure to the drug. Inappropriate treatment regimens, due to missed resistance, leads to amplification of resistance, and transmission. We put these results into the context of current WHO treatment regimens, underscoring the risks of treatment without knowledge of the full drug resistance profile.

Antibiotic resistance prediction for Mycobacterium tuberculosis from genome sequence data with Mykrobe

Two billion people are infected with Mycobacterium tuberculosis, leading to 10 million new cases of active tuberculosis and 1.5 million deaths annually. Universal access to drug susceptibility testing (DST) has become a World Health Organization priority. We previously developed a software tool, Mykrobe predictor, which provided offline species identification and drug resistance predictions for M. tuberculosis from whole genome sequencing (WGS) data. Performance was insufficient to support the use of WGS as an alternative to conventional phenotype-based DST, due to mutation catalogue limitations.  Here we present a new tool, Mykrobe, which provides the same functionality based on a new software implementation. Improvements include i) an updated mutation catalogue giving greater sensitivity to detect pyrazinamide resistance, ii) support for user-defined resistance catalogues, iii) improved identification of non-tuberculous mycobacterial species, and iv) an updated statistical model for Oxford Nanopore Technologies sequencing data. Mykrobe is released under MIT license at https://github.com/mykrobe-tools/mykrobe. We incorporate mutation catalogues from the CRyPTIC consortium et al. (2018) and from Walker et al. (2015), and make improvements based on performance on an initial set of 3206 and an independent set of 5845 M. tuberculosis Illumina sequences. To give estimates of error rates, we use a prospectively collected dataset of 4362 M. tuberculosis isolates. Using culture based DST as the reference, we estimate Mykrobe to be 100%, 95%, 82%, 99% sensitive and 99%, 100%, 99%, 99% specific for rifampicin, isoniazid, pyrazinamide and ethambutol resistance prediction respectively. We benchmark against four other tools on 10207 (=5845+4362) samples, and also show that Mykrobe gives concordant results with nanopore data.  We measure the ability of Mykrobe-based DST to guide personalized therapeutic regimen design in the context of complex drug susceptibility profiles, showing 94% concordance of implied regimen with that driven by phenotypic DST, higher than all other benchmarked tools.

Accuracy of whole-genome sequencing to determine recent tuberculosis transmission: an 11-year population-based study in Hamburg, Germany

Controlling human-to-human tuberculosis (TB) transmission is key for achieving the targets of the End TB Strategy set by the World Health Organization (WHO) [1, 2]. Stopping TB transmission, in large cities especially, is a challenging top priority worldwide [3]. Metropolitan areas have higher TB case notification rates than the rest of a country, as they concentrate high-risk groups, such as homeless people, drug users and migrants often from (other) high TB incidence settings. Opportunities for transmission are amplified by population density and complex social interactions, regularly leading to large, temporally extended transmission networks [3]. Targeted interventions to interrupt transmission require the combination of effective genotyping of TB strains with enhanced epidemiological investigation. While classic IS6110 DNA fingerprinting and 24-locus MIRU–VNTR (mycobacterial interspersed repetitive units–variable number of tandem repeats) typing provide standardised and easily computable typing results with an online nomenclature system, several studies have now demonstrated that whole-genome sequencing (WGS) has a superior discriminatory power, allowing for an unparalleled resolution of outbreak strains [4–10]. However, predictivity of WGS for detecting transmission in metropolitan areas has not yet been quantified versus most deterministic references, i.e. tangible epidemiological links identified by ad hoc investigation, at extended time and population scales.

WGS typing with a five-SNP c:d:ut-off delineates recent transmission chains with highest accuracy and also provides high-resolution resistance patterns, thus enabling direct clinical benefitshttp://bit.ly/2Pk37Wo

Population structure of drug-resistant Mycobacterium tuberculosis in Central Asia

Background

Drug-resistant tuberculosis (TB) is a major public health concern threathing the success of TB control efforts, and this is particularily problematic in Central Asia. Here, we present the first analysis of the population structure of Mycobacterium tuberculosis complex isolates in the Central Asian republics Uzbekistan, Tajikistan, and Kyrgyzstan.

Methods

The study set consisted of 607 isolates with 235 from Uzbekistan, 206 from Tajikistan, and 166 from Kyrgyzstan. 24-loci MIRU-VNTR (Mycobacterial Interspersed Repetitive Units - Variable Number of Tandem Repeats) typing and spoligotyping were combined for genotyping. In addition, phenotypic drug suceptibility was performed.

Results

The population structure mainly comprises strains of the Beijing lineage (411/607). 349 of the 411 Beijing isolates formed clusters, compared to only 33 of the 196 isolates from other clades. Beijing 94–32 (n = 145) and 100–32 (n = 70) formed the largest clusters. Beijing isolates were more frequently multidrug-resistant, pre-extensively resistant (pre-XDR)- or XDR-TB than other genotypes.

Conclusions

Beijing clusters 94–32 and 100–32 are the dominant MTB genotypes in Central Asia. The relative size of 100–32 compared to previous studies in Kazakhstan and its unequal geographic distribution support the hypothesis of its more recent emergence in Central Asia. The data also demonstrate that clonal spread of resistant TB strains, particularly of the Beijing lineage, is a root of the so far uncontroled MDR-TB epidemic in Central Asia.

Longitudinal Outbreak of Multidrug-Resistant Tuberculosis in a Hospital Setting, Serbia

A retrospective population-based molecular epidemiologic study of multidrug-resistant Mycobacterium tuberculosis complex strains in Serbia (2008–2014) revealed an outbreak of TUR genotype strains in a psychiatric hospital starting around 1990. Drug unavailability, poor infection control, and schizophrenia likely fueled acquisition of additional resistance and bacterial fitness–related mutations over 2 decades.

The role of IS6110 in micro- and macroevolution of Mycobacterium tuberculosis lineage 2

The insertion sequence 6110 (IS6110) is the most studied transposable element in the Mycobacterium tuberculosis complex species. The element plays a significant role in genome plasticity of this important human pathogen, but still many causes and consequences of its transposition have not been fully studied. Here, we analyzed insertion sites for 902 Mycobacterium tuberculosis lineage 2 strains using whole-genome sequencing data. In total, 17,972 insertions were found, corresponding to 827 independent positions in the genome of the reference strain H37Rv. To trace the history of IS6110 expansion since proto-Beijing strains up to modern sublineages, we looked at the distribution of IS6110 across the genome-wide SNP-based phylogenetic tree. This analysis demonstrated a stepwise transposition of IS6110 that occurs by «copy-and-paste» mechanism. Additionally, we detected evolutionary-scale and sublineage-specific integration sites, which can be used for typing and for understanding the reasons for the success of the lineage. A significant part of such insertions affected the genes that are essential for the pathogen. Finally, we identified and confirmed deletions that occurred between differently oriented elements, which is uncommon for this family of insertion elements and appears to be another mechanism of genome variability.

Emergence of Low-level Delamanid and Bedaquiline Resistance During Extremely Drug-resistant Tuberculosis Treatment

Two new drugs, delamanid and bedaquiline, have recently been approved for treatment of multidrug-resistant and extensively drug-resistant (XDR) tuberculosis. Here, we report a case of clofazimine, bedaquiline, and low-level delamanid resistances acquired during treatment of a patient with XDR tuberculosis.

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.

Drug-resistance profiling and transmission dynamics of multidrug-resistant Mycobacterium tuberculosis in Saudi Arabia revealed by whole genome sequencing

Background

In Saudi Arabia, cross-border transmission of multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains might be particularly fostered by high immigration rates. Herein, we aimed to elucidate the transmission dynamics of MDR-MTBC strains and reveal a detailed prediction of all resistance-conferring mutations for the first- and second-line drugs.

Methods

We investigated all MDR-MTBC strains collected between 2015 and 2017 from provincial mycobacteria referral laboratories and compared demographic and clinical parameters to a cohort of non-MDR-TB patients using a whole genome sequencing approach. Clusters were defined based on a maximum strain-to-strain genetic distance of five single-nucleotide polymorphisms (SNPs) as surrogate marker for recent transmission, and then investigated molecular drug-resistance markers (37 genes).

Results

Forty-eight (67.6%) MDR-MTBC strains were grouped in 14 different clusters, ranging in size from two to six strains; 22.5% (16/71) of all MDR-MTBC isolates were predicted to be fully resistant to all five first-line drugs, and five strains (7.0%) exhibited fluoroquinolone resistance. Moreover, we revealed the presence of 12 compensatory mutations as well as 26 non-synonymous SNPs in the rpoC gene and non-hotspot region in rpoB, respectively.

Conclusion

Optimized TB molecular surveillance, diagnosis, and patient management are urgently needed to contain MDR-MTBC transmission and prevent the development of additional drug resistance.

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.

Compensatory evolution drives multidrug-resistant tuberculosis in Central Asia

Bacterial factors favoring the unprecedented multidrug-resistant tuberculosis (MDR-TB) epidemic in the former Soviet Union remain unclear. We utilized whole genome sequencing and Bayesian statistics to analyze the evolutionary history, temporal emergence of resistance and transmission networks of MDR Mycobacterium tuberculosis complex isolates from Karakalpakstan, Uzbekistan (2001-2006). One clade (termed Central Asian outbreak, CAO) dating back to 1974 (95% HPD 1969-1982) subsequently acquired resistance mediating mutations to eight anti-TB drugs. Introduction of standardized WHO-endorsed directly observed treatment, short-course in Karakalpakstan in 1998 likely selected for CAO-strains, comprising 75% of sampled MDR-TB isolates in 2005/2006. CAO-isolates were also identified in a published cohort from Russia (2008-2010). Similarly, the presence of mutations supposed to compensate bacterial fitness deficits was associated with transmission success and higher drug resistance rates. The genetic make-up of these MDR-strains threatens the success of both empirical and standardized MDR-TB therapies, including the newly WHO-endorsed short MDR-TB regimen in Uzbekistan.