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

The complex evolution of antibiotic resistance in Mycobacterium tuberculosis

Multidrug-resistant and extensively drug-resistant tuberculosis (TB) represent a major threat to the control of the disease worldwide. The mechanisms and pathways that result in the emergence and subsequent fixation of resistant strains of Mycobacterium tuberculosis are not fully understood and recent studies suggest that they are much more complex than initially thought. In this review, we highlight the exciting new areas of research within TB resistance that are beginning to fill these gaps in our understanding, whilst also raising new questions and providing future directions.

Precision methylome characterization of Mycobacterium tuberculosis complex (MTBC) using PacBio single-molecule real-time (SMRT) technology

Tuberculosis (TB) remains one of the most common infectious diseases caused by Mycobacterium tuberculosis complex (MTBC). To panoramically analyze MTBC's genomic methylation, we completed the genomes of 12 MTBC strains (Mycobacterium bovis; M. bovis BCG; M. microti; M. africanum; M. tuberculosis H37Rv; H37Ra; and 6 M. tuberculosis clinical isolates) belonging to different lineages and characterized their methylomes using single-molecule real-time (SMRT) technology. We identified three (m6)A sequence motifs and their corresponding methyltransferase (MTase) genes, including the reported mamA, hsdM and a newly discovered mamB. We also experimentally verified the methylated motifs and functions of HsdM and MamB. Our analysis indicated the MTase activities varied between 12 strains due to mutations/deletions. Furthermore, through measuring 'the methylated-motif-site ratio' and 'the methylated-read ratio', we explored the methylation status of each modified site and sequence-read to obtain the 'precision methylome' of the MTBC strains, which enabled intricate analysis of MTase activity at whole-genome scale. Most unmodified sites overlapped with transcription-factor binding-regions, which might protect these sites from methylation. Overall, our findings show enormous potential for the SMRT platform to investigate the precise character of methylome, and significantly enhance our understanding of the function of DNA MTase.

Within-host evolution of bacterial pathogens

Whole-genome sequencing has opened the way for investigating the dynamics and genomic evolution of bacterial pathogens during the colonization and infection of humans. The application of this technology to the longitudinal study of adaptation in an infected host--in particular, the evolution of drug resistance and host adaptation in patients who are chronically infected with opportunistic pathogens--has revealed remarkable patterns of convergent evolution, suggestive of an inherent repeatability of evolution. In this Review, we describe how these studies have advanced our understanding of the mechanisms and principles of within-host genome evolution, and we consider the consequences of findings such as a potent adaptive potential for pathogenicity. Finally, we discuss the possibility that genomics may be used in the future to predict the clinical progression of bacterial infections and to suggest the best option for treatment.

Acquired Drug Resistance in Mycobacterium tuberculosis and Poor Outcomes among Patients with Multidrug-Resistant Tuberculosis

Acquired drug resistance is common and an impediment to successful treatment outcomes.

Rates and risk factors for acquired drug resistance and association with outcomes among patients with multidrug-resistant tuberculosis (MDR TB) are not well defined. In an MDR TB cohort from the country of Georgia, drug susceptibility testing for second-line drugs (SLDs) was performed at baseline and every third month. Acquired resistance was defined as any SLD whose status changed from susceptible at baseline to resistant at follow-up. Among 141 patients, acquired resistance in Mycobacterium tuberculosis was observed in 19 (14%); prevalence was 9.1% for ofloxacin and 9.8% for capreomycin or kanamycin. Baseline cavitary disease and resistance to >6 drugs were associated with acquired resistance. Patients with M. tuberculosis that had acquired resistance were at significantly increased risk for poor treatment outcome compared with patients without these isolates (89% vs. 36%; p<0.01). Acquired resistance occurs commonly among patients with MDR TB and impedes successful treatment outcomes.

Characterization of mutations in streptomycin-resistant Mycobacterium tuberculosis isolates in Sichuan, China and the association between Beijing-lineage and dual-mutation in gidB

Mutations in rpsL, rrs, and gidB are well linked to streptomycin (STR) resistance, some of which are suggested to be potentially associated with Mycobacterium tuberculosis genotypic lineages in certain geographic regions. In this study, we aimed to investigate the mutation characteristics of streptomycin resistance and the relationship between the polymorphism of drug-resistant genes and the lineage of M. tuberculosis isolates in Sichuan, China. A total of 227 M. tuberculosis clinical isolates, including 180 STR-resistant and 47 pan-susceptible isolates, were analyzed for presence of mutations in the rpsL, rrs and gidB loci. Mutation K43R in rpsL was strongly associated with high-level streptomycin resistance (P < 0.01), while mutations in rrs and gidB potentially contributed to low-level resistance (P < 0.05). No general association was exhibited between STR resistance and Beijing genotype, however, in STR-resistant strains, Beijing genotype was significantly correlated with high-level STR resistance, as well as the rpsL mutation K43R (P < 0.01), indicating that Beijing genotype has an evolutionary advantage under streptomycin pressure. Notably, in all isolates of Beijing genotype, a dual mutation E92D (a276c) and A205A (a615g) in gidB was detected, suggesting a highly significant association between this dual mutation and Beijing genotype.

Genome sequence comparisons of serial multi-drug-resistant Mycobacterium tuberculosis isolates over 21 years of infection in a single patient

We report a case of chronic pulmonary multi-drug-resistant tuberculosis. Despite 14 years of treatment, Mycobacterium tuberculosis was persistently isolated from sputum. Following treatment cessation the patient remained well, although M. tuberculosis was isolated from sputum for a further 8 years. Genome sequencing of eight serial M. tuberculosis isolates cultured between 1991 and 2011 revealed 17 mutations (0.8 mutations per genome year^- 1). Eight of these were persisting mutations and only two mutations were detected in the 7 years following cessation of treatment in 2004. In four isolates there were mixed alleles, suggesting the likely presence of bacterial subpopulations. The initial 1991 isolate demonstrated genotypic resistance to isoniazid (katG W91R), rifampicin (rpoB S531L), ethambutol (embB M306V), streptomycin (gidB L16R), quinolones (gyrA S95T) and P-aminosalicylic acid (thyA T202A). Subsequent resistance mutations developed for pyrazinamide (pncA I31F) and ethionamide (ethA frameshift). Such information might have been instructive when developing a treatment regimen. In retrospect and with the benefit of high-resolution genomic hindsight we were able to determine that the patient received only one or two active anti-tuberculous agents for most of their treatment. Additionally, mutations in bacA and Rv2326c were detected, which may have contributed to the persistent but mild disease course. BacA is likely to be associated with maintenance of chronic infection and Rv2326c with a decreased bacterial metabolic state. These results expand our understanding of M. tuberculosis evolution during human infection and underline the link between antibiotic resistance and clinical persistence.

Diversity of Mycobacterium tuberculosis across Evolutionary Scales

