Key experimental evidence of chromosomal DNA transfer among selected tuberculosis-causing mycobacteria

Interred mechanisms of resistance and host immune evasion revealed through network-connectivity analysis of M. tuberculosis complex graph pangenome

Mycobacterium tuberculosis complex successfully adapts to environmental pressures through mechanisms of rapid adaptation which remain poorly understood despite knowledge gained through decades of research. In this study, we used 110 reference-quality, complete de novo assembled, long-read sequenced clinical genomes to study patterns of structural adaptation through a graph-based pangenome analysis, elucidating rarely studied mechanisms that enable enhanced clinical phenotypes offering a novel perspective to the species' adaptation. Across isolates, we identified a pangenome of 4,325 genes (3,767 core and 558 accessory), revealing 290 novel genes, and a substantially more complete account of difficult-to-sequence esx/pe/pgrs/ppe genes. Seventy-four percent of core genes were deemed non-essential in vitro, 38% of which support the pathogen's survival in vivo, suggesting a need to broaden current perspectives on essentiality. Through information-theoretic analysis, we reveal the ppe genes that contribute most to the species' diversity-several with known consequences for antigenic variation and immune evasion. Construction of a graph pangenome revealed topological variations that implicate genes known to modulate host immunity (Rv0071-73, Rv2817c, cas2), defense against phages/viruses (cas2, csm6, and Rv2817c-2821c), and others associated with host tissue colonization. Here, the prominent trehalose transport pathway stands out for its involvement in caseous granuloma catabolism and the development of post-primary disease. We show paralogous duplications of genes implicated in bedaquiline (mmpL5 in all L1 isolates) and ethambutol (embC-A) resistance, with a paralogous duplication of its regulator (embR) in 96 isolates. We provide hypotheses for novel mechanisms of immune evasion and antibiotic resistance through gene dosing that can escape detection by molecular diagnostics.IMPORTANCEM. tuberculosis complex (MTBC) has killed over a billion people in the past 200 years alone and continues to kill nearly 1.5 million annually. The pathogen has a versatile ability to diversify under immune and drug pressure and survive, even becoming antibiotic persistent or resistant in the face of harsh chemotherapy. For proper diagnosis and design of an appropriate treatment regimen, a full understanding of this diversification and its clinical consequences is desperately needed. A mechanism of diversification that is rarely studied systematically is MTBC's ability to structurally change its genome. In this article, we have de novo assembled 110 clinical genomes (the largest de novo assembled set to date) and performed a pangenomic analysis. Our pangenome provides structural variation-based hypotheses for novel mechanisms of immune evasion and antibiotic resistance through gene dosing that can compromise molecular diagnostics and lead to further emergence of antibiotic resistance.

Distribution of lineages and type II toxin-antitoxin systems among rifampin-resistant Mycobacterium Tuberculosis Isolates

Type II toxin-antitoxin systems such as mazEF3, vapBC3, and relJK play a role in antibiotic resistance and tolerance. Among the different known TA systems, mazEF3, vapBC3, and relJK, which are type II systems, have specific roles in drug resistance. Therefore, the aim of this study was to investigate the mutations in these genes in sensitive and resistant isolates of Mycobacterium tuberculosis. Thirty-two rifampin-resistant and 121 rifampin-sensitive M. tuberculosis isolates were collected from various regions of Iran. Lineage typing was performed using the ASO-PCR method. Mutations in the rpoB gene were analyzed in all isolates by MAS-PCR. Furthermore, mutations in the mazEF3, relJK, and vapBC3 genes of the type II toxin system were assessed through PCR sequencing. These sequences were analyzed using COBALT and SnapGene 2017, and submitted to the GenBank database. Among the 153 M. tuberculosis samples, lineages 4, 3 and 2 were the most common. Lineage 2 had the highest rate of rifampin resistance. Mutations in rpoB531 were the most frequent in resistant isolates. Examination of the toxin-antitoxin system showed that rifampin-resistant isolates belonging to lineage 3 had mutations in either the toxin or antitoxin parts of all three TA systems. A mutation in nucleotide 195 (codon 65) of mazF3 leading to an amino acid change from threonine to isoleucine was detected in all rifampin-resistant isolates. M. tuberculosis isolates belonging to lineage 2 exhibited the highest rifampin resistance in our study. Identifying the mutation in mazF3 in all rifampin-resistant isolates can highlight the significance of this mutation in the development of drug resistance in M. tuberculosis. Expanding the sample size in future studies can help develop a new method for identifying resistant isolates.

A nucleoid-associated protein is involved in the emergence of antibiotic resistance by promoting the frequent exchange of the replicative DNA polymerase in Mycobacterium smegmatis

Antibiotic resistance in Mycobacterium tuberculosis exclusively originates from chromosomal mutations, either during normal DNA replication or under stress, when the expression of error-prone DNA polymerases increases to repair damaged DNA. To bypass DNA lesions and catalyze error-prone DNA synthesis, translesion polymerases must be able to access the DNA, temporarily replacing the high-fidelity replicative polymerase. The mechanisms that govern polymerase exchange are not well understood, especially in mycobacteria. Here, using a suite of quantitative fluorescence imaging techniques, we discover that in Mycobacterium smegmatis, as in other bacterial species, the replicative polymerase, DnaE1, exchanges at a timescale much faster than that of DNA replication. Interestingly, this fast exchange rate depends on an actinobacteria-specific nucleoid-associated protein (NAP), Lsr2. In cells missing lsr2, DnaE1 exchanges less frequently, and the chromosome is replicated more faithfully. Additionally, in conditions that damage DNA, cells lacking lsr2 load the complex needed to bypass DNA lesions less effectively and, consistently, replicate with higher fidelity but exhibit growth defects. Together, our results show that Lsr2 promotes dynamic flexibility of the mycobacterial replisome, which is critical for robust cell growth and lesion repair in conditions that damage DNA.

IMPORTANCE

Unlike many other pathogens, Mycobacterium tuberculosis has limited ability for horizontal gene transfer, a major mechanism for developing antibiotic resistance. Thus, the mechanisms that facilitate chromosomal mutagenesis are of particular importance in mycobacteria. Here, we show that Lsr2, a nucleoid-associated protein, has a novel role in DNA replication and mutagenesis in the model mycobacterium Mycobacterium smegmatis. We find that Lsr2 promotes the fast exchange rate of the replicative DNA polymerase, DnaE1, at the replication fork and is important for the effective loading of the DnaE2-ImuA'-ImuB translesion complex. Without lsr2, M. smegmatis replicates its chromosome more faithfully and acquires resistance to rifampin at a lower rate, but at the cost of impaired survival to DNA damaging agents. Together, our work establishes Lsr2 as a potential factor in the emergence of mycobacterial antibiotic resistance.

