Mycobacterial Evolution Research Articles

Genome wide screening to discover novel toxin-antitoxin modules in Mycobacterium indicus pranii; perspective on gene acquisition during mycobacterial evolution.

Mycobacterium indicus pranii (MIP), a benign saprophyte with potent immunomodulatory attributes, holds a pivotal position in mycobacterial evolution, potentially serving as the precursor to the pathogenic Mycobacterium avium complex (MAC). Despite its established immunotherapeutic efficacy against leprosy and notable outcomes in gram-negative sepsis and COVID-19 cases, the genomic and biochemical features of MIP remain largely elusive. This study explores the uncharted territory of toxin-antitoxin (TA) systems within MIP, hypothesizing their role in mycobacterial pathogenicity regulation. Genome-wide screening, employing diverse databases, unveils putative TA modules in MIP, setting the stage for a comparative analysis with known modules in Mycobacterium tuberculosis, Mycobacterium smegmatis, Escherichia coli, and Vibrio cholerae. The study further delves into the TA network of MAC and Mycobacterium intracellulare, unraveling interactive properties and family characteristics of identified TA modules in MIP. This comprehensive exploration seeks to illuminate the contribution of TA modules in regulating virulence, habitat diversification, and the evolutionary pathogenicity of mycobacteria. The insights garnered from this investigation not only enhance our understanding of MIP's potential as a vaccine candidate but also hold promise in optimizing tuberculosis drug regimens for expedited recovery.

Cold Cas: reevaluating the occurrence of CRISPR/Cas systems in Mycobacteriaceae.

Bacterial CRISPR/Cas systems target foreign genetic elements such as phages and regulate gene expression by some pathogens, even in the host. The system is a marker for evolutionary history and has been used for inferences in Mycobacterium tuberculosis for 30 years. However, knowledge about mycobacterial CRISPR/Cas systems remains limited. It is believed that Type III-A Cas systems are exclusive to Mycobacterium canettii and the M. tuberculosis complex (MTBC) of organisms and that very few of the >200 diverse species of non-tuberculous mycobacteria (NTM) possess any CRISPR/Cas system. This study sought unreported CRISPR/Cas loci across NTM to better understand mycobacterial evolution, particularly in species phylogenetically near the MTBC. An analysis of available mycobacterial genomes revealed that Cas systems are widespread across Mycobacteriaceae and that some species contain multiple types. The phylogeny of Cas loci shows scattered presence in many NTM, with variation even within species, suggesting gains/losses of these loci occur frequently. Cas Type III-A systems were identified in pathogenic Mycobacterium heckeshornense and the geological environmental isolate Mycobacterium SM1. In summary, mycobacterial CRISPR/Cas systems are numerous, Type III-A systems are unreliable as markers for MTBC evolution, and mycobacterial horizontal gene transfer appears to be a frequent source of genetic variation.

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

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

Differential Sensitivity of Mycobacteria to Isoniazid Is Related to Differences in KatG-Mediated Enzymatic Activation of the Drug.

Isoniazid (INH) is a cornerstone of antitubercular therapy. Mycobacterium tuberculosis complex bacteria are the only mycobacteria sensitive to clinically relevant concentrations of INH. All other mycobacteria, including M. marinum and M. avium subsp. paratuberculosis are resistant. INH requires activation by bacterial KatG to inhibit mycobacterial growth. We tested the role of the differences between M. tuberculosis KatG and that of other mycobacteria in INH sensitivity. We cloned the M. boviskatG gene into M. marinum and M. avium subsp. paratuberculosis and measured the MIC of INH. We recombinantly expressed KatG of these mycobacteria and tested in vitro binding to, and activation of, INH. Introduction of katG from M. bovis into M. marinum and M. avium subsp. paratuberculosis rendered them 20 to 30 times more sensitive to INH. Analysis of different katG sequences across the genus found KatG evolution diverged from RNA polymerase-defined mycobacterial evolution. Biophysical and biochemical tests of M. bovis and nontuberculous mycobacteria (NTM) KatG proteins showed lower affinity to INH and substantially lower enzymatic capacity for the conversion of INH into the active form in NTM. The KatG proteins of M. marinum and M. avium subsp. paratuberculosis are substantially less effective in INH activation than that of M. tuberculosis, explaining the relative INH insensitivity of these microbes. These data indicate that the M. tuberculosis complex KatG is divergent from the KatG of NTM, with a reciprocal relationship between resistance to host defenses and INH resistance. Studies of bacteria where KatG is functionally active but does not activate INH may aid in understanding M. tuberculosis INH-resistance mechanisms, and suggest paths to overcome them.

