Compensatory evolution drives multidrug-resistant tuberculosis in Central Asia

Delineating the evolutionary pathway to multidrug-resistant outbreaks of a Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage

OBJECTIVES

We sought to capture the evolutionary itinerary of the Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, where it caused a major multidrug-resistant (MDR) tuberculosis outbreak in a context strictly negative for HIV infection.

METHODS

We combined whole genome sequencing and Bayesian approaches using a representative collection of drug-susceptible and drug-resistant L4.1.2.1/Haarlem clinical strains (n = 121) recovered from the outbreak region over 16 years.

RESULTS

In the absence of drug resistance, the L4.1.2.1/Haarlem sublineage showed a propensity for rapid transmission as witnessed by the high clustering (44.6%) and recent transmission rates (25%), as well as the reduced mean distance between genome pairs. The entire pool of L4.1.2.1/Haarlem MDR strains was found to be linked to either the aforementioned major outbreak (68 individuals, 2001-2016) or to a minor, newly uncovered outbreak (six cases, 2001-2011). Strikingly, the two outbreaks descended independently from a common ancestor that can be dated back to 1886.

CONCLUSIONS

Our data point to the intrinsic propensity for rapid transmission of the M. tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, linking the overall MDR tuberculosis epidemic to a single ancestor. These findings bring out the important role of the bacillus' genetic background in the emergence of successful MDR M. tuberculosis clones.

The epidemiology, transmission, diagnosis, and management of drug-resistant tuberculosis-lessons from the South African experience

Drug-resistant tuberculosis (DR-TB) threatens to derail tuberculosis control efforts, particularly in Africa where the disease remains out of control. The dogma that DR-TB epidemics are fueled by unchecked rates of acquired resistance in inadequately treated or non-adherent individuals is no longer valid in most high DR-TB burden settings, where community transmission is now widespread. A large burden of DR-TB in Africa remains undiagnosed due to inadequate access to diagnostic tools that simultaneously detect tuberculosis and screen for resistance. Furthermore, acquisition of drug resistance to new and repurposed drugs, for which diagnostic solutions are not yet available, presents a major challenge for the implementation of novel, all-oral, shortened (6-9 months) treatment. Structural challenges including poverty, stigma, and social distress disrupt engagement in care, promote poor treatment outcomes, and reduce the quality of life for people with DR-TB. We reflect on the lessons learnt from the South African experience in implementing state-of-the-art advances in diagnostic solutions, deploying recent innovations in pharmacotherapeutic approaches for rapid cure, understanding local transmission dynamics and implementing interventions to curtail DR-TB transmission, and in mitigating the catastrophic socioeconomic costs of DR-TB. We also highlight globally relevant and locally responsive research priorities for achieving DR-TB control in South Africa.

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.

Transmission dynamics and phylogeography of Mycobacterium tuberculosis in China based on whole-genome phylogenetic analysis

OBJECTIVES

This study aimed to describe the lineage-specific transmissibility and epidemiological migration of Mycobacterium tuberculosis in China.

METHODS

We curated a large set of whole-genome sequences from 3204 M. tuberculosis isolates, including thousands of newly sequenced genomes, and applied a series of metrics to compare the transmissibility of M. tuberculosis strains between lineages and sublineages. The countrywide transmission patterns of major lineages were explored.

RESULTS

We found that lineage 2 (L2) was the most prevalent lineage in China (85.7%), with the major sublineage 2.2.1 (80.9%), followed by lineage 4 (L4) (13.8%), which comprises major sublineages 4.2 (1.5%), 4.4 (6.2%) and 4.5 (5.8%). We showed evidence for frequent cross-regional spread and large cluster formation of L2.2.1 strains, whereas L4 strains were relatively geographically restricted in China. Next, we applied a series of genomic indices to evaluate M. tuberculosis strain transmissibility and uncovered higher transmissibility of L2.2.1 compared with the L2.2.2 and L4 sublineages. Phylogeographic analysis showed that southern, eastern, and northern China were highly connected regions for countrywide L2.2.1 strain spread.

CONCLUSIONS

The present study provides insights into the different transmission and migration patterns of the major M. tuberculosis lineages in China and highlights that transmissible L2.2.1 is a threat to tuberculosis control.

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.

Compensatory mutations are associated with increased in vitro growth in resistant clinical samples of Mycobacterium tuberculosis

Mutations in Mycobacterium tuberculosis associated with resistance to antibiotics often come with a fitness cost for the bacteria. Resistance to the first-line drug rifampicin leads to lower competitive fitness of M. tuberculosis populations when compared to susceptible populations. This fitness cost, introduced by resistance mutations in the RNA polymerase, can be alleviated by compensatory mutations (CMs) in other regions of the affected protein. CMs are of particular interest clinically since they could lock in resistance mutations, encouraging the spread of resistant strains worldwide. Here, we report the statistical inference of a comprehensive set of CMs in the RNA polymerase of M. tuberculosis, using over 70 000 M. tuberculosis genomes that were collated as part of the CRyPTIC project. The unprecedented size of this data set gave the statistical tests more power to investigate the association of putative CMs with resistance-conferring mutations. Overall, we propose 51 high-confidence CMs by means of statistical association testing and suggest hypotheses for how they exert their compensatory mechanism by mapping them onto the protein structure. In addition, we were able to show an association of CMs with higher in vitro growth densities, and hence presumably with higher fitness, in resistant samples in the more virulent M. tuberculosis lineage 2. Our results suggest the association of CM presence with significantly higher in vitro growth than for wild-type samples, although this association is confounded with lineage and sub-lineage affiliation. Our findings emphasize the integral role of CMs and lineage affiliation in resistance spread and increases the urgency of antibiotic stewardship, which implies accurate, cheap and widely accessible diagnostics for M. tuberculosis infections to not only improve patient outcomes but also prevent the spread of resistant strains.

Novel and reported compensatory mutations in rpoABC genes found in drug resistant tuberculosis outbreaks

Background

Rifampicin (RIF) is a key first-line drug used to treat tuberculosis, a primarily pulmonary disease caused by Mycobacterium tuberculosis. RIF resistance is caused by mutations in rpoB, at the cost of slower growth and reduced transcription efficiency. Antibiotic resistance to RIF is prevalent despite this fitness cost. Compensatory mutations in rpoABC genes have been shown to alleviate the fitness cost of rpoB:S450L, explaining how RIF resistant strains harbor this mutation can spread so rapidly. Unfortunately, the full set of RIF compensatory mutations is still unknown, particularly those compensating for rarer RIF resistance mutations.

Objectives

We performed an association study on a globally representative set of 4,309 whole genome sequenced clinical M. tuberculosis isolates to identify novel putative compensatory mutations, determine the prevalence of known and previously reported putative compensatory mutations, and determine which RIF resistance markers associate with these compensatory mutations.

Results and conclusions

Of the 1,079 RIF resistant isolates, 638 carried previously reported putative and high-probability compensatory mutations. Our strict criteria identified 46 additional mutations in rpoABC for which no strong prior evidence of their compensatory role exists. Of these, 35 have previously been reported. As such, our independent corroboration adds to the mounting evidence that these 35 also carry a compensatory role. The remaining 11 are novel putative compensatory markers, reported here for the first time. Six of these 11 novel putative compensatory mutations had two or more mutation events. Most compensatory mutations appear to be specifically compensating for the fitness loss due to rpoB:S450L. However, an outbreak of 22 closely related isolates each carried three rpoB mutations, the rare RIFR markers D435G and L452P and the putative compensatory mutation I1106T. This suggests compensation may require specific combinations of rpoABC mutations. Here, we report only mutations that met our very strict criteria. It is highly likely that many additional rpoABC mutations compensate for rare resistance-causing mutations and therefore did not carry the statistical power to be reported here. These findings aid in the identification of RIF resistant M. tuberculosis strains with restored fitness, which pose a greater risk of causing resistant outbreaks.

Genomic Characterization of Drug-Resistant Mycobacterium tuberculosis L2/Beijing Isolates from Astana, Kazakhstan