Tuberculosis (TB) is a global public health emergency. Increasingly drug resistant strains of Mycobacterium tuberculosis (M.tb) continue to emerge and spread, highlighting adaptability of this pathogen. Most studies of M.tb evolution have relied on ‘between-host’ samples, in which each person with TB is represented by a single M.tb isolate. However, individuals with TB commonly harbor populations of M.tb numbering in the billions. Here, we use analyses of M.tb genomic data from within and between hosts to gain insight into influences shaping genetic diversity of this pathogen. We find that the amount of M.tb genetic diversity harbored by individuals with TB can vary dramatically, likely as a function of disease severity. Surprisingly, we did not find an appreciable impact of TB treatment on M.tb diversity. In examining genomic data from M.tb samples within and between hosts with TB, we find that genes involved in the regulation, synthesis, and transportation of immunomodulatory cell envelope lipids appear repeatedly in the extremes of various statistical measures of diversity. Many of these genes have been identified as possible targets of selection in other studies employing different methods and data sets. Taken together, these observations suggest that M.tb cell envelope lipids are targets of selection within hosts. Many of these lipids are specific to pathogenic mycobacteria and, in some cases, human-pathogenic mycobacteria. We speculate that rapid adaptation of cell envelope lipids is facilitated by functional redundancy, flexibility in their metabolism, and their roles mediating interactions with the host.

Tuberculosis (TB) is a grave threat to global public health and is the second leading cause of death due to infectious disease. The causative agent, Mycobacterium tuberculosis (M.tb), has emerged in increasingly drug resistant forms that hamper our efforts to control TB. We need a better understanding of M.tb adaptation to guide development of more effective TB treatment and control strategies. The goal of this study was to gain insight into M.tb evolution within individual patients with TB. We found that TB patients harbor a diverse population of M.tb. We further found evidence to suggest that the bacterial population evolves measurably in response to selection pressures imposed by the environment within hosts. Changes were particularly notable in M.tb genes involved in the regulation, synthesis, and transportation of lipids and glycolipids of the bacterial cell envelope. These findings have important implications for drug and vaccine development, and provide insight into TB host pathogen interactions.

Genetic Structuration, Demography and Evolutionary History of Mycobacterium tuberculosis LAM9 Sublineage in the Americas as Two Distinct Subpopulations Revealed by Bayesian Analyses

Tuberculosis (TB) remains broadly present in the Americas despite intense global efforts for its control and elimination. Starting from a large dataset comprising spoligotyping (n = 21183 isolates) and 12-loci MIRU-VNTRs data (n = 4022 isolates) from a total of 31 countries of the Americas (data extracted from the SITVIT2 database), this study aimed to get an overview of lineages circulating in the Americas. A total of 17119 (80.8%) strains belonged to the Euro-American lineage 4, among which the most predominant genotypic family belonged to the Latin American and Mediterranean (LAM) lineage (n = 6386, 30.1% of strains). By combining classical phylogenetic analyses and Bayesian approaches, this study revealed for the first time a clear genetic structuration of LAM9 sublineage into two subpopulations named LAM9C1 and LAM9C2, with distinct genetic characteristics. LAM9C1 was predominant in Chile, Colombia and USA, while LAM9C2 was predominant in Brazil, Dominican Republic, Guadeloupe and French Guiana. Globally, LAM9C2 was characterized by higher allelic richness as compared to LAM9C1 isolates. Moreover, LAM9C2 sublineage appeared to expand close to twenty times more than LAM9C1 and showed older traces of expansion. Interestingly, a significant proportion of LAM9C2 isolates presented typical signature of ancestral LAM-RD^Rio MIRU-VNTR type (224226153321). Further studies based on Whole Genome Sequencing of LAM strains will provide the needed resolution to decipher the biogeographical structure and evolutionary history of this successful family.

Risk Factors for Acquired Rifamycin and Isoniazid Resistance: A Systematic Review and Meta-Analysis

BACKGROUND

Studies looking at acquired drug resistance (ADR) are diverse with respect to geographical distribution, HIV co-infection rates, retreatment status and programmatic factors such as regimens administered and directly observed therapy. Our objective was to examine and consolidate evidence from clinical studies of the multifactorial aetiology of acquired rifamycin and/or isoniazid resistance within the scope of a single systematic review. This is important to inform policy and identify key areas for further studies.

METHODS

Case-control and cohort studies and randomised controlled trials that reported ADR as an outcome during antitubercular treatment regimens including a rifamycin and examined the association of at least 1 risk factor were included. Post hoc, we carried out random effects Mantel-Haenszel weighted meta-analyses of the impact of 2 key risk factors 1) HIV and 2) baseline drug resistance on the binary outcome of ADR. Heterogeneity was assessed used I2 statistic. As a secondary outcome, we calculated median cumulative incidence of ADR, weighted by the sample size of the studies.

RESULTS

Meta-analysis of 15 studies showed increased risk of ADR with baseline mono- or polyresistance (RR 4.85 95% CI 3.26 to 7.23, heterogeneity I2 58%, 95% CI 26 to 76%). Meta-analysis of 8 studies showed that HIV co-infection was associated with increased risk of ADR (RR 3.02, 95% CI 1.28 to 7.11); there was considerable heterogeneity amongst these studies (I2 81%, 95% CI 64 to 90%). Non-adherence, extrapulmonary/disseminated disease and advanced immunosuppression in HIV co-infection were other risk factors noted. The weighted median cumulative incidence of acquired multi drug resistance calculated in 24 studies (assuming whole cohort as denominator, regardless of follow up DST) was 0.1% (5th to 95th percentile 0.07 to 3.2%).

CONCLUSION

Baseline drug resistance and HIV co-infection were significant risk factors for ADR. There was a trend of positive association with non-adherence which is likely to contribute to the outcome of ADR. The multifactorial aetiology of ADR in a programmatic setting should be further evaluated via appropriately designed studies.

Evolution of Extensively Drug-Resistant Tuberculosis over Four Decades: Whole Genome Sequencing and Dating Analysis of Mycobacterium tuberculosis Isolates from KwaZulu-Natal

BACKGROUND

The continued advance of antibiotic resistance threatens the treatment and control of many infectious diseases. This is exemplified by the largest global outbreak of extensively drug-resistant (XDR) tuberculosis (TB) identified in Tugela Ferry, KwaZulu-Natal, South Africa, in 2005 that continues today. It is unclear whether the emergence of XDR-TB in KwaZulu-Natal was due to recent inadequacies in TB control in conjunction with HIV or other factors. Understanding the origins of drug resistance in this fatal outbreak of XDR will inform the control and prevention of drug-resistant TB in other settings. In this study, we used whole genome sequencing and dating analysis to determine if XDR-TB had emerged recently or had ancient antecedents.

METHODS AND FINDINGS

