- Department of Microbiology, University of Washington, Seattle, WA, United States
Tuberculosis, caused by Mycobacterium tuberculosis (Mtb) is an ancient disease that has remained a leading cause of infectious death. Mtb has evolved drug resistance to every antibiotic regimen ever introduced, greatly complicating treatment, lowering rates of cure and menacing TB control in parts of the world. As technology has advanced, our understanding of antimicrobial resistance has improved, and our models of the phenomenon have evolved. In this review, we focus on recent research progress that supports an updated model for the evolution of drug resistance in Mtb. We highlight the contribution of drug tolerance on the path to resistance, and the influence of heterogeneity on tolerance. Resistance is likely to remain an issue for as long as drugs are needed to treat TB. However, with technology driving new insights and careful management of newly developed resources, antimicrobial resistance need not continue to threaten global progress against TB, as it has done for decades.
Introduction
Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), has remained a centrally important cause of morbidity and mortality for centuries, but our understanding of TB disease, what causes it and how to combat it has evolved substantially over that time. In some of the earliest written records of the disease, the ancient Greeks recognized that TB was especially deadly, with the venerable “Father of Medicine” Hippocrates warning other physicians against treating advanced cases because the inevitable bad outcomes would damage the doctor’s reputation (Herzog, 1998). Eventually, the 19th century discovery of the TB bacillus and the 20th century introduction of effective chemotherapies seemed to promise a new era in which TB was tamed if not eliminated (Daniel, 2006; Barberis et al., 2017). However, the emergence of drug-resistant isolates was noted in the very first TB chemotherapy trials, and resistance has appeared whenever a new anti-TB agent is introduced (Crofton and Mitchison, 1948; Gillespie, 2002). With an estimated 1.3 million deaths in 2020 (WHO, 2021b), and given the dramatic change in world population, TB may claim as many total lives today as in the years before TB chemotherapy was available. This dismal situation has several causes, including co-morbidities like HIV and diabetes, and gaps in timely diagnosis and treatment, but drug resistance stands widely recognized as one of the major challenges to effective TB control worldwide.
The classical model of drug treatment and resistance
Just as ideas about TB have changed, thinking about drug resistance has undergone a significant evolution over time. Antibiotics were initially hailed as “magic bullets”, capable of stopping even lethal infectious diseases in their tracks, but the emergence of resistance-fueled treatment failures led to deeper investigations into the biology of drug response. Driven by powerful advances in genetics and the emergence of molecular biology in the latter half of the 20th century, a concise model of antibiotic action and resistance developed over decades. Briefly, this model proposes that antibiotics work by inhibiting some essential target, generally an enzyme, in the pathogen. Mutations occur at random and exist in each population prior to antibiotic pressure. Resistance emerges when pre-existing mutations promote growth or survival in the presence of the drug (Figure 1, left side). A corollary of this model is that resistance is a numbers game. Any pathogen population of sufficient size will harbor at least one mutation conferring resistance to each agent that can be selected by drug exposure, so adding additional drugs to a regimen serves to reduce the rate at which resistance emerges.
Figure 1 Visual representation of the classical and updated models for development of antibiotic resistance in Mtb. The top portion shows a population of Mtb in the absence of drug selection (- DRUG). In the classical model, pre-existing resistant bacilli (red) in the population grow despite the addition of drug (+ DRUG) and susceptible bacteria (gray) are eliminated over time until only resistant cells remain. In the updated model, pre-existing resistance may still be present, but tolerant bacilli (yellow) may also be present and/or be induced by drug exposure. With drug selection, susceptible bacteria are again eventually eliminated, with resistance dominating the population over time. Figure created with BioRender.com.
Grounded in Darwinian natural selection, this model has long provided a useful framework for exploring antibiotic action and resistance. However, the available data have never all fit comfortably within this paradigm. For example, while current treatment for drug-susceptible TB is ~95% effective, most treatment failures are due to relapses with drug sensitive bacilli (Colangeli et al., 2018). These relapses are hard to explain with the classical model of drug treatment and resistance, since any drug-susceptible bacterium exposed to antibiotics should be eliminated. This phenomenon is not unique to TB – studies with many different bacteria have long revealed curious phenotypes of enhanced survival in the absence of any heritable resistance that are not readily explained by the classical model (Bigger, 1944; Levin and Rozen, 2006). Further, the rate of spontaneous resistance to a three-drug combination has been estimated at greater than 1 in 1018 bacteria (Gillespie, 2007), many orders of magnitude greater than the number of TB bacilli in any one patient, possibly more than in all current patients combined. Even considering how some drugs are excluded from some niches in vivo (Strydom et al., 2019), it is difficult to reconcile the not-infrequent emergence of widespread resistance to 3 or more agents using the classical model alone. Now, driven again by advances in molecular genetics and especially sequencing, thinking about the evolution of resistance in TB has itself evolved. Here we explore how new observations and insights are promoting an updated model of TB drug resistance, with the promise of new ways to combat the resistance problem.
Classical resistance determinants in Mtb
Antibiotic resistance is defined as a heritable trait that enables a bacterial population to both survive and replicate in the presence of an otherwise inhibitory antibiotic treatment (Box 1). Bacteria can evolve resistance through several well-known mechanisms including alteration of sequence or expression of the target (Ince and Hooper, 2003), alteration of drug modification enzymes (Robicsek et al., 2006), drug efflux (Nikaido, 2009), or gene amplification (Andersson and Hughes, 2009; Sandegren and Andersson, 2009). However, unlike in other bacteria, genetic resistance in Mtb is not known to be associated with horizontal gene transfer and instead is commonly the result of single nucleotide polymorphisms and insertions or deletions (Namouchi et al., 2012; Dookie et al., 2018). This lack of promiscuous gene transfer somewhat simplifies the process of uncovering genetic determinants of drug resistance in Mtb, however complexities remain at every level.
Box 1. Definitions
Antibiotic resistance is a stably heritable trait that enables a bacterial population to both survive and replicate in the presence of an otherwise inhibitory antibiotic concentration. Mechanisms of antibiotic resistance are tightly associated with mutations in drug targets, activating enzymes, efflux systems, or membrane porins. Of these, mutations in targets and activators are by far most common in Mtb.
Antibiotic tolerance. We refer to bacteria that survive high or prolonged concentrations of antibiotics in the absence of stably heritable mechanisms of resistance as drug tolerant persisters. Progeny of tolerant cells exhibit a drug susceptibility profile similar to that of the parental strain. Of note, the nomenclature of drug tolerance can be challenging, as some researchers employ different terms (eg. – phenotypic resistance; non-heritable resistance, antibiotic indifference), and others define tolerance and persistence in more limited ways e.g. (Brauner et al., 2016; Balaban et al., 2019; Urbaniec et al., 2022). In general, these nomenclature distinctions serve to highlight particular subsets of tolerance. Mechanisms of drug tolerance are varied, including reduced metabolic activities, low ATP levels, toxin-antitoxin systems, and stringent response. Note that the definition we employ for drug tolerant persisters is agnostic as to form of tolerance or mechanism.
Resistance to first- and second-line anti-TB drugs is generally associated with known mutations at particular loci. These mutations have been reviewed previously (Almeida Da Silva and Palomino, 2011; Cohen et al., 2019) and are summarized in Table 1. Nearly all of these mutations confer resistance through the alteration of a drug target or activator. For example, the target of rifampin (RIF) is the β subunit of bacterial RNA polymerase (RNAP), which is encoded by the rpoB gene (Goldstein, 2014). Mutations in rpoB have been used to predict resistance to RIF with such a high degree of success (Goldstein, 2014) that an 81-bp region of the rpoB gene is now designated the RIF resistance-determining region (Ohno et al., 1996; Ramaswamy and Musser, 1998). However, not all cases of RIF resistance are so straightforward. Distinctions have been drawn between high- and low-level RIF resistance, phenotypes caused by mutations within and outside of the known RIF resistance-determining region of rpoB (Shea et al., 2021). Similarly, the presence of mutations in rpoA and rpoC (Andersson, 2006) as well as secondary site mutations in rpoB (Brandis and Hughes, 2013; Meftahi et al., 2016) are now known to play a role in the Mtb RIF resistance profile. Additionally, mutation, and therefore resistance, can be induced under drug pressure. For example, the mutation frequency to RIF was found to increase more than a thousand-fold during 14-days of monotherapy (Kayigire et al., 2017).
The classical model can shed light on most drug resistant Mtb strains circulating today. Indeed, identifying point mutations in specific loci is the basis of highly successful PCR-based tests for Mtb drug resistance (Stevens et al., 2017). Recently, the WHO has catalogued whole genome sequences and drug resistance profiles of 38,215 Mtb clinical strains (WHO, 2021a). This catalogue makes clear the value of the classical resistance model, while also revealing many mutations of unknown mechanism are linked to resistance. Indeed, for every anti-TB agent, there are resistant strains that continue to elude molecular genetic characterization. In addition, the pre-existing mutation model sheds little light on how drug resistance evolves in Mtb. However, in recent years an updated model has emerged that seeks to incorporate older, seemingly anomalous observations with newer, technology-driven insights to explain more completely the global Mtb drug resistance landscape. To take the 18th century writer and polymath Johann Wolfgang von Goethe badly out of context, “tolerance comes of age”.
The updated model: Tolerance on the pathway to resistance
In 1944, Joseph Bigger described a subpopulation of Staphylococci that survived exposure to penicillin without generating heritable resistance. When those cells were cultured in fresh media and then re-exposed to penicillin, they retained the parent strain’s level of susceptibility (Bigger, 1944). Since then, non-heritable survival in the face of antibiotics has been noted in a variety of bacteria exposed to different agents, including Mtb and the phenomenon has been given many names (McCune and Tompsett, 1956; Levin and Rozen, 2006). Here we refer to bacteria that do not harbor stably heritable resistance and yet survive significant antibiotic exposure as drug tolerant persisters (see Box 1). Note that our definition is agnostic as to the mechanism(s) by which tolerance occurs.
A possible link between tolerance and the evolution of resistance immediately suggests itself – other things being equal, the longer bacteria survive, the greater the opportunity to mutate to a stably antibiotic resistant state. This reasoning underpins mathematical models that link the two phenomena (Levin-Reisman et al., 2019), and explains why we discuss tolerance in a review about drug resistance. However, it is important to bear in mind that experimental validation for this link is so far limited to a very small number of examples (Levin-Reisman et al., 2017; Sebastian et al., 2017).
Tolerance and heterogeneity
It was recognized early on that most drugs are less effective on Mtb during infection than they are in vitro (McCune and Tompsett, 1956). One important reason is that most antibiotics work best on actively dividing cells, but a robust immune response is one of several mechanisms to slow the Mtb replication rate (Gill et al., 2009; Ford et al., 2011; Colangeli et al., 2014). More recently it has become evident that Mtb has evolved multiple strategies to generate subpopulations of phenotypically distinct bacteria, each with separate growth rates and levels of drug tolerance (Aldridge et al., 2012; Sarathy et al., 2018; Richards et al., 2019). In any given mycobacterial population, variations in replication, DNA repair, transcription, translation, metabolism, and efflux all promote bacterial heterogeneity and are also linked to drug tolerance (Szumowski et al., 2013; Chung et al., 2022). A similar phenomenon is evident within the human lung, where some lesions can support Mtb growth and expand during drug treatment even as other lesions shrink and the patient improves overall (Akira et al., 2000; Xie et al., 2021). Finally, heterogeneity exists among Mtb strains worldwide, driving differences in the accumulation and spread of drug resistance. Recently appreciation has been growing for how widely different mechanisms that promote and maintain bacterial heterogeneity are linked to drug tolerance, treatment failure, and ultimately the evolution of resistance. While this review makes no attempt to be comprehensive, some relevant examples of these mechanisms are provided below (Figure 2).
Figure 2 Examples of intracellular mechanisms to generate heterogeneity and tolerance in Mtb. Shown are various cellular processes (blue) at which heterogeneity and tolerance can be generated. Examples of known mechanisms are given at each level with examples in red. One expected outcome of these mechanisms is shown at the bottom, with a sensitive population (gray) becoming differentially tolerant (shades of yellow) after induction of any of these systems. Figure created with BioRender.com.
DNA replication, growth and division
Unlike other bacterial species, mycobacteria elongate from the cell poles, providing a straightforward opportunity for asymmetric growth and division leading to daughter cells with different sizes and growth rates (Aldridge et al., 2012; Santi et al., 2013; Chung et al., 2022). A key member of mycobacterial divisome complex, LamA has been identified as a mediator of asymmetric growth by inhibiting cell-wall synthesis at the new poles (Rego et al., 2017). In the absence of lamA, daughter cells are less heterogeneous in size and also display increased sensitivity to RIF and vancomycin compared to wild-type cells (Rego et al., 2017). During in vitro growth, Mtb cells shorter in length were found to be more susceptible to RIF, along with being more sensitive to oxidative and nitrite stress (Vijay et al., 2017). Clinically, Mtb with an increased cell length has been associated with host stresses such as iron deficiency and oxidative stress along with an increase in severe pulmonary disease (Vijay et al., 2017). Multidrug-resistant (MDR)-Mtb strains were also found to exhibit increased heterogeneity in cell length compared to drug-sensitive strains during intramacrophage growth and during RIF exposure, which when combined, was further associated with an increase in cell length (Vijay et al., 2017).
In a separate genome-wide association study of clinical Mtb strains, variants in the essential DNA replication initiation factor, dnaA, were found to be associated with drug resistance (Hicks et al., 2020). These dnaA variants enhance Mtb survival during isoniazid (INH) treatment by interacting with the Rv0010c-Rv0011c intergenic region and reduced expression of katG, the activator of INH. However, the connection between Rv0010c-Rv0011c and altered katG expression is not currently understood (Hicks et al., 2020).
Phase variation
Phase variation is an adaptive mechanism that mediates reversible switching of a gene by genotypic changes, which in turn can lead to reversible or transient drug resistance. Phase variation results from insertions or deletions in a homopolymeric tract (HT) located within the promoter or open reading frame of a gene. Reversible frameshift mutations in HTs are thought to be a result of slipped-strand mispairing errors during replication. Inactivating transient frameshift mutations in the HT of 7 cytosines in the glpK gene, which encodes a glycerol-3-kinase, were found to confer heritable drug resistance to INH, RIF, ethambutol (EMB), pyrazinamide (PZA) and moxifloxacin (MOX), but these changes were reversible with additional insertions or deletions in the same HT (Bellerose et al., 2019; Safi et al., 2019). These frameshift mutations were also identified in Mtb-infected mice and in clinical Mtb strains (Bellerose et al., 2019; Safi et al., 2019).
In another example, reversible frameshift mutations in the Mtb orn gene in combination with low-level EMB-resistant double embB-aftA mutant produced a small colony variant with a significantly higher MIC and further culture led to a rapid reversion of the orn frameshift (Safi et al., 2020).
Genetic regulators mediating tolerance and resistance
Heterogeneity in levels or activity of proteins and flux of metabolites is often facilitated by heterogeneity in gene expression. Recent work has provided examples of transcription factors and regulatory circuits that directly impact drug susceptibility. For instance, bedaquiline (BDQ) is a newer antitubercular drug that inhibits the F1F0-ATP synthase of Mtb. Network analysis of the Mtb transcriptional response to BDQ suggested that regulators Rv0324 and Rv0880 played important roles in the response to this drug (Peterson et al., 2016). Subsequent work showed that deletion of either Rv0324 or Rv0880 led to improved killing by BDQ but not other drugs. Exposure to capreomycin and MOX also induced Rv0324 suggesting potential antagonism with BDQ, while exposure to pretomanid decreased expression of the Rv0880 BDQ response regulon (Peterson et al., 2016). The addition of pretomanid to novel BDQ-containing regimens was found to improve clearance and reduce relapse in several murine models of TB (Xu et al., 2019), consistent with the synergistic effect of BDQ and pretomanid predicted by the transcriptional response and network analysis.
In another example, a network-based genetic screening approach, the transcriptional regulator-induced phenotype (TRIP) screen, was used to identify regulators that alter susceptibility to INH. Expression of mce3R, a TetR family transcription factor, was found to mediate heightened sensitivity to INH (Ma et al., 2021), which was linked to repression of ctpD, a gene encoding a metal binding protein (Raimunda et al., 2014; Patel et al., 2016; Ma et al., 2021) not previously known to play a role in INH susceptibility. Disruption of ctpD conferred hypersusceptibilty to INH, with increased intracellular accumulation of INH and INH-NAD adduct.
Drug induced tolerance
Drug exposure has also been shown to induce transcriptional changes conferring tolerance. When exposed to low RIF concentrations, transcription from one promoter is inhibited, allowing increased rpoB expression from a second promoter and after a delay, growth can resume despite ongoing drug exposure (Zhu et al., 2018). Subsequent studies have found that the absence of a functional lepA, a translation-associated elongation factor, increased RIF tolerance by the upregulation of rpoB expression (Wang et al., 2020). Mutations in lepA found in clinical Mtb strains were also found to have increased tolerance to RIF (Wang et al., 2020).
DNA repair
Environmental stresses such as hypoxia, nutrient deprivation, and host immune effectors can generate genome-wide mutations in subsets of cells, leading to heritable genetic resistance (Sebastian et al., 2017; Ragheb et al., 2019; Hicks et al., 2020; Jakkala et al., 2020; Swaminath et al., 2020). The activity and efficiency of different DNA repair pathways influences the rate at which these cells are a reservoir from which resistant isolates emerge.
A recently described link between DNA repair and drug resistance involved the DNA translocase protein Mfd. Mfd has long been thought to initiate nucleotide excision repair at sites of stalled RNAPs, though Mtb cells deficient in Mfd are not more sensitive to DNA-damaging agents. However, Mfd is found to associate with RNAP in the absence of DNA damage and has also been identified as an anti-backtracking factor for RNAP. Interactions of Mfd with the β-subunit of RNAP promote mutagenesis leading to drug resistance in several bacterial pathogens, including Mtb. Deletion of mfd led to a decrease in the number of spontaneous resistant mutants to RIF, EMB and ciprofloxacin (Ragheb et al., 2019). Recently, Rv1019, a transcriptional regulator of the TetR family, was found to negatively regulate mfd expression. Overexpression of Rv1019 leads to the downregulation of mfd and decreased Mtb survival under oxidative stress (Pushparajan et al., 2020). Since Rv1019 is differently expressed during hypoxia and reactivation (Schubert et al., 2015), it would be interesting to know if Rv1019 is the key regulator of Mfd-mediated changes in Mtb leading to drug resistance.
In an in vitro persistence model, Mtb was found to develop resistance to RIF or MOX at a much higher frequency than predicted when exposed to continuous lethal concentrations of RIF (Sebastian et al., 2017). Mtb in the RIF persistence-phase was found to carry elevated levels of hydroxyl radicals leading to genome-wide random mutagenesis, generating not only mutants in rpoB, but also in gyrA demonstrating that exposure to one antibiotic can generate genetic resistance to a different antibiotic (Sebastian et al., 2017). Similarly, exposure to lethal concentrations of MOX also generated high levels of hydroxyl radicals leading to resistance not only to MOX but also to EMB and INH (Swaminath et al., 2020).
Mistranslation during protein synthesis
Mistranslation happens when an error occurs during the protein synthesis yielding either incorporation of an incorrect amino acid or a truncated protein product. Generally deleterious, mistranslation can be adaptive in the presence of RIF. As noted above, most RIF resistance is due to mutation in a small region of rpoB, the RIF resistance determining region (Gagneux et al., 2006b). Increasing the mistranslation rate in M. smegmatis led to substitutions of glutamate for glutamine and aspartate for asparagine within the same region, which improved survival during RIF exposure (Javid et al., 2014). This effect is regulated by levels of the GatCAB enzyme complex, where reduced expression results in increased mistranslation and RIF tolerance (Su et al., 2016). Clinical strains with mutations in gatA, cause a partial loss of function of the complex along with increased mistranslation and increased RIF tolerance (Su et al., 2016; Li et al., 2021).
Metabolic regulation
Mtb can utilize a wide variety of carbon sources to support in vitro growth. However, Mtb resides in vivo within intracellular and extracellular niches where the nutrient composition is thought to be sparse and growth-limiting (Berney and Berney-Meyer, 2017; Sarathy et al., 2018). This slowed growth has traditionally been associated with drug tolerance, as most antibiotics target metabolically active Mtb (Schaefer, 1954). In Mtb, the regulation of several interconnected pathways that control carbon and lipid metabolism contribute to this metabolic slowdown. Rerouting pathways from energy-generation towards energy storage is associated with growth arrest and reduced drug susceptibility. Importantly, even stochastic differences in expression or activity of regulators and rate-limiting steps in these pathways results in bacterial heterogeneity and differential susceptibility to drugs.
For instance, under stress conditions (including drug pressure), Mtb can shift from the growth-promoting TCA cycle to carbon storage in fatty acids via triacylglycerol (TAG) synthesis via the upregulation of the triglyceride synthase tgs1 (Sirakova et al., 2006). In the absence of tgs1, drug tolerance induced during hypoxia is reversed and Mtb remains sensitive to INH, streptomycin (STR), fluoroquinolones (FQ) and EMB (Baek et al., 2011). The slowing of the TCA cycle alters the turnover of alpha-ketoglutarate, oxaloacetate and reducing agents such as NADH are diminished, resulting in reduced amino acid synthesis and protein translation. Additionally, enzymes such as isocitrate lysases, which are involved in the the glyoxylate bypass, are induced upon exposure to INH, RIF and STR conferring cross-tolerance (Nandakumar et al., 2014). Deficiency in the gluconeogenic enzyme phosphoenolpyruvate carboxykinase encoded by pckA, has been implicated in enhanced drug tolerance to INH and BDQ, with this tolerance associated with the overaccumulation of methylcitrate cycle (MCC) intermediates (Quinonez et al., 2022). Similarly, depletion of phosphoenolpyruvate during hypoxia confers tolerance to INH (Lim et al., 2021). In prpR mutants, defective regulation of MCC leads to an accumulation of propionyl-CoA conferring tolerance to INH, RIF and ofloxacin (OFX) (Hicks et al., 2018).
In another example, trehalose in Mtb serves as both a carbohydrate store as well as a component of the cell surface glycolipids trehalose monomycolate (TMM) and trehalose dimycolate (TDM). During hypoxia, Mtb has been shown to down-regulate TMM and TDM and channel trehalose into the biosynthesis of central carbon metabolism (CCM) intermediates. In a biofilm model, drug-tolerant persisters were shown to shift trehalose metabolism towards CCM intermediates (Lee et al., 2019). treS deletion mutants were unable to shift trehalose metabolism to CMM and showed a rapid depletion of ATP and were also significantly more susceptible to BDQ (Lee et al., 2019).
Efflux
The Mtb genome encodes a significant number of efflux pumps with a known role in intrinsic and acquired drug resistance, and many of these pumps are also induced during infection (Schnappinger et al., 2003; Rengarajan et al., 2005; Gupta et al., 2010; Szumowski et al., 2013; Pule et al., 2016). Since efflux pumps can have broad substrate specificities, their induction under varying environmental conditions or drug exposure may confer tolerance or resistance to multiple drugs, as is seen with RIF exposure resulting in tolerance to OXF (Louw et al., 2011). Exposure to INH has been shown to induce efflux pumps mediating its tolerance, which then promote the emergence of genetically resistant INH strains (Machado et al., 2012). Additionally, upon infection of macrophages, Mtb has been shown to induce tolerance to numerous antitubercular drugs including INH, RIF, MOX and BDQ that is not tied to reduced growth rate and is sensitive to efflux pump inhibitors such as verapamil (Adams et al., 2011; Adams et al., 2014). Macrophage-induced tolerance to RIF was shown to be mediated by Rv1258c/Tap, an efflux pump also important for intracellular growth (Adams et al., 2011). Further study revealed that strains from all tested global lineages developed macrophage-induced tolerance to RIF except lineage 2 Beijing isolates, which harbor a natural loss-of-function mutation in Rv1258c (Villellas et al., 2013; Adams et al., 2019). In addition, mutations in Rv1258c that were identified in clinical isolates have been linked with resistance to INH, PZA, and STR (Liu et al., 2019).
Within-host Mtb differences affect drug response
One limitation with the classical model of resistance is that it does not consider within-host variation during TB infection; however, recent studies are beginning to bring evidence of heterogeneity in the Mtb response during drug treatment to light (Borrell and Gagneux, 2009; McGrath et al., 2014). Clinical Mtb strains were collected from patients with delayed culture conversion and WGS was performed. Exposing these strains to RIF in vitro revealed drug tolerant variants undetected by bulk WGS-analysis (Genestet et al., 2021). One variant identified by RIF treatment was also enriched during macrophage infection and was found to have a mutation in the gene mas, altering its cell surface lipids (Genestet et al., 2021). Characterizing these tolerant sub-populations may help identify patients at risk for treatment failure and the evolution of resistance.
In another study, Mtb isolates from treatment-naïve patients were subjected to WGS to assess within-host diversity. The accumulation of identified mutations varied substantially between isolates from the same individual and were elevated in HIV-negative patients, suggesting that the host immune environment may influence mutation rates (Liu et al., 2020). These results argue that the risk of developing new drug resistance mutations in vivo may vary with the host immune environment (Liu et al., 2020). This idea is consistent with evidence from TB patients that host gene signatures exhibiting heightened inflammatory and immune gene expression correlate with longer times to cure and a reduced cure rate (DiNardo et al., 2022).
Lineage-specific Mtb differences and drug response
It has become clear in the last several years that global variations in Mtb strains have strong impacts on drug response and the evolution of drug resistance. Worldwide Mtb has been separated into seven lineages and many sub-lineages with distinct characteristics that co-evolved with the human populations in which they are present (Brites and Gagneux, 2015). Global lineages differ in their ability to respond to drugs and develop resistance. Members of the modern lineages (2, 3 and 4) are associated with greater disease burden and drug resistance than the ancient lineages, possibly due to an increased spontaneous mutation rate (Borrell and Gagneux, 2009; McGrath et al., 2014). Further, differences in vitro and in vivo, the genetic background of the strain and the nature of the specific resistance mutation both influence outcomes. For example, there is an association between strain lineage and the type of resistance mutation identified, suggesting that certain Mtb lineages may have characteristics that encourage different routes to resistance (Gagneux et al., 2006a). One study using TnSeq showed that clinically distinct strains have different requirements for in vitro growth, including katG and glcB (Carey et al., 2018). The differences in TnSeq phenotypes of these strains were found to predict their drug resistance rates (Carey et al., 2018).
Compensatory mutations
Antibiotics, by their nature, target important functions of the bacterial cell. Thus, any mutation that renders a strain resistant has good potential to also reduce the strain’s fitness. This observation once led to the hope that simply reducing the use of antibiotics would lead to fitter, susceptible strains outcompeting resistant ones. However, reduced fitness can also be addressed by compensatory evolution and genetic co-selection (Andersson and Levin, 1999; Andersson, 2006; Andersson and Hughes, 2010).
The fitness cost of a resistance mutation can be measured in vitro with isogenic strains serially passaged or grown continuously in chemostats. Such results do not always translate into a host setting, so it is important to also consider how virulence and pathogen transmission are affected. For example, in the 1950s Middlebrook and colleagues found that many INH-resistant Mtb strains were less virulent in Guinea pigs (Middlebrook, 1954; Widelock et al., 1955; Wolinsky et al., 1956). Later, it was revealed that INH-resistant Mtb lacking KatG catalase-peroxidase activity could compensate by overexpressing an alkyl hydroperoxidase (Sherman et al., 1996). Similarly, with regard to RIF resistance, it has been shown that prolonged treatment can result in multidrug resistant strains that have no measurable fitness defect (Gagneux et al., 2006b). These examples illustrate the complex relationship of drug resistance and fitness, where initial costs can be corrected by compensatory mutations that retain the resistance phenotype. It is an important consideration, as such low and no cost mutations have been seen in clinical isolates (Sander et al., 2002). Specific compensatory mutations are shown in Table 1 and have been reviewed elsewhere (Alame Emane et al., 2021).
One non-canonical form of compensatory mutation that was recently described in Mtb involves restoring fitness of a capreomycin (CAP)-resistant mutant. CAP binds to 16S rRNA and inhibits translation. CAP resistance is conferred by 16S rRNA mutation that also reduces translation efficiency. However, translation can be largely restored by increased expression of an enzyme that methylates a nearby site on the 16S rRNA, significantly reducing the fitness cost of CAP resistance (Freihofer et al., 2016). This is a striking example of a compensatory mutation that relies on changes in expression but acts through post-transcriptional modification. Evidence of these ‘multi-level’ mechanisms of fitness alterations are rare, but it seems likely that more will be discovered and shown to be relevant in other contexts.
Epistasis
Epistasis refers to a phenomenon where the phenotypic effect of a particular gene allele depends on its genomic background (Hughes and Andersson, 2017). In the context of antibiotic resistance, epistatic interactions between resistance-conferring mutations have a major influence on the fitness of the multidrug-resistant (MDR) isolates and hence their evolution. Epistasis can have either positive or negative outcomes depending on the net effect on bacterial fitness in the absence of antibiotic pressure. Positive epistasis occurs when the fitness cost associated with multiple resistance-conferring mutations is lower than the anticipated additive cost of these mutations if calculated independently. Positive epistasis is a common phenomenon in mycobacteria and in other bacteria as illustrated in numerous studies. For example, a study by Borrel et al. showed that MDR isolates harboring double mutations in rpoB H526Y and gyrA D94G, conferring resistance to RIF and ofloxacin respectively, were associated with enhanced fitness as opposed to their respective single mutants (Borrell et al., 2013). Similarly, Sun et al. reported a positive epistatic interaction in MDR isolates with double mutations in rpsL K43M and gyrA D94Y, which confer resistance to STR and fluoroquinolones respectively (Sun et al., 2018). Another example of positive epistasis was also reported by Li et al. where MDR isolates with dual mutations in rpoB C531T and katG 315C were associated with enhanced fitness (Li et al., 2017). Importantly, those MDR isolates where positive epistasis conferred fitness benefits were associated with better transmissibility and thus were frequently encountered in clinical settings, which supports the idea that positive epistasis plays an important role in the evolution of MDR isolates (Trindade et al., 2009; Borrell et al., 2013).
In addition, positive lineage-specific epistatic interactions were found to be associated with particular Mtb clades. One study identified two epistatic interactions that were exclusively observed in lineage 4 (Coll et al., 2018). Compensatory mutations in pnaB2 and thyX–hsdS.1 promoter were found to be associated with resistances to PZA and para-aminosalicylic acid (PAS), due to mutations in pncA and thyA, respectively.
On the other hand, negative epistasis occurs when the interaction between two or more resistance-conferring mutations aggravates the fitness cost associated with these mutations. For example, FQ resistant isolates with double mutations in gyrA and gyrB were associated with diminished fitness (Luo et al., 2017). Similarly, in the Borrel et al. study, negative epistasis was observed in MDR isolates with double mutations in rpoB H526P and gyrA G88C (Borrell et al., 2013). In contrast to positive epistasis, MDR isolates where epistatic interactions resulted in diminished fitness were associated with low transmission rates and were rarely encountered in clinics.
Epistasis and the evolution of resistance
Several studies have revealed a strong correlation between Mtb lineages and particular drug resistance-conferring mutations, highlighting the major influence of the genetic background on the evolution of drug resistance. For example, one study reported lineage-specific differences in the level of INH resistance due to mutations in katG and inhA (Fenner et al., 2012). The katG mutations were more prevalent in lineage 2 isolates and conferred a high level of INH resistance, whereas inhA mutations were more prevalent in lineage 1 and were associated with low levels of INH resistance. Another study found that fluoroquinolone resistance due to mutated gyrA occurred more frequently in lineages 2 and 4 (Castro et al., 2020). Similarly, katG and rpoB mutations occur more frequently in modern Beijing sublineages compared to the ancient strains (Li et al., 2017). These examples demonstrate the major influence the genetic background could have on the evolution of drug response and also may explain why some Mtb lineages, particularly Beijing isolates, are often associated with multidrug resistance (Fenner et al., 2012; Nieto Ramirez et al., 2020; Fursov et al., 2021).
Epistasis and the level of drug resistance
The classical resistance model relies on using specific genetic determinants to define drug resistance. However, a key limitation of this model is that it tends to ignore the effect of epistasis on the level of drug resistance. Several studies have recently shown that bacterial cells can epistatically exhibit enhanced drug susceptibility despite the presence of a resistance-conferring mutation. For example, a study showed that a loss of function mutation in the eis coding region was able to restore amikacin susceptibility in resistant isolates harboring eis C-14T mutation (Vargas et al., 2021). Moreover, the same study questioned the validity of mmpR mutations as a determinant of bedaquiline and clofazimine resistances if loss of function mutations in mmpS5 and mmpL5 were present concomitantly (Vargas et al., 2021).
Summary and conclusions
So where does drug resistance in Mtb come from, and where is it going? Historically, the classical model (Figure 1, left side), in which pre-existing mutations are selected by drug pressure, has proven a very useful framework for our evolving understanding of resistance. However, the updated model (Figure 1, right side), with non-obligatory steps through tolerance on the path to resistance, does a better job of describing the rates, types, and patterns of drug resistance within communities and around the world. It is clear that pre-existing mutations conferring resistance do exist in any population of sufficient size, and that resistance does not require a tolerant pre-step. In practice however, with so many different routes to a tolerant state, it is entirely possible that the majority of resistant isolates worldwide today emerged from drug-tolerant precursors.
As described above, tolerance can be stochastic or genetically programmed, and it is frequently linked with the formation and maintenance of heterogeneous sub-populations. Heterogeneity can be recognized at all levels of TB disease, including bacterial subpopulations within individual lesions and across lesions in a single patient, within a single patient over time, within communities, and in different lineages across the globe. In all these cases, we should expect that heterogeneity contributes both to treatment failure and evolution of resistance. Further, since the drug tolerance spawned by heterogeneity is adaptive, we might anticipate that future work will demonstrate that the production of heterogeneous sub-populations is itself under genetic selection and control. In fact, multiple recent reports already point in that direction (Bellerose et al., 2019; Safi et al., 2019; Safi et al., 2020; Ma et al., 2021; Carey et al., 2022; Martini et al., 2022).
The updated model has implications for how we track drug susceptibility and resistance. Increasingly around the world, slow and labor-intensive microbiological drug susceptibility testing is being replaced by faster DNA-based methods, either PCR or next-generation sequencing (NGS) (Dookie et al., 2022; Rowlinson and Musser, 2022). For example, a recent study reported the whole genome sequences of more than 12,000 Mtb clinical isolates, along with sensitivity data for 13 different drugs (The CRyPTIC Consortium, 2022). DNA-based methods offer important advantages in speed, throughput, and safety, as well as altogether novel insights into drug resistance mechanisms (Hicks et al., 2019; The CRyPTIC Consortium, 2022). Catalogues of sequencing results should be widely available and analyzed regularly for potential associations and emerging mutations of interest. However, it is important that these methods are implemented with stringent controls for DNA extraction, sequencing and data handling. NGS sequencing and analysis are not always straightforward (Villellas et al., 2017; Kaniga et al., 2022), and global standards for the application of NGS data to drug susceptibility testing and data reporting are still emerging. Further, there are cases of discordance between whole genome sequencing and drug sensitivity data (Dookie et al., 2022), though these are rare and the extent to which they are due to experimental errors is not yet clear. Finally, NGS generally provides only a snapshot of a bulk sample, without conveying the subtleties of the heterogeneous subpopulations described above. Technologies are in development (Box 2) for the identification and characterization of subpopulations, but these are not yet commonplace, and are certainly not yet employed for drug susceptibility testing. Altogether, DNA-based approaches are revolutionizing how we monitor drug susceptibility and resistance and show much promise for further advances but making good on that promise will require both the development of new tools and the rigorous application of those tools in the lab and the clinic.
Box 2. Technology and our understanding of antibiotic resistance
As is common in biology, technological advances have been critical in updating our concepts of antibiotic action and the evolution of resistance. Important advances have occurred in:
Visualization – Advanced visualization tools such as multiparameter confocal microscopy (Gern et al., 2021; Plumlee et al., 2021) and mass cytometry (Xu et al., 2021) are helping to uncover the complexity of the host response to TB infection. Positron emission tomography (PET) imaging has brought to light the heterogeneity of TB lesions in live animals and humans (Lenaerts et al., 2015). Microfluidics (Molloy et al., 2021) and time-lapse microscopy (Herricks et al., 2020) are revealing the complexity of Mtb populations in vitro, and reporter gene technology (Huang et al., 2019) is providing similar insights in vivo.
Next-generation sequencing – High throughput sequencing has revolutionized to study of drug resistance. With thousands of Mtb genomes sequenced, the diversity of the Mtb pan genome is now evident. Many novel mutations have been linked with resistance to particular drugs, either alone or in association with known resistance loci (Zhang et al., 2013; Zeng et al., 2018; The CRyPTIC Consortium, 2022). Each new mutation must then be studied to see if it truly confers resistance or compensates for fitness defects imposed by mutations at other sites. Single-cell RNA-seq (Pisu et al., 2021) and dual-seq that simultaneously captures transcriptomes of Mtb and host (Pisu et al., 2020) have become important tools to study rare cell types and sub-populations in vivo. In addition, next-gen sequencing is central to the updated genetic screens described below.
Molecular genetics – Updated approaches in molecular genetics are also shedding new light on antibiotic action and resistance. Tn-seq is not really new technology, but it is being used to gain new insights into resistance mechanisms (Carey et al., 2018). Similarly, CRISPRi screens are identifying new loci associated with resistance to different agents (Li et al., 2022). Also, network-based TRIP screens have been employed to identify novel regulons and effector genes linked to drug sensitivity, tolerance and resistance (Ma et al., 2021).
The updated model also suggests new approaches to combat the emergence of drug resistance. If resistance frequently emerges from tolerant cells, then strategies to kill drug tolerant persisters or restrict their formation should slow the emergence of resistance. While not the topic of this review, eliminating persisters may also shorten the course of current therapy (Chung et al., 2022), so efforts to develop anti-persister therapy should receive high priority. Assays that identify small molecules targeting specific persister subpopulations have been reported (Sukheja et al., 2017; Gold et al., 2021). Hits from these screens could be combined with recent work to identify synergistic drug combinations (Cokol et al., 2017; Katzir et al., 2019; Ma et al., 2019) that can target multiple subpopulations at once. It may also be possible to directly target the machinery that promotes tolerance and resistance. For example, small molecules that inhibit the action of mfd (Ragheb et al., 2019) or DNA repair enzymes (Reiche et al., 2017) should reduce the rate at which resistance to other drugs emerge. Also, since small expression changes can have substantial effects on drug tolerance (Ma et al., 2021) and treatment outcome (Colangeli et al., 2018), it should be possible to identify small molecules that specifically alter Mtb gene expression away from tolerance-promoting states. When combined with NGS to characterize individual strains and efforts to uncover host-directed therapies, it is possible to imagine these approaches promoting an era of personalized TB therapy to achieve both improved outcomes and diminished resistance.
In conclusion, the history of drug development argues that resistance will emerge following the introduction of virtually any new agent. However, as the field has evolved and a new model of resistance has emerged, new strategies to protect and preserve agents can be envisioned. With careful monitoring and thoughtful development, we need not tolerate the loss of new agents to TB drug resistance any longer.
Author contributions
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
Funding
This review was supported by grants to DS: U19 AI 135976, U19 AI 1625598, and R01 AI 146194. The funders had no role in the decision to publish or prepare this manuscript.
Acknowledgments
The authors wish to thank David Alland and Shuyi Ma for helpful discussions and suggestions, and Eleanor Lamont for critical reading of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Adams, K. N., Szumowski, J. D., Ramakrishnan, L. (2014). Verapamil, and its metabolite norverapamil, inhibit macrophage-induced, bacterial efflux pump-mediated tolerance to multiple anti-tubercular drugs. J. Infect. Dis. 210 (3), 456–466. doi: 10.1093/infdis/jiu095
Adams, K. N., Takaki, K., Connolly, L. E., Wiedenhoft, H., Winglee, K., Humbert, O., et al. (2011). Drug tolerance in replicating mycobacteria mediated by a macrophage-induced efflux mechanism. Cell 145 (1), 39–53. doi: 10.1016/j.cell.2011.02.022
Adams, K. N., Verma, A. K., Gopalaswamy, R., Adikesavalu, H., Singhal, D. K., Tripathy, S., et al. (2019). Diverse clinical isolates of Mycobacterium tuberculosis develop macrophage-induced rifampin tolerance. J. Infect. Dis. 219 (10), 1554–1558. doi: 10.1093/infdis/jiy710
Akira, M., Sakatani, M., Ishikawa, H. (2000). Transient radiographic progression during initial treatment of pulmonary tuberculosis: CT findings. J. Comput. Assist. Tomogr 24 (3), 426–431. doi: 10.1097/00004728-200005000-00014
Alame Emane, A. K., Guo, X., Takiff, H. E., Liu, S. (2021). Drug resistance, fitness and compensatory mutations in mycobacterium tuberculosis. Tuberculosis (Edinb) 129, 102091. doi: 10.1016/j.tube.2021.102091
Alcaide, F., Pfyffer, G. E., Telenti, A. (1997). Role of embB in natural and acquired resistance to ethambutol in mycobacteria. Antimicrob. Agents Chemother. 41 (10), 2270–2273. doi: 10.1128/AAC.41.10.2270
Aldridge, B. B., Fernandez-Suarez, M., Heller, D., Ambravaneswaran, V., Irimia, D., Toner, M., et al. (2012). Asymmetry and aging of mycobacterial cells lead to variable growth and antibiotic susceptibility. Science 335 (6064), 100–104. doi: 10.1126/science.1216166
Almeida Da Silva, P. E., Palomino, J. C. (2011). Molecular basis and mechanisms of drug resistance in Mycobacterium tuberculosis: classical and new drugs. J. Antimicrob. Chemother. 66 (7), 1417–1430. doi: 10.1093/jac/dkr173
Almeida, D., Ioerger, T., Tyagi, S., Li, S. Y., Mdluli, K., Andries, K., et al. (2016). Mutations in pepQ confer low-level resistance to bedaquiline and clofazimine in Mycobacterium tuberculosis. Antimicrob. Agents Chemother. 60 (8), 4590–4599. doi: 10.1128/AAC.00753-16
Andersson, D. I. (2006). The biological cost of mutational antibiotic resistance: any practical conclusions? Curr. Opin. Microbiol. 9 (5), 461–465. doi: 10.1016/j.mib.2006.07.002
Andersson, D. I., Hughes, D. (2009). Gene amplification and adaptive evolution in bacteria. Annu. Rev. Genet. 43, 167–195. doi: 10.1146/annurev-genet-102108-134805
Andersson, D. I., Hughes, D. (2010). Antibiotic resistance and its cost: is it possible to reverse resistance? Nat. Rev. Microbiol. 8 (4), 260–271. doi: 10.1038/nrmicro2319
Andersson, D. I., Levin, B. R. (1999). The biological cost of antibiotic resistance. Curr. Opin. Microbiol. 2 (5), 489–493. doi: 10.1016/S1369-5274(99)00005-3
Andries, K., Verhasselt, P., Guillemont, J., Gohlmann, H. W., Neefs, J. M., Winkler, H., et al. (2005). A diarylquinoline drug active on the ATP synthase of Mycobacterium tuberculosis. Science 307 (5707), 223–227. doi: 10.1126/science.1106753
Andries, K., Villellas, C., Coeck, N., Thys, K., Gevers, T., Vranckx, L., et al. (2014). Acquired resistance of Mycobacterium tuberculosis to bedaquiline. PloS One 9 (7), e102135. doi: 10.1371/journal.pone.0102135
Baek, S. H., Li, A. H., Sassetti, C. M. (2011). Metabolic regulation of mycobacterial growth and antibiotic sensitivity. PloS Biol. 9 (5), e1001065. doi: 10.1371/journal.pbio.1001065
Balaban, N. Q., Helaine, S., Lewis, K., Ackermann, M., Aldridge, B., Andersson, D. I., et al. (2019). Definitions and guidelines for research on antibiotic persistence. Nat. Rev. Microbiol. 17 (7), 441–448. doi: 10.1038/s41579-019-0196-3
Barberis, I., Bragazzi, N. L., Galluzzo, L., Martini, M. (2017). The history of tuberculosis: from the first historical records to the isolation of koch's bacillus. J. Prev. Med. Hyg 58 (1), E9–E12. doi: 10.15167/2421-4248/jpmh2017.58.1.728
Beckert, P., Hillemann, D., Kohl, T. A., Kalinowski, J., Richter, E., Niemann, S., et al. (2012). rplC T460C identified as a dominant mutation in linezolid-resistant Mycobacterium tuberculosis strains. Antimicrob. Agents Chemother. 56 (5), 2743–2745. doi: 10.1128/AAC.06227-11
Bellerose, M. M., Baek, S. H., Huang, C. C., Moss, C. E., Koh, E. I., Proulx, M. K., et al. (2019). Common variants in the glycerol kinase gene reduce tuberculosis drug efficacy. mBio 10 (4), 1–15. doi: 10.1128/mBio.00663-19
Berney, M., Berney-Meyer, L. (2017). Mycobacterium tuberculosis in the face of host-imposed nutrient limitation. Microbiol. Spectr. 5 (3), 1–17. doi: 10.1128/microbiolspec.TBTB2-0030-2016
Bigger, J. W. (1944). Treatment of staphylococcal infections with penicillin by intermittent sterilisation. Lancet 244 (6320), 497–500. doi: 10.1016/s0140-6736(00)74210-3
Borrell, S., Gagneux, S. (2009). Infectiousness, reproductive fitness and evolution of drug-resistant Mycobacterium tuberculosis. Int. J. tuberculosis Lung Dis. 13 (12), 1456–1466.
Borrell, S., Teo, Y., Giardina, F., Streicher, E. M., Klopper, M., Feldmann, J., et al. (2013). Epistasis between antibiotic resistance mutations drives the evolution of extensively drug-resistant tuberculosis. Evolution Medicine Public Health 2013 (1), 65–74. doi: 10.1093/emph/eot003
Brandis, G., Hughes, D. (2013). Genetic characterization of compensatory evolution in strains carrying rpoB Ser531Leu, the rifampicin resistance mutation most frequently found in clinical isolates. J. Antimicrob. Chemother. 68 (11), 2493–2497. doi: 10.1093/jac/dkt224
Brauner, A., Fridman, O., Gefen, O., Balaban, N. Q. (2016). Distinguishing between resistance, tolerance and persistence to antibiotic treatment. Nat. Rev. Microbiol. 14 (5), 320–330. doi: 10.1038/nrmicro.2016.34
Brites, D., Gagneux, S. (2015). Co-Evolution of Mycobacterium tuberculosis and Homo sapiens. Immunol. Rev. 264 (1), 6–24. doi: 10.1111/imr.12264
Carey, A. F., Rock, J. M., Krieger, I. V., Chase, M. R., Fernandez-Suarez, M., Gagneux, S., et al. (2018). TnSeq of Mycobacterium tuberculosis clinical isolates reveals strain-specific antibiotic liabilities. PloS Pathog. 14 (3), e1006939. doi: 10.1371/journal.ppat.1006939
Carey, A. F., Wang, X., Cicchetti, N., Spaulding, C. N., Liu, Q., Hopkins, F., et al. (2022). Multiplexed strain phenotyping defines consequences of genetic diversity in Mycobacterium tuberculosis for infection and vaccination outcomes. mSystems 7 (3), e0011022. doi: 10.1128/msystems.00110-22
Castro, R. A. D., Ross, A., Kamwela, L., Reinhard, M., Loiseau, C., Feldmann, J., et al. (2020). The genetic background modulates the evolution of fluoroquinolone-resistance in Mycobacterium tuberculosis. Mol. Biol. Evol. 37 (1), 195–207. doi: 10.1093/molbev/msz214
Chung, E. S., Johnson, W. C., Aldridge, B. B. (2022). Types and functions of heterogeneity in mycobacteria. Nat. Rev. Microbiol 20, 529–541. doi: 10.1038/s41579-022-00721-0
Cohen, K. A., Manson, A. L., Desjardins, C. A., Abeel, T., Earl, A. M. (2019). Deciphering drug resistance in Mycobacterium tuberculosis using whole-genome sequencing: progress, promise, and challenges. Genome Med. 11 (1), 45. doi: 10.1186/s13073-019-0660-8
Cokol, M., Kuru, N., Bicak, E., Larkins-Ford, J., Aldridge, B. B. (2017). Efficient measurement and factorization of high-order drug interactions in Mycobacterium tuberculosis. Sci. Adv. 3 (10), e1701881. doi: 10.1126/sciadv.1701881
Colangeli, R., Arcus, V. L., Cursons, R. T., Ruthe, A., Karalus, N., Coley, K., et al. (2014). Whole genome sequencing of Mycobacterium tuberculosis reveals slow growth and low mutation rates during latent infections in humans. PloS One 9 (3), e91024. doi: 10.1371/journal.pone.0091024
Colangeli, R., Jedrey, H., Kim, S., Connell, R., Ma, S., Venkata, U. D. C., et al. (2018). Bacterial factors that predict relapse after tuberculosis therapy. New Engl. J. Med. 379 (9), 823–833. doi: 10.1056/nejmoa1715849
Coll, F., Phelan, J., Hill-Cawthorne, G. A., Nair, M. B., Mallard, K., Ali, S., et al. (2018). Genome-wide analysis of multi- and extensively drug-resistant Mycobacterium tuberculosis. Nat. Genet. 50 (2), 307–316. doi: 10.1038/s41588-017-0029-0
Comas, I., Borrell, S., Roetzer, A., Rose, G., Malla, B., Kato-Maeda, M., et al. (2011). Whole-genome sequencing of rifampicin-resistant Mycobacterium tuberculosis strains identifies compensatory mutations in RNA polymerase genes. Nat. Genet. 44 (1), 106–110. doi: 10.1038/ng.1038
Crofton, J., Mitchison, D. A. (1948). Streptomycin resistance in pulmonary tuberculosis. Br. Med. J. 2 (4588), 1009–1015. doi: 10.1136/bmj.2.4588.1009
Daniel, T. M. (2006). The history of tuberculosis. Respir. Med. 100 (11), 1862–1870. doi: 10.1016/j.rmed.2006.08.006
de Jonge, M. R., Koymans, L. H., Guillemont, J. E., Koul, A., Andries, K. (2007). A computational model of the inhibition of Mycobacterium tuberculosis ATPase by a new drug candidate R207910. Proteins 67 (4), 971–980. doi: 10.1002/prot.21376
DiNardo, A. R., Gandhi, T., Heyckendorf, J., Grimm, S. L., Rajapakshe, K., Nishiguchi, T., et al. (2022). Gene expression signatures identify biologically and clinically distinct tuberculosis endotypes. Eur. Respir. J 60 (3), 1–15. doi: 10.1183/13993003.02263-2021
Dookie, N., Khan, A., Padayatchi, N., Naidoo, K. (2022). Application of next generation sequencing for diagnosis and clinical management of drug-resistant tuberculosis: Updates on recent developments in the field. Front. Microbiol. 13. doi: 10.3389/fmicb.2022.775030
Dookie, N., Naidoo, K., Padayatchi, N. (2018). Whole-genome sequencing to guide the selection of treatment for drug-resistant tuberculosis. Antimicrob. Agents Chemother. 62 (8), 1–2. doi: 10.1128/AAC.00574-18
Fenner, L., Egger, M., Bodmer, T., Altpeter, E., Zwahlen, M., Jaton, K., et al. (2012). Effect of mutation and genetic background on drug resistance in Mycobacterium tuberculosis. Antimicrob. Agents Chemother. 56 (6), 3047–3053. doi: 10.1128/AAC.06460-11
Ford, C. B., Lin, P. L., Chase, M. R., Shah, R. R., Iartchouk, O., Galagan, J., et al. (2011). Use of whole genome sequencing to estimate the mutation rate of Mycobacterium tuberculosis during latent infection. Nat. Genet. 43 (5), 482–486. doi: 10.1038/ng.811
Freihofer, P., Akbergenov, R., Teo, Y., Juskeviciene, R., Andersson, D. I., Böttger, E. C. (2016). Nonmutational compensation of the fitness cost of antibiotic resistance in mycobacteria by overexpression of tlyA rRNA methylase. RNA 22 (12), 1836–1843. doi: 10.1261/rna.057257.116
Fursov, M. V., Shitikov, E. A., Lagutkin, D. A., Fursova, A. D., Ganina, E. A., Kombarova, T. I., et al. (2021). MDR and pre-XDR clinical mycobacterium tuberculosis Beijing strains: Assessment of virulence and host cytokine response in mice infectious model. Microorganisms 9 (8), 1–17. doi: 10.3390/microorganisms9081792
Gagneux, S., Burgos, M. V., DeRiemer, K., Enciso, A., Muñoz, S., Hopewell, P. C., et al. (2006a). Impact of bacterial genetics on the transmission of isoniazid-resistant Mycobacterium tuberculosis. PloS Pathog. 2 (6), e61. doi: 10.1371/journal.ppat.0020061
Gagneux, S., Long, C. D., Small, P. M., Van, T., Schoolnik, G. K., Bohannan, B. J. M. (2006b). The competitive cost of antibiotic resistance in Mycobacterium tuberculosis. Science 312 (5782), 1944–1946. doi: 10.1126/science.1124410
Genestet, C., Hodille, E., Barbry, A., Berland, J. L., Hoffmann, J., Westeel, E., et al. (2021). Rifampicin exposure reveals within-host Mycobacterium tuberculosis diversity in patients with delayed culture conversion. PloS Pathog. 17 (6), e1009643. doi: 10.1371/journal.ppat.1009643
Gern, B. H., Adams, K. N., Plumlee, C. R., Stoltzfus, C. R., Shehata, L., Moguche, A. O., et al. (2021). TGFbeta restricts expansion, survival, and function of T cells within the tuberculous granuloma. Cell Host Microbe 29 (4), 594–606.e596. doi: 10.1016/j.chom.2021.02.005
Gillespie, S. H. (2002). Evolution of drug resistance in Mycobacterium tuberculosis: clinical and molecular perspective. Antimicrob. Agents Chemother. 46 (2), 267–274. doi: 10.1128/AAC.46.2.267-274.2002
Gillespie, S. H. (2007). Tuberculosis: evolution in millennia and minutes. Biochem. Soc. Trans. 35 (Pt 5), 1317–1320. doi: 10.1042/BST0351317
Gill, W. P., Harik, N. S., Whiddon, M. R., Liao, R. P., Mittler, J. E., Sherman, D. R. (2009). A replication clock for Mycobacterium tuberculosis. Nat. Med. 15 (2), 211–214. doi: 10.1038/nm.1915
Goldstein, B. P. (2014). Resistance to rifampicin: a review. J. Antibiot (Tokyo) 67 (9), 625–630. doi: 10.1038/ja.2014.107
Gold, B., Warrier, T., Nathan, C. (2021). A multistress model for high throughput screening against nonreplicating Mycobacterium tuberculosis. Methods Mol. Biol. 2314, 611–635. doi: 10.1007/978-1-0716-1460-0_27
Gomez-Gonzalez, P. J., Perdigao, J., Gomes, P., Puyen, Z. M., Santos-Lazaro, D., Napier, G., et al. (2021). Genetic diversity of candidate loci linked to Mycobacterium tuberculosis resistance to bedaquiline, delamanid and pretomanid. Sci. Rep. 11 (1), 19431. doi: 10.1038/s41598-021-98862-4
Gupta, A. K., Katoch, V. M., Chauhan, D. S., Sharma, R., Singh, M., Venkatesan, K., et al. (2010). Microarray analysis of efflux pump genes in multidrug-resistant Mycobacterium tuberculosis during stress induced by common anti-tuberculous drugs. Microb. Drug Resist. 16 (1), 21–28. doi: 10.1089/mdr.2009.0054
Haver, H. L., Chua, A., Ghode, P., Lakshminarayana, S. B., Singhal, A., Mathema, B., et al. (2015). Mutations in genes for the F420 biosynthetic pathway and a nitroreductase enzyme are the primary resistance determinants in spontaneous in vitro-selected PA-824-resistant mutants of Mycobacterium tuberculosis. Antimicrob. Agents Chemother. 59 (9), 5316–5323. doi: 10.1128/AAC.00308-15
Herricks, T., Donczew, M., Mast, F. D., Rustad, T., Morrison, R., Sterling, T. R., et al. (2020). ODELAM, rapid sequence-independent detection of drug resistance in isolates of Mycobacterium tuberculosis. Elife 9, 1–23. doi: 10.7554/eLife.56613
Heym, B., Alzari, P. M., Honore, N., Cole, S. T. (1995). Missense mutations in the catalase-peroxidase gene, katG, are associated with isoniazid resistance in Mycobacterium tuberculosis. Mol. Microbiol. 15 (2), 235–245. doi: 10.1111/j.1365-2958.1995.tb02238.x
Hicks, N. D., Carey, A. F., Yang, J., Zhao, Y., Fortune, S. M. (2019). Bacterial genome-wide association identifies novel factors that contribute to ethionamide and prothionamide susceptibility in mycobacterium tuberculosis. mBio 10 (2), 1–14. doi: 10.1128/mBio.00616-19
Hicks, N. D., Giffen, S. R., Culviner, P. H., Chao, M. C., Dulberger, C. L., Liu, Q., et al. (2020). Mutations in dnaA and a cryptic interaction site increase drug resistance in Mycobacterium tuberculosis. PloS Pathog. 16 (11), e1009063. doi: 10.1371/journal.ppat.1009063
Hicks, N. D., Yang, J., Zhang, X., Zhao, B., Grad, Y. H., Liu, L., et al. (2018). Clinically prevalent mutations in Mycobacterium tuberculosis alter propionate metabolism and mediate multidrug tolerance. Nat. Microbiol. 3 (9), 1032–1042. doi: 10.1038/s41564-018-0218-3
Hillemann, D., Rusch-Gerdes, S., Richter, E. (2008). In vitro-selected linezolid-resistant Mycobacterium tuberculosis mutants. Antimicrob. Agents Chemother. 52 (2), 800–801. doi: 10.1128/AAC.01189-07
Huang, L., Nazarova, E. V., Russell, D. G. (2019). Mycobacterium tuberculosis: Bacterial fitness within the host macrophage. Microbiol. Spectr. 7 (2), 1–12. doi: 10.1128/microbiolspec.BAI-0001-2019
Hughes, D., Andersson, D. I. (2017). Evolutionary trajectories to antibiotic resistance. Annu. Rev. Microbiol. 71, 579–596. doi: 10.1146/annurev-micro-090816-093813
Huitric, E., Verhasselt, P., Koul, A., Andries, K., Hoffner, S., Andersson, D. I. (2010). Rates and mechanisms of resistance development in Mycobacterium tuberculosis to a novel diarylquinoline ATP synthase inhibitor. Antimicrob. Agents Chemother. 54 (3), 1022–1028. doi: 10.1128/AAC.01611-09
Ince, D., Hooper, D. C. (2003). Quinolone resistance due to reduced target enzyme expression. J. Bacteriol 185 (23), 6883–6892. doi: 10.1128/JB.185.23.6883-6892.2003
Jakkala, K., Paul, A., Pradhan, A., Nair, R. R., Sharan, D., Swaminath, S., et al. (2020). Unique mode of cell division by the mycobacterial genetic resister clones emerging De novo from the antibiotic-surviving population. mSphere. 5 (6), 1–25. doi: 10.1128/mSphere.00994-20
Javid, B., Sorrentino, F., Toosky, M., Zheng, W., Pinkham, J. T., Jain, N., et al. (2014). Mycobacterial mistranslation is necessary and sufficient for rifampicin phenotypic resistance. Proc. Natl. Acad. Sci. U.S.A. 111 (3), 1132–1137. doi: 10.1073/pnas.1317580111
Kaniga, K., Lounis, N., Zhuo, S., Bakare, N., Andries, K. (2022). Impact of Rv0678 mutations on patients with drug-resistant TB treated with bedaquiline. Int. J. Tuberc Lung Dis. 26 (6), 571–573. doi: 10.5588/ijtld.21.0670
Katzir, I., Cokol, M., Aldridge, B. B., Alon, U. (2019). Prediction of ultra-high-order antibiotic combinations based on pairwise interactions. PloS Comput. Biol. 15 (1), e1006774. doi: 10.1371/journal.pcbi.1006774
Kayigire, X. A., Friedrich, S. O., van der Merwe, L., Diacon, A. H. (2017). Acquisition of rifampin resistance in pulmonary tuberculosis. Antimicrob. Agents Chemother. 61 (4), 1–5. doi: 10.1128/AAC.02220-16
Konno, K., Feldmann, F. M., McDermott, W. (1967). Pyrazinamide susceptibility and amidase activity of tubercle bacilli. Am. Rev. Respir. Dis. 95 (3), 461–469. doi: 10.1164/arrd.1967.95.3.461
Lee, J. J., Lee, S. K., Song, N., Nathan, T. O., Swarts, B. M., Eum, S. Y., et al. (2019). Transient drug-tolerance and permanent drug-resistance rely on the trehalose-catalytic shift in Mycobacterium tuberculosis. Nat. Commun. 10 (1), 2928. doi: 10.1038/s41467-019-10975-7
Lenaerts, A., Barry, C. E., 3rd, Dartois, V. (2015). Heterogeneity in tuberculosis pathology, microenvironments and therapeutic responses. Immunol. Rev. 264 (1), 288–307. doi: 10.1111/imr.12252
Levin-Reisman, I., Brauner, A., Ronin, I., Balaban, N. Q. (2019). Epistasis between antibiotic tolerance, persistence, and resistance mutations. Proc. Natl. Acad. Sci. U.S.A. 116 (29), 14734–14739. doi: 10.1073/pnas.1906169116
Levin-Reisman, I., Ronin, I., Gefen, O., Braniss, I., Shoresh, N., Balaban, N. Q. (2017). Antibiotic tolerance facilitates the evolution of resistance. Science 355 (6327), 826–830. doi: 10.1126/science.aaj2191
Levin, B. R., Rozen, D. E. (2006). Non-inherited antibiotic resistance. Nat. Rev. Microbiol. 4 (7), 556–562. doi: 10.1038/nrmicro1445
Li, Y. Y., Cai, R. J., Yang, J. Y., Hendrickson, T. L., Xiang, Y., Javid, B. (2021). Clinically relevant mutations of mycobacterial GatCAB inform regulation of translational fidelity. mBio 12 (4), e0110021. doi: 10.1128/mBio.01100-21
Li, Q.-j., Jiao, W.-w., Yin, Q.-q., Li, Y.-j., Li, J.-q., Xu, F., et al. (2017). Positive epistasis of major low-cost drug resistance mutations rpoB531-TTG and katG315-ACC depends on the phylogenetic background of Mycobacterium tuberculosis strains. Int. J. Antimicrobial Agents 49 (6), 757–762. doi: 10.1016/j.ijantimicag.2017.02.009
Lim, J., Lee, J. J., Lee, S. K., Kim, S., Eum, S. Y., Eoh, H. (2021). Phosphoenolpyruvate depletion mediates both growth arrest and drug tolerance of Mycobacterium tuberculosis in hypoxia. Proc. Natl. Acad. Sci. U.S.A. 118 (35), 1–12. doi: 10.1073/pnas.2105800118
Li, S., Poulton, N. C., Chang, J. S., Azadian, Z. A., DeJesus, M. A., Ruecker, N., et al. (2022). CRISPRi chemical genetics and comparative genomics identify genes mediating drug potency in Mycobacterium tuberculosis. Nat. Microbiol. 7 (6), 766–779. doi: 10.1038/s41564-022-01130-y
Liu, J., Shi, W., Zhang, S., Hao, X., Maslov, D. A., Shur, K. V., et al. (2019). Mutations in efflux pump Rv1258c (Tap) cause resistance to pyrazinamide, isoniazid, and streptomycin in M. tuberculosis. Front. Microbiol. 10. doi: 10.3389/fmicb.2019.00216
Liu, Q., Wei, J., Li, Y., Wang, M., Su, J., Lu, Y., et al. (2020). Mycobacterium tuberculosis clinical isolates carry mutational signatures of host immune environments. Sci. Adv. 6 (22), eaba4901. doi: 10.1126/sciadv.aba4901
Louw, G. E., Warren, R. M., Gey van Pittius, N. C., Leon, R., Jimenez, A., Hernandez-Pando, R., et al. (2011). Rifampicin reduces susceptibility to ofloxacin in rifampicin-resistant Mycobacterium tuberculosis through efflux. Am. J. Respir. Crit. Care Med. 184 (2), 269–276. doi: 10.1164/rccm.201011-1924OC
Luo, T., Yuan, J., Peng, X., Yang, G., Mi, Y., Sun, C., et al. (2017). Double mutation in DNA gyrase confers moxifloxacin resistance and decreased fitness of mycobacterium smegmatis. J. Antimicrobial Chemotherapy 72 (7), 1893–1900. doi: 10.1093/jac/dkx110
Machado, D., Couto, I., Perdigao, J., Rodrigues, L., Portugal, I., Baptista, P., et al. (2012). Contribution of efflux to the emergence of isoniazid and multidrug resistance in Mycobacterium tuberculosis. PloS One 7 (4), e34538. doi: 10.1371/journal.pone.0034538
Ma, S., Jaipalli, S., Larkins-Ford, J., Lohmiller, J., Aldridge, B. B., Sherman, D. R., et al. (2019). Transcriptomic signatures predict regulators of drug synergy and clinical regimen efficacy against tuberculosis. mBio 10 (6), 1–16. doi: 10.1128/mBio.02627-19
Ma, S., Morrison, R., Hobbs, S. J., Soni, V., Farrow-Johnson, J., Frando, A., et al. (2021). Transcriptional regulator-induced phenotype screen reveals drug potentiators in Mycobacterium tuberculosis. Nat. Microbiol. 6 (1), 44–50. doi: 10.1038/s41564-020-00810-x
Martini, M. C., Hicks, N. D., Xiao, J., Alonso, M. N., Barbier, T., Sixsmith, J., et al. (2022). Loss of RNase J leads to multi-drug tolerance and accumulation of highly structured mRNA fragments in Mycobacterium tuberculosis. PloS Pathog. 18 (7), e1010705. doi: 10.1371/journal.ppat.1010705
Maus, C. E., Plikaytis, B. B., Shinnick, T. M. (2005). Molecular analysis of cross-resistance to capreomycin, kanamycin, amikacin, and viomycin in Mycobacterium tuberculosis. Antimicrob. Agents Chemother. 49 (8), 3192–3197. doi: 10.1128/AAC.49.8.3192-3197.2005
McCune, R. M., Jr., Tompsett, R. (1956). Fate of Mycobacterium tuberculosis in mouse tissues as determined by the microbial enumeration technique. i. the persistence of drug-susceptible tubercle bacilli in the tissues despite prolonged antimicrobial therapy. J. Exp. Med. 104 (5), 737–762. doi: 10.1084/jem.104.5.737
McGrath, M., Pittius, N., Helden, P., Warren, R. M., Warner, D. F. (2014). Mutation rate and the emergence of drug resistance in Mycobacterium tuberculosis. J. Antimicrobial Chemotherapy. 69 (2), 292–302. doi: 10.1093/jac/dkt364
Meftahi, N., Namouchi, A., Mhenni, B., Brandis, G., Hughes, D., Mardassi, H. (2016). Evidence for the critical role of a secondary site rpoB mutation in the compensatory evolution and successful transmission of an MDR tuberculosis outbreak strain. J. Antimicrob. Chemother. 71 (2), 324–332. doi: 10.1093/jac/dkv345
Meier, A., Kirschner, P., Bange, F. C., Vogel, U., Bottger, E. C. (1994). Genetic alterations in streptomycin-resistant Mycobacterium tuberculosis: mapping of mutations conferring resistance. Antimicrob. Agents Chemother. 38 (2), 228–233. doi: 10.1128/AAC.38.2.228
Middlebrook, G. (1954). Isoniazid-resistance and catalase activity of tubercle bacilli; a preliminary report. Am. Rev. tuberculosis. 69 (3), 471–472. doi: 10.1164/art.1954.69.3.471
Molloy, A., Harrison, J., McGrath, J. S., Owen, Z., Smith, C., Liu, X., et al. (2021). Microfluidics as a novel technique for tuberculosis: From diagnostics to drug discovery. Microorganisms 9 (11), 1–24. doi: 10.3390/microorganisms9112330
Muzondiwa, D., Hlanze, H., Reva, O. N. (2021). The epistatic landscape of antibiotic resistance of different clades of Mycobacterium tuberculosis. Antibiotics (Basel) 10 (7), 1–14. doi: 10.3390/antibiotics10070857
Nair, J., Rouse, D. A., Bai, G. H., Morris, S. L. (1993). The rpsL gene and streptomycin resistance in single and multiple drug-resistant strains of Mycobacterium tuberculosis. Mol. Microbiol. 10 (3), 521–527. doi: 10.1111/j.1365-2958.1993.tb00924.x
Namouchi, A., Didelot, X., Schock, U., Gicquel, B., Rocha, E. P. (2012). After the bottleneck: Genome-wide diversification of the Mycobacterium tuberculosis complex by mutation, recombination, and natural selection. Genome Res. 22 (4), 721–734. doi: 10.1101/gr.129544.111
Nandakumar, M., Nathan, C., Rhee, K. Y. (2014). Isocitrate lyase mediates broad antibiotic tolerance in Mycobacterium tuberculosis. Nat. Commun. 5, 4306. doi: 10.1038/ncomms5306
Nieto Ramirez, L. M., Quintero Vargas, K., Diaz, G. (2020). Whole genome sequencing for the analysis of drug resistant strains of Mycobacterium tuberculosis: A systematic review for bedaquiline and delamanid. Antibiotics (Basel) 9 (3), 1–14. doi: 10.3390/antibiotics9030133
Nikaido, H. (2009). Multidrug resistance in bacteria. Biochemistry 78 (1), 119–146. doi: 10.1146/annurev.biochem.78.082907
Ohno, H., Koga, H., Kohno, S., Tashiro, T., Hara, K. (1996). Relationship between rifampin MICs for and rpoB mutations of Mycobacterium tuberculosis strains isolated in Japan. Antimicrob. Agents Chemother. 40 (4), 1053–1056. doi: 10.1128/AAC.40.4.1053
Patel, S. J., Lewis, B. E., Long, J. E., Nambi, S., Sassetti, C. M., Stemmler, T. L., et al. (2016). Fine-tuning of substrate affinity leads to alternative roles of Mycobacterium tuberculosis Fe2+-ATPases. J. Biol. Chem. 291 (22), 11529–11539. doi: 10.1074/jbc.M116.718239
Peterson, E. J. R., Ma, S., Sherman, D. R., Baliga, N. S. (2016). Network analysis identifies Rv0324 and Rv0880 as regulators of bedaquiline tolerance in Mycobacterium tuberculosis. Nat. Microbiol. 1 (8), nmicrobiol201678. doi: 10.1038/nmicrobiol.2016.78
Pi, R., Liu, Q., Takiff, H. E., Gao, Q. (2020). Fitness cost and compensatory evolution in levofloxacin-resistant Mycobacterium aurum. Antimicrob. Agents Chemother. 64 (8), 1–5. doi: 10.1128/AAC.00224-20
Pisu, D., Huang, L., Grenier, J. K., Russell, D. G. (2020). Dual RNA-seq of mtb-infected macrophages In vivo reveals ontologically distinct host-pathogen interactions. Cell Rep. 30 (2), 335–350.e334. doi: 10.1016/j.celrep.2019.12.033
Pisu, D., Huang, L., Narang, V., Theriault, M., Le-Bury, G., Lee, B., et al. (2021). Single cell analysis of M. tuberculosis phenotype and macrophage lineages in the infected lung. J. Exp. Med. 218 (9), 1–25. doi: 10.1084/jem.20210615
Plumlee, C. R., Duffy, F. J., Gern, B. H., Delahaye, J. L., Cohen, S. B., Stoltzfus, C. R., et al. (2021). Ultra-low dose aerosol infection of mice with Mycobacterium tuberculosis more closely models human tuberculosis. Cell Host Microbe 29 (1), 68–82.e65. doi: 10.1016/j.chom.2020.10.003
Pule, C. M., Sampson, S. L., Warren, R. M., Black, P. A., van Helden, P. D., Victor, T. C., et al. (2016). Efflux pump inhibitors: targeting mycobacterial efflux systems to enhance TB therapy. J. Antimicrob. Chemother. 71 (1), 17–26. doi: 10.1093/jac/dkv316
Pushparajan, A. R., Ramachandran, R., Gopi Reji, J., Ajay Kumar, R. (2020). Mycobacterium tuberculosis TetR family transcriptional regulator Rv1019 is a negative regulator of the mfd-mazG operon encoding DNA repair proteins. FEBS Lett. 594 (17), 2867–2880. doi: 10.1002/1873-3468.13861
Quinonez, C. G., Lee, J. J., Lim, J., Odell, M., Lawson, C. P., Anyogu, A., et al. (2022). The role of fatty acid metabolism in drug tolerance of Mycobacterium tuberculosis. mBio 13 (1), e0355921. doi: 10.1128/mbio.03559-21
Ragheb, M. N., Thomason, M. K., Hsu, C., Nugent, P., Gage, J., Samadpour, A. N., et al. (2019). Inhibiting the evolution of antibiotic resistance. Mol. Cell 73 (1), 157–165.e155. doi: 10.1016/j.molcel.2018.10.015
Raimunda, D., Long, J. E., Padilla-Benavides, T., Sassetti, C. M., Argüello, J. M. (2014). Differential roles for the Co(2+) /Ni(2+) transporting ATPases, CtpD and CtpJ, in Mycobacterium tuberculosis virulence. Mol. Microbiol. 91 (1), 185–197. doi: 10.1111/mmi.12454
Ramaswamy, S., Musser, J. M. (1998). Molecular genetic basis of antimicrobial agent resistance in Mycobacterium tuberculosis: 1998 update. Tuber Lung Dis. 79 (1), 3–29. doi: 10.1054/tuld.1998.0002
Rego, E. H., Audette, R. E., Rubin, E. J. (2017). Deletion of a mycobacterial divisome factor collapses single-cell phenotypic heterogeneity. Nature 546 (7656), 153–157. doi: 10.1038/nature22361
Reiche, M. A., Warner, D. F., Mizrahi, V. (2017). Targeting DNA replication and repair for the development of novel therapeutics against tuberculosis. Front. Mol. Biosci. 4. doi: 10.3389/fmolb.2017.00075
Rengarajan, J., Bloom, B. R., Rubin, E. J. (2005). Genome-wide requirements for Mycobacterium tuberculosis adaptation and survival in macrophages. Proc. Natl. Acad. Sci. U.S.A. 102 (23), 8327–8332. doi: 10.1073/pnas.0503272102
Rengarajan, J., Sassetti, C. M., Naroditskaya, V., Sloutsky, A., Bloom, B. R., Rubin, E. J. (2004). The folate pathway is a target for resistance to the drug para-aminosalicylic acid (PAS) in mycobacteria. Mol. Microbiol. 53 (1), 275–282. doi: 10.1111/j.1365-2958.2004.04120.x
Richards, J. P., Cai, W., Zill, N. A., Zhang, W., Ojha, A. K. (2019). Adaptation of Mycobacterium tuberculosis to biofilm growth is genetically linked to drug tolerance. Antimicrob. Agents Chemother. 63 (11), 1–17. doi: 10.1128/AAC.01213-19
Robicsek, A., Strahilevitz, J., Jacoby, G. A., Macielag, M., Abbanat, D., Park, C. H., et al. (2006). Fluoroquinolone-modifying enzyme: a new adaptation of a common aminoglycoside acetyltransferase. Nat. Med. 12 (1), 83–88. doi: 10.1038/nm1347
Rowlinson, M.-C., Musser, K. A. (2022). Current methods and role of next-generation sequencing in the diagnosis of antimicrobial resistance in tuberculosis. Clin. Microbiol. Newslett. 44 (1), 1–12. doi: 10.1016/j.clinmicnews.2021.12.001
Safi, H., Gopal, P., Lingaraju, S., Ma, S., Levine, C., Dartois, V., et al. (2019). Phase variation in Mycobacterium tuberculosis glpK produces transiently heritable drug tolerance. Proc. Natl. Acad. Sci. U.S.A. 116 (39), 19665–19674. doi: 10.1073/pnas.1907631116
Safi, H., Lingaraju, S., Amin, A., Kim, S., Jones, M., Holmes, M., et al. (2013). Evolution of high-level ethambutol-resistant tuberculosis through interacting mutations in decaprenylphosphoryl-beta-D-arabinose biosynthetic and utilization pathway genes. Nat. Genet. 45 (10), 1190–1197. doi: 10.1038/ng.2743
Safi, H., Lingaraju, S., Ma, S., Husain, S., Hoque, M., Soteropoulos, P., et al. (2020). Rapidly correcting frameshift mutations in the Mycobacterium tuberculosis orn gene produce reversible ethambutol resistance and small-Colony-Variant morphology. Antimicrob. Agents Chemother. 64 (9), 1–11. doi: 10.1128/AAC.00213-20
Sandegren, L., Andersson, D. I. (2009). Bacterial gene amplification: implications for the evolution of antibiotic resistance. Nat. Rev. Microbiol. 7 (8), 578–588. doi: 10.1038/nrmicro2174
Sander, P., Springer, B., Prammananan, T., Sturmfels, A., Kappler, M., Pletschette, M., et al. (2002). Fitness cost of chromosomal drug resistance-conferring mutations. Antimicrobial Agents Chemotherapy 46 (5), 1204–1211. doi: 10.1128/aac.46.5.1204-1211.2002
Santi, I., Dhar, N., Bousbaine, D., Wakamoto, Y., McKinney, J. D. (2013). Single-cell dynamics of the chromosome replication and cell division cycles in mycobacteria. Nat. Commun. 4, 2470. doi: 10.1038/ncomms3470
Sarathy, J. P., Via, L. E., Weiner, D., Blanc, L., Boshoff, H., Eugenin, E. A., et al. (2018). Extreme drug tolerance of Mycobacterium tuberculosis in caseum. Antimicrob. Agents Chemother. 62 (2), 1–11. doi: 10.1128/AAC.02266-17
Schaefer, W. B. (1954). The effect of isoniazid on growing and resting tubercle bacilli. Am. Rev. Tuberc 69 (1), 125–127. doi: 10.1164/art.1954.69.1.125
Schnappinger, D., Ehrt, S., Voskuil, M. I., Liu, Y., Mangan, J. A., Monahan, I. M., et al. (2003). Transcriptional adaptation of mycobacterium tuberculosis within macrophages: Insights into the phagosomal environment. J. Exp. Med. 198 (5), 693–704. doi: 10.1084/jem.20030846
Schubert, O. T., Ludwig, C., Kogadeeva, M., Zimmermann, M., Rosenberger, G., Gengenbacher, M., et al. (2015). Absolute proteome composition and dynamics during dormancy and resuscitation of Mycobacterium tuberculosis. Cell Host Microbe 18 (1), 96–108. doi: 10.1016/j.chom.2015.06.001
Scorpio, A., Zhang, Y. (1996). Mutations in pncA, a gene encoding pyrazinamidase/nicotinamidase, cause resistance to the antituberculous drug pyrazinamide in tubercle bacillus. Nat. Med. 2 (6), 662–667. doi: 10.1038/nm0696-662
Sebastian, J., Swaminath, S., Nair, R. R., Jakkala, K., Pradhan, A., Ajitkumar, P. (2017). De novo emergence of genetically resistant mutants of Mycobacterium tuberculosis from the persistence phase cells formed against antituberculosis drugs In vitro. Antimicrobial Agents Chemotherapy 61 (2), 1–25. doi: 10.1128/aac.01343-16
Shea, J., Halse, T. A., Kohlerschmidt, D., Lapierre, P., Modestil, H. A., Kearns, C. H., et al. (2021). Low-level rifampin resistance and rpoB mutations in Mycobacterium tuberculosis: an analysis of whole-genome sequencing and drug susceptibility test data in new York. J. Clin. Microbiol. 59 (4), 1–11. doi: 10.1128/JCM.01885-20
Sherman, D. R., Mdluli, K., Hickey, M. J., Arain, T. M., Morris, S. L., Iii, C. E. B., et al. (1996). Compensatory ahpC gene expression in isoniazid-resistant Mycobacterium tuberculosis. Science 272 (5268), 1641–1643. doi: 10.1126/science.272.5268.1641
Sirakova, T. D., Dubey, V. S., Deb, C., Daniel, J., Korotkova, T. A., Abomoelak, B., et al. (2006). Identification of a diacylglycerol acyltransferase gene involved in accumulation of triacylglycerol in Mycobacterium tuberculosis under stress. Microbiol. (Reading) 152 (Pt 9), 2717–2725. doi: 10.1099/mic.0.28993-0
Stevens, W. S., Scott, L., Noble, L., Gous, N., Dheda, K. (2017). Impact of the GeneXpert MTB/RIF technology on tuberculosis control. Microbiol. Spectr. 5 (1), 1–21. doi: 10.1128/microbiolspec.TBTB2-0040-2016
Strydom, N., Gupta, S. V., Fox, W. S., Via, L. E., Bang, H., Lee, M., et al. (2019). Tuberculosis drugs' distribution and emergence of resistance in patient's lung lesions: A mechanistic model and tool for regimen and dose optimization. PloS Med. 16 (4), e1002773. doi: 10.1371/journal.pmed.1002773
Sukheja, P., Kumar, P., Mittal, N., Li, S. G., Singleton, E., Russo, R., et al. (2017). A novel small-molecule inhibitor of the Mycobacterium tuberculosis demethylmenaquinone methyltransferase MenG is bactericidal to both growing and nutritionally deprived persister cells. mBio 8 (1), 1–21. doi: 10.1128/mBio.02022-16
Sun, H., Zeng, J., Li, S., Liang, P., Zheng, C., Liu, Y., et al. (2018). Interaction between rpsL and gyrA mutations affects the fitness and dual resistance of Mycobacterium tuberculosis clinical isolates against streptomycin and fluoroquinolones. Infection Drug Resistance 11, 431–440. doi: 10.2147/idr.s152335
Su, H. W., Zhu, J. H., Li, H., Cai, R. J., Ealand, C., Wang, X., et al. (2016). The essential mycobacterial amidotransferase GatCAB is a modulator of specific translational fidelity. Nat. Microbiol. 1 (11), 16147. doi: 10.1038/nmicrobiol.2016.147
Swaminath, S., Paul, A., Pradhan, A., Sebastian, J., Nair, R. R., Ajitkumar, P. (2020). Mycobacterium smegmatis moxifloxacin persister cells produce high levels of hydroxyl radical, generating genetic resisters selectable not only with moxifloxacin, but also with ethambutol and isoniazid. Microbiol. (Reading) 166 (2), 180–198. doi: 10.1099/mic.0.000874
Szumowski, J. D., Adams, K. N., Edelstein, P. H., Ramakrishnan, L. (2013). Antimicrobial efflux pumps and Mycobacterium tuberculosis drug tolerance: evolutionary considerations. Curr. Top. Microbiol. Immunol. 374, 81–108. doi: 10.1007/82_2012_300
Takiff, H. E., Salazar, L., Guerrero, C., Philipp, W., Huang, W. M., Kreiswirth, B., et al. (1994). Cloning and nucleotide sequence of Mycobacterium tuberculosis gyrA and gyrB genes and detection of quinolone resistance mutations. Antimicrob. Agents Chemother. 38 (4), 773–780. doi: 10.1128/AAC.38.4.773
Telenti, A., Imboden, P., Marchesi, F., Lowrie, D., Cole, S., Colston, M. J., et al. (1993). Detection of rifampicin-resistance mutations in Mycobacterium tuberculosis. Lancet 341 (8846), 647–650. doi: 10.1016/0140-6736(93)90417-f
Telenti, A., Philipp, W. J., Sreevatsan, S., Bernasconi, C., Stockbauer, K. E., Wieles, B., et al. (1997). The emb operon, a gene cluster of Mycobacterium tuberculosis involved in resistance to ethambutol. Nat. Med. 3 (5), 567–570. doi: 10.1038/nm0597-567
The CRyPTIC Consortium (2022). A data compendium associating the genomes of 12,289 Mycobacterium tuberculosis isolates with quantitative resistance phenotypes to 13 antibiotics. PloS Biol. 20 (8), e3001721. doi: 10.1371/journal.pbio.3001721
Trindade, S., Sousa, A., Xavier, K. B., Dionisio, F., Ferreira, M. G., Gordo, I. (2009). Positive epistasis drives the acquisition of multidrug resistance. PloS Genet. 5 (7), e1000578. doi: 10.1371/journal.pgen.1000578
Urbaniec, J., Xu, Y., Hu, Y., Hingley-Wilson, S., McFadden, J. (2022). Phenotypic heterogeneity in persisters: a novel 'hunker' theory of persistence. FEMS Microbiol. Rev. 46 (1), 1–16. doi: 10.1093/femsre/fuab042
Vargas, R., Freschi, L., Spitaleri, A., Tahseen, S., Barilar, I., Niemann, S., et al. (2021). Role of epistasis in amikacin, kanamycin, bedaquiline, and clofazimine resistance in Mycobacterium tuberculosis complex. Antimicrobial Agents Chemotherapy 65 (11), e01164-01121. doi: 10.1128/aac.01164-21
Vijay, S., Vinh, D. N., Hai, H. T., Ha, V. T. N., Dung, V. T. M., Dinh, T. D., et al. (2017). Influence of stress and antibiotic resistance on cell-length distribution in Mycobacterium tuberculosis clinical isolates. Front. Microbiol. 8. doi: 10.3389/fmicb.2017.02296
Villellas, C., Aristimuno, L., Vitoria, M. A., Prat, C., Blanco, S., Garcia de Viedma, D., et al. (2013). Analysis of mutations in streptomycin-resistant strains reveals a simple and reliable genetic marker for identification of the Mycobacterium tuberculosis Beijing genotype. J. Clin. Microbiol. 51 (7), 2124–2130. doi: 10.1128/JCM.01944-12
Villellas, C., Coeck, N., Meehan, C. J., Lounis, N., de Jong, B., Rigouts, L., et al. (2017). Unexpected high prevalence of resistance-associated Rv0678 variants in MDR-TB patients without documented prior use of clofazimine or bedaquiline. J. Antimicrob. Chemother. 72 (3), 684–690. doi: 10.1093/jac/dkw502
Wang, B. W., Zhu, J. H., Javid, B. (2020). Clinically relevant mutations in mycobacterial LepA cause rifampicin-specific phenotypic resistance. Sci. Rep. 10 (1), 8402. doi: 10.1038/s41598-020-65308-2
WHO (2021a). Catalogue of mutations in mycobacterium tuberculosis complex and their association with drug resistance (Geneva: World Health Organization). Licence: CC BY-NC-SA 3.0 IGO ed.
Widelock, D., Peizer, L. R., Klein, S. (1955). Public health significance of tubercle bacilli resistant to isoniazid. Am. J. Public Health Nations Health 45 (1), 79–83. doi: 10.2105/AJPH.45.1.79
Wolinsky, E., Smith, M. M., Steenken, W. (1956). Isoniazid susceptibility, catalase activity, and guinea pig virulence of recently isolated cultures of tubercle bacilli. Am. Rev. tuberculosis 73 (5), 768–772. doi: 10.1164/artpd.1956.73.5.768
Xie, Y. L., de Jager, V. R., Chen, R. Y., Dodd, L. E., Paripati, P., Via, L. E., et al. (2021). Fourteen-day PET/CT imaging to monitor drug combination activity in treated individuals with tuberculosis. Sci. Transl. Med. 13 (579), 1–11. doi: 10.1126/scitranslmed.abd7618
Xu, J., Li, S. Y., Almeida, D. V., Tasneen, R., Barnes-Boyle, K., Converse, P. J., et al. (2019). Contribution of pretomanid to novel regimens containing bedaquiline with either linezolid or moxifloxacin and pyrazinamide in murine models of tuberculosis. Antimicrob. Agents Chemother. 63 (5), 1–14. doi: 10.1128/AAC.00021-19
Xu, W., Snell, L. M., Guo, M., Boukhaled, G., Macleod, B. L., Li, M., et al. (2021). Early innate and adaptive immune perturbations determine long-term severity of chronic virus and Mycobacterium tuberculosis coinfection. Immunity 54 (3), 526–541.e527. doi: 10.1016/j.immuni.2021.01.003
Zeng, X., Kwok, J. S., Yang, K. Y., Leung, K. S., Shi, M., Yang, Z., et al. (2018). Whole genome sequencing data of 1110 Mycobacterium tuberculosis isolates identifies insertions and deletions associated with drug resistance. BMC Genomics 19 (1), 365. doi: 10.1186/s12864-018-4734-6
Zhang, Y., Heym, B., Allen, B., Young, D., Cole, S. (1992). The catalase-peroxidase gene and isoniazid resistance of Mycobacterium tuberculosis. Nature 358 (6387), 591–593. doi: 10.1038/358591a0
Zhang, H., Li, D., Zhao, L., Fleming, J., Lin, N., Wang, T., et al. (2013). Genome sequencing of 161 Mycobacterium tuberculosis isolates from China identifies genes and intergenic regions associated with drug resistance. Nat. Genet. 45 (10), 1255–1260. doi: 10.1038/ng.2735
Keywords: tuberculosis, mycobacteria, resistance, tolerance, heterogeneity, antibiotics
Citation: Jones RM, Adams KN, Eldesouky HE and Sherman DR (2022) The evolving biology of Mycobacterium tuberculosis drug resistance. Front. Cell. Infect. Microbiol. 12:1027394. doi: 10.3389/fcimb.2022.1027394
Received: 24 August 2022; Accepted: 20 September 2022;
Published: 05 October 2022.
Edited by:
Ben Gold, Weill Cornell Medicine, United StatesReviewed by:
Thomas Dick, Hackensack Meridian Health, United StatesAnna D. Tischler, University of Minnesota Twin Cities, United States
Alexandre Gouzy, NewYork-Presbyterian, United States
Copyright © 2022 Jones, Adams, Eldesouky and Sherman. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: David R. Sherman, ZHNoZXJtYW5AdXcuZWR1