Drug-susceptible Mycobacterium tuberculosis Beijing genotype does not develop mutation-conferred resistance to rifampin at an elevated rate

The Mycobacterium tuberculosis genome at 25 years: lessons and lingering questions

First achieved in 1998 by Cole et al., the complete genome sequence of Mycobacterium tuberculosis continues to provide an invaluable resource to understand tuberculosis (TB), the leading cause of global infectious disease mortality. At the 25-year anniversary of this accomplishment, we describe how insights gleaned from the M. tuberculosis genome have led to vital tools for TB research, epidemiology, and clinical practice. The increasing accessibility of whole-genome sequencing across research and clinical settings has improved our ability to predict antibacterial susceptibility, to track epidemics at the level of individual outbreaks and wider historical trends, to query the efficacy of the bacille Calmette-Guérin (BCG) vaccine, and to uncover targets for novel antitubercular therapeutics. Likewise, we discuss several recent efforts to extract further discoveries from this powerful resource.

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.

Estimating mutation rates in a Markov branching process using approximate Bayesian computation

Estimating microbial mutation rates is an essential task in evolutionary biology, with wide range applications in related fields such as virology, epidemiology, clinic and public health, and antibiotic research. Significant progress has been made on this research since 1943 when Luria-Delbrück fluctuation analysis was first introduced. However, existing estimators of mutation rates are heavily reliant on model assumptions in fluctuation analysis, and become less applicable to real microbial experiments which deviate from the model assumptions. To overcome this difficulty, we propose to model fluctuation experimental data by a two-type Markov branching process (MBP) and use approximate Bayesian computation (ABC) to estimate the mutation probability parameters. Such an ABC-based mutation rate estimator is based on intensive simulations from the mutation process, thereby taking advantage of modern computing power. Most importantly, its likelihood-free feature allows more complex and realistic setups of the mutation process, especially when the distribution of the number of mutants cannot be easily derived. To further improve computation efficiency, we use a Gaussian process surrogate to substitute the simulator in the ABC algorithm, and call the resulting estimator GPS-ABC. Simulation studies show that, when used to estimate constant mutation rate in MBP, ABC-based estimators generally outperform traditional moment or likelihood-based estimators. When mutations occur in two stages, i.e., in MBP with a piece-wise constant mutation rate function, traditional mutation rate estimators become not applicable, yet GPS-ABC still achieves reasonable estimates. Finally, the proposed GPS-ABC estimator is used to analyze real fluctuation experimental datasets for studying drug resistance.

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

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

Emergence of phenotypic and genotypic antimicrobial resistance in Mycobacterium tuberculosis

Concentration dependency of phenotypic and genotypic isoniazid-rifampicin resistance emergence was investigated to obtain a mechanistic understanding on how anti-mycobacterial drugs facilitate the emergence of bacterial populations that survive throughout treatment. Using static kill curve experiments, observing two evolution cycles, it was demonstrated that rifampicin resistance was the result of non-specific mechanisms and not associated with accumulation of drug resistance encoding SNPs. Whereas, part of isoniazid resistance could be accounted for by accumulation of specific SNPs, which was concentration dependent. Using a Hollow Fibre Infection Model it was demonstrated that emergence of resistance did not occur at concentration-time profiles mimicking the granuloma. This study showed that disentangling and quantifying concentration dependent emergence of resistance provides an improved rational for drug and dose selection although further work to understand the underlying mechanisms is needed to improve the drug development pipeline.

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.

Molecular characteristics of Mycobacterium tuberculosis drug-resistant isolates from HIV- and HIV+ tuberculosis patients in Russia

Background

High burden of drug-resistant (DR) tuberculosis (TB) is a significant threat to national TB control programs all over the world and in the Russian Federation. Different Mycobacterium tuberculosis (MTB) genotypes are hypothesized to have specific characteristics affecting TB control programs. For example, Beijing strains are supposed to have higher mutation rates compared to strains of other genotypes and subsequently higher capability to develop drug-resistance.

Results

Clinical MTB isolates from HIV- and HIV+ patients from four regions of Russia were analyzed for genotypes and mutations conferring resistance to Isoniazid, Rifampicin, Ethambutol, aminoglycosides, and fluoroquinolones. Analysis of genotypes and polymorphism of genomic loci according to the HIV status of the patients – sources of MTB isolates were performed. Studied MTB isolates from HIV- TB patients belonged to 15 genotypes and from HIV + TB patients – to 6 genotypes. Beijing clinical isolates dominated in HIV- (64,7%) and HIV+ (74,4%) groups. Other isolates were of LAM (including LAM1 and LAM9), Ural, and 4 minor groups of genotypes (including 5 subclones T). The spectrum of genotypes in the HIV- group was broader than in the HIV+ group. PR of B0/W148 Beijing was significantly lower than of other Beijing genotypes in susceptible and MDR-XDR isolates. Rates of isolates belonging to non-Beijing genotypes were higher than Beijing in susceptible isolates from HIV- patients.

Conclusions

Beijing genotype isolates prevailed in clinical isolates of all drug susceptibility profiles both from HIV- and HIV+ patients, although B0/W148 Beijing genotype did not dominate in this study. Genome loci and mutations polymorphisms were more pronounced in clinical isolates from HIV- patients, than from HIV+.

Genomic signatures of pre-resistance in Mycobacterium tuberculosis

