Use of Whole Genome Sequencing for Mycobacterium tuberculosis Complex Antimicrobial Susceptibility Testing: From Sequence Data to Resistance Profiles.

Whole genome sequencing of Mycobacterium tuberculosis complex (MTBC) isolates has been shown to provide accurate predictions for resistance and susceptibility for many first- and second-line anti-tuberculosis drugs. However, bioinformatic pipelines and mutation catalogs to predict antimicrobial resistances in MTBC isolates are often customized and detailed protocols are difficult to access. Here, we provide a step-by-step workflow for the processing and interpretation of short-read sequencing data and give an overview of available analysis pipelines.

Use of Whole Genome Sequencing for Mycobacterium tuberculosis Complex Antimicrobial Susceptibility Testing.

Whole genome sequencing (WGS) is becoming an important diagnostic tool for antimicrobial susceptibility testing of Mycobacterium tuberculosis complex (MTBC) isolates in many countries. WGS protocols usually start with the preparation of a DNA-library: the critical first step in the process. A DNA-library represents the genomic content of a DNA sample and consists of unique short DNA fragments. Although available DNA-library protocols come with manufacturer instructions, details of the entire process, including quality controls, instrument parameters, and run evaluations, often need to be developed and customized by each laboratory to implement WGS technology effectively. Here, we provide a detailed workflow for a DNA-library preparation based on an adapted Illumina protocol optimized for the reduction of reagent costs.

Phylogenomic and genomic analysis reveals unique and shared genetic signatures of Mycobacterium kansasii complex species.

Species belonging to the Mycobacterium kansasii complex (MKC) are frequently isolated from humans and the environment and can cause serious diseases. The most common MKC infections are caused by the species M. kansasii (sensu stricto), leading to tuberculosis-like disease. However, a broad spectrum of virulence, antimicrobial resistance and pathogenicity of these non-tuberculous mycobacteria (NTM) are observed across the MKC. Many genomic aspects of the MKC that relate to these broad phenotypes are not well elucidated. Here, we performed genomic analyses from a collection of 665 MKC strains, isolated from environmental, animal and human sources. We inferred the MKC pangenome, mobilome, resistome, virulome and defence systems and show that the MKC species harbours unique and shared genomic signatures. High frequency of presence of prophages and different types of defence systems were observed. We found that the M. kansasii species splits into four lineages, of which three are lowly represented and mainly in Brazil, while one lineage is dominant and globally spread. Moreover, we show that four sub-lineages of this most distributed M. kansasii lineage emerged during the twentieth century. Further analysis of the M. kansasii genomes revealed almost 300 regions of difference contributing to genomic diversity, as well as fixed mutations that may explain the M. kansasii's increased virulence and drug resistance.

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.

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.

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.

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.

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.