TbD1 deletion as a driver of the evolutionary success of modern epidemic Mycobacterium tuberculosis lineages

Analysis of complete genomes of Mycobacterium tuberculosis sublineage 2.1 (Proto-Beijing) revealed the presence of three pe_pgrs3-pe_pgrs4-like genes

Mycobacterium tuberculosis Complex (MTBC), the etiological agent of tuberculosis (TB), demonstrates considerable genotypic diversity with distinct geographic distributions and variable virulence profiles. The pe-ppe gene family is especially noteworthy for its extensive variability and roles in host immune response modulation and virulence enhancement. We sequenced an Mtb genotype L2.1 isolate from Chiangrai, Northern Thailand, using second and third-generation sequencing technologies. Comparative genomic analysis with two additional L2.1 isolates and two L2.2.AA3 (Asia Ancestral 3 Beijing) isolates revealed significant pe-ppe gene variations. Notably, all L2.1 isolates harbored three copies of pe_pgrs3-pe_pgrs4-like genes (pe_pgrs3*, pe_pgrs4*, and pe_pgrs4), different from L2.2.AA3 and H37Rv strains. Additionally, ppe53 was duplicated in all but H37Rv, and ppe50 was deleted in L2.1 isolates, contrasting with an extended ppe50 in an L2.2 isolate (Mtb 18b), which contains an additional SVP motif. Complete deletion of ppe66 and loss of wag22 were observed in L2.1 isolates. These findings highlight the high structural variability of the pe-ppe gene family, emphasizing its complex roles in Mtb-host immune interactions. This genetic complexity offers potentially critical insights into mycobacterial pathogenesis, with significant implications for vaccine development and diagnostics.

Comprehensive analysis of Mycobacterium tuberculosis genomes reveals genetic variations in bacterial virulence

Tuberculosis, a disease caused by Mycobacterium tuberculosis (Mtb), is a significant health problem worldwide. Here, we developed a method to detect large insertions and deletions (indels), which have been generally understudied. Leveraging this framework, we performed a comprehensive analysis of single nucleotide variants and small and large indels across 1,960 Mtb clinical isolates. Comparing the distribution of variants demonstrated that gene disruptive variants are underrepresented in genes essential for bacterial survival. An evolutionary analysis revealed that Mtb genomes are enriched in partially deleterious mutations. Genome-wide association studies identified small and large deletions in eccB2 significantly associated with patient prognosis. Additionally, we unveil significant associations with antibiotic resistance in 23 non-canonical genes. Among these, large indels are primarily found in genetic regions of Rv1216c, Rv1217c, fadD11, and ctpD. This study provides a comprehensive catalog of genetic variations and highlights their potential impact for the future treatment and risk prediction of tuberculosis.

Genomic Interactions Between Mycobacterium tuberculosis and Humans

Mycobacterium tuberculosis is considered by many to be the deadliest microbe, with the estimated annual cases numbering more than 10 million. The bacteria, including Mycobacterium africanum, are classified into nine major lineages and hundreds of sublineages, each with different geographical distributions and levels of virulence. The phylogeographic patterns can be a result of recent and early human migrations as well as coevolution between the bacteria and various human populations, which may explain why many studies on human genetic factors contributing to tuberculosis have not been replicable in different areas. Moreover, several studies have revealed the significance of interactions between human genetic variations and bacterial genotypes in determining the development of tuberculosis, suggesting coadaptation. The increased availability of whole-genome sequence data from both humans and bacteria has enabled a better understanding of these interactions, which can inform the development of vaccines and other control measures.

Incipient tuberculosis: a comprehensive overview

In the context of the evolving global health landscape shaped by the COVID-19 pandemic, tuberculosis (TB) is gaining renewed attention as a reemerging threat even in low-endemic countries. Immunological tests such as the tuberculin skin test (TST) and interferon-gamma release assay (IGRA) are pivotal in identifying tuberculosis infection (TBI). However, their inability to distinguish between past and ongoing infection poses a diagnostic challenge, possibly leading to the unnecessary treatment of a significant portion of the population with potential side effects. This review delves into the concept of incipient tuberculosis (ITB), a dynamic, presymptomatic stage characterized by heightened Mycobacterium tuberculosis complex (MTC) metabolic activity and replication that result in minimal radiological changes, signifying a transitional state between TBI and TB. Key focus areas include epidemiological factors, underlying pathogenesis, imaging findings, and the ongoing challenges in the identification of individuals with ITB through the development of new biomarkers and the use of whole-genome sequencing-based analyses to implement early treatment strategies.

Analysis of the limited M. tuberculosis accessory genome reveals potential pitfalls of pan-genome analysis approaches

Pan-genome analysis is a fundamental tool for studying bacterial genome evolution; however, the variety of methods used to define and measure the pan-genome poses challenges to the interpretation and reliability of results. To quantify sources of bias and error related to common pan-genome analysis approaches, we evaluated different approaches applied to curated collection of 151 Mycobacterium tuberculosis ( Mtb ) isolates. Mtb is characterized by its clonal evolution, absence of horizontal gene transfer, and limited accessory genome, making it an ideal test case for this study. Using a state-of-the-art graph-genome approach, we found that a majority of the structural variation observed in Mtb originates from rearrangement, deletion, and duplication of redundant nucleotide sequences. In contrast, we found that pan-genome analyses that focus on comparison of coding sequences (at the amino acid level) can yield surprisingly variable results, driven by differences in assembly quality and the softwares used. Upon closer inspection, we found that coding sequence annotation discrepancies were a major contributor to inflated Mtb accessory genome estimates. To address this, we developed panqc, a software that detects annotation discrepancies and collapses nucleotide redundancy in pan-genome estimates. When applied to Mtb and E. coli pan-genomes, panqc exposed distinct biases influenced by the genomic diversity of the population studied. Our findings underscore the need for careful methodological selection and quality control to accurately map the evolutionary dynamics of a bacterial species.

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.

The association between phylogenetic lineage and the subclinical phenotype of pulmonary tuberculosis: A retrospective 2-cohort study

BACKGROUND

Subclinical pulmonary tuberculosis (PTB) is an asymptomatic disease state between established TB infection and symptomatic (clinical) TB disease. It is present in 20-25% of PTB patients in high-income countries. Mycobacterium tuberculosis complex (MTBC) genetic heterogeneity, and differential host immunological responses, have been implicated in its pathogenesis.

METHODS

To determine the association between MTBC lineage and PTB disease phenotype, we used two retrospective cohorts of PTB patients in Canada and two independent lineage attribution methods (DNA fingerprinting and genome sequencing). The first cohort, Cohort 1, consisted of consecutively diagnosed PTB patients between 2014 and 2020. The second, Cohort 2, consisted of newly-arrived foreign-born PTB patients who either were or were not referred for post-landing medical surveillance between 2004 and 2017. Univariable and multivariable logistic regression models were sequentially fitted to both cohorts, adjusting for age, sex, disease type, drug resistance and HIV. Evolution of radiographic features was correlated to lineage in Cohort 2.

FINDINGS

Cohort 1 and 2 included 874 (209 subclinical) and 111 (44 subclinical) patients, respectively. In both cohorts, subclinical patients were more likely than clinical patients to have relapse/retreatment disease, be smear-negative, have longer times-to-culture positivity and to harbor an ancestral MTBC lineage (Indo-Oceanic or Mycobacterium africanum). Relapse/retreatment disease and ancestral MTBC lineage were independent predictors of subclinical disease (ORs and 95% CIs in Cohort 1, 1.85 [1.07,3.28], p < 0.029 and 2.30 [1.66,3.18], p < 0.001, respectively, and Cohort 2, 5.74 [1.37-24.06], p < 0.017 and 3.21 (1.29,7.97], p < 0.012, respectively). The geographic distribution of Indo-Oceanic strains causing subclinical disease was uneven. Non-progressive lung disease was more common in patients infected with ancestral than modern lineages in Cohort 2, 56.0% vs 25.4%, p < 0.005.

INTERPRETATION

MTBC lineage is a strong predictor of PTB disease phenotype. The genetic drivers of this association, and the relative contribution of other explanatory variables, are unknown.

