High functional diversity in Mycobacterium tuberculosis driven by genetic drift and human demography

Bacterial diversity dominates variable macrophage responses of tuberculosis patients in Tanzania

The Mycobacterium tuberculosis complex (MTBC) comprises nine human-adapted lineages that differ in their geographical distribution. Local adaptation of specific MTBC genotypes to the respective human host population has been invoked in this context. We aimed to assess if bacterial genetics governs MTBC pathogenesis or if local co-adaptation translates into differential susceptibility of human macrophages to infection by different MTBC genotypes. We generated macrophages from cryopreserved blood mononuclear cells of Tanzanian tuberculosis patients, from which the infecting MTBC strains had previously been phylogenetically characterized. We infected these macrophages ex vivo with a phylogenetically similar MTBC strain ("matched infection") or with strains representative of other MTBC lineages ("mismatched infection"). We found that L1 infections resulted in a significantly lower bacterial burden and that the intra-cellular replication rate of L2 strains was significantly higher compared the other MTBC lineages, irrespective of the MTBC lineage originally infecting the patients. Moreover, L4-infected macrophages released significantly greater amounts of TNF-α, IL-6, IL-10, MIP-1β, and IL-1β compared to macrophages infected by all other strains. While our results revealed no measurable effect of local adaptation, they further highlight the strong impact of MTBC phylogenetic diversity on the variable outcome of the host-pathogen interaction in human tuberculosis.

Evolution and emergence of Mycobacterium tuberculosis

Tuberculosis (TB) remains one of the deadliest infectious diseases in human history, prevailing even in the 21st century. The causative agents of TB are represented by a group of closely related bacteria belonging to the Mycobacterium tuberculosis complex (MTBC), which can be subdivided into several lineages of human- and animal-adapted strains, thought to have shared a last common ancestor emerged by clonal expansion from a pool of recombinogenic Mycobacterium canettii-like tubercle bacilli. A better understanding of how MTBC populations evolved from less virulent mycobacteria may allow for discovering improved TB control strategies and future epidemiologic trends. In this review, we highlight new insights into the evolution of mycobacteria at the genus level, describing different milestones in the evolution of mycobacteria, with a focus on the genomic events that have likely enabled the emergence and the dominance of the MTBC. We also review the recent literature describing the various MTBC lineages and highlight their particularities and differences with a focus on host preferences and geographic distribution. Finally, we discuss on putative mechanisms driving the evolution of tubercle bacilli and mycobacteria in general, by taking the mycobacteria-specific distributive conjugal transfer as an example.

Polyketide synthases mutation in tuberculosis transmission revealed by whole genomic sequence, China, 2011-2019

Introduction: Tuberculosis (TB) is an infectious disease caused by a bacterium called Mycobacterium tuberculosis (Mtb). Previous studies have primarily focused on the transmissibility of multidrug-resistant (MDR) or extensively drug-resistant (XDR) Mtb. However, variations in virulence across Mtb lineages may also account for differences in transmissibility. In Mtb, polyketide synthase (PKS) genes encode large multifunctional proteins which have been shown to be major mycobacterial virulence factors. Therefore, this study aimed to identify the role of PKS mutations in TB transmission and assess its risk and characteristics. Methods: Whole genome sequences (WGSs) data from 3,204 Mtb isolates was collected from 2011 to 2019 in China. Whole genome single nucleotide polymorphism (SNP) profiles were used for phylogenetic tree analysis. Putative transmission clusters (≤10 SNPs) were identified. To identify the role of PKS mutations in TB transmission, we compared SNPs in the PKS gene region between "clustered isolates" and "non-clustered isolates" in different lineages. Results: Cluster-associated mutations in ppsA, pks12, and pks13 were identified among different lineage isolates. They were statistically significant among clustered strains, indicating that they may enhance the transmissibility of Mtb. Conclusion: Overall, this study provides new insights into the function of PKS and its localization in M. tuberculosis. The study found that ppsA, pks12, and pks13 may contribute to disease progression and higher transmission of certain strains. We also discussed the prospective use of mutant ppsA, pks12, and pks13 genes as drug targets.

Probable cases of tuberculosis in Early Medieval pastoralists of Eastern Europe

Two anthropological collections of the 8-9th century AD from the forest-steppe area of the southern Eastern Europe has been the subject of analysis of the osteo-articular alterations concerning tuberculosis related lesions. According to archaeological data, the Mayaki and Dmitrievka groups are from the territory of Khazar Khaganate, which is well-known by historical data. It was a traditional population with settled pastoral farming. There were studied 292 adult human remains by macroscopic morphological and radiological methods. The completeness of the skeletons could have been better due to the intentional selection of their parts for the museum funds, so there were limitations for the differential diagnosis of TB by morphological criteria. In general, 31 individuals (20 men, 11 women) were marked as possible candidates with skeletal TB. The number of skeletons with possible specific infections from the two series is comparable, but according to sex distribution, there are some significant differences. If in the Dmitrievka group, there was approximately the same distribution of the number of cases of possible infected men and women; in the Mayaki group then there was a significant difference in the prevalence of the male sample. Various reasons can explain the observed differences. Firstly, we cannot rule out an error in the diagnosis of tuberculosis only by morphological methods since the preservation of the skeletons leaves much to be desired. Secondly, as it is known, only a few percent of tuberculosis patients show skeletal alterations, so the direct quantitative comparison cannot be adequate in paleopathological studies. The ambiguity of morphological criteria, particularly for destructive lesions of bone other than vertebrae, does not allow asserting that all identified cases result from TB. However, the results of differential diagnosis can underline that at least half of the skeletons have skeletal alterations most likely related to TB. The final analysis of the remains by aDNA methods will permit more certain confirmation of the diagnosis.

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.

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

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

Inhibitors of Aspartate Transcarbamoylase Inhibit Mycobacterium tuberculosis Growth

Aspartate transcarbamoylase (ATCase) plays a key role in the second step of de novo pyrimidine biosynthesis in eukaryotes and has been proposed to be a target to suppress cell proliferation in E. coli, human cells and the malarial parasite. We hypothesized that a library of ATCase inhibitors developed for malarial ATCase (PfATCase) may also contain inhibitors of the tubercular ATCase and provide a similar inhibition of cellular proliferation. Of the 70 compounds screened, 10 showed single-digit micromolar inhibition in an in vitro activity assay and were tested for their effect on M. tuberculosis cell growth in culture. The most promising compound demonstrated a MIC90 of 4 μM. A model of MtbATCase was generated using the experimental coordinates of PfATCase. In silico docking experiments showed this compound can occupy a similar allosteric pocket on MtbATCase to that seen on PfATCase, explaining the observed species selectivity seen for this compound series.

Investigating the composition and recruitment of the mycobacterial ImuA'-ImuB-DnaE2 mutasome

A DNA damage-inducible mutagenic gene cassette has been implicated in the emergence of drug resistance in Mycobacterium tuberculosis during anti-tuberculosis (TB) chemotherapy. However, the molecular composition and operation of the encoded 'mycobacterial mutasome' - minimally comprising DnaE2 polymerase and ImuA' and ImuB accessory proteins - remain elusive. Following exposure of mycobacteria to DNA damaging agents, we observe that DnaE2 and ImuB co-localize with the DNA polymerase III β subunit (β clamp) in distinct intracellular foci. Notably, genetic inactivation of the mutasome in an imuBAAAAGG mutant containing a disrupted β clamp-binding motif abolishes ImuB-β clamp focus formation, a phenotype recapitulated pharmacologically by treating bacilli with griselimycin and in biochemical assays in which this β clamp-binding antibiotic collapses pre-formed ImuB-β clamp complexes. These observations establish the essentiality of the ImuB-β clamp interaction for mutagenic DNA repair in mycobacteria, identifying the mutasome as target for adjunctive therapeutics designed to protect anti-TB drugs against emerging resistance.

Lack of methoxy-mycolates characterizes the geographically restricted lineage 7 of Mycobacterium tuberculosis complex

