Prison as a driver of recent transmissions of multidrug-resistant tuberculosis in Callao, Peru: a cross-sectional study.

Background: We sought to identify resistance patterns and key drivers of recent multidrug-resistant tuberculosis (MDR-TB) transmission in a TB-prevalent area in Peru. Methods: Cross-sectional study including MDR Mycobacterium tuberculosis complex (Mtbc) strains identified in Callao-Peru between April 2017 and February 2019. Mtbc DNA was extracted for whole genome sequencing which was used for phylogenetic inference, clustering, and resistance mutation analyses. Clusters indicative of recent transmission were defined based on a strain-to-strain distance of ≤5 (D5) single nucleotide polymorphisms (SNPs). Epidemiologic factors linked to MDR-TB clustering were analyzed using Poisson regression. Findings: 171 unique MDR-Mtbc strains were included; 22 (13%) had additional fluoroquinolone resistance and were classified as pre-XDR. Six strains (3.5%) harboured bedaquiline (BDQ) resistance mutations and were classified as MDR + BDQ. 158 (92%) Mtbc strains belonged to lineage 4 and 13 (8%) to lineage 2. Using a cluster threshold of ≤5 SNPs, 98 (57%) strains were grouped in one of the 17 D5 clusters indicative of recent transmission, ranging in size from 2 to the largest cluster formed by 53 4.3.3 strains (group_1). Lineage 4.3.3 strains showed the overall highest cluster rate (43%). In multivariate analyses, current or previous imprisonment was independently associated with being part of any MDR-TB transmission clusters (adjusted prevalence ratio [aPR], 1.45; 95% CI, 1.09-1.92). Interpretation: Pre-XDR-TB emerged in more than 10% of the MDR-TB strains investigated. Transmission of 4.3.3 Mtbc strains especially of the dominant group_1 clone is a major driver of the MDR-TB epidemic in Callao. Current or previous imprisonment was linked to recent MDR-TB transmissions, indicating an important role of prisons in driving the MDR-TB epidemic. Funding: This work was supported in part by the ERANet-LAC Network of the European Union, Latin America and the Caribbean Countries on Joint Innovation and Research Activities, and FONDECYT. Additional support was received from Leibniz Science Campus Evolutionary Medicine of the Lung, the Deutsche Forschungsgemeinschaft (German Research Foundation, under Germany's Excellence Strategy-EXC 2167 Precision Medicine in Inflammation), and the Research Training Group 2501 TransEvo.

Molecular determinants of multidrug-resistant tuberculosis in Sierra Leone.

Multidrug-resistant tuberculosis (MDR-TB) management has become a serious global health challenge. Understanding its epidemic determinants on the regional level is crucial for developing effective control measures. We used whole genome sequencing data of 238 of Mycobacterium tuberculosis complex (MTBC) strains to determine drug resistance profiles, phylogeny, and transmission dynamics of MDR/rifampicin-resistant (RR) MTBC strains from Sierra Leone. Forty-two strains were classified as RR, 196 as MDR, 5 were resistant to bedaquiline (BDQ) and clofazimine (CFZ), but none was found to be resistant to fluoroquinolones. Sixty-one (26%) strains were resistant to all first-line drugs, three of which had additional resistance to BDQ/CFZ. The strains were classified into six major MTBC lineages (L), with strains of L4 being the most prevalent, 62% (n = 147), followed by L6 (Mycobacterium africanum) strains, (21%, n = 50). The overall clustering rate (using ≤d12 single-nucleotide polymorphism threshold) was 44%, stratified into 31 clusters ranging from 2 to 16 strains. The largest cluster (n = 16) was formed by sublineage 2.2.1 Beijing Ancestral 3 strains, which developed MDR several times. Meanwhile, 10 of the L6 strains had a primary MDR transmission. We observed a high diversity of drug resistance mutations, including borderline resistance mutations to isoniazid and rifampicin, and mutations were not detected by commercial assays. In conclusion, one in five strains investigated was resistant to all first-line drugs, three of which had evidence of BDQ/CFZ resistance. Implementation of interventions such as rapid diagnostics that prevent further resistance development and stop MDR-TB transmission chains in the country is urgently needed. IMPORTANCE: A substantial proportion of MDR-TB strains in Sierra Leone were resistant against all first line drugs; however this makes the all-oral-six-month BPaLM regimen or other 6-9 months all oral regimens still viable, mainly because there was no FQ resistance.Resistance to BDQ was detected, as well as RR, due to mutations outside of the hotspot region. While the prevalence of those resistances was low, it is still cause for concern and needs to be closely monitored.

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.

Discovery of dual-active ethionamide boosters inhibiting the Mycobacterium tuberculosis ESX-1 secretion system.

Drug-resistant Mycobacterium tuberculosis (Mtb) remains a major public health concern requiring complementary approaches to standard anti-tuberculous regimens. Anti-virulence molecules or compounds that enhance the activity of antimicrobial prodrugs are promising alternatives to conventional antibiotics. Exploiting host cell-based drug discovery, we identified an oxadiazole compound (S3) that blocks the ESX-1 secretion system, a major virulence factor of Mtb. S3-treated mycobacteria showed impaired intracellular growth and a reduced ability to lyse macrophages. RNA sequencing experiments of drug-exposed bacteria revealed strong upregulation of a distinct set of genes including ethA, encoding a monooxygenase activating the anti-tuberculous prodrug ethionamide. Accordingly, we found a strong ethionamide boosting effect in S3-treated Mtb. Extensive structure-activity relationship experiments revealed that anti-virulence and ethionamide-boosting activity can be uncoupled by chemical modification of the primary hit molecule. To conclude, this series of dual-active oxadiazole compounds targets Mtb via two distinct mechanisms of action.

Emergence of bedaquiline-resistant tuberculosis and of multidrug-resistant and extensively drug-resistant Mycobacterium tuberculosis strains with rpoB Ile491Phe mutation not detected by Xpert MTB/RIF in Mozambique: a retrospective observational study.

BACKGROUND: In 2021, an estimated 4800 people developed rifampicin-resistant tuberculosis in Mozambique, 75% of which went undiagnosed. Detailed molecular data on rifampicin-resistant and multidrug-resistant (MDR) tuberculosis are not available. Here, we aimed at gaining precise data on the determinants of rifampicin-resistant and MDR tuberculosis in Mozambique. METHODS: In this retrospective observational study, we performed whole-genome sequencing of 704 rifampicin-resistant Mycobacterium tuberculosis complex (Mtbc) strains submitted to the National Tuberculosis Reference Laboratory (NTRL) in Maputo, Mozambique, between 2015 and 2021. Phylogenetic strain classification, genomic resistance prediction, and cluster analysis were performed. FINDINGS: Between Jan 1, 2015, and July 31, 2021, 2606 Mtbc isolates with an isoniazid or rifampicin resistance were identified in the NTRL biobank, of which, 1483 (56·9%) were from men, 1114 (42·7%) from women, and nine (0·4%) were unknown. Genome-based drug-resistant prediction classified 704 Mtbc strains as rifampicin resistant. 628 (89%) of the 704 Mtbc strains were classified MDR; of those, 146 (23%) were pre-extensively drug resistant (pre-XDR; additional fluoroquinolone resistance), and 24 (4%) extensively drug resistant (XDR; combined fluoroquinolone and bedaquiline resistance). Overall, 61 (9%) of 704 strains revealed resistance to bedaquiline: five (7%) of 76 rifampicin resistant plus bedaquiline resistant, 32 (7%) of 458 MDR plus bedaquiline resistant, and 24 (100%) of 24 XDR. Prevalence of bedaquiline resistance increased from 3% in 2016 to 14% in 2021. The cluster rate (12 single-nucleotide polymorphism threshold) was 42% for rifampicin-resistant strains, 78% for MDR strains, 94% for pre-XDR strains, and 96% for XDR Mtbc strains. 31 (4%) of 704 Mtbc strains, belonging to a diagnostic escape outbreak strain previously described in Eswatini (group_56), had an rpoB Ile491Phe mutation which is not detected by Xpert MTB/RIF (no other rpoB mutation). Of these, 23 (74%) showed additional resistance to bedaquiline, 13 (42%) had bedaquiline and fluoroquinolone resistance, and two (6%) were bedaquiline, fluoroquinolone, and delamanid resistant. INTERPRETATION: Pre-XDR resistance is highly prevalent among MDR Mtbc strains in Mozambique and so is bedaquiline resistance; and the frequency of bedaquiline resistance quadrupled over time and was found even in Mtbc strains without fluoroquinolone resistance. Importantly, strains with Ile491Phe mutation were frequent, accounting for 31% (n=10) of MDR plus bedaquiline-resistant strains and 54% (n=13) of XDR Mtbc strains. Given the current diagnostic algorithms and treatment regimens, both the emergence of rifampicin resistance due to Ile491Phe and bedaquiline resistance might jeopardise MDR tuberculosis prevention and care unless sequencing-based technology is rolled out. The potential cross border spread of diagnostic escape strains needs further investigation. FUNDING: The German Ministry of Health through the Seq_MDRTB-Net project, the Deutsche Forschungsgemeinschaft under Germany's Excellence Strategy Precision Medicine in Inflammation and the Research Training Group 2501 TransEvo, the Leibniz Science Campus Evolutionary Medicine of the Lung, and the German Ministry of Education and Research via the German Center for Infection Research.

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.

Implementation of targeted next-generation sequencing for the diagnosis of drug-resistant tuberculosis in low-resource settings: a programmatic model, challenges, and initial outcomes.

Targeted next-generation sequencing (tNGS) from clinical specimens has the potential to become a comprehensive tool for routine drug-resistance (DR) prediction of Mycobacterium tuberculosis complex strains (MTBC), the causative agent of tuberculosis (TB). However, TB mainly affects low- and middle-income countries, in which the implementation of new technologies have specific needs and challenges. We propose a model for programmatic implementation of tNGS in settings with no or low previous sequencing capacity/experience. We highlight the major challenges and considerations for a successful implementation. This model has been applied to build NGS capacity in Namibia, an upper middle-income country located in Southern Africa and suffering from a high-burden of TB and TB-HIV, and we describe herein the outcomes of this process.

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.

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 Sklodowska-Curie Actions.

Genomic epidemiology of SARS-CoV-2 during the first four waves in Mozambique.

Mozambique reported the first case of coronavirus disease 2019 (COVID-19) in March 2020 and it has since spread to all provinces in the country. To investigate the introductions and spread of SARS-CoV-2 in Mozambique, 1 142 whole genome sequences sampled within Mozambique were phylogenetically analyzed against a globally representative set, reflecting the first 25 months of the epidemic. The epidemic in the country was marked by four waves of infection, the first associated with B.1 ancestral lineages, while the Beta, Delta, and Omicron Variants of Concern (VOCs) were responsible for most infections and deaths during the second, third, and fourth waves. Large-scale viral exchanges occurred during the latter three waves and were largely attributed to southern African origins. Not only did the country remain vulnerable to the introductions of new variants but these variants continued to evolve within the borders of the country. Due to the Mozambican health system already under constraint, and paucity of data in Mozambique, there is a need to continue to strengthen and support genomic surveillance in the country as VOCs and Variants of interests (VOIs) are often reported from the southern African region.

Phylogroup-specific variation shapes the clustering of antimicrobial resistance genes and defence systems across regions of genome plasticity in Pseudomonas aeruginosa.

BACKGROUND: Pseudomonas aeruginosa is an opportunistic pathogen consisting of three phylogroups (hereafter named A, B, and C). Here, we assessed phylogroup-specific evolutionary dynamics across available and also new P. aeruginosa genomes. METHODS: In this genomic analysis, we first generated new genome assemblies for 18 strains of the major P. aeruginosa clone type (mPact) panel, comprising a phylogenetically diverse collection of clinical and environmental isolates for this species. Thereafter, we combined these new genomes with 1991 publicly available P. aeruginosa genomes for a phylogenomic and comparative analysis. We specifically explored to what extent antimicrobial resistance (AMR) genes, defence systems, and virulence genes vary in their distribution across regions of genome plasticity (RGPs) and "masked" (RGP-free) genomes, and to what extent this variation differs among the phylogroups. FINDINGS: We found that members of phylogroup B possess larger genomes, contribute a comparatively larger number of pangenome families, and show lower abundance of CRISPR-Cas systems. Furthermore, AMR and defence systems are pervasive in RGPs and integrative and conjugative/mobilizable elements (ICEs/IMEs) from phylogroups A and B, and the abundance of these cargo genes is often significantly correlated. Moreover, inter- and intra-phylogroup interactions occur at the accessory genome level, suggesting frequent recombination events. Finally, we provide here the mPact panel of diverse P. aeruginosa strains that may serve as a valuable reference for functional analyses. INTERPRETATION: Altogether, our results highlight distinct pangenome characteristics of the P. aeruginosa phylogroups, which are possibly influenced by variation in the abundance of CRISPR-Cas systems and are shaped by the differential distribution of other defence systems and AMR genes. FUNDING: German Science Foundation, Max-Planck Society, Leibniz ScienceCampus Evolutionary Medicine of the Lung, BMBF program Medical Infection Genomics, Kiel Life Science Postdoc Award.

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.

Origin and Global Expansion of Mycobacterium tuberculosis Complex Lineage 3.

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

Whole genome sequencing-based drug resistance predictions of multidrug-resistant Mycobacterium tuberculosis isolates from Tanzania.

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.

Sub-Lineage Specific Phenolic Glycolipid Patterns in the Mycobacterium tuberculosis Complex Lineage 1.

"Ancestral" Mycobacterium tuberculosis complex (MTBC) strains of Lineage 1 (L1, East African Indian) are a prominent tuberculosis (TB) cause in countries around the Indian Ocean. However, the pathobiology of L1 strains is insufficiently characterized. Here, we used whole genome sequencing (WGS) of 312 L1 strains from 43 countries to perform a characterization of the global L1 population structure and correlate this to the analysis of the synthesis of phenolic glycolipids (PGL) - known MTBC polyketide-derived virulence factors. Our results reveal the presence of eight major L1 sub-lineages, whose members have specific mutation signatures in PGL biosynthesis genes, e.g., pks15/1 or glycosyltransferases Rv2962c and/or Rv2958c. Sub-lineage specific PGL production was studied by NMR-based lipid profiling and strains with a completely abolished phenolphthiocerol dimycoserosate biosynthesis showed in average a more prominent growth in human macrophages. In conclusion, our results show a diverse population structure of L1 strains that is associated with the presence of specific PGL types. This includes the occurrence of mycoside B in one sub-lineage, representing the first description of a PGL in an M. tuberculosis lineage other than L2. Such differences may be important for the evolution of L1 strains, e.g., allowing adaption to different human populations.

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.

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 blaTEM-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.

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.