Back-to-Africa introductions of Mycobacterium tuberculosis as the main cause of tuberculosis in Dar es Salaam, Tanzania.

In settings with high tuberculosis (TB) endemicity, distinct genotypes of the Mycobacterium tuberculosis complex (MTBC) often differ in prevalence. However, the factors leading to these differences remain poorly understood. Here we studied the MTBC population in Dar es Salaam, Tanzania over a six-year period, using 1,082 unique patient-derived MTBC whole-genome sequences (WGS) and associated clinical data. We show that the TB epidemic in Dar es Salaam is dominated by multiple MTBC genotypes introduced to Tanzania from different parts of the world during the last 300 years. The most common MTBC genotypes deriving from these introductions exhibited differences in transmission rates and in the duration of the infectious period, but little differences in overall fitness, as measured by the effective reproductive number. Moreover, measures of disease severity and bacterial load indicated no differences in virulence between these genotypes during active TB. Instead, the combination of an early introduction and a high transmission rate accounted for the high prevalence of L3.1.1, the most dominant MTBC genotype in this setting. Yet, a longer co-existence with the host population did not always result in a higher transmission rate, suggesting that distinct life-history traits have evolved in the different MTBC genotypes. Taken together, our results point to bacterial factors as important determinants of the TB epidemic in Dar es Salaam.

Bedaquiline Resistance after Effective Treatment of Multidrug-Resistant Tuberculosis, Namibia.

Bedaquiline is currently a key drug for treating multidrug-resistant or rifampin-resistant tuberculosis. We report and discuss the unusual development of resistance to bedaquiline in a teenager in Namibia, despite an optimal background regimen and adherence. The report highlights the risk for bedaquiline resistance development and the need for rapid drug-resistance testing.

Challenging the gold standard: the limitations of molecular assays for detection of Mycobacterium tuberculosis heteroresistance.

OBJECTIVES: Heteroresistant infections are defined as infections in which a mixture of drug-resistant and drug-susceptible populations are present. In Mycobacterium tuberculosis (M. tb), heteroresistance poses a challenge in diagnosis and has been linked with poor treatment outcomes. We compared the analytical sensitivity of molecular methods, such as GeneXpert and whole genome sequencing (WGS) in detecting heteroresistance when compared with the 'gold standard' phenotypic assay: the agar proportion method (APM). METHODS: Using two rounds of proficiency surveys with defined monoresistant BCG strains and mixtures of susceptible/resistant M. tb, we determined the limit of detection (LOD) of known resistance associated mutations. RESULTS: The LOD for rifampin-R (RIF-R) detection was 1% using APM, 60% using GeneXpert MTB/RIF, 10% using GeneXpert MTB/RIF Ultra and 10% using WGS. While WGS could detect mutations beyond those associated with RIF resistance, the LOD for these other mutations was also 10%. Additionally, we observed instances where laboratories did not report resistance in the majority population, yet the mutations were present in the raw sequence data. CONCLUSION: The gold standard APM detects minority resistant populations at a lower proportion than molecular tests. Mycobacterium bovis BCG strains with defined resistance and extracted DNA from M. tb provided concordant results and can serve in quality control of laboratories offering molecular testing for resistance. Further research is required to determine whether the higher LOD of molecular tests is associated with negative treatment outcomes.

First detection of a Mycobacterium tuberculosis XDR clinical isolate harbouring an RpoB I491F mutation in a Ukrainian patient treated in Germany, October 2023.

This report documents the case of a Ukrainian patient infected with an extensively drug-resistant (XDR) lineage 2 Mycobacterium tuberculosis strain harbouring the rifampicin resistance mutation RpoB I491F. This mutation is not detected by routine molecular WHO-recommended rapid diagnostics, complicating the detection and treatment of these strains. The occurrence of such mutations underscores the need for enhanced diagnostic techniques and tailored treatment regimens, especially in eastern Europe where lineage 2 strains and XDR-tuberculosis are prevalent.

Tuberculosis Variant with Rifampin Resistance Undetectable by Xpert MTB/RIF, Botswana.