Lineage 7 (L7) emerged in the phylogeny of the Mycobacterium tuberculosis complex (MTBC) subsequent to the branching of 'ancient' lineage 1 and prior to the Eurasian dispersal of 'modern' lineages 2, 3 and 4. In contrast to the major MTBC lineages, the current epidemiology suggests that prevalence of L7 is highly confined to the Ethiopian population, or when identified outside of Ethiopia, it has mainly been in patients of Ethiopian origin. To search for microbiological factors that may contribute to its restricted distribution, we compared the genome of L7 to the genomes of globally dispersed MTBC lineages. The frequency of predicted functional mutations in L7 was similar to that documented in other lineages. These include mutations characteristic of modern lineages - such as constitutive expression of nitrate reductase - as well as mutations in the VirS locus that are commonly found in ancient lineages. We also identified and characterized multiple lineage-specific mutations in L7 in biosynthesis pathways of cell wall lipids, including confirmed deficiency of methoxy-mycolic acids due to a stop-gain mutation in the mmaA3 gene that encodes a methoxy-mycolic acid synthase. We show that the abolished biosynthesis of methoxy-mycolates of L7 alters the cell structure and colony morphology on selected growth media and impacts biofilm formation. The loss of these mycolic acid moieties may change the host-pathogen dynamic for L7 isolates, explaining the limited geographical distribution of L7 and contributing to further understanding the spread of MTBC lineages across the globe.

Liposomal Delivery of Saquinavir to Macrophages Overcomes Cathepsin Blockade by Mycobacterium tuberculosis and Helps Control the Phagosomal Replicative Niches

Mycobacterium tuberculosis is able to establish a chronic colonization of lung macrophages in a controlled replication manner, giving rise to a so-called latent infection. Conversely, when intracellular bacteria undergo actively uncontrolled replication rates, they provide the switch for the active infection called tuberculosis to occur. Our group found that the pathogen is able to manipulate the activity of endolysosomal enzymes, cathepsins, directly at the level of gene expression or indirectly by regulating their natural inhibitors, cystatins. To provide evidence for the crucial role of cathepsin manipulation for the success of tuberculosis bacilli in their intracellular survival, we used liposomal delivery of saquinavir. This protease inhibitor was previously found to be able to increase cathepsin proteolytic activity, overcoming the pathogen induced blockade. In this study, we demonstrate that incorporation in liposomes was able to increase the efficiency of saquinavir internalization in macrophages, reducing cytotoxicity at higher concentrations. Consequently, our results show a significant impact on the intracellular killing not only to reference and clinical strains susceptible to current antibiotic therapy but also to multidrug- and extensively drug-resistant (XDR) Mtb strains. Altogether, this indicates the manipulation of cathepsins as a fine-tuning strategy used by the pathogen to survive and replicate in host cells.

Bedaquiline resistant Mycobacterium tuberculosis clinical isolates with and without rv0678 mutations have similar growth patterns under varying BDQ drug pressure

Resistance associated mutations have been reported to alter the growth of Mycobacterium tuberculosis (MTB) isolates under drug pressure. However, there is little information on the growth characteristics of bedaquiline (BDQ) resistant isolates in the presence of BDQ. To further understand the role of rv0678, we aimed to study whether the presence of rv0678 variants in BDQ resistant isolates alters the killing effect of BDQ. We, therefore, selected BDQ resistant clinical MTB isolates with (n = 6) and without (n = 3) variants in rv0678 gene. Using time kill assays, growth inhibition; taken as the relative change in log average colony forming unit (CFU)/ml at selected time points (24-96 h), was studied at Minimum Inhibitory Concentrations (MICs): 0x, 1x, 2.5x, 5x, 7.5x, 10x for these isolates. Growth inhibition was then analyzed using Kruskal Wallis and Kolmogorov Smirnov tests in PRISM vr.9. During the 24-96 h lag phase isolates with and without variants in rv0678 showed a similar growth inhibition pattern. No difference was noted in growth inhibition between BDQ resistant isolates and H37Rv. These findings suggest that role of alternate mechanisms in contributing to BDQ tolerance needs to be explored.

Biocultural perspectives of infectious diseases and demographic evolution: Tuberculosis and its comorbidities through history

Anthropologists recognize the importance of conceptualizing health in the context of the mutually evolving nature of biology and culture through the biocultural approach, but biocultural anthropological perspectives of infectious diseases and their impacts on humans (and vice versa) through time are relatively underrepresented. Tuberculosis (TB) has been a constant companion of humans for thousands of years and has heavily influenced population health in almost every phase of cultural and demographic evolution. TB in human populations has been dramatically influenced by behavior, demographic and epidemiological shifts, and other comorbidities through history. This paper critically discusses TB and some of its major comorbidities through history within a biocultural framework to show how transitions in human demography and culture affected the disease-scape of TB. In doing so, I address the potential synthesis of biocultural and epidemiological transition theory to better comprehend the mutual evolution of infectious diseases and humans.

An ancestral mycobacterial effector promotes dissemination of infection

The human pathogen Mycobacterium tuberculosis typically causes lung disease but can also disseminate to other tissues. We identified a M. tuberculosis (Mtb) outbreak presenting with unusually high rates of extrapulmonary dissemination and bone disease. We found that the causal strain carried an ancestral full-length version of the type VII-secreted effector EsxM rather than the truncated version present in other modern Mtb lineages. The ancestral EsxM variant exacerbated dissemination through enhancement of macrophage motility, increased egress of macrophages from established granulomas, and alterations in macrophage actin dynamics. Reconstitution of the ancestral version of EsxM in an attenuated modern strain of Mtb altered the migratory mode of infected macrophages, enhancing their motility. In a zebrafish model, full-length EsxM promoted bone disease. The presence of a derived nonsense variant in EsxM throughout the major Mtb lineages 2, 3, and 4 is consistent with a role for EsxM in regulating the extent of dissemination.

The rate and role of pseudogenes of the Mycobacterium tuberculosis complex

Whole-genome sequence analyses have significantly contributed to the understanding of virulence and evolution of the Mycobacterium tuberculosis complex (MTBC), the causative pathogens of tuberculosis. Most MTBC evolutionary studies are focused on single nucleotide polymorphisms and deletions, but rare studies have evaluated gene content, whereas none has comprehensively evaluated pseudogenes. Accordingly, we describe an extensive study focused on quantifying and predicting possible functions of MTBC and Mycobacterium canettii pseudogenes. Using NCBI's PGAP-detected pseudogenes, we analysed 25 837 pseudogenes from 158 MTBC and M. canetii strains and combined transcriptomics and proteomics of M. tuberculosis H37Rv to gain insights about pseudogenes' expression. Our results indicate significant variability concerning rate and conservancy of in silico predicted pseudogenes among different ecotypes and lineages of tuberculous mycobacteria and pseudogenization of important virulence factors and genes of the metabolism and antimicrobial resistance/tolerance. We show that in silico predicted pseudogenes contribute considerably to MTBC genetic diversity at the population level. Moreover, the transcription machinery of M. tuberculosis can fully transcribe most pseudogenes, indicating intact promoters and recent pseudogene evolutionary emergence. Proteomics of M. tuberculosis and close evaluation of mutational lesions driving pseudogenization suggest that few in silico predicted pseudogenes are likely capable of neofunctionalization, nonsense mutation reversal, or phase variation, contradicting the classical definition of pseudogenes. Such findings indicate that genome annotation should be accompanied by proteomics and protein function assays to improve its accuracy. While indels and insertion sequences are the main drivers of the observed mutational lesions in these species, population bottlenecks and genetic drift are likely the evolutionary processes acting on pseudogenes' emergence over time. Our findings unveil a new perspective on MTBC's evolution and genetic diversity.

Evidence, Challenges, and Knowledge Gaps Regarding Latent Tuberculosis in Animals

Mycobacterium bovis and other Mycobacterium tuberculosis complex (MTBC) pathogens that cause domestic animal and wildlife tuberculosis have received considerably less attention than M. tuberculosis, the primary cause of human tuberculosis (TB). Human TB studies have shown that different stages of infection can exist, driven by host–pathogen interactions. This results in the emergence of heterogeneous subpopulations of mycobacteria in different phenotypic states, which range from actively replicating (AR) cells to viable but slowly or non-replicating (VBNR), viable but non-culturable (VBNC), and dormant mycobacteria. The VBNR, VBNC, and dormant subpopulations are believed to underlie latent tuberculosis (LTB) in humans; however, it is unclear if a similar phenomenon could be happening in animals. This review discusses the evidence, challenges, and knowledge gaps regarding LTB in animals, and possible host–pathogen differences in the MTBC strains M. tuberculosis and M. bovis during infection. We further consider models that might be adapted from human TB research to investigate how the different phenotypic states of bacteria could influence TB stages in animals. In addition, we explore potential host biomarkers and mycobacterial changes in the DosR regulon, transcriptional sigma factors, and resuscitation-promoting factors that may influence the development of LTB.

Whole Genome Sequencing Assessing Impact of Diabetes Mellitus on Tuberculosis Mutations and Type of Recurrence in India

BACKGROUND

