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

Clinical implication of novel drug resistance-conferring mutations in resistant tuberculosis

Evolving novel and/or unfamiliar mutations are revolutionizing the pathways of antibiotic resistance of clinical tuberculosis. The accumulation and interaction of these poorly characterized mutations augment the complexity of resistant pathogenic strains and raise public health concerns. This article reviews our current understanding of the genetic changes that characterize drug resistance in tuberculosis and highlights the imperative for further investigations focusing on the effects of an individual mutation and interacting mutations with detailed strain epidemiology, particularly as these pertain to technology-limited countries with high tuberculosis incidence rates. Concomitantly, there is a need for the development, testing, and uptake of new tools for studying the effects of these mutations in drug resistance and fitness cost of the pathogen. Such genetic data are critical for effective localized and global tuberculosis control interventions and for accurate epidemiological predictions.

A review of computational and mathematical modeling contributions to our understanding of Mycobacterium tuberculosis within-host infection and treatment

Tuberculosis (TB) is an ancient and deadly disease characterized by complex host-pathogen dynamics playing out over multiple time and length scales and physiological compartments. Computational modeling can be used to integrate various types of experimental data and suggest new hypotheses, mechanisms, and therapeutic approaches to TB. Here, we offer a first-time comprehensive review of work on within-host TB models that describe the immune response of the host to infection, including the formation of lung granulomas. The models include systems of ordinary and partial differential equations and agent-based models as well as hybrid and multi-scale models that are combinations of these. Many aspects of M. tuberculosis infection, including host dynamics in the lung (typical site of infection for TB), granuloma formation, roles of cytokine and chemokine dynamics, and bacterial nutrient availability have been explored. Finally, we survey applications of these within-host models to TB therapy and prevention and suggest future directions to impact this global disease.

What We Have Learned and What We Have Missed in Tuberculosis Pathophysiology for a New Vaccine Design: Searching for the "Pink Swan"

This is a call to encourage the search for a new vaccine to stop the progression of Mycobacterium tuberculosis infection to tuberculosis (TB) disease. TB is a highly discreet and stigmatized disease, with a massive impact on human health. It has killed 1.2 billion people in the last 200 years and still kills 1.5 million people per year. Over the last 20 years, the TB vaccine field has experienced spectacular developments, and we have learned about (1) the importance of the Th1 response in controlling infection, mainly against RD1 and Ag85 antigens; (2) the stability of the antigenic repertoire; (3) the dynamics of M. tuberculosis granulomas; or (4) the link between typical and atypical pulmonary TB and the immune status of the host. However, we still do not (1) know how to avoid M. tuberculosis infection and reinfection; (2) understand the major role of the increase in lesion size in progression from infection to disease; (3) the role of interlobular septa in encapsulating pulmonary lesions; or (4) the role of neutrophilic infiltration and an exaggerated inflammatory response in the development of TB disease. These are strong reasons to pursue new, imaginative proposals involving both the antibody response and a balanced, tolerant immune response that averts progression toward TB. So far, the scientific mindset has been quite monolithic and has mainly focused on the stimulation of conventional T cells. But this approach has failed. For that reason, we are seeking unconventional perspectives to find a “pink swan,” a more efficacious and safer vaccine candidate.

A Nonsynonymous SNP Catalog of Mycobacterium tuberculosis Virulence Genes and Its Use for Detecting New Potentially Virulent Sublineages

Mycobacterium tuberculosis is divided into several distinct lineages, and various genetic markers such as IS-elements, VNTR, and SNPs are used for lineage identification. We propose an M. tuberculosis classification approach based on functional polymorphisms in virulence genes. An M. tuberculosis virulence genes catalog has been established, including 319 genes from various protein groups, such as proteases, cell wall proteins, fatty acid and lipid metabolism proteins, sigma factors, toxin–antitoxin systems. Another catalog of 1,573 M. tuberculosis isolates of different lineages has been developed. The developed SNP-calling program has identified 3,563 nonsynonymous SNPs. The constructed SNP-based phylogeny reflected the evolutionary relationship between lineages and detected new sublineages. SNP analysis of sublineage F15/LAM4/KZN revealed four lineage-specific mutations in cyp125, mce3B, vapC25, and vapB34. The Ural lineage has been divided into two geographical clusters based on different SNPs in virulence genes. A new sublineage, B0/N-90, was detected inside the Beijing-B0/W-148 by SNPs in irtB, mce3F and vapC46. We have found 27 members of B0/N-90 among the 227 available genomes of the Beijing-B0/W-148 sublineage. Whole-genome sequencing of strain B9741, isolated from an HIV-positive patient, was demonstrated to belong to the new B0/N-90 group. A primer set for PCR detection of B0/N-90 lineage-specific mutations has been developed. The prospective use of mce3 mutant genes as genetically engineered vaccine is discussed.

Mycobacterium tuberculosis Complex Exhibits Lineage-Specific Variations Affecting Protein Ductility and Epitope Recognition

The advent of whole-genome sequencing has provided an unprecedented detail about the evolution and genetic significance of species-specific variations across the whole Mycobacterium tuberculosis Complex. However, little attention has been focused on understanding the functional roles of these variations in the protein coding sequences. In this work, we compare the coding sequences from 74 sequenced mycobacterial species including M. africanum, M. bovis, M. canettii, M. caprae, M. orygis, and M. tuberculosis. Results show that albeit protein variations affect all functional classes, those proteins involved in lipid and intermediary metabolism and respiration have accumulated mutations during evolution. To understand the impact of these mutations on protein functionality, we explored their implications on protein ductility/disorder, a yet unexplored feature of mycobacterial proteomes. In agreement with previous studies, we found that a Gly71Ile substitution in the PhoPR virulence system severely affects the ductility of its nearby region in M. africanum and animal-adapted species. In the same line of evidence, the SmtB transcriptional regulator shows amino acid variations specific to the Beijing lineage, which affects the flexibility of the N-terminal trans-activation domain. Furthermore, despite the fact that MTBC epitopes are evolutionary hyperconserved, we identify strain- and lineage-specific amino acid mutations affecting previously known T-cell epitopes such as EsxH and FbpA (Ag85A). Interestingly, in silico studies reveal that these variations result in differential interaction of epitopes with the main HLA haplogroups.

Positive epistasis of major low-cost drug resistance mutations rpoB531-TTG and katG315-ACC depends on the phylogenetic background of Mycobacterium tuberculosis strains