Genetic diversity, evolution and drug resistance of Mycobacterium tuberculosis lineage 2

Mycobacterium tuberculosis causes a chronic infectious disease called tuberculosis. Phylogenetic lineage 2 (L2) of M. tuberculosis, also known as the East Asian lineage, is associated with high virulence, increased transmissibility, and the spread of multidrug-resistant strains. This review article examines the genomic characteristics of the M. tuberculosis genome and M. tuberculosis lineage 2, such as the unique insertion sequence and spoligotype patterns, as well as MIRU-VNTR typing, and SNP-based barcoding. The review describes the geographical distribution of lineage 2 and its history of origin. In addition, the article discusses recent studies on drug resistance and compensatory mechanisms of M. tuberculosis lineage 2 and its impact on the pathogen's transmissibility and virulence. This review article discusses the importance of establishing a unified classification for lineage 2 to ensure consistency in terminology and criteria across different studies and settings.

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.

Large-scale Pan Genomic Analysis of Mycobacterium tuberculosis Reveals Key Insights Into Molecular Evolutionary Rate of Specific Processes and Functions

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), an infectious disease that is a major killer worldwide. Due to selection pressure caused by the use of antibacterial drugs, Mtb is characterised by mutational events that have given rise to multi drug resistant (MDR) and extensively drug resistant (XDR) phenotypes. The rate at which mutations occur is an important factor in the study of molecular evolution, and it helps understand gene evolution. Within the same species, different protein-coding genes evolve at different rates. To estimate the rates of molecular evolution of protein-coding genes, a commonly used parameter is the ratio dN/dS, where dN is the rate of non-synonymous substitutions and dS is the rate of synonymous substitutions. Here, we determined the estimated rates of molecular evolution of select biological processes and molecular functions across 264 strains of Mtb. We also investigated the molecular evolutionary rates of core genes of Mtb by computing the dN/dS values, and estimated the pan genome of the 264 strains of Mtb. Our results show that the cellular amino acid metabolic process and the kinase activity function evolve at a significantly higher rate, while the carbohydrate metabolic process evolves at a significantly lower rate for M. tuberculosis. These high rates of evolution correlate well with Mtb physiology and pathogenicity. We further propose that the core genome of M. tuberculosis likely experiences varying rates of molecular evolution which may drive an interplay between core genome and accessory genome during M. tuberculosis evolution.

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.

A nucleoid-associated protein is involved in the emergence of antibiotic resistance by promoting the frequent exchange of the replicative DNA polymerase in M. smegmatis

Antibiotic resistance in M. tuberculosis exclusively originates from chromosomal mutations, either during normal DNA replication or under stress, when the expression of error-prone DNA polymerases increases to repair damaged DNA. To bypass DNA lesions and catalyze error-prone DNA synthesis, translesion polymerases must be able to access the DNA, temporarily replacing the high-fidelity replicative polymerase. The mechanisms that govern polymerase exchange are not well understood, especially in mycobacteria. Here, using a suite of quantitative fluorescence imaging techniques, we discover that, as in other bacterial species, in M. smegmatis, the replicative polymerase, DnaE1, exchanges at a timescale much faster than that of DNA replication. Interestingly, this fast exchange rate depends on an actinobacteria-specific nucleoid-associated protein (NAP), Lsr2. In cells missing lsr2, DnaE1 exchanges less frequently, and the chromosome is replicated more faithfully. Additionally, in conditions that damage DNA, cells lacking lsr2 load the complex needed to bypass DNA lesions less effectively and, consistently, replicate with higher fidelity but exhibit growth defects. Together, our results show that Lsr2 promotes dynamic flexibility of the mycobacterial replisome, which is critical for robust cell growth and lesion repair in conditions that damage DNA.

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

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

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.

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.

Mechanisms of Linezolid Resistance in Mycobacteria

Mycobacteria form some of the most notorious and difficult-to-treat bacterial pathogens. As a group, they are intrinsically resistant to many commonly used antibiotics, such as tetracyclines and beta-lactams. In addition to intrinsic resistances, acquired multidrug resistance has also been observed and documented in Mycobacterium tuberculosis (MTB), Mycobacterium leprae and non-tuberculous mycobacteria (NTM). To combat multidrug resistant infections by these pathogens, innovative antimicrobials and treatment regimens are required. In this regard, linezolid, an oxazolidinone introduced for clinical use just two decades ago, was added to the therapeutic armamentarium for drug-resistant mycobacteria. It exhibits antibacterial activity by binding to the 50S ribosomal subunit and inhibiting protein synthesis. Unfortunately, linezolid resistance has now been documented in MTB and NTM, in many parts of the world. Most linezolid-resistant mycobacterial strains show mutations in the ribosome or related genes, such as in the rplC, rrl and tsnR genes. Non-ribosomal mechanisms appear to be rare. One such mechanism was associated with a mutation in fadD32, which encodes a protein that plays an important role in mycolic acid synthesis. Mycobacterial efflux proteins have also been implicated in linezolid resistance. This review summarises current knowledge of genetic determinants of linezolid resistance in mycobacteria, with the aim of contributing information that could facilitate the discovery of new therapeutic approaches to overcome, delay or avoid further developments of drug resistance among these important pathogens.

Mycobacterium smegmatis: The Vanguard of Mycobacterial Research

The genus Mycobacterium contains several slow-growing human pathogens, including Mycobacterium tuberculosis, Mycobacterium leprae, and Mycobacterium avium. Mycobacterium smegmatis is a nonpathogenic and fast growing species within this genus. In 1990, a mutant of M. smegmatis, designated mc2155, that could be transformed with episomal plasmids was isolated, elevating M. smegmatis to model status as the ideal surrogate for mycobacterial research. Classical bacterial models, such as Escherichia coli, were inadequate for mycobacteria research because they have low genetic conservation, different physiology, and lack the novel envelope structure that distinguishes the Mycobacterium genus. By contrast, M. smegmatis encodes thousands of conserved mycobacterial gene orthologs and has the same cell architecture and physiology. Dissection and characterization of conserved genes, structures, and processes in genetically tractable M. smegmatis mc2155 have since provided previously unattainable insights on these same features in its slow-growing relatives. Notably, tuberculosis (TB) drugs, including the first-line drugs isoniazid and ethambutol, are active against M. smegmatis, but not against E. coli, allowing the identification of their physiological targets. Furthermore, Bedaquiline, the first new TB drug in 40 years, was discovered through an M. smegmatis screen. M. smegmatis has become a model bacterium, not only for M. tuberculosis, but for all other Mycobacterium species and related genera. With a repertoire of bioinformatic and physical resources, including the recently established Mycobacterial Systems Resource, M. smegmatis will continue to accelerate mycobacterial research and advance the field of microbiology.