Kazakhstan ranks among the countries with the highest number of MDR-TB patients per 100,000 population worldwide. The successful transmission of local MDR strains of Mycobacterium tuberculosis (Mtb) poses a significant threat to disease control. In this study, we employed whole-genome sequencing to examine drug resistance, compensatory mutations, population structure, and transmission patterns in a sample of 24 clinical isolates of L2/Beijing Mtb collected in Astana, Kazakhstan between 2021 and 2022. The genotypic prediction of Mtb susceptibility to anti-TB agents was consistent with the phenotypic susceptibility, except for bedaquiline. An analysis of resistance-associated genes characterized most of the isolates as pre-extensively drug-resistant tuberculosis (pre-XDR-TB) (n = 15; 62.5%). The phylogenetic analysis grouped the isolates into four transmission clusters; the dominant cluster was assigned to the "aggressive" Central Asia outbreak (CAO) clade of L2/Beijing (n = 15; 62.5%). Thirteen mutations with putative compensatory effects were observed exclusively in Mtb isolates containing the rpoB S450L mutation. The putative compensatory mutations had a stabilizing effect on RpoABC protein stability and dynamics. The high prevalence of the CAO clade in the population structure of Mtb may explain the rapid spread of MDR-TB in Kazakhstan.

A multifaceted interplay between virulence, drug resistance, and the phylogeographic landscape of Mycobacterium tuberculosis

Latin-American Mediterranean (LAM) family is one of the most significant and global genotypes of Mycobacterium tuberculosis. Here, we used the murine model to study the virulence and lethality of the genetically and epidemiologically distinct LAM strains. The pathobiological characteristics of the four LAM strains (three drug resistant and one drug susceptible) and the susceptible reference strain H37Rv were studied in the C57BL/6 mouse model. The whole-genome sequencing was performed using the HiSeq Illumina platform, followed by bioinformatics and phylogenetic analysis. The susceptible strain H37Rv showed the highest virulence. Drug-susceptible LAM strain (spoligotype SIT264) was more virulent than three multidrug-resistant (MDR) strains (SIT252, SIT254, and SIT266). All three MDR isolates were low lethal, while the susceptible isolate and H37Rv were moderately/highly lethal. Putting the genomic, phenotypic, and virulence features of the LAM strains/spoligotypes in the context of their dynamic phylogeography over 20 years reveals three types of relationships between virulence, resistance, and transmission. First, the most virulent and more lethal drug-susceptible SIT264 increased its circulation in parts of Russia. Second, moderately virulent and pre-XDR SIT266 was prevalent in Belarus and continues to be visible in North-West Russia. Third, the low virulent and MDR strain SIT252 previously considered as emerging has disappeared from the population. These findings suggest that strain virulence impacts the transmission, irrespective of drug resistance properties. The increasing circulation of susceptible but more virulent and lethal strains implies that personalized TB treatment should consider not only resistance but also the virulence of the infecting M. tuberculosis strains. IMPORTANCE The study is multidisciplinary and investigates the epidemically/clinically important and global lineage of Mycobacterium tuberculosis, named Latin-American-Mediterranean (LAM), yet insufficiently studied with regard to its pathobiology. We studied different LAM strains (epidemic vs endemic and resistant vs susceptible) in the murine model and using whole-genome analysis. We also collected long-term, 20-year data on their prevalence in Eurasia. The findings are both expected and unexpected. (i) We observe that a drug-susceptible but highly virulent strain increased its prevalence. (ii) By contrast, the multidrug-resistant (MDR) but low-virulent, low-lethal strain (that we considered as emerging 15 years ago) has almost disappeared. (iii) Finally, an intermediate case is the MDR strain with moderate virulence that continues to circulate. We conclude that (i) the former and latter strains are the most hazardous and require close epidemiological monitoring, and (ii) personalized TB treatment should consider not only drug resistance but also the virulence of the infecting strains and development of anti-virulence drugs is warranted.

Host-pathogen relationship in retreated tuberculosis with major rifampicin resistance-conferring mutations

Introduction

It is assumed that host defense systems eliminating the pathogen and regulating tissue damage make a strong impact on the outcome of tuberculosis (TB) disease and that these processes are affected by rifampicin (RIF) resistance-conferring mutations of Mycobacterium tuberculosis (Mtb). However, the host responses to the pathogen harboring different mutations have not been studied comprehensively in clinical settings. We analyzed clinico-epidemiological factors and blood transcriptomic signatures associated with major rpoB mutations conferring RIF resistance in a cohort study.

Methods

Demographic data were collected from 295 active pulmonary TB patients with treatment history in Hanoi, Vietnam. When recruited, drug resistance-conferring mutations and lineage-specific variations were identified using whole-genome sequencing of clinical Mtb isolates. Before starting retreatment, total RNA was extracted from the whole blood of HIV-negative patients infected with Mtb that carried either the rpoB H445Y or rpoB S450L mutation, and the total RNA was subjected to RNA sequencing after age-gender matching. The individual RNA expression levels in the blood sample set were also measured using real-time RT-PCR. Logistic and linear regression models were used to assess possible associations.

Results

In our cohort, rpoB S450L and rpoB H445Y were major RIF resistance-conferring mutations [32/87 (36.8%) and 15/87 (17.2%), respectively]. H445Y was enriched in the ancient Beijing genotype and was associated with nonsynonymous mutations of Rv1830 that has been reported to regulate antibiotic resilience. H445Y was also more frequently observed in genetically clustered strains and in samples from patients who had received more than one TB treatment episode. According to the RNA sequencing, gene sets involved in the interferon-γ and-α pathways were downregulated in H445Y compared with S450L. The qRT-PCR analysis also confirmed the low expression levels of interferon-inducible genes, including BATF2 and SERPING1, in the H445Y group, particularly in patients with extensive lesions on chest X-ray.

Discussion

Our study results showed that rpoB mutations as well as Mtb sublineage with additional genetic variants may have significant effects on host response. These findings strengthen the rationale for investigation of host-pathogen interactions to develop countermeasures against epidemics of drug-resistant TB.

Effect of compensatory evolution in the emergence and transmission of rifampicin-resistant Mycobacterium tuberculosis in Cape Town, South Africa: a genomic epidemiology study

BACKGROUND

Experimental data show that drug-resistance-conferring mutations are often associated with a decrease in the replicative fitness of bacteria in vitro, and that this fitness cost can be mitigated by compensatory mutations; however, the role of compensatory evolution in clinical settings is less clear. We assessed whether compensatory evolution was associated with increased transmission of rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa.

METHODS

We did a genomic epidemiological study by analysing available M tuberculosis isolates and their associated clinical data from individuals routinely diagnosed with rifampicin-resistant tuberculosis in primary care and hospitals in Khayelitsha, Cape Town, South Africa. Isolates were collected as part of a previous study. All individuals diagnosed with rifampicin-resistant tuberculosis and with linked biobanked specimens were included in this study. We applied whole-genome sequencing, Bayesian reconstruction of transmission trees, and phylogenetic multivariable regression analysis to identify individual and bacterial factors associated with the transmission of rifampicin-resistant M tuberculosis strains.

FINDINGS

Between Jan 1, 2008, and Dec 31, 2017, 2161 individuals were diagnosed with multidrug-resistant or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa. Whole-genome sequences were available for 1168 (54%) unique individual M tuberculosis isolates. Compensatory evolution was associated with smear-positive pulmonary disease (adjusted odds ratio 1·49, 95% CI 1·08-2·06) and a higher number of drug-resistance-conferring mutations (incidence rate ratio 1·38, 95% CI 1·28-1·48). Compensatory evolution was also associated with increased transmission of rifampicin-resistant disease between individuals (adjusted odds ratio 1·55; 95% CI 1·13-2·12), independent of other patient and bacterial factors.

INTERPRETATION

Our findings suggest that compensatory evolution enhances the in vivo fitness of drug-resistant M tuberculosis genotypes, both within and between patients, and that the in vitro replicative fitness of rifampicin-resistant M tuberculosis measured in the laboratory correlates with the bacterial fitness measured in clinical settings. These results emphasise the importance of enhancing surveillance and monitoring efforts to prevent the emergence of highly transmissible clones capable of rapidly accumulating new drug resistance mutations. This concern becomes especially crucial at present, because treatment regimens incorporating novel drugs are being implemented.

FUNDING

Funding for this study was provided by a Swiss and South Africa joint research award (grant numbers 310030_188888, CRSII5_177163, and IZLSZ3_170834), the European Research Council (grant number 883582), and a Wellcome Trust fellowship (to HC; reference number 099818/Z/12/Z). ZS-D was funded through a PhD scholarship from the South African National Research Foundation and RMW was funded through the South African Medical Research Council.