GeneXpert MTB/RIF, a tool widely used for diagnosing tuberculosis, has limitations for detecting rifampin resistance in certain variants. We report transmission of a pre-extensively drug-resistant variant in Botswana that went undetected by GeneXpert. The public health impact of misdiagnosis emphasizes the need for comprehensive molecular testing to identify resistance and guide treatment.

Use of High-Resolution Geospatial and Genomic Data to Characterize Recent Tuberculosis Transmission, Botswana.

Combining genomic and geospatial data can be useful for understanding Mycobacterium tuberculosis transmission in high-burden tuberculosis (TB) settings. We performed whole-genome sequencing on M. tuberculosis DNA extracted from sputum cultures from a population-based TB study conducted in Gaborone, Botswana, during 2012-2016. We determined spatial distribution of cases on the basis of shared genotypes among isolates. We considered clusters of isolates with ≤5 single-nucleotide polymorphisms identified by whole-genome sequencing to indicate recent transmission and clusters of ≥10 persons to be outbreaks. We obtained both molecular and geospatial data for 946/1,449 (65%) participants with culture-confirmed TB; 62 persons belonged to 5 outbreaks of 10-19 persons each. We detected geospatial clustering in just 2 of those 5 outbreaks, suggesting heterogeneous spatial patterns. Our findings indicate that targeted interventions applied in smaller geographic areas of high-burden TB identified using integrated genomic and geospatial data might help interrupt TB transmission during outbreaks.

Characterization of Two Novel Inhibitors of the Mycobacterium tuberculosis Cytochrome bc1 Complex.

Drug-resistant tuberculosis is a global health care threat calling for novel effective treatment options. Here, we report on two novel cytochrome bc1 inhibitors (MJ-22 and B6) targeting the Mycobacterium tuberculosis respiratory chain with excellent intracellular activities in human macrophages. Both hit compounds revealed very low mutation frequencies and distinct cross-resistance patterns with other advanced cytochrome bc1 inhibitors.

Using genetic data to identify transmission risk factors: Statistical assessment and application to tuberculosis transmission.

Identifying host factors that influence infectious disease transmission is an important step toward developing interventions to reduce disease incidence. Recent advances in methods for reconstructing infectious disease transmission events using pathogen genomic and epidemiological data open the door for investigation of host factors that affect onward transmission. While most transmission reconstruction methods are designed to work with densely sampled outbreaks, these methods are making their way into surveillance studies, where the fraction of sampled cases with sequenced pathogens could be relatively low. Surveillance studies that use transmission event reconstruction then use the reconstructed events as response variables (i.e., infection source status of each sampled case) and use host characteristics as predictors (e.g., presence of HIV infection) in regression models. We use simulations to study estimation of the effect of a host factor on probability of being an infection source via this multi-step inferential procedure. Using TransPhylo-a widely-used method for Bayesian estimation of infectious disease transmission events-and logistic regression, we find that low sensitivity of identifying infection sources leads to dilution of the signal, biasing logistic regression coefficients toward zero. We show that increasing the proportion of sampled cases improves sensitivity and some, but not all properties of the logistic regression inference. Application of these approaches to real world data from a population-based TB study in Botswana fails to detect an association between HIV infection and probability of being a TB infection source. We conclude that application of a pipeline, where one first uses TransPhylo and sparsely sampled surveillance data to infer transmission events and then estimates effects of host characteristics on probabilities of these events, should be accompanied by a realistic simulation study to better understand biases stemming from imprecise transmission event inference.

Resistance patterns among drug-resistant tuberculosis patients and trends-over-time analysis of national surveillance data in Gabon, Central Africa.