Unraveling the mechanisms of intrinsic drug resistance in Mycobacterium tuberculosis

Tuberculosis (TB) is among the most difficult infections to treat, requiring several months of multidrug therapy to produce a durable cure. The reasons necessitating long treatment times are complex and multifactorial. However, one major difficulty of treating TB is the resistance of the infecting bacterium, Mycobacterium tuberculosis (Mtb), to many distinct classes of antimicrobials. This review will focus on the major gaps in our understanding of intrinsic drug resistance in Mtb and how functional and chemical-genetics can help close those gaps. A better understanding of intrinsic drug resistance will help lay the foundation for strategies to disarm and circumvent these mechanisms to develop more potent antitubercular therapies.

Evolution of Mycobacterium tuberculosis drug resistance in the genomic era

Mycobacterium tuberculosis has acquired drug resistance to all drugs that have been used against it, including those only recently introduced into clinical practice. Compared to other bacteria, it has a well conserved genome due to its role as an obligate human pathogen that has adapted to a niche over five to ten thousand years. These features facilitate reconstruction and dating of M. tuberculosis phylogenies, giving key insights into how resistance has been acquired and spread globally. Resistance to each new drug has occurred within five to ten years of clinical use and has occurred even more rapidly with recently introduced drugs. In most cases, resistance-conferring mutations come with a fitness cost, but this can be overcome by compensatory mutations which restore fitness to that of wild-type bacteria. It is likely that M. tuberculosis acquires drug resistance while maintaining limited genomic variability due the generation of low frequency within-host variation, combined with ongoing purifying selection causing loss of variants without a clear fitness advantage. However, variants that do confer an advantage, such as drug resistance, can increase in prevalence amongst all bacteria within a host and become the dominant clone. These resistant strains can then be transmitted leading to primary drug resistant infection in a new host. As many countries move towards genomic methods for diagnosis of M. tuberculosis infection and drug resistance, it is important to be aware of the implications for the evolution of resistance. Currently, understanding of resistance-conferring mutations is incomplete, and some targeted genetic diagnostics create their own selective pressures. We discuss an example where a rifampicin resistance-conferring mutation which was not routinely covered by standard testing became dominant. Finally, resistance to new drugs such as bedaquiline and delamanid is caused by individually rare mutations occurring across a large mutational genomic target that have been detected over a short time, and do not provide statistical power for genotype-phenotype correlation – in contrast to longer-established drugs that form the backbone of drug-sensitive antituberculosis therapy. Therefore, we need a different approach to identify resistance-conferring mutations of new drugs before their resistance becomes widespread, abrogating their usefulness.

The rate and role of pseudogenes of the Mycobacterium tuberculosis complex

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

Evolution of Tuberculosis Pathogenesis

Mycobacterium tuberculosis is a globally distributed, lethal pathogen of humans. The virulence armamentarium of M. tuberculosis appears to have been developed on a scaffold of antiphagocytic defenses found among diverse, mostly free-living species of Mycobacterium. Pathoadaptation was further aided by the modularity, flexibility, and interactivity characterizing mycobacterial effectors and their regulators. During emergence of M. tuberculosis, novel genetic material was acquired, created, and integrated with existing tools. The major mutational mechanisms underlying these adaptations are discussed in this review, with examples. During its evolution, M. tuberculosis lost the ability and/or opportunity to engage in lateral gene transfer, but despite this it has retained the adaptability that characterizes mycobacteria. M. tuberculosis exemplifies the evolutionary genomic mechanisms underlying adoption of the pathogenic niche, and studies of its evolution have uncovered a rich array of discoveries about how new pathogens are made. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

A critical evaluation of Mycobacterium bovis pangenomics, with reference to its utility in outbreak investigation

The increased accessibility of next generation sequencing has allowed enough genomes from a given bacterial species to be sequenced to describe the distribution of genes in the pangenome, without limiting analyses to genes present in reference strains. Although some taxa have thousands of whole genome sequences available on public databases, most genomes were sequenced with short read technology, resulting in incomplete assemblies. Studying pangenomes could lead to important insights into adaptation, pathogenicity, or molecular epidemiology, however given the known information loss inherent in analyzing contig-level assemblies, these inferences may be biased or inaccurate. In this study we describe the pangenome of a clonally evolving pathogen, Mycobacterium bovis , and examine the utility of gene content variation in M. bovis outbreak investigation. We constructed the M. bovis pangenome using 1463 de novo assembled genomes. We tested the assumption of strict clonal evolution by studying evidence of recombination in core genes and analyzing the distribution of accessory genes among core monophyletic groups. To determine if gene content variation could be utilized in outbreak investigation, we carefully examined accessory genes detected in a well described M. bovis outbreak in Minnesota. We found significant errors in accessory gene classification. After accounting for these errors, we show that M. bovis has a much smaller accessory genome than previously described and provide evidence supporting ongoing clonal evolution and a closed pangenome, with little gene content variation generated over outbreaks. We also identified frameshift mutations in multiple genes, including a mutation in glpK , which has recently been associated with antibiotic tolerance in Mycobacterium tuberculosis . A pangenomic approach enables a more comprehensive analysis of genome dynamics than is possible with reference-based approaches; however, without critical evaluation of accessory gene content, inferences of transmission patterns employing these loci could be misguided.

The impact of frequently neglected model violations on bacterial recombination rate estimation: a case study in Mycobacterium canettii and Mycobacterium tuberculosis

Mycobacterium canettii is a causative agent of tuberculosis in humans, along with the members of the Mycobacterium tuberculosis complex. Frequently used as an outgroup to the M. tuberculosis complex in phylogenetic analyses, M. canettii is thought to offer the best proxy for the progenitor species that gave rise to the complex. Here, we leverage whole-genome sequencing data and biologically relevant population genomic models to compare the evolutionary dynamics driving variation in the recombining M. canettii with that in the nonrecombining M. tuberculosis complex, and discuss differences in observed genomic diversity in the light of expected levels of Hill-Robertson interference. In doing so, we highlight the methodological challenges of estimating recombination rates through traditional population genetic approaches using sequences called from populations of microorganisms and evaluate the likely mis-inference that arises owing to a neglect of common model violations including purifying selection, background selection, progeny skew, and population size change. In addition, we compare performance when full within-host polymorphism data are utilized, versus the more common approach of basing analyses on within-host consensus sequences.