The fitness cost of horizontally transferred and mutational antimicrobial resistance in Escherichia coli

Antimicrobial resistance (AMR) in bacteria implies a tradeoff between the benefit of resistance under antimicrobial selection pressure and the incurred fitness cost in the absence of antimicrobials. The fitness cost of a resistance determinant is expected to depend on its genetic support, such as a chromosomal mutation or a plasmid acquisition, and on its impact on cell metabolism, such as an alteration in an essential metabolic pathway or the production of a new enzyme. To provide a global picture of the factors that influence AMR fitness cost, we conducted a systematic review and meta-analysis focused on a single species, Escherichia coli. By combining results from 46 high-quality studies in a multilevel meta-analysis framework, we find that the fitness cost of AMR is smaller when provided by horizontally transferable genes such as those encoding beta-lactamases, compared to mutations in core genes such as those involved in fluoroquinolone and rifampicin resistance. We observe that the accumulation of acquired AMR genes imposes a much smaller burden on the host cell than the accumulation of AMR mutations, and we provide quantitative estimates of the additional cost of a new gene or mutation. These findings highlight that gene acquisition is more efficient than the accumulation of mutations to evolve multidrug resistance, which can contribute to the observed dominance of horizontally transferred genes in the current AMR epidemic.

Insight into Population Structure and Drug Resistance of Pediatric Tuberculosis Strains from China and Russia Gained through Whole-Genome Sequencing

This study aimed to determine phenotypic and genotypic drug resistance patterns of Mycobacterium tuberculosis strains from children with tuberculosis (TB) in China and Russia, two high-burden countries for multi/extensively-drug resistant (MDR/XDR) TB. Whole-genome sequencing data of M. tuberculosis isolates from China (n = 137) and Russia (n = 60) were analyzed for phylogenetic markers and drug-resistance mutations, followed by comparison with phenotypic susceptibility data. The Beijing genotype was detected in 126 Chinese and 50 Russian isolates. The Euro-American lineage was detected in 10 Russian and 11 Chinese isolates. In the Russian collection, the Beijing genotype and Beijing B0/W148-cluster were dominated by MDR strains (68% and 94%, respectively). Ninety percent of B0/W148 strains were phenotypically pre-XDR. In the Chinese collection, neither of the Beijing sublineages was associated with MDR/pre-XDR status. MDR was mostly caused by low fitness cost mutations (rpoB S450L, katG S315T, rpsL K43R). Chinese rifampicin-resistant strains demonstrated a higher diversity of resistance mutations than Russian isolates (p = 0.003). The rifampicin and isoniazid resistance compensatory mutations were detected in some MDR strains, but they were not widespread. The molecular mechanisms of M. tuberculosis adaptation to anti-TB treatment are not unique to the pediatric strains, but they reflect the general situation with TB in Russia and China.

Vaccines against Tuberculosis: Where Are We Now?

Tuberculosis (TB) is among the top 10 leading causes of death in low-income countries. Statistically, TB kills more than 30,000 people each week and leads to more deaths than any other infectious disease, such as acquired immunodeficiency syndrome (AIDS) and malaria. TB treatment is largely dependent on BCG vaccination and impacted by the inefficacy of drugs, absence of advanced vaccines, misdiagnosis improper treatment, and social stigma. The BCG vaccine provides partial effectiveness in demographically distinct populations and the prevalence of multidrug-resistant (MDR) and extensively drug-resistant (XDR) TB incidences demands the design of novel TB vaccines. Various strategies have been employed to design vaccines against TB, such as: (a) The protein subunit vaccine; (b) The viral vector vaccine; (c) The inactivation of whole-cell vaccine, using related mycobacteria, (d) Recombinant BCG (rBCG) expressing Mycobacterium tuberculosis (M.tb) protein or some non-essential gene deleted BCG. There are, approximately, 19 vaccine candidates in different phases of clinical trials. In this article, we review the development of TB vaccines, their status and potential in the treatment of TB. Heterologous immune responses generated by advanced vaccines will contribute to long-lasting immunity and might protect us from both drug-sensitive and drug-resistant TB. Therefore, advanced vaccine candidates need to be identified and developed to boost the human immune system against TB.

The relative transmission fitness of multidrug-resistant Mycobacterium tuberculosis in a drug resistance hotspot

Multidrug-resistant tuberculosis (MDR-TB) is among the most frequent causes of death due to antimicrobial resistance. Although only 3% of global TB cases are MDR, geographical hotspots with up to 40% of MDR-TB have been observed in countries of the former Soviet Union. While the quality of TB control and patient-related factors are known contributors to such hotspots, the role of the pathogen remains unclear. Here we show that in the country of Georgia, a known hotspot of MDR-TB, MDR Mycobacterium tuberculosis strains of lineage 4 (L4) transmit less than their drug-susceptible counterparts, whereas most MDR strains of L2 suffer no such defect. Our findings further indicate that the high transmission fitness of these L2 strains results from epistatic interactions between the rifampicin resistance-conferring mutation RpoB S450L, compensatory mutations in the RNA polymerase, and other pre-existing genetic features of L2/Beijing clones that circulate in Georgia. We conclude that the transmission fitness of MDR M. tuberculosis strains is heterogeneous, but can be as high as drug-susceptible forms, and that such highly drug-resistant and transmissible strains contribute to the emergence and maintenance of hotspots of MDR-TB. As these strains successfully overcome the metabolic burden of drug resistance, and given the ongoing rollout of new treatment regimens against MDR-TB, proper surveillance should be implemented to prevent these strains from acquiring resistance to the additional drugs.

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.

Genetic Diversity and Primary Drug Resistance of Mycobacterium tuberculosis Beijing Genotype Strains in Northwestern Russia

The Beijing genotype is the main family of Mycobacterium tuberculosis in Russia. We analyzed its diversity and drug resistance in provinces across Northwestern Russia to identify the epidemiologically relevant Beijing strains. The study collection included 497 isolates from newly-diagnosed tuberculosis (TB) patients. Bacterial isolates were subjected to drug-susceptibility testing and genotyping. The Beijing genotype was detected in 57.5% (286/497); 50% of the Beijing strains were multidrug-resistant (MDR). Central Asian/Russian and B0/W148 groups included 176 and 77 isolates, respectively. MDR was more frequent among B0/W148 strains compared to Central Asian/Russian strains (85.7% vs. 40.3%, p < 0.0001). Typing of 24 minisatellite loci of Beijing strains revealed 82 profiles; 230 isolates were in 23 clusters. The largest Central Asian/Russian types were 94-32 (n = 75), 1065-32 (n = 17), and 95-32 (n = 12). B0/W148 types were 100-32 (n = 59) and 4737-32 (n = 5). MDR was more frequent in types 1065-32 (88.2%), 100-32 (83.1%), and 4737-32 (100%). In contrast, type 9391-32 (n = 9) included only drug-susceptible strains. To conclude, M. tuberculosis Beijing genotype is dominant in Northwestern Russia, and an active transmission of overwhelmingly MDR B0/W148 types explains the reported increase of MDR-TB. The presence of MDR-associated minor variants (type 1071-32/ancient Beijing and Central Asia Outbreak strain) in some of the studied provinces also requires attention.

Evolutionary history and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage, Tunisia

BACKGROUND

To infer the origin and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage in a Mediterranean country, Tunisia, where it predominates.

METHODS

We combined Bayesian (STRUCTURE) and maximum likelihood (MIGRAINE) estimation approaches based on a global 24-loci mycobacterial interspersed repetitive units-variable numbers of tandem repeats (MIRU-VNTR24) genotyping dataset consisting of 1573 L4.3/LAM clinical strains from four continents, including 252 isolates originating from Tunisia.

RESULTS

Phylogenetic analyses coupled with Bayesian estimations suggested that the most predominant L4.3/LAM subpopulation in Tunisia (65.07%), which is dominated by a single clonal complex, TUN4.3_CC1 (94.51%), has evolved from an ancestral pool that is restricted to Europe and Africa, contrasting with the remaining L4.3/LAM subpopulations whose ancestry was traced all over the word. Maximum likelihood analysis revealed that TUN4.3_CC1 has been undergoing a demographic expansion since 131 years ago (CI95%: 90.7-205), thus explaining its preponderance relative to the second most predominant CC, TUN4.3_CC2, whose population was found under contraction.