Recent advances in bacterial whole-genome sequencing have resulted in a comprehensive catalog of antibiotic resistance genomic signatures in Mycobacterium tuberculosis. With a view to pre-empt the emergence of resistance, we hypothesized that pre-existing polymorphisms in susceptible genotypes (pre-resistance mutations) could increase the risk of becoming resistant in the future. We sequenced whole genomes from 3135 isolates sampled over a 17-year period. After reconstructing ancestral genomes on time-calibrated phylogenetic trees, we developed and applied a genome-wide survival analysis to determine the hazard of resistance acquisition. We demonstrate that M. tuberculosis lineage 2 has a higher risk of acquiring resistance than lineage 4, and estimate a higher hazard of rifampicin resistance evolution following isoniazid mono-resistance. Furthermore, we describe loci and genomic polymorphisms associated with a higher risk of resistance acquisition. Identifying markers of future antibiotic resistance could enable targeted therapy to prevent resistance emergence in M. tuberculosis and other pathogens.

The within-host evolution of antimicrobial resistance in Mycobacterium tuberculosis

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

Control of Genome Stability by EndoMS/NucS-Mediated Non-Canonical Mismatch Repair

The DNA repair endonuclease EndoMS/NucS is highly conserved in Archaea and Actinobacteria. This enzyme is able to recognize and cleave dsDNA carrying a mismatched base pair, and its activity is enhanced by the interaction with the sliding clamp of the replisome. Today, EndoMS/NucS has been established as the key protein of a non-canonical mismatch repair (MMR) pathway, acting specifically in the repair of transitions and being essential for maintaining genome stability. Despite having some particularities, such as its lower activity on transversions and the inability to correct indels, EndoMS/NucS meets the main hallmarks of a MMR. Its absence leads to a hypermutator phenotype, a transition-biased mutational spectrum and an increase in homeologous recombination. Interestingly, polymorphic EndoMS/NucS variants with a possible effect in mutation rate have been detected in clinical isolates of the relevant actinobacterial pathogen Mycobacterium tuberculosis. Considering that MMR defects are often associated with the emergence of resistant bacteria, the existence of EndoMS/NucS-defective mutators could have an important role in the acquisition of antibiotic resistance in M. tuberculosis. Therefore, a further understanding of the EndoMS/NucS-mediated non-canonical MMR pathway may reveal new strategies to predict and fight drug resistance. This review is focused on the recent progress in NucS, with special emphasis on its effect on genome stability and evolvability in Actinobacteria.

Mechanisms and detection methods of Mycobacterium tuberculosis rifampicin resistance: The phenomenon of drug resistance is complex

Tuberculosis (TB) is an infectious disease that poses a serious threat to human health. Rifampin (RIF) is an important first-line anti-TB drug, and rifampin resistance (RIF-R) is a key factor in formulating treatment regimen and evaluating the prognosis of TB. Compared with other drugs resistance, the RIF-R mechanism of Mycobacterium tuberculosis (M. tuberculosis) is one of the clearest, which is mainly caused by RIF resistance-related mutations in the rpoB gene. This provides a convenient condition for developing rapid detection methods, and also an ideal object for studying the general drug resistance mechanisms of M. tuberculosis. This review focuses on the mechanisms that influence the RIF resistance of M. tuberculosis and related detection methods. Besides the mutations in rpoB, M. tuberculosis can decrease the amount of drugs entering the cells, enhance the drugs efflux, and be heterogeneous RIF susceptibility to resist drug pressure. Based on the results of current researches, many genes participate in influencing the susceptibility to RIF, which indicates the phenomenon of M. tuberculosis drug resistance is very complex.

New approaches to mutation rate fold change in Luria-Delbrück fluctuation experiments

For nearly eight decades the Luria-Delbrück protocol remains the preferred method for experimentally determining microbial mutation rates. However, earnest development and rigorous applications of statistical methods for mutation rate comparison using fluctuation assay data are a relatively recent phenomenon. While likelihood ratio tests tailored for the fluctuation protocol give investigators appropriate tools, researchers sometimes may prefer to view the comparison of two mutation rates through the lens of the ratio of the two mutation rates. The idea of using the bootstrap technique to construct intervals for mutation rate fold change was proposed nearly a decade ago, but it failed to gain traction partly due to a failure to incorporate likelihood-based estimation. In addition to extending the bootstrap method, this paper proposes two new methods of constructing intervals for mutation rate fold change: a profile likelihood method and a Bayesian method utilizing Monte Carlo Markov chain. All three methods are assessed by large-scale simulations.

The Genetic Background Modulates the Evolution of Fluoroquinolone-Resistance in Mycobacterium tuberculosis

Fluoroquinolones (FQ) form the backbone in experimental treatment regimens against drug-susceptible tuberculosis. However, little is known on whether the genetic variation present in natural populations of Mycobacterium tuberculosis (Mtb) affects the evolution of FQ-resistance (FQ-R). To investigate this question, we used nine genetically distinct drug-susceptible clinical isolates of Mtb and measured their frequency of resistance to the FQ ofloxacin (OFX) in vitro. We found that the Mtb genetic background led to differences in the frequency of OFX-resistance (OFX-R) that spanned two orders of magnitude and substantially modulated the observed mutational profiles for OFX-R. Further, in vitro assays showed that the genetic background also influenced the minimum inhibitory concentration and the fitness effect conferred by a given OFX-R mutation. To test the clinical relevance of our in vitro work, we surveyed the mutational profile for FQ-R in publicly available genomic sequences from clinical Mtb isolates, and found substantial Mtb lineage-dependent variability. Comparison of the clinical and the in vitro mutational profiles for FQ-R showed that 51% and 39% of the variability in the clinical frequency of FQ-R gyrA mutation events in Lineage 2 and Lineage 4 strains, respectively, can be attributed to how Mtb evolves FQ-R in vitro. As the Mtb genetic background strongly influenced the evolution of FQ-R in vitro, we conclude that the genetic background of Mtb also impacts the evolution of FQ-R in the clinic.