Evidence describing the impact of diabetes mellitus (DM) on the recurrence and mutation rate of Mycobacterium tuberculosis (Mtb) is limited.

METHODS

This study was nested in 3 cohort studies of tuberculosis (TB) patients with and without DM in India. Paired Mtb isolates recovered at baseline and treatment failure/recurrence underwent whole genome sequencing. We compared acquisition of single-nucleotide polymorphisms (SNPs), TB drug resistance mutations, and type of recurrence (endogenous reactivation [<8 SNPs] or exogenous reinfection [≥8 SNPs]) by DM status.

RESULTS

Of 1633 enrolled in the 3 parent cohorts, 236 (14.5%) had microbiologically confirmed TB treatment failure/recurrence; 76 Mtb isolate pairs were available for sequencing (22 in TB-DM and 54 in TB-only). The SNP acquisition rate was overall was 0.43 (95% confidence interval [CI], .25-.64) per 1 person-year (PY); 0.77 (95% CI, .40-1.35) per 1 PY, and 0.44 (95% CI, .19-.86) per 1 PY at treatment failure and recurrence, respectively. Significant difference in SNP rates by DM status was seen at recurrence (0.21 [95% CI, .04-.61]) per 1 PY for TB-only vs 1.28 (95% CI, .41-2.98) per 1 PY for TB-DM; P = .02). No significant difference in SNP rates by DM status was observed at treatment failure. Acquired TB drug resistance was seen in 4 of 18 (22%) in TB-DM vs 4 of 45 (9%) in TB-only (P = .21). Thirteen (17%) participants had exogenous reinfection; the reinfection rate at recurrence was 25% (3/12) for TB-DM vs 17% (4/24) in TB-only (P = .66).

CONCLUSIONS

Considerable intrahost Mtb mutation rates were present at recurrence among patients with DM in India. One-fourth of patients with DM had exogenous reinfection at recurrence.

Genomic analysis of diversity, biogeography, and drug resistance in Mycobacterium bovis

Mycobacterium bovis is the cause of bovine tuberculosis, and it can also cause disease in humans, with symptoms similar to those caused by Mycobacterium tuberculosis. However, our understanding of its genomic diversity, biogeography, and drug resistance remains incomplete. We performed a comparative and phylogenetic analysis of 3,228 M. bovis genomes from 24 countries. Following drug susceptibility testing, we applied a bacterial genome-wide association study to capture associations between genomic variation and drug resistance in 74 newly isolated strains from China. The data show that the cattle-adapted M. bovis were divided into six lineages with a strong phylogeographical population structure. Lineage 1 and Lineage 6 are the most widespread globally, while others show a strong geographical restriction. 17.39% of M. bovis isolates were resistant to at least one drug in China. Furthermore, we identify genomic variations associated with an increased risk of resistance acquisition. This study furthers our knowledge of M. bovis diversity, biogeography, and drug resistance and will facilitate more deeply informed genomic tracking and surveillance to minimize its threat to human health, as a cause of zoonotic tuberculosis. This article is protected by copyright. All rights reserved.

Tuberculosis caused by Mycobacterium africanum: Knowns and unknowns

Tuberculosis (TB), one of the deadliest threats to human health, is mainly caused by 2 highly related and human-adapted bacteria broadly known as Mycobacterium tuberculosis and Mycobacterium africanum. Whereas M. tuberculosis is widely spread, M. africanum is restricted to West Africa, where it remains a significant cause of tuberculosis. Although several differences have been identified between these 2 pathogens, M. africanum remains a lot less studied than M. tuberculosis. Here, we discuss the genetic, phenotypic, and clinical similarities and differences between strains of M. tuberculosis and M. africanum. We also discuss our current knowledge on the immune response to M. africanum and how it possibly articulates with distinct disease progression and with the geographical restriction attributed to this pathogen. Understanding the functional impact of the diversity existing in TB-causing bacteria, as well as incorporating this diversity in TB research, will contribute to the development of better, more specific approaches to tackle TB.

Emergence of multi-drug-resistant Mycobacterium tuberculosis in Niger: A snapshot based on whole-genome sequencing

Background

Among other West African countries experiencing the high endemicity of deadly tuberculosis, the situation in Niger is poorly evidenced by microbiological investigations.

Methodology/Principal findings

The study of 42 isolates of Mycobacterium tuberculosis from Niger by whole genome sequencing using Illumina iSeq technology yielded four M. tuberculosis lineages: Indo-Oceanic L1 (n = 1) (2.3%), East-Asian (n = 1) (2.3%), East-African Indian L3 (n = 2) (4.7%) and Euro-American L4 (n = 38) (90.4%). The sub-lineage L4.1.3 comprising 18 isolates (47.3%) was predominant, followed by the L4.6.2.2 sub-lineage (Cameroon genotype, n = 13 isolates) (34.2%). Investigating drug resistance profile for 12 antibiotics found 8/42 (19%) pan-susceptible isolates and 34/42 (81%) resistant isolates; with 40/42 (95.2%) isolates being susceptible to clofazimine-bedaquiline.

Conclusions/Significance

These unprecedented data from Niger highlight the dynamics of tuberculosis transmission and drug resistance in Niger and may assist tuberculosis control in this country which continues to support a high burden of tuberculosis.

Tuberculosis is a major public health problem in Niger, in West Africa. Niger has an emerging problem with multidrug-resistant tuberculosis (MDR-TB). Whole genome sequencing was used to understand the epidemiology of tuberculosis and genetics of multi-drug resistance among patients from the regions in Niger. In this study, most isolates of M. tuberculosis from this dataset belonged to the L4.6.2.2 sub-lineage and L4.1.3 sub-lineage within the Euro-American lineage. Thirty-four out of 42 (81%) isolates were detected as resistant isolates. Our study highlights the need for epidemiological surveillance and more concerted efforts to ensure that patients are put through treatment.

Gene evolutionary trajectories in Mycobacterium tuberculosis reveal temporal signs of selection

SignificancePrevious attempts to identify the action of natural selection in the Mycobacterium tuberculosis complex (MTBC) were limited by sample size and averaging across time and lineages. We investigate changes in selective pressures across time for every single gene of the MTBC. We developed a methodology to analyze temporal signals of selection in a large dataset (∼5,000 complete genomes) and showed that 1) almost half of the genes seem to have been under positive selection at some point in time; 2) experimentally confirmed epitopes tend to accumulate more mutations in deeper branches than in external branches; and 3) temporal signals identify genes that were conserved in the past but under positive selection in the present, suggesting ongoing adaptation to the host.

The TbD1 Locus Mediates a Hypoxia-Induced Copper Response in Mycobacterium bovis

The Mycobacterium tuberculosis complex (MTBC) contains the causative agents of tuberculosis (TB) in mammals. The archetypal members of the MTBC, Mycobacterium tuberculosis and Mycobacterium bovis, cause human tuberculosis and bovine tuberculosis, respectively. Although M. tuberculosis and M. bovis share over 99.9% genome identity, they show distinct host adaptation for humans and animals; hence, while the molecular basis of host adaptation is encoded in their genomes, the mechanistic basis of host tropism is still unclear. Exploration of the in vitro phenotypic consequences of known genetic difference between M. bovis and M. tuberculosis offers one route to explore genotype–phenotype links that may play a role in host adaptation. The TbD1 (“Mycobacterium tuberculosis deletion 1 region”) locus encompasses the mmpS6 and mmpL6 genes. TbD1 is absent in M. tuberculosis “modern” lineages (Lineages 2, 3, and 4) but present in “ancestral” M. tuberculosis (Lineages 1 and 7), Mycobacterium africanum lineages (Lineages 5 and 6), newly identified M. tuberculosis lineages (Lineages 8 and 9), and animal adapted strains, such as M. bovis. The function of TbD1 has previously been investigated in M. tuberculosis, where conflicting data has emerged on the role of TbD1 in sensitivity to oxidative stress, while the underlying mechanistic basis of such a phenotype is unclear. In this study, we aimed to shed further light on the role of the TbD1 locus by exploring its function in M. bovis. Toward this, we constructed an M. bovis TbD1 knockout (ΔTbD1) strain and conducted comparative transcriptomics to define global gene expression profiles of M. bovis wild-type (WT) and the ΔTbD1 strains under in vitro culture conditions (rolling and standing cultures). This analysis revealed differential induction of a hypoxia-driven copper response in WT and ΔTbD1 strains. In vitro phenotypic assays demonstrated that the deletion of TbD1 sensitized M. bovis to H2O2 and hypoxia-specific copper toxicity. Our study provides new information on the function of the TbD1 locus in M. bovis and its role in stress responses in the MTBC.