Gene evolutionary trajectories in Mycobacterium tuberculosis reveal temporal signs of selection

Genetic differences between different Mycobacterium tuberculosis complex (MTBC) strains determine their ability to transmit within different host populations, their latency times, and their drug resistance profiles. Said differences usually emerge through de novo mutations and are maintained or discarded by the balance of evolutionary forces. Using a dataset of ∼5,000 strains representing global MTBC diversity, we determined the past and present selective forces that have shaped the current variability observed in the pathogen population. We identified regions that have evolved under changing types of selection since the time of the MTBC common ancestor. Our approach highlighted striking differences in the genome regions relevant for host–pathogen interaction and, in particular, suggested an adaptive role for the sensor protein of two-component systems. In addition, we applied our approach to successfully identify potential determinants of resistance to drugs administered as second-line tuberculosis treatments.

A Polymorphic Gene within the Mycobacterium smegmatis esx1 Locus Determines Mycobacterial Self-Identity and Conjugal Compatibility

Mycobacteria mediate horizontal gene transfer (HGT) by a process called distributive conjugal transfer (DCT) that is mechanistically distinct from oriT-mediated plasmid transfer. The transfer of multiple, independent donor chromosome segments generates transconjugants with genomes that are mosaic blends of their parents. Previously, we had characterized contact-dependent conjugation between two independent isolates of Mycobacterium smegmatis. Here, we expand our analyses to include five independent isolates of M. smegmatis and establish that DCT is both active and prevalent among natural isolates of M. smegmatis. Two of these five strains were recipients but exhibited distinct conjugal compatibilities with donor strains, suggesting an ability to distinguish between potential donor partners. We determined that a single gene, Msmeg0070, was responsible for conferring mating compatibility using a combination of comparative DNA sequence analysis, bacterial genome-wide association studies (GWAS), and targeted mutagenesis. Msmeg0070 maps within the esx1 secretion locus, and we establish that it confers mycobacterial self-identity with parallels to kin recognition. Similar to other kin model systems, orthologs of Msmeg0070 are highly polymorphic. The identification of a kin recognition system in M. smegmatis reinforces the concept that communication between cells is an important checkpoint prior to DCT commitment and implies that there are likely to be other, unanticipated forms of social behaviors in mycobacteria. IMPORTANCE Conjugation, unlike other forms of HGT, requires direct interaction between two viable bacteria, which must be capable of distinguishing between mating types to allow successful DNA transfer from donor to recipient. We show that the conjugal compatibility of Mycobacterium smegmatis isolates is determined by a single, polymorphic gene located within the conserved esx1 secretion locus. This gene confers self-identity; the expression of identical Msmeg0070 proteins in both donor-recipient partners prevents DNA transfer. The presence of this polymorphic locus in many environmental mycobacteria suggests that kin identification is important in promoting beneficial gene flow between nonkin mycobacteria. Cell-cell communication, mediated by kin recognition and ESX secretion, is a key checkpoint in mycobacterial conjugation and likely plays a more global role in mycobacterial biology.

Genome reorganization during emergence of host-associated Mycobacterium abscessus

Mycobacterium abscessus is a rapid growing, free-living species of bacterium that also causes lung infections in humans. Human infections are usually acquired from the environment; however, dominant circulating clones (DCCs) have emerged recently in both M. abscessus subsp. massiliense and subsp. abscessus that appear to be transmitted among humans and are now globally distributed. These recently emerged clones are potentially informative about the ecological and evolutionary mechanisms of pathogen emergence and host adaptation. The geographical distribution of DCCs has been reported, but the genomic processes underlying their transition from environmental bacterium to human pathogen are not well characterized. To address this knowledge gap, we delineated the structure of M. abscessus subspecies abscessus and massiliense using genomic data from 200 clinical isolates of M. abscessus from seven geographical regions. We identified differences in overall patterns of lateral gene transfer (LGT) and barriers to LGT between subspecies and between environmental and host-adapted bacteria. We further characterized genome reorganization that accompanied bacterial host adaptation, inferring selection pressures acting at both genic and intergenic loci. We found that both subspecies encode an expansive pangenome with many genes at rare frequencies. Recombination appears more frequent in M. abscessus subsp. massiliense than in subsp. abscessus, consistent with prior reports. We found evidence suggesting that phage are exchanged between subspecies, despite genetic barriers evident elsewhere throughout the genome. Patterns of LGT differed according to niche, with less LGT observed among host-adapted DCCs versus environmental bacteria. We also found evidence suggesting that DCCs are under distinct selection pressures at both genic and intergenic sites. Our results indicate that host adaptation of M. abscessus was accompanied by major changes in genome evolution, including shifts in the apparent frequency of LGT and impacts of selection. Differences were evident among the DCCs as well, which varied in the degree of gene content remodelling, suggesting they were placed differently along the evolutionary trajectory toward host adaptation. These results provide insight into the evolutionary forces that reshape bacterial genomes as they emerge into the pathogenic niche.

A new nomenclature for the livestock-associated Mycobacterium tuberculosis complex based on phylogenomics

Background:  The bacteria that compose the Mycobacterium tuberculosis complex (MTBC) cause tuberculosis (TB) in humans and in different animals, including livestock. Much progress has been made in understanding the population structure of the human-adapted members of the MTBC by combining phylogenetics with genomics. Accompanying the discovery of new genetic diversity, a body of operational nomenclature has evolved to assist comparative and molecular epidemiological studies of human TB. By contrast, for the livestock-associated MTBC members, Mycobacterium bovis, M. caprae and M. orygis, there has been a lack of comprehensive nomenclature to accommodate new genetic diversity uncovered by emerging phylogenomic studies. We propose to fill this gap by putting forward a new nomenclature covering the main phylogenetic groups within M. bovis, M. caprae and M. orygis. Methods:  We gathered a total of 8,736 whole-genome sequences (WGS) from public sources and 39 newly sequenced strains, and selected a subset of 829 WGS, representative of the worldwide diversity of M. bovis, M. caprae and M. orygis. We used phylogenetics and genetic diversity patterns inferred from WGS to define groups. Results:  We propose to divide M. bovis, M. caprae and M. orygis in three main phylogenetic lineages, which we named La1, La2 and La3, respectively. Within La1, we identified several monophyletic groups, which we propose to classify into eight sublineages (La1.1-La1.8). These sublineages differed in geographic distribution, with some being geographically restricted and others globally widespread, suggesting different expansion abilities. To ease molecular characterization of these MTBC groups by the community, we provide phylogenetically informed, single nucleotide polymorphisms that can be used as barcodes for genotyping. These markers were implemented in KvarQ and TB-Profiler, which are platform-independent, open-source tools. Conclusions:  Our results contribute to an improved classification of the genetic diversity within the livestock-associated MTBC, which will benefit future molecular epidemiological and evolutionary studies.