CONCLUSIONS

The preponderance of L4.3/LAM in Tunisia stems from a 130-year expansion process of a locally evolved clone.

Whole-Genome Sequencing for Resistance Prediction and Transmission Analysis of Mycobacterium tuberculosis Complex Strains from Namibia

Namibia is among 30 countries with a high burden of tuberculosis (TB), with an estimated incidence of 460 per 100,000 population and around 800 new multidrug-resistant (MDR) TB cases per year. Still, data on the transmission and evolution of drug-resistant Mycobacterium tuberculosis complex (Mtbc) strains are not available. Whole-genome sequencing data of 136 rifampicin-resistant (RIFr) Mtbc strains obtained from 2016 to 2018 were used for phylogenetic classification, resistance prediction, and cluster analysis and linked with phenotypic drug susceptibility testing (pDST) data. Roughly 50% of the strains investigated were resistant to all first-line drugs. Furthermore, 13% of the MDR Mtbc strains were already pre-extensively drug resistant (pre-XDR). The cluster rates were high, at 74.6% among MDR and 85% among pre-XDR strains. A significant proportion of strains had borderline resistance-conferring mutations, e.g., inhA promoter mutations or rpoB L430P. Accordingly, 25% of the RIFr strains tested susceptible by pDST. Finally, we determined a potentially new bedaquiline resistance mutation (Rv0678 D88G) occurring in two independent clusters. High rates of resistance to first-line drugs in line with emerging pre-XDR and likely bedaquiline resistance linked with the ongoing recent transmission of MDR Mtbc clones underline the urgent need for the implementation of interventions that allow rapid diagnostics to break MDR TB transmission chains in the country. A borderline RIFr mutation in the dominant outbreak strain causing discrepancies between phenotypic and genotypic resistance testing results may require breakpoint adjustments but also may allow individualized regimens with high-dose treatment. IMPORTANCE The transmission of drug-resistant tuberculosis (TB) is a major problem for global TB control. Using genome sequencing, we showed that 13% of the multidrug-resistant (MDR) M. tuberculosis complex strains from Namibia are already pre-extensively drug resistant (pre-XDR), which is substantial in an African setting. Our data also indicate that the ongoing transmission of MDR and pre-XDR strains contributes significantly to the problem. In contrast to other settings with higher rates of drug resistance, we found a high proportion of strains having so-called borderline low-level resistance mutations, e.g., inhA promoter mutations or rpoB L430P. This led to the misclassification of 25% of the rifampicin-resistant strains as susceptible by phenotypic drug susceptibility testing. This observation potentially allows individualized regimens with high-dose treatment as a potential option for patients with few treatment options. We also found a potentially new bedaquiline resistance mutation in rv0678.

Evolution of Mycobacterium tuberculosis drug resistance in the genomic era

Mycobacterium tuberculosis has acquired drug resistance to all drugs that have been used against it, including those only recently introduced into clinical practice. Compared to other bacteria, it has a well conserved genome due to its role as an obligate human pathogen that has adapted to a niche over five to ten thousand years. These features facilitate reconstruction and dating of M. tuberculosis phylogenies, giving key insights into how resistance has been acquired and spread globally. Resistance to each new drug has occurred within five to ten years of clinical use and has occurred even more rapidly with recently introduced drugs. In most cases, resistance-conferring mutations come with a fitness cost, but this can be overcome by compensatory mutations which restore fitness to that of wild-type bacteria. It is likely that M. tuberculosis acquires drug resistance while maintaining limited genomic variability due the generation of low frequency within-host variation, combined with ongoing purifying selection causing loss of variants without a clear fitness advantage. However, variants that do confer an advantage, such as drug resistance, can increase in prevalence amongst all bacteria within a host and become the dominant clone. These resistant strains can then be transmitted leading to primary drug resistant infection in a new host. As many countries move towards genomic methods for diagnosis of M. tuberculosis infection and drug resistance, it is important to be aware of the implications for the evolution of resistance. Currently, understanding of resistance-conferring mutations is incomplete, and some targeted genetic diagnostics create their own selective pressures. We discuss an example where a rifampicin resistance-conferring mutation which was not routinely covered by standard testing became dominant. Finally, resistance to new drugs such as bedaquiline and delamanid is caused by individually rare mutations occurring across a large mutational genomic target that have been detected over a short time, and do not provide statistical power for genotype-phenotype correlation – in contrast to longer-established drugs that form the backbone of drug-sensitive antituberculosis therapy. Therefore, we need a different approach to identify resistance-conferring mutations of new drugs before their resistance becomes widespread, abrogating their usefulness.

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.

Mechanisms of antibiotic action shape the fitness landscapes of resistance mutations

Antibiotic-resistant pathogens are a major public health threat. A deeper understanding of how an antibiotic's mechanism of action influences the emergence of resistance would aid in the design of new drugs and help to preserve the effectiveness of existing ones. To this end, we developed a model that links bacterial population dynamics with antibiotic-target binding kinetics. Our approach allows us to derive mechanistic insights on drug activity from population-scale experimental data and to quantify the interplay between drug mechanism and resistance selection. We find that both bacteriostatic and bactericidal agents can be equally effective at suppressing the selection of resistant mutants, but that key determinants of resistance selection are the relationships between the number of drug-inactivated targets within a cell and the rates of cellular growth and death. We also show that heterogeneous drug-target binding within a population enables resistant bacteria to evolve fitness-improving secondary mutations even when drug doses remain above the resistant strain's minimum inhibitory concentration. Our work suggests that antibiotic doses beyond this "secondary mutation selection window" could safeguard against the emergence of high-fitness resistant strains during treatment.

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.

Genome-Wide Study of Drug Resistant Mycobacterium tuberculosis and Its Intra-Host Evolution during Treatment

The emergence of drug resistant Mycobacterium tuberculosis (MTB) strains has become a global public health problem, while, at the same time, there has been development of new antimicrobial agents. The main goals of this study were to determine new variants associated with drug resistance in MTB and to observe which polymorphisms emerge in MTB genomes after anti-tuberculosis treatment. We performed whole-genome sequencing of 152 MTB isolates including 70 isolates as 32 series of pre- and post-treatment MTB. Based on genotypes and phenotypic drug susceptibility, we conducted phylogenetic convergence-based genome-wide association study (GWAS) with streptomycin-, isoniazid-, rifampicin-, ethambutol-, fluoroquinolones-, and aminoglycosides-resistant MTB against susceptible ones. GWAS revealed statistically significant associations of SNPs within Rv2820c, cyp123 and indels in Rv1269c, Rv1907c, Rv1883c, Rv2407, Rv3785 genes with resistant MTB phenotypes. Comparisons of serial isolates showed that treatment induced different patterns of intra-host evolution. We found indels within Rv1435c and ppsA that were not lineage-specific. In addition, Beijing-specific polymorphisms within Rv0036c, Rv0678, Rv3433c, and dop genes were detected in post-treatment isolates. The appearance of Rv3785 frameshift insertion in 2 post-treatment strains compared to pre-treatment was also observed. We propose that the insertion within Rv3785, which was a GWAS hit, might affect cell wall biosynthesis and probably mediates a compensatory mechanism in response to treatment. These results may shed light on the mechanisms of MTB adaptation to chemotherapy and drug resistance formation.

Understanding drivers of phylogenetic clustering and terminal branch lengths distribution in epidemics of Mycobacterium tuberculosis

Detecting factors associated with transmission is important to understand disease epidemics, and to design effective public health measures. Clustering and terminal branch lengths (TBL) analyses are commonly applied to genomic data sets of Mycobacterium tuberculosis (MTB) to identify sub-populations with increased transmission. Here, I used a simulation-based approach to investigate what epidemiological processes influence the results of clustering and TBL analyses, and whether differences in transmission can be detected with these methods. I simulated MTB epidemics with different dynamics (latency, infectious period, transmission rate, basic reproductive number R0, sampling proportion, sampling period, and molecular clock), and found that all considered factors, except for the length of the infectious period, affect the results of clustering and TBL distributions. I show that standard interpretations of this type of analyses ignore two main caveats: (1) clustering results and TBL depend on many factors that have nothing to do with transmission, (2) clustering results and TBL do not tell anything about whether the epidemic is stable, growing, or shrinking, unless all the additional parameters that influence these metrics are known, or assumed identical between sub-populations. An important consequence is that the optimal SNP threshold for clustering depends on the epidemiological conditions, and that sub-populations with different epidemiological characteristics should not be analyzed with the same threshold. Finally, these results suggest that different clustering rates and TBL distributions, that are found consistently between different MTB lineages, are probably due to intrinsic bacterial factors, and do not indicate necessarily differences in transmission or evolutionary success.