Comparative genomics of Mycobacterium mucogenicum and Mycobacterium neoaurum clade members emphasizing tRNA and non-coding RNA

BackgroundMycobacteria occupy various ecological niches and can be isolated from soil, tap water and ground water. Several cause diseases in humans and animals. To get deeper insight into our understanding of mycobacterial evolution focusing on tRNA and non-coding (nc)RNA, we conducted a comparative genome analysis of Mycobacterium mucogenicum (Mmuc) and Mycobacterium neoaurum (Mneo) clade members.ResultsGenome sizes for Mmuc- and Mneo-clade members vary between 5.4 and 6.5 Mbps with the complete MmucT (type strain) genome encompassing 6.1 Mbp. The number of tRNA genes range between 46 and 79 (including one pseudo tRNA gene) with 39 tRNA genes common among the members of these clades, while additional tRNA genes were probably acquired through horizontal gene transfer. Selected tRNAs and ncRNAs (RNase P RNA, tmRNA, 4.5S RNA, Ms1 RNA and 6C RNA) are expressed, and the levels for several of these are higher in stationary phase compared to exponentially growing cells. The rare tRNAIleTAT isoacceptor and two for mycobacteria novel ncRNAs: the Lactobacillales-derived GOLLD RNA and a homolog to the antisense Salmonella typhimurium phage Sar RNA, were shown to be present and expressed in certain Mmuc-clade members.ConclusionsPhages, IS elements, horizontally transferred tRNA gene clusters, and phage-derived ncRNAs appears to have influenced the evolution of the Mmuc- and Mneo-clades. While the number of predicted coding sequences correlates with genome size, the number of tRNA coding genes does not. The majority of the tRNA genes in mycobacteria are transcribed mainly from single genes and the levels of certain ncRNAs, including RNase P RNA (essential for the processing of tRNAs), are higher at stationary phase compared to exponentially growing cells. We provide supporting evidence that Ms1 RNA represents a mycobacterial 6S RNA variant. The evolutionary routes for the ncRNAs RNase P RNA, tmRNA and Ms1 RNA are different from that of the core genes.

Mycobacterial Evolution Intersects With Host Tolerance.

Over the past 200 years, tuberculosis (TB) has caused more deaths than any other infectious disease, likely infecting more people than it has at any other time in human history. Mycobacterium tuberculosis (Mtb), the etiologic agent of TB, is an obligate human pathogen that has evolved through the millennia to become an archetypal human-adapted pathogen. This review focuses on the evolutionary framework by which Mtb emerged as a specialized human pathogen and applies this perspective to the emergence of specific lineages that drive global TB burden. We consider how evolutionary pressures, including transmission dynamics, host tolerance, and human population patterns, may have shaped the evolution of diverse mycobacterial genomes.

Epidemiological and genomic determinants of tuberculosis outbreaks in First Nations communities in Canada

