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

Differential rates of Mycobacterium tuberculosis transmission associate with host-pathogen sympatry

Several human-adapted Mycobacterium tuberculosis complex (Mtbc) lineages exhibit a restricted geographical distribution globally. These lineages are hypothesized to transmit more effectively among sympatric hosts, that is, those that share the same geographical area, though this is yet to be confirmed while controlling for exposure, social networks and disease risk after exposure. Using pathogen genomic and contact tracing data from 2,279 tuberculosis cases linked to 12,749 contacts from three low-incidence cities, we show that geographically restricted Mtbc lineages were less transmissible than lineages that have a widespread global distribution. Allopatric host-pathogen exposure, in which the restricted pathogen and host are from non-overlapping areas, had a 38% decrease in the odds of infection among contacts compared with sympatric exposures. We measure tenfold lower uptake of geographically restricted lineage 6 strains compared with widespread lineage 4 strains in allopatric macrophage infections. We conclude that Mtbc strain-human long-term coexistence has resulted in differential transmissibility of Mtbc lineages and that this differs by human population.

Spread of antibiotic resistance genes in drinking water reservoirs: Insights from a deep metagenomic study using a curated database

The exploration of antibiotic resistance genes (ARGs) in drinking water reservoirs is an emerging field. Using a curated database, we enhanced the ARG detection and conducted a comprehensive analysis using 2.2 Tb of deep metagenomic sequencing data to determine the distribution of ARGs across 16 drinking water reservoirs and associated environments. Our findings reveal a greater diversity of ARGs in sediments than in water, underscoring the importance of extensive background surveys. Crucial ARG carriers-specifically Acinetobacter, Pseudomonas, and Mycobacterium were identified in drinking water reservoirs. Extensive analysis of the data uncovered a considerable concern for drinking water safety, particularly in regions reliant on river sources. Mobile genetic elements have been found to contribute markedly to the propagation of ARGs. The results of this research suggest that the establishment of drinking water reservoirs for supplying raw water may be an effective strategy for alleviating the spread of water-mediated ARGs.

Large-scale statistical analysis of Mycobacterium tuberculosis genome sequences identifies compensatory mutations associated with multi-drug resistance

Tuberculosis (TB), caused by Mycobacterium tuberculosis, has a significant impact on global health worldwide. The development of multi-drug resistant strains that are resistant to the first-line drugs isoniazid and rifampicin threatens public health security. Rifampicin and isoniazid resistance are largely underpinned by mutations in rpoB and katG respectively and are associated with fitness costs. Compensatory mutations are considered to alleviate these fitness costs and have been observed in rpoC/rpoA (rifampicin) and oxyR'-ahpC (isoniazid). We developed a framework (CompMut-TB) to detect compensatory mutations from whole genome sequences from a large dataset comprised of 18,396 M. tuberculosis samples. We performed association analysis (Fisher's exact tests) to identify pairs of mutations that are associated with drug-resistance, followed by mediation analysis to identify complementary or full mediators of drug-resistance. The analyses revealed several potential mutations in rpoC (N = 47), rpoA (N = 4), and oxyR'-ahpC (N = 7) that were considered either 'highly likely' or 'likely' to confer compensatory effects on drug-resistance, including mutations that have previously been reported and validated. Overall, we have developed the CompMut-TB framework which can assist with identifying compensatory mutations which is important for more precise genome-based profiling of drug-resistant TB strains and to further understanding of the evolutionary mechanisms that underpin drug-resistance.

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.

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.

On the onset and dispersal of a major MDR TB clone among HIV-negative patients, Tunisia

BACKGROUND

To carry out a whole genome sequencing (WGS)-based investigation on the emergence and spread of the largest multidrug-resistant tuberculosis (MDR TB) outbreak that has been thriving among HIV-negative patients, Tunisia, since the early 2000s.

METHODS

We performed phylogeographic analyses and molecular dating based on a WGS dataset representing 68 unique Mycobacterium tuberculosis isolates, covering almost the entire MDR TB outbreak for the time period 2001-2016.

RESULTS