Mycobacterium tuberculosis Beijing genotype strains increasingly circulate in different world regions, either as historical endemic, e.g. in East Asia, or recently imported, e.g. in South America, and this family is regarded as the most successful lineage of the global tuberculosis (TB) epidemic. Here we analysed the transmission capacity of these strains in the context of their phylogenetic background and drug resistance mutations. The study collection included all multidrug resistant (MDR) strains of Beijing genotype isolated in Beijing Chest Hospital, the largest tertiary TB facility in North China, in 2011-2013 (n = 278). Strains were subjected to NTF/IS6110 and 24-loci MIRU-VNTR analysis. Drug resistance mutations were detected in rpoB, katG, inhA and oxyR-ahpC. A total of 58 and 220 strains were assigned to the ancient and modern Beijing sublineages, respectively. 24-MIRU-VNTR clustering was higher in modern versus ancient Beijing strains (35.9% vs. 12.1%; P <0.001). After taking into consideration the presence of rpoB and katG mutations, clustering decreased to 15.9% in modern and 0% in ancient strains. The most frequent combination of mutations (rpoB531-TTG and katG315-ACC) was more prevalent in clustered versus non-clustered isolates in the modern sublineage (23/35 vs. 47/185; P <0.0001). To conclude, a combination of the known low-fitness-cost rpoB531-TTG and katG315-ACC mutations likely facilitates the increased transmission ability of MDR strains of the modern but not ancient Beijing sublineage. Accordingly, positive epistasis of major low-cost drug resistance-conferring mutations is influenced by the phylogenetic background of M. tuberculosis strains.

The within-host population dynamics of Mycobacterium tuberculosis vary with treatment efficacy

Background

Combination therapy is one of the most effective tools for limiting the emergence of drug resistance in pathogens. Despite the widespread adoption of combination therapy across diseases, drug resistance rates continue to rise, leading to failing treatment regimens. The mechanisms underlying treatment failure are well studied, but the processes governing successful combination therapy are poorly understood. We address this question by studying the population dynamics of Mycobacterium tuberculosis within tuberculosis patients undergoing treatment with different combinations of antibiotics.

Results

By combining very deep whole genome sequencing (~1000-fold genome-wide coverage) with sequential sputum sampling, we were able to detect transient genetic diversity driven by the apparently continuous turnover of minor alleles, which could serve as the source of drug-resistant bacteria. However, we report that treatment efficacy has a clear impact on the population dynamics: sufficient drug pressure bears a clear signature of purifying selection leading to apparent genetic stability. In contrast, M. tuberculosis populations subject to less drug pressure show markedly different dynamics, including cases of acquisition of additional drug resistance.

Conclusions

Our findings show that for a pathogen like M. tuberculosis, which is well adapted to the human host, purifying selection constrains the evolutionary trajectory to resistance in effectively treated individuals. Nonetheless, we also report a continuous turnover of minor variants, which could give rise to the emergence of drug resistance in cases of drug pressure weakening. Monitoring bacterial population dynamics could therefore provide an informative metric for assessing the efficacy of novel drug combinations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1196-0) contains supplementary material, which is available to authorized users.

Acquisition of Rifampin Resistance in Pulmonary Tuberculosis

Mycobacterium tuberculosis strains with spontaneous mutations conferring resistance to rifampin (RIF) are exceedingly rare, and fixed drug combinations typically prevent augmentation of resistance to single drugs. Fourteen newly diagnosed tuberculosis patients were treated with RIF alone for 14 days, and bacterial loads, including mutation frequencies, were determined. A statistical model estimated that 1% of the remaining viable mycobacteria could be RIF resistant after 30 days of monotherapy. This indicates that temporal and spatial windows of RIF monotherapy due to uneven drug distribution within lung lesions could contribute to the acquisition of resistance to RIF.

The epidemiology, pathogenesis, transmission, diagnosis, and management of multidrug-resistant, extensively drug-resistant, and incurable tuberculosis

Global tuberculosis incidence has declined marginally over the past decade, and tuberculosis remains out of control in several parts of the world including Africa and Asia. Although tuberculosis control has been effective in some regions of the world, these gains are threatened by the increasing burden of multidrug-resistant (MDR) and extensively drug-resistant (XDR) tuberculosis. XDR tuberculosis has evolved in several tuberculosis-endemic countries to drug-incurable or programmatically incurable tuberculosis (totally drug-resistant tuberculosis). This poses several challenges similar to those encountered in the pre-chemotherapy era, including the inability to cure tuberculosis, high mortality, and the need for alternative methods to prevent disease transmission. This phenomenon mirrors the worldwide increase in antimicrobial resistance and the emergence of other MDR pathogens, such as malaria, HIV, and Gram-negative bacteria. MDR and XDR tuberculosis are associated with high morbidity and substantial mortality, are a threat to health-care workers, prohibitively expensive to treat, and are therefore a serious public health problem. In this Commission, we examine several aspects of drug-resistant tuberculosis. The traditional view that acquired resistance to antituberculous drugs is driven by poor compliance and programmatic failure is now being questioned, and several lines of evidence suggest that alternative mechanisms-including pharmacokinetic variability, induction of efflux pumps that transport the drug out of cells, and suboptimal drug penetration into tuberculosis lesions-are likely crucial to the pathogenesis of drug-resistant tuberculosis. These factors have implications for the design of new interventions, drug delivery and dosing mechanisms, and public health policy. We discuss epidemiology and transmission dynamics, including new insights into the fundamental biology of transmission, and we review the utility of newer diagnostic tools, including molecular tests and next-generation whole-genome sequencing, and their potential for clinical effectiveness. Relevant research priorities are highlighted, including optimal medical and surgical management, the role of newer and repurposed drugs (including bedaquiline, delamanid, and linezolid), pharmacokinetic and pharmacodynamic considerations, preventive strategies (such as prophylaxis in MDR and XDR contacts), palliative and patient-orientated care aspects, and medicolegal and ethical issues.

Quadruple-first line drug resistance in Mycobacterium tuberculosis in Vietnam: What can we learn from genes?

In Vietnam, a country with high tuberculosis (137/100.000 population) and multidrug-resistant (MDR)-TB burdens (7.8/100.000 population), little is known about the molecular signatures of drug resistance in general and more particularly of second line drug (SLD) resistance. This study is specifically focused on Mycobacterium tuberculosis isolates resistant to four first-line drugs (FLDs) that make TB much more difficult to treat. The aim is to determine the proportion of SLD resistance in these quadruple drug resistant isolates and the genetic determinants linked to drug resistance to better understand the genetic processes leading to quadruple and extremely drug resistance (XDR). 91 quadruple (rifampicin, isoniazid, ethambutol and streptomycin) FLD resistant and 55 susceptible isolates were included. Spoligotyping and 24-locus MIRU-VNTR techniques were performed and 9 genes and promoters linked to FLD and SLD resistance were sequenced. SLD susceptibility testing was carried out on a subsample of isolates. High proportion of quadruple-FLD resistant isolates was resistant to fluoroquinolones (27%) and second-line injectable drugs (30.2%) by drug susceptibility testing. The sequencing revealed high mutation diversity with prevailing mutations at positions katG315, inhA-15, rpoB531, embB306, rrs1401, rpsL43 and gyrA94. The sensitivity and specificity were high for most drug resistances (>86%), but the sensitivity was lower for injectable drug resistances (<69%). The mutation patterns revealed 23.1% of pre-XDR and 7.7% of XDR isolates, mostly belonging to Beijing family. The genotypic diversity and the variety of mutations reflect the existence of various evolutionary paths leading to FLD and SLD resistance. Nevertheless, particular mutation patterns linked to high-level resistance and low fitness costs seem to be favored.

A non-canonical mismatch repair pathway in prokaryotes