Origin and Global Expansion of Mycobacterium tuberculosis Complex Lineage 3

Simple Summary

Tuberculosis still causes 1.5 million deaths annually and is mainly caused by Mycobacterium tuberculosis complex strains belonging to three evolutionary modern lineages (Lineages 2–4). While Lineage 2 and Lineage 4 virtually conquered the world, Lineage 3 is particularly successful in Northern and Eastern Africa, as well as in Southern Asia, the suspected evolutionary origin of these strains. Here, we sought to understand how Lineage 3 strains came to the African continent. To this end, we performed routine genotyping to characterize over 2500 clinical isolates from 38 countries. We then selected a representative collection of 373 isolates for a whole-genome analysis and a modeling approach to infer the geographic origin of different sublineages. In fact, the origin of Lineage 3 could be located in India, and we found evidence for independent introductions of four distinct sublineages into North/East Africa, in line with known ancient exchanges and migrations between both world regions. Our study illustrates that the evolutionary history of humans and their pathogens are closely connected and further provides a systematic understanding of the genomic diversity of Lineage 3, which could be important for the development of new tuberculosis vaccines or new therapeutics.

Abstract

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.

Comprehensive genomic analysis of Mycobacterium tuberculosis reveals limited impact of high-fitness genotypes on MDR-TB transmission

OBJECTIVES

Environmental and host-related factors that contribute to the transmission of multidrug-resistant tuberculosis (MDR-TB) have become an increasing concern, but the impact of bacterial genetic factors associated with bacterial fitness on MDR-TB transmission is poorly understood. Here, we present a global view of the correlation between common fitness-related genotypes and MDR-TB transmission by analyzing a representative number of MDR-TB isolates.

METHODS

We assembled a global whole genome sequencing (WGS) dataset of MDR-TB strains collected through retrospective cohorts or population-based approaches using public databases and literature curation. WGS-based clusters were defined as groups of strains with genomic difference of ≤ 5 SNPs.

RESULTS

We curated high-quality WGS data of 4696 MDR-TB isolates from 17 countries with a mean clustering rate of 48% (range 0-100%). Correlational analysis showed that increased risk of MDR-TB strain clustering was not associated with compensatory mutations (OR 1.07, 95% CI 0.72-1.59), low-fitness cost drug-resistant mutations (katG S315T: OR 1.42, 95% CI 0.82-2.47; rpoB S450L: OR 1.26, 95% CI 0.87-1.83) or Lineage 2 (OR 1.50, 95% CI 0.95-2.39).

CONCLUSIONS

The factors most commonly thought to increase bacterial fitness were not significantly associated with increased MDR-TB transmission, and thus do not appear to be major contributors to the current epidemic of MDR-TB.

Multidrug-resistant tuberculosis control in Rwanda overcomes a successful clone that causes most disease over a quarter century

SUMMARY BACKGROUND

Multidrug-resistant (MDR) tuberculosis (TB) poses an important challenge in TB management and control. Rifampicin resistance (RR) is a solid surrogate marker of MDR-TB. We investigated the RR-TB clustering rates, bacterial population dynamics to infer transmission dynamics, and the impact of changes to patient management on these dynamics over 27 years in Rwanda.

METHODS

We analysed whole genome sequences of a longitudinal collection of nationwide RR-TB isolates. The collection covered three important periods: before programmatic management of MDR-TB (PMDT; 1991-2005), the early PMDT phase (2006-2013), in which rifampicin drug-susceptibility testing (DST) was offered to retreatment patients only, and the consolidated phase (2014-2018), in which all bacteriologically confirmed TB patients had rifampicin DST done mostly via Xpert MTB/RIF assay. We constructed clusters based on a 5 SNP cut-off and resistance conferring SNPs. We used Bayesian modelling for dating and population size estimations, TransPhylo to estimate the number of secondary cases infected by each patient, and multivariable logistic regression to assess predictors of being infected by the dominant clone.

RESULTS

Of 308 baseline RR-TB isolates considered for transmission analysis, the clustering analysis grouped 259 (84.1%) isolates into 13 clusters. Within these clusters, a single dominant clone was discovered containing 213 isolates (82.2% of clustered and 69.1% of all RR-TB), which we named the "Rwanda Rifampicin-Resistant clone" (R3clone). R3clone isolates belonged to Ugandan sub-lineage 4.6.1.2 and its rifampicin and isoniazid resistance were conferred by the Ser450Leu mutation in rpoB and Ser315Thr in katG genes, respectively. All R3clone isolates had Pro481Thr, a putative compensatory mutation in the rpoC gene that likely restored its fitness. The R3clone was estimated to first arise in 1987 and its population size increased exponentially through the 1990s', reaching maximum size (∼84%) in early 2000 s', with a declining trend since 2014. Indeed, the highest proportion of R3clone (129/157; 82·2%, 95%CI: 75·3-87·8%) occurred between 2000 and 13, declining to 64·4% (95%CI: 55·1-73·0%) from 2014 onward. We showed that patients with R3clone detected after an unsuccessful category 2 treatment were more likely to generate secondary cases than patients with R3clone detected after an unsuccessful category 1 treatment regimen.

CONCLUSIONS

RR-TB in Rwanda is largely transmitted. Xpert MTB/RIF assay as first diagnostic test avoids unnecessary rounds of rifampicin-based TB treatment, thus preventing ongoing transmission of the dominant R3clone. As PMDT was intensified and all TB patients accessed rifampicin-resistance testing, the nationwide R3clone burden declined. To our knowledge, our findings provide the first evidence supporting the impact of universal DST on the transmission of RR-TB.

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.

The host-pathogen-environment triad: Lessons learned through the study of the multidrug-resistant Mycobacterium tuberculosis M strain

Multidrug-resistant tuberculosis is one of the major obstacles that face the tuberculosis eradication efforts. Drug-resistant Mycobacterium tuberculosis clones were initially disregarded as a public health threat, because they were assumed to have paid a high fitness cost in exchange of resistance acquisition. However, some genotypes manage to overcome the impact of drug-resistance conferring mutations, retain transmissibility and cause large outbreaks. In Argentina, the HIV-AIDS epidemics fuelled the expansion of the so-called M strain in the early 1990s, which is responsible for the largest recorded multidrug-resistant tuberculosis cluster of Latin America. The aim of this work is to review the knowledge gathered after nearly three decades of multidisciplinary research on epidemiological, microbiological and immunological aspects of this highly successful strain. Collectively, our results indicate that the successful transmission of the M strain could be ascribed to its unaltered virulence, low Th1/Th17 response, a low fitness cost imposed by the resistance conferring mutations and a high resistance to host-related stress. In the early 2000s, the incident cases due to the M strain steadily declined and stabilized in the latest years. Improvements in the management, diagnosis and treatment of multidrug-resistant tuberculosis along with societal factors such as the low domestic and international mobility of the patients affected by this strain probably contributed to the outbreak containment. This stresses the importance of sustaining the public health interventions to avoid the resurgence of this conspicuous multidrug-resistant strain.

Molecular Epidemiology of Mycobacterium tuberculosis Complex Strains in Urban and Slum Settings of Nairobi, Kenya