Efflux pump Rv1258c activates novel functions of the oxidative stress and via the VII secretion system ESX-3-mediated iron metabolic pathway in Mycobacterium tuberculosis

Oxidative stress and iron metabolism are essential for Mycobacterium tuberculosis (M.tb) survival in host cells. The efflux pump Rv1258c belongs to the major facilitator superfamily (MFS) and can actively pump drugs to promote certain drug resistance in M.tb. In this study, we compared H37RvΔRv1258c with wild-type (WT) M.tb H37Rv. The qRT-PCR results suggested that Rv1258c is potentially involved in the iron metabolic pathway by regulating the expression of ESX-3, which is required for iron uptake. Protein-Protein Affinity Predictor (PPA-Pred2) and the artificial intelligence program AlphaFold 2 were used for prediction and showed that Rv1258c has direct interactions with PPE4 and EccD3 but weak interactions with EccB3. This was further determined via protein-protein interaction analysis of the yeast two-hybrid expression system. By comparing mutant H37RvΔRv1258c strains with WT strains, we discovered that the absence of Rv1258c led to elevated intracellular H+ potential and NAD+/NADH ratios in M.tb, thereby resulting in oxidative stress. We hypothesize that the efflux pump Rv1258c not only has the function of regulating drug resistance in M.tb but also has a novel function in activating oxidative stress and regulating ESX-3-associated iron metabolism in M.tb.

An imputed ancestral reference genome for the Mycobacterium tuberculosis complex better captures structural genomic diversity for reference-based alignment workflows

Reference-based alignment of short-reads is a widely used technique in genomic analysis of the Mycobacterium tuberculosis complex (MTBC) and the choice of reference sequence impacts the interpretation of analyses. The most widely used reference genomes include the ATCC type strain (H37Rv) and the putative MTBC ancestral sequence of Comas et al. both of which are based on a lineage 4 sequence. As such, these reference sequences do not capture all of the structural variation known to be present in the ancestor of the MTBC. To better represent the base of the MTBC, we generated an imputed ancestral genomic sequence, termed MTBC0 from reference-free alignments of closed MTBC genomes. When used as a reference sequence in alignment workflows, MTBC0 mapped more short sequencing reads and called more pairwise SNPs relative to the Comas et al. sequence while exhibiting minimal impact on the overall phylogeny of MTBC. The results also show that MTBC0 provides greater fidelity in capturing genomic variation and allows for the inclusion of regions absent from H37Rv in standard MTBC workflows without additional steps. The use of MTBC0 as an ancestral reference sequence in standard workflows modestly improved read mapping, SNP calling and intuitively facilitates the study of structural variation and evolution in MTBC.

Dynamics of tuberculosis infection in various populations during the 19th and 20th century: The impact of conservative and pharmaceutical treatments

Humans and Mycobacterium tuberculosis have co-evolved together for thousands of years. Many individuals are infected with the bacterium, but few show signs and symptoms of tuberculosis (TB). Pharmacotherapy to treat those who develop disease is useful, but drug resistance and non-adherence significantly impact the efficacy of these treatments. Prior to the introduction of antibiotic therapies, public health strategies were used to reduce TB mortality. This work shows how these strategies were able to reduce TB mortality in 19th and 20th century populations, compared with antibiotic treatments. Previously published mortality data from historical records for several populations (Switzerland, Germany, England and Wales, Scotland, USA, Japan, Brazil and South Africa) were used. Curvilinear regression was used to examine the reduction in mortality before and after the introduction of antibiotic treatments (1946). A strong decline in TB mortality was already occurring in Switzerland, Germany, England and Wales, Scotland and the USA prior to the introduction of antibiotic treatment. This occurred following many public health interventions including improved sanitation, compulsory reporting of TB cases, diagnostic techniques and sanatoria treatments. Following the introduction of antibiotics, mortality rates declined further, however, this had a smaller effect than the previously employed strategies. In Japan, Brazil and South Africa, reductions in mortality rates were largely driven by antibiotic treatments that caused rapid decline of mortality, with a smaller contribution from public health strategies. For the development of active disease, immune status is important. Individuals infected with the bacterium are more likely to develop signs and symptoms if their immune function is reduced. Effective strategies against TB can therefore include enhancing immune function of the population by improving nutrition, as well as reducing transmission by improving living conditions and public health. This has been effective in the past. Improving immunity may be an important strategy against drug resistant TB.

Evaluating the Sensitivity of Different Molecular Techniques for Detecting Mycobacterium tuberculosis Complex in Patients with Pulmonary Infection

This study aimed to evaluate the accuracy of detecting drug-resistant Mycobacterium tuberculosis complex (MTBC)-specific DNA in sputum specimens from 48 patients diagnosed with pulmonary tuberculosis. The presence of MTBC DNA in the specimens was validated using the GeneXpert MTB/RIF system and compared with a specific PCR assay targeting the IS6110 and the mtp40 gene sequence fragments. Additionally, the results obtained by multiplex PCR assays to detect the most frequently encountered rifampin, isoniazid, and ethambutol resistance-conferring mutations were matched with those obtained by GeneXpert and phenotypic culture-based drug susceptibility tests. Of the 48 sputum samples, 25 were positive for MTBC using the GeneXpert MTB/RIF test. Nevertheless, the IS6110 and mtp40 single-step PCR revealed the IS6110 in 27 of the 48 sputum samples, while the mtp40 gene fragment was found in only 17 of them. Furthermore, multiplex PCR assays detected drug-resistant conferring mutations in 21 (77.8%) of the 27 samples with confirmed MTBC DNA, 10 of which contained single drug-resistant conferring mutations towards ethambutol and two towards rifampin, and the remaining nine contained double-resistant mutations for ethambutol and rifampin. In contrast, only five sputum specimens (18.5%) contained drug-resistant MTBC isolates, and two contained mono-drug-resistant MTBC species toward ethambutol and rifampin, respectively, and the remaining three were designated as multi-drug resistant toward both drugs using GeneXpert and phenotypic culture-based drug susceptibility tests. Such discrepancies in the results emphasize the need to develop novel molecular tests that associate with phenotypic non-DNA-based assays to improve the detection of drug-resistant isolates in clinical specimens in future studies.

Tuberculosis Severity Predictive Model Using Mtb Variants and Serum Biomarkers in a Colombian Cohort of APTB Patients

Currently, tuberculosis (TB) is a bacterial infection caused by Mycobacterium tuberculosis (Mtb) that primarily affects the lungs. The severity of active pulmonary TB (APTB) is an important determinant of transmission, morbidity, mortality, disease experience, and treatment outcomes. Several publications have shown a high prevalence of disabling complications in individuals who have had severe APTB. Furthermore, certain strains of Mtb were associated with more severe disease outcomes. The use of biomarkers to predict severe APTB patients who are candidates for host-directed therapies, due to the high risk of developing post-tuberculous lung disease (PTLD), has not yet been implemented in the management of TB patients. We followed 108 individuals with APTB for 6 months using clinical tools, flow cytometry, and whole-genome sequencing (WGS). The median age of the study population was 26.5 years, and the frequency of women was 53.7%. In this study, we aimed to identify biomarkers that could help us to recognize individuals with APTB and improve our understanding of the immunopathology in these individuals. In this study, we conducted a follow-up on the treatment progress of 121 cases of APTB. The follow-up process commenced at the time of diagnosis (T0), continued with a control visit at 2 months (T2), and culminated in an exit appointment at 6 months following the completion of medical treatment (T6). People classified with severe APTB showed significantly higher levels of IL-6 (14.7 pg/mL; p < 0.05) compared to those with mild APTB (7.7 pg/mL) at T0. The AUCs for the ROC curves and the Matthews correlation coefficient values (MCC) demonstrate correlations ranging from moderate to very strong. We conducted WGS on 88 clinical isolates of Mtb, and our analysis revealed a total of 325 genes with insertions and deletions (Indels) within their coding regions when compared to the Mtb H37Rv reference genome. The pattern of association was found between serum levels of CHIT1 and the presence of Indels in Mtb isolates from patients with severe APTB. A key finding in our study was the high levels of CHIT1 in severe APTB patients. We identified a biomarker profile (IL-6, IFN-γ, IL-33, and CHIT1) that allows us to identify individuals with severe APTB, as well as the identification of a panel of polymorphisms (125) in clinical isolates of Mtb from individuals with severe APTB. Integrating these findings into a predictive model of severity would show promise for the management of APTB patients in the future, to guide host-directed therapy and reduce the prevalence of PTLD.