The data indicate that the ancestor of the MDR TB outbreak emerged in the region of Bizerte, as early as 1974 (95% CI 1951-1985), from where it spread to other regions by 1992 (95% CI 1980-1996). Analysis of a minimum spanning tree based on core genome Multilocus Sequence Typing (cgMLST) uncovered the early spill-over of the fitness-compensated MDR TB strain from the prison into the general population. Indeed, cases with history of incarceration were found to be directly or indirectly linked to up to 22 new outbreak cases (32.35%) among the non-imprisoned population. By around 2008, the MDR TB outbreak strain had acquired additional resistance, leading to an XDR phenotype.

CONCLUSIONS

WGS allowed refining our understanding of the emergence and evolution of the largest MDR TB outbreak in Tunisia, whose causative strain has been circulating silently for almost 26 years before. Our study lends further support to the critical role of prisons-related cases in the early spread of the outbreak among the general population. The shift to an XDR phenotype of such an epidemic clone prompts an urgent need to undertake drastic control measures.

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.

Signatures of transmission in within-host M. tuberculosis variation

Background

Because M. tuberculosis evolves slowly, transmission clusters often contain multiple individuals with identical consensus genomes, making it difficult to reconstruct transmission chains. Finding additional sources of shared M. tuberculosis variation could help overcome this problem. Previous studies have reported M. tuberculosis diversity within infected individuals; however, whether within-host variation improves transmission inferences remains unclear.

Methods

To evaluate the transmission information present in within-host M. tuberculosis variation, we re-analyzed publicly available sequence data from three household transmission studies, using household membership as a proxy for transmission linkage between donor-recipient pairs.

Findings

We found moderate levels of minority variation present in M. tuberculosis sequence data from cultured isolates that varied significantly across studies (mean: 6, 7, and 170 minority variants above a 1% minor allele frequency threshold, outside of PE/PPE genes). Isolates from household members shared more minority variants than did isolates from unlinked individuals in the three studies (mean 98 shared minority variants vs. 10; 0.8 vs. 0.2, and 0.7 vs. 0.2, respectively). Shared within-host variation was significantly associated with household membership (OR: 1.51 [1.30,1.71], for one standard deviation increase in shared minority variants). Models that included shared within-host variation improved the accuracy of predicting household membership in all three studies as compared to models without within-host variation (AUC: 0.95 versus 0.92, 0.99 versus 0.95, and 0.93 versus 0.91).

Interpretation

Within-host M. tuberculosis variation persists through culture and could enhance the resolution of transmission inferences. The substantial differences in minority variation recovered across studies highlights the need to optimize approaches to recover and incorporate within-host variation into automated phylogenetic and transmission inference.

Funding

NIAID: 5K01AI173385.

Regulatory fine-tuning of mcr-1 increases bacterial fitness and stabilises antibiotic resistance in agricultural settings

Antibiotic resistance tends to carry fitness costs, making it difficult to understand how resistance can be maintained in the absence of continual antibiotic exposure. Here we investigate this problem in the context of mcr-1, a globally disseminated gene that confers resistance to colistin, an agricultural antibiotic that is used as a last resort for the treatment of multi-drug resistant infections. Here we show that regulatory evolution has fine-tuned the expression of mcr-1, allowing E. coli to reduce the fitness cost of mcr-1 while simultaneously increasing colistin resistance. Conjugative plasmids have transferred low-cost/high-resistance mcr-1 alleles across an incredible diversity of E. coli strains, further stabilising mcr-1 at the species level. Regulatory mutations were associated with increased mcr-1 stability in pig farms following a ban on the use of colistin as a growth promoter that decreased colistin consumption by 90%. Our study shows how regulatory evolution and plasmid transfer can combine to stabilise resistance and limit the impact of reducing antibiotic consumption.

Prioritizing persons deprived of liberty in global guidelines for tuberculosis preventive treatment

In this Policy Forum piece, Aditya Narayan and colleagues discuss the challenges and opportunities for tuberculosis preventive treatment in carceral settings.

Insight into the drug-resistant characteristics and genetic diversity of multidrug-resistant Mycobacterium tuberculosis in China