Mismatch repair (MMR) is a near ubiquitous pathway, essential for the maintenance of genome stability. Members of the MutS and MutL protein families perform key steps in mismatch correction. Despite the major importance of this repair pathway, MutS-MutL are absent in almost all Actinobacteria and many Archaea. However, these organisms exhibit rates and spectra of spontaneous mutations similar to MMR-bearing species, suggesting the existence of an alternative to the canonical MutS-MutL-based MMR. Here we report that Mycobacterium smegmatis NucS/EndoMS, a putative endonuclease with no structural homology to known MMR factors, is required for mutation avoidance and anti-recombination, hallmarks of the canonical MMR. Furthermore, phenotypic analysis of naturally occurring polymorphic NucS in a M. smegmatis surrogate model, suggests the existence of M. tuberculosis mutator strains. The phylogenetic analysis of NucS indicates a complex evolutionary process leading to a disperse distribution pattern in prokaryotes. Together, these findings indicate that distinct pathways for MMR have evolved at least twice in nature.

Genomic analysis of globally diverse Mycobacterium tuberculosis strains provides insights into the emergence and spread of multidrug resistance

Multidrug-resistant tuberculosis (MDR-TB), caused by drug resistant strains of Mycobacterium tuberculosis, is an increasingly serious problem worldwide. In this study, we examined a dataset of 5,310 M. tuberculosis whole genome sequences from five continents. Despite great diversity with respect to geographic point of isolation, genetic background and drug resistance, patterns of drug resistance emergence were conserved globally. We have identified harbinger mutations that often precede MDR. In particular, the katG S315T mutation, conferring resistance to isoniazid, overwhelmingly arose before rifampicin resistance across all lineages, geographic regions, and time periods. Molecular diagnostics that include markers for rifampicin resistance alone will be insufficient to identify pre-MDR strains. Incorporating knowledge of pre-MDR polymorphisms, particularly katG S315, into molecular diagnostics will enable targeted treatment of patients with pre-MDR-TB to prevent further development of MDR-TB.

Metabolic Perspectives on Persistence

Accumulating evidence has left little doubt about the importance of persistence or metabolism in the biology and chemotherapy of tuberculosis. However, knowledge of the intersection between these two factors has only recently begun to emerge. Here, we provide a focused review of metabolic characteristics associated with Mycobacterium tuberculosis persistence. We focus on metabolism because it is the biochemical foundation of all physiologic processes and a distinguishing hallmark of M. tuberculosis physiology and pathogenicity. In addition, it serves as the chemical interface between host and pathogen. Existing knowledge, however, derives largely from physiologic contexts in which replication is the primary biochemical objective. The goal of this review is to reframe current knowledge of M. tuberculosis metabolism in the context of persistence, where quiescence is often a key distinguishing characteristic. Such a perspective may help ongoing efforts to develop more efficient cures and inform on novel strategies to break the cycle of transmission sustaining the pandemic.

Strain Diversity and the Evolution of Antibiotic Resistance

Drug resistance is best thought of as an ongoing biological process. Resistant bacteria must emerge, become established and ultimately transmit in order to be relevant to human health. In this context, genetic diversity can influence the rate and likelihood of resistance emerging; it can also modulate the net physiological impact of resistance and the propensity of an organism to improve any defects that arise from it. Combined, these effects can have an impact on a larger scale, with highly transmissible drug-resistant bacterial strains posing a formidable threat to global health. These considerations are pertinent to the future of tuberculosis control as well. In this chapter, we review our current understanding of the impact of genetic diversity in the broadest sense on the evolution of drug-resistant members of the Mycobacterium tuberculosis complex.

Unique Regulation of the DosR Regulon in the Beijing Lineage of Mycobacterium tuberculosis

The DosR regulon, a set of 48 genes normally expressed in Mycobacterium tuberculosis under conditions that inhibit aerobic respiration, is controlled via the DosR-DosS/DosT two-component system. While the regulon requires induction in most M. tuberculosis isolates, for members of the Beijing lineage, its expression is uncoupled from the need for signaling. In our attempts to understand the mechanistic basis for this uncoupling in the Beijing background, we previously reported the identification of two synonymous single-nucleotide polymorphisms (SNPs) within the adjacent Rv3134c gene. In the present study, we have interrogated the impact of these SNPs on dosR expression in wild-type strains, as well as a range of dosR-dosS-dosT mutants, for both Beijing and non-Beijing M. tuberculosis backgrounds. In this manner, we have unequivocally determined that the C601T dosR promoter SNP is the sole requirement for the dramatic shift in the pattern of DosR regulon expression seen in this globally important lineage. Interestingly, we also show that DosT is completely nonfunctional within these strains. Thus, a complex series of evolutionary steps has led to the present-day Beijing DosR phenotype that, in turn, potentially confers a fitness advantage in the face of some form of host-associated selective pressure.

IMPORTANCE

Mycobacterium tuberculosis strains of the Beijing lineage have been described as being of enhanced virulence compared to other lineages, and in certain regions, they are associated with the dramatic spread of multidrug-resistant tuberculosis (TB). In terms of trying to understand the functional basis for these broad epidemiological phenomena, it is interesting that, in contrast to the other major lineages, the Beijing strains all constitutively overexpress members of the DosR regulon. Here, we identify the mutational events that led to the evolution of this unique phenotype. In addition, our work highlights the fact that important phenotypic differences exist between distinct M. tuberculosis lineages, with the potential to impact the efficacy of diagnosis, vaccination, and treatment programs.

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.

Genomic epidemiology of Lineage 4 Mycobacterium tuberculosis subpopulations in New York City and New Jersey, 1999-2009

Background

Whole genome sequencing (WGS) has rapidly become an important research tool in tuberculosis epidemiology and is likely to replace many existing methods in public health microbiology in the near future. WGS-based methods may be particularly useful in areas with less diverse Mycobacterium tuberculosis populations, such as New York City, where conventional genotyping is often uninformative and field epidemiology often difficult. This study applies four candidate strategies for WGS-based identification of emerging M. tuberculosis subpopulations, employing both phylogenomic and population genetics methods.

Results

M. tuberculosis subpopulations in New York City and New Jersey can be distinguished via phylogenomic reconstruction, evidence of demographic expansion and subpopulation-specific signatures of selection, and by determination of subgroup-defining nucleotide substitutions. These methods identified known historical outbreak clusters and previously unidentified subpopulations within relatively monomorphic M. tuberculosis endemic clone groups. Neutrality statistics based on the site frequency spectrum were less useful for identifying M. tuberculosis subpopulations, likely due to the low levels of informative genetic variation in recently diverged isolate groups. In addition, we observed that isolates from New York City endemic clone groups have acquired multiple non-synonymous SNPs in virulence- and growth-associated pathways, and relatively few mutations in drug resistance-associated genes, suggesting that overall pathoadaptive fitness, rather than the acquisition of drug resistance mutations, has played a central role in the evolutionary history and epidemiology of M. tuberculosis subpopulations in New York City.

Conclusions

Our results demonstrate that some but not all WGS-based methods are useful for detection of emerging M. tuberculosis clone groups, and support the use of phylogenomic reconstruction in routine tuberculosis laboratory surveillance, particularly in areas with relatively less diverse M. tuberculosis populations. Our study also supports the use of wider-reaching phylogenomic and population genomic methods in tuberculosis public health practice, which can support tuberculosis control activities by identifying genetic polymorphisms contributing to epidemiological success in local M. tuberculosis populations and possibly explain why certain isolate groups are apparently more successful in specific host populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3298-6) contains supplementary material, which is available to authorized users.

