Whole genome sequencing of Mycobacterium tuberculosis reveals slow growth and low mutation rates during latent infections in humans

Mycobacterium tuberculosis is less likely to acquire pathogenic mutations during latent infection than during active disease

Most people infected with Mycobacterium tuberculosis (Mtb) are believed to be in a state of latent tuberculosis (TB) infection (LTBI). Although LTBI is asymptomatic and not infectious, there is a risk of developing active disease even decades after infection. Here, to characterize mutations acquired during LTBI, we collected and analyzed Mtb genomes from seven Japanese patient pairs, each pair consisting of two active TB patients whose starting dates of developing active disease were >3 years apart; one had a high suspicion of LTBI before developing active disease, whereas the other did not. Thereafter, we compared these genomes with those of longitudinal sample pairs within a host of chronic active TB infections combined with public data. The bacterial populations in patients with LTBI were genetically more homogeneous and accumulated single nucleotide polymorphisms (SNPs) slower than those from active disease. Moreover, the lower proportion of nonsynonymous SNPs indicated weaker selective pressures during LTBI than active disease. Finally, the different mutation spectrums indicated different mutators between LTBI and active disease. These results suggest that the likelihood of the acquisition of mutations responsible for antibiotic resistance and increased virulence was lower in the Mtb population from LTBI than active disease.IMPORTANCEControlling latent tuberculosis (TB) infection (LTBI) activation is an effective strategy for TB elimination, where understanding Mycobacterium tuberculosis (Mtb) dynamics within the host plays an important role. Previous studies on chronic active disease reported that Mtb accumulated genomic mutations within the host, possibly resulting in acquired drug resistance and increased virulence. However, several reports suggest that fewer mutations accumulate during LTBI than during the active disease, but the associated risk is largely unknown. Here, we analyzed the genomic dynamics of Mtb within the host during LTBI. Our results statistically suggest that Mtb accumulates mutations during LTBI, but most mutations are under low selective pressures, which induce mutations responsible for drug resistance and virulence. Thus, we propose that LTBI acts as a source for new TB disease rather than as a period for in-host genome evolution.

Unraveling tuberculosis patient cluster transmission chains: integrating WGS-based network with clinical and epidemiological insights

Background

Tuberculosis remains a global health threat, and the World Health Organization reports a limited reduction in disease incidence rates, including both new and relapse cases. Therefore, studies targeting tuberculosis transmission chains and recurrent episodes are crucial for developing the most effective control measures. Herein, multiple tuberculosis clusters were retrospectively investigated by integrating patients' epidemiological and clinical information with median-joining networks recreated based on whole genome sequencing (WGS) data of Mycobacterium tuberculosis isolates.

Methods

Epidemiologically linked tuberculosis patient clusters were identified during the source case investigation for pediatric tuberculosis patients. Only M. tuberculosis isolate DNA samples with previously determined spoligotypes identical within clusters were subjected to WGS and further median-joining network recreation. Relevant clinical and epidemiological data were obtained from patient medical records.

Results

We investigated 18 clusters comprising 100 active tuberculosis patients 29 of whom were children at the time of diagnosis; nine patients experienced recurrent episodes. M. tuberculosis isolates of studied clusters belonged to Lineages 2 (sub-lineage 2.2.1) and 4 (sub-lineages 4.3.3, 4.1.2.1, 4.8, and 4.2.1), while sub-lineage 4.3.3 (LAM) was the most abundant. Isolates of six clusters were drug-resistant. Within clusters, the maximum genetic distance between closely related isolates was only 5-11 single nucleotide variants (SNVs). Recreated median-joining networks, integrated with patients' diagnoses, specimen collection dates, sputum smear microscopy, and epidemiological investigation results indicated transmission directions within clusters and long periods of latent infection. It also facilitated the identification of potential infection sources for pediatric patients and recurrent active tuberculosis episodes refuting the reactivation possibility despite the small genetic distance of ≤5 SNVs between isolates. However, unidentified active tuberculosis cases within the cluster, the variable mycobacterial mutation rate in dormant and active states, and low M. tuberculosis genetic variability inferred precise transmission chain delineation. In some cases, heterozygous SNVs with an allelic frequency of 10-73% proved valuable in identifying direct transmission events.

Conclusion

The complex approach of integrating tuberculosis cluster WGS-data-based median-joining networks with relevant epidemiological and clinical data proved valuable in delineating epidemiologically linked patient transmission chains and deciphering causes of recurrent tuberculosis episodes within clusters.

Differential Selection for Translation Efficiency Shapes Translation Machineries in Bacterial Species

Different bacterial species have dramatically different generation times, from 20-30 min in Escherichia coli to about two weeks in Mycobacterium leprae. The translation machinery in a cell needs to synthesize all proteins for a new cell in each generation. The three subprocesses of translation, i.e., initiation, elongation, and termination, are expected to be under stronger selection pressure to optimize in short-generation bacteria (SGB) such as Vibrio natriegens than in the long-generation Mycobacterium leprae. The initiation efficiency depends on the start codon decoded by the initiation tRNA, the optimal Shine-Dalgarno (SD) decoded by the anti-SD (aSD) sequence on small subunit rRNA, and the secondary structure that may embed the initiation signals and prevent them from being decoded. The elongation efficiency depends on the tRNA pool and codon usage. The termination efficiency in bacteria depends mainly on the nature of the stop codon and the nucleotide immediately downstream of the stop codon. By contrasting SGB with long-generation bacteria (LGB), we predict (1) SGB to have more ribosome RNA operons to produce ribosomes, and more tRNA genes for carrying amino acids to ribosomes, (2) SGB to have a higher percentage of genes using AUG as the start codon and UAA as the stop codon than LGB, (3) SGB to exhibit better codon and anticodon adaptation than LGB, and (4) SGB to have a weaker secondary structure near the translation initiation signals than LGB. These differences between SGB and LGB should be more pronounced in highly expressed genes than the rest of the genes. We present empirical evidence in support of these predictions.

Parallel signatures of Mycobacterium tuberculosis and human Y-chromosome phylogeography support the Two Layer model of East Asian population history

The Two Layer hypothesis is fast becoming the favoured narrative describing East Asian population history. Under this model, hunter-gatherer groups who initially peopled East Asia via a route south of the Himalayas were assimilated by agriculturalist migrants who arrived via a northern route across Eurasia. A lack of ancient samples from tropical East Asia limits the resolution of this model. We consider insight afforded by patterns of variation within the human pathogen Mycobacterium tuberculosis (Mtb) by analysing its phylogeographic signatures jointly with the human Y-chromosome. We demonstrate the Y-chromosome lineages enriched in the traditionally hunter-gatherer groups associated with East Asia's first layer of peopling to display deep roots, low long-term effective population size, and diversity patterns consistent with a southern entry route. These characteristics mirror those of the evolutionarily ancient Mtb lineage 1. The remaining East Asian Y-chromosome lineage is almost entirely absent from traditionally hunter-gatherer groups and displays spatial and temporal characteristics which are incompatible with a southern entry route, and which link it to the development of agriculture in modern-day China. These characteristics mirror those of the evolutionarily modern Mtb lineage 2. This model paves the way for novel host-pathogen coevolutionary research hypotheses in East Asia.