Multidrug-resistant tuberculosis (MDR-TB) has a severe impact on public health. To investigate the drug-resistant profile, compensatory mutations and genetic variations among MDR-TB isolates, a total of 546 MDR-TB isolates from China underwent drug-susceptibility testing and whole genome sequencing for further analysis. The results showed that our isolates have a high rate of fluoroquinolone resistance (45.60%, 249/546) and a low proportion of conferring resistance to bedaquiline, clofazimine, linezolid, and delamanid. The majority of MDR-TB isolates (77.66%, 424/546) belong to Lineage 2.2.1, followed by Lineage 4.5 (6.41%, 35/546), and the Lineage 2 isolates have a strong association with pre-XDR/XDR-TB (P < 0.05) in our study. Epidemic success analysis using time-scaled haplotypic density (THD) showed that clustered isolates outperformed non-clustered isolates. Compensatory mutations happened in rpoA, rpoC, and non-RRDR of rpoB genes, which were found more frequently in clusters and were associated with the increase of THD index, suggesting that increased bacterial fitness was associated with MDR-TB transmission. In addition, the variants in resistance associated genes in MDR isolates are mainly focused on single nucleotide polymorphism mutations, and only a few genes have indel variants, such as katG, ethA. We also found some genes underwent indel variation correlated with the lineage and sub-lineage of isolates, suggesting the selective evolution of different lineage isolates. Thus, this analysis of the characterization and genetic diversity of MDR isolates would be helpful in developing effective strategies for treatment regimens and tailoring public interventions. IMPORTANCE Multidrug-resistant tuberculosis (MDR-TB) is a serious obstacle to tuberculosis prevention and control in China. This study provides insight into the drug-resistant characteristics of MDR combined with phenotypic drug-susceptibility testing and whole genome sequencing. The compensatory mutations and epidemic success analysis were analyzed by time-scaled haplotypic density (THD) method, suggesting clustered isolates and compensatory mutations are associated with MDR-TB transmission. In addition, the insertion and deletion variants happened in some genes, which are associated with the lineage and sub-lineage of isolates, such as the mpt64 gene. This study offered a valuable reference and increased understanding of MDR-TB in China, which could be crucial for achieving the objective of precision medicine in the prevention and treatment of MDR-TB.

Genotypic and spatial analysis of transmission dynamics of tuberculosis in Shanghai, China: a 10-year prospective population-based surveillance study

Background

With improved tuberculosis (TB) control programs, the incidence of TB in China declined dramatically over the past few decades, but recently the rate of decrease has slowed, especially in large cities such as Shanghai. To help formulate strategies to further reduce TB incidence, we performed a 10-year study in Songjiang, a district of Shanghai, to delineate the characteristics, transmission patterns, and dynamic changes of the local TB burden.

Methods

We conducted a population-based study of culture-positive pulmonary TB patients diagnosed in Songjiang during 2011-2020. Genomic clusters were defined with a threshold distance of 12-single-nucleotide-polymorphisms based on whole-genome sequencing, and risk factors for clustering were identified by logistic regression. Transmission inference was performed using phybreak. The distances between the residences of patients were compared to the genomic distances of their isolates. Spatial patient hotspots were defined with kernel density estimation.

Findings

Of 2212 enrolled patients, 74.7% (1652/2212) were internal migrants. The clustering rate (25.2%, 558/2212) and spatial concentrations of clustered and unclustered patients were unchanged over the study period. Migrants had significantly higher TB rates but less clustering than residents. Clustering was highest in male migrants, younger patients and both residents and migrants employed in physical labor. Only 22.1% of transmission events occurred between residents and migrants, with residents more likely to transmit to migrants. The clustering risk decreased rapidly with increasing distances between patient residences, but more than half of clustered patient pairs lived ≥5 km apart. Epidemiologic links were identified for only 15.6% of clustered patients, mostly in close contacts.

Interpretation

Although some of the TB in Songjiang's migrant population is caused by strains brought by infected migrants, local, recent transmission is an important driver of the TB burden. These results suggest that further reductions in TB incidence require novel strategies to detect TB early and interrupt urban transmission.

Funding

Shanghai Municipal Science and Technology Major Project (ZD2021CY001), National Natural Science Foundation of China (82272376), National Research Council of Science and Technology Major Project of China (2017ZX10201302-006).

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

Scaling up evidence-based approaches to tuberculosis screening in prisons