We performed whole genome sequencing and drug susceptibility testing on 337 clinical isolates of Mycobacterium tuberculosis collected in KwaZulu-Natal from 2008 to 2013, in addition to three historical isolates, collected from patients in the same province and including an isolate from the 2005 Tugela Ferry XDR outbreak, a multidrug-resistant (MDR) isolate from 1994, and a pansusceptible isolate from 1995. We utilized an array of whole genome comparative techniques to assess the relatedness among strains, to establish the order of acquisition of drug resistance mutations, including the timing of acquisitions leading to XDR-TB in the LAM4 spoligotype, and to calculate the number of independent evolutionary emergences of MDR and XDR. Our sequencing and analysis revealed a 50-member clone of XDR M. tuberculosis that was highly related to the Tugela Ferry XDR outbreak strain. We estimated that mutations conferring isoniazid and streptomycin resistance in this clone were acquired 50 y prior to the Tugela Ferry outbreak (katG S315T [isoniazid]; gidB 130 bp deletion [streptomycin]; 1957 [95% highest posterior density (HPD): 1937-1971]), with the subsequent emergence of MDR and XDR occurring 20 y (rpoB L452P [rifampicin]; pncA 1 bp insertion [pyrazinamide]; 1984 [95% HPD: 1974-1992]) and 10 y (rpoB D435G [rifampicin]; rrs 1400 [kanamycin]; gyrA A90V [ofloxacin]; 1995 [95% HPD: 1988-1999]) prior to the outbreak, respectively. We observed frequent de novo evolution of MDR and XDR, with 56 and nine independent evolutionary events, respectively. Isoniazid resistance evolved before rifampicin resistance 46 times, whereas rifampicin resistance evolved prior to isoniazid only twice. We identified additional putative compensatory mutations to rifampicin in this dataset. One major limitation of this study is that the conclusions with respect to ordering and timing of acquisition of mutations may not represent universal patterns of drug resistance emergence in other areas of the globe.

CONCLUSIONS

In the first whole genome-based analysis of the emergence of drug resistance among clinical isolates of M. tuberculosis, we show that the ancestral precursor of the LAM4 XDR outbreak strain in Tugela Ferry gained mutations to first-line drugs at the beginning of the antibiotic era. Subsequent accumulation of stepwise resistance mutations, occurring over decades and prior to the explosion of HIV in this region, yielded MDR and XDR, permitting the emergence of compensatory mutations. Our results suggest that drug-resistant strains circulating today reflect not only vulnerabilities of current TB control efforts but also those that date back 50 y. In drug-resistant TB, isoniazid resistance was overwhelmingly the initial resistance mutation to be acquired, which would not be detected by current rapid molecular diagnostics employed in South Africa that assess only rifampicin resistance.

Quantifying the impact of drug combination regimens on TB treatment efficacy and multidrug resistance probability

OBJECTIVES

TB patients' non-adherence to the multidrug treatment regimen is thought to be the main cause of the emergence of drug resistance. The purpose of this study was to quantify the impacts of two-drug combination regimens and non-adherence to these regimens on treatment efficacy and drug resistance probability.

METHODS

A drug treatment modelling strategy was developed by incorporating a pharmacokinetic/pharmacodynamic model into a bacterial population dynamic model to explore the dynamics of TB bacilli and evolution of resistance during multidrug combination therapy, with an emphasis on non-adherence. A Hill-equation-based pharmacodynamic model was used to assess the bactericidal efficacy of single drugs and to estimate drug interactions.

RESULTS

Non-adherence to the treatment regimen increased treatment duration by nearly 1.6- and 3.4-fold relative to compliance with treatment. Symptom-based intermittent treatment, a form of non-adherence, might lead to treatment failure and accelerated growth and evolution of resistant mutants, resulting in a dramatically higher probability of 4.17 × 10(-3) (95% CI 2.10 × 10(-4)-1.28 × 10(-2)) for the emergence of MDR TB. Overall, determination of the optimal treatment regimen depended on the different types of medication adherence.

CONCLUSIONS

Our model not only predicts evolutionary dynamics, but also quantifies treatment efficacy. More broadly, our model provides a quantitative framework for improving treatment protocols and establishing an emergence threshold of resistance that can be used to prevent drug resistance.

Beijing clades of Mycobacterium tuberculosis are associated with differential survival in HIV-negative Russian patients

We conducted a prospective study to establish factors associated with survival in tuberculosis patients in Russia including social, clinical and pathogen-related genetic parameters. Specifically we wished to determine whether different strains/clades of the Beijing lineage exerted a differential effect of survival. HIV-negative culture-confirmed cases were recruited during 2008-2010 across Samara Oblast and censored in December 2011. Molecular characterization was performed by a combination of spoligotyping, multilocus VNTR typing and whole genome sequencing (WGS). We analyzed 2602 strains and detected a high prevalence of Beijing family (n=1933; 74%) represented largely by two highly homogenous dominant clades A (n=794) and B (n=402) and non-A/non-B (n=737). Multivariable analysis of 1366 patients with full clinical and genotyping data showed that multi- and extensive drug resistance (HR=1.86; 95%CI: 1.52, 2.28 and HR=2.19; 95%CI: 1.55, 3.11) had the largest impact on survival. In addition older age, extensive lung damage, shortness of breath, treatment in the past and alcohol abuse reduced survival time. After adjustment for clinical and demographic predictors there was evidence that clades A and B combined were associated with poorer survival than other Beijing strains (HR=0.48; 95%CI 0.34, 0.67). All other pathogen-related factors (polymorphisms in genes plcA, plcB, plcC, lipR, dosT and pks15/1) had no effect on survival. In conclusion, drug resistance exerted the greatest effect on survival of TB patients. Nevertheless we provide evidence for the independent biological effect on survival of different Beijing family strains even within the same defined geographical population. Better understanding of the role of different strain factors in active disease and their influence on outcome is essential.

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.

Analysis of IS6110 insertion sites provide a glimpse into genome evolution of Mycobacterium tuberculosis

Insertion sequence (IS) 6110 is found at multiple sites in the Mycobacterium tuberculosis genome and displays a high degree of polymorphism with respect to copy number and insertion sites. Therefore, IS6110 is considered to be a useful molecular marker for diagnosis and strain typing of M. tuberculosis. Generally IS6110 elements are identified using experimental methods, useful for analysis of a limited number of isolates. Since short read genome sequences generated using next-generation sequencing (NGS) platforms are available for a large number of isolates, a computational pipeline for identification of IS6110 elements from these datasets was developed. This study shows results from analysis of NGS data of 1377 M. tuberculosis isolates. These isolates represent all seven major global lineages of M. tuberculosis. Lineage specific copy number patterns and preferential insertion regions were observed. Intra-lineage differences were further analyzed for identifying spoligotype specific variations. Copy number distribution and preferential locations of IS6110 in different lineages imply independent evolution of IS6110, governed mainly through ancestral insertion, fitness (gene truncation, promoter activity) and recombinational loss of some copies. A phylogenetic tree based on IS6110 insertion data of different isolates was constructed in order to understand genome level variations of different markers across different lineages.