Mutation tendency of mutator Plasmodium berghei with proofreading-deficient DNA polymerase δ

In this study, we investigated the mutation tendency of a mutator rodent malaria parasite, Plasmodium berghei, with proofreading-deficient DNA polymerase δ. Wild-type and mutator parasites were maintained in mice for over 24 weeks, and the genome-wide accumulated mutations were determined by high-throughput sequencing. The mutator P. berghei had a significant preference for C/G to A/T substitutions; thus, its genome had a trend towards a higher AT content. The mutation rate was influenced by the sequence context, and mutations were markedly elevated at TCT. Some genes mutated repeatedly in replicate passage lines. In particular, knockout mutations of the AP2-G gene were frequent, which conferred strong growth advantages on parasites during the blood stage but at the cost of losing the ability to form gametocytes. This is the first report to demonstrate a biased mutation tendency in malaria parasites, and its results help to promote our basic understanding of Plasmodium genetics.

Genomic diversity in autopsy samples reveals within-host dissemination of HIV-associated Mycobacterium tuberculosis

Mycobacterium tuberculosis remains a leading cause of death worldwide, especially among individuals infected with HIV. Whereas phylogenetic analysis has revealed M. tuberculosis spread throughout history and in local outbreaks, much less is understood about its dissemination within the body. Here we report genomic analysis of 2,693 samples collected post mortem from lung and extrapulmonary biopsies of 44 subjects in KwaZulu-Natal, South Africa, who received minimal antitubercular treatment and most of whom were HIV seropositive. We found that purifying selection occurred within individual patients, without the need for patient-to-patient transmission. Despite negative selection, mycobacteria diversified within individuals to form sublineages that co-existed for years. These sublineages, as well as distinct strains from mixed infections, were differentially distributed throughout the lung, suggesting temporary barriers to pathogen migration. As a consequence, samples taken from the upper airway often captured only a fraction of the population diversity, challenging current methods of outbreak tracing and resistance diagnostics. Phylogenetic analysis indicated that dissemination from the lungs to extrapulmonary sites was as frequent as between lung sites, supporting the idea of similar migration routes within and between organs, at least in subjects with HIV. Genomic diversity therefore provides a record of pathogen diversification and repeated dissemination across the body.

Whole genome sequencing for M/XDR tuberculosis surveillance and for resistance testing

Whole genome sequencing (WGS) can help to relate Mycobacterium tuberculosis genomes to one another to assess genetic relatedness and infer the likelihood of transmission between cases. The same sequence data are now increasingly being used to predict drug resistance and susceptibility. Controlling the spread of tuberculosis and providing patients with the correct treatment are central to the World Health Organization's target to 'End TB' by 2035, for which the global prevalence of drug-resistant tuberculosis remains one of the main obstacles to success. So far, WGS has been applied largely to drug-susceptible strains for the purposes of understanding transmission, leaving a number of analytical considerations before transferring what has been learnt from drug-susceptible disease to drug-resistant tuberculosis. We discuss these potential problems here, alongside some of the challenges to characterizing the Mycobacterium tuberculosis 'resistome'-the optimal knowledge-base required for WGS-based assays to successfully direct individualized treatment regimens through the prediction of drug resistance and susceptibility in the future.

Turkish and Japanese Mycobacterium tuberculosis sublineages share a remote common ancestor

Two geographically distant M. tuberculosis sublineages, Tur from Turkey and T3-Osaka from Japan, exhibit partially identical genotypic signatures (identical 12-loci MIRU-VNTR profiles, distinct spoligotyping patterns). We investigated T3-Osaka and Tur sublineages characteristics and potential genetic relatedness, first using MIRU-VNTR locus analysis on 21 and 25 samples of each sublineage respectively, and second comparing Whole Genome Sequences of 8 new samples to public data from 45 samples uncovering human tuberculosis diversity. We then tried to date their Most Recent Common Ancestor (MRCA) using three calibrations of SNP accumulation rate (long-term=0.03SNP/genome/year, derived from a tuberculosis ancestor of around 70,000years old; intermediate=0.2SNP/genome/year derived from a Peruvian mummy; short-term=0.5SNP/genome/year). To disentangle between these scenarios, we confronted the corresponding divergence times with major human history events and knowledge on human genetic divergence. We identified relatively high intrasublineage diversity for both T3-Osaka and Tur. We definitively proved their monophyly; the corresponding super-sublineage (referred to as "T3-Osa-Tur") shares a common ancestor with T3-Ethiopia and Ural sublineages but is only remotely related to other Euro-American sublineages such as X, LAM, Haarlem and S. The evolutionary scenario based on long-term evolution rate being valid until T3-Osa-Tur MRCA was not supported by Japanese fossil data. The evolutionary scenario relying on short-term evolution rate since T3-Osa-Tur MRCA was contradicted by human history and potential traces of past epidemics. T3-Osaka and Tur sublineages were found likely to have diverged between 800y and 2000years ago, potentially at the time of Mongol Empire. Altogether, this study definitively proves a strong genetic link between Turkish and Japanese tuberculosis. It provides a first hypothesis for calibrating TB Euro-American lineage molecular clock; additional studies are needed to reliably date events corresponding to intermediate depths in tuberculosis phylogeny.

Whole Genome Sequencing Demonstrates Limited Transmission within Identified Mycobacterium tuberculosis Clusters in New South Wales, Australia

Australia has a low tuberculosis incidence rate with most cases occurring among recent immigrants. Given suboptimal cluster resolution achieved with 24-locus mycobacterium interspersed repetitive unit (MIRU-24) genotyping, the added value of whole genome sequencing was explored. MIRU-24 profiles of all Mycobacterium tuberculosis culture-confirmed tuberculosis cases diagnosed between 2009 and 2013 in New South Wales (NSW), Australia, were examined and clusters identified. The relatedness of cases within the largest MIRU-24 clusters was assessed using whole genome sequencing and phylogenetic analyses. Of 1841 culture-confirmed TB cases, 91.9% (1692/1841) had complete demographic and genotyping data. East-African Indian (474; 28.0%) and Beijing (470; 27.8%) lineage strains predominated. The overall rate of MIRU-24 clustering was 20.1% (340/1692) and was highest among Beijing lineage strains (35.7%; 168/470). One Beijing and three East-African Indian (EAI) clonal complexes were responsible for the majority of observed clusters. Whole genome sequencing of the 4 largest clusters (30 isolates) demonstrated diverse single nucleotide polymorphisms (SNPs) within identified clusters. All sequenced EAI strains and 70% of Beijing lineage strains clustered by MIRU-24 typing demonstrated distinct SNP profiles. The superior resolution provided by whole genome sequencing demonstrated limited M. tuberculosis transmission within NSW, even within identified MIRU-24 clusters. Routine whole genome sequencing could provide valuable public health guidance in low burden settings.

Drug resistance mechanisms and novel drug targets for tuberculosis therapy