People deprived of liberty have among the highest rates of tuberculosis globally. The incidence of tuberculosis is ten times greater than the incidence of tuberculosis in the general population. In 2021, WHO updated its guidance to strongly recommend systematic screening for tuberculosis in prisons and penitentiary systems. Which case-finding strategies should be adopted, and how to effectively implement these strategies in these settings, will be crucial questions facing ministries of health and justice. In this Viewpoint, we review the evidence base for tuberculosis screening and diagnostic strategies in prisons, highlighting promising approaches and knowledge gaps. Drawing upon past experiences of implementing active case-finding and care programmes in settings with a high tuberculosis burden, we discuss challenges and opportunities for improving the tuberculosis diagnosis and treatment cascade in these settings. We argue that improved transparency in reporting of tuberculosis notifications and outcomes in prisons and renewed focus and resourcing from WHO and other stakeholders will be crucial for building the commitment and investments needed from countries to address the continued crisis of tuberculosis in prisons.

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.

Discrepancy in the transmissibility of multidrug-resistant Mycobacterium tuberculosis in urban and rural areas in China

The fitness of multidrug-resistant tuberculosis (MDR-TB) is thought to be an important determinant of a strain's ability to be transmitted and cause outbreaks. Studies in the laboratory have demonstrated that MDR-TB strains have reduced fitness but the relative transmissibility of MDR-TB versus drug-susceptible (DS) TB strains in human populations remains unresolved. We used data on genomic clustering from our previous molecular epidemiological study in Songjiang (2011-2020) and Wusheng (2009-2020), China, to compare the relative transmissibility of MDR-TB versus DS-TB. Genomic clusters were defined with a threshold distance of 12-single-nucleotide-polymorphisms and the risk for MDR-TB clustering was analyzed by logistic regression. In total, 2212 culture-positive pulmonary TB patients were enrolled in Songjiang and 1289 in Wusheng. The clustering rates of MDR-TB and DS-TB strains were 19.4% (20/103) and 26.3% (509/1936), respectively in Songjiang, and 43.9% (29/66) and 26.0% (293/1128) in Wusheng. The risk of MDR-TB clustering was 2.34 (95% CI 1.38-3.94) times higher than DS-TB clustering in Wusheng and 0.64 (95% CI 0.38-1.06) times lower in Songjiang. Neither lineage 2, compensatory mutations nor rpoB S450L were significantly associated with MDR-TB transmission, and katG S315T increased MDR-TB transmission only in Wusheng (OR 5.28, 95% CI 1.42-19.21). MDR-TB was not more transmissible than DS-TB in either Songjiang or Wusheng. It appears that the different transmissibility of MDR-TB in Songjiang and Wusheng is likely due to differences in the quality of the local TB control programs. These results suggest that the most effective way to control MDR-TB is by improving local TB control programs.

The Evolution and Transmission Dynamics of Multidrug-Resistant Tuberculosis in an Isolated High-Plateau Population of Tibet, China

On the Tibetan Plateau, most tuberculosis is caused by indigenous Mycobacterium tuberculosis strains with a monophyletic structure and high-level drug resistance. This study investigated the emergence, evolution, and transmission dynamics of multidrug-resistant tuberculosis (MDR-TB) in Tibet. The whole-genome sequences of 576 clinical strains from Tibet were analyzed with the TB-profiler tool to identify drug-resistance mutations. The evolution of the drug resistance was then inferred based on maximum-likelihood phylogeny and dated trees that traced the serial acquisition of mutations conferring resistance to different drugs. Among the 576 clinical M. tuberculosis strains, 346 (60.1%) carried at least 1 resistance-conferring mutation and 231 (40.1%) were MDR-TB. Using a pairwise distance of 50 single nucleotide polymorphisms (SNPs), most strains (89.9%, 518/576) were phylogenetically separated into 50 long-term transmission clusters. Eleven large drug-resistant clusters contained 76.1% (176/231) of the local multidrug-resistant strains. A total of 85.2% of the isoniazid-resistant strains were highly transmitted with an average of 6.6 cases per cluster, of which most shared the mutation KatG Ser315Thr. A lower proportion (71.6%) of multidrug-resistant strains were transmitted, with an average cluster size of 2.9 cases. The isoniazid-resistant clusters appear to have undergone substantial bacterial population growth in the 1970s to 1990s and then subsequently accumulated multiple rifampicin-resistance mutations and caused the current local MDR-TB burden. These findings highlight the importance of detecting and curing isoniazid-resistant strains to prevent the emergence of endemic MDR-TB. IMPORTANCE Emerging isoniazid resistance in the 1970s allowed M. tuberculosis strains to spread and form into large multidrug-resistant tuberculosis clusters in the isolated plateau of Tibet, China. The epidemic was driven by the high risk of transmission as well as the potential of acquiring further drug resistance from isoniazid-resistant strains. Eleven large drug-resistant clusters consisted of the majority of local multidrug-resistant cases. Among the clusters, isoniazid resistance overwhelmingly evolved before all the other resistance types. A large bacterial population growth of isoniazid-resistant clusters occurred between 1970s and 1990s, which subsequently accumulated rifampicin-resistance-conferring mutations in parallel and accounted for the local multidrug-resistant tuberculosis burden. The results of our study indicate that it may be possible to restrict MDR-TB evolution and dissemination by prioritizing screening for isoniazid (INH)-resistant TB strains before they become MDR-TB and by adopting measures that can limit their transmission.

