A Comprehensive Map of Mycobacterium tuberculosis Complex Regions of Difference

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.

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.

Diversification of gene content in the Mycobacterium tuberculosis complex is determined by phylogenetic and ecological signatures

We analyzed the pan-genome and gene content modulation of the most diverse genome data set of the Mycobacterium tuberculosis complex (MTBC) gathered to date. The closed pan-genome of the MTBC was characterized by reduced accessory and strain-specific genomes, compatible with its clonal nature. However, significantly fewer gene families were shared between MTBC genomes as their phylogenetic distance increased. This effect was only observed in inter-species comparisons, not within-species, which suggests that species-specific ecological characteristics are associated with changes in gene content. Gene loss, resulting from genomic deletions and pseudogenization, was found to drive the variation in gene content. This gene erosion differed among MTBC species and lineages, even within M. tuberculosis, where L2 showed more gene loss than L4. We also show that phylogenetic proximity is not always a good proxy for gene content relatedness in the MTBC, as the gene repertoire of Mycobacterium africanum L6 deviated from its expected phylogenetic niche conservatism. Gene disruptions of virulence factors, represented by pseudogene annotations, are mostly not conserved, being poor predictors of MTBC ecotypes. Each MTBC ecotype carries its own accessory genome, likely influenced by distinct selective pressures such as host and geography. It is important to investigate how gene loss confer new adaptive traits to MTBC strains; the detected heterogeneous gene loss poses a significant challenge in elucidating genetic factors responsible for the diverse phenotypes observed in the MTBC. By detailing specific gene losses, our study serves as a resource for researchers studying the MTBC phenotypes and their immune evasion strategies.IMPORTANCEIn this study, we analyzed the gene content of different ecotypes of the Mycobacterium tuberculosis complex (MTBC), the pathogens of tuberculosis. We found that changes in their gene content are associated with their ecological features, such as host preference. Gene loss was identified as the primary driver of these changes, which can vary even among different strains of the same ecotype. Our study also revealed that the gene content relatedness of these bacteria does not always mirror their evolutionary relationships. In addition, some genes of virulence can be variably lost among strains of the same MTBC ecotype, likely helping them to evade the immune system. Overall, our study highlights the importance of understanding how gene loss can lead to new adaptations in these bacteria and how different selective pressures may influence their genetic makeup.

Whole genome sequencing and comparative genomics of Mycobacterium orygis isolated from different animal hosts to identify specific diagnostic markers

Introduction

Mycobacterium orygis, a member of MTBC has been identified in higher numbers in the recent years from animals of South Asia. Comparative genomics of this important zoonotic pathogen is not available which can provide data on the molecular difference between other MTBC members. Hence, the present study was carried out to isolate, whole genome sequence M. orygis from different animal species (cattle, buffalo and deer) and to identify molecular marker for the differentiation of M. orygis from other MTBC members.

Methods

Isolation and whole genome sequencing of M. orygis was carried out for 9 samples (4 cattle, 4 deer and 1 buffalo) died due to tuberculosis. Comparative genomics employing 53 genomes (44 from database and 9 newly sequenced) was performed to identify SNPs, spoligotype, pangenome structure, and region of difference.

Results

M. orygis was isolated from water buffalo and sambar deer which is the first of its kind report worldwide. Comparative pangenomics of all M. orygis strains worldwide (n= 53) showed a closed pangenome structure which is also reported for the first time. Pairwise SNP between TANUVAS_2, TANUVAS_4, TANUVAS_5, TANUVAS_7 and NIRTAH144 was less than 15 indicating that the same M. orygis strain may be the cause for infection. Region of difference prediction showed absence of RD7, RD8, RD9, RD10, RD12, RD301, RD315 in all the M. orygis analyzed. SNPs in virulence gene, PE35 was found to be unique to M. orygis which can be used as marker for identification.

Conclusion

The present study is yet another supportive evidence that M. orygis is more prevalent among animals in South Asia and the zoonotic potential of this organism needs to be evaluated.