Whole genome sequencing-based identification of human tuberculosis caused by animal-lineage Mycobacterium orygis

A recently described member of the Mycobacterium tuberculosis complex (MTBC) is Mycobacterium orygis, which can cause disease primarily in animals but also in humans. Although M. orygis has been reported from different geographic regions around the world, due to a lack of proper identification techniques, the contribution of this emerging pathogen to the global burden of zoonotic tuberculosis is not fully understood. In the present work, we report single nucleotide polymorphism (SNP) analysis using whole genome sequencing (WGS) that can accurately identify M. orygis and differentiate it from other members of the MTBC species. WGS-based SNP analysis was performed for 61 isolates from different provinces in Canada that were identified as M. orygis. A total of 56 M. orygis sequences from the public databases were also included in the analysis. Several unique SNPs in the gyrB, PPE55, Rv2042c, leuS, mmpL6, and mmpS6 genes were used to determine their effectiveness as genetic markers for the identification of M. orygis. To the best of our knowledge, five of these SNPs, viz., gyrB 277 (A→G), gyrB 1478 (T→C), leuS 1064 (A→T), mmpL6 486 (T→C), and mmpS6 334 (C→G), are reported for the first time in this study. Our results also revealed several SNPs specific to other species within MTBC. The phylogenetic analysis shows that the studied genomes were genetically diverse and clustered with M. orygis sequences of human and animal origin reported from different geographic locations. Therefore, the present study provides a new insight into the high-confidence identification of M. orygis from MTBC species based on WGS data, which can be useful for reference and diagnostic laboratories.

Bradyrhizobium ottawaense efficiently reduces nitrous oxide through high nosZ gene expression

N2O is an important greenhouse gas influencing global warming, and agricultural land is the predominant (anthropogenic) source of N2O emissions. Here, we report the high N2O-reducing activity of Bradyrhizobium ottawaense, suggesting the potential for efficiently mitigating N2O emission from agricultural lands. Among the 15 B. ottawaense isolates examined, the N2O-reducing activities of most (13) strains were approximately five-fold higher than that of Bradyrhizobium diazoefficiens USDA110T under anaerobic conditions. This robust N2O-reducing activity of B. ottawaense was confirmed by N2O reductase (NosZ) protein levels and by mitigation of N2O emitted by nodule decomposition in laboratory system. While the NosZ of B. ottawaense and B. diazoefficiens showed high homology, nosZ gene expression in B. ottawaense was over 150-fold higher than that in B. diazoefficiens USDA110T, suggesting the high N2O-reducing activity of B. ottawaense is achieved by high nos expression. Furthermore, we examined the nos operon transcription start sites and found that, unlike B. diazoefficiens, B. ottawaense has two transcription start sites under N2O-respiring conditions, which may contribute to the high nosZ expression. Our study indicates the potential of B. ottawaense for effective N2O reduction and unique regulation of nos gene expression towards the high performance of N2O mitigation in the soil.

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.

Molecular investigations of Mycobacterium tuberculosis genotypes among baseline and follow-up strains circulating in four regions of Eswatini

BACKGROUND

The tuberculosis (TB) epidemic remains a major global health problem and Eswatini is not excluded. Our study investigated the circulating genotypes in Eswatini and compared them at baseline (start of treatment) and follow-up during TB treatment.

METHODS

Three hundred and ninety (n = 390) participants were prospectively enrolled from referral clinics and patients who met the inclusion criteria, were included in the study. A total of 103 participants provided specimens at baseline and follow-up within six months. Mycobacterium tuberculosis (M.tb) strains were detected by GeneXpert® MTB/RIF assay (Cephied, USA) and Ziehl -Neelsen (ZN) microscopy respectively at baseline and follow-up time-points respectively. The 206 collected specimens were decontaminated and cultured on BACTEC™ MGIT™ 960 Mycobacteria Culture System (Becton Dickinson, USA). Drug sensitivity testing was performed at both baseline and follow-up time points. Spoligotyping was performed on both baseline and follow-up strains after DNA extraction.

RESULTS

Resistance to at least one first line drug was detected higher at baseline compared to follow-up specimens with most of them developing into multidrug-resistant (MDR)-TB. A total of four lineages and twenty genotypes were detected. The distribution of the lineages varied among the different regions in Eswatini. The Euro-American lineage was the most prevalent with 46.12% (95/206) followed by the East Asian with 24.27% (50/206); Indo-Oceanic at 9.71% (20/206) and Central Asian at 1.94% (4/206). Furthermore, a high proportion of the Beijing genotype at 24.27% (50/206) and S genotype at 16.50% (34/206) were detected. The Beijing genotype was predominant in follow-up specimens collected from the Manzini region with 48.9% (23/47) (p = 0.001). A significant proportion of follow-up specimens developed MDR-TB (p = 0.001) with Beijing being the major genotype in most follow-up specimens (p < 0.000).

CONCLUSION

Eswatini has a high M.tb genotypic diversity. A significant proportion of the TB infected participants had the Beijing genotype associated with MDR-TB in follow-up specimens and thus indicate community wide transmission.

Ion mobility mass spectrometry for the study of mycobacterial mycolic acids

Lipids are highly structurally diverse molecules involved in a wide variety of biological processes. The involvement of lipids is even more pronounced in mycobacteria, including the human pathogen Mycobacterium tuberculosis, which produces a highly complex and diverse set of lipids in the cell envelope. These lipids include mycolic acids, which are among the longest fatty acids in nature and can contain up to 90 carbon atoms. Mycolic acids are ubiquitously found in mycobacteria and are alpha branched and beta hydroxylated lipids. Discrete modifications, such as alpha, alpha', epoxy, methoxy, keto, and carboxy, characterize mycolic acids at the species level. Here, we used high precision ion mobility-mass spectrometry to build a database including 206 mass-resolved collision cross sections (CCSs) of mycolic acids originating from the strict human pathogen M. tuberculosis, the opportunistic strains M. abscessus, M. marinum and M. avium, and the nonpathogenic strain M. smegmatis. Primary differences between the mycolic acid profiles could be observed between mycobacterial species. Acyl tail length and modifications were the primary structural descriptors determining CCS magnitude. As a resource for researchers, this work provides a detailed catalogue of the mass-resolved collision cross sections for mycolic acids along with a workflow to generate and analyse the dataset generated.

Genomic changes underpinning the emergence of a successful Mycobacterium tuberculosis Latin American and Mediterranean clonal complex

Introduction

The Latin American and Mediterranean sublineage (L4.3/LAM) is the most common generalist sublineage of Mycobacterium tuberculosis lineage 4 (L4), yet certain L4.3/LAM genotypes appear to be confined to particular geographic regions. This is typically the case of a L4.3/LAM clonal complex (CC), TUN4.3_CC1, which is the most preponderant in Tunisia (61.5% of L4.3/LAM).

Methods

Here, we used whole-genome sequencing data of 346 globally distributed L4 clinical strains, including 278 L4.3/LAM isolates, to reconstruct the evolutionary history of TUN4.3_CC1 and delineate critical genomic changes underpinning its success.

Results and Discussion

Phylogenomic coupled to phylogeographic analyses indicated that TUN4.3_CC1 has evolved locally, being confined mainly to North Africa. Maximum likelihood analyses using the site and branch-site models of the PAML package disclosed strong evidence of positive selection in the gene category "cell wall and cell processes" of TUN4.3_CC1. Collectively, the data indicate that TUN4.3_CC1 has inherited several mutations, which could have potentially contributed to its evolutionary success. Of particular interest are amino acid replacements at the esxK and eccC2 genes of the ESX/Type VII secretion system, which were found to be specific to TUN4.3_CC1, being common to almost all isolates. Because of its homoplastic nature, the esxK mutation could potentially have endowed TUN4.3_CC1 with a selective advantage. Moreover, we noticed the occurrence of additional, previously described homoplasic nonsense mutations in ponA1 and Rv0197. The mutation in the latter gene, a putative oxido-reductase, has previously been shown to be correlated with enhanced transmissibility in vivo. In sum, our findings unveiled several features underpinning the success of a locally evolved L4.3/LAM clonal complex, lending further support to the critical role of genes encoded by the ESX/type VII secretion system.