Hypobiosis of Mycobacteria: Biochemical Aspects

Under suboptimal growth conditions, bacteria can transit to the dormant forms characterized by a significantly reduced metabolic activity, resistance to various stress factors, and absence of cell proliferation. Traditionally, the dormant state is associated with the formation of highly differentiated cysts and spores. However, non-spore-forming bacteria can transfer to the dormant-like hypobiotic state with the generation of less differentiated cyst-like forms (which are different from spores). This review focuses on morphological and biochemical changes occurred during formation of dormant forms of mycobacteria in particular pathogenic M. tuberculosis (Mtb) caused latent forms of tuberculosis. These forms are characterized by the low metabolic activity, the absence of cell division, resistance to some antibiotics, marked morphological changes, and loss of ability to grow on standard solid media ("non-culturable" state). Being produced in vitro, dormant Mtb retained ability to maintain latent infection in mice. After a long period of dormancy, mycobacteria retain a number of stable proteins with a potential enzymatic activity which could participate in maintaining of low-level metabolic activity in period of dormancy. Indeed, the metabolomic analysis showed significant levels of metabolites in the dormant cells even after a long period of dormancy, which may be indicative of residual metabolism in dormant mycobacteria. Special role may play intracellularly accumulated trehalose in dormant mycobacteria. Trehalose appears to stabilize dormant cells, as evidenced by the direct correlation between the trehalose content and cell viability during the long-term dormancy. In addition, trehalose can be considered as a reserve energy substrate consumed during reactivation of dormant mycobacteria due to the ATP-dependent conversion of trehalase from the latent to the active state. Another feature of dormant mycobacteria is a high representation of proteins participating in the enzymatic defense against stress factors and of low-molecular-weight compounds protecting cells in the absence of replication. Dormant mycobacteria contain a large number of hydrolyzing enzymes, which, on the one hand, ensure inactivation of biomolecules damaged by stress. On the other hand, the products of these enzymatic reactions can be used for the maintenance of energy state and vital activity of bacterial cells during their long-term survival in the dormant state, i.e., for creating a situation that we propose to refer to as the "catabolic survival". In general, dormant non-replicating mycobacterial cells can be described as morphologically altered forms that contain principal macromolecules and are stabilized and protected from the damaging factors by an arsenal of proteins and low-molecular-weight compounds. Because of the presumable occurrence of metabolic reactions in such cells, this form of survival should be referred to as hypobiosis.

Tuberkulose im Erwachsenenalter

Die Tuberkulose ist in Deutschland eine seltene, überwiegend gut behandelbare Erkrankung. Weltweit ist sie eine der häufigsten Infektionserkrankungen mit ca. 10 Millionen Neuerkrankungen/Jahr. Auch bei einer niedrigen Inzidenz in Deutschland bleibt Tuberkulose insbesondere aufgrund der internationalen Entwicklungen und Migrationsbewegungen eine wichtige Differenzialdiagnose. In Deutschland besteht, aufgrund der niedrigen Prävalenz der Erkrankung und der damit verbundenen abnehmenden klinischen Erfahrung, ein Informationsbedarf zu allen Aspekten der Tuberkulose und ihrer Kontrolle. Diese Leitlinie umfasst die mikrobiologische Diagnostik, die Grundprinzipien der Standardtherapie, die Behandlung verschiedener Organmanifestationen, den Umgang mit typischen unerwünschten Arzneimittelwirkungen, die Besonderheiten in der Diagnostik und Therapie resistenter Tuberkulose sowie die Behandlung bei TB-HIV-Koinfektion. Sie geht darüber hinaus auf Versorgungsaspekte und gesetzliche Regelungen wie auch auf die Diagnosestellung und präventive Therapie einer latenten tuberkulösen Infektion ein. Es wird ausgeführt, wann es der Behandlung durch spezialisierte Zentren bedarf.

Die Aktualisierung der S2k-Leitlinie „Tuberkulose im Erwachsenenalter“ soll allen in der Tuberkuloseversorgung Tätigen als Richtschnur für die Prävention, die Diagnose und die Therapie der Tuberkulose dienen und helfen, den heutigen Herausforderungen im Umgang mit Tuberkulose in Deutschland gewachsen zu sein.

Estimation of the mutation rate of Mycobacterium tuberculosis in cases with recurrent tuberculosis using whole genome sequencing

The study of tuberculosis latency is problematic due to the difficulty of isolating the bacteria in the dormancy state. Despite this, several in vivo approaches have been taken to mimic the latency process. Our group has studied the evolution of the bacteria in 18 cases of recurrent tuberculosis. We found that HIV positive patients develop recurrent tuberculosis earlier, generally in the first two years (p value = 0.041). The genome of the 36 Mycobacterium tuberculosis paired isolates (first and relapsed isolates) showed that none of the SNPs found within each pair was observed more than once, indicating that they were not directly related to the recurrence process. Moreover, some IS6110 movements were found in the paired isolates, indicating the presence of different clones within the patient. Finally, our results suggest that the mutation rate remains constant during all the period as no correlation was found between the number of SNPs and the time to relapse.

The evolving biology of Mycobacterium tuberculosis drug resistance

Tuberculosis, caused by Mycobacterium tuberculosis (Mtb) is an ancient disease that has remained a leading cause of infectious death. Mtb has evolved drug resistance to every antibiotic regimen ever introduced, greatly complicating treatment, lowering rates of cure and menacing TB control in parts of the world. As technology has advanced, our understanding of antimicrobial resistance has improved, and our models of the phenomenon have evolved. In this review, we focus on recent research progress that supports an updated model for the evolution of drug resistance in Mtb. We highlight the contribution of drug tolerance on the path to resistance, and the influence of heterogeneity on tolerance. Resistance is likely to remain an issue for as long as drugs are needed to treat TB. However, with technology driving new insights and careful management of newly developed resources, antimicrobial resistance need not continue to threaten global progress against TB, as it has done for decades.

Feasibility of novel approaches to detect viable Mycobacterium tuberculosis within the spectrum of the tuberculosis disease

Tuberculosis (TB) is a global disease caused by Mycobacterium tuberculosis (Mtb) and is manifested as a continuum spectrum of infectious states. Both, the most common and clinically asymptomatic latent tuberculosis infection (LTBI), and the symptomatic disease, active tuberculosis (TB), are at opposite ends of the spectrum. Such binary classification is insufficient to describe the existing clinical heterogeneity, which includes incipient and subclinical TB. The absence of clinically TB-related symptoms and the extremely low bacterial burden are features shared by LTBI, incipient and subclinical TB states. In addition, diagnosis relies on cytokine release after antigenic T cell stimulation, yet several studies have shown that a high proportion of individuals with immunoreactivity never developed disease, suggesting that they were no longer infected. LTBI is estimated to affect to approximately one fourth of the human population and, according to WHO data, reactivation of LTBI is the main responsible of TB cases in developed countries. Assuming the drawbacks associated to the current diagnostic tests at this part of the disease spectrum, properly assessing individuals at real risk of developing TB is a major need. Further, it would help to efficiently design preventive treatment. This quest would be achievable if information about bacterial viability during human silent Mtb infection could be determined. Here, we have evaluated the feasibility of new approaches to detect viable bacilli across the full spectrum of TB disease. We focused on methods that specifically can measure host-independent parameters relying on the viability of Mtb either by its direct or indirect detection.

