Ecology and evolution of Mycobacterium tuberculosis

Newly Identified Mycobacterium africanum Lineage 10, Central Africa

Analysis of genome sequencing data from >100,000 genomes of Mycobacterium tuberculosis complex using TB-Annotator software revealed a previously unknown lineage, proposed name L10, in central Africa. Phylogenetic reconstruction suggests L10 could represent a missing link in the evolutionary and geographic migration histories of M. africanum.

An additional proofreader contributes to DNA replication fidelity in mycobacteria

The removal of mis-incorporated nucleotides by proofreading activity ensures DNA replication fidelity. Whereas the ε-exonuclease DnaQ is a well-established proofreader in the model organism Escherichia coli, it has been shown that proofreading in a majority of bacteria relies on the polymerase and histidinol phosphatase (PHP) domain of replicative polymerase, despite the presence of a DnaQ homolog that is structurally and functionally distinct from E. coli DnaQ. However, the biological functions of this type of noncanonical DnaQ remain unclear. Here, we provide independent evidence that noncanonical DnaQ functions as an additional proofreader for mycobacteria. Using the mutation accumulation assay in combination with whole-genome sequencing, we showed that depletion of DnaQ in Mycolicibacterium smegmatis leads to an increased mutation rate, resulting in AT-biased mutagenesis and increased insertions/deletions in the homopolymer tract. Our results showed that mycobacterial DnaQ binds to the β clamp and functions synergistically with the PHP domain proofreader to correct replication errors. Furthermore, the loss of dnaQ results in replication fork dysfunction, leading to attenuated growth and increased mutagenesis on subinhibitory fluoroquinolones potentially due to increased vulnerability to fork collapse. By analyzing the sequence polymorphism of dnaQ in clinical isolates of Mycobacterium tuberculosis (Mtb), we demonstrated that a naturally evolved DnaQ variant prevalent in Mtb lineage 4.3 may enable hypermutability and is associated with drug resistance. These results establish a coproofreading model and suggest a division of labor between DnaQ and PHP domain proofreader. This study also provides real-world evidence that a mutator-driven evolutionary pathway may exist during the adaptation of Mtb.

Automated Diagnosis and Phenotyping of Tuberculosis Using Serum Metabolic Fingerprints

Tuberculosis (TB) stands as the second most fatal infectious disease after COVID-19, the effective treatment of which depends on accurate diagnosis and phenotyping. Metabolomics provides valuable insights into the identification of differential metabolites for disease diagnosis and phenotyping. However, TB diagnosis and phenotyping remain great challenges due to the lack of a satisfactory metabolic approach. Here, a metabolomics-based diagnostic method for rapid TB detection is reported. Serum metabolic fingerprints are examined via an automated nanoparticle-enhanced laser desorption/ionization mass spectrometry platform outstanding by its rapid detection speed (measured in seconds), minimal sample consumption (in nanoliters), and cost-effectiveness (approximately $3). A panel of 14 m z-1 features is identified as biomarkers for TB diagnosis and a panel of 4 m z-1 features for TB phenotyping. Based on the acquired biomarkers, TB metabolic models are constructed through advanced machine learning algorithms. The robust metabolic model yields a 97.8% (95% confidence interval (CI), 0.964-0.986) area under the curve (AUC) in TB diagnosis and an 85.7% (95% CI, 0.806-0.891) AUC in phenotyping. In this study, serum metabolic biomarker panels are revealed and develop an accurate metabolic tool with desirable diagnostic performance for TB diagnosis and phenotyping, which may expedite the effective implementation of the end-TB strategy.

Whole-genome sequencing-based genetic diversity, transmission dynamics, and drug-resistant mutations in Mycobacterium tuberculosis isolated from extrapulmonary tuberculosis patients in western Ethiopia

Background

Extrapulmonary tuberculosis (EPTB) refers to a form of Tuberculosis (TB) where the infection occurs outside the lungs. Despite EPTB being a devastating disease of public health concern, it is frequently overlooked as a public health problem. This study aimed to investigate genetic diversity, identify drug-resistance mutations, and trace ongoing transmission chains.

Methods

A cross-sectional study was undertaken on individuals with EPTB in western Ethiopia. In this study, whole-genome sequencing (WGS) was employed to analyze Mycobacterium tuberculosis (MTB) samples obtained from EPTB patients. Out of the 96 genomes initially sequenced, 89 met the required quality standards for genetic diversity, and drug-resistant mutations analysis. The data were processed using robust bioinformatics tools.

Results

Our analysis reveals that the majority (87.64%) of the isolates can be attributed to Lineage-4 (L4), with L4.6.3 and L4.2.2.2 emerging as the predominant sub-lineages, constituting 34.62% and 26.92%, respectively. The overall clustering rate and recent transmission index (RTI) were 30 and 17.24%, respectively. Notably, 7.87% of the isolates demonstrated resistance to at least one anti-TB drug, although multi-drug resistance (MDR) was observed in only 1.12% of the isolates.

Conclusions

The genetic diversity of MTBC strains in western Ethiopia was found to have low inter-lineage diversity, with L4 predominating and exhibiting high intra-lineage diversity. The notably high clustering rate in the region implies a pressing need for enhanced TB infection control measures to effectively disrupt the transmission chain. It's noteworthy that 68.75% of resistance-conferring mutations went undetected by both GeneXpert MTB/RIF and the line probe assay (LPA) in western Ethiopia. The identification of resistance mutations undetected by both GeneXpert and LPA, along with the detection of mixed infections through WGS, emphasizes the value of adopting WGS as a high-resolution approach for TB diagnosis and molecular epidemiological surveillance.

LILRB1-HLA-G axis defines a checkpoint driving natural killer cell exhaustion in tuberculosis

Chronic infections, including Mycobacterium tuberculosis (Mtb)-caused tuberculosis (TB), can induce host immune exhaustion. However, the key checkpoint molecules involved in this process and the underlying regulatory mechanisms remain largely undefined, which impede the application of checkpoint-based immunotherapy in infectious diseases. Here, through adopting time-of-flight mass cytometry and transcriptional profiling to systematically analyze natural killer (NK) cell surface receptors, we identify leukocyte immunoglobulin like receptor B1 (LILRB1) as a critical checkpoint receptor that defines a TB-associated cell subset (LILRB1+ NK cells) and drives NK cell exhaustion in TB. Mechanistically, Mtb-infected macrophages display high expression of human leukocyte antigen-G (HLA-G), which upregulates and activates LILRB1 on NK cells to impair their functions by inhibiting mitogen-activated protein kinase (MAPK) signaling via tyrosine phosphatases SHP1/2. Furthermore, LILRB1 blockade restores NK cell-dependent anti-Mtb immunity in immuno-humanized mice. Thus, LILRB1-HLA-G axis constitutes a NK cell immune checkpoint in TB and serves as a promising immunotherapy target.

Pan-lineage Mycobacterium tuberculosis reference genome for enhanced molecular diagnosis

In Mycobacterium tuberculosis (MTB) control, whole genome sequencing-based molecular drug susceptibility testing (molDST-WGS) has emerged as a pivotal tool. However, the current reliance on a single-strain reference limits molDST-WGS's true potential. To address this, we introduce a new pan-lineage reference genome, 'MtbRf'. We assembled 'unmapped' reads from 3,614 MTB genomes (751 L1; 881 L2; 1,700 L3; and 282 L4) into 35 shared, annotated contigs (54 coding sequences [CDSs]). We constructed MtbRf through: (1) searching for contig homologues among genome database that precipitate results uniquely within Mycobacteria genus; (2) comparing genomes with H37Rv ('lift-over') to define 18 insertions; and (3) filling gaps in H37Rv with insertions. MtbRf adds 1.18% sequences to H37rv, salvaging >60% of previously unmapped reads. Transcriptomics confirmed gene expression of new CDSs. The new variants provided a moderate DST predictive value (AUROC 0.60-0.75). MtbRf thus unveils previously hidden genomic information and lays the foundation for lineage-specific molDST-WGS.

Transcriptome analysis and molecular characterization of novel small RNAs in Mycobacterium tuberculosis Lineage 1