Pangenome Reconstruction of Mycobacterium tuberculosis as a Guide to Reveal Genomic Features Associated with Strain Clinical Phenotype

Tuberculosis (TB) is one of the leading causes of human deaths worldwide caused by infectious diseases. TB infection by Mycobacterium tuberculosis can occur in the lungs, causing pulmonary tuberculosis (PTB), or in any other organ of the body, resulting in extrapulmonary tuberculosis (EPTB). There is no consensus on the genetic determinants of this pathogen that may contribute to EPTB. In this study, we constructed the M. tuberculosis pangenome and used it as a tool to seek genomic signatures associated with the clinical presentation of TB based on its accessory genome differences. The analysis carried out in the present study includes the raw reads of 490 M. tuberculosis genomes (PTB n = 245, EPTB n = 245) retrieved from public databases that were assembled, as well as ten genomes from Mexican strains (PTB n = 5, EPTB n = 5) that were sequenced and assembled. All genomes were annotated and then used to construct the pangenome with Roary and Panaroo. The pangenome obtained using Roary consisted of 2231 core genes and 3729 accessory genes. On the other hand, the pangenome resulting from Panaroo consisted of 2130 core genes and 5598 accessory genes. Associations between the distribution of accessory genes and the PTB/EPTB phenotypes were examined using the Scoary and Pyseer tools. Both tools found a significant association between the hspR, plcD, Rv2550c, pe_pgrs5, pe_pgrs25, and pe_pgrs57 genes and the PTB genotype. In contrast, the deletion of the aceA, esxR, plcA, and ppe50 genes was significantly associated with the EPTB phenotype. Rv1759c and Rv3740 were found to be associated with the PTB phenotype according to Scoary; however, these associations were not observed when using Pyseer. The robustness of the constructed pangenome and the gene-phenotype associations is supported by several factors, including the analysis of a large number of genomes, the inclusion of the same number of PTB/EPTB genomes, and the reproducibility of results thanks to the different bioinformatic tools used. Such characteristics surpass most of previous M. tuberculosis pangenomes. Thus, it can be inferred that the deletion of these genes can lead to changes in the processes involved in stress response and fatty acid metabolism, conferring phenotypic advantages associated with pulmonary or extrapulmonary presentation of TB. This study represents the first attempt to use the pangenome to seek gene-phenotype associations in M. tuberculosis.

Understanding the Phage-Host Interaction Mechanism toward Improving the Efficacy of Current Antibiotics in Mycobacterium abscessus

Pulmonary infections caused by Mycobacterium abscessus (MAB) have been increasing in incidence in recent years, leading to chronic and many times fatal infections due to MAB's natural resistance to most available antimicrobials. The use of bacteriophages (phages) in clinics is emerging as a novel treatment strategy to save the lives of patients suffering from drug-resistant, chronic, and disseminated infections. The substantial research indicates that phage-antibiotic combination therapy can display synergy and be clinically more effective than phage therapy alone. However, there is limited knowledge in the understanding of the molecular mechanisms in phage-mycobacteria interaction and the synergism of phage-antibiotic combinations. We generated the lytic mycobacteriophage library and studied phage specificity and the host range in MAB clinical isolates and characterized the phage's ability to lyse the pathogen under various environmental and mammalian host stress conditions. Our results indicate that phage lytic efficiency is altered by environmental conditions, especially in conditions of biofilm and intracellular states of MAB. By utilizing the MAB gene knockout mutants of the MAB_0937c/MmpL10 drug efflux pump and MAB_0939/pks polyketide synthase enzyme, we discovered the surface glycolipid diacyltrehalose/polyacyltrehalose (DAT/PAT) as one of the major primary phage receptors in mycobacteria. We also established a set of phages that alter the MmpL10 multidrug efflux pump function in MAB through an evolutionary trade-off mechanism. The combination of these phages with antibiotics significantly decreases the number of viable bacteria when compared to phage or antibiotic-alone treatments. This study deepens our understanding of phage-mycobacteria interaction mechanisms and identifies therapeutic phages that can lower bacterial fitness by impairing an antibiotic efflux function and attenuating the MAB intrinsic resistance mechanism via targeted therapy.

Performance of lipid fingerprint by routine matrix-assisted laser desorption/ionization time of flight for the diagnosis of Mycobacterium tuberculosis complex species

OBJECTIVES

Rapid detection of bacterial pathogens at species and sub-species levels is crucial for appropriate treatment, infection control, and public health management. Currently, one of the challenges in clinical microbiology is the discrimination of mycobacterial sub-species within the M. tuberculosis complex (MTBC). Our objective was to evaluate the ability of a biosafe mycobacterial lipid-based approach to identify MTBC cultures and sub-species.

METHODS

A blinded study was conducted using 90 mycobacterial clinical isolate strains comprising MTBC strains sub-cultured in Middlebrook 7H11 medium supplemented with 10% oleic-acid, dextrose, catalase growth supplement and incubated for up to 6 weeks at 37°C and using the following seven reference strains (M. tuberculosis H37Rv, M canettii, M. africanum, M. pinnipedii, M. caprae, M. bovis, and M. bovis BCG) grown under the same conditions, to set the reference lipid database and test it against the 90 MTBC clinical isolates. Cultured mycobacteria were heat-inactivated and loaded onto the matrix-assisted laser desorption/ionization target followed by the addition of the matrix. Acquisition of the data was performed using the positive ion mode.

RESULTS

Based on the identification of clear and defined lipid signatures from the seven reference strains, the method that we developed was fast (<10 minutes) and produced interpretable profiles for all but four isolates, caused by poor ionization giving an n = 86 with interpretable spectra. The sensitivity and specificity of the matrix-assisted laser desorption/ionization time of flight were 94.4 (95% CI, 86.4-98.5) and 94.4 (95% CI, 72.7-99.9), respectively.

CONCLUSIONS

Mycobacterial lipid profiling provides a means of rapid, safe, and accurate discrimination of species within the MTBC.

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.

Inhibition of the Niemann-Pick C1 protein is a conserved feature of multiple strains of pathogenic mycobacteria

Mycobacterium tuberculosis (Mtb) survives and replicates within host macrophages (MΦ) and subverts multiple antimicrobial defense mechanisms. Previously, we reported that lipids shed by pathogenic mycobacteria inhibit NPC1, the lysosomal membrane protein deficient in the lysosomal storage disorder Niemann-Pick disease type C (NPC). Inhibition of NPC1 leads to a drop in lysosomal calcium levels, blocking phagosome-lysosome fusion leading to mycobacterial survival. We speculated that the production of specific cell wall lipid(s) that inhibit NPC1 could have been a critical step in the evolution of pathogenicity. We therefore investigated whether lipid extracts from clinical Mtb strains from multiple Mtb lineages, Mtb complex (MTBC) members and non-tubercular mycobacteria (NTM) inhibit the NPC pathway. We report that inhibition of the NPC pathway was present in all clinical isolates from Mtb lineages 1, 2, 3 and 4, Mycobacterium bovis and the NTM, Mycobacterium abscessus and Mycobacterium avium. However, lipid extract from Mycobacterium canettii, which is considered to resemble the common ancestor of the MTBC did not inhibit the NPC1 pathway. We conclude that the evolution of NPC1 inhibitory mycobacterial cell wall lipids evolved early and post divergence from Mycobacterium canettii-related mycobacteria and that this activity contributes significantly to the promotion of disease.

Identification of region of difference and H37Rv-related deletion in Mycobacterium tuberculosis complex by structural variant detection and genome assembly