Type VII secretion systems: structure, functions and transport models

Type VII secretion systems (T7SSs) have a key role in the secretion of effector proteins in non-pathogenic mycobacteria and pathogenic mycobacteria such as Mycobacterium tuberculosis, the main causative agent of tuberculosis. Tuberculosis-causing mycobacteria, still accounting for 1.4 million deaths annually, rely on paralogous T7SSs to survive in the host and efficiently evade its immune response. Although it is still unknown how effector proteins of T7SSs cross the outer membrane of the diderm mycobacterial cell envelope, recent advances in the structural characterization of these secretion systems have revealed the intricate network of interactions of conserved components in the plasma membrane. This structural information, added to recent advances in the molecular biology and regulation of mycobacterial T7SSs as well as progress in our understanding of their secreted effector proteins, is shedding light on the inner working of the T7SS machinery. In this Review, we highlight the implications of these studies and the derived transport models, which provide new scenarios for targeting the deathly human pathogen M. tuberculosis.

The within-host evolution of antimicrobial resistance in Mycobacterium tuberculosis

Tuberculosis (TB) has been responsible for the greatest number of human deaths due to an infectious disease in general, and due to antimicrobial resistance (AMR) in particular. The etiological agents of human TB are a closely-related group of human-adapted bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Understanding how MTBC populations evolve within-host may allow for improved TB treatment and control strategies. In this review, we highlight recent works that have shed light on how AMR evolves in MTBC populations within individual patients. We discuss the role of heteroresistance in AMR evolution, and review the bacterial, patient and environmental factors that likely modulate the magnitude of heteroresistance within-host. We further highlight recent works on the dynamics of MTBC genetic diversity within-host, and discuss how spatial substructures in patients' lungs, spatiotemporal heterogeneity in antimicrobial concentrations and phenotypic drug tolerance likely modulates the dynamics of MTBC genetic diversity in patients during treatment. We note the general characteristics that are shared between how the MTBC and other bacterial pathogens evolve in humans, and highlight the characteristics unique to the MTBC.

Conserved and specialized functions of Type VII secretion systems in non-tuberculous mycobacteria

Non-tuberculous mycobacteria (NTM) are a large group of micro-organisms comprising more than 200 individual species. Most NTM are saprophytic organisms and are found mainly in terrestrial and aquatic environments. In recent years, NTM have been increasingly associated with infections in both immunocompetent and immunocompromised individuals, prompting significant efforts to understand the diverse pathogenic and signalling traits of these emerging pathogens. Since the discovery of Type VII secretion systems (T7SS), there have been significant developments regarding the role of these complex systems in mycobacteria. These specialised systems, also known as Early Antigenic Secretion (ESX) systems, are employed to secrete proteins across the inner membrane. They also play an essential role in virulence, nutrient uptake and conjugation. Our understanding of T7SS in mycobacteria has significantly benefited over the last few years, from the resolution of ESX-3 structure in Mycobacterium smegmatis, to ESX-5 structures in Mycobacterium xenopi and Mycobacterium tuberculosis. In addition, ESX-4, considered until recently as a non-functional system in both pathogenic and non-pathogenic mycobacteria, has been proposed to play an important role in the virulence of Mycobacterium abscessus; an increasingly recognized opportunistic NTM causing severe lung diseases. These major findings have led to important new insights into the functional mechanisms of these biological systems, their implication in virulence, nutrient acquisitions and cell wall shaping, and will be discussed in this review.

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

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

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

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

Measurable genomic changes in Mycobacterium avium subsp. hominissuis after long-term adaptation in Acanthamoeba lenticulata and reduced persistence in macrophages

Free-living amoebae are ubiquitous in aquatic environments and act as environmental reservoirs for nontuberculous mycobacteria. Mycobacterium avium subsp. hominissuis recovered from Acanthamoeba has been demonstrated to be more virulent in both human and murine models. Here, we investigate the persistence of M. avium subsp. hominissuis after short-term (2 weeks) and long-term (42 weeks) co-culture in Acanthamoeba lenticulata We hypothesize that A. lenticulata-adapted M. avium subsp. hominissuis demonstrate phenotypic and genomic changes facilitating intracellular persistence in naïve Acanthamoeba and human macrophages. M. avium subsp. hominissuis CFU in co-culture with A. lenticulata were recorded every 2 weeks up to 60 weeks. While A. lenticulata-associated M. avium subsp. hominissuis CFU did not significantly change across 60 weeks of co-culture, longer adaptation time in amoebae reduced colony size. Isolates recovered after 2 or 42 weeks of amoebae co-culture were referred as "early-adapted" and "late-adapted" M. avium subsp. hominissuis, respectively. Whole genome sequencing was performed on amoebae-adapted isolates with pan-genome comparisons to the original M. avium subsp. hominissuis isolate. Next, amoebae-adapted isolates were assessed for their persistence in A. lenticulata, A. castellanii, and human THP-1 macrophages. Multiplex cytokine/chemokine analyses were conducted on THP-1 culture supernatants.  Compared to the original isolate, counts of late-adapted M. avium subsp. hominissuis were reduced in Acanthamoeba and contrary to expectations, lower counts were also observed in THP-1 macrophages with concomitant decrease in TNFa, IL-6, and MIP-1b suggesting that host adaptation may influence the inflammatory properties of M. avium IMPORTANCE Short-term interaction between Acanthamoeba and M. avium has been demonstrated to increase infectivity in human and murine models of infection, establishing the paradigm that amoebae "train" M. avium in the environment by selecting for phenotypes capable of enduring in human cells. We investigate this phenomenon further by determining the consequence of long-term amoebae adaptation on M. avium subsp. hominissuis persistence in host cells. We monitored genomic changes across long-term Acanthamoeba co-culture and report significant changes to the M. avium subsp. hominissuis genome in response to amoebae-adaptation and reduced colony size. Furthermore, we examined isolates co-cultured with A. lenticulata for 2 or 42 weeks and provide biological evidence that long-term co-culture in amoebae reduces M. avium persistence in human macrophages.

Pathogenic Determinants of the Mycobacterium kansasii Complex: An Unsuspected Role for Distributive Conjugal Transfer