Kenya is a country with a high tuberculosis (TB) burden. However, knowledge on the genetic diversity of Mycobacterium tuberculosis complex (MTBC) strains and their transmission dynamics is sparsely available. Hence, we used whole-genome sequencing (WGS) to depict the genetic diversity, molecular markers of drug resistance, and possible transmission clusters among MTBC strains in urban and slum settings of Nairobi. We analyzed 385 clinical MTBC isolates collected between 2010 and 2015 in combination with patients’ demographics. We showed that the MTBC population mainly comprises strains of four lineages (L1–L4). The two dominating lineages were L4 with 55.8% (n = 215) and L3 with 25.7% (n = 99) of all strains, respectively. Genome-based cluster analysis showed that 30.4% (117/385) of the strains were clustered using a ≤5 single-nucleotide polymorphism (SNP) threshold as a surrogate marker for direct patient-to-patient MTBC transmission. Moreover, 5.2% (20/385) of the strains were multidrug-resistant (MDR), and 50.0% (n = 10) were part of a genome-based cluster (i.e., direct MDR MTBC transmission). Notably, 30.0% (6/20) of the MDR strains were resistant to all first-line drugs and are part of one molecular cluster. Moreover, TB patients in urban living setting had 3.8 times the odds of being infected with a drug-resistant strain as compared to patients from slums (p-value = 0.002). Our results show that L4 strains are the main causative agent of TB in Nairobi and MDR strain transmission is an emerging concern in urban settings. This emphasizes the need for more focused infection control measures and contact tracing of patients with MDR TB to break the transmission chains.

Phylogeography and transmission of M. tuberculosis in Moldova: A prospective genomic analysis

BACKGROUND

The incidence of multidrug-resistant tuberculosis (MDR-TB) remains critically high in countries of the former Soviet Union, where >20% of new cases and >50% of previously treated cases have resistance to rifampin and isoniazid. Transmission of resistant strains, as opposed to resistance selected through inadequate treatment of drug-susceptible tuberculosis (TB), is the main driver of incident MDR-TB in these countries.

METHODS AND FINDINGS

We conducted a prospective, genomic analysis of all culture-positive TB cases diagnosed in 2018 and 2019 in the Republic of Moldova. We used phylogenetic methods to identify putative transmission clusters; spatial and demographic data were analyzed to further describe local transmission of Mycobacterium tuberculosis. Of 2,236 participants, 779 (36%) had MDR-TB, of whom 386 (50%) had never been treated previously for TB. Moreover, 92% of multidrug-resistant M. tuberculosis strains belonged to putative transmission clusters. Phylogenetic reconstruction identified 3 large clades that were comprised nearly uniformly of MDR-TB: 2 of these clades were of Beijing lineage, and 1 of Ural lineage, and each had additional distinct clade-specific second-line drug resistance mutations and geographic distributions. Spatial and temporal proximity between pairs of cases within a cluster was associated with greater genomic similarity. Our study lasted for only 2 years, a relatively short duration compared with the natural history of TB, and, thus, the ability to infer the full extent of transmission is limited.

CONCLUSIONS

The MDR-TB epidemic in Moldova is associated with the local transmission of multiple M. tuberculosis strains, including distinct clades of highly drug-resistant M. tuberculosis with varying geographic distributions and drug resistance profiles. This study demonstrates the role of comprehensive genomic surveillance for understanding the transmission of M. tuberculosis and highlights the urgency of interventions to interrupt transmission of highly drug-resistant M. tuberculosis.

Alterations in molecular response of Mycobacterium tuberculosis against anti-tuberculosis drugs

BACKGROUND

Tuberculosis (TB), an infectious disease caused by Mycobacterium tuberculosis, has plagued humans since the early middle-ages. More than one million deaths are recorded annually due to TB, even in present times. These deaths are primarily attributed to the constant appearance of resistant TB strains. Even with the advent of new therapeutics and diagnostics techniques, tuberculosis remains challenging to control due to resistant M. tuberculosis strains. Aided by various molecular changes, these strains adapt to stress created by anti-tuberculosis drugs.

MATERIALS AND METHODS

The review thus is an overview of ongoing research in the genome and transcriptome of antibiotic-resistant TB. It explores omics-based research to identify mutation and utilization of differential gene expression.

CONCLUSIONS

This study shows several mutations distinctive in the first- and second-line drug-resistant M. tuberculosis strains. It also explores the expressional differences of genes involved in the fundamental process of the cells and how they help in drug resistance. With the development of transcriptomics-based studies, a new insight has developed to inquire about gene expression changes in drug resistance. This information on expressional pattern changes can be utilized to design the basic platform of anti-TB treatments and therapeutic approaches. These novel insights can be instrumental in disease diagnosis and global containment of resistant TB.

Connection between two historical tuberculosis outbreak sites in Japan, Honshu, by a new ancestral Mycobacterium tuberculosis L2 sublineage

By gathering 680 publicly available Sequence Read Archives from isolates of Mycobacterium tuberculosis complex (MTBC) including 190 belonging to the lineage 2 Beijing , and using an in-house bioinformatical pipeline, the TB-Annotator , that analyses more than 50 000 characters, we describe herein a new L2 sublineage from 20 isolates found in the Tochigi province, (Japan), that we designate as asia ancestral 5 (AAnc5). These isolates harbour a number of specific criteria (42 SNPs) and their intra-cluster pairwise distance suggests historical and not epidemiological transmission. These isolates harbour a mutation in rpoC , and do not fulfil, any of the modern Beijing lineage criteria, nor any of the other ancestral Beijing lineages described so far. Asia ancestral 5 isolates do not possess mutT2 58 and ogt 12 characteristics of modern Beijing , but possess ancestral Beijing SNPs characteristics. By looking into the literature, we found a reference isolate ID381, described in Kobe and Osaka belonging to the ‘G3’ group, sharing 36 out of the 42 specific SNPs found in AAnc5. We also assessed the intermediate position of the asia ancestral 4 (AAnc4) sublineage recently described in Thailand and propose an improved classification of the L2 that now includes AAnc4 and AAnc5. By increasing the recruitment into TB-Annotator to around 3000 genomes (including 642 belonging to L2), we confirmed our results and discovered additional historical ancestral L2 branches that remain to be investigated in more detail. We also present, in addition, some anthropological and historical data from Chinese and Japan history of tuberculosis, as well as from Korea, that could support our results on L2 evolution. This study shows that the reconstruction of the early history of tuberculosis in Asia is likely to reveal complex patterns since its emergence.

Revised nomenclature and SNP barcode for Mycobacterium tuberculosis lineage 2

Mycobacterium tuberculosis (Mtb) lineage 2 (L2) strains are present globally, contributing to a widespread tuberculosis (TB) burden, particularly in Asia where both prevalence of TB and numbers of drug resistant TB are highest. The increasing availability of whole-genome sequencing (WGS) data worldwide provides an opportunity to improve our understanding of the global genetic diversity of Mtb L2 and its association with the disease epidemiology and pathogenesis. However, existing L2 sublineage classification schemes leave >20 % of the Modern Beijing isolates unclassified. Here, we present a revised SNP-based classification scheme of L2 in a genomic framework based on phylogenetic analysis of >4000 L2 isolates from 34 countries in Asia, Eastern Europe, Oceania and Africa. Our scheme consists of over 30 genotypes, many of which have not been described before. In particular, we propose six main genotypes of Modern Beijing strains, denoted L2.2.M1–L2.2.M6. We also provide SNP markers for genotyping L2 strains from WGS data. This fine-scale genotyping scheme, which can classify >98 % of the studied isolates, serves as a basis for more effective monitoring and reporting of transmission and outbreaks, as well as improving genotype-phenotype associations such as disease severity and drug resistance. This article contains data hosted by Microreact.

Bacterial evolution during human infection: Adapt and live or adapt and die

Microbes are constantly evolving. Laboratory studies of bacterial evolution increase our understanding of evolutionary dynamics, identify adaptive changes, and answer important questions that impact human health. During bacterial infections in humans, however, the evolutionary parameters acting on infecting populations are likely to be much more complex than those that can be tested in the laboratory. Nonetheless, human infections can be thought of as naturally occurring in vivo bacterial evolution experiments, which can teach us about antibiotic resistance, pathogenesis, and transmission. Here, we review recent advances in the study of within-host bacterial evolution during human infection and discuss practical considerations for conducting such studies. We focus on 2 possible outcomes for de novo adaptive mutations, which we have termed "adapt-and-live" and "adapt-and-die." In the adapt-and-live scenario, a mutation is long lived, enabling its transmission on to other individuals, or the establishment of chronic infection. In the adapt-and-die scenario, a mutation is rapidly extinguished, either because it carries a substantial fitness cost, it arises within tissues that block transmission to new hosts, it is outcompeted by more fit clones, or the infection resolves. Adapt-and-die mutations can provide rich information about selection pressures in vivo, yet they can easily elude detection because they are short lived, may be more difficult to sample, or could be maladaptive in the long term. Understanding how bacteria adapt under each of these scenarios can reveal new insights about the basic biology of pathogenic microbes and could aid in the design of new translational approaches to combat bacterial infections.

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.