Mycobacterium tuberculosis DNA repair helicase UvrD1 is activated by redox-dependent dimerization via a 2B domain cysteine

Mycobacterium tuberculosis (Mtb) causes tuberculosis and, during infection, is exposed to reactive oxygen species and reactive nitrogen intermediates from the host immune response that can cause DNA damage. UvrD-like proteins are involved in DNA repair and replication and belong to the SF1 family of DNA helicases that use ATP hydrolysis to catalyze DNA unwinding. In Mtb, there are two UvrD-like enzymes, where UvrD1 is most closely related to other family members. Previous studies have suggested that UvrD1 is exclusively monomeric; however, it is well known that Escherichia coli UvrD and other UvrD family members exhibit monomer-dimer equilibria and unwind as dimers in the absence of accessory factors. Here, we reconcile these incongruent studies by showing that Mtb UvrD1 exists in monomer, dimer, and higher-order oligomeric forms, where dimerization is regulated by redox potential. We identify a 2B domain cysteine, conserved in many Actinobacteria, that underlies this effect. We also show that UvrD1 DNA-unwinding activity correlates specifically with the dimer population and is thus titrated directly via increasing positive (i.e., oxidative) redox potential. Consistent with the regulatory role of the 2B domain and the dimerization-based activation of DNA unwinding in UvrD family helicases, these results suggest that UvrD1 is activated under oxidizing conditions when it may be needed to respond to DNA damage during infection.

Comparative genomic analyses of multi-drug resistant Mycobacterium tuberculosis isolates from Nepal and other geographical locations

Nepal exhibits a tuberculosis (TB) incidence rate that is comparable to neighbouring high TB incidence countries. In addition, it records >500 cases of multi-drug resistant (MDR) TB each year. The objective of this study was to perform whole-genome bioinformatic analysis on MDR-TB isolates from Nepal (n = 19) to identify the specific mutations underlying their phenotypic resistance. In addition, we examined the dominant genotype among the Nepal MDR-TB isolates, the East-Asian Beijing sub-lineage, to determine its relatedness to a panel of 1274 genomes of international strains available from public databases. These analyses provided evidence that the XDR-TB isolates in our collection were not derived from importation of primary XDR-TB to Nepal but were more likely the result of acquisition of second-line drug resistance in Nepal. Resistance to fluoroquinolones was detected among a high proportion of the Nepal isolates. This has implications for the management of TB, including appropriate antimicrobial stewardship and susceptibility testing for fluoroquinolones and other second-line TB drugs, to minimise the development of XDR-TB among Nepal TB cases.

Host-pathogen genetic interactions underlie tuberculosis susceptibility in genetically diverse mice

The outcome of an encounter with Mycobacterium tuberculosis (Mtb) depends on the pathogen’s ability to adapt to the variable immune pressures exerted by the host. Understanding this interplay has proven difficult, largely because experimentally tractable animal models do not recapitulate the heterogeneity of tuberculosis disease. We leveraged the genetically diverse Collaborative Cross (CC) mouse panel in conjunction with a library of Mtb mutants to create a resource for associating bacterial genetic requirements with host genetics and immunity. We report that CC strains vary dramatically in their susceptibility to infection and produce qualitatively distinct immune states. Global analysis of Mtb transposon mutant fitness (TnSeq) across the CC panel revealed that many virulence pathways are only required in specific host microenvironments, identifying a large fraction of the pathogen’s genome that has been maintained to ensure fitness in a diverse population. Both immunological and bacterial traits can be associated with genetic variants distributed across the mouse genome, making the CC a unique population for identifying specific host-pathogen genetic interactions that influence pathogenesis.

Era of TB elimination: Growing need to understand diversities of Mycobacterium tuberculosis lineages!

INTRODUCTION

The mycobacterium tuberculosis complex (MTBC) has highly clonal population structure which made the organism spread globally mirroring human migration out of Africa and resulted in the formation of seven lineages. We conducted this study to determine the proportion of spoligotype lineages and drug susceptibility profile of Mycobacterium tuberculosis isolates among smear positive TB patients attending a tertiary care hospital in Mysore, Karnataka, India.

METHODS

It is a descriptive study conducted at JSS Hospital a tertiary care centre at Mysore, India during 2018-19. The sputum smear positive samples were subjected to solid culture and drug susceptibility testing and spoligotyping for identification of lineages.

RESULTS

Of the 100 samples which were culture positive, 94 isolates were clustered into five spoligotype international types with SIT-126 (EAI-5) being the largest cluster of 46 (46%) isolates, followed by SIT-62 (H1) with 24 (24%), SIT -26 (CAS 1-DELHI) with 20 (20%), SIT-53 (T1) with 03 (3%) and SIT-482 (BOV-1) with 01 (1%). Among the remaining six isolates, two had unique Cameroon spoligotypes and four were orphans CONCLUSION: The study finding reveals that a diverse pattern of genotypes is circulating in the region of which EAI-5, Harleem (H1) and CAS-DELHI pattern forms the majority (88%). It is evident that there is a wide range of MTB genetic lineages in circulation and further research is needed to understand the diversity across the country.

DNA Methyltransferase HsdM Induce Drug Resistance on Mycobacterium tuberculosis via Multiple Effects

Besides the genomic variants, epigenetic mechanisms such as DNA methylation also have an effect on drug resistance. This study aimed to investigate the methylomes of totally/extensively drug-resistant M. tuberculosis clinical isolates using the PacBio single-molecule real-time technology. The results showed they were almost the same as the pan-susceptible ones. Genetics and bioinformatics analysis confirmed three DNA methyltransferases-MamA, MamB, and HsdM. Moreover, anti-tuberculosis drug treatment did not change the methylomes. In addition, the knockout of the DNA methyltransferase hsdM gene in the extensively drug-resistant clinical isolate 11826 revealed that the motifs of GTAYN4ATC modified by HsdM were completely demethylated. Furthermore, the results of the methylated DNA target analysis found that HsdM was mainly involved in redox-related pathways, especially the prodrug isoniazid active protein KatG. HsdM also targeted three drug-targeted genes, eis, embB, and gyrA, and three drug transporters (Rv0194, Rv1410, and Rv1877), which mildly affected the drug susceptibility. The overexpression of HsdM in M. smegmatis increased the basal mutation rate. Our results suggested that DNA methyltransferase HsdM affected the drug resistance of M. tuberculosis by modulating the gene expression of redox, drug targets and transporters, and gene mutation.

Badania kopalnego DNA – możliwości i ograniczenia

Ostatnie cztery dekady przyniosły znaczący rozwój archeologii molekularnej i badania nad kopalnym DNA (aDNA). Nowatorskie metody uwzględniają szeroki zakres badań, począwszy od sekwencjonowania niewielkich fragmentów mitochondrialnego DNA po wielkoskalowe badania całych populacji, łączące sekwencjonowanie genomów mitochondrialnych, genów podlegających doborowi naturalnemu, jak i całych genomów jądrowych. Postęp, zwłaszcza w dziedzinie technologii sekwencjonowania DNA, umożliwił pozyskanie informacji ze szczątków paleontologicznych i materiału archeologicznego, umożliwiając zbadanie związków filogenetycznych między wymarłymi i współczesnymi gatunkami. Dzięki zastosowaniu technologii sekwencjonowania nowej generacji możliwe stało się poznanie sekwencji DNA nie tylko bezpośrednio ze szczątków ludzkich lub zwierzęcych, ale także z osadów sedymentacyjnych z głębin jezior oraz jaskiń. W artykule przedstawiono możliwości i ograniczenia występujące w badaniach nad kopalnym DNA ludzi, zwierząt czy bakterii z podkreśleniem wkładu polskich badaczy w rozwój tej dziedziny nauki.

Structural implications of lipoarabinomannan glycans from global clinical isolates in diagnosis of Mycobacterium tuberculosis infection