OBJECTIVE: Routinely generated surveillance data are important for monitoring the effectiveness of MDR-TB control strategies. Incidence of rifampicin-resistant tuberculosis (RR-TB) is a key indicator for monitoring MDR-TB. METHODS: In a longitudinal nationwide retrospective study, 8 years (2014-2021) of sputum samples from presumptively drug-resistant tuberculosis patients from all regions of Gabon were referred to the national tuberculosis reference laboratory. Samples were analysed using GeneXpert MTB/RIF and Genotype MTBDRsl version 2/Line Probe Assay. RESULTS: Of 3057 sputum samples from presumptive tuberculosis patients, both from local hospital and from referral patients, 334 were RR-TB. The median patient age was 33 years (interquartile range 26-43); one third was newly diagnosed drug-resistant tuberculosis patients; one-third was HIV-positive. The proportion of men with RR-TB was significantly higher than that of women (55% vs 45%; p < 0.0001). Patients aged 25-35 years were most affected (32%; 108/334). The cumulative incidence of RR-TB was 17 (95% CI 15-19)/100,000 population over 8 years. The highest incidences were observed in 2020 and 2021. A total of 281 samples were analysed for second-line drug resistance. The proportions of study participants with MDR-TB, pre-XDR-TB and XDR-TB were 90.7% (255/281), 9% (25/281) and 0.3% (1/281), respectively. The most-common mutations in fluoroquinolones resistance isolates was gyrA double mutation gyrA MUT3B and MUT3C (23%; 4/17). Most (64%; 6/8) second-line injectable drugs resistance isolates were characterised by missing both rrs WT2 and MUT2 banding. CONCLUSION: The increasing incidence of MDR-TB infection in Gabon is alarming. It is highest in the 25-35 years age category. The incidence of MDR-TB infection in treatment-naïve patients calls for case finding and contact tracing strategy improvement.

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.

Ancient and recent differences in the intrinsic susceptibility of Mycobacterium tuberculosis complex to pretomanid.

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.

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.

Loss of Mobile Genomic Islands in Metal-Resistant, Hydrogen-Oxidizing Cupriavidus metallidurans.

The genome of the metal-resistant, hydrogen-oxidizing bacterium Cupriavidus metallidurans strain CH34 contains horizontally acquired plasmids and genomic islands. Metal-resistance determinants on the two plasmids may exert genetic dominance over other related determinants. To investigate whether these recessive determinants can be activated in the absence of the dominant ones, the transcriptome of the highly zinc-sensitive deletion mutant Δe4 (ΔcadA ΔzntA ΔdmeF ΔfieF) of the plasmid-free parent AE104 was characterized using gene arrays. As a consequence of some unexpected results, close examination by PCR and genomic resequencing of strains CH34, AE104, Δe4, and others revealed that the genomic islands CMGI2, 3, 4, D, and E, but no other islands or recessive determinants, were deleted in some of these strains. Provided that wild-type CH34 was kept under alternating zinc and nickel selection pressure, no comparable deletions occurred. All current data suggest that genes were actually deleted and were not, as surmised previously, silenced in the respective strain. As a consequence, a cured database was compiled from the newly generated and previously published gene array data. An analysis of data from this database indicated that some genes of recessive, no longer needed determinants were nevertheless expressed and upregulated. Their products may interact with those of the dominant determinants to mediate a mosaic phenotype. The ability to contribute to such a mosaic phenotype may prevent deletion of the recessive determinant. The data suggest that the bacterium actively modifies its genome to deal with metal stress and at the same time ensures metal homeostasis. IMPORTANCE In their natural environment, bacteria continually acquire genes by horizontal gene transfer, and newly acquired determinants may become dominant over related ones already present in the host genome. When a bacterium is taken into laboratory culture, it is isolated from the horizontal gene transfer network. It can no longer gain genes but instead may lose them. This phenomenon was indeed observed in Cupriavidus metallidurans for the loss key metal resistance determinants when no selection pressure was kept continuously. However, some recessive metal resistance determinants were maintained in the genome. It is proposed that they might contribute some accessory genes to related dominant resistance determinants, for instance periplasmic metal-binding proteins or two-component regulatory systems. Alternatively, they may remain in the genome only because their DNA serves as a scaffold for the nucleoid. Using C. metallidurans as an example, this study sheds light on the fate and function of horizontally acquired genes in bacteria.

A multi-phenotype genome-wide association study of clades causing tuberculosis in a Ghanaian- and South African cohort.