Mutation of Mycobacterium tuberculosis and Implications for Using Whole-Genome Sequencing for Investigating Recent Tuberculosis Transmission

Tuberculosis (TB) control programs use whole-genome sequencing (WGS) of Mycobacterium tuberculosis (Mtb) for detecting and investigating TB case clusters. Existence of few genomic differences between Mtb isolates might indicate TB cases are the result of recent transmission. However, the variable and sometimes long duration of latent infection, combined with uncertainty in the Mtb mutation rate during latency, can complicate interpretation of WGS results. To estimate the association between infection duration and single nucleotide polymorphism (SNP) accumulation in the Mtb genome, we first analyzed pairwise SNP differences among TB cases from Los Angeles County, California, with strong epidemiologic links. We found that SNP distance alone was insufficient for concluding that cases are linked through recent transmission. Second, we describe a well-characterized cluster of TB cases in California to illustrate the role of genomic data in conclusions regarding recent transmission. Longer presumed latent periods were inconsistently associated with larger SNP differences. Our analyses suggest that WGS alone cannot be used to definitively determine that a case is attributable to recent transmission. Methods for integrating clinical, epidemiologic, and genomic data can guide conclusions regarding the likelihood of recent transmission, providing local public health practitioners with better tools for monitoring and investigating TB transmission.

Treatment in Latent Tuberculosis Uveitis-Is Immunosuppression Effective or Is Conventional 3- or 4-Drug Antituberculosis Therapy Mandatory?

Background/Aims: Controversy exists regarding 3- or 4 drug antituberculosis therapy (conventional ATT) in uveitis patients having latent tuberculosis (LTB), especially while initiating therapy with corticosteroids and/or other immunosuppressants. Methods: We performed a monocentral retrospective analysis of posterior uveitis patients with latent TB. Latent TB was diagnosed, in case of a positive QuantiFERON®-TB-Gold test and normal chest imaging, after ruling out other causes of infectious and noninfectious uveitis. Patients with active TB were excluded. From 2016 to 2020 we included 17 patients. Ophthalmological evaluation consisted of Best corrected visual acuity (BCVA), slit lamp examination, fundoscopy, OCT, and fluorescein- and indocyaningreen- angiography before and at months 3, 6, 12, 24, and the last follow-up after treatment. Results: Initially, all patients had active posterior uveitis with occlusive (n = 5 patients) and nonocclusive retinal vasculitis (n = 12 patients). Mean follow up was 28 ± 15 months. Therapy was started with systemic corticosteroids (mean prednisolone equivalent 71.3 mg/d) and already after 3 months it could be tapered to a mean maintenance dosage of 8.63 mg/d. Additional immunosuppressive treatment with cs- or bDMARDs was initiated in 14 patients (82%) due to recurrences of uveitis while tapering the corticosteroids <10 mg per/day or because of severe inflammation at the initial visit. While being on immunosuppression, best corrected visual acuity increased from 0.56 logMAR to 0.32 logMAR during follow-up and only three patients had one uveitis relapse, which was followed by switch of immunosuppressive treatment. As recommended, TB prophylaxis with 300 mg/d isoniazid was administered in 11 patients for at least 9 months while being on TNF-alpha-blocking agents. No patient developed active tuberculosis during immunosuppressive therapy. Conclusion: Mainly conventional ATT is strongly recommended—as monotherapy or in combination with immunosuppressives—for effective treatment in patients with uveitis due to latent TB. Although in our patient group no conventional ATT was initiated, immunosuppression alone occurred as an efficient treatment. Nevertheless, due to possible activation of TB, isoniazid prophylaxis is mandatory in latent TB patients while being on TNF-alpha blocking agents.

The Biological and Clinical Aspects of a Latent Tuberculosis Infection

Tuberculosis (TB), caused by bacilli from the Mycobacterium tuberculosis complex, remains a serious global public health problem, representing one of the main causes of death from infectious diseases. About one quarter of the world’s population is infected with Mtb and has a latent TB infection (LTBI). According to the World Health Organization (WHO), an LTBI is characterized by a lasting immune response to Mtb antigens without any TB symptoms. Current LTBI diagnoses and treatments are based on this simplified definition, although an LTBI involves a broad range of conditions, including when Mtb remains in the body in a persistent form and the immune response cannot be detected. The study of LTBIs has progressed in recent years; however, many biological and medical aspects of an LTBI are still under discussion. This review focuses on an LTBI as a broad spectrum of states, both of the human body, and of Mtb cells. The problems of phenotypic insusceptibility, diagnoses, chemoprophylaxis, and the necessity of treatment are discussed. We emphasize the complexity of an LTBI diagnosis and its treatment due to its ambiguous nature. We consider alternative ways of differentiating an LTBI from active TB, as well as predicting TB reactivation based on using mycobacterial “latency antigens” for interferon gamma release assay (IGRA) tests and the transcriptomic analysis of human blood cells.

Targeting Non-Replicating Mycobacterium tuberculosis and Latent Infection: Alternatives and Perspectives (Mini-Review)

Latent tuberculosis infection (LTBI) represents a major challenge to curing TB disease. Current guidelines for LTBI management include only three older drugs and their combinations-isoniazid and rifamycins (rifampicin and rifapentine). These available control strategies have little impact on latent TB elimination, and new specific therapeutics are urgently needed. In the present mini-review, we highlight some of the alternatives that may potentially be included in LTBI treatment recommendations and a list of early-stage prospective small molecules that act on drug targets specific for Mycobacterium tuberculosis latency.

The within-host evolution of antimicrobial resistance in Mycobacterium tuberculosis

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

Accounting for the Biological Complexity of Pathogenic Fungi in Phylogenetic Dating

In the study of pathogen evolution, temporal dating of phylogenies provides information on when species and lineages may have diverged in the past. When combined with spatial and epidemiological data in phylodynamic models, these dated phylogenies can also help infer where and when outbreaks occurred, how pathogens may have spread to new geographic locations and/or niches, and how virulence or drug resistance has developed over time. Although widely applied to viruses and, increasingly, to bacterial pathogen outbreaks, phylogenetic dating is yet to be widely used in the study of pathogenic fungi. Fungi are complex organisms with several biological processes that could present issues with appropriate inference of phylogenies, clock rates, and divergence times, including high levels of recombination and slower mutation rates although with potentially high levels of mutation rate variation. Here, we discuss some of the key methodological challenges in accurate phylogeny reconstruction for fungi in the context of the temporal analyses conducted to date and make recommendations for future dating studies to aid development of a best practices roadmap in light of the increasing threat of fungal outbreaks and antifungal drug resistance worldwide.