Mycobacterium tuberculosis (Mtb), the tuberculosis-causing agent, exhibits diverse genetic lineages, with known links to virulence. While genomic and transcriptomic variations between modern and ancient Mtb lineages have been explored, the role of small non-coding RNA (sRNA) in post-translational gene regulation remains largely uncharted. In this study, Mtb Lineage 1 (L1) Sabahan strains (n = 3) underwent sRNA sequencing, revealing 351 sRNAs, including 23 known sRNAs and 328 novel ones identified using ANNOgesic. Thirteen sRNAs were selected based on the best average cut-off value of 300, with RT-qPCR revealing significant expression differences for sRNA 1 (p = 0.0132) and sRNA 29 (p = 0.0012) between Mtb L1 and other lineages (L2 and L4, n = 3) (p > 0.05). Further characterization using RACE (rapid amplification of cDNA ends), followed by target prediction with TargetRNA3 unveils that sRNA 1 (55 base pairs) targets Rv0506, Rv0697, and Rv3590c, and sRNA 29 (86 base pairs) targets Rv33859c, Rv3345c, Rv0755c, Rv0107c, Rv1817, Rv2950c, Rv1181, Rv3610c, and Rv3296. Functional characterization with Mycobrowser reveals these targets involved in regulating intermediary metabolism and respiration, cell wall and cell processes, lipid metabolism, information pathways, and PE/PPE. In summary, two novel sRNAs, sRNA 1 and sRNA 29, exhibited differential expression between L1 and other lineages, with predicted roles in essential Mtb functions. These findings offer insights into Mtb regulatory mechanisms, holding promise for the development of improved tuberculosis treatment strategies in the future.

Integrative and comparative genomic analyses of mammalian macrophage responses to intracellular mycobacterial pathogens

Mycobacterium tuberculosis, the causative agent of human tuberculosis (hTB), is a close evolutionary relative of Mycobacterium bovis, which causes bovine tuberculosis (bTB), one of the most damaging infectious diseases to livestock agriculture. Previous studies have shown that the pathogenesis of bTB disease is comparable to hTB disease, and that the bovine and human alveolar macrophage (bAM and hAM, respectively) transcriptomes are extensively reprogrammed in response to infection with these intracellular mycobacterial pathogens. In this study, a multi-omics integrative approach was applied with functional genomics and GWAS data sets across the two primary hosts (Bos taurus and Homo sapiens) and both pathogens (M. bovis and M. tuberculosis). Four different experimental infection groups were used: 1) bAM infected with M. bovis, 2) bAM infected with M. tuberculosis, 3) hAM infected with M. tuberculosis, and 4) human monocyte-derived macrophages (hMDM) infected with M. tuberculosis. RNA-seq data from these experiments 24 h post-infection (24 hpi) was analysed using three computational pipelines: 1) differentially expressed genes, 2) differential gene expression interaction networks, and 3) combined pathway analysis. The results were integrated with high-resolution bovine and human GWAS data sets to detect novel quantitative trait loci (QTLs) for resistance to mycobacterial infection and resilience to disease. This revealed common and unique response macrophage pathways for both pathogens and identified 32 genes (12 bovine and 20 human) significantly enriched for SNPs associated with disease resistance, the majority of which encode key components of the NF-κB signalling pathway and that also drive formation of the granuloma.

Mechanisms of host adaptation by bacterial pathogens

The emergence of new infectious diseases poses a major threat to humans, animals, and broader ecosystems. Defining factors that govern the ability of pathogens to adapt to new host species is therefore a crucial research imperative. Pathogenic bacteria are of particular concern, given dwindling treatment options amid the continued expansion of antimicrobial resistance. In this review, we summarize recent advancements in the understanding of bacterial host species adaptation, with an emphasis on pathogens of humans and related mammals. We focus particularly on molecular mechanisms underlying key steps of bacterial host adaptation including colonization, nutrient acquisition, and immune evasion, as well as suggest key areas for future investigation. By developing a greater understanding of the mechanisms of host adaptation in pathogenic bacteria, we may uncover new strategies to target these microbes for the treatment and prevention of infectious diseases in humans, animals, and the broader environment.

Regulatory role of Mycobacterium tuberculosis MtrA on dormancy/resuscitation revealed by a novel target gene-mining strategy

Introduction

The unique dormancy of Mycobacterium tuberculosis plays a significant role in the major clinical treatment challenge of tuberculosis, such as its long treatment cycle, antibiotic resistance, immune escape, and high latent infection rate.

Methods

To determine the function of MtrA, the only essential response regulator, one strategy was developed to establish its regulatory network according to high-quality genome-wide binding sites.

Results and discussion

The complex modulation mechanisms were implied by the strong bias distribution of MtrA binding sites in the noncoding regions, and 32.7% of the binding sites were located inside the target genes. The functions of 288 potential MtrA target genes predicted according to 294 confirmed binding sites were highly diverse, and DNA replication and damage repair, lipid metabolism, cell wall component biosynthesis, cell wall assembly, and cell division were the predominant pathways. Among the 53 pathways shared between dormancy/resuscitation and persistence, which accounted for 81.5% and 93.0% of the total number of pathways, respectively, MtrA regulatory genes were identified not only in 73.6% of their mutual pathways, but also in 75.4% of the pathways related to dormancy/resuscitation and persistence respectively. These results suggested the pivotal roles of MtrA in regulating dormancy/resuscitation and the apparent relationship between dormancy/resuscitation and persistence. Furthermore, the finding that 32.6% of the MtrA regulons were essential in vivo and/or in vitro for M. tuberculosis provided new insight into its indispensability. The findings mentioned above indicated that MtrA is a novel promising therapeutic target for tuberculosis treatment since the crucial function of MtrA may be a point of weakness for M. tuberculosis.

Structural and functional characterization of mycobacterial PhoH2 and identification of potential inhibitor of its enzymatic activity

Mycobacterium tuberculosis is composed of a cumbersome signaling and protein network which partakes in bacterial survival and augments its pathogenesis. Mycobacterial PhoH2 (Mt-PhoH2) is a signaling element and a predictive phosphate starvation protein that works in an ATP-dependent manner. Here, we elaborated the characterization of Mt-PhoH2 through biophysical, biochemical, and computational methods. In addition to its intrinsic ATPase activity, the biochemical experiments revealed its GTPase activity and both activities are metal ion dependent. Magnesium, manganese, copper, iron, nickel, zinc, cesium, calcium, and lithium were examined for their effect on activity, and the optimum activity was found with 10 mM of Mg2+ ions. The kinetic parameters of 3 µM Mt-PhoH2 were observed as Km 4.873 ± 0.44 µM, Vmax 12.3817 ± 0.084 µM/min/mg, Kcat 0.0075 ± 0.00005 s-1, and Kcat/Km 0.0015 ± 0.000001 µM-1 s-1 with GTP. In the case of GTP as a substrate, a 20% decrease in enzymatic activity and a 50% increase in binding affinity of Mt-PhoH2 were observed. The substrates ADP and GDP inhibit the ATPase and GTPase activity of Mt-PhoH2. CD spectroscopy showed the dominance of alpha helix in the secondary structure of Mt-PhoH2, and this structural pattern was altered upon addition of ATP and GTP. In silico inhibitor screening revealed ML141 and NAV_2729 as two potential inhibitors of the catalytic activity of Mt-PhoH2. Mt-PhoH2 is essential for mycobacterial growth as its knockdown strain showed a decreased growth effect. Overall, the present article emphasizes the factors essential for the proper functioning of Mt-PhoH2 which is a participant in the toxin-antitoxin machinery and may also play an important role in phosphate starvation.

Origin and dispersal of the Mycobacterium tuberculosis Haarlem genotype: Clues from its phylogeographic landscape and human migration