Mycobacterium tuberculosis complex (MTBC), the main cause of TB in humans and animals, is an extreme example of genetic homogeneity, whereas it is still nevertheless separated into various lineages by numerous typing methods, which differ in phenotype, virulence, geographic distribution, and host preference. The large sequence polymorphism (LSP), incorporating region of difference (RD) and H37Rv-related deletion (RvD), is considered to be a powerful means of constructing phylogenetic relationships within MTBC. Although there have been many studies on LSP already, focusing on the distribution of RDs in MTBC and their impact on MTB phenotypes, a crumb of new lineages or sub-lineages have been excluded and RvDs have received less attention. We, therefore, sampled a dataset of 1,495 strains, containing 113 lineages from the laboratory collection, to screen for RDs and RvDs by structural variant detection and genome assembly, and examined the distribution of RvDs in MTBC, including RvD2, RvD5, and cobF region. Consistent with genealogical delineation by single nucleotide polymorphism (SNP), we identified 125 RDs and 5 RvDs at the species, lineage, or sub-lineage levels. The specificities of RDs and RvDs were further investigated in the remaining 10,218 strains, suggesting that most of them were highly specific to distinct phylogenetic groups, could be used as stable genetic markers in genotyping. More importantly, we identified 34 new lineage or evolutionary branch specific RDs and 2 RvDs, also demonstrated the distribution of known RDs and RvDs in MTBC. This study provides novel details about deletion events that have occurred in distinct phylogenetic groups and may help to understand the genealogical differentiation.

Understanding drivers of phylogenetic clustering and terminal branch lengths distribution in epidemics of Mycobacterium tuberculosis

Detecting factors associated with transmission is important to understand disease epidemics, and to design effective public health measures. Clustering and terminal branch lengths (TBL) analyses are commonly applied to genomic data sets of Mycobacterium tuberculosis (MTB) to identify sub-populations with increased transmission. Here, I used a simulation-based approach to investigate what epidemiological processes influence the results of clustering and TBL analyses, and whether differences in transmission can be detected with these methods. I simulated MTB epidemics with different dynamics (latency, infectious period, transmission rate, basic reproductive number R0, sampling proportion, sampling period, and molecular clock), and found that all considered factors, except for the length of the infectious period, affect the results of clustering and TBL distributions. I show that standard interpretations of this type of analyses ignore two main caveats: (1) clustering results and TBL depend on many factors that have nothing to do with transmission, (2) clustering results and TBL do not tell anything about whether the epidemic is stable, growing, or shrinking, unless all the additional parameters that influence these metrics are known, or assumed identical between sub-populations. An important consequence is that the optimal SNP threshold for clustering depends on the epidemiological conditions, and that sub-populations with different epidemiological characteristics should not be analyzed with the same threshold. Finally, these results suggest that different clustering rates and TBL distributions, that are found consistently between different MTB lineages, are probably due to intrinsic bacterial factors, and do not indicate necessarily differences in transmission or evolutionary success.

Origin and Global Expansion of Mycobacterium tuberculosis Complex Lineage 3

Simple Summary

Tuberculosis still causes 1.5 million deaths annually and is mainly caused by Mycobacterium tuberculosis complex strains belonging to three evolutionary modern lineages (Lineages 2–4). While Lineage 2 and Lineage 4 virtually conquered the world, Lineage 3 is particularly successful in Northern and Eastern Africa, as well as in Southern Asia, the suspected evolutionary origin of these strains. Here, we sought to understand how Lineage 3 strains came to the African continent. To this end, we performed routine genotyping to characterize over 2500 clinical isolates from 38 countries. We then selected a representative collection of 373 isolates for a whole-genome analysis and a modeling approach to infer the geographic origin of different sublineages. In fact, the origin of Lineage 3 could be located in India, and we found evidence for independent introductions of four distinct sublineages into North/East Africa, in line with known ancient exchanges and migrations between both world regions. Our study illustrates that the evolutionary history of humans and their pathogens are closely connected and further provides a systematic understanding of the genomic diversity of Lineage 3, which could be important for the development of new tuberculosis vaccines or new therapeutics.

Abstract

Mycobacterium tuberculosis complex (MTBC) Lineage 3 (L3) strains are abundant in world regions with the highest tuberculosis burden. To investigate the population structure and the global diversity of this major lineage, we analyzed a dataset comprising 2682 L3 strains from 38 countries over 5 continents, by employing 24-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats genotyping (MIRU-VNTR) and drug susceptibility testing. We further combined whole-genome sequencing (WGS) and phylogeographic analysis for 373 strains representing the global L3 genetic diversity. Ancestral state reconstruction confirmed that the origin of L3 strains is located in Southern Asia and further revealed multiple independent introduction events into North-East and East Africa. This study provides a systematic understanding of the global diversity of L3 strains and reports phylogenetic variations that could inform clinical trials which evaluate the effectivity of new drugs/regimens or vaccine candidates.

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.

Rapid Identification of Mycobacterium tuberculosis Complex Using Mass Spectrometry: A Proof of Concept

Mycobacteria that form the Mycobacterium tuberculosis complex are responsible for deadly tuberculosis in animals and patients. Identification of these pathogens at the species level is of primary importance for treatment and source tracing and currently relies on DNA analysis, including whole genome sequencing (WGS), which requires a whole day. In this study, we report the unprecedented discrimination of M. tuberculosis complex species using matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF-MS), with WGS as the comparative reference standard. In the first step, optimized peptide extraction applied to 36 isolates otherwise identified in five of the 11 M. tuberculosis complex variants by WGS yielded 139 MALDI-TOF spectra, which were used to identify biomarkers of interest that facilitate differentiation between variants. In a second step, 70/80 (88%) other isolates were correctly classified by an algorithm based on specific peaks. This study is the first to report a MALDI-TOF-MS method for discriminating M. tuberculosis complex mycobacteria that is easily implemented in clinical microbiology laboratories.

Sub-Lineage Specific Phenolic Glycolipid Patterns in the Mycobacterium tuberculosis Complex Lineage 1

“Ancestral” Mycobacterium tuberculosis complex (MTBC) strains of Lineage 1 (L1, East African Indian) are a prominent tuberculosis (TB) cause in countries around the Indian Ocean. However, the pathobiology of L1 strains is insufficiently characterized. Here, we used whole genome sequencing (WGS) of 312 L1 strains from 43 countries to perform a characterization of the global L1 population structure and correlate this to the analysis of the synthesis of phenolic glycolipids (PGL) – known MTBC polyketide-derived virulence factors. Our results reveal the presence of eight major L1 sub-lineages, whose members have specific mutation signatures in PGL biosynthesis genes, e.g., pks15/1 or glycosyltransferases Rv2962c and/or Rv2958c. Sub-lineage specific PGL production was studied by NMR-based lipid profiling and strains with a completely abolished phenolphthiocerol dimycoserosate biosynthesis showed in average a more prominent growth in human macrophages. In conclusion, our results show a diverse population structure of L1 strains that is associated with the presence of specific PGL types. This includes the occurrence of mycoside B in one sub-lineage, representing the first description of a PGL in an M. tuberculosis lineage other than L2. Such differences may be important for the evolution of L1 strains, e.g., allowing adaption to different human populations.

Dual RNA Sequencing of Mycobacterium tuberculosis-Infected Human Splenic Macrophages Reveals a Strain-Dependent Host-Pathogen Response to Infection

Tuberculosis (TB) is caused by Mycobacterium tuberculosis (Mtb), leading to pulmonary and extrapulmonary TB, whereby Mtb is disseminated to many other organs and tissues. Dissemination occurs early during the disease, and bacteria can be found first in the lymph nodes adjacent to the lungs and then later in the extrapulmonary organs, including the spleen. The early global gene expression response of human tissue macrophages and intracellular clinical isolates of Mtb has been poorly studied. Using dual RNA-seq, we have explored the mRNA profiles of two closely related clinical strains of the Latin American and Mediterranean (LAM) family of Mtb in infected human splenic macrophages (hSMs). This work shows that these pathogens mediate a distinct host response despite their genetic similarity. Using a genome-scale host–pathogen metabolic reconstruction to analyze the data further, we highlight that the infecting Mtb strain also determines the metabolic response of both the host and pathogen. Thus, macrophage ontogeny and the genetic-derived program of Mtb direct the host–pathogen interaction.

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.