Identifying likely transmissions in Mycobacterium bovis infected populations of cattle and badgers using the Kolmogorov Forward Equations

Established methods for whole-genome-sequencing (WGS) technology allow for the detection of single-nucleotide polymorphisms (SNPs) in the pathogen genomes sourced from host samples. The information obtained can be used to track the pathogen’s evolution in time and potentially identify ‘who-infected-whom’ with unprecedented accuracy. Successful methods include ‘phylodynamic approaches’ that integrate evolutionary and epidemiological data. However, they are typically computationally intensive, require extensive data, and are best applied when there is a strong molecular clock signal and substantial pathogen diversity. To determine how much transmission information can be inferred when pathogen genetic diversity is low and metadata limited, we propose an analytical approach that combines pathogen WGS data and sampling times from infected hosts. It accounts for ‘between-scale’ processes, in particular within-host pathogen evolution and between-host transmission. We applied this to a well-characterised population with an endemic Mycobacterium bovis (the causative agent of bovine/zoonotic tuberculosis, bTB) infection. Our results show that, even with such limited data and low diversity, the computation of the transmission probability between host pairs can help discriminate between likely and unlikely infection pathways and therefore help to identify potential transmission networks. However, the method can be sensitive to assumptions about within-host evolution.

Higher genome mutation rates of Beijing lineage of Mycobacterium tuberculosis during human infection

Mycobacterium tuberculosis (Mtb) strains of Beijing lineage have caused great concern because of their rapid emergence of drug resistance and worldwide spread. DNA mutation rates that reflect evolutional adaptation to host responses and the appearance of drug resistance have not been elucidated in human-infected Beijing strains. We tracked and obtained an original Mtb isolate of Beijing lineage from the 1999 tuberculosis outbreak in Japan, as well as five other isolates that spread in humans, and two isolates from the patient caused recurrence. Three isolates were from patients who developed TB within one year after infection (rapid-progressor, RP), and the other three isolates were from those who developed TB more than one year after infection (slow-progressor, SP). We sequenced genomes of these isolates and analyzed the propensity and rate of genomic mutations. Generation time versus mutation rate curves were significantly higher for RP. The ratio of oxidative versus non-oxidation damages induced mutations was higher in SP than RP, suggesting that persistent Mtb are exposed to oxidative stress in the latent state. Our data thus demonstrates that higher mutation rates of Mtb Beijing strains during human infection is likely to account for the higher adaptability and an emergence ratio of drug resistance.

Mycobacterium tuberculosis progresses through two phases of latent infection in humans

Little is known about the physiology of latent Mycobacterium tuberculosis infection. We studied the mutational rates of 24 index tuberculosis (TB) cases and their latently infected household contacts who developed active TB up to 5.25 years later, as an indication of bacterial physiological state and possible generation times during latent TB infection in humans. Here we report that the rate of new mutations in the M. tuberculosis genome decline dramatically after two years of latent infection (two-sided p < 0.001, assuming an 18 h generation time equal to log phase M. tuberculosis, with latency period modeled as a continuous variable). Alternatively, assuming a fixed mutation rate, the generation time increases over the latency duration. Mutations indicative of oxidative stress do not increase with increasing latency duration suggesting a lack of host or bacterial derived mutational stress. These results suggest that M. tuberculosis enters a quiescent state during latency, decreasing the risk for mutational drug resistance and increasing generation time, but potentially increasing bacterial tolerance to drugs that target actively growing bacteria.

Evaluating the contributions of purifying selection and progeny-skew in dictating within-host Mycobacterium tuberculosis evolution

The within-host evolutionary dynamics of tuberculosis (TB) remain unclear, and underlying biological characteristics render standard population genetic approaches based upon the Wright-Fisher model largely inappropriate. In addition, the compact genome combined with an absence of recombination is expected to result in strong purifying selection effects. Thus, it is imperative to establish a biologically relevant evolutionary framework incorporating these factors in order to enable an accurate study of this important human pathogen. Further, such a model is critical for inferring fundamental evolutionary parameters related to patient treatment, including mutation rates and the severity of infection bottlenecks. We here implement such a model and infer the underlying evolutionary parameters governing within-patient evolutionary dynamics. Results demonstrate that the progeny skew associated with the clonal nature of TB severely reduces genetic diversity and that the neglect of this parameter in previous studies has led to significant mis-inference of mutation rates. As such, our results suggest an underlying de novo mutation rate that is considerably faster than previously inferred, and a progeny distribution differing significantly from Wright-Fisher assumptions. This inference represents a more appropriate evolutionary null model, against which the periodic effects of positive selection, associated with drug-resistance for example, may be better assessed.

Application of whole-genome sequencing in a case study of renal tuberculosis in a child

BACKGROUND

Urogenital tuberculosis (TB) is rare in children and usually develops due to reactivation of the foci in the genitourinary tract after the latency period following initial infection. Urogenital TB in children has no pathognomonic clinical features that can result in overlooking or misdiagnosing this clinical entity. Here, we report important findings regarding the pathogenesis and transmission of TB by using genotyping and whole-genome sequencing (WGS) in a study of renal TB case in a child.

CASE PRESENTATION

A 13-year-old boy was admitted to the hospital because of high fever, severe dry cough, flank pain and painful urination. Abdominal ultrasonography and CT revealed an 8 mm calculus in the kidney, and clinical findings were initially interpreted as nephrolithiasis. Nevertheless, due to the atypical clinical presentation of kidney stone disease, additional investigations for possible TB were performed. The QuantiFERON®-TB Gold Plus test was positive, and the Mantoux test resulted in 15 mm of induration, confirming infection with Mycobacterium tuberculosis (Mtb). Chest X-ray was normal. Chest CT revealed calcified intrathoracic lymph nodes. The urine sample tested positive for acid-fast bacilli, and Mtb cultures were obtained from urine and bronchial aspirate samples, resulting in a final diagnosis of intrathoracic lymph node and renal TB. Contact investigation revealed that the child's father was diagnosed with TB when the child was 1 year old. Genotyping and WGS analysis of Mtb isolates of the child and his father confirmed the epidemiological link and pointed to the latency of infection in the child.

CONCLUSIONS

This case report confirmed the development of active TB from calcified lesions in adolescent after 12 years of exposure, demonstrated the absence of microevolutionary changes in the Mtb genome during the period of latency, and proved the importance of appropriate evaluation and management to prevent the progression of TB infection to active TB disease. The use of WGS provided the ultimate resolution for the detection of TB transmission and reactivation events.

Comparison of routine field epidemiology and whole genome sequencing to identify tuberculosis transmission in a remote setting