BackgroundIn Canada, tuberculosis disproportionately affects foreign-born and First Nations populations. Within First Nations’ peoples, a high proportion of cases occur in association with outbreaks. Tuberculosis transmission in the context of outbreaks is thought to result from the convergence of several factors including characteristics of the cases, contacts, the environment, and the pathogen.MethodsWe examined the epidemiological and genomic determinants of two well-characterized tuberculosis outbreaks attributed to two super-spreaders among First Nations in the province of Alberta. These outbreaks were associated with two distinct DNA fingerprints (restriction fragment-length polymorphisms or RFLPs 0.0142 and 0.0728). We compared outbreak isolates with endemic isolates not spatio-temporarily linked to outbreak cases. We extracted epidemiological variables pertaining to tuberculosis cases and contacts from individual public health records and the provincial tuberculosis registry. We conducted group analyses using parametric and non-parametric statistical tests. We carried out whole-genome sequencing and bioinformatic analysis using validated protocols.ResultsWe observed differences between outbreak and endemic groups in the mean number of total and child-aged contacts and the number of contacts with new positive and converted tuberculin skin tests in all group comparisons (p < 0.05). Differences were also detected in the proportion of cases with cavitation on a chest radiograph and the mean number of close contacts in selected group comparisons (p < 0.02). A phylogenetic network analysis of whole-genome sequencing data indicated that most outbreak and endemic strains were closely related to the source case for the 0.0142 fingerprint. For the 0.0728 fingerprint, the source case haplotype was circulating among endemic cases prior to the outbreak. Genetic and temporal distances were not correlated for either RFLP 0.0142 (r2 = − 0.05) or RFLP 0.0728 (r2 = 0.09) when all isolates were analyzed.ConclusionsWe found no evidence that endemic strains acquired mutations resulting in their emergence in outbreak form. We conclude that the propagation of these outbreaks was likely driven by the combination of characteristics of the source cases, contacts, and the environment. The role of whole-genome sequencing in understanding mycobacterial evolution and in assisting public health authorities in conducting contact investigations and managing outbreaks is important and expected to grow in the future.

ESX secretion systems: mycobacterial evolution to counter host immunity.

Mycobacterium tuberculosis uses sophisticated secretion systems, named 6 kDa early secretory antigenic target (ESAT6) protein family secretion (ESX) systems (also known as type VII secretion systems), to export a set of effector proteins that helps the pathogen to resist or evade the host immune response. Since the discovery of the esx loci during the M. tuberculosis H37Rv genome project, structural biology, cell biology and evolutionary analyses have advanced our knowledge of the function of these systems. In this Review, we highlight the intriguing roles that these studies have revealed for ESX systems in bacterial survival and pathogenicity during infection with M. tuberculosis. Furthermore, we discuss the diversity of ESX systems that has been described among mycobacteria and selected non-mycobacterial species. Finally, we consider how our knowledge of ESX systems might be applied to the development of novel strategies for the treatment and prevention of disease.

Identification of a Non-Pentapeptide Region Associated with Rapid Mycobacterial Evolution

A large portion of the coding capacity of Mycobacterium tuberculosis is devoted to the production of proteins containing several copies of the pentapeptide-2 repeat, namely the PE/PPE_MPTR proteins. Protein domain repeats have a variety of binding properties and are involved in protein-protein interactions as well as binding to other ligands such as DNA and RNA. They are not as common in prokaryotes, compared to eukaryotes, but the enrichment of pentapeptide-2 repeats in Mycobacteria constitutes an exception to that rule. The genes encoding the PE/PPE_MPTR proteins have undergone many rearrangements and here we have identified the expansion patterns across the Mycobacteria. We have performed a reclassification of the PE/PPE_MPTR proteins using cohesive regions rather than sparse domain architectures. It is clear that these proteins have undergone large insertions of several pentapeptide-2 domains appearing adjacent to one another in a repetitive pattern. Further, we have identified a non-pentapeptide motif associated with rapid mycobacterial evolution. The sequence composition of this region suggests a different structure compared to pentapeptide-2 repeats. By studying the evolution of the PE/PPE_MPTR proteins, we have distinguished features pertaining to tuberculosis-inducing species. Further studies of the non-pentapeptide region associated with repeat expansions promises to shed light on the pathogenicity of Mycobacterium tuberculosis.