Multifidelity Analysis for Predicting Rare Events in Stochastic Computational Models of Complex Biological Systems

Rare events such as genetic mutations or cell-cell interactions are important contributors to dynamics in complex biological systems, eg, in drug-resistant infections. Computational approaches can help analyze rare events that are difficult to study experimentally. However, analyzing the frequency and dynamics of rare events in computational models can also be challenging due to high computational resource demands, especially for high-fidelity stochastic computational models. To facilitate analysis of rare events in complex biological systems, we present a multifidelity analysis approach that uses medium-fidelity analysis (Monte Carlo simulations) and/or low-fidelity analysis (Markov chain models) to analyze high-fidelity stochastic model results. Medium-fidelity analysis can produce large numbers of possible rare event trajectories for a single high-fidelity model simulation. This allows prediction of both rare event dynamics and probability distributions at much lower frequencies than high-fidelity models. Low-fidelity analysis can calculate probability distributions for rare events over time for any frequency by updating the probabilities of the rare event state space after each discrete event of the high-fidelity model. To validate the approach, we apply multifidelity analysis to a high-fidelity model of tuberculosis disease. We validate the method against high-fidelity model results and illustrate the application of multifidelity analysis in predicting rare event trajectories, performing sensitivity analyses and extrapolating predictions to very low frequencies in complex systems. We believe that our approach will complement ongoing efforts to enable accurate prediction of rare event dynamics in high-fidelity computational models.

Emergence and selection of isoniazid and rifampin resistance in tuberculosis granulomas

Drug resistant tuberculosis is increasing world-wide. Resistance against isoniazid (INH), rifampicin (RIF), or both (multi-drug resistant TB, MDR-TB) is of particular concern, since INH and RIF form part of the standard regimen for TB disease. While it is known that suboptimal treatment can lead to resistance, it remains unclear how host immune responses and antibiotic dynamics within granulomas (sites of infection) affect emergence and selection of drug-resistant bacteria. We take a systems pharmacology approach to explore resistance dynamics within granulomas. We integrate spatio-temporal host immunity, INH and RIF dynamics, and bacterial dynamics (including fitness costs and compensatory mutations) in a computational framework. We simulate resistance emergence in the absence of treatment, as well as resistance selection during INH and/or RIF treatment. There are four main findings. First, in the absence of treatment, the percentage of granulomas containing resistant bacteria mirrors the non-monotonic bacterial dynamics within granulomas. Second, drug-resistant bacteria are less frequently found in non-replicating states in caseum, compared to drug-sensitive bacteria. Third, due to a steeper dose response curve and faster plasma clearance of INH compared to RIF, INH-resistant bacteria have a stronger influence on treatment outcomes than RIF-resistant bacteria. Finally, under combination therapy with INH and RIF, few MDR bacteria are able to significantly affect treatment outcomes. Overall, our approach allows drug-specific prediction of drug resistance emergence and selection in the complex granuloma context. Since our predictions are based on pre-clinical data, our approach can be implemented relatively early in the treatment development process, thereby enabling pro-active rather than reactive responses to emerging drug resistance for new drugs. Furthermore, this quantitative and drug-specific approach can help identify drug-specific properties that influence resistance and use this information to design treatment regimens that minimize resistance selection and expand the useful life-span of new antibiotics.

DNA Replication Fidelity in the Mycobacterium tuberculosis Complex

Mycobacterium tuberculosis is genetically isolated, with no evidence for horizontal gene transfer or the acquisition of episomal genetic information in the modern evolution of strains of the Mycobacterium tuberculosis complex. When considered in the context of the specific features of the disease M. tuberculosis causes (e.g., transmission via cough aerosol, replication within professional phagocytes, subclinical persistence, and stimulation of a destructive immune pathology), this implies that to understand the mechanisms ensuring preservation of genomic integrity in infecting mycobacterial populations is to understand the source of genetic variation, including the emergence of microdiverse sub-populations that may be linked to the acquisition of drug resistance. In this chapter, we focus on mechanisms involved in maintaining DNA replication fidelity in M. tuberculosis, and consider the potential to target components of the DNA replication machinery as part of novel therapeutic regimens designed to curb the emerging threat of drug-resistance.

Antimicrobial resistance in Mycobacterium tuberculosis: mechanistic and evolutionary perspectives

Antibiotic-resistant Mycobacterium tuberculosis strains are threatening progress in containing the global tuberculosis epidemic. Mycobacterium tuberculosis is intrinsically resistant to many antibiotics, limiting the number of compounds available for treatment. This intrinsic resistance is due to a number of mechanisms including a thick, waxy, hydrophobic cell envelope and the presence of drug degrading and modifying enzymes. Resistance to the drugs which are active against M. tuberculosis is, in the absence of horizontally transferred resistance determinants, conferred by chromosomal mutations. These chromosomal mutations may confer drug resistance via modification or overexpression of the drug target, as well as by prevention of prodrug activation. Drug resistance mutations may have pleiotropic effects leading to a reduction in the bacterium's fitness, quantifiable e.g. by a reduction in the in vitro growth rate. Secondary so-called compensatory mutations, not involved in conferring resistance, can ameliorate the fitness cost by interacting epistatically with the resistance mutation. Although the genetic diversity of M. tuberculosis is low compared to other pathogenic bacteria, the strain genetic background has been demonstrated to influence multiple aspects in the evolution of drug resistance. The rate of resistance evolution and the fitness costs of drug resistance mutations may vary as a function of the genetic background.