Yukon Territory (YT) is a remote region in northern Canada with ongoing spread of tuberculosis (TB). To explore the utility of whole genome sequencing (WGS) for TB surveillance and monitoring in a setting with detailed contact tracing and interview data, we used a mixed-methods approach. Our analysis included all culture-confirmed cases in YT (2005-2014) and incorporated data from 24-locus Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) genotyping, WGS and contact tracing. We compared field-based (contact investigation (CI) data + MIRU-VNTR) and genomic-based (WGS + MIRU-VNTR + basic case data) investigations to identify the most likely source of each person's TB and assessed the knowledge, attitudes and practices of programme personnel around genotyping and genomics using online, multiple-choice surveys (n = 4) and an in-person group interview (n = 5). Field- and genomics-based approaches agreed for 26 of 32 (81%) cases on likely location of TB acquisition. There was less agreement in the identification of specific source cases (13/22 or 59% of cases). Single-locus MIRU-VNTR variants and limited genetic diversity complicated the analysis. Qualitative data indicated that participants viewed genomic epidemiology as a useful tool to streamline investigations, particularly in differentiating latent TB reactivation from the recent transmission. Based on this, genomic data could be used to enhance CIs, focus resources, target interventions and aid in TB programme evaluation.

Focusing on DNA Repair and Damage Tolerance Mechanisms in Mycobacterium tuberculosis: An Emerging Therapeutic Theme

Tuberculosis (TB) is one such disease that has become a nuisance in the world scenario and one of the most deadly diseases of the current times. The etiological agent of tuberculosis, Mycobacterium tuberculosis (M. tb) kills millions of people each year. Not only 1.7 million people worldwide are estimated to harbor M. tb in the latent form but also 5 to 15 percent of which are expected to acquire an infection during a lifetime. Though curable, a long duration of drug regimen and expense leads to low patient adherence. The emergence of multi-, extensive- and total- drug-resistant strains of M. tb further complicates the situation. Owing to high TB burden, scientists worldwide are trying to design novel therapeutics to combat this disease. Therefore, to identify new drug targets, there is a growing interest in targeting DNA repair pathways to fight this infection. Thus, this review aims to explore DNA repair and damage tolerance as an efficient target for drug development by understanding M. tb DNA repair and tolerance machinery and its regulation, its role in pathogenesis and survival, mutagenesis, and consequently, in the development of drug resistance.

Dispersal of Mycobacterium tuberculosis Driven by Historical European Trade in the South Pacific

Mycobacterium tuberculosis (Mtb) is a globally distributed bacterial pathogen whose population structure has largely been shaped by the activities of its obligate human host. Oceania was the last major global region to be reached by Europeans and is the last region for which the dispersal and evolution of Mtb remains largely unexplored. Here, we investigated the evolutionary history of the Euro-American L4.4 sublineage and its dispersal to the South Pacific. Using a phylodynamics approach and a dataset of 236 global Mtb L4.4 genomes we have traced the origins and dispersal of L4.4 strains to New Zealand. These strains are predominantly found in indigenous Māori and Pacific people and we identify a clade of European, likely French, origin that is prevalent in indigenous populations in both New Zealand and Canada. Molecular dating suggests the expansion of European trade networks in the early 19th century drove the dispersal of this clade to the South Pacific. We also identify historical and social factors within the region that have contributed to the local spread and expansion of these strains, including recent Pacific migrations to New Zealand and the rapid urbanization of Māori in the 20th century. Our results offer new insight into the expansion and dispersal of Mtb in the South Pacific and provide a striking example of the role of historical European migrations in the global dispersal of Mtb.

Investigating Evolutionary Rate Variation in Bacteria

Rates of molecular evolution are known to vary between species and across all kingdoms of life. Here, we explore variation in the rate at which bacteria accumulate mutations (accumulation rates) in their natural environments over short periods of time. We have compiled estimates of the accumulation rate for over 34 species of bacteria, the majority of which are pathogens evolving either within an individual host or during outbreaks. Across species, we find that accumulation rates vary by over 3700-fold. We investigate whether accumulation rates are associated to a number potential correlates including genome size, GC content, measures of the natural selection and the time frame over which the accumulation rates were estimated. After controlling for phylogenetic non-independence, we find that the accumulation rate is not significantly correlated to any factor. Furthermore, contrary to previous results, we find that it is not impacted by the time frame of which the estimate was made. However, our study, with only 34 species, is likely to lack power to detect anything but large effects. We suggest that much of the rate variation may be explained by differences between species in the generation time in the wild.

Beyond the SNP Threshold: Identifying Outbreak Clusters Using Inferred Transmissions

Whole-genome sequencing (WGS) is increasingly used to aid the understanding of pathogen transmission. A first step in analyzing WGS data is usually to define "transmission clusters," sets of cases that are potentially linked by direct transmission. This is often done by including two cases in the same cluster if they are separated by fewer single-nucleotide polymorphisms (SNPs) than a specified threshold. However, there is little agreement as to what an appropriate threshold should be. We propose a probabilistic alternative, suggesting that the key inferential target for transmission clusters is the number of transmissions separating cases. We characterize this by combining the number of SNP differences and the length of time over which those differences have accumulated, using information about case timing, molecular clock, and transmission processes. Our framework has the advantage of allowing for variable mutation rates across the genome and can incorporate other epidemiological data. We use two tuberculosis studies to illustrate the impact of our approach: with British Columbia data by using spatial divisions; with Republic of Moldova data by incorporating antibiotic resistance. Simulation results indicate that our transmission-based method is better in identifying direct transmissions than a SNP threshold, with dissimilarity between clusterings of on average 0.27 bits compared with 0.37 bits for the SNP-threshold method and 0.84 bits for randomly permuted data. These results show that it is likely to outperform the SNP-threshold method where clock rates are variable and sample collection times are spread out. We implement the method in the R package transcluster.

DNA Replication Fidelity in the Mycobacterium tuberculosis Complex

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

DNA Replication in Mycobacterium tuberculosis

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

Differential carriage of virulence-associated loci in the New Zealand Rangipo outbreak strain of Mycobacterium tuberculosis

BACKGROUND

The Rangipo strain of Mycobacterium tuberculosis achieved notoriety in New Zealand due to its role in several tuberculosis (TB) outbreaks. Why this strain should be the source of relatively large clusters of the disease is unknown. In this work, we performed an in-depth analysis of the genome of the Rangipo strain to determine whether it offers clues to understanding its prevalence.

METHODS

Next-generation sequencing was performed on nine isolates which matched the Rangipo genotypic profile. Sequence reads were assembled against the H37Rv reference genome and single-locus variants identified. Unmapped reads were compared against the genome sequences of other M. tuberculosis strains, in particular CDC1551, Haarlem and Erdman.

RESULTS

Across the nine Rangipo strains, a total of 727 single-locus variants were identified with respect to H37Rv, of which 700 were common to all Rangipo strains sequenced. Within the common variants, 386 were non-synonymous, with 12 occurring in genes associated with M. tuberculosis virulence. Next-generation and Sanger sequencing determined the presence of three genes in the Rangipo isolates, which are absent in H37Rv, but which have been reported to be important for the pathogenicity of M. tuberculosis. The differentially encoded Rangipo genes consisted of transcriptional regulator EmbR2, and molybdopterin cofactor biosynthesis proteins A and B. The Rangipo strain also harbours an extended DNA helicase and an additional adenylate cyclase.