New insights on tuberculosis transmission dynamics and drug susceptibility profiles among the prison population in Southern Brazil based on whole-genome sequencing

BACKGROUND

The rate of tuberculosis (TB) infection among the prison population (PP) in Brazil is 28 times higher than that in the general population, and prison environment favors the spread of TB.

OBJECTIVE

To describe TB transmission dynamics and drug resistance profiles in PP using whole-genome sequencing (WGS).

METHODS

This was a retrospective study of Mycobacterium tuberculosis cultivated from people incarcerated in 55 prisons between 2016 and 2019; only one isolate per prisoner was included. Information about movement from one prison to another was tracked. Clinical information was collected, and WGS was performed on isolates obtained at the time of TB diagnosis.

RESULTS

Among 134 prisoners included in the study, we detected 16 clusters with a total of 58 (43%) cases of M. tuberculosis. Clusters ranged from two to seven isolates with five or fewer single nucleotide polymorphism (SNP) differences, suggesting a recent transmission. Six (4.4%) isolates were resistant to at least one anti-TB drug. Two of these clustered together and showed resistance to rifampicin, isoniazid, and fluoroquinolones, with 100% concordance between WGS and phenotypic drug-susceptibility testing. Prisoners with clustered isolates had a high amount of movement between prisons (two to eight moves) during the study period.

CONCLUSIONS

WGS demonstrated the recent transmission of TB within prisons in Brazil. The high movement among prisoners seems to be related to the transmission of the same M. tuberculosis strain within the prison system. Screening for TB before and after the movement of prisoners using rapid molecular tests could play a role in reducing transmission.

Phylogeography and transmission of Mycobacterium tuberculosis spanning prisons and surrounding communities in Paraguay

Recent rises in incident tuberculosis (TB) cases in Paraguay and the increasing concentration of TB within prisons highlight the urgency of targeting strategies to interrupt transmission and prevent new infections. However, whether specific cities or carceral institutions play a disproportionate role in transmission remains unknown. We conducted prospective genomic surveillance, sequencing 471 Mycobacterium tuberculosis complex genomes, from inside and outside prisons in Paraguay's two largest urban areas, Asunción and Ciudad del Este, from 2016 to 2021. We found genomic evidence of frequent recent transmission within prisons and transmission linkages spanning prisons and surrounding populations. We identified a signal of frequent M. tuberculosis spread between urban areas and marked recent population size expansion of the three largest genomic transmission clusters. Together, our findings highlight the urgency of strengthening TB control programs to reduce transmission risk within prisons in Paraguay, where incidence was 70 times that outside prisons in 2021.

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.

Population-based sequencing of Mycobacterium tuberculosis reveals how current population dynamics are shaped by past epidemics