rSalvador: An R Package for the Fluctuation Experiment

The past few years have seen a surge of novel applications of the Luria-Delbrück fluctuation assay protocol in bacterial research. Appropriate analysis of fluctuation assay data often requires computational methods that are unavailable in the popular web tool FALCOR. This paper introduces an R package named rSalvador to bring improvements to the field. The paper focuses on rSalvador’s capabilities to alleviate three kinds of problems found in recent investigations: (i) resorting to partial plating without properly accounting for the effects of partial plating; (ii) conducting attendant fitness assays without incorporating mutants’ relative fitness in subsequent data analysis; and (iii) comparing mutation rates using methods that are in general inapplicable to fluctuation assay data. In addition, the paper touches on rSalvador’s capabilities to estimate sample size and the difficulties related to parameter nonidentifiability.

Impact of Genetic Diversity on the Biology of Mycobacterium tuberculosis Complex Strains

Tuberculosis (TB) remains the most deadly bacterial infectious disease worldwide. Its treatment and control are threatened by increasing numbers of multidrug-resistant (MDR) or nearly untreatable extensively drug-resistant (XDR) strains. New concepts are therefore urgently needed to understand the factors driving the TB epidemics and the spread of different strain populations, especially in association with drug resistance. Classical genotyping and, more recently, whole-genome sequencing (WGS) revealed that the world population of tubercle bacilli is more diverse than previously thought. Several major phylogenetic lineages can be distinguished, which are associated with their sympatric host population. Distinct clonal (sub)populations can even coexist within infected patients. WGS is now used as the ultimate approach for differentiating clinical isolates and for linking phenotypic to genomic variation from lineage to strain levels. Multiple lines of evidence indicate that the genetic diversity of TB strains translates into pathobiological consequences, and key molecular mechanisms probably involved in differential pathoadaptation of some main lineages have recently been identified. Evidence also accumulates on molecular mechanisms putatively fostering the emergence and rapid expansion of particular MDR and XDR strain groups in some world regions. However, further integrative studies will be needed for complete elucidation of the mechanisms that allow the pathogen to infect its host, acquire multidrug resistance, and transmit so efficiently. Such knowledge will be key for the development of the most effective new diagnostics, drugs, and vaccination strategies.

DNA Replication in Mycobacterium tuberculosis

Faithful replication and maintenance of the genome are essential to the ability of any organism to survive and propagate. For an obligate pathogen such as Mycobacterium tuberculosis that has to complete successive cycles of transmission, infection, and disease in order to retain a foothold in the human population, this requires that genome replication and maintenance must be accomplished under the metabolic, immune, and antibiotic stresses encountered during passage through variable host environments. Comparative genomic analyses have established that chromosomal mutations enable M. tuberculosis to adapt to these stresses: the emergence of drug-resistant isolates provides direct evidence of this capacity, so too the well-documented genetic diversity among M. tuberculosis lineages across geographic loci, as well as the microvariation within individual patients that is increasingly observed as whole-genome sequencing methodologies are applied to clinical samples and tuberculosis (TB) disease models. However, the precise mutagenic mechanisms responsible for M. tuberculosis evolution and adaptation are poorly understood. Here, we summarize current knowledge of the machinery responsible for DNA replication in M. tuberculosis, and discuss the potential contribution of the expanded complement of mycobacterial DNA polymerases to mutagenesis. We also consider briefly the possible role of DNA replication-in particular, its regulation and coordination with cell division-in the ability of M. tuberculosis to withstand antibacterial stresses, including host immune effectors and antibiotics, through the generation at the population level of a tolerant state, or through the formation of a subpopulation of persister bacilli-both of which might be relevant to the emergence and fixation of genetic drug resistance.

Armed conflict and population displacement as drivers of the evolution and dispersal of Mycobacterium tuberculosis

The "Beijing" Mycobacterium tuberculosis (Mtb) lineage 2 (L2) is spreading globally and has been associated with accelerated disease progression and increased antibiotic resistance. Here we performed a phylodynamic reconstruction of one of the L2 sublineages, the central Asian clade (CAC), which has recently spread to western Europe. We find that recent historical events have contributed to the evolution and dispersal of the CAC. Our timing estimates indicate that the clade was likely introduced to Afghanistan during the 1979-1989 Soviet-Afghan war and spread further after population displacement in the wake of the American invasion in 2001. We also find that drug resistance mutations accumulated on a massive scale in Mtb isolates from former Soviet republics after the fall of the Soviet Union, a pattern that was not observed in CAC isolates from Afghanistan. Our results underscore the detrimental effects of political instability and population displacement on tuberculosis control and demonstrate the power of phylodynamic methods in exploring bacterial evolution in space and time.

A second look at the final number of cells in a fluctuation experiment