The Mycobacterium kansasii species comprises six subtypes that were recently classified into six closely related species; Mycobacterium kansasii (formerly M. kansasii subtype 1), Mycobacterium persicum (subtype 2), Mycobacterium pseudokansasii (subtype 3), Mycobacterium ostraviense (subtype 4), Mycobacterium innocens (subtype 5) and Mycobacterium attenuatum (subtype 6). Together with Mycobacterium gastri, they form the M. kansasii complex. M. kansasii is the most frequent and most pathogenic species of the complex. M. persicum is classically associated with diseases in immunosuppressed patients, and the other species are mostly colonizers, and are only very rarely reported in ill patients. Comparative genomics was used to assess the genetic determinants leading to the pathogenicity of members of the M. kansasii complex. The genomes of 51 isolates collected from patients with and without disease were sequenced and compared with 24 publicly available genomes. The pathogenicity of each isolate was determined based on the clinical records or public metadata. A comparative genomic analysis showed that all M. persicum, M. ostraviense, M innocens and M. gastri isolates lacked the ESX-1-associated EspACD locus that is thought to play a crucial role in the pathogenicity of M. tuberculosis and other non-tuberculous mycobacteria. Furthermore, M. kansasii was the only species exhibiting a 25-Kb-large genomic island encoding for 17 type-VII secretion system-associated proteins. Finally, a genome-wide association analysis revealed that two consecutive genes encoding a hemerythrin-like protein and a nitroreductase-like protein were significantly associated with pathogenicity. These two genes may be involved in the resistance to reactive oxygen and nitrogen species, a required mechanism for the intracellular survival of bacteria. Three non-pathogenic M. kansasii lacked these genes likely due to two distinct distributive conjugal transfers (DCTs) between M. attenuatum and M. kansasii, and one DCT between M. persicum and M. kansasii. To our knowledge, this is the first study linking DCT to reduced pathogenicity.

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

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

Direct PCR on Tissue Samples To Detect Mycobacterium tuberculosis Complex: an Alternative to the Bacteriological Culture

Bovine tuberculosis (bTB) is an ongoing issue in several countries within the European Union. Microbiological culture is the official confirmation technique for the presence of Mycobacterium tuberculosis complex (MTBC) members in bovine tissues, but several methodological issues, such as moderate sensitivity and long incubation times, require the development of more sensitive and rapid techniques. This study evaluates the analytical and diagnostic performance, comparative to culture, of a real-time PCR targeting the MTBC-specific IS6110 transposon using a panel of bovine tissue samples sourced from the Spanish bTB eradication campaign. Robustness and repeatability were evaluated in an interlaboratory trial between European Union National Reference Laboratories. The limit of detection with 95% confidence was established at 65 fg/reaction of purified genomic equivalents. Diagnostic sensitivity (Se) and specificity (Sp) were, respectively, 96.45% and 93.66%, and the overall agreement (κ) was 0.88. Cross-reactivity was detected against two mycobacterial isolates identified as Mycobacterium marinum and "Mycobacterium avium subsp. hominissuis," and whole-genome sequencing (WGS) analysis of the latter isolate revealed an IS6110-like sequence with 83% identity. An identical IS-like element was found in other Mycobacterium avium complex species in the NCBI nucleotide and WGS databases. Despite this finding, this methodology is considered a valuable alternative to culture, and the strategy of use should be defined depending on the control or eradication programs.

A fine-scale map of genome-wide recombination in divergent Escherichia coli population

Recombination is one of the most important molecular mechanisms of prokaryotic genome evolution, but its exact roles are still in debate. Here we try to infer genome-wide recombination within a species, utilizing a dataset of 149 complete genomes of Escherichia coli from diverse animal hosts and geographic origins, including 45 in-house sequenced with the single-molecular real-time platform. Two major clades identified based on physiological, clinical and ecological characteristics form distinct genetic lineages based on scarcity of interclade gene exchanges. By defining gene-based syntenies for genomic segments within and between the two clades, we build a fine-scale recombination map for this representative global E. coli population. The map suggests extensive within-clade recombination that often breaks physical linkages among individual genes but seldom interrupts the structure of genome organizational frameworks as well as primary metabolic portfolios supported by the framework integrity, possibly due to strong natural selection for both physiological compatibility and ecological fitness. In contrast, the between-clade recombination declines drastically when phylogenetic distance increases to the extent where a 10-fold reduction can be observed, establishing a firm genetic barrier between clades. Our empirical data suggest a critical role for such recombination events in the early stage of speciation where recombination rate is associated with phylogenetic distance in addition to sequence and gene variations. The extensive intraclade recombination binds sister strains into a quasisexual group and optimizes genes or alleles to streamline physiological activities, whereas the sharply declined interclade recombination split the population into clades adaptive to divergent ecological niches.

Deep learning for HGT insertion sites recognition

Background

Horizontal Gene Transfer (HGT) refers to the sharing of genetic materials between distant species that are not in a parent-offspring relationship. The HGT insertion sites are important to understand the HGT mechanisms. Recent studies in main agents of HGT, such as transposon and plasmid, demonstrate that insertion sites usually hold specific sequence features. This motivates us to find a method to infer HGT insertion sites according to sequence features.

Results

In this paper, we propose a deep residual network, DeepHGT, to recognize HGT insertion sites. To train DeepHGT, we extracted about 1.55 million sequence segments as training instances from 262 metagenomic samples, where the ratio between positive instances and negative instances is about 1:1. These segments are randomly partitioned into three subsets: 80% of them as the training set, 10% as the validation set, and the remaining 10% as the test set. The training loss of DeepHGT is 0.4163 and the validation loss is 0.423. On the test set, DeepHGT has achieved the area under curve (AUC) value of 0.8782. Furthermore, in order to further evaluate the generalization of DeepHGT, we constructed an independent test set containing 689,312 sequence segments from another 147 gut metagenomic samples. DeepHGT has achieved the AUC value of 0.8428, which approaches the previous test AUC value. As a comparison, the gradient boosting classifier model implemented in PyFeat achieve an AUC value of 0.694 and 0.686 on the above two test sets, respectively. Furthermore, DeepHGT could learn discriminant sequence features; for example, DeepHGT has learned a sequence pattern of palindromic subsequences as a significantly (P-value=0.0182) local feature. Hence, DeepHGT is a reliable model to recognize the HGT insertion site.

Conclusion

DeepHGT is the first deep learning model that can accurately recognize HGT insertion sites on genomes according to the sequence pattern.

Heterologous Production of 1-Tuberculosinyladenosine in Mycobacterium kansasii Models Pathoevolution towards the Transcellular Lifestyle of Mycobacterium tuberculosis