In Mycobacterium tuberculosis (Mtb), surface-exposed Lipoarabinomannan (LAM) is a key determinant of immunogenicity, yet its intrinsic heterogeneity confounds typical structure–function analysis. Recently, LAM gained a strong foothold as a validated marker for active tuberculosis (TB) infection and has shown great potential in new diagnostic efforts. However, no efforts have yet been made to model or evaluate the impact of mixed polyclonal Mtb infections (infection with multiple strains) on TB diagnostic procedures other than antibiotic susceptibility testing. Here, we selected three TB clinical isolates (HN878, EAI, and IO) and purified LAM from these strains to present an integrated analytical approach of one-dimensional and two-dimensional Nuclear Magnetic Resonance (NMR) spectroscopy, as well as enzymatic digestion and site-specific mass spectrometry (MS) to probe LAM structure and behavior at multiple levels. Overall, we found that the glycan was similar in all LAM preparations, albeit with subtle variations. Succinates, lactates, hydroxybutyrate, acetate, and the hallmark of Mtb LAM-methylthioxylose (MTX), adorned the nonreducing terminal arabinan of these LAM species. Newly identified acetoxy/hydroxybutyrate was present only in LAM from EAI and IO Mtb strains. Notably, detailed LC/MS-MS unambiguously showed that all acyl modifications and the lactyl ether in LAM are at the 3-OH position of the 2-linked arabinofuranose adjacent to the terminal β-arabinofuranose. Finally, after sequential enzymatic deglycosylation of LAM, the residual glycan that has ∼50% of α−arabinofuranose -(1→5) linked did not bind to monoclonal antibody CS35. These data clearly indicate the importance of the arabinan termini arrangements for the antigenicity of LAM.

Genome-wide estimation of recombination, mutation and positive selection enlightens diversification drivers of Mycobacterium bovis

Genome sequencing has reinvigorated the infectious disease research field, shedding light on disease epidemiology, pathogenesis, host-pathogen interactions and also evolutionary processes exerted upon pathogens. Mycobacterium tuberculosis complex (MTBC), enclosing M. bovis as one of its animal-adapted members causing tuberculosis (TB) in terrestrial mammals, is a paradigmatic model of bacterial evolution. As other MTBC members, M. bovis is postulated as a strictly clonal, slowly evolving pathogen, with apparently no signs of recombination or horizontal gene transfer. In this work, we applied comparative genomics to a whole genome sequence (WGS) dataset composed by 70 M. bovis from different lineages (European and African) to gain insights into the evolutionary forces that shape genetic diversification in M. bovis. Three distinct approaches were used to estimate signs of recombination. Globally, a small number of recombinant events was identified and confirmed by two independent methods with solid support. Still, recombination reveals a weaker effect on M. bovis diversity compared with mutation (overall r/m = 0.037). The differential r/m average values obtained across the clonal complexes of M. bovis in our dataset are consistent with the general notion that the extent of recombination may vary widely among lineages assigned to the same taxonomical species. Based on this work, recombination in M. bovis cannot be excluded and should thus be a topic of further effort in future comparative genomics studies for which WGS of large datasets from different epidemiological scenarios across the world is crucial. A smaller M. bovis dataset (n = 42) from a multi-host TB endemic scenario was then subjected to additional analyses, with the identification of more than 1,800 sites wherein at least one strain showed a single nucleotide polymorphism (SNP). The majority (87.1%) was located in coding regions, with the global ratio of non-synonymous upon synonymous alterations (dN/dS) exceeding 1.5, suggesting that positive selection is an important evolutionary force exerted upon M. bovis. A higher percentage of SNPs was detected in genes enriched into "lipid metabolism", "cell wall and cell processes" and "intermediary metabolism and respiration" functional categories, revealing their underlying importance in M. bovis biology and evolution. A closer look on genes prone to horizontal gene transfer in the MTBC ancestor and included in the 3R (DNA repair, replication and recombination) system revealed a global average negative value for Taijima's D neutrality test, suggesting that past selective sweeps and population expansion after a recent bottleneck remain as major evolutionary drivers of the obligatory pathogen M. bovis in its struggle with the host.

Whole genome sequencing of clinical samples reveals extensively drug resistant tuberculosis (XDR TB) strains from the Beijing lineage in Nigeria, West Africa

Multi-drug (MDR) and extensively drug-resistant (XDR) tuberculosis (TB) continues to be a global public health problem especially in high TB burden countries like Nigeria. Many of these cases are undetected and go on to infect high risk individuals. Clinical samples from positive rifampicin resistant Xpert®MTB/Rif assay were subjected to direct whole genome sequencing and bioinformatics analysis to identify the full antibiotics resistance and lineage profile. We report two (2) XDR TB samples also belonging to the East-Asian/Beijing family of lineage 2 Mycobacterium tuberculosis complex from clinical samples in Nigeria. Our findings further reveal the presence of mutations that confer resistance to first-line drugs (rifampicin, isoniazid, ethambutol and pyrazanimide), second-line injectables (capreomycin, streptomycin, kanamycin and/or amikacin) and at least one of the fluoroquinolones (ofloxacin, moxifloxacin, levofloxacin and/or ciprofloxacin) in both samples. The genomic sequence data from this study not only provide the first evidence of XDR TB in Nigeria and West Africa, but also emphasize the importance of WGS in accurately detecting MDR and XDR TB, to ensure adequate and proper management treatment regimens for affected individuals. This will greatly aid in preventing the spread of drug resistance TB in high burden countries.

Cathepsins and Their Endogenous Inhibitors in Host Defense During Mycobacterium tuberculosis and HIV Infection

The moment a very old bacterial pathogen met a young virus from the 80's defined the beginning of a tragic syndemic for humanity. Such is the case for the causative agent of tuberculosis and the human immunodeficiency virus (HIV). Syndemic is by definition a convergence of more than one disease resulting in magnification of their burden. Both pathogens work synergistically contributing to speed up the replication of each other. Mycobacterium tuberculosis (Mtb) and HIV infections are in the 21st century among the leaders of morbidity and mortality of humankind. There is an urgent need for development of new approaches for prevention, better diagnosis, and new therapies for both infections. Moreover, these approaches should consider Mtb and HIV as a co-infection, rather than just as separate problems, to prevent further aggravation of the HIV-TB syndemic. Both pathogens manipulate the host immune responses to establish chronic infections in intracellular niches of their host cells. This includes manipulation of host relevant antimicrobial proteases such as cathepsins or their endogenous inhibitors. Here we discuss recent understanding on how Mtb and HIV interact with cathepsins and their inhibitors in their multifactorial functions during the pathogenesis of both infections. Particularly we will address the role on pathogen transmission, during establishment of intracellular chronic niches and in granuloma clinical outcome and tuberculosis diagnosis. This area of research will open new avenues for the design of innovative therapies and diagnostic interventions so urgently needed to fight this threat to humanity.

Genetic composition and evolution of the prevalent Mycobacterium tuberculosis lineages 2 and 4 in the Chinese and Zhejiang Province populations

Background

There are seven human-adaptation lineages of Mycobacterium tuberculosis (Mtb). Tuberculosis (TB) dissemination is strongly influenced by human movements and host genetics. The detailed lineage distribution evolution of Mtb in Zhejiang Province is unknown. We aim to determine how different sub-lineages are transmitted and distributed within China and Zhejiang Province.

Methods

We analysed whole-genome sequencing data for a worldwide collection of 1154 isolates and a provincial collection of 1296 isolates, constructed the best-scoring maximum likelihood phylogenetic tree. Bayesian evolutionary analysis was used to calculate the latest common ancestor of lineages 2 and 4. The antigenic diversity of human T cell epitopes was evaluated by calculating the pairwise dN/dS ratios.

Results

Of the Zhejiang isolates, 964 (74.38%) belonged to lineage 2 and 332 (25.62%) belonged to lineage 4. The distributions of the sub-lineages varied across the geographic regions of Zhejiang Province. L2.2 is the most ancient sub-lineage in Zhejiang, first appearing approximately 6897 years ago (95% highest posterior density interval (HDI): 6513–7298). L4.4 is the most modern sub-lineage, first appearing approximately 2217 years ago (95% HDI: 1864–2581). The dN/dS ratios showed that the epitope and non-epitope regions of lineage 2 strains were significantly (P < 0.001) more conserved than those of lineage 4.

Conclusions

An increase in the frequency of lineage 4 may reflect its successful transmission over the last 20 years. The recent common ancestors of the sub-lineages and their transmission routes are relevant to the entry of humans into China and Zhejiang Province. Diversity in T cell epitopes may prevent Mycobacterium tuberculosis from being recognized by the immune system.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-021-00673-7.

The Relevance of Genomic Epidemiology for Control of Tuberculosis in West Africa