Mycobacterial Pan-Genome Analysis Suggests Important Role of Plasmids in the Radiation of Type VII Secretion Systems.

In mycobacteria, various type VII secretion systems corresponding to different ESX (ESAT-6 secretory) types, are contributing to pathogenicity, iron acquisition, and/or conjugation. In addition to the known chromosomal ESX loci, the existence of plasmid-encoded ESX systems was recently reported. To investigate the potential role of ESX-encoding plasmids on mycobacterial evolution, we analyzed a large representative collection of mycobacterial genomes, including both chromosomal and plasmid-borne sequences. Data obtained for chromosomal ESX loci confirmed the previous five classical ESX types and identified a novel mycobacterial ESX-4-like type, termed ESX-4-bis. Moreover, analysis of the plasmid-encoded ESX loci showed extensive diversification, with at least seven new ESX profiles, identified. Three of them (ESX-P clusters 1–3) were found in multiple plasmids, while four corresponded to singletons. Our phylogenetic and gene-order-analyses revealed two main groups of ESX types: 1) ancestral types, including ESX-4 and ESX-4-like systems from mycobacterial and non-mycobacterial actinobacteria and 2) mycobacteria-specific ESX systems, including ESX-1-2-3-5 systems and the plasmid-encoded ESX types. Synteny analysis revealed that ESX-P systems are part of phylogenetic groups that derived from a common ancestor, which diversified and resulted in the different ESX types through extensive gene rearrangements. A converging body of evidence, derived from composition bias-, phylogenetic-, and synteny analyses points to a scenario in which ESX-encoding plasmids have been a major driving force for acquisition and diversification of type VII systems in mycobacteria, which likely played (and possibly still play) important roles in the adaptation to new environments and hosts during evolution of mycobacterial pathogenesis.

Unraveling Mycobacterium tuberculosis genomic diversity and evolution in Lisbon, Portugal, a highly drug resistant setting.

BackgroundMultidrug- (MDR) and extensively drug resistant (XDR) tuberculosis (TB) presents a challenge to disease control and elimination goals. In Lisbon, Portugal, specific and successful XDR-TB strains have been found in circulation for almost two decades.ResultsIn the present study we have genotyped and sequenced the genomes of 56 Mycobacterium tuberculosis isolates recovered mostly from Lisbon. The genotyping data revealed three major clusters associated with MDR-TB, two of which are associated with XDR-TB. Whilst the genomic data contributed to elucidate the phylogenetic positioning of circulating MDR-TB strains, showing a high predominance of a single SNP cluster group 5. Furthermore, a genome-wide phylogeny analysis from these strains, together with 19 publicly available genomes of Mycobacterium tuberculosis clinical isolates, revealed two major clades responsible for M/XDR-TB in the region: Lisboa3 and Q1 (LAM).The data presented by this study yielded insights on microevolution and identification of novel compensatory mutations associated with rifampicin resistance in rpoB and rpoC. The screening for other structural variations revealed putative clade-defining variants. One deletion in PPE41, found among Lisboa3 isolates, is proposed to contribute to immune evasion and as a selective advantage. Insertion sequence (IS) mapping has also demonstrated the role of IS6110 as a major driver in mycobacterial evolution by affecting gene integrity and regulation.ConclusionsGlobally, this study contributes with novel genome-wide phylogenetic data and has led to the identification of new genomic variants that support the notion of a growing genomic diversity facing both setting and host adaptation.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-991) contains supplementary material, which is available to authorized users.

Distributive Conjugal Transfer: New Insights into Horizontal Gene Transfer and Genetic Exchange in Mycobacteria.