Mycobacterium kansasii is an environmental nontuberculous mycobacterium that causes opportunistic tuberculosis-like disease. It is one of the most closely related species to the Mycobacterium tuberculosis complex. Using M. kansasii as a proxy for the M. kansasii-M. tuberculosis common ancestor, we asked whether introducing the M. tuberculosis-specific gene pair Rv3377c-Rv3378c into M. kansasii affects the course of experimental infection. Expression of these genes resulted in the production of an adenosine-linked lipid species, known as 1-tuberculosinyladenosine (1-TbAd), but did not alter growth in vitro under standard conditions. Production of 1-TbAd enhanced growth of M. kansasii under acidic conditions through a bacterial cell-intrinsic mechanism independent of controlling pH in the bulk extracellular and intracellular spaces. Production of 1-TbAd led to greater burden of M. kansasii in the lungs of C57BL/6 mice during the first 24 h after infection, and ex vivo infections of alveolar macrophages recapitulated this phenotype within the same time frame. However, in long-term infections, production of 1-TbAd resulted in impaired bacterial survival in both C57BL/6 mice and Ccr2^-/- mice. We have demonstrated that M. kansasii is a valid surrogate of M. tuberculosis to study virulence factors acquired by the latter organism, yet shown the challenge inherent to studying the complex evolution of mycobacterial pathogenicity with isolated gene complementation.IMPORTANCE This work sheds light on the role of the lipid 1-tuberculosinyladenosine in the evolution of an environmental ancestor to M. tuberculosis On a larger scale, it reinforces the importance of horizontal gene transfer in bacterial evolution and examines novel models and methods to provide a better understanding of the subtle effects of individual M. tuberculosis-specific virulence factors in infection settings that are relevant to the pathogen.

Insertion and deletion evolution reflects antibiotics selection pressure in a Mycobacterium tuberculosis outbreak

In genome evolution, genetic variants are the source of diversity, which natural selection acts upon. Treatment of human tuberculosis (TB) induces a strong selection pressure for the emergence of antibiotic resistance-conferring variants in the infecting Mycobacterium tuberculosis (MTB) strains. MTB evolution in response to treatment has been intensively studied and mainly attributed to point substitutions. However, the frequency and contribution of insertions and deletions (indels) to MTB genome evolution remains poorly understood. Here, we analyzed a multi-drug resistant MTB outbreak for the presence of high-quality indels and substitutions. We find that indels are significantly enriched in genes conferring antibiotic resistance. Furthermore, we show that indels are inherited during the outbreak and follow a molecular clock with an evolutionary rate of 5.37e-9 indels/site/year, which is 23 times lower than the substitution rate. Inherited indels may co-occur with substitutions in genes along related biological pathways; examples are iron storage and resistance to second-line antibiotics. This suggests that epistatic interactions between indels and substitutions affect antibiotic resistance and compensatory evolution in MTB.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Delving Into the Functional Meaning of Phenotypic Variation in Mycobacterial Persistence: Who Benefits the Most From Programmed Death of Individual Cells?

The lengthy tuberculosis therapy is emblematic of how hard drug-persistent infections are to eradicate. Phenotypic variation within clonal bacterial communities contributes to drug evasion and has major implications for the treatment of drug-persistent infections. We reported that single mycobacterial cells exhibit differential drug susceptibility, contingent on their inherent phenotypic variation in DNA damage response. Individual cells experiencing severe DNA damage massively induce the SOS response and exhibit signs of programmed cell death (PCD), such as unbalanced growth, chromosomal fragmentation, autolysis, and release of the intracellular content. Toxin-antitoxin systems are known to contribute to PCD in model microorganisms by targeting essential cellular processes, and they might function similarly in mycobacteria. We have found that the toxin MazF and a Clp protease, possibly responsible for degrading the MazF cognate antitoxin MazE, are induced during harsh conditions in a model organism for tuberculosis, and that cells that are about to lyse from drug exposure display a buildup of toxin. Deeper analysis of PCD in mycobacteria may reveal whether this process belongs to a broader strategy for the community's survival. Finally, disrupting the balance between survival and PCD may prove useful to tackle drug evasion in mycobacterial persistent subpopulations.

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

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

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

Shared Pathogenomic Patterns Characterize a New Phylotype, Revealing Transition toward Host-Adaptation Long before Speciation of Mycobacterium tuberculosis

Tuberculosis remains one of the deadliest infectious diseases of humanity. To better understand the evolutionary history of host-adaptation of tubercle bacilli (MTB), we sought for mycobacterial species that were more closely related to MTB than the previously used comparator species Mycobacterium marinum and Mycobacterium kansasii. Our phylogenomic approach revealed some recently sequenced opportunistic mycobacterial pathogens, Mycobacterium decipiens, Mycobacterium lacus, Mycobacterium riyadhense, and Mycobacterium shinjukuense, to constitute a common clade with MTB, hereafter called MTB-associated phylotype (MTBAP), from which MTB have emerged. Multivariate and clustering analyses of genomic functional content revealed that the MTBAP lineage forms a clearly distinct cluster of species that share common genomic characteristics, such as loss of core genes, shift in dN/dS ratios, and massive expansion of toxin–antitoxin systems. Consistently, analysis of predicted horizontal gene transfer regions suggests that putative functions acquired by MTBAP members were markedly associated with changes in microbial ecology, for example adaption to intracellular stress resistance. Our study thus considerably deepens our view on MTB evolutionary history, unveiling a decisive shift that promoted conversion to host-adaptation among ancestral founders of the MTBAP lineage long before Mycobacterium tuberculosis has adapted to the human host.

The Genetic Background Modulates the Evolution of Fluoroquinolone-Resistance in Mycobacterium tuberculosis

Fluoroquinolones (FQ) form the backbone in experimental treatment regimens against drug-susceptible tuberculosis. However, little is known on whether the genetic variation present in natural populations of Mycobacterium tuberculosis (Mtb) affects the evolution of FQ-resistance (FQ-R). To investigate this question, we used nine genetically distinct drug-susceptible clinical isolates of Mtb and measured their frequency of resistance to the FQ ofloxacin (OFX) in vitro. We found that the Mtb genetic background led to differences in the frequency of OFX-resistance (OFX-R) that spanned two orders of magnitude and substantially modulated the observed mutational profiles for OFX-R. Further, in vitro assays showed that the genetic background also influenced the minimum inhibitory concentration and the fitness effect conferred by a given OFX-R mutation. To test the clinical relevance of our in vitro work, we surveyed the mutational profile for FQ-R in publicly available genomic sequences from clinical Mtb isolates, and found substantial Mtb lineage-dependent variability. Comparison of the clinical and the in vitro mutational profiles for FQ-R showed that 51% and 39% of the variability in the clinical frequency of FQ-R gyrA mutation events in Lineage 2 and Lineage 4 strains, respectively, can be attributed to how Mtb evolves FQ-R in vitro. As the Mtb genetic background strongly influenced the evolution of FQ-R in vitro, we conclude that the genetic background of Mtb also impacts the evolution of FQ-R in the clinic.