Tuberculosis (TB), an airborne infectious disease caused by Mycobacterium tuberculosis complex (MTBC), remains a global health problem. West Africa has a unique epidemiology of TB that is characterized by medium- to high-prevalence. Moreover, the geographical restriction of M. africanum to the sub-region makes West Africa have an extra burden to deal with a two-in-one pathogen. The region is also burdened with low case detection, late reporting, poor treatment adherence leading to development of drug resistance and relapse. Sporadic studies conducted within the subregion report higher burden of drug resistant TB (DRTB) than previously thought. The need for more sensitive and robust tools for routine surveillance as well as to understand the mechanisms of DRTB and transmission dynamics for the design of effective control tools, cannot be overemphasized. The advancement in molecular biology tools including traditional fingerprinting and next generation sequencing (NGS) technologies offer reliable tools for genomic epidemiology. Genomic epidemiology provides in-depth insight of the nature of pathogens, circulating strains and their spread as well as prompt detection of the emergence of new strains. It also offers the opportunity to monitor treatment and evaluate interventions. Furthermore, genomic epidemiology can be used to understand potential emergence and spread of drug resistant strains and resistance mechanisms allowing the design of simple but rapid tools. In this review, we will describe the local epidemiology of MTBC, highlight past and current investigations toward understanding their biology and spread as well as discuss the relevance of genomic epidemiology studies to TB control in West Africa.

Molecular epidemiology and drug susceptibility profiles of Mycobacterium tuberculosis complex isolates from Northern Ghana

OBJECTIVE

We conducted a cross-sectional study in the five administrative regions of Northern Ghana to determine the diversity of Mycobacterium tuberculosis complex (MTBC) sub/lineages and their susceptibility to isoniazid (INH) and rifampicin (RIF).

METHODS

Sputum specimens were collected and cultured from 566 pulmonary tuberculosis patients reporting to 17 health facilities from 2015 to 2019. Mycobacterial isolates obtained from solid cultures were confirmed as members of the MTBC by PCR amplification of IS6110 and rpoß and assigned lineages and sub-lineages using spoligotyping.

RESULTS

Of 294 mycobacterial isolates recovered, MTBC species identified were: M. tuberculosis sensu stricto (Mtbss) 241 (82.0%), M. africanum 41 (13.9%) and M. bovis four (1.4%) with eight (2.7%) unidentified. The human-adapted lineages (L) identified (N=279) were L1 (8/279, 2.9%), L2 (15/279, 5.4%), L3 (7/279, 2.5%), L4 (208/279, 74.5%), L5 (13/279, 4.7%) and L6 (28/279, 10.0%) with three unidentified lineages. Among the 208 L4, the dominant sub-lineages in the region were the Cameroon 120/208 (57.7%) and Ghana 50/208 (24.0%). We found 4.4% (13/294) and 0.7% (2/294) of the patients infected with MTBC isolates resistant to INH only and RIF only, respectively, with 2.4% (7/294) being infected with MDR strains. Whereas L6 was associated with the elderly, we identified that the Ghana sub-lineage of L4 was associated with both INH and MDR (p<0.05), making them important TB pathogens in Northern Ghana and a growing public health concern.

Multiple Merger Genealogies in Outbreaks of Mycobacterium tuberculosis

The Kingman coalescent and its developments are often considered among the most important advances in population genetics of the last decades. Demographic inference based on coalescent theory has been used to reconstruct the population dynamics and evolutionary history of several species, including Mycobacterium tuberculosis (MTB), an important human pathogen causing tuberculosis. One key assumption of the Kingman coalescent is that the number of descendants of different individuals does not vary strongly, and violating this assumption could lead to severe biases caused by model misspecification. Individual lineages of MTB are expected to vary strongly in reproductive success because 1) MTB is potentially under constant selection due to the pressure of the host immune system and of antibiotic treatment, 2) MTB undergoes repeated population bottlenecks when it transmits from one host to the next, and 3) some hosts show much higher transmission rates compared with the average (superspreaders).

Here, we used an approximate Bayesian computation approach to test whether multiple-merger coalescents (MMC), a class of models that allow for large variation in reproductive success among lineages, are more appropriate models to study MTB populations. We considered 11 publicly available whole-genome sequence data sets sampled from local MTB populations and outbreaks and found that MMC had a better fit compared with the Kingman coalescent for 10 of the 11 data sets. These results indicate that the null model for analyzing MTB outbreaks should be reassessed and that past findings based on the Kingman coalescent need to be revisited.

Phylogenomics of Mycobacterium africanum reveals a new lineage and a complex evolutionary history

Human tuberculosis (TB) is caused by members of the Mycobacterium tuberculosis complex (MTBC). The MTBC comprises several human-adapted lineages known as M. tuberculosis sensu stricto, as well as two lineages (L5 and L6) traditionally referred to as Mycobacterium africanum. Strains of L5 and L6 are largely limited to West Africa for reasons unknown, and little is known of their genomic diversity, phylogeography and evolution. Here, we analysed the genomes of 350 L5 and 320 L6 strains, isolated from patients from 21 African countries, plus 5 related genomes that had not been classified into any of the known MTBC lineages. Our population genomic and phylogeographical analyses showed that the unclassified genomes belonged to a new group that we propose to name MTBC lineage 9 (L9). While the most likely ancestral distribution of L9 was predicted to be East Africa, the most likely ancestral distribution for both L5 and L6 was the Eastern part of West Africa. Moreover, we found important differences between L5 and L6 strains with respect to their phylogeographical substructure and genetic diversity. Finally, we could not confirm the previous association of drug-resistance markers with lineage and sublineages. Instead, our results indicate that the association of drug resistance with lineage is most likely driven by sample bias or geography. In conclusion, our study sheds new light onto the genomic diversity and evolutionary history of M. africanum, and highlights the need to consider the particularities of each MTBC lineage for understanding the ecology and epidemiology of TB in Africa and globally.

Human tuberculosis and Mycobacterium tuberculosis complex: A review on genetic diversity, pathogenesis and omics approaches in host biomarkers discovery

Mycobacterium tuberculosis complex (MTBC) refers to a group of mycobacteria encompassing nine members of closely related species that causes tuberculosis in animals and humans. Among the nine members, Mycobacterium tuberculosis (M. tuberculosis) remains the main causative agent for human tuberculosis that results in high mortality and morbidity globally. In general, MTBC species are low in diversity but exhibit distinctive biological differences and phenotypes among different MTBC lineages. MTBC species are likely to have evolved from a common ancestor through insertions/deletions processes resulting in species speciation with different degrees of pathogenicity. The pathogenesis of human tuberculosis is complex and remains poorly understood. It involves multi-interactions or evolutionary co-options between host factors and bacterial determinants for survival of the MTBC. Granuloma formation as a protection or survival mechanism in hosts by MTBC remains controversial. Additionally, MTBC species are capable of modulating host immune response and have adopted several mechanisms to evade from host immune attack in order to survive in humans. On the other hand, current diagnostic tools for human tuberculosis are inadequate and have several shortcomings. Numerous studies have suggested the potential of host biomarkers in early diagnosis of tuberculosis, in disease differentiation and in treatment monitoring. "Multi-omics" approaches provide holistic views to dissect the association of MTBC species with humans and offer great advantages in host biomarkers discovery. Thus, in this review, we seek to understand how the genetic variations in MTBC lead to species speciation with different pathogenicity. Furthermore, we also discuss how the host and bacterial players contribute to the pathogenesis of human tuberculosis. Lastly, we provide an overview of the journey of "omics" approaches in host biomarkers discovery in human tuberculosis and provide some interesting insights on the challenges and directions of "omics" approaches in host biomarkers innovation and clinical implementation.

Variation in Mycobacterium tuberculosis population structure in Iran: a systemic review and meta-analysis

Background

Acquiring comprehensive insight into the dynamics of Mycobacterium tuberculosis (Mtb) population structure is an essential step to adopt effective tuberculosis (TB) control strategies and improve therapeutic methods and vaccines. Accordingly, we performed this systematic review and meta-analysis to determine the overall prevalence of Mtb genotypes/ sublineages in Iran.

Methods

We carried out a comprehensive literature search using the international databases of MEDLINE and Scopus as well as Iranian databases. Articles published until April 2020 were selected based on the PRISMA flow diagram. The overall prevalence of the Mtb genotypes/sublineage in Iran was determined using the random effects or fixed effect model. The metafor R package and MedCalc software were employed for performing this meta-analysis.

Results

We identified 34 studies for inclusion in this study, containing 8329 clinical samples. Based on the pooled prevalence of the Mtb genotypes, NEW1 (21.94, 95% CI: 16.41–28.05%), CAS (19.21, 95% CI: 14.95–23.86%), EAI (12.95, 95% CI: 7.58–19.47%), and T (12.16, 95% CI: 9.18–15.50%) were characterized as the dominant circulating genotypes in Iran. West African (L 5/6), Cameroon, TUR and H37Rv were identified as genotypes with the lowest prevalence in Iran (< 2%). The highest pooled prevalence rates of multidrug-resistant strains were related to Beijing (2.52, 95% CI) and CAS (1.21, 95% CI).