The last decade has seen an explosion in the application of genomic tools across all biological disciplines. This is also true for mycobacteria, where whole genome sequences are now available for pathogens and non-pathogens alike. Genomes within the Mycobacterium tuberculosis Complex (MTBC) bear the hallmarks of horizontal gene transfer (HGT). Conjugation is the form of HGT with the highest potential capacity and evolutionary influence. Donor and recipient strains of Mycobacterium smegmatis actively conjugate upon co-culturing in biofilms and on solid media. Whole genome sequencing of the transconjugant progeny demonstrated the incredible scale and range of genomic variation that conjugation generates. Transconjugant genomes are complex mosaics of the parental strains. Some transconjugant genomes are up to one-quarter donor-derived, distributed over 30 segments. Transferred segments range from ~50 bp to ~225,000 bp in length, and are exchanged with their recipient orthologs all around the genome. This unpredictable genome-wide infusion of DNA sequences is called Distributive Conjugal Transfer (DCT), to distinguish it from traditional oriT-based conjugation. The mosaicism generated in a single transfer event resembles that seen from meiotic recombination in sexually reproducing organisms, and contrasts with traditional models of HGT. This similarity allowed the application of a GWAS-like approach to map the donor genes that confer a donor mating identity phenotype. The mating identity genes map to the esx1 locus, expanding the central role of ESX-1 function in conjugation. The potential for DCT to instantaneously blend genomes will affect how we view mycobacterial evolution, and provide new tools for the facile manipulation of mycobacterial genomes.

Mapping of Genotype–Phenotype Diversity among Clinical Isolates of Mycobacterium tuberculosis by Sequence-Based Transcriptional Profiling

Genome sequencing has identified an extensive repertoire of single nucleotide polymorphisms among clinical isolates of Mycobacterium tuberculosis, but the extent to which these differences influence phenotypic properties of the bacteria remains to be elucidated. To determine whether these polymorphisms give rise to phenotypic diversity, we have integrated genome data sets with RNA sequencing to assess their impact on the comparative transcriptome profiles of strains belonging to M. tuberculosis Lineages 1 and 2. We observed clear correlations between genotype and transcriptional phenotype. These arose by three mechanisms. First, lineage-specific changes in amino acid sequence of transcriptional regulators were associated with alterations in their ability to control gene expression. Second, changes in nucleotide sequence were associated with alteration of promoter activity and generation of novel transcriptional start sites in intergenic regions and within coding sequences. We show that in some cases this mechanism is expected to generate functionally active truncated proteins involved in innate immune recognition. Finally, genes showing lineage-specific patterns of differential expression not linked directly to primary mutations were characterized by a striking overrepresentation of toxin–antitoxin pairs. Taken together, these findings advance our understanding of mycobacterial evolution, contribute to a systems level understanding of this important human pathogen, and more broadly demonstrate the application of state-of-the-art techniques to provide novel insight into mechanisms by which intergenic and silent mutations contribute to diversity.

Comparative genomic and phylogenetic approaches to characterize the role of genetic recombination in mycobacterial evolution.

The genus Mycobacterium encompasses over one hundred named species of environmental and pathogenic organisms, including the causative agents of devastating human diseases such as tuberculosis and leprosy. The success of these human pathogens is due in part to their ability to rapidly adapt to their changing environment and host. Recombination is the fastest way for bacterial genomes to acquire genetic material, but conflicting results about the extent of recombination in the genus Mycobacterium have been reported. We examined a data set comprising 18 distinct strains from 13 named species for evidence of recombination. Genomic regions common to all strains (accounting for 10% to 22% of the full genomes of all examined species) were aligned and concatenated in the chromosomal order of one mycobacterial reference species. The concatenated sequence was screened for evidence of recombination using a variety of statistical methods, with each proposed event evaluated by comparing maximum-likelihood phylogenies of the recombinant section with the non-recombinant portion of the dataset. Incongruent phylogenies were identified by comparing the site-wise log-likelihoods of each tree using multiple tests. We also used a phylogenomic approach to identify genes that may have been acquired through horizontal transfer from non-mycobacterial sources. The most frequent associated lineages (and potential gene transfer partners) in the Mycobacterium lineage-restricted gene trees are other members of suborder Corynebacterinae, but more-distant partners were identified as well. In two examined cases of potentially frequent and habitat-directed transfer (M. abscessus to Segniliparus and M. smegmatis to Streptomyces), observed sequence distances were small and consistent with a hypothesis of transfer, while in a third case (M. vanbaalenii to Streptomyces) distances were larger. The analyses described here indicate that whereas evidence of recombination in core regions within the genus is relatively sparse, the acquisition of genes from non-mycobacterial lineages is a significant feature of mycobacterial evolution.