The Haarlem family belongs to the Euro-American phylogenetic lineage of Mycobacterium tuberculosis and is one of the globally spread genotypes of this important human pathogen. In spite of the sporadic observations on drug resistance and peculiar virulence profile, Haarlem remains in the shade of other M. tuberculosis genotypes. I analyzed genotyping data of the Haarlem genotype in light of its pathogenic properties and relevant human migration, to gain insight into its origin, evolutionary history, and current spread. Central Europe is marked with a very high prevalence of both major Haarlem subclades ancestral H3/SIT50 and derived H1, jointly making 33-41% in Czechia, Austria, and Hungary. There is a declining gradient of Haarlem beyond central Europe with 10-18% in Italy, France, Belgium, 10-13% in the Balkan countries and Turkey. Placing the available genetic diversity and ancient DNA data within the historical context, I hypothesize that M. tuberculosis Haarlem genotype likely originated in Central Europe and its primary long-term circulation occurred within the area of the former Austria/Austria-Hungary Empire in the 14th-19th centuries. The genotype is not highly transmissible and its spread was driven by long-term human migration. The European colonial expansion (when accompanied by a sufficient volume of migration) was a vehicle of its secondary dissemination. I conclude that human migration and its lack thereof (but not strain pathobiology) was a major driving force that shaped the population structure of this global lineage of M. tuberculosis. At the same time, Haarlem strains appear over-represented in some ethnic groups which warrants in-depth experimental research.

Drug resistance and epidemiological success of modern Mycobacterium tuberculosis lineages in western India

BACKGROUND

Transmission is contributing to the slow decline of tuberculosis (TB) incidence globally. Drivers of TB transmission in India, the country estimated to carry a quarter of the World's burden, are not well studied. We conducted a genomic epidemiology study to compare epidemiological success, host factors and drug resistance (DR) among the four major Mycobacterium tuberculosis (Mtb) lineages (L1-4) circulating in Pune, India.

METHODS

We performed whole-genome sequencing (WGS) of Mtb sputum culture-positive isolates from participants in two prospective cohort studies and predicted genotypic susceptibility using a validated random forest model. We used maximum likelihood estimation to build phylogenies. We compared lineage specific phylogenetic and time-scaled metrics to assess epidemiological success.

RESULTS