Conclusions

This systematic review showed that Mtb populations are genetically diverse in Iran, and further studies are needed to gain a better insight into the national diversity of Mtb populations and their drug resistance pattern.

Characterization of DNA methylation in Malawian Mycobacterium tuberculosis clinical isolates

Background

Although Mycobacterium tuberculosis (Mtb) strains exhibit genomic homology of >99%, there is considerable variation in the phenotype. The underlying mechanisms of phenotypic heterogeneity in Mtb are not well understood but epigenetic variation is thought to contribute. At present the methylome of Mtb has not been completely characterized.

Methods

We completed methylomes of 18 Mycobacterium tuberculosis (Mtb) clinical isolates from Malawi representing the largest number of Mtb genomes to be completed in a single study using Single Molecule Real Time (SMRT) sequencing to date.

Results

We replicate and confirm four methylation disrupting mutations in 4 lineages of Mtb. For the first time we report complete loss of methylation courtesy of C758T (S253L) mutation in the MamB gene of Indo-oceanic lineage of Mtb. Additionally, we report a novel missense mutation G454A (G152S) in the MamA gene of the Euro-American lineage which could potentially be attributed to total disruption of methylation in the CCCAG motif but partial loss in a partner motif. Through a genomic and methylome comparative analysis with a global sample of sixteen, we report previously unknown mutations affecting the pks15/1 locus in L6 isolates. We confirm that methylation in Mtb is lineage specific although some unresolved issues still remain.

Unexpected diversity of CRISPR unveils some evolutionary patterns of repeated sequences in Mycobacterium tuberculosis

BACKGROUND

Diversity of the CRISPR locus of Mycobacterium tuberculosis complex has been studied since 1997 for molecular epidemiology purposes. By targeting solely the 43 spacers present in the two first sequenced genomes (H37Rv and BCG), it gave a biased idea of CRISPR diversity and ignored diversity in the neighbouring cas-genes.

RESULTS

We set up tailored pipelines to explore the diversity of CRISPR-cas locus in Short Reads. We analyzed data from a representative set of 198 clinical isolates as evidenced by well-characterized SNPs. We found a relatively low diversity in terms of spacers: we recovered only the 68 spacers that had been described in 2000. We found no partial or global inversions in the sequences, letting always the Direct Variant Repeats (DVR) in the same order. In contrast, we found an unexpected diversity in the form of: SNPs in spacers and in Direct Repeats, duplications of various length, and insertions at various locations of the IS6110 insertion sequence, as well as blocks of DVR deletions. The diversity was in part specific to lineages. When reconstructing evolutionary steps of the locus, we found no evidence for SNP reversal. DVR deletions were linked to recombination between IS6110 insertions or between Direct Repeats.

CONCLUSION

This work definitively shows that CRISPR locus of M. tuberculosis did not evolve by classical CRISPR adaptation (incorporation of new spacers) since the last most recent common ancestor of virulent lineages. The evolutionary mechanisms that we discovered could be involved in bacterial adaptation but in a way that remains to be identified.

Investigating extradomiciliary transmission of tuberculosis: An exploratory approach using social network patterns of TB cases and controls and the genotyping of Mycobacterium tuberculosis

Extradomiciliary contacts have been overlooked in the study of TB transmission due to difficulties in identifying actual contacts in large populations. Complex network analysis provides a framework to model the structure of contacts, specially extradomiciliary ones. We conducted a study of incident sputum-positive TB cases and healthy controls occurring in a moderate TB burden city. Cases and controls were interviewed to obtain data regarding the usual locations of residence, work, study, and leisure. Mycobacterium tuberculosis isolated from sputum was genotyped. The collected data were used to build networks based on a framework of putative social interactions indicating possible TB transmission. A user-friendly open source environment (GraphTube) was setup to extract information from the collected data. Networks based on the likelihood of patient-patient, patient-healthy, and healthy-healthy contacts were setup, depending on a constraint of geographical distance of places attended by the volunteers. Using a threshold for the geographical distance of 300 m, the differences between TB cases and controls are revealed. Several clusters formed by social network nodes with high genotypic similarity were characterized. The developed framework provided consistent results and can be used to support the targeted search of potentially infected individuals and to help to understand the TB transmission.

The population genomics of within-host Mycobacterium tuberculosis

Recent progress in genomic sequencing from patient samples has allowed for the first detailed insight into the within-host genetic diversity of Mycobacterium tuberculosis (M.TB), revealing remarkably low levels of variation. While this has often been attributed to low mutation rates, other factors have been described, including resistance evolution (i.e., selective sweeps), widespread purifying and background selection, and, more recently, progeny skew. Here we review recent findings pertaining to the processes governing the evolutionary dynamics of M.TB, discuss their implications for improving our understanding of this important human pathogen, and make recommendations for future work. Significantly, this emerging evolutionary framework involving the joint estimation of demographic, selective, and reproductive processes is forming a new paradigm for the study of within-host pathogen evolution that will be widely applicable across organisms.

'TB or not TB': the conundrum of pre-European contact tuberculosis in the Pacific

Tuberculosis (TB) is a major global health threat, infecting one-third of the world's population. Despite this prominence, the age, origin and spread of the disease have been topics of contentious debate. Molecular studies suggest that Mycobacterium tuberculosis 'sensu stricto', the most common strain of TB infecting humans today, originated in Africa and from there spread into Europe and Asia. The M. tuberculosis strains most commonly found across the Pacific and the Americas today are most closely related to European strains, supporting a hypothesis that the disease only reached these regions relatively recently via European sailors or settlers. However, this hypothesis is inconsistent with palaeopathological evidence of TB-like lesions in human remains from across the Pacific that predate European contact. Similarly, genetic evidence from pre-European South American mummies challenges the notion of a European introduction of the disease into the Pacific. Here, we review the complex evidence for the age and origin of TB in the Pacific, and discuss key gaps in our knowledge and how these may be addressed. This article is part of the theme issue 'Insights into health and disease from ancient biomolecules'.

Cryptic speciation of a pelagic Roseobacter population varying at a few thousand nucleotide sites

A drop of seawater contains numerous microspatial niches at the scale relevant to microbial activities. Examples are abiotic niches such as detrital particles that show different sizes and organic contents, and biotic niches resulting from bacteria-phage and bacteria-phytoplankton interactions. A common practice to investigate the impact of microenvironments on bacterial evolution is to separate the microenvironments physically and compare the bacterial inhabitants from each. It remains poorly understood, however, which microenvironment primarily drives bacterioplankton evolution in the pelagic ocean. By applying a dilution cultivation approach to an undisturbed coastal water sample, we isolate a bacterial population affiliated with the globally dominant Roseobacter group. Although varying at just a few thousand nucleotide sites across the whole genomes, members of this clonal population are diverging into two genetically separated subspecies. Genes underlying speciation are not unique to subspecies but instead clustered at the shared regions that represent ~6% of the genomic DNA. They are primarily involved in vitamin synthesis, motility, oxidative defense, carbohydrate, and amino acid utilization, consistent with the known strategies that roseobacters take to interact with phytoplankton and particles. Physiological assays corroborate that one subspecies outcompetes the other in these traits. Our results indicate that the microenvironments in the pelagic ocean represented by phytoplankton and organic particles are likely important niches that drive the cryptic speciation of the Roseobacter population, though microhabitats contributed by other less abundant pelagic hosts cannot be ruled out.

A seventeenth-century Mycobacterium tuberculosis genome supports a Neolithic emergence of the Mycobacterium tuberculosis complex

BACKGROUND

Although tuberculosis accounts for the highest mortality from a bacterial infection on a global scale, questions persist regarding its origin. One hypothesis based on modern Mycobacterium tuberculosis complex (MTBC) genomes suggests their most recent common ancestor followed human migrations out of Africa approximately 70,000 years before present. However, studies using ancient genomes as calibration points have yielded much younger dates of less than 6000 years. Here, we aim to address this discrepancy through the analysis of the highest-coverage and highest-quality ancient MTBC genome available to date, reconstructed from a calcified lung nodule of Bishop Peder Winstrup of Lund (b. 1605-d. 1679).

RESULTS

A metagenomic approach for taxonomic classification of whole DNA content permitted the identification of abundant DNA belonging to the human host and the MTBC, with few non-TB bacterial taxa comprising the background. Genomic enrichment enabled the reconstruction of a 141-fold coverage M. tuberculosis genome. In utilizing this high-quality, high-coverage seventeenth-century genome as a calibration point for dating the MTBC, we employed multiple Bayesian tree models, including birth-death models, which allowed us to model pathogen population dynamics and data sampling strategies more realistically than those based on the coalescent.