In a fluctuation experiment, the number of cells existing in a culture immediately before plating (commonly known as Nt) varies across the parallel cultures. However, most existing mathematical models for fluctuation assay data do not recognize the variation in Nt. Despite repeated attempts in the past to integrate this source of variability in the estimation of microbial mutation rates, several questions of practical importance remain unanswered. The present investigation finds that the variation needs accounting for only when the coefficient of variation for Nt is large, and experimental data suggest that the coefficient of variation is often moderate or small. Moreover, an increase in the inoculum size can reduce the coefficient of variation. Through extensive simulation, several existing methods that accommodate the variation in Nt are compared. It was found that a newly devised likelihood method based on the existing gamma mixture model outperforms other existing methods. The investigation focuses on the estimation of mutation rates using the Lea-Coulson model, under which mutation is selectively neutral; however, the paper also explores the major findings' implications for the comparison of mutation rates using the likelihood ratio test, and for the estimation of mutation rates using the Mandelbrot-Koch model that allows for non-neutral mutations.

Antimicrobial Resistance in Mycobacterium tuberculosis: The Odd One Out

Antimicrobial resistance (AMR) threats are typically represented by bacteria capable of extensive horizontal gene transfer (HGT). One clear exception is Mycobacterium tuberculosis (Mtb). It is an obligate human pathogen with limited genetic diversity and a low mutation rate which lacks any evidence for HGT. Such features should, in principle, reduce its ability to rapidly evolve AMR. We identify key features in its biology and epidemiology that allow it to overcome its low adaptive potential. We focus in particular on its innate resistance to drugs, its unusual life cycle, including an often extensive latent phase, and its ability to shelter from exposure to antimicrobial drugs within cavities it induces in the lungs.

Genetic features of Mycobacterium tuberculosis modern Beijing sublineage

Mycobacterium tuberculosis (MTB) Beijing strains have caused a great concern because of their rapid emergence and increasing prevalence in worldwide regions. Great efforts have been made to investigate the pathogenic characteristics of Beijing strains such as hypervirulence, drug resistance and favoring transmission. Phylogenetically, MTB Beijing family was divided into modern and ancient sublineages. Modern Beijing strains displayed enhanced virulence and higher prevalence when compared with ancient Beijing strains, but the genetic basis for this difference remains unclear. In this study, by analyzing previously published sequencing data of 1082 MTB Beijing isolates, we determined the genetic changes that were commonly present in modern Beijing strains but absent in ancient Beijing strains. These changes include 44 single-nucleotide polymorphisms (SNPs) and two short genomic deletions. Through bioinformatics analysis, we demonstrated that these genetic changes had high probability of functional effects. For example, 4 genes were frameshifted due to premature stop mutation or genomic deletions, 19 nonsynonymous SNPs located in conservative codons, and there is a significant enrichment in regulatory network for all nonsynonymous mutations. Besides, three SNPs located in promoter regions were verified to alter downstream gene expressions. Our study precisely defined the genetic features of modern Beijing strains and provided interesting clues for future researches to elucidate the mechanisms that underlie this sublineage's successful expansion. These findings from the analysis of the modern Beijing sublineage could provide us a model to understand the dynamics of pathogenicity of MTB.

A new practical guide to the Luria-Delbrück protocol

Since 2000 several review papers have been published about the analysis of experimental data obtained using the Luria-Delbrück protocol. These timely papers cleared much of the confusion surrounding various methods for estimating or comparing mutation rates. As a result, today the fluctuation test is more widely applied with much improved accuracy. The present paper provides guidelines on a few remaining problems that continue to baffle mutation researchers. Among the issues addressed are incomplete plating, relative fitness, and comparison of experiments where average final cell population sizes differ. It also offers a fresh view on the estimation methods that are based on the sample median.

Clinical value of whole-genome sequencing of Mycobacterium tuberculosis

Whole-genome sequencing (WGS) is now common as a result of new technologies that can rapidly sequence a complete bacterial genome for US$500 or less. Many studies have addressed questions about tuberculosis with WGS, and knowing the sequence of the entire genome, rather than only a few fragments, has greatly increased the precision of molecular epidemiology and contact tracing. Additionally, topics such as the mutation rate, drug resistance, the target of new drugs, and the phylogeny and evolution of the Mycobacterium tuberculosis complex bacteria have been elucidated by WGS. Nonetheless, WGS has not explained differences in transmissibility between strains, or why some strains are more virulent than others or more prone to development of multidrug resistance. With advances in technology, WGS of clinical specimens could become routine in high-income countries; however, its relevance will probably depend on easy to use software to efficiently process the sequences produced and accessible genomic databases that can be mined in future studies.

Clinical application of whole-genome sequencing to inform treatment for multidrug-resistant tuberculosis cases

The treatment of drug-resistant tuberculosis cases is challenging, as drug options are limited, and the existing diagnostics are inadequate. Whole-genome sequencing (WGS) has been used in a clinical setting to investigate six cases of suspected extensively drug-resistant Mycobacterium tuberculosis (XDR-TB) encountered at a London teaching hospital between 2008 and 2014. Sixteen isolates from six suspected XDR-TB cases were sequenced; five cases were analyzed in a clinically relevant time frame, with one case sequenced retrospectively. WGS identified mutations in the M. tuberculosis genes associated with antibiotic resistance that are likely to be responsible for the phenotypic resistance. Thus, an evidence base was developed to inform the clinical decisions made around antibiotic treatment over prolonged periods. All strains in this study belonged to the East Asian (Beijing) lineage, and the strain relatedness was consistent with the expectations from the case histories, confirming one contact transmission event. We demonstrate that WGS data can be produced in a clinically relevant time scale some weeks before drug sensitivity testing (DST) data are available, and they actively help clinical decision-making through the assessment of whether an isolate (i) has a particular resistance mutation where there are absent or contradictory DST results, (ii) has no further resistance markers and therefore is unlikely to be XDR, or (iii) is identical to an isolate of known resistance (i.e., a likely transmission event). A small number of discrepancies between the genotypic predictions and phenotypic DST results are discussed in the wider context of the interpretation and reporting of WGS results.