CONCLUSIONS

Our study provides new insights into the genomic content of the New Zealand Rangipo strain of M. tuberculosis and highlights the presence of additional virulence-related loci not found in H37Rv.

Using whole genome sequencing to investigate transmission in a multi-host system: bovine tuberculosis in New Zealand

Background

Bovine tuberculosis (bTB), caused by Mycobacterium bovis, is an important livestock disease raising public health and economic concerns around the world. In New Zealand, a number of wildlife species are implicated in the spread and persistence of bTB in cattle populations, most notably the brushtail possum (Trichosurus vulpecula). Whole Genome Sequenced (WGS) M. bovis isolates sourced from infected cattle and wildlife across New Zealand were analysed. Bayesian phylogenetic analyses were conducted to estimate the substitution rate of the sampled population and investigate the role of wildlife. In addition, the utility of WGS was examined with a view to these methods being incorporated into routine bTB surveillance.

Results

A high rate of exchange was evident between the sampled wildlife and cattle populations but directional estimates of inter-species transmission were sensitive to the sampling strategy employed. A relatively high substitution rate was estimated, this, in combination with a strong spatial signature and a good agreement to previous typing methods, acts to endorse WGS as a typing tool.

Conclusions

In agreement with the current knowledge of bTB in New Zealand, transmission of M. bovis between cattle and wildlife was evident. Without direction, these estimates are less informative but taken in conjunction with the low prevalence of bTB in New Zealand’s cattle population it is likely that, currently, wildlife populations are acting as the main bTB reservoir. Wildlife should therefore continue to be targeted if bTB is to be eradicated from New Zealand. WGS will be a considerable aid to bTB eradication by greatly improving the discriminatory power of molecular typing data. The substitution rates estimated here will be an important part of epidemiological investigations using WGS data.

Electronic supplementary material

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

The Evolution of Strain Typing in the Mycobacterium tuberculosis Complex

Tuberculosis (TB) is a contagious disease with a complex epidemiology. Therefore, molecular typing (genotyping) of Mycobacterium tuberculosis complex (MTBC) strains is of primary importance to effectively guide outbreak investigations, define transmission dynamics and assist global epidemiological surveillance of the disease. Large-scale genotyping is also needed to get better insights into the biological diversity and the evolution of the pathogen. Thanks to its shorter turnaround and simple numerical nomenclature system, mycobacterial interspersed repetitive unit-variable-number tandem repeat (MIRU-VNTR) typing, based on 24 standardized plus 4 hypervariable loci, optionally combined with spoligotyping, has replaced IS6110 DNA fingerprinting over the last decade as a gold standard among classical strain typing methods for many applications. With the continuous progress and decreasing costs of next-generation sequencing (NGS) technologies, typing based on whole genome sequencing (WGS) is now increasingly performed for near complete exploitation of the available genetic information. However, some important challenges remain such as the lack of standardization of WGS analysis pipelines, the need of databases for sharing WGS data at a global level, and a better understanding of the relevant genomic distances for defining clusters of recent TB transmission in different epidemiological contexts. This chapter provides an overview of the evolution of genotyping methods over the last three decades, which culminated with the development of WGS-based methods. It addresses the relative advantages and limitations of these techniques, indicates current challenges and potential directions for facilitating standardization of WGS-based typing, and provides suggestions on what method to use depending on the specific research question.

The tuberculocidal activity of polyaniline and functionalised polyanilines

Tuberculosis is considered a leading cause of death worldwide. More than 95% of cases and deaths occur in low- and middle-income countries. In resource-limited countries, hospitals often lack adequate facilities to manage and isolate patients with infectious tuberculosis (TB), relying instead on personal protective equipment, such as facemasks, to reduce nosocomial transmission of the disease. Facemasks impregnated with an antimicrobial agent may be a cost-effective way of adding an extra level of protection against the spread of TB by reducing the risk of disease transmission. Conducting polymers, such as polyaniline (PANI), and their functionalised derivatives are a novel class of antimicrobial agents with potential as non-leaching additives to provide contamination resistant surfaces. We have investigated the antimicrobial action of PANI and a functionalised derivative, poly-3-aminobenzoic acid (P3ABA), against mycobacteria and have determined the optimal treatment time and concentration to achieve significant knockdown of Mycobacterium smegmatis and Mycobacterium tuberculosis on an agar surface. Results indicated that P3ABA is a potential candidate for use as an anti-tuberculoid agent in facemasks to reduce TB transmission.

Genomic Analysis of the Evolution of Fluoroquinolone Resistance in Mycobacterium tuberculosis Prior to Tuberculosis Diagnosis

Fluoroquinolones (FQs) are effective second-line drugs for treating antibiotic-resistant tuberculosis (TB) and are being considered for use as first-line agents. Because FQs are used to treat a range of infections, in a setting of undiagnosed TB, there is potential to select for drug-resistant Mycobacterium tuberculosis mutants during FQ-based treatment of other infections, including pneumonia. Here we present a detailed characterization of ofloxacin-resistant M. tuberculosis samples isolated directly from patients in Taiwan, which demonstrates that selection for FQ resistance can occur within patients who have not received FQs for the treatment of TB. Several of these samples showed no mutations in gyrA or gyrB based on PCR-based molecular assays, but genome-wide next-generation sequencing (NGS) revealed minority populations of gyrA and/or gyrB mutants. In other samples with PCR-detectable gyrA mutations, NGS revealed subpopulations containing alternative resistance-associated genotypes. Isolation of individual clones from these apparently heterogeneous samples confirmed the presence of the minority drug-resistant variants suggested by the NGS data. Further NGS of these purified clones established evolutionary links between FQ-sensitive and -resistant clones derived from the same patient, suggesting de novo emergence of FQ-resistant TB. Importantly, most of these samples were isolated from patients without a history of FQ treatment for TB. Thus, selective pressure applied by FQ monotherapy in the setting of undiagnosed TB infection appears to be able to drive the full or partial emergence of FQ-resistant M. tuberculosis, which has the potential to confound diagnostic tests for antibiotic susceptibility and limit the effectiveness of FQs in TB treatment.

Methodological and Clinical Aspects of the Molecular Epidemiology of Mycobacterium tuberculosis and Other Mycobacteria

Molecular typing has revolutionized epidemiological studies of infectious diseases, including those of a mycobacterial etiology. With the advent of fingerprinting techniques, many traditional concepts regarding transmission, infectivity, or pathogenicity of mycobacterial bacilli have been revisited, and their conventional interpretations have been challenged. Since the mid-1990s, when the first typing methods were introduced, a plethora of other modalities have been proposed. So-called molecular epidemiology has become an essential subdiscipline of modern mycobacteriology. It serves as a resource for understanding the key issues in the epidemiology of tuberculosis and other mycobacterial diseases. Among these issues are disclosing sources of infection, quantifying recent transmission, identifying transmission links, discerning reinfection from relapse, tracking the geographic distribution and clonal expansion of specific strains, and exploring the genetic mechanisms underlying specific phenotypic traits, including virulence, organ tropism, transmissibility, or drug resistance. Since genotyping continues to unravel the biology of mycobacteria, it offers enormous promise in the fight against and prevention of the diseases caused by these pathogens. In this review, molecular typing methods for Mycobacterium tuberculosis and nontuberculous mycobacteria elaborated over the last 2 decades are summarized. The relevance of these methods to the epidemiological investigation, diagnosis, evolution, and control of mycobacterial diseases is discussed.