Of the 642 isolates that underwent WGS, 612 met sequence quality criteria. Most isolates belonged to L3 (44.6%). The majority (61.1%) of multidrug-resistant isolates belonged to L2 (P < 0.001). In molecular dating, L2 demonstrated a higher rate and more recent resistance acquisition. We measured higher clustering, and time-scaled haplotypic density (THD) for L4 and L2 compared to L3 and/or L1 suggesting higher epidemiological success. L4 demonstrated higher THD and clustering (OR 5.1 (95% CI 2.3-12.3) in multivariate models controlling for host factors and DR.

CONCLUSION

L2 shows a higher frequency of DR and both L2 and L4 demonstrate evidence of higher epidemiological success than L3 or L1 in the study setting. Our findings highlight the need for contact tracing around TB cases, and heightened surveillance of TB DR in India.

High-throughput nanopore targeted sequencing for efficient drug resistance assay of Mycobacterium tuberculosis

Drug-resistant tuberculosis (TB), especially multidrug-resistant tuberculosis (MDR-TB) and extensively drug-resistant tuberculosis (XDR-TB), is one of the urgent clinical problems and public health challenges. Culture-based phenotypic drug susceptibility testing (pDST) is time-consuming, and PCR-based assays are limited to hotspot mutations. In this study, we developed and validated a convenient and efficient approach based on high-throughput nanopore sequencing technology combined with multiplex PCR, namely nanopore targeted sequencing (NTS), to simultaneously sequence 18 genes associated with antibiotic resistance in Mycobacterium tuberculosis (MTB). The analytical performance of NTS was evaluated, and 99 clinical samples were collected to assess its clinical performance. The NTS results showed that MTB and its drug resistance were successfully identified in approximately 7.5 h. Furthermore, compared to the pDST and Xpert MTB/RIF assays, NTS provided much more drug resistance information, covering 14 anti-TB drugs, and it identified 20 clinical cases of drug-resistant MTB. The mutations underlying these drug-resistant cases were all verified using Sanger sequencing. Our approach for this TB drug resistance assay offers several advantages, including being culture-free, efficient, high-throughput, and highly accurate, which would be very helpful for clinical patient management and TB infection control.

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.

Machine learning-based drug design for identification of thymidylate kinase inhibitors as a potential anti-Mycobacterium tuberculosis

The rise of antibiotic-resistant Mycobacterium tuberculosis (Mtb) has reduced the availability of medications for tuberculosis therapy, resulting in increased morbidity and mortality globally. Tuberculosis spreads from the lungs to other parts of the body, including the brain and spine. Developing a single drug can take several decades, making drug discovery costly and time-consuming. Machine learning algorithms like support vector machines (SVM), k-nearest neighbor (k-NN), random forest (RF) and Gaussian naive base (GNB) are fast and effective and are commonly used in drug discovery. These algorithms are ideal for the virtual screening of large compound libraries to classify molecules as active or inactive. For the training of the models, a dataset of 307 was downloaded from BindingDB. Among 307 compounds, 85 compounds were labeled as active, having an IC50 below 58 mM, while 222 compounds were labeled inactive against thymidylate kinase, with 87.2% accuracy. The developed models were subjected to an external ZINC dataset of 136,564 compounds. Furthermore, we performed the 100-ns dynamic simulation and post trajectories analysis of compounds having good interaction and score in molecular docking. As compared to the standard reference compound, the top three hits revealed greater stability and compactness. In conclusion, our predicted hits can inhibit thymidylate kinase overexpression to combat Mycobacterium tuberculosis.Communicated by Ramaswamy H. Sarma.

Identification of positively selected genes in Mycobacterium tuberculosis from southern Xinjiang Uygur autonomous region of China

Background

Tuberculosis (TB), mainly caused by Mycobacterium tuberculosis (Mtb), remains a serious public health problem. Increasing evidence supports that selective evolution is an important force affecting genomic determinants of Mtb phenotypes. It is necessary to further understand the Mtb selective evolution and identify the positively selected genes that probably drive the phenotype of Mtb.

Methods

This study mainly focused on the positive selection of 807 Mtb strains from Southern Xinjiang of China using whole genome sequencing (WGS). PAML software was used for identifying the genes and sites under positive selection in 807 Mtb strains.

Results

Lineage 2 (62.70%) strains were the dominant strains in this area, followed by lineage 3 (19.45%) and lineage 4 (17.84%) strains. There were 239 codons in 47 genes under positive selection, and the genes were majorly associated with the functions of transcription, defense mechanisms, and cell wall/membrane/envelope biogenesis. There were 28 codons (43 mutations) in eight genes (gyrA, rpoB, rpoC, katG, pncA, embB, gid, and cut1) under positive selection in multi-drug resistance (MDR) strains but not in drug-susceptible (DS) strains, in which 27 mutations were drug-resistant loci, 9 mutations were non-drug-resistant loci but were in drug-resistant genes, 2 mutations were compensatory mutations, and 5 mutations were in unknown drug-resistant gene of cut1. There was a codon in Rv0336 under positive selection in L3 strains but not in L2 and L4 strains. The epitopes of T and B cells were both hyper-conserved, particularly in the T-cell epitopes.

Conclusion

This study revealed the ongoing selective evolution of Mtb. We found some special genes and sites under positive selection which may contribute to the advantage of MDR and L3 strains. It is necessary to further study these mutations to understand their impact on phenotypes for providing more useful information to develop new TB interventions.

Pharmacogenetic Study of Drugs Affecting Mycobacterium tuberculosis

BACKGROUND

Pharmacogenetic research has led to significant progress in understanding how genetic factors influence drug response in tuberculosis (TB) treatment. One ongoing challenge is the variable occurrence of adverse drug reactions in some TB patients. Previous studies have indicated that genetic variations in the N-acetyltransferase 2 (NAT2) and solute carrier organic anion transporter family member 1B1 (SLCO1B1) genes can impact the blood concentrations of the first-line anti-TB drugs isoniazid (INH) and rifampicin (RIF), respectively. This study aimed to investigate the influence of pharmacogenetic markers in the NAT2 and SLCO1B1 genes on TB treatment outcomes using whole-exome sequencing (WES) analysis.

METHODS

DNA samples were collected from 30 healthy Iranian adults aged 18-40 years. The allelic frequencies of single-nucleotide polymorphisms (SNPs) in the NAT2 and SLCO1B1 genes were determined through WES.

RESULTS

Seven frequent SNPs were identified in the NAT2 gene (rs1041983, rs1801280, rs1799929, rs1799930, rs1208, rs1799931, rs2552), along with 16 frequent SNPs in the SLCO1B1 gene (rs2306283, rs11045818, rs11045819, rs4149056, rs4149057, rs2291075, rs201722521, rs11045852, rs11045854, rs756393362, rs11045859, rs74064211, rs201556175, rs34671512, rs71581985, rs4149085).

CONCLUSION

Genetic variations in NAT2 and SLCO1B1 can affect the metabolism of INH and RIF, respectively. A better understanding of the pharmacogenetic profile in the study population may facilitate the design of more personalized and effective TB treatment strategies. Further research is needed to directly correlate these genetic markers with clinical outcomes in TB patients.

Important role of DNA methylation hints at significant potential in tuberculosis

Tuberculosis (TB), an infectious disease caused by Mycobacterium tuberculosis (Mtb) infection, has persisted as a major global public health threat for millennia. Until now, TB continues to challenge efforts aimed at controlling it, with drug resistance and latent infections being the two main factors hindering treatment efficacy. The scientific community is still striving to understand the underlying mechanisms behind Mtb's drug resistance and latent infection. DNA methylation, a critical epigenetic modification occurring throughout an individual's growth and development, has gained attention following advances in high-throughput sequencing technologies. Researchers have observed abnormal DNA methylation patterns in the host genome during Mtb infection. Given the escalating issue of drug-resistant Mtb, delving into the role of DNA methylation in TB's development is crucial. This review article explores DNA methylation's significance in human growth, development and disease, and its role in regulating Mtb's evolution and infection processes. Additionally, it discusses potential applications of DNA methylation research in tuberculosis.

Refined understanding of the impact of the Mycobacterium tuberculosis complex diversity on the intrinsic susceptibility to pretomanid

Previous work reported unprecedented differences in the intrinsic in vitro susceptibility of the Mycobacterium tuberculosis complex (MTBC) to pretomanid (Pa) using the Mycobacteria Growth Indicator Tube (MGIT) system. We tested 125 phylogenetically diverse strains from all known MTBC lineages (1-9) without known Pa resistance mutations and four strains with known resistance mutations as controls. This confirmed that MTBC, unlike most bacteria-antimicrobial combinations, displayed substantial differences in the intrinsic susceptibility relative to the technical variation of Pa MIC testing. This was also the case for the Middlebrook 7H11 (7H11) medium, demonstrating that these differences were not specific to MGIT. Notably, lineage 1 was confirmed to have intrinsically elevated MICs compared with lineages 2, 3, 4, and 7 (L2-4/7), underlining the urgent need for WHO to publish its decision of whether lineage 1 should be deemed treatable by BPaL(M), the now preferred all-oral regimen for treating rifampin-resistant tuberculosis. Lineages 5 and 6, which are most frequent in West Africa, responded differently to Pa, with lineage 5 being more similar to L2-4/7 and lineage 6 being more susceptible. More data are needed to determine whether 7H11 MICs are systematically lower than those in MGIT.

IMPORTANCE

This study confirmed that the Mycobacterium tuberculosis complex lineage 1, responsible for 28% of global tuberculosis cases, is less susceptible to pretomanid (Pa). It also refined the understanding of the intrinsic susceptibilities of lineages 5 and 6, most frequent in West Africa, and lineages 8 and 9. Regulators must review whether these in vitro differences affect the clinical efficacy of the WHO-recommended BPaL(M) regimen and set breakpoints for antimicrobial susceptibility testing accordingly. Notably, regulators should provide detailed justifications for their decisions to facilitate public scrutiny.

Transmission dynamics and phylogeography of Mycobacterium tuberculosis in China based on whole-genome phylogenetic analysis

OBJECTIVES

This study aimed to describe the lineage-specific transmissibility and epidemiological migration of Mycobacterium tuberculosis in China.

METHODS

We curated a large set of whole-genome sequences from 3204 M. tuberculosis isolates, including thousands of newly sequenced genomes, and applied a series of metrics to compare the transmissibility of M. tuberculosis strains between lineages and sublineages. The countrywide transmission patterns of major lineages were explored.

RESULTS

We found that lineage 2 (L2) was the most prevalent lineage in China (85.7%), with the major sublineage 2.2.1 (80.9%), followed by lineage 4 (L4) (13.8%), which comprises major sublineages 4.2 (1.5%), 4.4 (6.2%) and 4.5 (5.8%). We showed evidence for frequent cross-regional spread and large cluster formation of L2.2.1 strains, whereas L4 strains were relatively geographically restricted in China. Next, we applied a series of genomic indices to evaluate M. tuberculosis strain transmissibility and uncovered higher transmissibility of L2.2.1 compared with the L2.2.2 and L4 sublineages. Phylogeographic analysis showed that southern, eastern, and northern China were highly connected regions for countrywide L2.2.1 strain spread.

CONCLUSIONS

The present study provides insights into the different transmission and migration patterns of the major M. tuberculosis lineages in China and highlights that transmissible L2.2.1 is a threat to tuberculosis control.

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.

Metabolites from Streptomyces aureus (VTCC43181) and Their Inhibition of Mycobacterium tuberculosis ClpC1 Protein

Tuberculosis is one of the most common infectious diseases in the world, caused by Mycobacterium tuberculosis. The outbreak of multiple drug-resistant tuberculosis has become a major challenge to prevent this disease worldwide. ClpC1 is a Clp ATPase protein of Mycobacterium tuberculosis, functioning as a chaperon when combined with the Clp complex. ClpC1 has emerged as a new target to discover anti-tuberculosis drugs. This study aimed to explore the ClpC1 inhibitors from actinomycetes, which have been known to provide abundant sources of antibiotics. Two cyclic peptides, including nocardamin (1), halolitoralin A (3), and a lactone pleurone (2), were isolated from the culture of Streptomyces aureus (VTCC43181). The structures of these compounds were determined based on the detailed analysis of their spectral data and comparison with references. This is the first time these compounds have been isolated from S. aureus. Compounds 1-3 were evaluated for their affection of ATPase activity of the recombinant ClpC1 protein. Of these compounds, halolitoralin A (1), a macrocyclic peptide, was effective for the ATPase hydrolysis of the ClpC1 protein.

Ongoing evolution of the Mycobacterium tuberculosis lactate dehydrogenase reveals the pleiotropic effects of bacterial adaption to host pressure

The bacterial determinants that facilitate Mycobacterium tuberculosis (Mtb) adaptation to the human host environment are poorly characterized. We have sought to decipher the pressures facing the bacterium in vivo by assessing Mtb genes that are under positive selection in clinical isolates. One of the strongest targets of selection in the Mtb genome is lldD2, which encodes a quinone-dependent L-lactate dehydrogenase (LldD2) that catalyzes the oxidation of lactate to pyruvate. Lactate accumulation is a salient feature of the intracellular environment during infection and lldD2 is essential for Mtb growth in macrophages. We determined the extent of lldD2 variation across a set of global clinical isolates and defined how prevalent mutations modulate Mtb fitness. We show the stepwise nature of lldD2 evolution that occurs as a result of ongoing lldD2 selection in the background of ancestral lineage-defining mutations and demonstrate that the genetic evolution of lldD2 additively augments Mtb growth in lactate. Using quinone-dependent antibiotic susceptibility as a functional reporter, we also find that the evolved lldD2 mutations functionally increase the quinone-dependent activity of LldD2. Using 13C-lactate metabolic flux tracing, we find that lldD2 is necessary for robust incorporation of lactate into central carbon metabolism. In the absence of lldD2, label preferentially accumulates in dihydroxyacetone phosphate (DHAP) and glyceraldehyde-3-phosphate (G3P) and is associated with a discernible growth defect, providing experimental evidence for accrued lactate toxicity via the deleterious buildup of sugar phosphates. The evolved lldD2 variants increase lactate incorporation to pyruvate while altering triose phosphate flux, suggesting both an anaplerotic and detoxification benefit to lldD2 evolution. We further show that the mycobacterial cell is transcriptionally sensitive to the changes associated with altered lldD2 activity which affect the expression of genes involved in cell wall lipid metabolism and the ESX- 1 virulence system. Together, these data illustrate a multifunctional role of LldD2 that provides context for the selective advantage of lldD2 mutations in adapting to host stress.

Phylogenetic Survival Analysis

Going back in time through a phylogenetic tree makes it possible to evaluate ancestral genomes and assess their potential to acquire key polymorphisms of interest over evolutionary time. Knowledge of this kind may allow for the emergence of key traits to be predicted and pre-empted from currently circulating strains in the future. Here, we present a novel genome-wide survival analysis and use the emergence of drug resistance in Mycobacterium tuberculosis as an example to demonstrate the potential and utility of the technique.

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.

Towards the reconstruction of a global TB history using a new pipeline "TB-Annotator"

Mycobacterium tuberculosis complex (MTBC) has a population structure consisting of 9 human and animal lineages. The genomic diversity within these lineages is a pathogenesis factor that affects virulence, transmissibility, host response, and antibiotic resistance. Hence it is important to develop improved information systems for tracking and understanding the spreading and evolution of genomes. We present results obtained thanks to a new informatics platform for computational biology of MTBC, that uses a convenience sample from public/private SRAs, designated as TB-Annotator. Version 1 was a first interactive graphic-based web tool based on 15,901 representative genomes. Version 2, still interactive, is a more sophisticated database, developed using the Snakemake Workflow Management System (WMS) that allows an unsupervised global and scalable analysis of the content of the USA National Center for Biotechnology Information Short Read Archives database. This platform analyzes nucleotide variants, the presence/absence of genes, known regions of difference and detect new deletions, the insertion sites of mobile genetic elements, and allows phylogenetic trees to be built, imported in a graphical interface and interactively analyzed between the data and the tree. The objective of TB-Annotator is triple: detect recent epidemiological links, reconstruct distant phylogeographical histories as well as perform more complex phenotypic/genotypic Genome-Wide Association Studies (GWAS). In this paper, we compare the various taxonomic SNPs-based labels and hierarchies previously described in recent reference papers for L1, and present a comparative analysis that allows identification of alias and thus provides the basis of a future unifying naming scheme for L1 sublineages. We present a global phylogenetic tree built with RAxML-NG, and one on L2; at the time of writing, we characterized about 200 sublineages, with many new ones; a detail tree for Modern L2 and a hierarchical scheme allowing to facilitate L2 lineage assignment are also presented.

Tuberculosis in mummies - New findings, perspectives and limitations

The molecular analysis of ancient pathogen DNA represents a unique opportunity for the study of infectious diseases in ancient human remains. Among other diseases, paleogenetic studies have been successful in detecting tuberculous DNA in ancient human remains. In the beginning of ancient DNA (aDNA) studies, the presence of tuberculosis (TB) DNA was assessed using a PCR-based assay targeting specific regions of the Mycobacterium tuberculosis (MTB) complex, such as the repetitive element IS6110. The advent of high-throughput sequencing has enabled the reconstruction of full ancient TB genomes in the field of paleomicrobiology. However, despite the numerous paleopathological and PCR-based studies on the presence of tuberculosis in historic human remains, full genome wide reconstructions are still limited to well-preserved specimens with low environmental contamination and connected with extensive screening efforts. This has led to some controversies regarding the evolutionary history of its causative agent Mycobacterium tuberculosis. In this context, mummies have been shown to be a good source for the detection of MTB complex DNA due to a low exposure to environmental influences and the overall good state of preservation of hard and soft tissues in the human remains. Here, we present the major findings on the presence of TB infections in the 18th century naturally mummified human remains from Vác, Hungary and the current status of the detection of MTB complex DNA in mummified human remains. The future perspectives of detecting tuberculosis in mummies will be discussed in the light of methodological aspects, as well as ethical and curational challenges.

Infectious disease in the Pleistocene: Old friends or old foes?

The impact of endemic and epidemic disease on humans has traditionally been seen as a comparatively recent historical phenomenon associated with the Neolithisation of human groups, an increase in population size led by sedentarism, and increasing contact with domesticated animals as well as species occupying opportunistic symbiotic and ectosymbiotic relationships with humans. The orthodox approach is that Neolithisation created the conditions for increasing population size able to support a reservoir of infectious disease sufficient to act as selective pressure. This orthodoxy is the result of an overly simplistic reliance on skeletal data assuming that no skeletal lesions equated to a healthy individual, underpinned by the assumption that hunter-gatherer groups were inherently healthy while agricultural groups acted as infectious disease reservoirs. The work of van Blerkom, Am. J. Phys. Anthropol., vol. suppl 37 (2003), Wolfe et al., Nature, vol. 447 (2007) and Houldcroft and Underdown, Am. J. Phys. Anthropol., vol. 160, (2016) has changed this landscape by arguing that humans and pathogens have long been fellow travelers. The package of infectious diseases experienced by our ancient ancestors may not be as dissimilar to modern infectious diseases as was once believed. The importance of DNA, from ancient and modern sources, to the study of the antiquity of infectious disease, and its role as a selective pressure cannot be overstated. Here we consider evidence of ancient epidemic and endemic infectious diseases with inferences from modern and ancient human and hominin DNA, and from circulating and extinct pathogen genomes. We argue that the pandemics of the past are a vital tool to unlock the weapons needed to fight pandemics of the future.

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.

New Meroterpenoid Derivatives from the Pomegranate-Derived Endophytic Fungus Talaromyces purpureogenus

In this study, we report the isolation of two new meroterpenoids, miniolutelide D (1) and miniolutelide E (13-epi-miniolutelide C) (2), along with two meroterpenoidal analogues (3 and 4) and two phenolic compounds (5 and 6) from the endophytic fungus Talaromyces purpureogenus derived from Punica granatum fruits. Their structures were elucidated using extensive MS, 1D, and 2D NMR spectroscopic analyses as well as by comparing with data in the literature. The absolute configurations of 1 and 2 were determined using TDDFT-ECD calculations. Antimicrobial activity was evaluated. Compound 5 displayed significant activity against methicillin-resistant Staphylococcus aureus strain ATCC 700699 and moderate activity against S. aureus strain ATCC 29213.

A smooth tubercle bacillus from Ethiopia phylogenetically close to the Mycobacterium tuberculosis complex

The Mycobacterium tuberculosis complex (MTBC) includes several human- and animal-adapted pathogens. It is thought to have originated in East Africa from a recombinogenic Mycobacterium canettii-like ancestral pool. Here, we describe the discovery of a clinical tuberculosis strain isolated in Ethiopia that shares archetypal phenotypic and genomic features of M. canettii strains, but represents a phylogenetic branch much closer to the MTBC clade than to the M. canettii strains. Analysis of genomic traces of horizontal gene transfer in this isolate and previously identified M. canettii strains indicates a persistent albeit decreased recombinogenic lifestyle near the emergence of the MTBC. Our findings support that the MTBC emergence from its putative free-living M. canettii-like progenitor is evolutionarily very recent, and suggest the existence of a continuum of further extant derivatives from ancestral stages, close to the root of the MTBC, along the Great Rift Valley.

Prolonged survival of a patient with active MDR-TB HIV co-morbidity: insights from a Mycobacterium tuberculosis strain with a unique genomic deletion

Coinfection of HIV and multidrug-resistant tuberculosis (MDR-TB) presents significant challenges in terms of the treatment and prognosis of tuberculosis, leading to complexities in managing the disease and impacting the overall outcome for TB patients. This study presents a remarkable case of a patient with MDR-TB and HIV coinfection who survived for over 8 years, despite poor treatment adherence and comorbidities. Whole genome sequencing (WGS) of the infecting Mycobacterium tuberculosis (Mtb) strain revealed a unique genomic deletion, spanning 18 genes, including key genes involved in hypoxia response, intracellular survival, immunodominant antigens, and dormancy. This deletion, that we have called "Del-X," potentially exerts a profound influence on the bacterial physiology and its virulence. Only few similar deletions were detected in other non-related Mtb genomes worldwide. In vivo evolution analysis identified drug resistance and metabolic adaptation mutations and their temporal dynamics during the patient's treatment course.

Mycobacterial DnaQ is an Alternative Proofreader Ensuring DNA Replication Fidelity

Remove of mis-incorporated nucleotides ensures replicative fidelity. Although the ε-exonuclease DnaQ is a well-established proofreader in the model organism Escherichia coli, proofreading in mycobacteria relies on the polymerase and histidinol phosphatase (PHP) domain of replicative polymerase despite the presence of an alternative DnaQ homolog. Here, we show that depletion of DnaQ in Mycolicibacterium smegmatis results in increased mutation rate, leading to AT-biased mutagenesis and elevated insertions/deletions in homopolymer tract. We demonstrated that mycobacterial DnaQ binds to the b-clamp and functions synergistically with the PHP domain to correct replication errors. Further, we found that the mycobacterial DnaQ sustains replicative fidelity upon chromosome topological stress. Intriguingly, we showed that a naturally evolved DnaQ variant prevalent in clinical Mycobacterium tuberculosis isolates enables hypermutability and is associated with extensive drug resistance. These results collectively establish that the alternative DnaQ functions in proofreading, and thus reveal that mycobacteria deploy two proofreaders to maintain replicative fidelity.

Mycobacterium tuberculosis Sub-Lineage 4.2.2/SIT149 as Dominant Drug-Resistant Clade in Northwest Ethiopia 2020-2022: In-silico Whole-Genome Sequence Analysis

Introduction

Drug resistance (DR) in Mycobacterium tuberculosis complex (MTBC) is mainly associated with certain lineages and varies across regions and countries. The Beijing genotype is the leading resistant lineage in Asia and western countries. M. tuberculosis (Mtb) (sub) lineages responsible for most drug resistance in Ethiopia are not well described. Hence, this study aimed to identify the leading drug resistance sub-lineages and characterize first-line anti-tuberculosis drug resistance-associated single nucleotide polymorphisms (SNPs).

Methods

A facility-based cross-sectional study was conducted in 2020-2022 among new and presumptive multidrug resistant-TB (MDR-TB) cases in Northwest Ethiopia. Whole-genome sequencing (WGS) was performed on 161 isolates using Illumina NovaSeq 6000 technology. The SNP mutations associated with drug resistance were identified using MtbSeq and TB profiler Bioinformatics softwares.

Results

Of the 146 Mtb isolates that were successfully genotyped, 20 (13.7%) harbored one or more resistance-associated SNPs. L4.2.2.ETH was the leading drug-resistant sub-lineage, accounting for 10/20 (50%) of the resistant Mtb. MDR-TB isolates showed extensive mutations against first-line anti-TB drugs. Ser450Leu/(tcg/tTg) for Rifampicin (RIF), Ser315Thr/(agc/aCc) for Isoniazid (INH), Met306Ile/(atg/atA(C)) for Ethambutol (EMB), and Gly69Asp for Streptomycin (STR) were the leading resistance associated mutations which accounted for 56.5%, 89.5%, 47%, and 29.4%, respectively. The presence of both clustered and non-clustered drug resistance (DR) isolates indicated that the epidemics is driven by both new DR development and acquired resistance.

Conclusion

The high prevalence of drug-resistant TB due to geographically restricted sub-lineages (L4.2.2.ETH) indicates the ongoing local micro epidemics. The Mtb drug resistance surveillance system must be improved. Further evolutionary analysis of L4.2.2.ETH strain is highly desirable to understand evolutionary forces that leads L4.2.2.ETH in to high level DR and transmissible sub-lineage.

The heterogenous and diverse population of prophages in Mycobacterium genomes

IMPORTANCE

Mycobacterium species include several human pathogens and mycobacteriophages show potential for therapeutic use to control Mycobacterium infections. However, phage infection profiles vary greatly among Mycobacterium abscessus clinical isolates and phage therapies must be personalized for individual patients. Mycobacterium phage susceptibility is likely determined primarily by accessory parts of bacterial genomes, and we have identified the prophage and phage-related genomic regions across sequenced Mycobacterium strains. The prophages are numerous and diverse, especially in M. abscessus genomes, and provide a potentially rich reservoir of new viruses that can be propagated lytically and used to expand the repertoire of therapeutically useful phages.

Recent progress in the genotyping of bovine tuberculosis and its rapid diagnosis via nanoparticle-based electrochemical biosensors

Bovine tuberculosis (bTB) is considered a worldwide infectious zoonotic disease. Mycobacterium bovis causes bTB disease. It is one of the Mycobacterium tuberculosis complex (MTBC) members. MTBC is a clonal complex of close relatives with approximately 99.95% similarity. M. bovis is a spillover pathogen that can transmit from animals to humans and rarely from humans to animals with contact. Genotyping techniques are important to discriminate and differentiate between MTBC species. Spoligotyping and mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) are widely used but they have some limitations. As an alternative, whole genome sequencing approaches have been utilized due to their high-resolution power. They are employed in typing M. bovis and explain the evolutionary and phylogenetic relationships between isolates. The control of bTB disease has attracted a large amount of attention. Rapid and proper diagnosis is necessary for monitoring the disease as an initial step for its control and treatment. Nanotechnology has a potential impact on the rapid diagnosis and treatment of bTB through the use of nanocarrier and metal nanoparticles (NPs). Special attention has been paid to voltammetric and impedimetric electrochemical strategies as facile, sensitive, and selective methods for the efficient detection of tuberculosis. The efficacy of these sensors is enhanced in the presence of NPs, which act as recognition and/or redox probes. Gold, silver, copper, cobalt, graphene, and magnetic NPs, as well as polypyrrole nanowires and multiwalled carbon nanotubes have been employed for detecting tuberculosis. Overall, NP-based electrochemical sensors represent a promising tool for the diagnosis of bTB.