Drug-resistant tuberculosis (TB) poses a significant challenge to the successful treatment and control of TB worldwide. Resistance to anti-TB drugs has existed since the beginning of the chemotherapy era. New insights into the resistant mechanisms of anti-TB drugs have been provided. Better understanding of drug resistance mechanisms helps in the development of new tools for the rapid diagnosis of drug-resistant TB. There is also a pressing need in the development of new drugs with novel targets to improve the current treatment of TB and to prevent the emergence of drug resistance in Mycobacterium tuberculosis. This review summarizes the anti-TB drug resistance mechanisms, furnishes some possible novel drug targets in the development of new agents for TB therapy and discusses the usefulness using known targets to develop new anti-TB drugs. Whole genome sequencing is currently an advanced technology to uncover drug resistance mechanisms in M. tuberculosis. However, further research is required to unravel the significance of some newly discovered gene mutations in their contribution to drug resistance.

Multiple major disease-associated clones of Legionella pneumophila have emerged recently and independently

Legionella pneumophila is an environmental bacterium and the leading cause of Legionnaires' disease. Just five sequence types (ST), from more than 2000 currently described, cause nearly half of disease cases in northwest Europe. Here, we report the sequence and analyses of 364 L. pneumophila genomes, including 337 from the five disease-associated STs and 27 representative of the species diversity. Phylogenetic analyses revealed that the five STs have independent origins within a highly diverse species. The number of de novo mutations is extremely low with maximum pairwise single-nucleotide polymorphisms (SNPs) ranging from 19 (ST47) to 127 (ST1), which suggests emergences within the last century. Isolates sampled geographically far apart differ by only a few SNPs, demonstrating rapid dissemination. These five STs have been recombining recently, leading to a shared pool of allelic variants potentially contributing to their increased disease propensity. The oldest clone, ST1, has spread globally; between 1940 and 2000, four new clones have emerged in Europe, which show long-distance, rapid dispersal. That a large proportion of clinical cases is caused by recently emerged and internationally dispersed clones, linked by convergent evolution, is surprising for an environmental bacterium traditionally considered to be an opportunistic pathogen. To simultaneously explain recent emergence, rapid spread and increased disease association, we hypothesize that these STs have adapted to new man-made environmental niches, which may be linked by human infection and transmission.

Sequence element enrichment analysis to determine the genetic basis of bacterial phenotypes

Bacterial genomes vary extensively in terms of both gene content and gene sequence. This plasticity hampers the use of traditional SNP-based methods for identifying all genetic associations with phenotypic variation. Here we introduce a computationally scalable and widely applicable statistical method (SEER) for the identification of sequence elements that are significantly enriched in a phenotype of interest. SEER is applicable to tens of thousands of genomes by counting variable-length k-mers using a distributed string-mining algorithm. Robust options are provided for association analysis that also correct for the clonal population structure of bacteria. Using large collections of genomes of the major human pathogens Streptococcus pneumoniae and Streptococcus pyogenes, SEER identifies relevant previously characterized resistance determinants for several antibiotics and discovers potential novel factors related to the invasiveness of S. pyogenes. We thus demonstrate that our method can answer important biologically and medically relevant questions.

Plasticity and clonal population structure in bacterial genomes can hinder traditional SNP-based genetic association studies. Here, Corander and colleagues present a method to identify variable-length sequence elements enriched in a phenotype of interest, and demonstrate its use in human pathogens.

Parallel reaction monitoring of clinical Mycobacterium tuberculosis lineages reveals pre-existent markers of rifampicin tolerance in the emerging Beijing lineage

The spread of multidrug resistant Mycobacterium tuberculosis is one of the major challenges in tuberculosis control. In Eurasia, the spread of multidrug resistant tuberculosis is driven by the M. tuberculosis Beijing genotype. In this study, we examined whether selective advantages are present in the proteome of Beijing isolates that contribute to the emergence of this genotype. To this end, we compared the proteome of M. tuberculosis Beijing to that of M. tuberculosis H37Rv, both in the presence and absence of the first-line antibiotic rifampicin. During rifampicin exposure, both M. tuberculosis genotypes express proteins belonging to the DosR dormancy regulon, which induces a metabolically hypoactive-, drug tolerant phenotype. However, these markers of rifampicin tolerance were already more abundant in the M. tuberculosis Beijing isolate prior to drug exposure. To determine whether the a priori high abundance of specific proteins contribute to the formation of antibiotic resistance in M. tuberculosis Beijing, we quantified the abundance of 33 selected proteins in 27 clinical isolates from the five most common M. tuberculosis lineages using parallel reaction monitoring. The observed pre-existing high abundance of dormancy proteins in Beijing strains provides an evolutionary advantage that allows these strains to persist for prolonged periods during rifampicin treatment.

SIGNIFICANCE

M. tuberculosis is the leading cause of death by a bacterial infection worldwide. Treatment-regimen to eradicate this pathogen make use of the first-line antibiotic rifampicin, which is considered to be the cornerstone of modern day anti-tuberculosis treatment. Despite the potency of rifampicin, there is an increasing occurrence of rifampicin resistant mutants in a specific cluster of M. tuberculosis, the Beijing genotype. Using both a data dependent acquisition and a targeted proteomic approach we identified markers of rifampicin tolerance to be high abundant in members of the M. tuberculosis Beijing genotype, already prior drug exposure. The identification of this M. tuberculosis Beijing specific trait will contribute to improved diagnostics and treatment of M. tuberculosis.

Drug-resistant tuberculosis viewed from bacterial and host genomes

The outcome of infection with Mycobacterium tuberculosis (MTB) is largely influenced by the host-pathogen interaction in which both the human host and the MTB genetic backgrounds play an important role. Whether this interaction also influences the selection and expansion of drug-resistant MTB strains is the primary focus of this review. We first outline the main and recent findings regarding MTB determinants implicated in the development of drug resistance. Second, we examine data regarding human genetic factors that may play a role in TB drug resistance. We highlight interesting openings for TB research and therapy.

Whole-Genome Sequencing Analysis of Serially Isolated Multi-Drug and Extensively Drug Resistant Mycobacterium tuberculosis from Thai Patients

Multi-drug and extensively drug-resistant tuberculosis (MDR and XDR-TB) are problems that threaten public health worldwide. Only some genetic markers associated with drug-resistant TB are known. Whole-genome sequencing (WGS) is a promising tool for distinguishing between re-infection and persistent infection in isolates taken at different times from a single patient, but has not yet been applied in MDR and XDR-TB. We aim to detect genetic markers associated with drug resistance and distinguish between reinfection and persistent infection from MDR and XDR-TB patients based on WGS analysis. Samples of Mycobacterium tuberculosis (n = 7), serially isolated from 2 MDR cases and 1 XDR-TB case, were retrieved from Siriraj Hospital, Bangkok. The WGS analysis used an Illumina Miseq sequencer. In cases of persistent infection, MDR-TB isolates differed at an average of 2 SNPs across the span of 2–9 months whereas in the case of reinfection, isolates differed at 61 SNPs across 2 years. Known genetic markers associated with resistance were detected from strains susceptible to streptomycin (2/7 isolates), p-aminosalicylic acid (3/7 isolates) and fluoroquinolone drugs. Among fluoroquinolone drugs, ofloxacin had the highest phenotype-genotype concordance (6/7 isolates), whereas gatifloxcain had the lowest (3/7 isolates). A putative candidate SNP in Rv2477c associated with kanamycin and amikacin resistance was suggested for further validation. WGS provided comprehensive results regarding molecular epidemiology, distinguishing between persistent infection and reinfection in M/XDR-TB and potentially can be used for detection of novel mutations associated with drug resistance.