A hybridization target enrichment approach for pathogen genomics

IMPORTANCE

Emerging infectious diseases require continuous pathogen monitoring. Rapid clinical diagnosis by nucleic acid amplification is limited to a small number of targets and may miss target detection due to new mutations in clinical isolates. Whole-genome sequencing (WGS) identifies genome-wide variations that may be used to determine a pathogen's drug resistance patterns and phylogenetically characterize isolates to track disease origin and transmission. WGS is typically performed using DNA isolated from cultured clinical isolates. Culturing clinical specimens increases turn-around time and may not be possible for fastidious bacteria. To overcome some of these limitations, direct sequencing of clinical specimens has been attempted using expensive capture probes to enrich the entire genomes of target pathogens. We present a method to produce a cost-effective, time-efficient, and large-scale synthesis of probes for whole-genome enrichment. We envision that our method can be used for direct clinical sequencing of a wide range of microbial pathogens for genomic epidemiology.

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.

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.

Mapping the phylogeny and lineage history of geographically distinct BCG vaccine strains

The bacillus Calmette-Guérin (BCG) vaccine has been in use for prevention of tuberculosis for over a century. It remains the only widely available tuberculosis vaccine and its protective efficacy has varied across geographical regions. Since it was developed, the BCG vaccine strain has been shared across different laboratories around the world, where use of differing culture methods has resulted in genetically distinct strains over time. Whilst differing BCG vaccine efficacy around the world is well documented, and the reasons for this may be multifactorial, it has been hypothesized that genetic differences in BCG vaccine strains contribute to this variation. Isolates from an historic archive of lyophilized BCG strains were regrown, DNA was extracted and then whole-genome sequenced using Oxford Nanopore Technologies. The resulting whole-genome data were plotted on a phylogenetic tree and analysed to identify the presence or absence of regions of difference (RDs) and single-nucleotide polymorphisms (SNPs) relating to virulence, growth and cell wall structure. Of 50 strains available, 36 were revived in culture and 39 were sequenced. Morphology differed between the strains distributed before and after 1934. There was phylogenetic association amongst certain geographically classified strains, most notably BCG-Russia, BCG-Japan and BCG-Danish. RD2, RD171 and RD713 deletions were associated with late strains (seeded after 1927). When mapped to BCG-Pasteur 1172, the SNPs in sigK, plaA, mmaA3 and eccC5 were associated with early strains. Whilst BCG-Russia, BCG-Japan and BCG-Danish showed strong geographical isolate clustering, the late strains, including BCG-Pasteur, showed more variation. A wide range of SNPs were seen within geographically classified strains, and as much intra-strain variation as between-strain variation was seen. The date of distribution from the original Pasteur laboratory (early pre-1927 or late post-1927) gave the strongest association with genetic differences in regions of difference and virulence-related SNPs, which agrees with the previous literature.

Insight into pathogenomics and phylogeography of hypervirulent and highly-lethal Mycobacterium tuberculosis strain cluster

BACKGROUND

. The Mycobacterium tuberculosis Beijing genotype is globally spread lineage with important medical properties that however vary among its subtypes. M. tuberculosis Beijing 14717-15-cluster was recently discovered as both multidrug-resistant, hypervirulent, and highly-lethal strain circulating in the Far Eastern region of Russia. Here, we aimed to analyze its pathogenomic features and phylogeographic pattern.

RESULTS

. The study collection included M. tuberculosis DNA collected between 1996 and 2020 in different world regions. The bacterial DNA was subjected to genotyping and whole genome sequencing followed by bioinformatics and phylogenetic analysis. The PCR-based assay to detect specific SNPs of the Beijing 14717-15-cluster was developed and used for its screening in the global collections. Phylogenomic and phylogeographic analysis confirmed endemic prevalence of the Beijing 14717-15-cluster in the Asian part of Russia, and distant common ancestor with isolates from Korea (> 115 SNPs). The Beijing 14717-15-cluster isolates had two common resistance mutations RpsL Lys88Arg and KatG Ser315Thr and belonged to spoligotype SIT269. The Russian isolates of this cluster were from the Asian Russia while 4 isolates were from the Netherlands and Spain. The cluster-specific SNPs that significantly affect the protein function were identified in silico in genes within different categories (lipid metabolism, regulatory proteins, intermediary metabolism and respiration, PE/PPE, cell wall and cell processes).