Diversity and evolution of Mycobacterium tuberculosis: moving to whole-genome-based approaches

Genotyping of clinical Mycobacterium tuberculosis complex (MTBC) strains has become a standard tool for epidemiological tracing and for the investigation of the local and global strain population structure. Of special importance is the analysis of the expansion of multidrug (MDR) and extensively drug-resistant (XDR) strains. Classical genotyping and, more recently, whole-genome sequencing have revealed that the strains of the MTBC are more diverse than previously anticipated. Globally, several phylogenetic lineages can be distinguished whose geographical distribution is markedly variable. Strains of particular (sub)lineages, such as Beijing, seem to be more virulent and associated with enhanced resistance levels and fitness, likely fueling their spread in certain world regions. The upcoming generalization of whole-genome sequencing approaches will expectedly provide more comprehensive insights into the molecular and epidemiological mechanisms involved and lead to better diagnostic and therapeutic tools.

Unbiased estimation of mutation rates under fluctuating final counts

Estimation methods for mutation rates (or probabilities) in Luria-Delbrück fluctuation analysis usually assume that the final number of cells remains constant from one culture to another. We show that this leads to systematically underestimate the mutation rate. Two levels of information on final numbers are considered: either the coefficient of variation has been independently estimated, or the final number of cells in each culture is known. In both cases, unbiased estimation methods are proposed. Their statistical properties are assessed both theoretically and through Monte-Carlo simulation. As an application, the data from two well known fluctuation analysis studies on Mycobacterium tuberculosis are reexamined.

Whole genome sequencing reveals complex evolution patterns of multidrug-resistant Mycobacterium tuberculosis Beijing strains in patients

Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains represent a major threat for tuberculosis (TB) control. Treatment of MDR-TB patients is long and less effective, resulting in a significant number of treatment failures. The development of further resistances leads to extensively drug-resistant (XDR) variants. However, data on the individual reasons for treatment failure, e.g. an induced mutational burst, and on the evolution of bacteria in the patient are only sparsely available. To address this question, we investigated the intra-patient evolution of serial MTBC isolates obtained from three MDR-TB patients undergoing longitudinal treatment, finally leading to XDR-TB. Sequential isolates displayed identical IS6110 fingerprint patterns, suggesting the absence of exogenous re-infection. We utilized whole genome sequencing (WGS) to screen for variations in three isolates from Patient A and four isolates from Patient B and C, respectively. Acquired polymorphisms were subsequently validated in up to 15 serial isolates by Sanger sequencing. We determined eight (Patient A) and nine (Patient B) polymorphisms, which occurred in a stepwise manner during the course of the therapy and were linked to resistance or a potential compensatory mechanism. For both patients, our analysis revealed the long-term co-existence of clonal subpopulations that displayed different drug resistance allele combinations. Out of these, the most resistant clone was fixed in the population. In contrast, baseline and follow-up isolates of Patient C were distinguished each by eleven unique polymorphisms, indicating an exogenous re-infection with an XDR strain not detected by IS6110 RFLP typing. Our study demonstrates that intra-patient microevolution of MDR-MTBC strains under longitudinal treatment is more complex than previously anticipated. However, a mutator phenotype was not detected. The presence of different subpopulations might confound phenotypic and molecular drug resistance tests. Furthermore, high resolution WGS analysis is necessary to accurately detect exogenous re-infection as classical genotyping lacks discriminatory power in high incidence settings.

Mycobacterium tuberculosis mutation rate estimates from different lineages predict substantial differences in the emergence of drug-resistant tuberculosis

A critical question in tuberculosis control is why some strains of Mycobacterium tuberculosis are preferentially associated with multiple drug resistances. We demonstrate that M. tuberculosis strains from Lineage 2 (East Asian lineage and Beijing sublineage) acquire drug resistances in vitro more rapidly than M. tuberculosis strains from Lineage 4 (Euro-American lineage) and that this higher rate can be attributed to a higher mutation rate. Moreover, the in vitro mutation rate correlates well with the bacterial mutation rate in humans as determined by whole genome sequencing of clinical isolates. Finally, using a stochastic mathematical model, we demonstrate that the observed differences in mutation rate predict a substantially higher probability that patients infected with a drug susceptible Lineage 2 strain will harbor multidrug resistant bacteria at the time of diagnosis. These data suggest that interventions to prevent the emergence of drug resistant tuberculosis should target bacterial as well as treatment-related risk factors.

DNA metabolism in mycobacterial pathogenesis

Fundamental aspects of the lifestyle of Mycobacterium tuberculosis implicate DNA metabolism in bacillary survival and adaptive evolution. The environments encountered by M. tuberculosis during successive cycles of infection and transmission are genotoxic. Moreover, as an obligate pathogen, M. tuberculosis has the ability to persist for extended periods in a subclinical state, suggesting that active DNA repair is critical to maintain genome integrity and bacterial viability during prolonged infection. In this chapter, we provide an overview of the major DNA metabolic pathways identified in M. tuberculosis, and situate key recent findings within the context of mycobacterial pathogenesis. Unlike many other bacterial pathogens, M. tuberculosis is genetically secluded, and appears to rely solely on chromosomal mutagenesis to drive its microevolution within the human host. In turn, this implies that a balance between high versus relaxed fidelity mechanisms of DNA metabolism ensures the maintenance of genome integrity, while accommodating the evolutionary imperative to adapt to hostile and fluctuating environments. The inferred relationship between mycobacterial DNA repair and genome dynamics is considered in the light of emerging data from whole-genome sequencing studies of clinical M. tuberculosis isolates which have revealed the potential for considerable heterogeneity within and between different bacterial and host populations.