MDR and Pre-XDR Clinical Mycobacterium tuberculosis Beijing Strains: Assessment of Virulence and Host Cytokine Response in Mice Infectious Model

Mycobacterium tuberculosis Beijing genotype associated with drug resistance is a growing public health problem worldwide. The aim of this study was the assessment of virulence for C57BL/6 mice after infection by clinical M. tuberculosis strains 267/47 and 120/26, which belong to the modern sublineages B0/W148 and Central Asia outbreak of the Beijing genotype, respectively. The sublineages were identified by the analysis of the strains' whole-genomes. The strains 267/47 and 120/26 were characterized as agents of pre-extensively drug-resistant (pre-XDR) and multidrug-resistant (MDR) tuberculosis, respectively. Both clinical strains were slow-growing in 7H9 broth compared to the M. tuberculosis H37Rv strain. The survival rates of C57BL/6 mice infected by 267/47, 120/26, and H37Rv on the 150th day postinfection were 10%, 40%, and 70%, respectively. Mycobacterial load in the lungs, spleen, and liver was higher and histopathological changes were more expressed for mice infected by the 267/47 strain compared to those infected by the 120/26 and H37Rv strains. The cytokine response in the lungs of C57BL/6 mice after infection with the 267/47, 120/26, and H37Rv strains was different. Notably, proinflammatory cytokine genes Il-1α, Il-6, Il-7, and Il-17, as well as anti-inflammatory genes Il-6 and Il-13, were downregulated after an infection caused by the 267/47 strain compared to those after infection with the H37Rv strain.

The within-host evolution of antimicrobial resistance in Mycobacterium tuberculosis

Tuberculosis (TB) has been responsible for the greatest number of human deaths due to an infectious disease in general, and due to antimicrobial resistance (AMR) in particular. The etiological agents of human TB are a closely-related group of human-adapted bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Understanding how MTBC populations evolve within-host may allow for improved TB treatment and control strategies. In this review, we highlight recent works that have shed light on how AMR evolves in MTBC populations within individual patients. We discuss the role of heteroresistance in AMR evolution, and review the bacterial, patient and environmental factors that likely modulate the magnitude of heteroresistance within-host. We further highlight recent works on the dynamics of MTBC genetic diversity within-host, and discuss how spatial substructures in patients' lungs, spatiotemporal heterogeneity in antimicrobial concentrations and phenotypic drug tolerance likely modulates the dynamics of MTBC genetic diversity in patients during treatment. We note the general characteristics that are shared between how the MTBC and other bacterial pathogens evolve in humans, and highlight the characteristics unique to the MTBC.

Expression Dysregulation as a Mediator of Fitness Costs in Antibiotic Resistance

Antimicrobial resistance (AMR) poses a threat to global health and the economy. Rifampicin-resistant Mycobacterium tuberculosis accounts for a third of the global AMR burden. Gaining the upper hand on AMR requires a deeper understanding of the physiology of resistance. AMR often results in a fitness cost in the absence of drug. Identifying the molecular mechanisms underpinning this cost could help strengthen future treatment regimens. Here, we used a collection of M. tuberculosis strains that provide an evolutionary and phylogenetic snapshot of rifampicin resistance and subjected them to genome-wide transcriptomic and proteomic profiling to identify key perturbations of normal physiology. We found that the clinically most common rifampicin resistance-conferring mutation, RpoB Ser450Leu, imparts considerable gene expression changes, many of which are mitigated by the compensatory mutation in RpoC Leu516Pro. However, our data also provide evidence for pervasive epistasis—the same resistance mutation imposed a different fitness cost and functionally distinct changes to gene expression in genetically unrelated clinical strains. Finally, we report a likely posttranscriptional modulation of gene expression that is shared in most of the tested strains carrying RpoB Ser450Leu, resulting in an increased abundance of proteins involved in central carbon metabolism. These changes contribute to a more general trend in which the disruption of the composition of the proteome correlates with the fitness cost of the RpoB Ser450Leu mutation in different strains.

Compensatory effects of M. tuberculosis rpoB mutations outside the rifampicin resistance-determining region

Mycobacterium tuberculosis has been observed to develop resistance to the frontline anti-tuberculosis drug rifampicin, primarily through mutations in the rifampicin resistance-determining region (RRDR) of rpoB. While these mutations have been determined to confer a fitness cost, compensatory mutations in rpoA and rpoC that may enhance the fitness of resistant strains have been demonstrated. Recent genomic studies identified several rpoB non-RRDR mutations that co-occurred with RRDR mutations in clinical isolates without rpoA/rpoC mutations and may confer fitness compensation. In this study, we identified 33 evolutionarily convergent rpoB non-RRDR mutations through phylogenomic analysis of public genomic data for clinical M. tuberculosis isolates. We found that none of these mutations, except V170F and I491F, can cause rifampin resistance in Mycolicibacterium smegmatis. The compensatory effects of five representative mutations across rpoB were evaluated by an in vitro competition assay, through which we observed that each of these mutations can significantly improve the relative fitness of the initial S450L mutant (0.97–1.08 vs 0.87). Furthermore, we observed that the decreased RNAP transcription efficiency introduced by S450L was significantly alleviated by each of the five mutations. Structural analysis indicated that the fitness compensation observed for the non-RRDR mutations might be achieved by modification of the RpoB active centre or by changes in interactions between RNAP subunits. Our results provide experimental evidence supporting that compensatory effects are exerted by several rpoB non-RRDR mutations, which could be utilized as additional molecular markers for predicting the fitness of clinical rifampin-resistant M. tuberculosis strains.

The Epistatic Landscape of Antibiotic Resistance of Different Clades of Mycobacterium tuberculosis

Drug resistance (DR) remains a global challenge in tuberculosis (TB) control. In order to develop molecular-based diagnostic methods to replace the traditional culture-based diagnostics, there is a need for a thorough understanding of the processes that govern TB drug resistance. The use of whole-genome sequencing coupled with statistical and computational methods has shown great potential in unraveling the complexity of the evolution of DR-TB. In this study, we took an innovative approach that sought to determine nonrandom associations between polymorphic sites in Mycobacterium tuberculosis (Mtb) genomes. Attributable risk statistics were applied to identify the epistatic determinants of DR in different clades of Mtb and the possible evolutionary pathways of DR development. It was found that different lineages of Mtb exploited different evolutionary trajectories towards multidrug resistance and compensatory evolution to reduce the DR-associated fitness cost. Epistasis of DR acquisition is a new area of research that will aid in the better understanding of evolutionary biological processes and allow predicting upcoming multidrug-resistant pathogens before a new outbreak strikes humanity.

Prisons as ecological drivers of fitness-compensated multidrug-resistant Mycobacterium tuberculosis

Multidrug-resistant tuberculosis (MDR-TB) accounts for one third of the annual deaths due to antimicrobial resistance¹. Drug resistance-conferring mutations frequently cause fitness costs in bacteria^2-5. Experimental work indicates that these drug resistance-related fitness costs might be mitigated by compensatory mutations^6-10. However, the clinical relevance of compensatory evolution remains poorly understood. Here we show that, in the country of Georgia, during a 6-year nationwide study, 63% of MDR-TB was due to patient-to-patient transmission. Compensatory mutations and patient incarceration were independently associated with transmission. Furthermore, compensatory mutations were overrepresented among isolates from incarcerated individuals that also frequently spilled over into the non-incarcerated population. As a result, up to 31% of MDR-TB in Georgia was directly or indirectly linked to prisons. We conclude that prisons fuel the epidemic of MDR-TB in Georgia by acting as ecological drivers of fitness-compensated strains with high transmission potential.

Multiple Merger Genealogies in Outbreaks of Mycobacterium tuberculosis