Synonymous Genetic Variation in Natural Isolates of Escherichia coli Does Not Predict Where Synonymous Substitutions Occur in a Long-Term Experiment

Synonymous genetic differences vary by more than 20-fold among genes in natural isolates of Escherichia coli. One hypothesis to explain this heterogeneity is that genes with high levels of synonymous variation mutate at higher rates than genes with low synonymous variation. If so, then one would expect to observe similar mutational patterns in evolution experiments. In fact, however, the pattern of synonymous substitutions in a long-term evolution experiment with E. coli does not support this hypothesis. In particular, the extent of synonymous variation across genes in that experiment does not reflect the variation observed in natural isolates of E. coli. Instead, gene length alone predicts with high accuracy the prevalence of synonymous changes in the experimental populations. We hypothesize that patterns of synonymous variation in natural E. coli populations are instead caused by differences across genomic regions in their effective population size that, in turn, reflect different histories of recombination, horizontal gene transfer, selection, and population structure.

Genomic epidemiology of multidrug-resistant Mycobacterium tuberculosis during transcontinental spread

The transcontinental spread of multidrug-resistant (MDR) tuberculosis is poorly characterized in molecular epidemiologic studies. We used genomic sequencing to understand the establishment and dispersion of MDR Mycobacterium tuberculosis within a group of immigrants to the United States. We used a genomic epidemiology approach to study a genotypically matched (by spoligotype, IS6110 restriction fragment length polymorphism, and mycobacterial interspersed repetitive units-variable number of tandem repeat signature) lineage 2/Beijing MDR strain implicated in an outbreak of tuberculosis among refugees in Thailand and consecutive cases within California. All 46 MDR M. tuberculosis genomes from both Thailand and California were highly related, with a median difference of 10 single-nucleotide polymorphisms (SNPs). The Wat Tham Krabok (WTK) strain is a new sequence type distinguished from all known Beijing strains by 55 SNPs and a genomic deletion (Rv1267c) associated with increased fitness. Sequence data revealed a highly prevalent MDR strain that included several closely related but distinct allelic variants within Thailand, rather than the occurrence of a single outbreak. In California, sequencing data supported multiple independent introductions of WTK with subsequent transmission and reactivation within the state, as well as a potential super spreader with a prolonged infectious period. Twenty-seven drug resistance-conferring mutations and 4 putative compensatory mutations were found within WTK strains. Genomic sequencing has substantial epidemiologic value in both low- and high-burden settings in understanding transmission chains of highly prevalent MDR strains.

Tackling tuberculosis: Insights from an international TB Summit in London

Tuberculosis (TB) poses a grave predicament to the world as it is not merely a scientific challenge but a socio-economic burden as well. A prime cause of mortality in human due to an infectious disease; the malady and its cause, Mycobacterium tuberculosis have remained an enigma with many questions that remain unanswered. The ability of the pathogen to survive and switch between varied physiological states necessitates a protracted therapeutic regimen that exerts an excessive strain on low-resource countries. To complicate things further, there has been a significant rise of antimicrobial resistance. Existing control measures, including treatment regimens have remained fairly uniform globally for at least half a century and require reinvention. Overcoming the societal and scientific challenges requires an increase in dialog to identify key regions that need attention and effective partners with whom successful collaborations can be fostered. In this report, we explore the discussions held at the International TB Summit 2015 hosted by EuroSciCon, which served as an excellent platform for researchers to share their recent findings. Ground-breaking results require outreach to affect policy design, governance and control of the disease. Hence, we feel it is important that meetings such as these reach a wider, global audience.

Southern East Asian origin and coexpansion of Mycobacterium tuberculosis Beijing family with Han Chinese

The Beijing family is the most successful genotype of Mycobacterium tuberculosis and responsible for more than a quarter of the global tuberculosis epidemic. As the predominant genotype in East Asia, the Beijing family has been emerging in various areas of the world and is often associated with disease outbreaks and antibiotic resistance. Revealing the origin and historical dissemination of this strain family is important for understanding its current global success. Here we characterized the global diversity of this family based on whole-genome sequences of 358 Beijing strains. We show that the Beijing strains endemic in East Asia are genetically diverse, whereas the globally emerging strains mostly belong to a more homogenous subtype known as "modern" Beijing. Phylogeographic and coalescent analyses indicate that the Beijing family most likely emerged around 30,000 y ago in southern East Asia, and accompanied the early colonization by modern humans in this area. By combining the genomic data and genotyping result of 1,793 strains from across China, we found the "modern" Beijing sublineage experienced massive expansions in northern China during the Neolithic era and subsequently spread to other regions following the migration of Han Chinese. Our results support a parallel evolution of the Beijing family and modern humans in East Asia. The dominance of the "modern" Beijing sublineage in East Asia and its recent global emergence are most likely driven by its hypervirulence, which might reflect adaption to increased human population densities linked to the agricultural transition in northern China.

Differential transcriptomic and metabolic profiles of M. africanum- and M. tuberculosis-infected patients after, but not before, drug treatment

The epidemiology of Mycobacterium tuberculosis (Mtb) and M. africanum (Maf) suggest differences in their virulence, but the host immune profile to better understand the pathogenesis of tuberculosis (TB) have not been studied. We compared the transcriptomic and metabolic profiles between Mtb and Maf-infected TB cases to identify host biomarkers associated with lineages-specific pathogenesis and response to anti-TB chemotherapy. Venous blood samples from Mtb- and Maf-infected patients obtained before and after anti-TB treatment were analysed for cell composition, gene expression and metabolic profiles. Prior to treatment, similar transcriptomic profiles were seen in Maf- and Mtb-infected patients. In contrast, post-treatment, over 1600 genes related to immune responses and metabolic diseases were differentially expressed between the groups. Notably, the upstream regulator hepatocyte nuclear factor 4-alpha (HNF4α), which regulated 15% of these genes, was markedly enriched. Serum metabolic profiles were similar in both group pre-treatment, but the decline in pro-inflammatory metabolites post-treatment were most pronounced in Mtb-infected patients. Together, the differences in both peripheral blood transcriptomic and serum metabolic profiles between Maf- and Mtb-infected patients observed over the treatment period, might be indicative of intrinsic host factors related to susceptibility to TB and/or differential efficacy of the standard anti-TB treatment on the two lineages.

Evolution of drug resistance in tuberculosis: recent progress and implications for diagnosis and therapy