Transmission is a driver of tuberculosis (TB) epidemics in high-burden regions, with assumed negligible impact in low-burden areas. However, we still lack a full characterization of transmission dynamics in settings with similar and different burdens. Genomic epidemiology can greatly help to quantify transmission, but the lack of whole genome sequencing population-based studies has hampered its application. Here, we generate a population-based dataset from Valencia region and compare it with available datasets from different TB-burden settings to reveal transmission dynamics heterogeneity and its public health implications. We sequenced the whole genome of 785 Mycobacterium tuberculosis strains and linked genomes to patient epidemiological data. We use a pairwise distance clustering approach and phylodynamic methods to characterize transmission events over the last 150 years, in different TB-burden regions. Our results underscore significant differences in transmission between low-burden TB settings, i.e., clustering in Valencia region is higher (47.4%) than in Oxfordshire (27%), and similar to a high-burden area as Malawi (49.8%). By modeling times of the transmission links, we observed that settings with high transmission rate are associated with decades of uninterrupted transmission, irrespective of burden. Together, our results reveal that burden and transmission are not necessarily linked due to the role of past epidemics in the ongoing TB incidence, and highlight the need for in-depth characterization of transmission dynamics and specifically tailored TB control strategies.

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.

The role of prisons in disseminating tuberculosis in Brazil: A genomic epidemiology study

Background

Globally, prisons are high-incidence settings for tuberculosis. Yet the role of prisons as reservoirs of M. tuberculosis, propagating epidemics through spillover to surrounding communities, has been difficult to measure directly.

Methods

To quantify the role of prisons in driving wider community M. tuberculosis transmission, we conducted prospective genomic surveillance in Central West Brazil from 2014 to 2019. We whole genome sequenced 1152 M. tuberculosis isolates collected during active and passive surveillance inside and outside prisons and linked genomes to detailed incarceration histories. We applied multiple phylogenetic and genomic clustering approaches and inferred timed transmission trees.

Findings

M. tuberculosis sequences from incarcerated and non-incarcerated people were closely related in a maximum likelihood phylogeny. The majority (70.8%; 46/65) of genomic clusters including people with no incarceration history also included individuals with a recent history of incarceration. Among cases in individuals with no incarceration history, 50.6% (162/320) were in clusters that included individuals with recent incarceration history, suggesting that transmission chains often span prisons and communities. We identified a minimum of 18 highly probable spillover events, M. tuberculosis transmission from people with a recent incarceration history to people with no prior history of incarceration, occurring in the state’s four largest cities and across sampling years. We additionally found that frequent transfers of people between the state’s prisons creates a highly connected prison network that likely disseminates M. tuberculosis across the state.

Interpretation

We developed a framework for measuring spillover from high-incidence environments to surrounding communities by integrating genomic and spatial information. Our findings indicate that, in this setting, prisons serve not only as disease reservoirs, but also disseminate M. tuberculosis across highly connected prison networks, both amplifying and propagating M. tuberculosis risk in surrounding communities.

Funding

Brazil’s National Council for Scientific and Technological Development and US National Institutes of Health.

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.

Collateral Sensitivity Interactions between Antibiotics Depend on Local Abiotic Conditions

Mutations conferring resistance to one antibiotic can increase (cross-resistance) or decrease (collateral sensitivity) resistance to others. Antibiotic combinations displaying collateral sensitivity could be used in treatments that slow resistance evolution. However, lab-to-clinic translation requires understanding whether collateral effects are robust across different environmental conditions. Here, we isolated and characterized resistant mutants of Escherichia coli using five antibiotics, before measuring collateral effects on resistance to other paired antibiotics. During both isolation and phenotyping, we varied conditions in ways relevant in nature (pH, temperature, and bile). This revealed that local abiotic conditions modified expression of resistance against both the antibiotic used during isolation and other antibiotics. Consequently, local conditions influenced collateral sensitivity in two ways: by favoring different sets of mutants (with different collateral sensitivities) and by modifying expression of collateral effects for individual mutants. These results place collateral sensitivity in the context of environmental variation, with important implications for translation to real-world applications. IMPORTANCE When bacteria become resistant to an antibiotic, the genetic changes involved sometimes increase (cross-resistance) or decrease (collateral sensitivity) their resistance to other antibiotics. Antibiotic combinations showing repeatable collateral sensitivity could be used in treatment to slow resistance evolution. However, collateral sensitivity interactions may depend on the local environmental conditions that bacteria experience, potentially reducing repeatability and clinical application. Here, we show that variation in local conditions (pH, temperature, and bile salts) can influence collateral sensitivity in two ways: by favoring different sets of mutants during bacterial resistance evolution (with different collateral sensitivities to other antibiotics) and by modifying expression of collateral effects for individual mutants. This suggests that translation from the lab to the clinic of new approaches exploiting collateral sensitivity will be influenced by local abiotic conditions.

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.