The Evolutionary History, Demography, and Spread of the Mycobacterium tuberculosis Complex

With the advent of next-generation sequencing technology, the genotyping of clinical Mycobacterium tuberculosis strains went through a major breakup that dramatically improved the field of molecular epidemiology but also revolutionized our deep understanding of the M. tuberculosis complex evolutionary history. The intricate paths of the pathogen and its human host are reflected by a common geographical origin in Africa and strong biogeographical associations that largely reflect the past migration waves out of Africa. This long coevolutionary history is cardinal for our understanding of the host-pathogen dynamic, including past and ongoing demographic components, strains’ genetic background, as well as the immune system genetic architecture of the host. Coalescent- and Bayesian-based analyses allowed us to reconstruct population size changes of M. tuberculosis through time, to date the most recent common ancestor and the several phylogenetic lineages. This information will ultimately help us to understand the spread of the Beijing lineage, the rise of multidrug-resistant sublineages, or the fall of others in the light of socioeconomic events, antibiotic programs, or host population densities. If we leave the present and go through the looking glass, thanks to our ability to handle small degraded molecules combined with targeted capture, paleomicrobiology covering the Pleistocene era will possibly unravel lineage replacements, dig out extinct ones, and eventually ask for major revisions of the current model.

The Rate and Spectrum of Spontaneous Mutations in Mycobacterium smegmatis, a Bacterium Naturally Devoid of the Postreplicative Mismatch Repair Pathway

Mycobacterium smegmatis is a bacterium that is naturally devoid of known postreplicative DNA mismatch repair (MMR) homologs, mutS and mutL, providing an opportunity to investigate how the mutation rate and spectrum has evolved in the absence of a highly conserved primary repair pathway. Mutation accumulation experiments of M. smegmatis yielded a base-substitution mutation rate of 5.27 × 10⁻¹⁰ per site per generation, or 0.0036 per genome per generation, which is surprisingly similar to the mutation rate in MMR-functional unicellular organisms. Transitions were found more frequently than transversions, with the A:T→G:C transition rate significantly higher than the G:C→A:T transition rate, opposite to what is observed in most studied bacteria. We also found that the transition-mutation rate of M. smegmatis is significantly lower than that of other naturally MMR-devoid or MMR-knockout organisms. Two possible candidates that could be responsible for maintaining high DNA fidelity in this MMR-deficient organism are the ancestral-like DNA polymerase DnaE1, which contains a highly efficient DNA proofreading histidinol phosphatase (PHP) domain, and/or the existence of a uracil-DNA glycosylase B (UdgB) homolog that might protect the GC-rich M. smegmatis genome against DNA damage arising from oxidation or deamination. Our results suggest that M. smegmatis has a noncanonical Dam (DNA adenine methylase) methylation system, with target motifs differing from those previously reported. The mutation features of M. smegmatis provide further evidence that genomes harbor alternative routes for improving replication fidelity, even in the absence of major repair pathways.

Substantial molecular evolution and mutation rates in prolonged latent Mycobacterium tuberculosis infection in humans

The genome of Mycobacterium tuberculosis (Mtb) of latently infected individuals may hold the key to understanding the processes that lead to reactivation and progression to clinical disease. We report here analysis of pairs of Mtb isolates from putative prolonged latent TB cases. We identified two confirmed cases, and used whole genome sequencing to investigate the mutational processes that occur over decades in latent Mtb. We found an estimated mutation rate between 0.2 and 0.3 over 33 years, suggesting that latent Mtb accumulates mutations at rates similar to observations from cases of active disease.

Host Immune Responses Differ between M. africanum- and M. tuberculosis-Infected Patients following Standard Anti-tuberculosis Treatment

Epidemiological differences exist between Mycobacterium africanum (Maf)- and Mycobacterium tuberculosis (Mtb)-infected patients, but to date, contributing host factors have not been characterised. We analysed clinical outcomes, as well as soluble markers and gene expression profiles in unstimulated, and ESAT6/CFP-10-, whole-Maf- and Mtb-stimulated blood samples of 26 Maf- and 49 Mtb-HIV-negative tuberculosis patients before, and after 2 and 6 months of anti-tuberculosis therapy. Before treatment, both groups had similar clinical parameters, but differed in few cytokines concentration and gene expression profiles. Following treatment the body mass index, skinfold thickness and chest X-ray scores showed greater improvement in the Mtb- compared to Maf-infected patients, after adjusting for age, sex and ethnicity (p = 0.02; 0.04 and 0.007, respectively). In addition, in unstimulated blood, IL-12p70, IL12A and TLR9 were significantly higher in Maf-infected patients, while IL-15, IL-8 and MIP-1α were higher in Mtb-infected patients. Overnight stimulation with ESAT-6/CFP-10 induced significantly higher levels of IFN-γ and TNF-α production, as well as gene expression of CCL4, IL1B and TLR4 in Mtb- compared to Maf-infected patients. Our study confirms differences in clinical features and immune genes expression and concentration of proteins associated with inflammatory processes between Mtb- and Maf-infected patients following anti-tuberculosis treatment These findings have public health implications for treatment regimens, and biomarkers for tuberculosis diagnosis and susceptibility.

Rifampin Resistance Mutations Are Associated with Broad Chemical Remodeling of Mycobacterium tuberculosis

Global control of tuberculosis has become increasingly complicated with the emergence of multidrug-resistant strains of Mycobacterium tuberculosis First-line treatments are anchored by two antibiotics, rifampin and isoniazid. Most rifampin resistance occurs through the acquisition of missense mutations in the rifampin resistance-determining region, an 81-base pair region encoding the rifampin binding site on the β subunit of RNA polymerase (rpoB). Although these mutations confer a survival advantage in the presence of rifampin, they may alter the normal process of transcription, thereby imposing significant fitness costs. Because the downstream biochemical consequences of the rpoB mutations are unknown, we used an organism-wide screen to identify the number and types of lipids changed after rpoB mutation. A new mass spectrometry-based profiling platform systematically compared ∼10,000 cell wall lipids in a panel of rifampin-resistant mutants within two genetically distinct strains, CDC1551and W-Beijing. This unbiased lipidomic survey detected quantitative alterations (>2-fold, p < 0.05) in more than 100 lipids in each mutant. By focusing on molecular events that change among most mutants and in both genetic backgrounds, we found that rifampin resistance mutations lead to altered concentrations of mycobactin siderophores and acylated sulfoglycolipids. These findings validate a new organism-wide lipidomic analysis platform for drug-resistant mycobacteria and provide direct evidence for characteristic remodeling of cell wall lipids in rifampin-resistant strains of M. tuberculosis The specific links between rifampin resistance and named lipid factors provide diagnostic and therapeutic targets that may be exploited to address the problem of drug resistance.