Genetic diversity in Mycobacterium tuberculosis

Recent years have witnessed an increased appreciation of the extent and relevance of strain-to-strain variation in Mycobacterium tuberculosis. This paradigm shift can largely be attributed to an improved understanding of the global population structure of this organism, and to the realisation that the various members of the M. tuberculosis complex (MTBC) harbour more genetic diversity than previously realised. Moreover, many studies using experimental models of infection have demonstrated that MTBC diversity translates into significant differences in immunogenecity and virulence . However, linking these experimental phenotypes to relevant clinical phenotypes has been difficult, and to date, largely unsuccessful. Nevertheless, emerging high-throughput technologies, in particular next-generation sequencing , offer new opportunities, and have already lead to important new insights. Given the complexity of the host-pathogen interaction in tuberculosis, systems approaches will be key to define the role of MTBC diversity in the fight against one of humankind's most important pathogens.

The heterogeneous evolution of multidrug-resistant Mycobacterium tuberculosis

Recent surveillance data of multidrug-resistant tuberculosis (MDR-TB) reported the highest rates of resistance ever documented. As further amplification of resistance in MDR strains of Mycobacterium tuberculosis occurs, extensively drug-resistant (XDR) and totally drug-resistant (TDR) TB are beginning to emerge. Although for the most part, the epidemiological factors involved in the spread of MDR-TB are understood, insights into the bacterial drivers of MDR-TB have been gained only recently, largely owing to novel technologies and research in other organisms. Herein, we review recent findings on how bacterial factors, such as persistence, hypermutation, the complex interrelation between drug resistance and fitness, compensatory evolution, and epistasis affect the evolution of multidrug resistance in M. tuberculosis. Improved knowledge of these factors will help better predict the future trajectory of MDR-TB, and contribute to the development of new tools and strategies to combat this growing public health threat.

Mycobacterium tuberculosis Beijing strains favor transmission but not drug resistance in China

BACKGROUND

The Mycobacterium tuberculosis Beijing strains are widespread globally. We aimed to determine whether Beijing strains in China are more likely than other strains to spread, and whether they are more likely to become drug resistant. We also sought to determine whether different Beijing sublineages have distinct phenotypic characteristics.

METHODS

We conducted a population-based molecular epidemiologic study in 6 provinces in China from 2009 to 2010. We analyzed data and specimens from culture-confirmed pulmonary tuberculosis patients. Each patient's isolate was genotyped using 16-loci variable number of tandem repeats and 6 single-nucleotide polymorphisms.

RESULTS

By genotyping, 75.0% (1031/1375) of the strains of M. tuberculosis were Beijing strains. Beijing strains were more likely than non-Beijing strains to be in a genotypic cluster (odds ratio, 2.40, P < .001), and were significantly associated with younger age (P(trend) < .05). There was no significant difference in the proportion of Beijing strains and non-Beijing strains that were drug resistant, even when stratified by new vs retreatment patients. We identified 6 sublineages of Beijing strains in the study population. The modern sublineage of Beijing strains were more likely than the ancient sublineages to be clustered (odds ratio, 2.27, P < .001).

CONCLUSIONS

Beijing strains of M. tuberculosis were significantly associated with genotypic clustering, reflecting recent transmission, and younger age, but were not associated with drug resistance. Future studies of Beijing family strains should avoid assuming and attributing characteristics to the entire family and should assess strains of specific sublineages and/or settings.

Drug susceptibility of Mycobacterium tuberculosis Beijing genotype and association with MDR TB

To determine differences in the ability of Mycobacterium tuberculosis strains to withstand antituberculosis drug treatment, we compared the activity of antituberculosis drugs against susceptible Beijing and East-African/Indian genotype M. tuberculosis strains. Beijing genotype strains showed high rates of mutation within a wide range of drug concentrations, possibly explaining this genotype’s association with multidrug-resistant tuberculosis.

Translating basic science insight into public health action for multidrug- and extensively drug-resistant tuberculosis

Multidrug (MDR)- and extensively drug-resistant (XDR) tuberculosis (TB) impose a heavy toll of human suffering and social costs. Controlling drug-resistant TB is a complex global public health challenge. Basic science advances including elucidation of the genetic basis of resistance have enabled development of new assays that are transforming the diagnosis of MDR-TB. Molecular epidemiological approaches have provided new insights into the natural history of TB with important implications for drug resistance. In the future, progress in understanding Mycobacterium tuberculosis strain-specific human immune responses, integration of systems biology approaches with traditional epidemiology and insight into the biology of mycobacterial persistence have potential to be translated into new tools for diagnosis and treatment of MDR- and XDR-TB. We review recent basic sciences developments that have contributed or may contribute to improved public health response.

A population-based cohort study of Mycobacterium tuberculosis Beijing strains: an emerging public health threat in an immigrant-receiving country?

Introduction

Mycobacterium tuberculosis Beijing strains are frequently associated with tuberculosis outbreaks and drug resistance. However, contradictory evidence and limited study generalizability make it difficult to foresee if the emergence of Beijing strains in high-income immigrant-receiving countries poses an increased public health threat. The purpose of this study was to determine if Beijing strains are associated with high risk disease presentations relative to other strains within Canada.