Drug-resistant tuberculosis is a growing threat to global public health. Recent efforts to understand the evolution of drug resistance have shown that changes in drug–target interactions are only the first step in a longer adaptive process. The emergence of transmissible drug-resistant Mycobacterium tuberculosis is the result of a multitude of additional genetic mutations, many of which interact, a phenomenon known as epistasis. The varied effects of these epistatic interactions include compensating for the reduction of the biological cost associated with the development of drug resistance, increasing the level of resistance, and possibly accommodating broader changes in the physiology of resistant bacteria. Knowledge of these processes and our ability to detect them as they happen informs the development of diagnostic tools and better control strategies. In particular, the use of whole genome sequencing combined with surveillance efforts in the field could provide a powerful instrument to prevent future epidemics of drug-resistant tuberculosis.

Host-Directed Therapies for Tuberculosis

Host-directed therapies are a relatively new and promising approach to treatment of tuberculosis. Modulation of specific host immune pathways, including those that impact inflammation and immunopathology, can limit mycobacterial infection and pathology, both in cell culture and in animal models. This review explores a range of host pathways and drugs, some already approved for clinical use that have the potential to provide new adjunctive therapies for tuberculosis. Drugs targeting host processes may largely avoid the development of bacterial antibiotic resistance, a major public health concern for tuberculosis. However, these drugs may also have generally increased risk for side effects on the host. Understanding the specific mechanisms by which these drugs act and the relationship of these mechanisms to Mycobacterium tuberculosis pathogenesis will be critical in selecting appropriate host-directed therapy. Overall, these host-directed compounds provide a novel strategy for antituberculosis therapy.

The Association between Mycobacterium Tuberculosis Genotype and Drug Resistance in Peru

Background

The comparison of Mycobacterium tuberculosis bacterial genotypes with phenotypic, demographic, geospatial and clinical data improves our understanding of how strain lineage influences the development of drug-resistance and the spread of tuberculosis.

Methods

To investigate the association of Mycobacterium tuberculosis bacterial genotype with drug-resistance. Drug susceptibility testing together with genotyping using both 15-loci MIRU-typing and spoligotyping, was performed on 2,139 culture positive isolates, each from a different patient in Lima, Peru. Demographic, geospatial and socio-economic data were collected using questionnaires, global positioning equipment and the latest national census.

Results