Quantifying and Characterizing Distributive Conjugal Transfer in Mycobacterium smegmatis

Horizontal gene transfer (HGT) in prokaryotes disseminates genetic information throughout a population and can facilitate adaptation and evolution of the species. Mycobacteria utilize an atypical method of conjugation called distributive conjugal transfer (DCT), which results in mosaic genomes and the potential for accelerated evolution beyond that enabled by the more classical oriT-mediated conjugation. The following is a description of the basic DCT protocol, some possible variations of the assay, and examples of downstream applications to better understand mycobacterial functions.

Three-reaction high-resolution melting assay for rapid differentiation of Mycobacterium tuberculosis complex members

The possibility of introducing a reliable assay for a quick identification and differentiation of the main species of Mycobacterium tuberculosis complex (MTBC) supports the improvement of efficient tuberculosis combating strategies worldwide. Commercially available assays are often based on cultured samples; however, due to the long cultivation time of mycobacteria, results are delayed. Developed PCR approaches have been published previously, though, when testing intricate veterinary samples, the complex composition of multiplex qPCRs frequently leads to assay failure. In order to overcome those limits, a paradigm of a three-reaction high-resolution melting (HRM) assay for the simultaneous identification and differentiation of the main members of MTBC was established. The assay is based on single nucleotide polymorphisms within gyrB and gyrA, which have been used as target for the establishment of two highly specific HRM assays (HRM assays 1 and 2) discriminating M. tuberculosis/ Mycobacterium canetti, Mycobacterium bovis/M. bovis BCG, Mycobacterium caprae/rare M. caprae/M. bovis ecotypes, Mycobacterium africanum/Mycobacterium orygis/ Mycobacterium pinnipedii/Clade A1, Mycobacterium microti, and a rare subtype of M. canettii followed by a third HRM assay (HRM assay 3) allowing a further differentiation of M. bovis, M. bovis BCG, and a rare subtype of M. caprae/M. bovis, which is considered to be a novel ecotype. High-resolution melting assay 1 is described in a previously published report. High-resolution melting assay 2 showed 100% correlation of all 39 examined isolates with the results of a commercial identification kit. 96% of the clinical samples tested demonstrated concordant results. High-resolution melting assay 3 showed an accordance of 100% with the results of the commercially available identification kit of all 22 samples analyzed. The proposed strategy of the three-reaction HRM assay can be used for an accurate differentiation of up to seven groups of MTBC and potentially to identify a rare subtype of M. canettii either on isolates or on clinical samples.

ESX/Type VII Secretion Systems-An Important Way Out for Mycobacterial Proteins

The causative agent of human tuberculosis, Mycobacterium tuberculosis, has a complex lipid-rich diderm envelope, which acts as a major barrier protecting the bacterium against the hostile environment inside the host cells. For the transfer of diverse molecules across this complex cell envelope, M. tuberculosis has a series of general and specialized protein secretion systems, characterized by the SecA general secretion pathway, the twin-arginine translocation pathway, and five specific ESX type VII secretion systems. In this review, we focus on the latter systems, known as ESX-1 to ESX-5, which were first discovered almost 20 years ago during the in silico analysis of the genome sequence of M. tuberculosis H37Rv. Since then, these systems have been the subject of highly dynamic research due to their involvement in several key biological processes and host-pathogen interactions of the tubercle bacilli.

Protein Export into and across the Atypical Diderm Cell Envelope of Mycobacteria

Mycobacteria, including the infamous pathogen Mycobacterium tuberculosis, are high-GC Gram-positive bacteria with a distinctive cell envelope. Although there is a typical inner membrane, the mycobacterial cell envelope is unusual in having its peptidoglycan layer connected to a polymer of arabinogalactan, which in turn is covalently attached to long-chain mycolic acids that help form a highly impermeable mycobacterial outer membrane. This complex double-membrane, or diderm, cell envelope imparts mycobacteria with unique requirements for protein export into and across the cell envelope for secretion into the extracellular environment. In this article, we review the four protein export pathways known to exist in mycobacteria: two conserved systems that exist in all types of bacteria (the Sec and Tat pathways) and two specialized systems that exist in mycobacteria, corynebacteria, and a subset of low-GC Gram-positive bacteria (the SecA2 and type VII secretion pathways). We describe the progress made over the past 15 years in understanding each of these mycobacterial export pathways, and we highlight the need for research to understand the specific steps of protein export across the mycobacterial outer membrane.

Comparing the Metabolic Capabilities of Bacteria in the Mycobacterium tuberculosis Complex

Pathogenic mycobacteria are known for their ability to maintain persistent infections in various mammals. The canonical pathogen in this genus is Mycobacterium tuberculosis and this bacterium is particularly successful at surviving and replicating within macrophages. Here, we will highlight the metabolic processes that M. tuberculosis employs during infection in macrophages and compare these findings with what is understood for other pathogens in the M. tuberculosis complex.

Genomic determinants of speciation and spread of the Mycobacterium tuberculosis complex

Models on how bacterial lineages differentiate increase our understanding of early bacterial speciation events and the genetic loci involved. Here, we analyze the population genomics events leading to the emergence of the tuberculosis pathogen. The emergence is characterized by a combination of recombination events involving core pathogenesis functions and purifying selection on early diverging loci. We identify the phoR gene, the sensor kinase of a two-component system involved in virulence, as a key functional player subject to pervasive positive selection after the divergence of the Mycobacterium tuberculosis complex from its ancestor. Previous evidence showed that phoR mutations played a central role in the adaptation of the pathogen to different host species. Now, we show that phoR mutations have been under selection during the early spread of human tuberculosis, during later expansions, and in ongoing transmission events. Our results show that linking pathogen evolution across evolutionary and epidemiological time scales points to past and present virulence determinants.

Revised Interpretation of the Hain Lifescience GenoType MTBC To Differentiate Mycobacterium canettii and Members of the Mycobacterium tuberculosis Complex

Using 894 phylogenetically diverse genomes of the Mycobacterium tuberculosis complex (MTBC), we simulated in silico the ability of the Hain Lifescience GenoType MTBC assay to differentiate the causative agents of tuberculosis. Here, we propose a revised interpretation of this assay to reflect its strengths (e.g., it can distinguish some strains of Mycobacterium canettii and variants of Mycobacterium bovis that are not intrinsically resistant to pyrazinamide) and limitations (e.g., Mycobacterium orygis cannot be differentiated from Mycobacterium africanum).