CONCLUSIONS

. We developed a simple method based on real-time PCR to detect clinically significant MDR and hypervirulent Beijing 14717-15-cluster. Most of the identified cluster-specific mutations were previously unreported and could potentially be associated with increased pathogenic properties of this hypervirulent M. tuberculosis strain. Further experimental study to assess the pathobiological role of these mutations is warranted.

Large genomic deletions delineate Mycobacterium tuberculosis L4 sublineages in South American countries

Mycobacterium tuberculosis (Mtb) is still one of the primary pathogens of humans causing tuberculosis (TB) disease. Mtb embraces nine well-defined phylogenetic lineages with biological and geographical disparities. The lineage L4 is the most globally widespread of all lineages and was introduced to America with European colonization. Taking advantage of many genome projects available in public repositories, we undertake an evolutionary and comparative genomic analysis of 522 L4 Latin American Mtb genomes. Initially, we performed careful quality control of public read datasets and applied several thresholds to filter out low-quality data. Using a genome de novo assembly strategy and phylogenomic methods, we spotted novel south American clades that have not been revealed yet. Additionally, we describe genomic deletion profiles of these strains from an evolutionary perspective and report Mycobacterium tuberculosis L4 sublineages signature-like gene deletions, some of the novel. One is a specific deletion of 6.5 kbp that is only present in sublineage 4.1.2.1. This deletion affects a complex group of 10 genes with putative products annotated, among others, as a lipoprotein, transmembrane protein, and toxin/antitoxin system proteins. The second novel deletion spans for 4.9 kbp and specific of a particular clade of the 4.8 sublineage and affects 7 genes. The last novel deletion affects 4 genes, extends for 4.8 kbp., and is specific to some strains within the 4.1.2.1 sublineage that are present in Colombia, Peru and Brasil.

Features of Mycobacterium bovis Complete Genomes Belonging to 5 Different Lineages

Mammalian tuberculosis (TB) is a zoonotic disease mainly due to Mycobacterium bovis (M. bovis). A current challenge for its eradication is understanding its transmission within multi-host systems. Improvements in long-read sequencing technologies have made it possible to obtain complete bacterial genomes that provide a comprehensive view of species-specific genomic features. In the context of TB, new genomic references based on complete genomes genetically close to field strains are also essential to perform precise field molecular epidemiological studies. A total of 10 M. bovis strains representing each genetic lineage identified in France and in other countries were selected for performing complete assembly of their genomes. Pangenome analysis revealed a "closed" pangenome composed of 3900 core genes and only 96 accessory genes. Whole genomes-based alignment using progressive Mauve showed remarkable conservation of the genomic synteny except that the genomes have a variable number of copies of IS6110. Characteristic genomic traits of each lineage were identified through the discovery of specific indels. Altogether, these results provide new genetic features that improve the description of M. bovis lineages. The availability of new complete representative genomes of M. bovis will be useful to epidemiological studies and better understand the transmission of this clonal-evolving pathogen.

Application of metagenomic next-generation sequencing in cutaneous tuberculosis

Tuberculous infection in a skin wound is a rare but well-known condition. This study describes a child infected with tuberculosis after being wounded. Because of swelling and pain in his wrist tissue, he was admitted to the Affiliated Hospital of Jining Medical University of Shandong Province on 16 October 2021. His medical history only included a wound. He was discharged after debridement. The laboratory data were normal. Two months after surgery, his wound was still swollen and painful. Secretions from the wound were sent for metagenomic next-generation sequencing (mNGS), which revealed three reads related to the Mycobacterium tuberculosis complex group (MTBC). A diagnosis of cutaneous tuberculosis (TB) was made. The wound disappeared after anti-TB drugs were administered. This case demonstrates that, while TB presenting as a severe cutaneous wound is rare, it should be considered in the clinical diagnosis. Clinicians should also pay attention to extrapulmonary infection with MTBC in patients, particularly in some long-suffering patients, and identify the specific pathogen as soon as possible. mNGS could help to identify pathogens and facilitate early treatment, thereby improving the prognosis.

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.