The Kingman coalescent and its developments are often considered among the most important advances in population genetics of the last decades. Demographic inference based on coalescent theory has been used to reconstruct the population dynamics and evolutionary history of several species, including Mycobacterium tuberculosis (MTB), an important human pathogen causing tuberculosis. One key assumption of the Kingman coalescent is that the number of descendants of different individuals does not vary strongly, and violating this assumption could lead to severe biases caused by model misspecification. Individual lineages of MTB are expected to vary strongly in reproductive success because 1) MTB is potentially under constant selection due to the pressure of the host immune system and of antibiotic treatment, 2) MTB undergoes repeated population bottlenecks when it transmits from one host to the next, and 3) some hosts show much higher transmission rates compared with the average (superspreaders).

Here, we used an approximate Bayesian computation approach to test whether multiple-merger coalescents (MMC), a class of models that allow for large variation in reproductive success among lineages, are more appropriate models to study MTB populations. We considered 11 publicly available whole-genome sequence data sets sampled from local MTB populations and outbreaks and found that MMC had a better fit compared with the Kingman coalescent for 10 of the 11 data sets. These results indicate that the null model for analyzing MTB outbreaks should be reassessed and that past findings based on the Kingman coalescent need to be revisited.

Impact of bedaquiline on treatment outcomes of multidrug-resistant tuberculosis in a high-burden country

BACKGROUND

Evaluation of novel anti-tuberculosis (TB) drugs for the treatment of multidrug-resistant (MDR)-TB continues to be of high interest on the TB research agenda. We assessed treatment outcomes in patients with pulmonary MDR-TB who received bedaquiline-containing treatment regimens in the Republic of Moldova, a high-burden MDR-TB country.

METHOD

We systematically analysed the SIMETB national electronic TB database and performed a retrospective propensity score-matched comparison of treatment outcomes in a cohort of patients with MDR-TB who started treatment during 2016-2018 with a bedaquiline-containing regimen (bedaquiline cohort) and a cohort of patients treated without bedaquiline (non-bedaquiline cohort).

RESULTS

Following propensity score matching, 114 patients were assigned to each cohort of MDR-TB patients. Patients in the bedaquiline cohort had a higher 6-month sputum culture conversion rate than those in the non-bedaquiline cohort (66.7% versus 40.3%; p<0.001). Patients under bedaquiline-containing regimens had a higher cure rate assessed by both World Health Organization (WHO) and TBnet definitions (55.3% versus 24.6%; p=0.001 and 43.5% versus 19.6%; p=0.004, respectively), as well as a lower mortality rate (8.8% versus 20.2%; p<0.001 and 10.9% versus 25.2%; p=0.01, respectively). In patients who previously failed on MDR-TB treatment, >40% of patients achieved a cure with a bedaquiline-containing regimen.

CONCLUSIONS

Bedaquiline-based MDR-TB treatment regimens result in better disease resolution when compared with bedaquiline-sparing MDR-TB treatment regimens under programmatic conditions in a country with a high burden of MDR-TB.

Highly transmitted M. tuberculosis strains are more likely to evolve MDR/XDR and cause outbreaks, but what makes them highly transmitted?

Multi-Drug-Resistant strains of Mycobacterium tuberculosis (MDR-TB) are a serious obstacle to global TB eradication. While most MDR-TB strains are infrequently transmitted, a few cause large transmission clusters that contribute substantially to local MDR-TB burdens. Here we examine whether the known mutations in these strains can explain their success. Drug resistance mutations differ in fitness costs and strains can also acquire compensatory mutations (CM) to restore fitness, but some highly transmitted MDR strains have no CM. The acquisition of resistance mutations that maintain high transmissibility seems to occur by chance and are more likely in strains that are intrinsically highly transmitted and cause many cases. Modern Beijing lineage strains have caused several large outbreaks, but MDR outbreaks are also caused by ancient Beijing and lineage 4 strains, suggesting the lineage is less important than the characteristics of the individual strain. The development of fluoroquinolone resistance appears to represent another level of selection, in which strains must surmount unknown fitness costs of gyrA mutations. The genetic determinants of high transmission are poorly defined but may involve genes encoding proteins involved in molybdenum acquisition and the Esx systems. In addition, strains eliciting lower cytokine responses and producing more caseating granulomas may have advantages for transmission. Successful MDR/XDR strains generally evolve from highly transmitted drug sensitive parent strains due to selection pressures from deficiencies in local TB control programs. Until TB incidence is considerably reduced, there will likely be highly transmitted strains that develop resistance to any new antibiotic.

Quantifying transmission fitness costs of multi-drug resistant tuberculosis

As multi-drug resistant tuberculosis (MDR-TB) continues to spread, investigating the transmission potential of different drug-resistant strains becomes an ever more pressing topic in public health. While phylogenetic and transmission tree inferences provide valuable insight into possible transmission chains, phylodynamic inference combines evolutionary and epidemiological analyses to estimate the parameters of the underlying epidemiological processes, allowing us to describe the overall dynamics of disease spread in the population. In this study, we introduce an approach to Mycobacterium tuberculosis (M. tuberculosis) phylodynamic analysis employing an existing computationally efficient model to quantify the transmission fitness costs of drug resistance with respect to drug-sensitive strains. To determine the accuracy and precision of our approach, we first perform a simulation study, mimicking the simultaneous spread of drug-sensitive and drug-resistant tuberculosis (TB) strains. We analyse the simulated transmission trees using the phylodynamic multi-type birth-death model (MTBD, (Kühnert et al., 2016)) within the BEAST2 framework and show that this model can estimate the parameters of the epidemic well, despite the simplifying assumptions that MTBD makes compared to the complex TB transmission dynamics used for simulation. We then apply the MTBD model to an M. tuberculosis lineage 4 dataset that primarily consists of MDR sequences. Some of the MDR strains additionally exhibit resistance to pyrazinamide - an important first-line anti-tuberculosis drug. Our results support the previously proposed hypothesis that pyrazinamide resistance confers a transmission fitness cost to the bacterium, which we quantify for the given dataset. Importantly, our sensitivity analyses show that the estimates are robust to different prior distributions on the resistance acquisition rate, but are affected by the size of the dataset - i.e. we estimate a higher fitness cost when using fewer sequences for analysis. Overall, we propose that MTBD can be used to quantify the transmission fitness cost for a wide range of pathogens where the strains can be appropriately divided into two or more categories with distinct properties.

Thin-Layer-Agar-Based Direct Phenotypic Drug Susceptibility Testing on Sputum in Eswatini Rapidly Detects Mycobacterium tuberculosis Growth and Rifampicin Resistance Otherwise Missed by WHO-Endorsed Diagnostic Tests

Xpert MTB/RIF rapidly detects resistance to rifampicin (RR); however, this test misses I491F-RR conferring rpoB mutation, common in southern Africa. In addition, Xpert MTB/RIF does not distinguish between viable and dead Mycobacterium tuberculosis (MTB). We aimed to investigate the ability of thin-layer agar (TLA) direct drug-susceptibility testing (DST) to detect MTB and its drug-resistance profiles in field conditions in Eswatini. Consecutive samples were tested in parallel with Xpert MTB/RIF and TLA for rifampicin (1.0 μg/ml) and ofloxacin (2.0 μg/ml). TLA results were compared at the Reference Laboratory in Antwerp with indirect-DST on Löwenstein-Jensen or 7H11 solid media and additional phenotypic and genotypic testing to resolve discordance. TLA showed a positivity rate for MTB detection of 7.1% versus 10.0% for Xpert MTB/RIF. Of a total of 4,547 samples included in the study, 200 isolates were available for comparison to the composite reference. Within a median of 18.4 days, TLA detected RR with 93.0% sensitivity (95% confidence interval [CI], 77.4 to 98.0) and 99.4% specificity (95% CI, 96.7 to 99.9) versus 62.5% (95% CI, 42.7 to 78.8) and 99.3% (95% CI, 96.2 to 99.9) for Xpert MTB/RIF. Eight isolates, 28.6% of all RR-confirmed isolates, carried the I491F mutation, all detected by TLA. TLA also correctly identified 183 of the 184 ofloxacin-susceptible isolates (99.5% specificity; 95% CI, 97.0 to 99.9). In field conditions, TLA rapidly detects RR, and in this specific setting, it contributed to detection of additional RR patients over Xpert MTB/RIF, mainly but not exclusively due to I491F. TLA also accurately excluded fluoroquinolone resistance.