Recent developments in genomics, bioinformatics and drug discovery to combat emerging drug-resistant tuberculosis

Emergence of drug-resistant tuberculosis (DR-TB) is a big challenge in TB control. The delay in diagnosis of DR-TB leads to its increased transmission, and therefore prevalence. Recent developments in genomics have enabled whole genome sequencing (WGS) of Mycobacterium tuberculosis (M. tuberculosis) from 3-day-old liquid culture and directly from uncultured sputa, while new bioinformatics tools facilitate to determine DR mutations rapidly from the resulting sequences. The present drug discovery and development pipeline is filled with candidate drugs which have shown efficacy against DR-TB. Furthermore, some of the FDA-approved drugs are being evaluated for repurposing, and this approach appears promising as several drugs are reported to enhance efficacy of the standard TB drugs, reduce drug tolerance, or modulate the host immune response to control the growth of intracellular M. tuberculosis. Recent developments in genomics and bioinformatics along with new drug discovery collectively have the potential to result in synergistic impact leading to the development of a rapid protocol to determine the drug resistance profile of the infecting strain so as to provide personalized medicine. Hence, in this review, we discuss recent developments in WGS, bioinformatics and drug discovery to perceive how they would transform the management of tuberculosis in a timely manner.

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.

Multiscale Model of Mycobacterium tuberculosis Infection Maps Metabolite and Gene Perturbations to Granuloma Sterilization Predictions

Granulomas are a hallmark of tuberculosis. Inside granulomas, the pathogen Mycobacterium tuberculosis may enter a metabolically inactive state that is less susceptible to antibiotics. Understanding M. tuberculosis metabolism within granulomas could contribute to reducing the lengthy treatment required for tuberculosis and provide additional targets for new drugs. Two key adaptations of M. tuberculosis are a nonreplicating phenotype and accumulation of lipid inclusions in response to hypoxic conditions. To explore how these adaptations influence granuloma-scale outcomes in vivo, we present a multiscale in silico model of granuloma formation in tuberculosis. The model comprises host immunity, M. tuberculosis metabolism, M. tuberculosis growth adaptation to hypoxia, and nutrient diffusion. We calibrated our model to in vivo data from nonhuman primates and rabbits and apply the model to predict M. tuberculosis population dynamics and heterogeneity within granulomas. We found that bacterial populations are highly dynamic throughout infection in response to changing oxygen levels and host immunity pressures. Our results indicate that a nonreplicating phenotype, but not lipid inclusion formation, is important for long-term M. tuberculosis survival in granulomas. We used virtual M. tuberculosis knockouts to predict the impact of both metabolic enzyme inhibitors and metabolic pathways exploited to overcome inhibition. Results indicate that knockouts whose growth rates are below ∼66% of the wild-type growth rate in a culture medium featuring lipid as the only carbon source are unable to sustain infections in granulomas. By mapping metabolite- and gene-scale perturbations to granuloma-scale outcomes and predicting mechanisms of sterilization, our method provides a powerful tool for hypothesis testing and guiding experimental searches for novel antituberculosis interventions.

Elucidating population-wide mycobacterial replication dynamics at the single-cell level

Mycobacterium tuberculosis infections result in a spectrum of clinical outcomes, and frequently the infection persists in a latent, clinically asymptomatic state. The within-host bacterial population is likely to be heterogeneous, and it is thought that persistent mycobacteria arise from a small population of viable, but non-replicating (VBNR) cells. These are likely to be antibiotic tolerant and necessitate prolonged treatment. Little is known about these persistent mycobacteria, since they are very difficult to isolate. To address this, we have successfully developed a replication reporter system for use in M. tuberculosis. This approach, termed fluorescence dilution, exploits two fluorescent reporters; a constitutive reporter allows the tracking of bacteria, while an inducible reporter enables the measurement of bacterial replication. The application of fluorescence single-cell analysis to characterize intracellular M. tuberculosis identified a distinct subpopulation of non-growing mycobacteria in murine macrophages. The presence of VBNR and actively replicating mycobacteria was observed within the same macrophage after 48 h of infection. Furthermore, our results suggest that macrophage uptake resulted in enrichment of non- or slowly replicating bacteria (as revealed by d-cycloserine treatment); this population is likely to be highly enriched for persisters, based on its drug-tolerant phenotype. These results demonstrate the successful application of the novel dual fluorescence reporter system both in vitro and in macrophage infection models to provide a window into mycobacterial population heterogeneity.

Whole genome analysis of Mycobacterium tuberculosis isolates from recurrent episodes of tuberculosis, Finland, 1995-2013

Recurrent tuberculosis (TB) is caused by an endogenous re-activation of the same strain of Mycobacterium tuberculosis (relapse) or exogenous infection with a new strain (re-infection). Recurrence of TB in Finland was analysed in a population-based, 19-year study, and genotyping was used to define relapse and re-infection. The M. tuberculosis isolates from patients with suspected relapse were further analysed by whole genome sequencing (WGS) to determine the number and type of mutations occurring in the bacterial genome between the first and second disease episodes. In addition, publicly available tools (PhyResSE and SpolPred) were used to predict drug resistance and spoligotype profile from the WGS data. Of the 8299 notified TB cases, 48 (0.6%) patients had episodes classified as recurrent. Forty-two patients had more than one culture-confirmed TB episode, and isolates from two episodes in 21 patients were available for genotyping. In 18 patients, the M. tuberculosis isolates obtained from the first and second TB episodes had identical spoligotypes. The WGS analysis of the 36 M. tuberculosis isolates from the 18 suspected relapse patients (average time between isolates 2.8 years) revealed 0 to 38 single nucleotide polymorphisms (median 1, mean 3.78) between the first and second isolate. There seemed to be no direct relation between the number of years between the two isolates, or treatment outcome, and the number of single nucleotide polymorphisms. The results suggest that the mutation rate may depend on multiple host-, strain- and treatment-related factors.

Interpreting whole genome sequencing for investigating tuberculosis transmission: a systematic review

BACKGROUND

Whole genome sequencing (WGS) is becoming an important part of epidemiological investigations of infectious diseases due to greater resolution and cost reductions compared to traditional typing approaches. Many public health and clinical teams will increasingly use WGS to investigate clusters of potential pathogen transmission, making it crucial to understand the benefits and assumptions of the analytical methods for investigating the data. We aimed to understand how different approaches affect inferences of transmission dynamics and outline limitations of the methods.

METHODS