Unraveling the metabolism of Mycobacterium caprae using comparative genomics

In our laboratory, Mycobacterium caprae has poor growth in standard medium (SM) 7H9-OADC supplemented with pyruvate and Tween-80. Our objectives were to identify mutations affecting M. caprae metabolism and use this information to design a culture medium to improve its growth. We selected 77 M. caprae genomes and sequenced M. caprae NLA000201913 used in our experiments. Mutations present in >95% of the strains compared to Mycobacterium tuberculosis H37Rv were analyzed in silico for their deleterious effects on proteins of metabolic pathways. Apart from the known defect in the pyruvate kinase, M. caprae has important lesions in enzymes of the TCA cycle, methylmalonyl cycle, B₁₂ metabolism, and electron-transport chain. We provide evidence of enzymatic redundancy elimination and epistatic mutations, and possible production of toxic metabolites hindering M. caprae growth in vitro. A newly designed SM supplemented with l-glutamate allowed faster growth and increased final microbial mass of M. caprae. However, possible accumulation of metabolic waste-products and/or nutritional limitations halted M. caprae growth prior to a M. tuberculosis-like stationary phase. Our findings suggest that M. caprae relies on GABA and/or glyoxylate shunts for in vitro growth in routine media. The newly developed medium will improve experiments with this bacterium by allowing faster growth in vitro.

Mycobacterium tuberculosis sensu stricto in African Apes, What Is Its True Health Impact?

Since the Symposium on Mycobacterial Infections of Zoo Animals held at the National Zoological Park, Smithsonian Institution in 1976, our understanding of tuberculosis (TB) in non-domestic animals has greatly expanded. Throughout the past decades, this knowledge has resulted in improved zoo-habitats and facilities design, stricter biosecurity measures, and advanced diagnostic methods, including molecular techniques, that have significantly decreased the number of clinical disease caused by Mycobacterium tuberculosis in apes under human care settings. In the other hand, exponential growth of human populations has led to human encroachment in wildlife habitat which has resulted in increased inter-species contact and recurrent conflict between humans and wild animals. Although it is widely accepted that non-human primates are susceptible to M. tb infection, opinions differ with regard to the susceptibility to develop disease amongst different taxa. Specifically, some authors suggest that African apes are less susceptible to clinical tuberculosis than other species of primates. The aim of this review article is to evaluate the current scientific literature to determine the actual health impact of disease caused by Mycobacterium tuberculosis and more specifically Mycobacterium tuberculosis sensu stricto in African apes. The literature review included literature databases: Web of Science, Pubmed, Scopus, Wiley, Springer and Science direct, without temporal limit and proceedings of annual conferences in the field of wildlife health. Our general inclusion criteria included information about serological, molecular, pathological (macroscopic and/or microscopic), and clinical evidence of TB in African apes; while our, our more stringent inclusion selection criteria required that in addition to a gross pathology, a molecular test confirmed Mycobacterium tuberculosis sensu stricto as the cause of disease or death. We identified eleven reports of tuberculosis in African apes; of those, only four reports met the more stringent selection criteria that confirmed M. tb sensu stricto in six individuals. All reports that confirmed M. tb sensu stricto originated from zoological collections. Our review suggests that there is little evidence of disease or mortality caused by M. tb in the different species of African apes both under human care and free ranging populations. Additional studies are needed in free-ranging, semi-captive populations (sanctuaries) and animals under human care (zoos and rescue centers) to definitely conclude that this mycobacteria has a limited health effect in African ape species.