MDR-TB Outbreak among HIV-Negative Tunisian Patients followed during 11 Years

Background

Multidrug-resistant tuberculosis (MDR-TB) outbreaks that evolve, from the outset, in a context strictly negative for HIV infection deserve special consideration since they reflect the true intrinsic epidemic potential of the causative strain. To our knowledge, the long-term evolution of such exceptional outbreaks and the treatment outcomes for the involved patients has never been reported hitherto. Here we provide a thorough description, over an 11-year period, of an MDR-TB outbreak that emerged and expanded in an HIV-negative context, Northern Tunisia.

Methodology/Principal Findings

From October 2001 to June 2011, the MDR-TB outbreak involved 48 HIV-negative individuals that are mainly young (mean age 31.09 yrs; 89.6% male) and noninstitutionalized. Drug susceptibility testing coupled to mutational analysis revealed that initial transmission involved an isolate that was simultaneously resistant to isoniazid, rifampicin, ethambutol, and streptomycin. The causative Haarlem3-ST50 outbreak strain expanded mainly as an 11-banded IS6110 RFLP profile (77.1%), from which a 12-banded subclone evolved. After undergoing a 2-year treatment with second-line drugs, 22 (45.8%) patients were cured and 3 (6.2%) completed treatment, thus yielding an overall treatment success rate of 52.1%. Among the patients that experienced unfavorable treatment outcomes, 10 (20.8%) failed treatment, 3 (6.2%) were lost to follow-up, 5 (10.4%) died, and 5 (10.4%) could not be evaluated. Poor adherence to treatment was found to be the main independent predictor of unfavorable outcomes (HR: 9.15; 95% CI 1.72–48.73; P = 0.014). Intriguingly, the evolved 12-banded subclone proved significantly associated with unfavorable outcomes (HR: 4.90; 95% CI 1.04–23.04, P = 0.044). High rate of fatality and relapse was further demonstrated at the long-term, since 70% of those whose treatment failed have died, and 24% among those deemed successfully treated have relapsed.

Conclusions/Significance

Taken together, the data obtained in this study indicate that MDR-TB clinical isolates could become fit enough to cause large and severe outbreaks in an HIV-negative context. Such MDR-TB outbreaks are characterized by low treatment success rates and could evolve towards increased severity, thus calling for early detection of cases and the necessity to raise the bar of surveillance throughout and beyond the treatment period.

Genetic signatures of Mycobacterium tuberculosis Nonthaburi genotype revealed by whole genome analysis of isolates from tuberculous meningitis patients in Thailand

Genome sequencing plays a key role in understanding the genetic diversity of Mycobacterium tuberculosis (M.tb). The genotype-specific character of M. tb contributes to tuberculosis severity and emergence of drug resistance. Strains of M. tb complex can be classified into seven lineages. The Nonthaburi (NB) genotype, belonging to the Indo-Oceanic lineage (lineage 1), has a unique spoligotype and IS6110-RFLP pattern but has not previously undergone a detailed whole genome analysis. In addition, there is not much information available on the whole genome analysis of M. tb isolates from tuberculous meningitis (TBM) patients in public databases. Isolates CSF3053, 46-5069 and 43-13838 of NB genotype were obtained from the cerebrospinal fluids of TBM Thai patients in Siriraj Hospital, Bangkok. The whole genomes were subjected to high throughput sequencing. The sequence data of each isolate were assembled into draft genome. The sequences were also aligned to reference genome, to determine genomic variations. Single nucleotide polymorphisms (SNPs) were obtained and grouped according to the functions of the genes containing them. They were compared with SNPs from 1,601 genomes, representing the seven lineages of M. tb complex, to determine the uniqueness of NB genotype. Susceptibility to first-line, second-line and other antituberculosis drugs were determined and related to the SNPs previously reported in drug-resistant related genes. The assembled genomes have an average size of 4,364,461 bp, 4,154 genes, 48 RNAs and 64 pseudogenes. A 500 base pairs deletion, which includes ppe50, was found in all isolates. RD239, specific for members of Indo Oceanic lineage, and RD147c were identified. A total of 2,202 SNPs were common to the isolates and used to classify the NB strains as members of sublineage 1.2.1. Compared with 1,601 genomes from the seven lineages of M. tb complex, mutation G2342203C was found novel to the isolates in this study. Three mutations (T28910C, C1180580T and C152178T) were found only in Thai NB isolates, including isolates from previous study. Although drug susceptibility tests indicated pan-susceptibility, non-synonymous SNPs previously reported to be associated with resistance to anti-tuberculous drugs; isoniazid, ethambutol, and ethionamide were identified in all the isolates. Non-synonymous SNPs were found in virulence genes such as the genes playing roles in apoptosis inhibition and phagosome arrest. We also report polymorphisms in essential genes, efflux pumps associated genes and genes with known epitopes. The analysis of the TBM isolates and the availability of the variations obtained will provide additional resources for global comparison of isolates from pulmonary tuberculosis and TBM. It will also contribute to the richness of genomic databases towards the prediction of antibiotic resistance, level of virulence and of origin of infection.

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.

Effect of bacillus Calmette-Guérin vaccination on CD4+Foxp3+ T cells during acquired immune response to Mycobacterium tuberculosis infection

Increasing information has shown that many newly emerging strains of Mycobacterium tuberculosis, including the highly prevalent and troublesome Beijing family of strains, can potently induce the emergence of Foxp3(+)CD4 Tregs Although the significance of this is still not fully understood, we have previously provided evidence that the emergence of this population can significantly ablate the protective effect of BCG vaccination, causing progressive fatal disease in the mouse model. However, whether the purpose of this response is to control inflammation or to directly dampen the acquired immune response is still unclear. In the present study, we have shown, using both cell depletion and adoptive transfer strategies, that Tregs can have either properties. Cell depletion resulted in a rapid, but transient, decrease in the lung bacterial load, suggesting release or temporary re-expansion of effector immunity. Transfer of Tregs into Rag2(-/-)or marked congenic mice worsened the disease course and depressed cellular influx of effector T cells into the lungs. Tregs from infected donors seemed to preferentially depress the inflammatory response and granulocytic influx. In contrast, those from BCG-vaccinated and then challenged donors seemed more focused on depression of acquired immunity. These qualitative differences might be related to increasing knowledge reflecting the plasticity of the Treg response.

Mechanisms of Phenotypic Rifampicin Tolerance in Mycobacterium tuberculosis Beijing Genotype Strain B0/W148 Revealed by Proteomics

The "successful" Russian clone B0/W148 of Mycobacterium tuberculosis Beijing is well-known for its capacity to develop antibiotic resistance. During treatment, resistant mutants can occur that have inheritable resistance to specific antibiotics. Next to mutations, M. tuberculosis has several mechanisms that increase their tolerance to a variety of antibiotics. Insights in the phenotypic mechanisms that contribute to drug tolerance will increase our understanding of how antibiotic resistance develops in M. tuberculosis. In this study, we examined the (phospho)proteome dynamics in M. tuberculosis Beijing strain B0/W148 when exposed to a high dose of rifampicin; one of the most potent first-line antibiotics. A total of 2,534 proteins and 191 phosphorylation sites were identified, and revealed the differential regulation of DosR regulon proteins, which are necessary for the development of a dormant phenotype that is less susceptible to antibiotics. By examining independent phenotypic markers of dormancy, we show that persisters of in vitro rifampicin exposure entered a metabolically hypoactive state, which yields rifampicin and other antibiotics largely ineffective. These new insights in the role of protein regulation and post-translational modifications during the initial phase of rifampicin treatment reveal a shortcoming in the antituberculosis regimen that is administered to 8-9 million individuals annually.