Genotypic and phenotypic diversity of Mycobacterium tuberculosis complex genotypes prevalent in West Africa

Findings from previous comparative genomics studies of the Mycobacterium tuberculosis complex (MTBC) suggest genomic variation among the genotypes may have phenotypic implications. We investigated the diversity in the phenotypic profiles of the main prevalent MTBC genotypes in West Africa. Thirty-six whole genome sequenced drug susceptible MTBC isolates belonging to lineages 4, 5 and 6 were included in this study. The isolates were phenotypically characterized for urease activity, tween hydrolysis, Thiophen-2-Carboxylic Acid Hydrazide (TCH) susceptibility, nitric oxide production, and growth rate in both liquid (7H9) and solid media (7H11 and Löwenstein-Jensen (L-J)). Lineage 4 isolates showed the highest growth rate in both liquid (p = 0.0003) and on solid (L-J) media supplemented with glycerol (p<0.001) or pyruvate (p = 0.005). L6 isolates optimally utilized pyruvate compared to glycerol (p<0.001), whereas L5 isolates grew similarly on both media (p = 0.05). Lineage 4 isolates showed the lowest average time to positivity (TTP) (p = 0.01; Average TTP: L4 = 15days, L5 = 16.7days, L6 = 29.7days) and the highest logCFU/mL (p = 0.04; average logCFU/mL L4 = 5.9, L5 = 5.0, L6 = 4.4) on 7H11 supplemented with glycerol, but there was no significant difference in growth on 7H11 supplemented with pyruvate (p = 0.23). The highest release of nitrite was recorded for L5 isolates, followed by L4 and L6 isolates. However, the reverse was observed in the urease activity for the lineages. All isolates tested were resistant to TCH except for one L6 isolate. Comparative genomic analyses revealed several mutations that might explain the diverse phenotypic profiles of these isolates. Our findings showed significant phenotypic diversity among the MTBC lineages used for this study.

Mycobacterium tuberculosis complex lineage 5 exhibits high levels of within-lineage genomic diversity and differing gene content compared to the type strain H37Rv

Pathogens of the Mycobacterium tuberculosis complex (MTBC) are considered to be monomorphic, with little gene content variation between strains. Nevertheless, several genotypic and phenotypic factors separate strains of the different MTBC lineages (L), especially L5 and L6 (traditionally termed Mycobacterium africanum) strains, from each other. However, this genome variability and gene content, especially of L5 strains, has not been fully explored and may be important for pathobiology and current approaches for genomic analysis of MTBC strains, including transmission studies. By comparing the genomes of 355 L5 clinical strains (including 3 complete genomes and 352 Illumina whole-genome sequenced isolates) to each other and to H37Rv, we identified multiple genes that were differentially present or absent between H37Rv and L5 strains. Additionally, considerable gene content variability was found across L5 strains, including a split in the L5.3 sub-lineage into L5.3.1 and L5.3.2. These gene content differences had a small knock-on effect on transmission cluster estimation, with clustering rates influenced by the selected reference genome, and with potential overestimation of recent transmission when using H37Rv as the reference genome. We conclude that full capture of the gene diversity, especially high-resolution outbreak analysis, requires a variation of the single H37Rv-centric reference genome mapping approach currently used in most whole-genome sequencing data analysis pipelines. Moreover, the high within-lineage gene content variability suggests that the pan-genome of M. tuberculosis is at least several kilobases larger than previously thought, implying that a concatenated or reference-free genome assembly (de novo) approach may be needed for particular questions.

Heterogeneous Host-Pathogen Encounters Coordinate Antibiotic Resilience in Mycobacterium tuberculosis

Successful treatment of tuberculosis (TB) depends on the eradication of its causative agent Mycobacterium tuberculosis (Mtb) in the host. However, the emergence of phenotypically drug-resistant Mtb in the host environment tempers the ability of antibiotics to cure disease. Host immunity produces diverse microenvironmental niches that are exploited by Mtb to mobilize adaptation programs. Such differential interactions amplify pre-existing heterogeneity in the host-pathogen milieu to influence disease pathology and therapy outcome. Therefore, comprehending the intricacies of phenotypic heterogeneity can be an empirical step forward in potentiating drug action. With this goal, we review the interconnectedness of the lesional, cellular, and bacterial heterogeneity underlying phenotypic drug resistance. Based on this information, we anticipate the development of new therapeutic strategies targeting host-pathogen heterogeneity to cure TB.

Underwhelming or Misunderstood? Genetic Variability of Pattern Recognition Receptors in Immune Responses and Resistance to Mycobacterium tuberculosis

Human genetic control is thought to affect a considerable part of the outcome of infection with Mycobacterium tuberculosis (Mtb). Most of us deal with the pathogen by containment (associated with clinical "latency") or sterilization, but tragically millions each year do not. After decades of studies on host genetic susceptibility to Mtb infection, genetic variation has been discovered to play a role in tuberculous immunoreactivity and tuberculosis (TB) disease. Genes encoding pattern recognition receptors (PRRs) enable a consistent, molecularly direct interaction between humans and Mtb which suggests the potential for co-evolution. In this review, we explore the roles ascribed to PRRs during Mtb infection and ask whether such a longstanding and intimate interface between our immune system and this pathogen plays a critical role in determining the outcome of Mtb infection. The scientific evidence to date suggests that PRR variation is clearly implicated in altered immunity to Mtb but has a more subtle role in limiting the pathogen and pathogenesis. In contrast to 'effectors' like IFN-γ, IL-12, Nitric Oxide and TNF that are critical for Mtb control, 'sensors' like PRRs are less critical for the outcome of Mtb infection. This is potentially due to redundancy of the numerous PRRs in the innate arsenal, such that Mtb rarely goes unnoticed. Genetic association studies investigating PRRs during Mtb infection should therefore be designed to investigate endophenotypes of infection - such as immunological or clinical variation - rather than just TB disease, if we hope to understand the molecular interface between innate immunity and Mtb.

Highly transmitted M. tuberculosis strains are more likely to evolve MDR/XDR and cause outbreaks, but what makes them highly transmitted?

Multi-Drug-Resistant strains of Mycobacterium tuberculosis (MDR-TB) are a serious obstacle to global TB eradication. While most MDR-TB strains are infrequently transmitted, a few cause large transmission clusters that contribute substantially to local MDR-TB burdens. Here we examine whether the known mutations in these strains can explain their success. Drug resistance mutations differ in fitness costs and strains can also acquire compensatory mutations (CM) to restore fitness, but some highly transmitted MDR strains have no CM. The acquisition of resistance mutations that maintain high transmissibility seems to occur by chance and are more likely in strains that are intrinsically highly transmitted and cause many cases. Modern Beijing lineage strains have caused several large outbreaks, but MDR outbreaks are also caused by ancient Beijing and lineage 4 strains, suggesting the lineage is less important than the characteristics of the individual strain. The development of fluoroquinolone resistance appears to represent another level of selection, in which strains must surmount unknown fitness costs of gyrA mutations. The genetic determinants of high transmission are poorly defined but may involve genes encoding proteins involved in molybdenum acquisition and the Esx systems. In addition, strains eliciting lower cytokine responses and producing more caseating granulomas may have advantages for transmission. Successful MDR/XDR strains generally evolve from highly transmitted drug sensitive parent strains due to selection pressures from deficiencies in local TB control programs. Until TB incidence is considerably reduced, there will likely be highly transmitted strains that develop resistance to any new antibiotic.

ESX-1-Independent Horizontal Gene Transfer by Mycobacterium tuberculosis Complex Strains