CONCLUSIONS

The results of our metagenomic analysis demonstrate the unique preservation environment calcified nodules provide for DNA. Importantly, we estimate a most recent common ancestor date for the MTBC of between 2190 and 4501 before present and for Lineage 4 of between 929 and 2084 before present using multiple models, confirming a Neolithic emergence for the MTBC.

Machine learning reveals that Mycobacterium tuberculosis genotypes and anatomic disease site impacts drug resistance and disease transmission among patients with proven extra-pulmonary tuberculosis

Background

There is a general dearth of information on extrapulmonary tuberculosis (EPTB). Here, we investigated Mycobacterium tuberculosis (Mtb) drug resistance and transmission patterns in EPTB patients treated in the Tshwane metropolitan area, in South Africa.

Methods

Consecutive Mtb culture-positive non-pulmonary samples from unique EPTB patients underwent mycobacterial genotyping and were assigned to phylogenetic lineages and transmission clusters based on spoligotypes. MTBDRplus assay was used to search mutations for isoniazid and rifampin resistance. Machine learning algorithms were used to identify clinically meaningful patterns in data. We computed odds ratio (OR), attributable risk (AR) and corresponding 95% confidence intervals (CI).

Results

Of the 70 isolates examined, the largest cluster comprised 25 (36%) Mtb strains that belonged to the East Asian lineage. East Asian lineage was significantly more likely to occur within chains of transmission when compared to the Euro-American and East-African Indian lineages: OR = 10.11 (95% CI: 1.56–116). Lymphadenitis, meningitis and cutaneous TB, were significantly more likely to be associated with drug resistance: OR = 12.69 (95% CI: 1.82–141.60) and AR = 0.25 (95% CI: 0.06–0.43) when compared with other EPTB sites, which suggests that poor rifampin penetration might be a contributing factor.

Conclusions

The majority of Mtb strains circulating in the Tshwane metropolis belongs to East Asian, Euro-American and East-African Indian lineages. Each of these are likely to be clustered, suggesting on-going EPTB transmission. Since 25% of the drug resistance was attributable to sanctuary EPTB sites notorious for poor rifampin penetration, we hypothesize that poor anti-tuberculosis drug dosing might have a role in the development of resistance.

The past, present and future of ancient bacterial DNA

Groundbreaking studies conducted in the mid-1980s demonstrated the possibility of sequencing ancient DNA (aDNA), which has allowed us to answer fundamental questions about the human past. Microbiologists were thus given a powerful tool to glimpse directly into inscrutable bacterial history, hitherto inaccessible due to a poor fossil record. Initially plagued by concerns regarding contamination, the field has grown alongside technical progress, with the advent of high-throughput sequencing being a breakthrough in sequence output and authentication. Albeit burdened with challenges unique to the analysis of bacteria, a growing number of viable sources for aDNA has opened multiple avenues of microbial research. Ancient pathogens have been extracted from bones, dental pulp, mummies and historical medical specimens and have answered focal historical questions such as identifying the aetiological agent of the black death as Yersinia pestis. Furthermore, ancient human microbiomes from fossilized faeces, mummies and dental plaque have shown shifts in human commensals through the Neolithic demographic transition and industrial revolution, whereas environmental isolates stemming from permafrost samples have revealed signs of ancient antimicrobial resistance. Culminating in an ever-growing repertoire of ancient genomes, the quickly expanding body of bacterial aDNA studies has also enabled comparisons of ancient genomes to their extant counterparts, illuminating the evolutionary history of bacteria. In this review we summarize the present avenues of research and contextualize them in the past of the field whilst also pointing towards questions still to be answered.

The geno-spatio analysis of Mycobacterium tuberculosis complex in hot and cold spots of Guangxi, China

BACKGROUND

At present, there are few studies on polymorphism of Mycobacterium tuberculosis (Mtb) gene and how it affects the TB epidemic. This study aimed to document the differences of polymorphisms between tuberculosis hot and cold spot areas of Guangxi Zhuang Autonomous Region, China.

METHODS

The cold and hot spot areas, each with 3 counties, had been pre-identified by TB incidence for 5 years from the surveillance database. Whole genome sequencing analysis was performed on all sputum Mtb isolates from the detected cases during January and June 2018. Single nucleotide polymorphism (SNP) of each isolate compared to the H37Rv strain were called and used for lineage and sub-lineage identification. Pairwise SNP differences between every pair of isolates were computed. Analyses of Molecular Variance (AMOVA) across counties of the same hot or cold spot area and between the two areas were performed.

RESULTS

As a whole, 59.8% (57.7% sub-lineage 2.2 and 2.1% sub-lineage 2.1) and 39.8% (17.8% sub-lineage 4.4, 6.5% sub-lineage 4.2 and 15.5% sub-lineage 4.5) of the Mtb strains were Lineage 2 and Lineage 4 respectively. The percentages of sub-lineage 2.2 (Beijing family strains) are significantly higher in hot spots. Through the MDS dimension reduction, the genomic population structure in the three hot spot counties is significantly different from those three cold spot counties (T-test p = 0.05). The median of SNPs distances among Mtb isolates in cold spots was greater than that in hot spots (897 vs 746, Rank-sum test p < 0.001). Three genomic clusters, each with genomic distance ≤12 SNPs, were identified with 2, 3 and 4 consanguineous strains. Two clusters were from hot spots and one was from cold spots.

CONCLUSION

Narrower genotype diversity in the hot area may indicate higher transmissibility of the Mtb strains in the area compared to those in the cold spot area.

A sister lineage of the Mycobacterium tuberculosis complex discovered in the African Great Lakes region

The human- and animal-adapted lineages of the Mycobacterium tuberculosis complex (MTBC) are thought to have expanded from a common progenitor in Africa. However, the molecular events that accompanied this emergence remain largely unknown. Here, we describe two MTBC strains isolated from patients with multidrug resistant tuberculosis, representing an as-yet-unknown lineage, named Lineage 8 (L8), seemingly restricted to the African Great Lakes region. Using genome-based phylogenetic reconstruction, we show that L8 is a sister clade to the known MTBC lineages. Comparison with other complete mycobacterial genomes indicate that the divergence of L8 preceded the loss of the cobF genome region - involved in the cobalamin/vitamin B12 synthesis - and gene interruptions in a subsequent common ancestor shared by all other known MTBC lineages. This discovery further supports an East African origin for the MTBC and provides additional molecular clues on the ancestral genome reduction associated with adaptation to a pathogenic lifestyle.

The human- and animal-adapted lineages of the Mycobacterium tuberculosis complex (MTBC) are thought to be evolved from a common progenitor in Africa. Here, the authors identify two MTBC strains isolated from patients with multidrug-resistant tuberculosis, representing an as-yet-unknown lineage further supporting an East African origin for the MTBC.

Homoplastic single nucleotide polymorphisms contributed to phenotypic diversity in Mycobacterium tuberculosis

Homoplastic mutations are mutations independently occurring in different clades of an organism. The homoplastic changes may be a result of convergence evolution due to selective pressures. Reports on the analysis of homoplastic mutations in Mycobacterium tuberculosis have been limited. Here we characterized the distribution of homoplastic single nucleotide polymorphisms (SNPs) among genomes of 1,170 clinical M. tuberculosis isolates. They were present in all functional categories of genes, with pe/ppe gene family having the highest ratio of homoplastic SNPs compared to the total SNPs identified in the same functional category. Among the pe/ppe genes, the homoplastic SNPs were common in a relatively small number of homologous genes, including ppe18, the protein of which is a component of a promising candidate vaccine, M72/AS01E. The homoplastic SNPs in ppe18 were particularly common among M. tuberculosis Lineage 1 isolates, suggesting the need for caution in extrapolating the results of the vaccine trial to the population where L1 is endemic in Asia. As expected, homoplastic SNPs strongly associated with drug resistance. Most of these mutations are already well known. However, a number of novel mutations associated with streptomycin resistance were identified, which warrants further investigation. A SNP in the intergenic region upstream of Rv0079 (DATIN) was experimentally shown to increase transcriptional activity of the downstream gene, suggesting that intergenic homoplastic SNPs should have effects on the physiology of the bacterial cells. Our study highlights the potential of homoplastic mutations to produce phenotypic changes. Under selective pressure and during interaction with the host, homoplastic mutations may confer advantages to M. tuberculosis and deserve further characterization.