The Latin American Mediterranean (LAM) clade (OR 2.4, p<0.001) was significantly associated with drug-resistance and alone accounted for more than half of all drug resistance in the region. Previously treated patients, prisoners and genetically clustered cases were also significantly associated with drug-resistance (OR's 2.5, 2.4 and 1.8, p<0.001, p<0.05, p<0.001 respectively).

Conclusions

Tuberculosis disease caused by the LAM clade was more likely to be drug resistant independent of important clinical, genetic and socio-economic confounding factors. Explanations for this include; the preferential co-evolution of LAM strains in a Latin American population, a LAM strain bacterial genetic background that favors drug-resistance or the "founder effect" from pre-existing LAM strains disproportionately exposed to drugs.

PhyTB: Phylogenetic tree visualisation and sample positioning for M. tuberculosis

Background

Phylogenetic-based classification of M. tuberculosis and other bacterial genomes is a core analysis for studying evolutionary hypotheses, disease outbreaks and transmission events. Whole genome sequencing is providing new insights into the genomic variation underlying intra- and inter-strain diversity, thereby assisting with the classification and molecular barcoding of the bacteria. One roadblock to strain investigation is the lack of user-interactive solutions to interrogate and visualise variation within a phylogenetic tree setting.

Results

We have developed a web-based tool called PhyTB (http://pathogenseq.lshtm.ac.uk/phytblive/index.php) to assist phylogenetic tree visualisation and identification of M. tuberculosis clade-informative polymorphism. Variant Call Format files can be uploaded to determine a sample position within the tree. A map view summarises the geographical distribution of alleles and strain-types. The utility of the PhyTB is demonstrated on sequence data from 1,601 M. tuberculosis isolates.

Conclusion

PhyTB contextualises M. tuberculosis genomic variation within epidemiological, geographical and phylogenic settings. Further tool utility is possible by incorporating large variants and phenotypic data (e.g. drug-resistance profiles), and an assessment of genotype-phenotype associations. Source code is available to develop similar websites for other organisms (http://sourceforge.net/projects/phylotrack).

A microbiological revolution meets an ancient disease: improving the management of tuberculosis with genomics

Tuberculosis (TB) is an ancient disease with an enormous global impact. Despite declining global incidence, the diagnosis, phenotyping, and epidemiological investigation of TB require significant clinical microbiology laboratory resources. Current methods for the detection and characterization of Mycobacterium tuberculosis consist of a series of laboratory tests varying in speed and performance, each of which yields incremental information about the disease. Since the sequencing of the first M. tuberculosis genome in 1998, genomic tools have aided in the diagnosis, treatment, and control of TB. Here we summarize genomics-based methods that are positioned to be introduced in the modern clinical TB laboratory, and we highlight how recent advances in genomics will improve the detection of antibiotic resistance-conferring mutations and the understanding of M. tuberculosis transmission dynamics and epidemiology. We imagine the future TB clinic as one that relies heavily on genomic interrogation of the M. tuberculosis isolate, allowing for more rapid diagnosis of TB and real-time monitoring of outbreak emergence.

DNA replication fidelity in Mycobacterium tuberculosis is mediated by an ancestral prokaryotic proofreader

The DNA replication machinery is an important target for antibiotic development for increasingly drug resistant bacteria including Mycobacterium tuberculosis¹. While blocking DNA replication leads to cell death, disrupting the processes used to ensure replication fidelity can accelerate mutation and the evolution of drug resistance. In E. coli, the proofreading subunit of the replisome, the ε-exonuclease, is essential for high fidelity DNA replication²; however, we find that it is completely dispensable in M. tuberculosis. Rather, the mycobacterial replicative polymerase, DnaE1, encodes a novel editing function that proofreads DNA replication, mediated by an intrinsic 3′-5′ exonuclease activity within its PHP domain. Inactivation of the DnaE1 PHP domain increases the mutation rate by greater than 3,000 fold. Moreover, phylogenetic analysis of DNA replication proofreading in the bacterial kingdom suggests that E. coli is a phylogenetic outlier and that PHP-domain mediated proofreading is widely conserved and indeed may be the ancestral prokaryotic proofreader.

Totally drug-resistant tuberculosis and adjunct therapies

The first cases of totally drug-resistant (TDR) tuberculosis (TB) were reported in Italy 10 years ago; more recently, cases have also been reported in Iran, India and South Africa. Although there is no consensus on terminology, it is most commonly described as 'resistance to all first- and second-line drugs used to treat TB'. Mycobacterium tuberculosis (M.tb) acquires drug resistance mutations in a sequential fashion under suboptimal drug pressure due to monotherapy, inadequate dosing, treatment interruptions and drug interactions. The treatment of TDR-TB includes antibiotics with disputed or minimal effectiveness against M.tb, and the fatality rate is high. Comorbidities such as diabetes and infection with human immunodeficiency virus further impact on TB treatment options and survival rates. Several new drug candidates with novel modes of action are under late-stage clinical evaluation (e.g., delamanid, bedaquiline, SQ109 and sutezolid). 'Repurposed' antibiotics have also recently been included in the treatment of extensively drug resistant TB. However, because of mutations in M.tb, drugs will not provide a cure for TB in the long term. Adjunct TB therapies, including therapeutic vaccines, vitamin supplementation and/or repurposing of drugs targeting biologically and clinically relevant molecular pathways, may achieve better clinical outcomes in combination with standard chemotherapy. Here, we review broader perspectives of drug resistance in TB and potential adjunct treatment options.

Construction of a virtual Mycobacterium tuberculosis consensus genome and its application to data from a next generation sequencer

BACKGROUND

Although Mycobacterium tuberculosis isolates are consisted of several different lineages and the epidemiology analyses are usually assessed relative to a particular reference genome, M. tuberculosis H37Rv, which might introduce some biased results. Those analyses are essentially based genome sequence information of M. tuberculosis and could be performed in sillico in theory, with whole genome sequence (WGS) data available in the databases and obtained by next generation sequencers (NGSs). As an approach to establish higher resolution methods for such analyses, whole genome sequences of the M. tuberculosis complexes (MTBCs) strains available on databases were aligned to construct virtual reference genome sequences called the consensus sequence (CS), and evaluated its feasibility in in sillico epidemiological analyses.

RESULTS

The consensus sequence (CS) was successfully constructed and utilized to perform phylogenetic analysis, evaluation of read mapping efficacy, which is crucial for detecting single nucleotide polymorphisms (SNPs), and various MTBC typing methods virtually including spoligotyping, VNTR, Long sequence polymorphism and Beijing typing. SNPs detected based on CS, in comparison with H37Rv, were utilized in concatemer-based phylogenetic analysis to determine their reliability relative to a phylogenetic tree based on whole genome alignment as the gold standard. Statistical comparison of phylogenic trees based on CS with that of H37Rv indicated the former showed always better results that that of later. SNP detection and concatenation with CS was advantageous because the frequency of crucial SNPs distinguishing among strain lineages was higher than those of H37Rv. The number of SNPs detected was lower with the consensus than with the H37Rv sequence, resulting in a significant reduction in computational time. Performance of each virtual typing was satisfactory and accorded with those published when those are available.

CONCLUSIONS

These results indicated that virtual CS constructed from genome sequence data is an ideal approach as a reference for MTBC studies.

Large-scale whole genome sequencing of M. tuberculosis provides insights into transmission in a high prevalence area

To improve understanding of the factors influencing tuberculosis transmission and the role of pathogen variation, we sequenced all available specimens from patients diagnosed over 15 years in a whole district in Malawi. Mycobacterium tuberculosis lineages were assigned and transmission networks constructed, allowing ≤10 single nucleotide polymorphisms (SNPs) difference. We defined disease as due to recent infection if the network-determined source was within 5 years, and assessed transmissibility from forward transmissions resulting in disease. High-quality sequences were available for 1687 disease episodes (72% of all culture-positive episodes): 66% of patients linked to at least one other patient. The between-patient mutation rate was 0.26 SNPs/year (95% CI 0.21–0.31). We showed striking differences by lineage in the proportion of disease due to recent transmission and in transmissibility (highest for lineage-2 and lowest for lineage-1) that were not confounded by immigration, HIV status or drug resistance. Transmissions resulting in disease decreased markedly over time.

DOI:http://dx.doi.org/10.7554/eLife.05166.001

Tuberculosis is an important public health threat around the globe and is particularly common in developing countries. It is difficult to control the spread of the disease because the bacteria that cause it can spread when an infected individual coughs or sneezes. It may take years for an infected individual to develop symptoms of tuberculosis so it can be hard to trace the source of an outbreak, and people infected with HIV are particularly susceptible to the disease.

The bacterium that causes the majority of cases of tuberculosis is called Mycobacterium tuberculosis. There are several different varieties or ‘lineages’ of M. tuberculosis, and it is thought that they may vary in their ability to spread and cause disease. However, the results of previous studies have been inconsistent and there also seems to be a lot of variation between strains within the same lineage.

In this study, Guerra-Assunção et al. used an approach called whole genome sequencing alongside more traditional methods to study the spread of tuberculosis in Malawi. They sequenced the genomes of every available sample of M. tuberculosis collected from patients in the Karonga district of Malawi over a 15-year period. This produced high-quality DNA sequence data about the bacteria responsible for almost 1700 cases of disease.

Using this massive amount of data, Guerra-Assunção et al. constructed networks that showed how the bacteria had spread in the community. This revealed that there were differences between the ability of the various M. tuberculosis lineages to cause disease and to spread in communities. For example, lineage 1 was less likely than the other lineages to cause disease soon after infecting an individual and was less able to spread.

The data also show that the proportion of cases of disease due to recent infection declined substantially during the 15-year period. This indicates that the tuberculosis and HIV control programmes in the area have been successful.

Guerra-Assunção et al.'s findings show that it is possible to understand how tuberculosis is transmitted on a large scale. The next challenge is to understand why the lineages differ in their ability to cause disease and spread between individuals.

DOI:http://dx.doi.org/10.7554/eLife.05166.002

Origin and proliferation of multiple-drug resistance in bacterial pathogens

SUMMARY

Many studies report the high prevalence of multiply drug-resistant (MDR) strains. Because MDR infections are often significantly harder and more expensive to treat, they represent a growing public health threat. However, for different pathogens, different underlying mechanisms are traditionally used to explain these observations, and it is unclear whether each bacterial taxon has its own mechanism(s) for multidrug resistance or whether there are common mechanisms between distantly related pathogens. In this review, we provide a systematic overview of the causes of the excess of MDR infections and define testable predictions made by each hypothetical mechanism, including experimental, epidemiological, population genomic, and other tests of these hypotheses. Better understanding the cause(s) of the excess of MDR is the first step to rational design of more effective interventions to prevent the origin and/or proliferation of MDR.

Modelling challenges in context: lessons from malaria, HIV, and tuberculosis

Malaria, HIV, and tuberculosis (TB) collectively account for several million deaths each year, with all three ranking among the top ten killers in low-income countries. Despite being caused by very different organisms, malaria, HIV, and TB present a suite of challenges for mathematical modellers that are particularly pronounced in these infections, but represent general problems in infectious disease modelling, and highlight many of the challenges described throughout this issue. Here, we describe some of the unifying challenges that arise in modelling malaria, HIV, and TB, including variation in dynamics within the host, diversity in the pathogen, and heterogeneity in human contact networks and behaviour. Through the lens of these three pathogens, we provide specific examples of the other challenges in this issue and discuss their implications for informing public health efforts.

Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage

Mycobacterium tuberculosis strains of the Beijing lineage are globally distributed and are associated with the massive spread of multidrug-resistant (MDR) tuberculosis in Eurasia. Here we reconstructed the biogeographical structure and evolutionary history of this lineage by genetic analysis of 4,987 isolates from 99 countries and whole-genome sequencing of 110 representative isolates. We show that this lineage initially originated in the Far East, from where it radiated worldwide in several waves. We detected successive increases in population size for this pathogen over the last 200 years, practically coinciding with the Industrial Revolution, the First World War and HIV epidemics. Two MDR clones of this lineage started to spread throughout central Asia and Russia concomitantly with the collapse of the public health system in the former Soviet Union. Mutations identified in genes putatively under positive selection and associated with virulence might have favored the expansion of the most successful branches of the lineage.

Single-bacterium nanomechanics in biomedicine: unravelling the dynamics of bacterial cells

The use of the atomic force microscope (AFM) in microbiology has progressed significantly throughout the years since its first application as a high-resolution imaging instrument. Modern AFM setups are capable of characterizing the nanomechanical behaviour of bacterial cells at both the cellular and molecular levels, where elastic properties and adhesion forces of single bacterium cells can be examined under different experimental conditions. Considering that bacterial and biofilm-mediated infections continue to challenge the biomedical field, it is important to understand the biophysical events leading towards bacterial adhesion and colonization on both biological and non-biological substrates. The purpose of this review is to present the latest findings concerning the field of single-bacterium nanomechanics, and discuss future trends and applications of nanoindentation and single-cell force spectroscopy techniques in biomedicine.

Human and tuberculosis co-evolution: An integrative view

Tuberculosis (TB) ranks as the second cause of death from an infectious disease worldwide after HIV. Archaeogenetics and evolutionary scenario for the Mycobacterium tuberculosis complex (MTBC) are in favor of a long-term interaction between tuberculosis and humans, predating the Neolithic period, contrary to the traditional belief. If tuberculosis evolved as a human pathogen in Africa and has spread outside Africa about more than ten-thousand years ago, its life history traits have been shaped by the immune system. Numerous studies described a variety of human susceptibility factors to TB, suggesting that MTBC strains have evolved different ways to overcome this system. However, the results of these studies reveal some inconsistencies even within populations. The temporally varying history of epidemics and ever-varying genetic diversity of pathogens and strains could easily contribute to blur out signal of selection in our human genome. Palaeomicrobiology gives the opportunity to genotype ancient TB strains circulating in past populations. Accessing ancient human pathogens allows us to a better understanding of infectious agents over a longer time scale and confrontation with the dynamic of modern TB strains. Nevertheless, we have to consider tuberculosis as a multifactorial disorder in which environmental factors interact tightly with human and pathogen genetic.

A first insight on the population structure of Mycobacterium tuberculosis complex as studied by spoligotyping and MIRU-VNTRs in Santiago, Chile

Tuberculosis (TB) remains a significant public health problem worldwide, but the ecology of the prevalent mycobacterial strains, and their transmission, can vary depending on country and region. Chile is a country with low incidence of TB, that has a geographically isolated location in relation to the rest of South American countries due to the Andes Mountains, but recent migration from neighboring countries has changed this situation. We aimed to assess the genotypic diversity of Mycobacterium tuberculosis complex (MTBC) strains in Santiago, Chile, and compare with reports from other Latin-American countries. We analyzed MTBC isolates from pulmonary tuberculosis cases collected between years 2008 and 2013 in Central Santiago, using two genotyping methods: spoligotyping and 12-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTRs). Data obtained were analyzed and compared to the SITVIT2 database. Mean age of the patients was 47.5 years and 61% were male; 11.6% were migrants. Of 103 strains (1 isolate/patient) included, there were 56 distinct spoligotype patterns. Of these, 16 strains (15.5%) corresponded to orphan strains in the SITVIT2 database, not previously reported. Latin American and Mediterranean (LAM) (34%) and T (33%) lineages were the most prevalent strains, followed by Haarlem lineage (16.5%). Beijing family was scarcely represented with only two cases (1.9%), one of them isolated from a Peruvian migrant. The most frequent clustered spoligotypes were SIT33/LAM3 (10.7%), SIT53/T1 (8.7%), SIT50/H3 (7.8%), and SIT37/T3 (6.8%). We conclude that LAM and T genotypes are the most prevalent genotypes of MTBC in Santiago, Chile, and together correspond to almost two thirds of analyzed strains, which is similar to strain distribution reported from other countries of Latin America. Nevertheless, the high proportion of SIT37/T3, which was rarely found in other Latin American countries, may underline a specific history or demographics of Chile related to probable human migrations and evolutions.

Reemergence and amplification of tuberculosis in the Canadian arctic

BACKGROUND

Between November 2011 and November 2012, a Canadian village of 933 persons had 50 culture-positive cases of tuberculosis, with 49 sharing the same genotype.

METHODS

We performed Illumina-based whole-genome sequencing on Mycobacterium tuberculosis isolates from this village, during and before the outbreak. Phylogenetic trees were generated using the maximum likelihood method.

RESULTS

Three distinct genotypes were identified. Strain I (n = 7) was isolated in 1991-1996. Strain II (n = 8) was isolated in 1996-2004. Strain III (n = 62) first appeared in 2007 and did not arise from strain I or II. Within strain III, there were 3 related but distinct clusters: IIIA, IIIB, and IIIC. Between 2007 and 2010, cluster IIIA predominated (11 of 22 vs 2 of 40; P < .001), whereas in 2011-2012 clusters IIIB (n = 18) and IIIC (n = 20) predominated over cluster IIIA (n = 11). Combined evolutionary and epidemiologic analysis of strain III cases revealed that the outbreak in 2011-2012 was the result of ≥6 temporally staggered events, spanning from 1 reactivation case to a point-source outbreak of 20 cases.

CONCLUSIONS

After the disappearance of 2 strains of M. tuberculosis in this village, its reemergence in 2007 was followed by an epidemiologic amplification, affecting >5% of the population.

A phylogeny-based sampling strategy and power calculator informs genome-wide associations study design for microbial pathogens

Whole genome sequencing is increasingly used to study phenotypic variation among infectious pathogens and to evaluate their relative transmissibility, virulence, and immunogenicity. To date, relatively little has been published on how and how many pathogen strains should be selected for studies associating phenotype and genotype. There are specific challenges when identifying genetic associations in bacteria which often comprise highly structured populations. Here we consider general methodological questions related to sampling and analysis focusing on clonal to moderately recombining pathogens. We propose that a matched sampling scheme constitutes an efficient study design, and provide a power calculator based on phylogenetic convergence. We demonstrate this approach by applying it to genomic datasets for two microbial pathogens: Mycobacterium tuberculosis and Campylobacter species.

Evolution of extensively drug-resistant Mycobacterium tuberculosis from a susceptible ancestor in a single patient

Background

Mycobacterium tuberculosis is characterized by a low mutation rate and a lack of genetic recombination. Yet, the rise of extensively resistant strains paints a picture of a microbe with an impressive adaptive potential. Here we describe the first documented case of extensively drug-resistant tuberculosis evolved from a susceptible ancestor within a single patient.

Results

Genome sequences of nine serial M. tuberculosis isolates from the same patient uncovered a dramatic turnover of competing lineages driven by the emergence, and subsequent fixation or loss of single nucleotide polymorphisms. For most drugs, resistance arose through independent emergence of mutations in more than one clone, of which only one ultimately prevailed as the clone carrying it expanded, displacing the other clones in the process. The vast majority of mutations identified over 3.5 years were either involved in drug resistance or hitchhiking in the genetic background of these. Additionally, RNA-sequencing of isolates grown in the absence of drug challenge revealed that the efflux-associated iniBAC operon was up-regulated over time, whereas down-regulated genes include those involved in mycolic acid synthesis.

Conclusions

We observed both rapid acquisitions of resistance to antimicrobial compounds mediated by individual mutations as well as a gradual increase in fitness in the presence of antibiotics, likely driven by stable gene expression reprogramming. The rapid turnover of resistance mutations and hitchhiking neutral mutations has major implications for inferring tuberculosis transmission events in situations where drug resistance evolves within transmission chains.

Electronic supplementary material

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

Lipidomics and genomics of Mycobacterium tuberculosis reveal lineage-specific trends in mycolic acid biosynthesis

Mycolic acids (MAs) are α-alkyl, β-hydroxy long-chain fatty acids found in abundance in the cell envelope of the Mycobacterium tuberculosis complex (MTBC). MAs form an efficient permeability barrier, modulate host innate immune responses, and are the targets of several anti-tuberculosis drugs. Using mass spectrometry, we measured the relative abundance of 80 MA species across 36 clinical isolates of MTBC covering four major phylogenetic lineages. We found significant variations in the MA patterns between different MTBC strains and lineages. MA patterns of “ancient” lineages contrasted those from “modern” lineages, with a lower representation of alpha-mycolates among Lineage 6 strains and an inversion of the methoxy: keto-mycolates ratio in Lineage 1 strains. By interrogating the whole genome sequences of these MTBC strains, we identified relevant single-nucleotide polymorphisms that may sustain the lineage-specific MA patterns. Our results show that the strain genetic background influences MA metabolism and suggests that strain diversity should be considered in the development of new anti-tuberculosis drugs that target MA synthesis.

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.

A robust SNP barcode for typing Mycobacterium tuberculosis complex strains

Strain-specific genomic diversity in the Mycobacterium tuberculosis complex (MTBC) is an important factor in pathogenesis that may affect virulence, transmissibility, host response and emergence of drug resistance. Several systems have been proposed to classify MTBC strains into distinct lineages and families. Here, we investigate single-nucleotide polymorphisms (SNPs) as robust (stable) markers of genetic variation for phylogenetic analysis. We identify ~92k SNP across a global collection of 1,601 genomes. The SNP-based phylogeny is consistent with the gold-standard regions of difference (RD) classification system. Of the ~7k strain-specific SNPs identified, 62 markers are proposed to discriminate known circulating strains. This SNP-based barcode is the first to cover all main lineages, and classifies a greater number of sublineages than current alternatives. It may be used to classify clinical isolates to evaluate tools to control the disease, including therapeutics and vaccines whose effectiveness may vary by strain type.

Genetic variation in Mycobacterium tuberculosis complex (MTBC) bacteria is responsible for differences in factors such as virulence and transmissibility. Here, the authors analyse the genomes of 1,601 MTBC isolates from diverse geographic locations and identify 62 SNPs that may be used to resolve lineages and sublineages of these strains.

Whole-genome sequencing to control antimicrobial resistance

Following recent improvements in sequencing technologies, whole-genome sequencing (WGS) is positioned to become an essential tool in the control of antibiotic resistance, a major threat in modern healthcare. WGS has already found numerous applications in this area, ranging from the development of novel antibiotics and diagnostic tests through to antibiotic stewardship of currently available drugs via surveillance and the elucidation of the factors that allow the emergence and persistence of resistance. Numerous proof-of-principle studies have also highlighted the value of WGS as a tool for day-to-day infection control and, for some pathogens, as a primary diagnostic tool to detect antibiotic resistance. However, appropriate data analysis platforms will need to be developed before routine WGS can be introduced on a large scale.

TB or not TB: that is no longer the question

Tuberculosis (TB) remains a devastating infectious disease and, with the emergence of multidrug-resistant forms, represents a major global threat. Much of our understanding of pathogenic and immunologic mechanisms in TB has derived from studies in experimental animals. However, it is becoming increasingly clear in TB as well as in other inflammatory diseases that there are substantial differences in immunological responses of humans not found or predicted by animal studies. Thus, it is critically important to understand mechanisms of pathogenesis and immunological protection in humans. In this review, we will address the key immunological question: What are the necessary and sufficient immune responses required for protection against TB infection and disease in people-specifically protection against infection, protection against the establishment of latency or persistence, and protection against transitioning from latent infection to active disease.

Disclosure of selective advantages in the "modern" sublineage of the Mycobacterium tuberculosis Beijing genotype family by quantitative proteomics

The Mycobacterium tuberculosis Beijing genotype, consisting of the more ancient (atypical) and modern (typical) emerging sublineage, is one of the most prevalent and genetically conserved genotype families and has often been associated with multidrug resistance. In this study, we employed a 2D-LC-FTICR MS approach, combined with dimethylation of tryptic peptides, to systematically compare protein abundance levels of ancient and modern Beijing strains and identify differences that could be associated with successful spread of the modern sublineage. The data is available via ProteomeXchange using the identifier PXD000931. Despite the highly uniform protein abundance ratios in both sublineages, we identified four proteins as differentially regulated between both sublineages, which could explain the apparent increased adaptation of the modern Beijing strains. These proteins are; Rv0450c/MmpL4, Rv1269c, Rv3137, and Rv3283/sseA. Transcriptional and functional analysis of these proteins in a large cohort of 29 Beijing strains showed that the mRNA levels of Rv0450c/MmpL4 are significantly higher in modern Beijing strains, whereas we also provide evidence that Rv3283/sseA is less abundant in the modern Beijing sublineage. Our findings provide a possible explanation for the increased virulence and success of the modern Beijing sublineage. In addition, in the established dataset of 1817 proteins, we demonstrate the pre-existence of several, possibly unique, antibiotic efflux pumps in the proteome of the Beijing strains. This may reflect an increased ability of Beijing strains to escape exposure to antituberculosis drugs.

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.