Current models of horizontal gene transfer (HGT) in mycobacteria are based on “distributive conjugal transfer” (DCT), an HGT type described in the fast-growing, saprophytic model organism Mycobacterium smegmatis, which creates genome mosaicism in resulting strains and depends on an ESX-1 type VII secretion system. In contrast, only few data on interstrain DNA transfer are available for tuberculosis-causing mycobacteria, for which chromosomal DNA transfer between two Mycobacterium canettii strains was reported, a process which, however, was not observed for Mycobacterium tuberculosis strains. Here, we have studied a wide range of human- and animal-adapted members of the Mycobacterium tuberculosis complex (MTBC) using an optimized filter-based mating assay together with three selected strains of M. canettii that acted as DNA recipients. Unlike in previous approaches, we obtained a high yield of thousands of recombinants containing transferred chromosomal DNA fragments from various MTBC donor strains, as confirmed by whole-genome sequence analysis of 38 randomly selected clones. While the genome organizations of the obtained recombinants showed mosaicisms of donor DNA fragments randomly integrated into a recipient genome backbone, reminiscent of those described as being the result of ESX-1-mediated DCT in M. smegmatis, we observed similar transfer efficiencies when ESX-1-deficient donor and/or recipient mutants were used, arguing that in tubercle bacilli, HGT is an ESX-1-independent process. These findings provide new insights into the genetic events driving the pathoevolution of M. tuberculosis and radically change our perception of HGT in mycobacteria, particularly for those species that show recombinogenic population structures despite the natural absence of ESX-1 secretion systems.

Breaching the phagosome, the case of the tuberculosis agent

The interactions between microbes and their hosts are among the most complex biological phenomena known today. The interaction may reach from overall beneficial interaction, as observed for most microbiome/microbiota related interactions to interaction with virulent pathogens, against which host cells have evolved sophisticated defence strategies. Among the latter, the confinement of invading pathogens in a phagosome plays a key role, which often results in the destruction of the invader, whereas some pathogens may counteract phagosomal arrest and survive by gaining access to the cytosol of the host cell. In the current review, we will discuss recent insights into this dynamic process of host-pathogen interaction, using Mycobacterium tuberculosis and related pathogenic mycobacteria as main examples.

Pan and Core Genome Analysis of 183 Mycobacterium tuberculosis Strains Revealed a High Inter-Species Diversity among the Human Adapted Strains

Tuberculosis (TB) is an airborne communicable disease with high morbidity and mortality rates, especially in developing countries. The causal agents of TB belong to the complex Mycobacterium tuberculosis (MTBc), which is composed of different human and animal TB associated species. Some animal associated species have zoonotic potential and add to the burden of TB management. The BCG ("Bacillus Calmette-Guérin") vaccine is widely used for the prevention against TB, but its use is limited in immunocompromised patients and animals due to the adverse effects and disseminated life-threatening complications. In this study, we aimed to carry out a comparative genome analysis between the human adapted species including BCG vaccine strains to identify and pinpoint the conserved genes related to the virulence across all the species, which could add a new value for vaccine development. For this purpose, the sequences of 183 Mycobacterium tuberculosis (MTB) strains were retrieved from the freely available WGS dataset at NCBI. The species included: 168 sensu stricto MTB species with other human MTB complex associated strains: M. tuberculosis var. africanum (3), M. tuberculosis var. bovis (2 draft genomes) and 10 BCG species, which enabled the analysis of core genome which contains the conserved genes and some virulence factor determinants. Further, a phylogenetic tree was constructed including the genomes of human (183); animals MTB adapted strains (6) and the environmental Mycobacterium strain "M. canettii". Our results showed that the core genome consists of 1166 conserved genes among these species, which represents a small portion of the pangenome (7036 genes). The remaining genes in the pangenome (5870) are accessory genes, adding a high inter-species diversity. Further, the core genome includes several virulence-associated genes and this could explain the rare infectiousness potential of some attenuated vaccine strains in some patients. This study reveals that low number of conserved genes in human adapted MTBc species and high inter-species diversity of the pan-genome could be considered for vaccine candidate development.

Lipoarabinomannan as a Point-of-Care Assay for Diagnosis of Tuberculosis: How Far Are We to Use It?

Tuberculosis (TB) is still a severe public health problem; the current diagnostic tests have limitations that delay treatment onset. Lipoarabinomannan (LAM) is a glycolipid that is a component of the cell wall of the bacillus Mycobacterium tuberculosis, the etiologic agent of TB. This glycolipid is excreted as a soluble form in urine. The World Health Organization has established that the design of new TB diagnostic methods is one of the priorities within the EndTB Strategy. LAM has been suggested as a biomarker to develop diagnostic tests based on its identification in urine, and it is one of the most prominent candidates to develop point-of-care diagnostic test because urine samples can be easily collected. Moreover, LAM can regulate the immune response in the host and can be found in the serum of TB patients, where it probably affects a wide variety of host cell populations, consequently influencing the quality of both innate and adaptive immune responses during TB infection. Here, we revised the evidence that supports that LAM could be used as a tool for the development of new point-of-care tests for TB diagnosis, and we discussed the mechanisms that could contribute to the low sensitivity of diagnostic testing.

One Size Fits All? Not in In Vivo Modeling of Tuberculosis Chemotherapeutics

Tuberculosis (TB) remains a global health problem despite almost universal efforts to provide patients with highly effective chemotherapy, in part, because many infected individuals are not diagnosed and treated, others do not complete treatment, and a small proportion harbor Mycobacterium tuberculosis (Mtb) strains that have become resistant to drugs in the standard regimen. Development and approval of new drugs for TB have accelerated in the last 10 years, but more drugs are needed due to both Mtb's development of resistance and the desire to shorten therapy to 4 months or less. The drug development process needs predictive animal models that recapitulate the complex pathology and bacterial burden distribution of human disease. The human host response to pulmonary infection with Mtb is granulomatous inflammation usually resulting in contained lesions and limited bacterial replication. In those who develop progressive or active disease, regions of necrosis and cavitation can develop leading to lasting lung damage and possible death. This review describes the major vertebrate animal models used in evaluating compound activity against Mtb and the disease presentation that develops. Each of the models, including the zebrafish, various mice, guinea pigs, rabbits, and non-human primates provides data on number of Mtb bacteria and pathology resolution. The models where individual lesions can be dissected from the tissue or sampled can also provide data on lesion-specific bacterial loads and lesion-specific drug concentrations. With the inclusion of medical imaging, a compound's effect on resolution of pathology within individual lesions and animals can also be determined over time. Incorporation of measurement of drug exposure and drug distribution within animals and their tissues is important for choosing the best compounds to push toward the clinic and to the development of better regimens. We review the practical aspects of each model and the advantages and limitations of each in order to promote choosing a rational combination of them for a compound's development.

Different PPD-stimulated cytokine responses from patients infected with genetically distinct Mycobacterium tuberculosis complex lineages

OBJECTIVES

The genetic diversity of Mycobacterium tuberculosis complex (MTBC) influences the immune response of the host, which may affect the immunodiagnostic tests and biomarker assessment studies used for tuberculosis (TB). This study aimed to determine whether the mycobacterial-antigen-stimulated cytokine responses vary with the genotype of the MTBC infecting the patient.

METHODS

Eighty-one patients with confirmed active pulmonary TB were recruited, and MTBC clinical strains were isolated from their sputum for bacterial lineage single-nucleotide polymorphism typing. Whole blood was drawn from the patients to measure the purified protein derivative (PPD)-stimulated cytokine responses (GM-CSF, IFN-γ, IL-1β, IL-2, IL-4, IL-5, IL-6, IL-7, IL-8, IL-10, TNF-α, IFN-α, IL-12, eotaxin, IL-13, IL-15, IL-17, MIP1-α, MIP1-β, MCP1, IL1RA, IP10, IL2R, MIG) with the Luminex multiplex immunoassay.

RESULTS

Of the 24 cytokines studied, three were produced differentially in whole blood dependent on the infecting lineage of MTBC. Decreased production of IL-17 was observed in patients infected with modern lineages compared with patients infected with ancestral lineages (P < 0.01), and production of IFN-γ and IL-2 was significantly decreased in patients infected with lineage 4 strains compared with patients infected with lineage 3 strains (P < 0.05).

CONCLUSION

MTBC strains belonging to lineage 4 induced a decreased whole-blood PPD-stimulated pro-inflammatory cytokine response.

Local adaptation in populations of Mycobacterium tuberculosis endemic to the Indian Ocean Rim