Methods

This was a retrospective population-based study of culture-confirmed active TB cases in a major immigrant-receiving province of Canada in 1991 through 2007. Of 1,852 eligible cases, 1,826 (99%) were successfully genotyped. Demographic, clinical, and mycobacteriologic surveillance data were combined with molecular diagnostic data. The main outcome measures were site of disease, lung cavitation, sputum smear positivity, bacillary load, and first-line antituberculosis drug resistance.

Results

A total of 350 (19%) patients had Beijing strains; 298 (85%) of these were born in the Western Pacific. Compared to non-Beijing strains, Beijing strains were significantly more likely to be associated with polyresistance (aOR 1.8; 95% CI 1.0–3.3; p = 0.046) and multidrug-resistance (aOR 3.4; 1.0–11.3; p = 0.049). Conversely, Beijing strains were no more likely than non-Beijing strains to be associated with respiratory disease (aOR 1.3; 1.0–1.8; p = 0.053), high bacillary load (aOR 1.2; 0.6–2.7), lung cavitation (aOR 1.0; 0.7–1.5), immediately life-threatening forms of tuberculosis (aOR 0.8; 0.5–1.6), and monoresistance (aOR 0.9; 0.6–1.3). In subgroup analyses, Beijing strains only had a significant association with multidrug-resistant tuberculosis (aOR 6.1; 1.2–30.4), and an association of borderline significance with polyresistant tuberculosis (aOR 1.8; 1.0–3.5; p = 0.062), among individuals born in the Western Pacific.

Conclusion

Other than an increased risk of polyresistant or multidrug-resistant tuberculosis, Beijing strains appear to pose no more of a public health threat than non-Beijing strains within a high-income immigrant-receiving country.

New mutations in the mycobacterial ATP synthase: new insights into the binding of the diarylquinoline TMC207 to the ATP synthase C-ring structure

TMC207 is a new antituberculous drug belonging to the diarylquinoline class which very efficiently inhibits the ATP synthase of mycobacteria such as Mycobacterium tuberculosis, one of the most important pathogens in the world. In order to map the amino acid residues involved in the binding of the drug, we have selected in vitro TMC207-resistant mutants from M. tuberculosis and diverse atypical mycobacteria. Six distinct mutations, Asp28 → Gly, Asp28 → Ala, Leu59 → Val, Glu61 → Asp, Ala63 → Pro, and Ile66 → Met, have been identified in the subunit c forming a C ring in the ATP synthase. They were studied by evaluating the levels of resistance that they confer in the selected clones and by using an isogenic complementation system in Mycobacterium smegmatis. The rates of increase of TMC207 MIC values (8- to 133-fold) were interpreted by constructing by homology modeling a structure of the mycobacterial C ring which was used for docking simulations with TMC207. Our results suggest that the residues found to be mutated in the resistant clones, together with a tyrosine specifically conserved at position 64 in mycobacteria, define a cleft located between two adjacent c subunits in the C ring. This cleft, which encompasses the proton-binding site (Glu61), is well fitted to bind TMC207 at the level of the bromoquinoline moiety, with the drug being anchored by several ionic, hydrogen, and halogen bonds with residues Glu61, Tyr64, and Asp28, respectively. These data shed light on the molecular interactions allowing TMC207 to bind specifically and efficiently at the level of the proton-binding site of the mycobacterial C ring.

Pre-existing isoniazid resistance, but not the genotype of Mycobacterium tuberculosis drives rifampicin resistance codon preference in vitro

Both the probability of a mutation occurring and the ability of the mutant to persist will influence the distribution of mutants that arise in a population. We studied the interaction of these factors for the in vitro selection of rifampicin (RIF)-resistant mutants of Mycobacterium tuberculosis. We characterised two series of spontaneous RIF-resistant in vitro mutants from isoniazid (INH)-sensitive and -resistant laboratory strains and clinical isolates, representing various M. tuberculosis genotypes. The first series were selected from multiple parallel 1 ml cultures and the second from single 10 ml cultures. RIF-resistant mutants were screened by Multiplex Ligation-dependent Probe Amplification (MLPA) or by sequencing the rpoB gene. For all strains the mutation rate for RIF resistance was determined with a fluctuation assay. The most striking observation was a shift towards rpoB-S531L (TCG→TTG) mutations in a panel of laboratory-generated INH-resistant mutants selected from the 10-ml cultures (p<0.001). All tested strains showed similar mutation rates (1.33×10⁻⁸ to 2.49×10⁻⁷) except one of the laboratory-generated INH mutants with a mutation rate measured at 5.71×10⁻⁷, more than 10 times higher than that of the INH susceptible parental strain (5.46–7.44×10⁻⁸). No significant, systematic difference in the spectrum of rpoB-mutations between strains of different genotypes was observed. The dramatic shift towards rpoB-S531L in our INH-resistant laboratory mutants suggests that the relative fitness of resistant mutants can dramatically impact the distribution of (subsequent) mutations that accumulate in a M. tuberculosis population, at least in vitro. We conclude that, against specific genetic backgrounds, certain resistance mutations are particularly likely to spread. Molecular screening for these (combinations of) mutations in clinical isolates could rapidly identify these particular pathogenic strains. We therefore recommend that isolates are screened for the distribution of resistance mutations, especially in regions that are highly endemic for (multi)drug resistant tuberculosis.