Plasmid-mediated drug resistance in mycobacteria: the tip of the iceberg?

Macrolides, such as clarithromycin, are crucial in the treatment of nontuberculous mycobacteria (NTM). NTM are notoriously innately drug resistant, which has made the dependence on macrolides for their treatment even more important. Not surprisingly, resistance to macrolides has been documented in some NTM, including Mycobacterium avium and Mycobacterium abscessus, which are the two NTM species most often identified in clinical isolates. Resistance is mediated by point mutations in the 23S ribosomal RNA or by methylation of the rRNA by a methylase (encoded by an erm gene). Chromosomally encoded erm genes have been identified in many of the macrolide-resistant isolates, but not in Mycobacterium chelonae. Now, Brown-Elliott et al. (J Clin Microbiol 61:e00428-23, 2023, https://doi.org/10.1128/JCM.00428-23) describe the identification of a new erm variant, erm(55), which was found either on the chromosome or on a plasmid in highly macrolide-resistant clinical isolates of M. chelonae. The chromosomal erm(55) gene appears to be associated with mobile elements; one gene is within a putative transposon and the second is in a large (37 kb) insertion/deletion. The plasmid carrying erm(55) also encodes type IV and type VII secretion systems, which are often linked on large mycobacterial plasmids and are hypothesized to mediate plasmid transfer. While the conjugative transfer of the erm(55)-containing plasmid between NTM has yet to be demonstrated, the inferences are clear, as evidenced by the dissemination of plasmid-mediated drug resistance in other medically important bacteria. Here, we discuss the findings of Brown-Elliott et al., and the potential ramifications on treatment of NTM infections.

Genomic, phenotypic and demographic characterization of Mycobacterium tuberculosis in Israel in 2021

According to World Health Organization WHO, Tuberculosis (TB) is the second cause of death from infectious disease worldwide. During 2021, 10.6 million people were infected with TB, and 1.6 million people died. TB is caused by pathogens belonging to the Mycobacterium tuberculosis complex (MTBC), mainly by Mycobacterium tuberculosis (M.tb). Members of this complex are acid-fast bacilli, which can cause intrapulmonary and extra pulmonary TB, and can be divided into various lineages, based on genomic markers. The main public health threat comes from drug resistant M.tb strains, which are responsible for about 25% of TB death and treatment failure worldwide. Treating drug resistant TB patients significantly raises the costs of TB treatment. This study provides an overview of the demographic and drug susceptibility characteristics of newly diagnosed TB patients in Israel in 2021. The State of Israel has a very low level of TB endemicity and is at a pre-elimination phase. Notably, only 11.7% of the newly diagnosed TB patients were born in Israel. In this report, of the 154 new laboratory-confirmed TB patients, 66.7% had pulmonary TB, while 16% had extrapulmonary TB. Males accounted for 52% of the patients, with the most prevalent age group being 21-40. Most patients were citizens of Israel (53.9%), while 37.7% had no Israeli citizenship. Among non-citizens, there was a predominance of males and patients aged 21-40. The susceptibility profile showed a high resistance rate to streptomycin (18.2%) and to a lower extent to isoniazid (13.6%), pyrazinamide (8.4%), rifampicin (7.8%), and ethambutol (3.2%). Only 2 cases of XDR-TB and 10 MDR-TB strains were detected in Israel in 2021, with both XDR strains and 5 out of 10 MDR strains belonging to the Beijing lineage. Most of Beijing isolates were resistant to at least one tested drug. Genomic sequencing of 134 out of 156 strains and bioinformatics analysis using the MTBseq program and WHO mutation catalogue shows a good match with only 9 discrepancies between phenotypic and genotypic susceptibility profiles in first line drugs. The most common lineage is Delhi-Cas (23%) followed by the Beijing lineage (17%). Most patients from the Delhi-Cas lineage were born in Africa, while patients with Beijing isolates were born in different countries. Minimum spanning tree analysis identified 15 clusters. The study highlights the need for ongoing surveillance of TB using molecular and phenotypic tools to further decreasing the spreading level of the disease and develop effective treatment strategies.

Phylogenetic distribution of malonate semialdehyde decarboxylase (MSAD) genes among strains within the genus Mycobacterium: evidence of MSAD gene loss in the evolution of pathogenic mycobacteria

Despite the great diversity of malonate semialdehyde decarboxylases (MSADs), one of five subgroups of the tautomerase superfamily (TSF) found throughout the biosphere, their distribution among strains within the genus Mycobacterium remains unknown. In this study, we sought to investigate the phylogenetic distribution of MSAD genes of mycobacterial species via genome analysis of 192 different reference Mycobacterium species or subspecies retrieved from NCBI databases. We found that in a total of 87 of 192 strains (45.3%), MSAD-1 and MSAD-2 were distributed in an exclusive manner among Mycobacterium species except for 12 strains, including Mycobacterium chelonae members, with both in their genome. Of note, Mycobacterium strains better adapted to the host and of high virulence potential, such as the Mycobacterium tuberculosis complex, Mycobacterium leprae, Mycobacterium marinum, Mycobacterium ulcerans, and Mycobacterium avium subsp. paratuberculosis, had no orthologs of MSAD in their genome, suggesting MSAD loss during species differentiation in pathogenic slow-growing Mycobacterium. To investigate the MSAD distribution among strains of M. avium subspecies, the genome sequences of a total of 255 reference strains from the four subspecies of M. avium (43 of subspecies avium, 162 of subspecies hominissuis, 49 of subspecies paratuberculosis, and 1 of subspecies silvaticum) were further analyzed. We found that only 121 of 255 strains (47.4%) had MSADs in their genome, with none of the 49 M. avium subsp. paratuberculosis strains having MSAD genes. Even in 13 of 121 M. avium strains with the MSAD-1 gene in their genome, deletion mutations in the 98th codon causing premature termination of MSAD were found, further highlighting the occurrence of MSAD pseudogenization during species or subspecies differentiation of M. avium. In conclusion, our data indicated that there are two distinct types of MSADs, MSAD-1 and MSAD-2, among strains in the Mycobacterium genus, but more than half of the strains, including pathogenic mycobacteria, M. tuberculosis and M. leprae, have no orthologs in their genome, suggesting MSAD loss during host adaptation of pathogenic mycobacteria. In the future, the role of two distinct MSADs, MSAD-1 and MSAD-2, in mycobacterial pathogenesis or evolution should be investigated.

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.

Ongoing evolution of the Mycobacterium tuberculosis lactate dehydrogenase reveals the pleiotropic effects of bacterial adaption to host pressure

The bacterial determinants that facilitate Mycobacterium tuberculosis (Mtb) adaptation to the human host environment are poorly characterized. We have sought to decipher the pressures facing the bacterium in vivo by assessing Mtb genes that are under positive selection in clinical isolates. One of the strongest targets of selection in the Mtb genome is lldD2 , which encodes a quinone-dependent L-lactate dehydrogenase (LldD2) that catalyzes the oxidation of lactate to pyruvate. Lactate accumulation is a salient feature of the intracellular environment during infection and lldD2 is essential for Mtb growth in macrophages. We determined the extent of lldD2 variation across a set of global clinical isolates and defined how prevalent mutations modulates Mtb fitness. We show the stepwise nature of lldD2 evolution that occurs as a result of ongoing lldD2 selection in the background of ancestral lineage defining mutations and demonstrate that the genetic evolution of lldD2 additively augments Mtb growth in lactate. Using quinone-dependent antibiotic susceptibility as a functional reporter, we also find that the evolved lldD2 mutations functionally increase the quinone-dependent activity of LldD2. Using 13 C-lactate metabolic flux tracing, we find that lldD2 is necessary for robust incorporation of lactate into central carbon metabolism. In the absence of lldD2 , label preferentially accumulates in methylglyoxal precursors dihydroxyacetone phosphate (DHAP) and glyceraldehyde-3-phosphate (G3P) and is associated with a discernible growth defect, providing experimental evidence for accumulated lactate toxicity via a methylglyoxal pathway that has been proposed previously. The evolved lldD2 variants increase lactate incorporation to pyruvate but also alter flux in the methylglyoxal pathway, suggesting both an anaplerotic and detoxification benefit to lldD2 evolution. We further show that the mycobacterial cell is transcriptionally sensitive to the changes associated with altered lldD2 activity which affect the expression of genes involved in cell wall lipid metabolism and the ESX-1 virulence system. Together, these data illustrate a multifunctional role of LldD2 that provide context for the selective advantage of lldD2 mutations in adapting to host stress.

Mycobacteriaceae Phenome Atlas (MPA): A Standardized Atlas for the Mycobacteriaceae Phenome Based on Heterogeneous Sources

The bacterial family Mycobacteriaceae includes pathogenic and nonpathogenic bacteria, and systematic research on their genome and phenome can give comprehensive perspectives for exploring their disease mechanism. In this study, the phenotypes of Mycobacteriaceae were inferred from available phenomic data, and 82 microbial phenotypic traits were recruited as data elements of the microbial phenome. This Mycobacteriaceae phenome contains five categories and 20 subcategories of polyphasic phenotypes, and three categories and eight subcategories of functional phenotypes, all of which are complementary to the existing data standards of microbial phenotypes. The phenomic data of Mycobacteriaceae strains were compiled by literature mining, third-party database integration, and bioinformatics annotation. The phenotypes were searchable and comparable from the website of the Mycobacteriaceae Phenome Atlas (MPA, https://www.biosino.org/mpa/). A topological data analysis of MPA revealed the co-evolution between Mycobacterium tuberculosis and virulence factors, and uncovered potential pathogenicity-associated phenotypes. Two hundred and sixty potential pathogen-enriched pathways were found by Fisher's exact test. The application of MPA may provide novel insights into the pathogenicity mechanism and antimicrobial targets of Mycobacteriaceae.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-023-00101-5.

Mycobacterium tuberculosis: Pathogenesis and therapeutic targets

Tuberculosis (TB) remains a significant public health concern in the 21st century, especially due to drug resistance, coinfection with diseases like immunodeficiency syndrome (AIDS) and coronavirus disease 2019, and the lengthy and costly treatment protocols. In this review, we summarize the pathogenesis of TB infection, therapeutic targets, and corresponding modulators, including first-line medications, current clinical trial drugs and molecules in preclinical assessment. Understanding the mechanisms of Mycobacterium tuberculosis (Mtb) infection and important biological targets can lead to innovative treatments. While most antitubercular agents target pathogen-related processes, host-directed therapy (HDT) modalities addressing immune defense, survival mechanisms, and immunopathology also hold promise. Mtb's adaptation to the human host involves manipulating host cellular mechanisms, and HDT aims to disrupt this manipulation to enhance treatment effectiveness. Our review provides valuable insights for future anti-TB drug development efforts.

Harnessing host-pathogen interactions for innovative drug discovery and host-directed therapeutics to tackle tuberculosis

Tuberculosis (TB) is a highly contagious bacterial infection caused by Mycobacterium tuberculosis (Mtb), which has been ranked as the second leading cause of death worldwide from a single infectious agent. As an intracellular pathogen, Mtb has well adapted to the phagocytic host microenvironment, influencing diverse host processes such as gene expression, trafficking, metabolism, and signaling pathways of the host to its advantage. These responses are the result of dynamic interactions of the bacteria with the host cell signaling pathways, whereby the bacteria attenuate the host cellular processes for their survival. Specific host genes and the mechanisms involved in the entry and subsequent stabilization of M. tuberculosis intracellularly have been identified in various genetic and chemical screens recently. The present understanding of the co-evolution of Mtb and macrophage system presented us the new possibilities for exploring host-directed therapeutics (HDT). Here, we discuss the host-pathogen interaction for Mtb, including the pathways adapted by Mtb to escape immunity. The review sheds light on different host-directed therapies (HDTs) such as repurposed drugs and vitamins, along with their targets such as granuloma, autophagy, extracellular matrix, lipids, and cytokines, among others. The article also examines the available clinical data on these drug molecules. In conclusion, the review presents a perspective on the current knowledge in the field of HDTs and the need for additional research to overcome the challenges associated HDTs.

A global atlas of marine antibiotic resistance genes and their expression

Oceans serve as global reservoirs of antibiotic-resistant bacteria and antibiotic resistance genes (ARGs). However, little is known about the traits and expression of ARGs in response to environmental factors. We analyzed 347 metagenomes and 182 metatranscriptomes to determine the distribution, hosts, and expression of ARGs in oceans. Our study found that the diversity and abundance of ARGs varied with latitude and depth. The core marine resistome mainly conferred glycopeptide and multidrug resistance. The hosts of this resistome were mainly limited to the core marine microbiome, with phylogenetic barriers to the horizontal transfer of ARGs, transfers being more frequent within species than between species. Sixty-five percent of the marine ARGs identified were expressed. More than 90% of high-risk ARGs were more likely to be expressed. Anthropogenic activity might affect the expression of ARGs by altering nitrate and phosphate concentrations and ocean temperature. Machine-learning models predict >97% of marine ARGs will change expression by 2100. High-risk ARGs will shift to low latitudes and regions with high anthropogenic activity, such as the Pacific and Atlantic Oceans. Certain ARGs serve a dual role in antibiotic resistance and potentially participate in element cycling, along with other unknown functions. Determining whether changes in ARG expression are beneficial to ecosystems and human health is challenging without comprehensive understanding of their functions. Our study identified a core resistome in the oceans and quantified the expression of ARGs for the development of future control strategies under global change.

Resistance patterns and transmission of mono- and polyresistant TB: clinical impact of WGS

Objectives

Rapidly diagnosing drug-resistant TB is crucial for improving treatment and transmission control. WGS is becoming increasingly accessible and has added value to the diagnosis and treatment of TB. The aim of the study was to perform WGS to determine the rate of false-positive results of phenotypic drug susceptibility testing (pDST) and characterize the molecular mechanisms of resistance and transmission of mono- and polyresistant Mycobacterium (M.) tuberculosis.

Methods

WGS was performed on 53 monoresistant and 25 polyresistant M. tuberculosis isolates characterized by pDST. Sequencing data were bioinformatically processed to infer mutations encoding resistance and determine the origin of resistance and phylogenetic relationship between isolates studied.

Results

The data showed the variable sensitivity and specificity of WGS in comparison with pDST as the gold standard: isoniazid 92.7% and 92.3%; streptomycin 41.9% and 100.0%; pyrazinamide 15% and 94.8%; and ethambutol 75.0% and 98.6%, respectively. We found novel mutations encoding resistance to streptomycin (in gidB) and pyrazinamide (in kefB). Most isolates belonged to lineage 4 (80.1%) and the overall clustering rate was 11.5%. We observed lineage-specific gene variations encoding resistance to streptomycin and pyrazinamide.

Conclusions

This study highlights the clinical potential of WGS in ruling out false-positive drug resistance following phenotypic or genetic drug testing, and recommend this technology together with the WHO catalogue in designing an optimal individualized treatment regimen and preventing the development of MDR TB. Our results suggest that resistance is primarily developed through spontaneous mutations or selective pressure.

Mycobacterium tuberculosis Rv3435c regulates inflammatory cytokines and promotes the intracellular survival of recombinant Mycobacteria

Mycobacterium tuberculosis is a pathogenic bacterium that is parasitic in macrophages and show high adaptation to the host's immune response. It can also trigger a complex immune response in the host. This relies on proteins encoded by a series of M. tuberculosis-encoded virulence genes. We found that the M. tuberculosis Rv3435c gene is highly conserved among pathogenic mycobacteria, and might be a virulence gene. To explore the gene function of Rv3435c, we used Mycobacterium smegmatis to construct a recombinant mycobacterium expressing Rv3435c heterologously. The results that Rv3435c is a cell wall-related protein that changes bacterial and colony morphology, inhibits the growth rate of recombinant mycobacteria, and enhances their resistance to various stresses. We also found that the fatty acid levels of the recombinant strain changed. Simultaneously, Rv3435c can inhibit the expression and secretion of inflammatory factors and host cell apoptosis, and enhance the survival of recombinant bacteria in macrophages. Experimental data indicated that Rv3435c might play an important role in Mycobacterium tuberculosis infection.

Host-directed therapy for bacterial infections -Modulation of the phagolysosome pathway

Bacterial infections still impose a significant burden on humanity, even though antimicrobial agents have long since been developed. In addition to individual severe infections, the f fatality rate of sepsis remains high, and the threat of antimicrobial-resistant bacteria grows with time, putting us at inferiority. Although tremendous resources have been devoted to the development of antimicrobial agents, we have yet to recover from the lost ground we have been driven into. Looking back at the evolution of treatment for cancer, which, like infectious diseases, has the similarity that host immunity eliminates the lesion, the development of drugs to eliminate the tumor itself has shifted from a single-minded focus on drug development to the establishment of a treatment strategy in which the de-suppression of host immunity is another pillar of treatment. In infectious diseases, on the other hand, the development of therapies that strengthen and support the immune system has only just begun. Among innate immunity, the first line of defense that bacteria encounter after invading the host, the molecular mechanisms of the phagolysosome pathway, which begins with phagocytosis to fusion with lysosome, have been elucidated in detail. Bacteria have a large number of strategies to escape and survive the pathway. Although the full picture is still unfathomable, the molecular mechanisms have been elucidated for some of them, providing sufficient clues for intervention. In this article, we review the host defense mechanisms and bacterial evasion mechanisms and discuss the possibility of host-directed therapy for bacterial infection by intervening in the phagolysosome pathway.

In silico studies of alkaloids and their derivatives against N-acetyltransferase EIS protein from Mycobacterium tuberculosis

Antibiotic resistance against Mycobacterium tuberculosis (M.tb.) has been a significant cause of death worldwide. The Enhanced intracellular survival (EIS) protein of the bacteria is an acetyltransferase that multiacetylates aminoglycoside antibiotics, preventing them from binding to the bacterial ribosome. To overcome the EIS-mediated antibiotics resistance of M.tb., we compiled 888 alkaloids and derivatives from five different databases and virtually screened them against the EIS receptor. The compound library was filtered down to 87 compounds, which underwent additional analysis and filtration. Moreover, the top 15 most prominent phytocompounds were obtained after the drug-likeness prediction and ADMET screening. Out of 15, nine compounds confirmed the maximum number of hydrogen bond interactions and reliable binding energies during molecular docking. Additionally, the Molecular dynamics (MD) simulation of nine compounds showed the three most stable complexes, further verified by re-docking with mutated protein. The density functional theory (DFT) calculation was performed to identify the HOMO-LUMO energy gaps of the selected three potential compounds. Finally, our selected top lead compounds i.e., Alkaloid AQC2 (PubChem85634496), Nobilisitine A (ChEbi68116), and N-methylcheilanthifoline (ChEbi140673) demonstrated more favourable outcomes when compared with reference compounds (i.e., 39b and 2i) in all parameters used in this study. Therefore, we anticipate that our findings will help to explore and develop natural compound therapy against multi and extensively drug-resistant strains of M.tb.Communicated by Ramaswamy H. Sarma.