We comprehensively searched electronic databases for studies that presented methods used to interpret WGS data for investigating tuberculosis (TB) transmission. Two authors independently selected studies for inclusion and extracted data. Due to considerable methodological heterogeneity between studies, we present summary data with accompanying narrative synthesis rather than pooled analyses.

RESULTS

Twenty-five studies met our inclusion criteria. Despite the range of interpretation tools, the usefulness of WGS data in understanding TB transmission often depends on the amount of genetic diversity in the setting. Where diversity is small, distinguishing re-infections from relapses may be impossible; interpretation may be aided by the use of epidemiological data, examining minor variants and deep sequencing. Conversely, when within-host diversity is large, due to genetic hitchhiking or co-infection of two dissimilar strains, it is critical to understand how it arose. Greater understanding of microevolution and mixed infection will enhance interpretation of WGS data.

CONCLUSIONS

As sequencing studies have sampled more intensely and integrated multiple sources of information, the understanding of TB transmission and diversity has grown, but there is still much to be learnt about the origins of diversity that will affect inferences from these data. Public health teams and researchers should combine epidemiological, clinical and WGS data to strengthen investigations of transmission.

The transmission of Mycobacterium tuberculosis in high burden settings

Unacceptable levels of Mycobacterium tuberculosis transmission are noted in high burden settings and a renewed focus on reducing person-to-person transmission in these communities is needed. We review recent developments in the understanding of airborne transmission. We outline approaches to measure transmission in populations and trials and describe the Wells-Riley equation, which is used to estimate transmission risk in indoor spaces. Present research priorities include the identification of effective strategies for tuberculosis infection control, improved understanding of where transmission occurs and the transmissibility of drug-resistant strains, and estimates of the effect of HIV and antiretroviral therapy on transmission dynamics. When research is planned and interventions are designed to interrupt transmission, resource constraints that are common in high burden settings-including shortages of health-care workers-must be considered.

Whole genome sequencing reveals genomic heterogeneity and antibiotic purification in Mycobacterium tuberculosis isolates

BACKGROUND

Whole genome sequencing has revolutionised the interrogation of mycobacterial genomes. Recent studies have reported conflicting findings on the genomic stability of Mycobacterium tuberculosis during the evolution of drug resistance. In an age where whole genome sequencing is increasingly relied upon for defining the structure of bacterial genomes, it is important to investigate the reliability of next generation sequencing to identify clonal variants present in a minor percentage of the population. This study aimed to define a reliable cut-off for identification of low frequency sequence variants and to subsequently investigate genetic heterogeneity and the evolution of drug resistance in M. tuberculosis.

METHODS

Genomic DNA was isolated from single colonies from 14 rifampicin mono-resistant M. tuberculosis isolates, as well as the primary cultures and follow up MDR cultures from two of these patients. The whole genomes of the M. tuberculosis isolates were sequenced using either the Illumina MiSeq or Illumina HiSeq platforms. Sequences were analysed with an in-house pipeline.

RESULTS

Using next-generation sequencing in combination with Sanger sequencing and statistical analysis we defined a read frequency cut-off of 30% to identify low frequency M. tuberculosis variants with high confidence. Using this cut-off we demonstrated a high rate of genetic diversity between single colonies isolated from one population, showing that by using the current sequencing technology, single colonies are not a true reflection of the genetic diversity within a whole population and vice versa. We further showed that numerous heterogeneous variants emerge and then disappear during the evolution of isoniazid resistance within individual patients. Our findings allowed us to formulate a model for the selective bottleneck which occurs during the course of infection, acting as a genomic purification event.

CONCLUSIONS

Our study demonstrated true levels of genetic diversity within an M. tuberculosis population and showed that genetic diversity may be re-defined when a selective pressure, such as drug exposure, is imposed on M. tuberculosis populations during the course of infection. This suggests that the genome of M. tuberculosis is more dynamic than previously thought, suggesting preparedness to respond to a changing environment.

Use of bacterial whole-genome sequencing to investigate local persistence and spread in bovine tuberculosis

•

We performed whole genome sequencing (WGS) of Mycobacterium bovis for a single molecular (VNTR) type.

•

Under-sampling of one lineage was caused by switching between VNTR-types.

•

Pairwise SNP distances showed a weak genetic isolation by distance pattern.

•

Bayesian phylogeographic inference was feasible despite a low substitution rate.

•

WGS studies of M. bovis need to account for slow evolution and molecular type switching.

Mycobacterium bovis is the causal agent of bovine tuberculosis, one of the most important diseases currently facing the UK cattle industry. Here, we use high-density whole genome sequencing (WGS) in a defined sub-population of M. bovis in 145 cattle across 66 herd breakdowns to gain insights into local spread and persistence. We show that despite low divergence among isolates, WGS can in principle expose contributions of under-sampled host populations to M. bovis transmission. However, we demonstrate that in our data such a signal is due to molecular type switching, which had been previously undocumented for M. bovis. Isolates from farms with a known history of direct cattle movement between them did not show a statistical signal of higher genetic similarity. Despite an overall signal of genetic isolation by distance, genetic distances also showed no apparent relationship with spatial distance among affected farms over distances <5 km. Using simulations, we find that even over the brief evolutionary timescale covered by our data, Bayesian phylogeographic approaches are feasible. Applying such approaches showed that M. bovis dispersal in this system is heterogeneous but slow overall, averaging 2 km/year. These results confirm that widespread application of WGS to M. bovis will bring novel and important insights into the dynamics of M. bovis spread and persistence, but that the current questions most pertinent to control will be best addressed using approaches that more directly integrate WGS with additional epidemiological data.

Mycobacterium tuberculosis Is Resistant to Isoniazid at a Slow Growth Rate by Single Nucleotide Polymorphisms in katG Codon Ser315

An important aim for improving TB treatment is to shorten the period of antibiotic therapy without increasing relapse rates or encouraging the development of antibiotic-resistant strains. In any M. tuberculosis population there is a proportion of bacteria that are drug-tolerant; this might be because of pre-existing populations of slow growing/non replicating bacteria that are protected from antibiotic action due to the expression of a phenotype that limits drug activity. We addressed this question by observing populations of either slow growing (constant 69.3h mean generation time) or fast growing bacilli (constant 23.1h mean generation time) in their response to the effects of isoniazid exposure, using controlled and defined growth in chemostats. Phenotypic differences were detected between the populations at the two growth rates including expression of efflux mechanisms and the involvement of antisense RNA/small RNA in the regulation of a drug-tolerant phenotype, which has not been explored previously for M. tuberculosis. Genotypic analyses showed that slow growing bacilli develop resistance to isoniazid through mutations specifically in katG codon Ser³¹⁵ which are present in approximately 50–90% of all isoniazid-resistant clinical isolates. The fast growing bacilli persisted as a mixed population with katG mutations distributed throughout the gene. Mutations in katG codon Ser³¹⁵ appear to have a fitness cost in vitro and particularly in fast growing cultures. Our results suggest a requirement for functional katG-encoded catalase-peroxide in the slow growers but not the fast-growing bacteria, which may explain why katG codon Ser³¹⁵ mutations are favoured in the slow growing cultures.

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.

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.