Background: Lineage 1 (L1) and 3 (L3) are two lineages of the Mycobacterium tuberculosis complex (MTBC) causing tuberculosis (TB) in humans. L1 and L3 are prevalent around the rim of the Indian Ocean, the region that accounts for most of the world's new TB cases. Despite their relevance for this region, L1 and L3 remain understudied. Methods: We analyzed 2,938 L1 and 2,030 L3 whole genome sequences originating from 69 countries. We reconstructed the evolutionary history of these two lineages and identified genes under positive selection. Results: We found a strongly asymmetric pattern of migration from South Asia toward neighboring regions, highlighting the historical role of South Asia in the dispersion of L1 and L3. Moreover, we found that several genes were under positive selection, including genes involved in virulence and resistance to antibiotics. For L1 we identified signatures of local adaptation at the esxH locus, a gene coding for a secreted effector that targets the human endosomal sorting complex, and is included in several vaccine candidates. Conclusions: Our study highlights the importance of genetic diversity in the MTBC, and sheds new light on two of the most important MTBC lineages affecting humans.

Survival of hypoxia-induced dormancy is not a common feature of all strains of the Mycobacterium tuberculosis complex

While persistence in a dormant state is crucial for the life cycle of Mycobacterium tuberculosis, no investigation regarding dormancy survival of different strains across different lineages was performed so far. We analyzed responses to oxygen starvation and recovery in terms of growth, metabolism, and transcription. All different strains belonging to the Euro-American lineage (L4) showed similar survival and resuscitation characteristics. Different clinical isolates from the Beijing (L2), East African-Indian (L3), and Delhi/Central Asian (L1) lineage did not survive oxygen starvation. We show that dormancy survival is lineage-dependent. Recovery from O₂ starvation was only observed in strains belonging to the Euro-American (L4) lineage but not in strains belonging to different lineages (L1, L2, L3). Thus, resuscitation from dormancy after oxygen starvation is not a general feature of all M. tuberculosis strains as thought before. Our findings are of key importance to understand infection dynamics of non-Euro-American vs Euro-American strains and to develop drugs targeting the dormant state.

The genetic proteome: Using genetics to inform the proteome of mycobacterial pathogens

Mycobacterial pathogens pose a sustained threat to human health. There is a critical need for new diagnostics, therapeutics, and vaccines targeting both tuberculous and nontuberculous mycobacterial species. Understanding the basic mechanisms used by diverse mycobacterial species to cause disease will facilitate efforts to design new approaches toward detection, treatment, and prevention of mycobacterial disease. Molecular, genetic, and biochemical approaches have been widely employed to define fundamental aspects of mycobacterial physiology and virulence. The recent expansion of genetic tools in mycobacteria has further increased the accessibility of forward genetic approaches. Proteomics has also emerged as a powerful approach to further our understanding of diverse mycobacterial species. Detection of large numbers of proteins and their modifications from complex mixtures of mycobacterial proteins is now routine, with efforts of quantification of these datasets becoming more robust. In this review, we discuss the “genetic proteome,” how the power of genetics, molecular biology, and biochemistry informs and amplifies the quality of subsequent analytical approaches and maximizes the potential of hypothesis-driven mycobacterial research. Published proteomics datasets can be used for hypothesis generation and effective post hoc supplementation to experimental data. Overall, we highlight how the integration of proteomics, genetic, molecular, and biochemical approaches can be employed successfully to define fundamental aspects of mycobacterial pathobiology.

Tuberculosis diagnostics: overcoming ancient challenges with modern solutions

Rapid, sensitive, accurate and portable diagnostics are a mainstay of modern medicine. Tuberculosis is a disease that has been with us since time immemorial and, despite the fact that it can be treated and cured, it still remains the world's biggest infectious killer, taking the lives of millions annually. There have been important developments in the diagnostic devices for tuberculosis however, these are often prone to error, expensive, lack the necessary sensitivity or accuracy and, crucially, not sufficiently portable and thus not applicable in the remote, rural areas, where they are most needed. Modern solutions have been emerging in the past decade, seeking to overcome many of the inhibiting issues in this field by utilising recent advances in molecular biology, genetics and sequencing or even completely ‘reinventing the wheel’, by developing novel and unprecedented diagnostic techniques. In this mini review, the issues and challenges arising from the historical methods of diagnosing tuberculosis are discussed, followed by outlaying their particular lack of appropriateness for regions of the world where tuberculosis still remains endemic. Subsequently, more recent developments of new methods and technological advancements as ‘modern weapons’ in the battle to defeat this disease and associated challenges are reviewed, and finally an outlook is presented, highlighting the future of the modern solutions under development, which are envisioned to lay the platform for improvements in delivering timely intervention, reduce immense expense and burden on healthcare systems worldwide, while saving millions of lives and eventually, may enable the eradication of this ancient disease.

6,11-Dioxobenzo[f]pyrido[1,2-a]indoles Kill Mycobacterium tuberculosis by Targeting Iron-Sulfur Protein Rv0338c (IspQ), A Putative Redox Sensor

Screening of a diversity-oriented compound library led to the identification of two 6,11-dioxobenzo[f]pyrido[1,2-a]indoles (DBPI) that displayed low micromolar bactericidal activity against the Erdman strain of Mycobacterium tuberculosis in vitro. The activity of these hit compounds was limited to tubercle bacilli, including the nonreplicating form, and to Mycobacterium marinum. On hit expansion and investigation of the structure activity relationship, selected modifications to the dioxo moiety of the DBPI scaffold were either neutral or led to reduction or abolition of antimycobacterial activity. To find the target, DBPI-resistant mutants of M. tuberculosis Erdman were raised and characterized first microbiologically and then by whole genome sequencing. Four different mutations, all affecting highly conserved residues, were uncovered in the essential gene rv0338c (ispQ) that encodes a membrane-bound protein, named IspQ, with 2Fe-2S and 4Fe-4S centers and putative iron-sulfur-binding reductase activity. With the help of a structural model, two of the mutations were localized close to the 2Fe-2S domain in IspQ and another in transmembrane segment 3. The mutant genes were recessive to the wild type in complementation experiments and further confirmation of the hit-target relationship was obtained using a conditional knockdown mutant of rv0338c in M. tuberculosis H37Rv. More mechanistic insight was obtained from transcriptome analysis, following exposure of M. tuberculosis to two different DBPI; this revealed strong upregulation of the redox-sensitive SigK regulon and genes induced by oxidative and thiol-stress. The findings of this investigation pharmacologically validate a novel target in tubercle bacilli and open a new vista for tuberculosis drug discovery.

Lineage-Specific Proteomic Signatures in the Mycobacterium tuberculosis Complex Reveal Differential Abundance of Proteins Involved in Virulence, DNA Repair, CRISPR-Cas, Bioenergetics and Lipid Metabolism

Despite the discovery of the tubercle bacillus more than 130 years ago, its physiology and the mechanisms of virulence are still not fully understood. A comprehensive analysis of the proteomes of members of the human-adapted Mycobacterium tuberculosis complex (MTBC) lineages 3, 4, 5, and 7 was conducted to better understand the evolution of virulence and other physiological characteristics. Unique and shared proteomic signatures in these modern, pre-modern and ancient MTBC lineages, as deduced from quantitative bioinformatics analyses of high-resolution mass spectrometry data, were delineated. The main proteomic findings were verified by using immunoblotting. In addition, analysis of multiple genome alignment of members of the same lineages was performed. Label-free peptide quantification of whole cells from MTBC lineages 3, 4, 5, and 7 yielded a total of 38,346 unique peptides derived from 3092 proteins, representing 77% coverage of the predicted proteome. MTBC lineage-specific differential expression was observed for 539 proteins. Lineage 7 exhibited a markedly reduced abundance of proteins involved in DNA repair, type VII ESX-3 and ESX-1 secretion systems, lipid metabolism and inorganic phosphate uptake, and an increased abundance of proteins involved in alternative pathways of the TCA cycle and the CRISPR-Cas system as compared to the other lineages. Lineages 3 and 4 exhibited a higher abundance of proteins involved in virulence, DNA repair, drug resistance and other metabolic pathways. The high throughput analysis of the MTBC proteome by super-resolution mass spectrometry provided an insight into the differential expression of proteins between MTBC lineages 3, 4, 5, and 7 that may explain the slow growth and reduced virulence, metabolic flexibility, and the ability to survive under adverse growth conditions of lineage 7.

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.