Massive gene acquisitions in Mycobacterium indicus pranii provide a perspective on mycobacterial evolution

Understanding the evolutionary and genomic mechanisms responsible for turning the soil-derived saprophytic mycobacteria into lethal intracellular pathogens is a critical step towards the development of strategies for the control of mycobacterial diseases. In this context, Mycobacterium indicus pranii (MIP) is of specific interest because of its unique immunological and evolutionary significance. Evolutionarily, it is the progenitor of opportunistic pathogens belonging to M. avium complex and is endowed with features that place it between saprophytic and pathogenic species. Herein, we have sequenced the complete MIP genome to understand its unique life style, basis of immunomodulation and habitat diversification in mycobacteria. As a case of massive gene acquisitions, 50.5% of MIP open reading frames (ORFs) are laterally acquired. We show, for the first time for Mycobacterium, that MIP genome has mosaic architecture. These gene acquisitions have led to the enrichment of selected gene families critical to MIP physiology. Comparative genomic analysis indicates a higher antigenic potential of MIP imparting it a unique ability for immunomodulation. Besides, it also suggests an important role of genomic fluidity in habitat diversification within mycobacteria and provides a unique view of evolutionary divergence and putative bottlenecks that might have eventually led to intracellular survival and pathogenic attributes in mycobacteria.

Genome-Wide Sequence Variation among Mycobacterium avium Subspecies paratuberculosis Isolates: A Better Understanding of Johne?s Disease Transmission Dynamics

Mycobacterium avium subspecies paratuberculosis (M. ap), the causative agent of Johne’s disease, infects many farmed ruminants, wild-life animals, and recently isolated from humans. To better understand the molecular pathogenesis of these infections, we analyzed the whole-genome sequences of several M. ap and M. avium subspecies avium (M. avium) isolates to gain insights into genomic diversity associated with variable hosts and environments. Using Next-generation sequencing technology, all six M. ap isolates showed a high percentage of similarity (98%) to the reference genome sequence of M. ap K-10 isolated from cattle. However, two M. avium isolates (DT 78 and Env 77) showed significant sequence diversity (only 87 and 40% similarity, respectively) compared to the reference strain M. avium 104, a reflection of the wide environmental niches of this group of mycobacteria. Within the M. ap isolates, genomic rearrangements (insertions/deletions) were not detected, and only unique single nucleotide polymorphisms (SNPs) were observed among M. ap isolates. While more of the SNPs (~100) in M. ap genomes were non-synonymous, a total of ~6,000 SNPs were detected among M. avium genomes, most of them were synonymous suggesting a differential selective pressure between M. ap and M. avium isolates. In addition, SNPs-based phylo-genomics had a enough discriminatory power to differentiate between isolates from different hosts but yet suggesting a bovine source of infection to other animals examined in this study. Interestingly, the human isolate (M. ap 4B) was closely related to a M. ap isolate from a dairy facility, suggesting a common source of infection. Overall, the identified phylo-genomes further supported the idea of a common ancestor to both M. ap and M. avium isolates. Genome-wide analysis described here could provide a strong foundation for a population genetic structure that could be useful for the analysis of mycobacterial evolution and for the tracking of Johne’s disease transmission among animals.