Despite decades of research and advancements in diagnostics and treatment, tuberculosis remains a major public health concern. New computational methods are needed to interrogate the intersection of host- and bacterial genomes. Paired host genotype datum and infecting bacterial isolate information were analysed for associations using a multinomial logistic regression framework implemented in SNPTest. A cohort of 853 admixed South African participants and a Ghanaian cohort of 1359 participants were included. Two directly genotyped variants, namely rs529920 and rs41472447, were identified in the Ghanaian cohort as being statistically significantly associated with risk for infection with strains of different members of the MTBC. Thus, a multinomial logistic regression using paired host-pathogen data may prove valuable for investigating the complex relationships driving infectious disease.

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 determined whether 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 1 January 2015 and 30 April 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 minimal inhibitory drug concentrations (MICs) 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 the presence of polymorphisms with yet unknown implications for drug resistance. Importantly, 549 of 561 (97.9%) predictions of drug susceptibility were phenotypically confirmed in MGIT, and 27 of 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 phenotypic DST. CONCLUSIONS: NGS-based genomic DST can reliably guide the design of effective MDR/RR-TB treatment regimens.

Secretome characterization of clinical isolates from the Mycobacterium abscessus complex provides insight into antigenic differences.

BACKGROUND: Mycobacterium abscessus (MAB) is a widely disseminated pathogenic non-tuberculous mycobacterium (NTM). Like with the M. tuberculosis complex (MTBC), excreted / secreted (ES) proteins play an essential role for its virulence and survival inside the host. Here, we used a robust bioinformatics pipeline to predict the secretome of the M. abscessus ATCC 19977 reference strain and 15 clinical isolates belonging to all three MAB subspecies, M. abscessus subsp. abscessus, M. abscessus subsp. bolletii, and M. abscessus subsp. massiliense. RESULTS: We found that ~ 18% of the proteins encoded in the MAB genomes were predicted as secreted and that the three MAB subspecies shared > 85% of the predicted secretomes. MAB isolates with a rough (R) colony morphotype showed larger predicted secretomes than isolates with a smooth (S) morphotype. Additionally, proteins exclusive to the secretomes of MAB R variants had higher antigenic densities than those exclusive to S variants, independent of the subspecies. For all investigated isolates, ES proteins had a significantly higher antigenic density than non-ES proteins. We identified 337 MAB ES proteins with homologues in previously investigated M. tuberculosis secretomes. Among these, 222 have previous experimental support of secretion, and some proteins showed homology with protein drug targets reported in the DrugBank database. The predicted MAB secretomes showed a higher abundance of proteins related to quorum-sensing and Mce domains as compared to MTBC indicating the importance of these pathways for MAB pathogenicity and virulence. Comparison of the predicted secretome of M. abscessus ATCC 19977 with the list of essential genes revealed that 99 secreted proteins corresponded to essential proteins required for in vitro growth. CONCLUSIONS: This study represents the first systematic prediction and in silico characterization of the MAB secretome. Our study demonstrates that bioinformatics strategies can help to broadly explore mycobacterial secretomes including those of clinical isolates and to tailor subsequent, complex and time-consuming experimental approaches accordingly. This approach can support systematic investigation exploring candidate proteins for new vaccines and diagnostic markers to distinguish between colonization and infection. All predicted secretomes were deposited in the Secret-AAR web-server ( http://microbiomics.ibt.unam.mx/tools/aar/index.php ).

Role of Epistasis in Amikacin, Kanamycin, Bedaquiline, and Clofazimine Resistance in Mycobacterium tuberculosis Complex.

Antibiotic resistance among bacterial pathogens poses a major global health threat. Mycobacterium tuberculosis complex (MTBC) is estimated to have the highest resistance rates of any pathogen globally. Given the low growth rate and the need for a biosafety level 3 laboratory, the only realistic avenue to scale up drug susceptibility testing (DST) for this pathogen is to rely on genotypic techniques. This raises the fundamental question of whether a mutation is a reliable surrogate for phenotypic resistance or whether the presence of a second mutation can completely counteract its effect, resulting in major diagnostic errors (i.e., systematic false resistance results). To date, such epistatic interactions have only been reported for streptomycin that is now rarely used. By analyzing more than 31,000 MTBC genomes, we demonstrated that the eis C-14T promoter mutation, which is interrogated by several genotypic DST assays endorsed by the World Health Organization, cannot confer resistance to amikacin and kanamycin if it coincides with loss-of-function (LoF) mutations in the coding region of eis. To our knowledge, this represents the first definitive example of antibiotic reversion in MTBC. Moreover, we raise the possibility that mmpR (Rv0678) mutations are not valid markers of resistance to bedaquiline and clofazimine if these coincide with an LoF mutation in the efflux pump encoded by mmpS5 (Rv0677c) and mmpL5 (Rv0676c).

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.

Use of a whole genome sequencing-based approach for Mycobacterium tuberculosis surveillance in Europe in 2017-2019: an ECDC pilot study.

Whole genome sequencing (WGS) can be used for molecular typing and characterisation of Mycobacterium tuberculosis complex (MTBC) strains. We evaluated the systematic use of a WGS-based approach for MTBC surveillance involving all European Union/European Economic Area (EU/EEA) countries and highlight the challenges and lessons learnt to be considered for the future development of a WGS-based surveillance system.WGS and epidemiological data of patients with rifampicin-resistant (RR) and multidrug-resistant (MDR) tuberculosis (TB) were collected from EU/EEA countries between January 2017 and December 2019. WGS-based genetic relatedness analysis was performed using a standardised approach including both core genome multilocus sequence typing (cgMLST) and single nucleotide polymorphism (SNP)-based calculation of distances on all WGS data that fulfilled minimum quality criteria to ensure data comparability.A total of 2218 RR/MDR-MTBC isolates were collected from 25 countries. Among these, 56 cross-border clusters with increased likelihood of recent transmission (≤5 SNPs distance) comprising 316 RR/MDR-MTBC isolates were identified. The cross-border clusters included between two and 30 resistant isolates from two to six countries, demonstrating different RR/MDR-TB transmission patterns in Western and Eastern EU countries.This pilot study shows that a WGS-based surveillance system is not only feasible but can efficiently elucidate the dynamics of in-country and cross-border RR/MDR-TB transmission across EU/EEA countries. Lessons learnt from this study highlight that the establishment of an EU/EEA centralised WGS-based surveillance system for TB will require strengthening of national integrated systems performing prospective WGS surveillance and the development of clear procedures to facilitate international collaboration for the investigation of cross-border clusters.

Deep amplicon sequencing for culture-free prediction of susceptibility or resistance to 13 anti-tuberculous drugs.

Conventional molecular tests for detecting Mycobacterium tuberculosis complex (MTBC) drug resistance on clinical samples cover a limited set of mutations. Whole-genome sequencing (WGS) typically requires culture.Here, we evaluated the Deeplex Myc-TB targeted deep-sequencing assay for prediction of resistance to 13 anti-tuberculous drugs/drug classes, directly applicable on sputum.With MTBC DNA tests, the limit of detection was 100-1000 genome copies for fixed resistance mutations. Deeplex Myc-TB captured in silico 97.1-99.3% of resistance phenotypes correctly predicted by WGS from 3651 MTBC genomes. On 429 isolates, the assay predicted 92.2% of 2369 first- and second-line phenotypes, with a sensitivity of 95.3% and a specificity of 97.4%. 56 out of 69 (81.2%) residual discrepancies with phenotypic results involved pyrazinamide, ethambutol and ethionamide, and low-level rifampicin or isoniazid resistance mutations, all notoriously prone to phenotypic testing variability. Only two out of 91 (2.2%) resistance phenotypes undetected by Deeplex Myc-TB had known resistance-associated mutations by WGS analysis outside Deeplex Myc-TB targets. Phenotype predictions from Deeplex Myc-TB analysis directly on 109 sputa from a Djibouti survey matched those of MTBSeq/PhyResSE/Mykrobe, fed with WGS data from subsequent cultures, with a sensitivity of 93.5/98.5/93.1% and a specificity of 98.5/97.2/95.3%, respectively. Most residual discordances involved gene deletions/indels and 3-12% heteroresistant calls undetected by WGS analysis or natural pyrazinamide resistance of globally rare "Mycobacterium canettii" strains then unreported by Deeplex Myc-TB. On 1494 arduous sputa from a Democratic Republic of the Congo survey, 14 902 out of 19 422 (76.7%) possible susceptible or resistance phenotypes could be predicted culture-free.Deeplex Myc-TB may enable fast, tailored tuberculosis treatment.