Identifying Likely Transmission Pathways within a 10-Year Community Outbreak of Tuberculosis by High-Depth Whole Genome Sequencing

Background

Improved tuberculosis control and the need to contain the spread of drug-resistant strains provide a strong rationale for exploring tuberculosis transmission dynamics at the population level. Whole-genome sequencing provides optimal strain resolution, facilitating detailed mapping of potential transmission pathways.

Methods

We sequenced 22 isolates from a Mycobacterium tuberculosis cluster in New South Wales, Australia, identified during routine 24-locus mycobacterial interspersed repetitive unit typing. Following high-depth paired-end sequencing using the Illumina HiSeq 2000 platform, two independent pipelines were employed for analysis, both employing read mapping onto reference genomes as well as de novo assembly, to control biases in variant detection. In addition to single-nucleotide polymorphisms, the analyses also sought to identify insertions, deletions and structural variants.

Results

Isolates were highly similar, with a distance of 13 variants between the most distant members of the cluster. The most sensitive analysis classified the 22 isolates into 18 groups. Four of the isolates did not appear to share a recent common ancestor with the largest clade; another four isolates had an uncertain ancestral relationship with the largest clade.

Conclusion

Whole genome sequencing, with analysis of single-nucleotide polymorphisms, insertions, deletions, structural variants and subpopulations, enabled the highest possible level of discrimination between cluster members, clarifying likely transmission pathways and exposing the complexity of strain origin. The analysis provides a basis for targeted public health intervention and enhanced classification of future isolates linked to the cluster.

The Role of Adherence and Retreatment in De Novo Emergence of MDR-TB

Treatment failure after therapy of pulmonary tuberculosis (TB) infections is an important challenge, especially when it coincides with de novo emergence of multi-drug-resistant TB (MDR-TB). We seek to explore possible causes why MDR-TB has been found to occur much more often in patients with a history of previous treatment. We develop a mathematical model of the replication of Mycobacterium tuberculosis within a patient reflecting the compartments of macrophages, granulomas, and open cavities as well as parameterizing the effects of drugs on the pathogen dynamics in these compartments. We use this model to study the influence of patient adherence to therapy and of common retreatment regimens on treatment outcome. As expected, the simulations show that treatment success increases with increasing adherence. However, treatment occasionally fails even under perfect adherence due to interpatient variability in pharmacological parameters. The risk of generating MDR de novo is highest between 40% and 80% adherence. Importantly, our simulations highlight the double-edged effect of retreatment: On the one hand, the recommended retreatment regimen increases the overall success rate compared to re-treating with the initial regimen. On the other hand, it increases the probability to accumulate more resistant genotypes. We conclude that treatment adherence is a key factor for a positive outcome, and that screening for resistant strains is advisable after treatment failure or relapse.

Simultaneous drug resistance detection and genotyping of Mycobacterium tuberculosis using a low-density hydrogel microarray

BACKGROUND

Nucleic acid amplification tests are widely used in TB diagnostics. Priority tasks in their development consist of increasing the specificity and sensitivity of the detection of resistance to a wide spectrum of anti-TB drugs.

METHODS

We developed a multiplexed assay allowing the detection of 116 drug resistance-determining mutations in the rpoB, katG, inhA, ahpC, gyrA, gyrB, rrs, eis and embB genes in the Mycobacterium tuberculosis complex genome and six SNPs to identify the main lineages circulating in Russia. The assay is based on the amplification of 17 fragments of the genome using the universal primer adapter technique and heat pulses at the elongation step, followed by hybridization on a microarray.

RESULTS

The method was evaluated using 264 pairs of clinical samples and corresponding isolates. A significant proportion (25%) of smear-negative samples were correctly analysed by microarray analysis in addition to 96% of smear-positive samples. The sensitivity and specificity of the assay exceeded 90% for rifampicin, isoniazid, ofloxacin and second-line injection drugs. In agreement with previous studies, the specificity of ethambutol resistance was as low as 57%, while the sensitivity was 89.9%. Strong association of the Beijing lineage with a resistant phenotype was observed. Euro-American lineage strains, excluding Ural and LAM, were predominantly associated with the susceptible phenotype.

CONCLUSIONS

The developed test has a high sensitivity and specificity and can be directly applied to clinical samples. The combination of mutation-based drug resistance profiling and basic genotyping could be useful for clinical microbiology studies and epidemiological surveillance of the M. tuberculosis complex.

Longitudinal whole genome analysis of pre and post drug treatment Mycobacterium tuberculosis isolates reveals progressive steps to drug resistance

Tuberculosis (TB) is one of the leading causes of death due to an infectious disease in the world. Understanding the mechanisms of drug resistance has become pivotal in the detection and treatment of newly emerging resistant TB cases. We have analyzed three pairs of Mycobacterium tuberculosis strains pre- and post-drug treatment to identify mutations involved in the progression of resistance to the drugs rifampicin and isoniazid. In the rifampicin resistant strain, we confirmed a mutation in rpoB (S450L) that is known to confer resistance to rifampicin. We discovered a novel L101R mutation in the katG gene of an isoniazid resistant strain, which may directly contribute to isoniazid resistance due to the proximity of the mutation to the katG isoniazid-activating site. Another isoniazid resistant strain had a rare mutation in the start codon of katG. We also identified a number of mutations in each longitudinal pair, such as toxin-antitoxin mutations that may influence the progression towards resistance or may play a role in compensatory fitness. These findings improve our knowledge of drug resistance progression during therapy and provide a methodology to monitor longitudinal strains using whole genome sequencing, polymorphism comparison, and functional annotation.

Double-Strand DNA Break Repair in Mycobacteria

Discontinuity of both strands of the chromosome is a lethal event in all living organisms because it compromises chromosome replication. As such, a diversity of DNA repair systems has evolved to repair double-strand DNA breaks (DSBs). In part, this diversity of DSB repair systems has evolved to repair breaks that arise in diverse physiologic circumstances or sequence contexts, including cellular states of nonreplication or breaks that arise between repeats. Mycobacteria elaborate a set of three genetically distinct DNA repair pathways: homologous recombination, nonhomologous end joining, and single-strand annealing. As such, mycobacterial DSB repair diverges substantially from the standard model of prokaryotic DSB repair and represents an attractive new model system. In addition, the presence in mycobacteria of a DSB repair system that can repair DSBs in nonreplicating cells (nonhomologous end joining) or when DSBs arise between repeats (single-strand annealing) has clear potential relevance to Mycobacterium tuberculosis pathogenesis, although the exact role of these systems in M. tuberculosis pathogenesis is still being elucidated. In this article we will review the genetics of mycobacterial DSB repair systems, focusing on recent insights.