TASmania: A bacterial Toxin-Antitoxin Systems database

Novel type II toxin-antitoxin systems with VapD-like proteins

Type II toxin-antitoxin (TA) systems are widespread in prokaryotes. They consist of neighboring genes encoding two small proteins: a toxin that inhibits a critical cellular process and an antitoxin that binds to and neutralizes the toxin. The VapD nuclease and the VapX antitoxin comprise a type II TA system that contributes to the virulence of the human pathogen Haemophilus influenzae. We analyzed the diversity and evolution of VapD-like proteins. By examining loci adjacent to genes coding for VapD-like proteins, we identified two novel families of antitoxins, which we named VapY and VapW. VapD toxins cognate to novel antitoxins induce the SOS response when overproduced, suggesting they target cellular processes related to genomic DNA integrity, maintenance, or replication. Though VapY has no sequence similarity to VapX, they share the same SH3 fold characterized by the five anti-parallel β sheets that form a barrel. VapW is a homolog of VapD without conserved catalytic residues required for nuclease activity. The crystal structure of the VapD-VapW complex reveals that VapW lacks the dimerization interface essential for the catalytic activity of VapD but retains the second interaction interface that enables VapD hexamerization. This allows VapW to bind VapD in the same manner that VapD dimers bind to each other in hexamers. Thus, though the VapD catalytic cleft remains accessible in the VapD-VapW complex, VapW may disrupt VapD oligomerization. To our knowledge, VapWD provides a unique example of TA systems evolution when a toxin loses its activity and becomes an antitoxin to itself.

IMPORTANCE

Genes encoding virulence-associated protein D (VapD) homologs are found in many pathogens such as Helicobacter pylori, Haemophilus influenzae, and Xylella fastidiosa. There are many indications that VapD proteins contribute to virulence, even though the exact mechanism is not known. VapD proteins are either encoded by stand-alone genes or form toxin-antitoxin pairs with VapX. We performed a comprehensive census of vapD-like genes and found two new antitoxins, VapW and VapY. The VapW antitoxins are catalytically inactivated variants of VapD, revealing a new evolutionary mechanism for the appearance of toxin-antitoxin pairs.

Carbapenem-Resistant Pseudomonas aeruginosa's Resistome: Pan-Genomic Plasticity, the Impact of Transposable Elements and Jumping Genes

Pseudomonas aeruginosa, a Gram-negative, motile bacterium, may cause significant infections in both community and hospital settings, leading to substantial morbidity and mortality. This opportunistic pathogen can thrive in various environments, making it a public health concern worldwide. P. aeruginosa's genomic pool is highly dynamic and diverse, with a pan-genome size ranging from 5.5 to 7.76 Mbp. This versatility arises from its ability to acquire genes through horizontal gene transfer (HGT) via different genetic elements (GEs), such as mobile genetic elements (MGEs). These MGEs, collectively known as the mobilome, facilitate the spread of genes encoding resistance to antimicrobials (ARGs), resistance to heavy metals (HMRGs), virulence (VGs), and metabolic functions (MGs). Of particular concern are the acquired carbapenemase genes (ACGs) and other β-lactamase genes, such as classes A, B [metallo-β-lactamases (MBLs)], and D carbapenemases, which can lead to increased antimicrobial resistance. This review emphasizes the importance of the mobilome in understanding antimicrobial resistance in P. aeruginosa.

Understanding the physiological role and cross-interaction network of VapBC35 toxin-antitoxin system from Mycobacterium tuberculosis

The VapBC toxin-antitoxin (TA) system, composed of VapC toxin and VapB antitoxin, has gained attention due to its relative abundance in members of the M. tuberculosis complex. Here, we have functionally characterised VapBC35 TA system from M. tuberculosis. We show that ectopic expression of VapC35 inhibits M. smegmatis growth in a bacteriostatic manner. Also, an increase in the VapB35 antitoxin to VapC35 toxin ratio results in a stronger binding affinity of the complex with the promoter-operator DNA. We show that VapBC35 is necessary for M. tuberculosis adaptation in oxidative stress conditions but is dispensable for M. tuberculosis growth in guinea pigs. Further, using a combination of co-expression studies and biophysical methods, we report that VapC35 also interacts with non-cognate antitoxin VapB3. Taken together, the present study advances our understanding of cross-interaction networks among VapBC TA systems from M. tuberculosis.

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.

Insight into the environmental cues modulating the expression of bacterial toxin-antitoxin systems

Bacteria require sophisticated sensing mechanisms to adjust their metabolism in response to stressful conditions and survive in hostile environments. Among them, toxin-antitoxin (TA) systems play a crucial role in bacterial adaptation to environmental challenges. TA systems are considered as stress-responsive elements, consisting of both toxin and antitoxin genes, typically organized in operons or encoded on complementary DNA strands. A decrease in the antitoxin-toxin ratio, often triggered by specific stress conditions, leads to toxin excess, disrupting essential cellular processes and inhibiting bacterial growth. These systems are categorized into eight types based on the nature of the antitoxin (RNA or protein) and the mode of action of toxin inhibition. While the well-established biological roles of TA systems include phage inhibition and the maintenance of genetic elements, the environmental cues regulating their expression remain insufficiently documented. In this review, we highlight the diversity and complexity of the environmental cues influencing TA systems expression. A comprehensive understanding of how these genetic modules are regulated could provide deeper insights into their functions and support the development of innovative antimicrobial strategies.

Genomic profiling and molecular dynamics analysis of parDEPa toxin-antitoxin homologs targeting DNA gyrase in Pseudomonas aeruginosa: insights from computational investigations

In the realm of hospital-acquired and chronic infections, Pseudomonas aeruginosa stands out, demonstrating significant associations with increased morbidity, mortality, and antibiotic resistance. Antibiotic-resistant strains are believed to contribute to thousands of deaths each year. Chronic and latent infections are associated with the bacterial toxin-antitoxin (TA) system, although the mechanisms involved are poorly understood. This study focuses on a novel type II TA system, parDEPa, identified in the genome of P. aeruginosa ATCC 27853. We explored its structural features, functional relationships, and genetic configurations. Our research identified parDEPa homologs in P. aeruginosa, clarified their interactions, and highlighted connections to essential cellular metabolic processes. Notably, homologs of the ParDPa antitoxin were found to be more conserved than the ParEPa toxin. Structural models of the ParEPa toxin and ParDPa antitoxin confirmed their integrity. Through docking and molecular dynamics simulations, we showed that the ParEPa toxin binds to DNA gyrase, inhibiting replication. The stability of the ParDPa-ParEPa complex is primarily driven by hydrophobic interactions (-1763.2 kcal/mol), while the ParEPa-GyrAPa interaction is sustained by strong electrostatic forces (-1294.9 kcal/mol). The RMSD scores indicated greater stability for the ParDPa-ParEPa complex (1.11 Å) than the ParEPa-GyrAPa complex (1.16 Å). RMSF analysis identified key residues involved in the ParDPa-ParEPa complex (Leu59, Gly60, Arg115, Asn116, Arg117) and the ParEPa-GyrAPa complex (Pro48, Gln49, Ser55, Asp94, Gln95). These findings significantly enhance our understanding of the structural and metabolic roles of the chromosomally encoded parDEPa TA module in P. aeruginosa.

Type I toxin-antitoxin systems in bacteria: from regulation to biological functions

Toxin-antitoxin systems are ubiquitous in the prokaryotic world and widely distributed among chromosomes and mobile genetic elements. Several different toxin-antitoxin system types exist, but what they all have in common is that toxin activity is prevented by the cognate antitoxin. In type I toxin-antitoxin systems, toxin production is controlled by an RNA antitoxin and by structural features inherent to the toxin messenger RNA. Most type I toxins are small membrane proteins that display a variety of cellular effects. While originally discovered as modules that stabilize plasmids, chromosomal type I toxin-antitoxin systems may also stabilize prophages, or serve important functions upon certain stress conditions and contribute to population-wide survival strategies. Here, we will describe the intricate RNA-based regulation of type I toxin-antitoxin systems and discuss their potential biological functions.

Draft genome sequence data of Serratia marcescens strain harboring blaNDM-7 from Dhaka, Bangladesh

Here, the draft genome sequence of a multi-drug resistant (MDR) Serratia marcescens strain BMD28, isolated from a clinical source from Dhaka, Bangladesh, has been reported. The sequence raw read files were generated using Illumina sequencing technology utilizing genomic DNA from the pure culture of this strain. The strain has a genome size of around 5.4 million base pairs, a GC content of 59.70 %, and 5,141 coding sequences. We conducted genomic studies using several bioinformatics tools focusing on resistance genes, virulence factors, toxin-antitoxin systems, and pangenome analysis. Strain BMD28 harbored the blaNDM-7 gene in an IncX3 plasmid. A phylogenomic study with S. marcescens strains isolated worldwide revealed that our strain is in the same clade as other strains reported in Bangladesh. The data can be used primarily to understand the genomic content, epidemiology, and evolution of S. marcescens in Bangladesh. The genome sequence data of BMD28 has been deposited in the NCBI database under BioSample accession number SAMN41260295.

Draft genome sequencing data of multidrug-resistant Staphylococcus haemolyticus from Bangladeshi hospitals

Here, we report the draft genome sequence data of two multidrug-resistant (MDR) Staphylococcus haemolyticus strains, SAC2 and SAC7, isolated from clinical samples from Dhaka, Bangladesh. The sequence raw read files were generated using Ion Torrent Sequencing Technology using the genomic DNA from the pure culture of the strains. These two Bangladeshi S. haemolyticus strains had an average genome size of 2.49 million base pairs with a GC content of 32.6 % and an average of 1783 coding sequences. We conducted genomic studies using bioinformatics tools focusing on resistance genes, virulence factors, and toxin-antitoxin systems. A phylogenomic study with S. haemolyticus strains isolated worldwide revealed that these two Bangladeshi strains are in different nodes but clustered together. The data can be used as a starting point for understanding the genomic content, epidemiology, and evolution of S. haemolyticus in Bangladesh. The genome sequence data of SAC2 and SAC7 strains have been deposited in the NCBI database under BioSample accession numbers SAMN35731443 and SAMN35731649, respectively.

Nucleotidyltransferase toxin MenT extends aminoacyl acceptor ends of serine tRNAs to control Mycobacterium tuberculosis growth

Toxins of toxin-antitoxin systems use diverse mechanisms to inhibit bacterial growth. In this study, we characterize the translation inhibitor toxin MenT3 of Mycobacterium tuberculosis, the bacterium responsible for tuberculosis in humans. We show that MenT3 is a robust cytidine specific tRNA nucleotidyltransferase in vitro, capable of modifying the aminoacyl acceptor ends of most tRNA but with a marked preference for tRNASer, to which long stretches of cytidines are added. Furthermore, transcriptomic-wide analysis of MenT3 targets in M. tuberculosis identifies tRNASer as the sole target of MenT3 and reveals significant detoxification attempts by the essential CCA-adding enzyme PcnA in response to MenT3. Finally, under physiological conditions, only in the presence the native menAT3 operon, an active pool of endogenous MenT3 targeting tRNASer in M. tuberculosis is detected, likely reflecting the importance of MenT3 during infection.

Chromosomal Type II Toxin-Antitoxin Systems May Enhance Bacterial Fitness of a Hybrid Pathogenic Escherichia coli Strain Under Stress Conditions

The functions of bacterial plasmid-encoded toxin-antitoxin (TA) systems are unambiguous in the sense of controlling cells that fail to inherit a plasmid copy. However, its role in chromosomal copies is contradictory, including stress-response-promoting fitness and antibiotic treatment survival. A hybrid pathogenic Escherichia coli strain may have the ability to colonize distinct host niches, facing contrasting stress environments. Herein, we determined the influence of multiple environmental stress factors on the bacterial growth dynamic and expression profile of previously described TA systems present in the chromosome of a hybrid atypical enteropathogenic and extraintestinal E. coli strain. Genomic analysis revealed 26 TA loci and the presence of five type II TA systems in the chromosome. Among the tested stress conditions, osmotic and acid stress significantly altered the growth dynamics of the hybrid strain, enhancing the necessary time to reach the stationary phase. Using qPCR analyses, 80% of the studied TA systems were differentially expressed in at least one of the tested conditions, either in the log or in the stationary phase. These data indicate that type II TA systems may contribute to the physiology of pathogenic hybrid strains, enabling their adaptation to different milieus.

Deciphering the role of VapBC13 and VapBC26 toxin antitoxin systems in the pathophysiology of Mycobacterium tuberculosis

The expansion of VapBC TA systems in M. tuberculosis has been linked with its fitness and survival upon exposure to stress conditions. Here, we have functionally characterized VapBC13 and VapBC26 TA modules of M. tuberculosis. We report that overexpression of VapC13 and VapC26 toxins in M. tuberculosis results in growth inhibition and transcriptional reprogramming. We have also identified various regulatory proteins as hub nodes in the top response network of VapC13 and VapC26 overexpression strains. Further, analysis of RNA protection ratios revealed potential tRNA targets for VapC13 and VapC26. Using in vitro ribonuclease assays, we demonstrate that VapC13 and VapC26 degrade serT and leuW tRNA, respectively. However, no significant changes in rRNA cleavage profiles were observed upon overexpression of VapC13 and VapC26 in M. tuberculosis. In order to delineate the role of these TA systems in M. tuberculosis physiology, various mutant strains were constructed. We show that in comparison to the parental strain, ΔvapBC13 and ΔvapBC26 strains were mildly susceptible to oxidative stress. Surprisingly, the growth patterns of parental and mutant strains were comparable in aerosol-infected guinea pigs. These observations imply that significant functional redundancy exists for some TA systems from M. tuberculosis.

Toxin-Antitoxin Systems Reflect Community Interactions Through Horizontal Gene Transfer

Bacterial evolution through horizontal gene transfer (HGT) reflects their community interactions. In this way, HGT networks do well at mapping community interactions, but offer little toward controlling them-an important step in the translation of synthetic strains into natural contexts. Toxin-antitoxin (TA) systems serve as ubiquitous and diverse agents of selection; however, their utility is limited by their erratic distribution in hosts. Here we examine the heterogeneous distribution of TAs as a consequence of their mobility. By systematically mapping TA systems across a 10,000 plasmid network, we find HGT communities have unique and predictable TA signatures. We propose these TA signatures arise from plasmid competition and have further potential to signal the degree to which plasmids, hosts, and phage interact. To emphasize these relationships, we construct an HGT network based solely on TA similarity, framing specific selection markers in the broader context of bacterial communities. This work both clarifies the evolution of TA systems and unlocks a common framework for manipulating community interactions through TA compatibility.

Deciphering the genomic character of the multidrug-resistant Staphylococcus aureus from Dhaka, Bangladesh

Staphylococcus aureus is one of the leading agents of nosocomial and community-acquired infections. In this study, we explored the genomic characterization of eight methicillin-resistant clinical isolates of S. aureus from Dhaka, Bangladesh. Notably, all strains were resistant to penicillin, cephalosporins, and monobactams, with partial susceptibility to meropenem and complete susceptibility to amikacin, vancomycin, and tigecycline antibiotics. The strains were found to have an average genome size of 2.73 Mbp and an average of 32.64% GC content. Multi-locus sequence typing analysis characterized the most predominant sequence type as ST361, which belongs to the clonal complex CC361. All isolates harbored the mecA gene, often linked to SCCmec_type IV variants. Multidrug resistance was attributed to efflux pumps NorA, NorC, SdrM, and LmrS alongside genes encoding beta-lactamase BlaZ and factors like ErmC and MepA. Additionally, virulence factors including adsA, sdrC, cap8D, harA, esaA, essC, isdB, geh, and lip were commonly identified. Furthermore, genes associated with heme uptake and clumping were present, highlighting their roles in S. aureus colonization and pathogenesis. Nine secondary metabolite biosynthetic gene clusters were found, of which six were common in all the strains. Numerous toxin-antitoxin systems were predicted, with ParE and ParB-like nuclease domains found to be the most prevalent toxin and antitoxin, respectively. Pan-genome analysis revealed 2007 core genes and 229 unique genes in the studied strains. Finally, the phylogenomic analysis showed that most Bangladeshi strains were grouped into two unique clades. This study provides a genomic and comparative insight into the multidrug resistance and pathogenicity of S. aureus strains, which will play a crucial role in the future antibiotic stewardship of Bangladesh.

Decoding the TAome and computational insights into parDE toxin-antitoxin systems in Pseudomonas aeruginosa

Toxin-antitoxin (TA) modules are widely found in the genomes of pathogenic bacteria. They regulate vital cellular functions like transcription, translation, and DNA replication, and are therefore essential to the survival of bacteria under stress. With a focus on the type II parDE modules, this study thoroughly examines TAome in Pseudomonas aeruginosa, a bacterium well-known for its adaptability and antibiotic resistance. We explored the TAome in three P. aeruginosa strains: ATCC 27,853, PAO1, and PA14, and found 15 type II TAs in ATCC 27,853, 12 in PAO1, and 13 in PA14, with significant variation in the associated mobile genetic elements. Five different parDE homologs were found by further TAome analysis in ATCC 27,853, and their relationships were confirmed by sequence alignments and precise genomic positions. After comparing these ParDE modules' sequences to those of other pathogenic bacteria, it was discovered that they were conserved throughout many taxa, especially Proteobacteria. Nucleic acids were predicted as potential ligands for ParD antitoxins, whereas ParE toxins interacted with a wide range of small molecules, indicating a diverse functional repertoire. The interaction interfaces between ParDE TAs were clarified by protein-protein interaction networks and docking studies, which also highlighted important residues involved in binding. This thorough examination improves our understanding of the diversity, evolutionary dynamics, and functional significance of TA systems in P. aeruginosa, providing insights into their roles in bacterial physiology and pathogenicity.

Mycobacterium tuberculosis strain with deletions in menT3 and menT4 is attenuated and confers protection in mice and guinea pigs

The genome of Mycobacterium tuberculosis encodes for a large repertoire of toxin-antitoxin systems. In the present study, MenT3 and MenT4 toxins belonging to MenAT subfamily of TA systems have been functionally characterized. We demonstrate that ectopic expression of these toxins inhibits bacterial growth and this is rescued upon co-expression of their cognate antitoxins. Here, we show that simultaneous deletion of menT3 and menT4 results in enhanced susceptibility of M. tuberculosis upon exposure to oxidative stress and attenuated growth in guinea pigs and mice. We observed reduced expression of transcripts encoding for proteins that are essential or required for intracellular growth in mid-log phase cultures of ΔmenT4ΔT3 compared to parental strain. Further, the transcript levels of proteins involved in efficient bacterial clearance were increased in lung tissues of ΔmenT4ΔT3 infected mice relative to parental strain infected mice. We show that immunization of mice and guinea pigs with ΔmenT4ΔT3 confers significant protection against M. tuberculosis infection. Remarkably, immunization of mice with ΔmenT4ΔT3 results in increased antigen-specific TH1 bias and activated memory T cell response. We conclude that MenT3 and MenT4 are important for M. tuberculosis pathogenicity and strains lacking menT3 and menT4 have the potential to be explored further as vaccine candidates.

HigA2 (Rv2021c) Is a Transcriptional Regulator with Multiple Regulatory Targets in Mycobacterium tuberculosis

Toxin-antitoxin (TA) systems are the major mechanism for persister formation in Mycobacterium tuberculosis (Mtb). Previous studies found that HigBA2 (Rv2022c-Rv2021c), a predicted type II TA system of Mtb, could be activated for transcription in response to multiple stresses such as anti-tuberculosis drugs, nutrient starvation, endure hypoxia, acidic pH, etc. In this study, we determined the binding site of HigA2 (Rv2021c), which is located in the coding region of the upstream gene higB2 (Rv2022c), and the conserved recognition motif of HigA2 was characterized via oligonucleotide mutation. Eight binding sites of HigA2 were further found in the Mtb genome according to the conserved motif. RT-PCR showed that HigA2 can regulate the transcription level of all eight of these genes and three adjacent downstream genes. DNA pull-down experiments showed that twelve functional regulators sense external regulatory signals and may regulate the transcription of the HigBA2 system. Of these, Rv0903c, Rv0744c, Rv0474, Rv3124, Rv2603c, and Rv3583c may be involved in the regulation of external stress signals. In general, we identified the downstream target genes and possible upstream regulatory genes of HigA2, which paved the way for the illustration of the persistence establishment mechanism in Mtb.

Metabolic engineering of Geobacillus thermoglucosidasius for polymer-grade lactic acid production at high temperature

The production and application of biodegradable polylactic acid are still severely hindered by the cost of its polymer-grade lactic acid monomers. High-temperature biomanufacturing has emerged as an increasingly attractive approach to enable low-cost and high-efficiency bulk chemical production. In this study, thermophilic Geobacillus thermoglucosidasius was reprogrammed to obtain optically pure l-lactic acid- and d-lactic acid-producing strains, G. thermoglucosidasius GTD17 and GTD7, by using rational metabolic engineering strategies including pathway construction, by-product elimination, and production enhancing. Moreover, semi-rational adaptive evolution was carried out to further improve their lactic acid synthesis performance. The final strains GTD17-55 and GTD7-144 produce 151.1 g/L of l-lactic acid and 153.1 g/L of d-lactic acid at 60 °C, respectively. In consideration of the high temperature, productive performance of these strains is superior compared to the state-of-the-art industrial strains. This study lays the foundation for the low-cost and efficient production of biodegradable plastic polylactic acid.

Unraveling the evolutionary dynamics of toxin-antitoxin systems in diverse genetic lineages of Escherichia coli including the high-risk clonal complexes

IMPORTANCE

Large-scale genomic studies of E. coli provide an invaluable opportunity to understand how genomic fine-tuning contributes to the transition of bacterial lifestyle from being commensals to mutualists or pathogens. Within this context, through machine learning-based studies, it appears that TA systems play an important role in the classification of high-risk clonal lineages and could be attributed to their epidemiological success. Due to these profound indications and assumptions, we attempted to provide unique insights into the ordered world of TA systems at the population level by investigating the diversity and evolutionary patterns of TA genes across 19 different STs of E. coli. Further in-depth analysis of ST-specific TA structures and associated genetic coordinates holds the potential to elucidate the functional implications of TA systems in bacterial cell survival and persistence, by and large.

Cassette recombination dynamics within chromosomal integrons are regulated by toxin-antitoxin systems

Integrons are adaptive bacterial devices that rearrange promoter-less gene cassettes into variable ordered arrays under stress conditions, thereby sampling combinatorial phenotypic diversity. Chromosomal integrons often carry hundreds of silent gene cassettes, with integrase-mediated recombination leading to rampant DNA excision and integration, posing a potential threat to genome integrity. How this activity is regulated and controlled, particularly through selective pressures, to maintain such large cassette arrays is unknown. Here, we show a key role of promoter-containing toxin-antitoxin (TA) cassettes as systems that kill the cell when the overall cassette excision rate is too high. These results highlight the importance of TA cassettes regulating the cassette recombination dynamics and provide insight into the evolution and success of integrons in bacterial genomes.

Editorial for Special Issue "Bacterial Toxin-Antitoxin Systems"

Toxin antitoxin systems (TAS) are widely distributed in bacterial chromosomes as well as on mobile genetic elements [...].

TADB 3.0: an updated database of bacterial toxin-antitoxin loci and associated mobile genetic elements

TADB 3.0 (https://bioinfo-mml.sjtu.edu.cn/TADB3/) is an updated database that provides comprehensive information on bacterial types I to VIII toxin-antitoxin (TA) loci. Compared with the previous version, three major improvements are introduced: First, with the aid of text mining and manual curation, it records the details of 536 TA loci with experimental support, including 102, 403, 8, 14, 1, 1, 3 and 4 TA loci of types I to VIII, respectively; Second, by leveraging the upgraded TA prediction tool TAfinder 2.0 with a stringent strategy, TADB 3.0 collects 211 697 putative types I to VIII TA loci predicted in 34 789 completely sequenced prokaryotic genomes, providing researchers with a large-scale dataset for further follow-up analysis and characterization; Third, based on their genomic locations, relationships of 69 019 TA loci and 60 898 mobile genetic elements (MGEs) are visualized by interactive networks accessible through the user-friendly web page. With the recent updates, TADB 3.0 may provide improved in silico support for comprehending the biological roles of TA pairs in prokaryotes and their functional associations with MGEs.

The Benefits of Toxicity: M. smegmatis VapBC TA Module Is Induced by Tetracycline Exposure and Promotes Survival

Toxin-antitoxin (TA) systems are widely present in bacterial genomes. Mycolicibacterium smegmatis, a common model organism for studying Mycobacterium tuberculosis physiology, has eight TA loci, including mazEF and vapBC. This study aims to investigate the physiological significance of these TA systems. Proteomic profiling was conducted on a culture overexpressing the VapC toxin, and the involvement of VapC in M. smegmatis stress responses to heat shock and antibiotic treatment was examined. While deciphering the underlying mechanisms of the altered stress resistance, we assessed the antibiotic susceptibility of vapBC, mazEF, and double vapBC-mazEF deletion mutants. Additionally, the mRNA levels of vapC and mazF were measured following tetracycline supplementation. The results reveal changes in the abundance of metabolic enzymes and stress response proteins associated with VapC overexpression. This activation of the general stress response leads to reduced thermosensitivity in M. smegmatis, but does not affect susceptibility to ciprofloxacin and isoniazid. Under tetracycline treatment, both vapC and mazF expression levels are increased, and the fate of the cell depends on the interaction between the corresponding TA systems.

A toxin-antidote selfish element increases fitness of its host

Selfish genetic elements can promote their transmission at the expense of individual survival, creating conflict between the element and the rest of the genome. Recently, a large number of toxin-antidote (TA) post-segregation distorters have been identified in non-obligate outcrossing nematodes. Their origin and the evolutionary forces that keep them at intermediate population frequencies are poorly understood. Here, we study a TA element in Caenorhabditis elegans called zeel-1;peel-1. Two major haplotypes of this locus, with and without the selfish element, segregate in C. elegans. We evaluate the fitness consequences of the zeel-1;peel-1 element outside of its role in gene drive in non-outcrossing animals and demonstrate that loss of the toxin peel-1 decreased fitness of hermaphrodites and resulted in reductions in fecundity and body size. These findings suggest a biological role for peel-1 beyond toxin lethality. This work demonstrates that a TA element can provide a fitness benefit to its hosts either during their initial evolution or by being co-opted by the animals following their selfish spread. These findings guide our understanding on how TA elements can remain in a population where gene drive is minimized, helping resolve the mystery of prevalent TA elements in selfing animals.

Bacterial toxin-antitoxin systems: Novel insights on toxin activation across populations and experimental shortcomings

The alarming rise in hard-to-treat bacterial infections is of great concern to human health. Thus, the identification of molecular mechanisms that enable the survival and growth of pathogens is of utmost urgency for the development of more efficient antimicrobial therapies. In challenging environments, such as presence of antibiotics, or during host infection, metabolic adjustments are essential for microorganism survival and competitiveness. Toxin-antitoxin systems (TASs) consisting of a toxin with metabolic modulating activity and a cognate antitoxin that antagonizes that toxin are important elements in the arsenal of bacterial stress defense. However, the exact physiological function of TA systems is highly debatable and with the exception of stabilization of mobile genetic elements and phage inhibition, other proposed biological functions lack a broad consensus. This review aims at gaining new insights into the physiological effects of TASs in bacteria and exploring the experimental shortcomings that lead to discrepant results in TAS research. Distinct control mechanisms ensure that only subsets of cells within isogenic cultures transiently develop moderate levels of toxin activity. As a result, TASs cause phenotypic growth heterogeneity rather than cell stasis in the entire population. It is this feature that allows bacteria to thrive in diverse environments through the creation of subpopulations with different metabolic rates and stress tolerance programs.

A Comprehensive Bioinformatics Resource Guide for Genome-Based Antimicrobial Resistance Studies

The use of high-throughput sequencing technologies and bioinformatic tools has greatly transformed microbial genome research. With the help of sophisticated computational tools, it has become easier to perform whole genome assembly, identify and compare different species based on their genomes, and predict the presence of genes responsible for proteins, antimicrobial resistance, and toxins. These bioinformatics resources are likely to continuously improve in quality, become more user-friendly to analyze the multiple genomic data, efficient in generating information and translating it into meaningful knowledge, and enhance our understanding of the genetic mechanism of AMR. In this manuscript, we provide an essential guide for selecting the popular resources for microbial research, such as genome assembly and annotation, antibiotic resistance gene profiling, identification of virulence factors, and drug interaction studies. In addition, we discuss the best practices in computer-oriented microbial genome research, emerging trends in microbial genomic data analysis, integration of multi-omics data, the appropriate use of machine-learning algorithms, and open-source bioinformatics resources for genome data analytics.

Focused Overview of Mycobacterium tuberculosis VapBC Toxin-Antitoxin Systems Regarding Their Structural and Functional Aspects: Including Insights on Biomimetic Peptides

Tuberculosis, caused by Mycobacterium tuberculosis, is a lethal infectious disease of significant public health concern. The rise of multidrug-resistant and drug-tolerant strains has necessitated novel approaches to combat the disease. Toxin-antitoxin (TA) systems, key players in bacterial adaptive responses, are prevalent in prokaryotic genomes and have been linked to tuberculosis. The genome of M. tuberculosis strains harbors an unusually high number of TA systems, prompting questions about their biological roles. The VapBC family, a representative type II TA system, is characterized by the VapC toxin, featuring a PilT N-terminal domain with nuclease activity. Its counterpart, VapB, functions as an antitoxin, inhibiting VapC's activity. Additionally, we explore peptide mimics designed to replicate protein helical structures in this review. Investigating these synthetic peptides offers fresh insights into molecular interactions, potentially leading to therapeutic applications. These synthetic peptides show promise as versatile tools for modulating cellular processes and protein-protein interactions. We examine the rational design strategies employed to mimic helical motifs, their biophysical properties, and potential applications in drug development and bioengineering. This review aims to provide an in-depth understanding of TA systems by introducing known complex structures, with a focus on both structural aspects and functional and molecular details associated with each system.

Type II bacterial toxin-antitoxins: hypotheses, facts, and the newfound plethora of the PezAT system

Toxin-antitoxin (TA) systems are entities found in the prokaryotic genomes, with eight reported types. Type II, the best characterized, is comprised of two genes organized as an operon. Whereas toxins impair growth, the cognate antitoxin neutralizes its activity. TAs appeared to be involved in plasmid maintenance, persistence, virulence, and defence against bacteriophages. Most Type II toxins target the bacterial translational machinery. They seem to be antecessors of Higher Eukaryotes and Prokaryotes Nucleotide-binding (HEPN) RNases, minimal nucleotidyltransferase domains, or CRISPR-Cas systems. A total of four TAs encoded by Streptococcus pneumoniae, RelBE, YefMYoeB, Phd-Doc, and HicAB, belong to HEPN-RNases. The fifth is represented by PezAT/Epsilon-Zeta. PezT/Zeta toxins phosphorylate the peptidoglycan precursors, thereby blocking cell wall synthesis. We explore the body of knowledge (facts) and hypotheses procured for Type II TAs and analyse the data accumulated on the PezAT family. Bioinformatics analyses showed that homologues of PezT/Zeta toxin are abundantly distributed among 14 bacterial phyla mostly in Proteobacteria (48%), Firmicutes (27%), and Actinobacteria (18%), showing the widespread distribution of this TA. The pezAT locus was found to be mainly chromosomally encoded whereas its homologue, the tripartite omega-epsilon-zeta locus, was found mostly on plasmids. We found several orphan pezT/zeta toxins, unaccompanied by a cognate antitoxin.

Genomic assortment and interactive insights of the chromosomal encoded control of cell death (ccd) toxin-antitoxin (TA) module in Xenorhabdus nematophila

In the present circumstances, toxin-antitoxin (TA) modules have a great consideration due to their elusive role in bacterial physiology. TA modules consist of a toxic part and a counteracting antitoxin part and these are abundant genetic loci harbored on bacterial plasmids and chromosomes. The control of cell death (ccd) TA locus was the first identified TA module and its unitary function (such as plasmid maintenance) has been described, however, the function of its chromosomal counterparts is still ambiguous. Here, we are exploring the genomic assortment, structural and functional association of chromosomally encoded ccdAB TA homolog (ccdABXn1) in the genome of an entomopathogenic bacterium Xenorhabdus nematophila. This bacterium is a symbiotic model with the nematode Steinernema carpocapsae that infects and kills the host insect. By genomic assortment analysis, our observations suggested that CcdA antitoxin homologs are not more closely related than CcdB toxin homologs. Further results suggest that the ccdABXn1 TA homolog has sulphonamide (such as 4C6, for CcdA homolog) and peptide (such as gyrase, for CcdB homolog) ligand partners with a typical TA interaction network that may affect essential cellular metabolism of the X. nematophila. Collectively, our results improve the knowledge and conception of the metabolic interactive role of ccdAB TA homologs in X. nematophila physiology.Communicated by Ramaswamy H. Sarma.

MenT nucleotidyltransferase toxins extend tRNA acceptor stems and can be inhibited by asymmetrical antitoxin binding

Mycobacterium tuberculosis, the bacterium responsible for human tuberculosis, has a genome encoding a remarkably high number of toxin-antitoxin systems of largely unknown function. We have recently shown that the M. tuberculosis genome encodes four of a widespread, MenAT family of nucleotidyltransferase toxin-antitoxin systems. In this study we characterize MenAT1, using tRNA sequencing to demonstrate MenT1 tRNA modification activity. MenT1 activity is blocked by MenA1, a short protein antitoxin unrelated to the MenA3 kinase. X-ray crystallographic analysis shows blockage of the conserved MenT fold by asymmetric binding of MenA1 across two MenT1 protomers, forming a heterotrimeric toxin-antitoxin complex. Finally, we also demonstrate tRNA modification by toxin MenT4, indicating conserved activity across the MenT family. Our study highlights variation in tRNA target preferences by MenT toxins, selective use of nucleotide substrates, and diverse modes of MenA antitoxin activity.

The structural basis of hyperpromiscuity in a core combinatorial network of type II toxin-antitoxin and related phage defense systems

Toxin-antitoxin (TA) systems are a large group of small genetic modules found in prokaryotes and their mobile genetic elements. Type II TAs are encoded as bicistronic (two-gene) operons that encode two proteins: a toxin and a neutralizing antitoxin. Using our tool NetFlax (standing for Network-FlaGs for toxins and antitoxins), we have performed a large-scale bioinformatic analysis of proteinaceous TAs, revealing interconnected clusters constituting a core network of TA-like gene pairs. To understand the structural basis of toxin neutralization by antitoxins, we have predicted the structures of 3,419 complexes with AlphaFold2. Together with mutagenesis and functional assays, our structural predictions provide insights into the neutralizing mechanism of the hyperpromiscuous Panacea antitoxin domain. In antitoxins composed of standalone Panacea, the domain mediates direct toxin neutralization, while in multidomain antitoxins the neutralization is mediated by other domains, such as PAD1, Phd-C, and ZFD. We hypothesize that Panacea acts as a sensor that regulates TA activation. We have experimentally validated 16 NetFlax TA systems and used domain annotations and metabolic labeling assays to predict their potential mechanisms of toxicity (such as membrane disruption, and inhibition of cell division or protein synthesis) as well as biological functions (such as antiphage defense). We have validated the antiphage activity of a RosmerTA system encoded by Gordonia phage Kita, and used fluorescence microscopy to confirm its predicted membrane-depolarizing activity. The interactive version of the NetFlax TA network that includes structural predictions can be accessed at http://netflax.webflags.se/.

Bacterial Toxin-Antitoxin Systems' Cross-Interactions-Implications for Practical Use in Medicine and Biotechnology

Toxin-antitoxin (TA) systems are widely present in bacterial genomes. They consist of stable toxins and unstable antitoxins that are classified into distinct groups based on their structure and biological activity. TA systems are mostly related to mobile genetic elements and can be easily acquired through horizontal gene transfer. The ubiquity of different homologous and non-homologous TA systems within a single bacterial genome raises questions about their potential cross-interactions. Unspecific cross-talk between toxins and antitoxins of non-cognate modules may unbalance the ratio of the interacting partners and cause an increase in the free toxin level, which can be deleterious to the cell. Moreover, TA systems can be involved in broadly understood molecular networks as transcriptional regulators of other genes' expression or modulators of cellular mRNA stability. In nature, multiple copies of highly similar or identical TA systems are rather infrequent and probably represent a transition stage during evolution to complete insulation or decay of one of them. Nevertheless, several types of cross-interactions have been described in the literature to date. This implies a question of the possibility and consequences of the TA system cross-interactions, especially in the context of the practical application of the TA-based biotechnological and medical strategies, in which such TAs will be used outside their natural context, will be artificially introduced and induced in the new hosts. Thus, in this review, we discuss the prospective challenges of system cross-talks in the safety and effectiveness of TA system usage.

Enterococcal Linear Plasmids Adapt to Enterococcus faecium and Spread within Multidrug-Resistant Clades

Antimicrobial resistance (AMR) of bacterial pathogens, including enterococci, is a global concern, and plasmids are crucial for spreading and maintaining AMR genes. Plasmids with linear topology were identified recently in clinical multidrug-resistant enterococci. The enterococcal linear-form plasmids, such as pELF1, confer resistance to clinically important antimicrobials, including vancomycin; however, little information exists about their epidemiological and physiological effects. In this study, we identified several lineages of enterococcal linear plasmids that are structurally conserved and occur globally. pELF1-like linear plasmids show plasticity in acquiring and maintaining AMR genes, often via transposition with the mobile genetic element IS1216E. This linear plasmid family has several characteristics enabling long-term persistence in the bacterial population, including high horizontal self-transmissibility, low-level transcription of plasmid-carried genes, and a moderate effect on the Enterococcus faecium genome alleviating fitness cost and promoting vertical inheritance. Combining all of these factors, the linear plasmid is an important factor in the spread and maintenance of AMR genes among enterococci.

The DarT/DarG Toxin-Antitoxin ADP-Ribosylation System as a Novel Target for a Rational Design of Innovative Antimicrobial Strategies

The chemical modification of cellular macromolecules by the transfer of ADP-ribose unit(s), known as ADP-ribosylation, is an ancient homeostatic and stress response control system. Highly conserved across the evolution, ADP-ribosyltransferases and ADP-ribosylhydrolases control ADP-ribosylation signalling and cellular responses. In addition to proteins, both prokaryotic and eukaryotic transferases can covalently link ADP-ribosylation to different conformations of nucleic acids, thus highlighting the evolutionary conservation of archaic stress response mechanisms. Here, we report several structural and functional aspects of DNA ADP-ribosylation modification controlled by the prototype DarT and DarG pair, which show ADP-ribosyltransferase and hydrolase activity, respectively. DarT/DarG is a toxin-antitoxin system conserved in many bacterial pathogens, for example in Mycobacterium tuberculosis, which regulates two clinically important processes for human health, namely, growth control and the anti-phage response. The chemical modulation of the DarT/DarG system by selective inhibitors may thus represent an exciting strategy to tackle resistance to current antimicrobial therapies.

Toxin-Antitoxin Systems Alter Adaptation of Mycobacterium smegmatis to Environmental Stress

Toxin-antitoxin (TA) systems are ubiquitous genetic elements in prokaryotes, but their biological importance is poorly understood. Mycobacterium smegmatis contains eight putative TA systems. Previously, seven TAs have been studied, with five of them being verified as functional. Here, we show that Ms0251-0252 is a novel TA system in that expression of the toxin Ms0251 leads to growth inhibition that can be rescued by the antitoxin Ms0252. To investigate the functional roles of TA systems in M. smegmatis, we deleted the eight putative TA loci and assayed the mutants for resistance to various stresses. Deletion of all eight TA loci resulted in decreased survival under starvation conditions and altered fitness when exposed to environmental stresses. Furthermore, we showed that deletion of the eight TA loci decreased resistance to phage infection in Sauton medium compared with the results using 7H10 medium, suggesting that TA systems might have different contributions depending on the nutrient environment. Furthermore, we found that MazEF specifically played a dominant role in resistance to phage infection. Finally, transcriptome analysis revealed that MazEF overexpression led to differential expression of multiple genes, including those related to iron acquisition. Altogether, we demonstrate that TA systems coordinately function to allow M. smegmatis to adapt to changing environmental conditions. IMPORTANCE Toxin-antitoxin (TA) systems are mechanisms for rapid adaptation of bacteria to environmental changes. Mycobacterium smegmatis, a model bacterium for studying Mycobacterium tuberculosis, encodes eight putative TA systems. Here, we constructed an M. smegmatis mutant with deletions of all eight TA-encoding genes and evaluated the resistance of these mutants to environmental stresses. Our results showed that different TA systems have overlapping and, in some cases, opposing functions in adaptation to various stresses. We suggest that complementary TA modules may function together to regulate the bacterial stress response, enabling adaptation to changing environments. Together, this study provides key insights into the roles of TA systems in resistance to various environmental stresses, drug tolerance, and defense against phage infection.

Disruption of MenT2 toxin impairs the growth of Mycobacterium tuberculosis in guinea pigs

Toxin-antitoxin (TA) systems are abundantly present in the genomes of various bacterial pathogens. TA systems have been implicated in either plasmid maintenance or protection against phage infection, stress adaptation or disease pathogenesis. The genome of Mycobacterium tuberculosis encodes for more than 90 TA systems and 4 of these belong to the type IV subfamily (MenAT family). The toxins and antitoxins belonging to type IV TA systems share sequence homology with the AbiEii family of nucleotidyl transferases and the AbiEi family of putative transcriptional regulators, respectively. Here, we have performed experiments to understand the role of MenT2, a toxin from the type IV TA system, in mycobacterial physiology and disease pathogenesis. The ectopic expression of MenT2 using inducible vectors does not inhibit bacterial growth in liquid cultures. Bioinformatic and molecular modelling analysis suggested that the M. tuberculosis genome has an alternative start site upstream of the annotated menT2 gene. The overexpression of the reannotated MenT2 resulted in moderate growth inhibition of Mycobacterium smegmatis. We show that both menT2 and menA2 transcript levels are increased when M. tuberculosis is exposed to nitrosative stress, in vitro. When compared to the survival of the wild-type and the complemented strain, the ΔmenT2 mutant strain of M. tuberculosis was more resistant to being killed by nitrosative stress. However, the survival of both the ΔmenT2 mutant and the wild-type strain was similar in macrophages and when exposed to other stress conditions. Here, we show that MenT2 is required for the establishment of disease in guinea pigs. Gross pathology and histopathology analysis of lung tissues from guinea pigs infected with the ∆menT2 strain revealed significantly reduced tissue damage and inflammation. In summary, these results provide new insights into the role of MenT2 in mycobacterial pathogenesis.

PAT: a comprehensive database of prokaryotic antimicrobial toxins

Antimicrobial toxins help prokaryotes win competitive advantages in intraspecific or interspecific conflicts and are also a critical factor affecting the pathogenicity of many pathogens that threaten human health. Although many studies have revealed that antagonism based on antimicrobial toxins plays a central role in prokaryotic life, a database on antimicrobial toxins remains lacking. Here, we present the prokaryotic antimicrobial toxin database (PAT, http://bioinfo.qd.sdu.edu.cn/PAT/), a comprehensive data resource collection on experimentally validated antimicrobial toxins. PAT has organized information, derived from the reported literature, on antimicrobial toxins, as well as the corresponding immunity proteins, delivery mechanisms, toxin activities, structural characteristics, sequences, etc. Moreover, we also predict potential antimicrobial toxins in prokaryotic reference genomes and show the taxonomic information and environmental distribution of typical antimicrobial toxins. These details have been fully incorporated into the PAT database, where users can browse, search, download, analyse and view informative statistics and detailed information. PAT resources have already been used in our prediction and identification of prokaryotic antimicrobial toxins and may contribute to promoting the efficient investigation of antimicrobial toxin functions, the discovery of novel antimicrobial toxins, and an improved understanding of the biological roles and significance of these toxins.

Genomic Characterization and Antimicrobial Susceptibility of Dromedary-Associated Staphylococcaceae from the Horn of Africa

Members of the Staphylococcaceae family, particularly those of the genus Staphylococcus, encompass important human and animal pathogens. We collected and characterized Staphylococcaceae strains from apparently healthy and diseased camels (n = 84) and cattle (n = 7) in Somalia and Kenya. We phenotypically characterized the strains, including their antimicrobial inhibitory concentrations. Then, we sequenced their genomes using long-read sequencing, closed their genomes, and subsequently compared and mapped their virulence- and resistance-associated gene pools. Genome-based phylogenetics revealed 13 known Staphylococcaceae and at least two novel species. East African strains of different species encompassed novel sequence types and phylogenetically distant clades. About one-third of the strains had non-wild-type MICs. They were resistant to at least one of the following antimicrobials: tetracycline, benzylpenicillin, oxacillin, erythromycin, clindamycin, trimethoprim, gentamicin, or streptomycin, encoded by tet(K), blaZ/bla_ARL, mecA/mecA1, msrA/mphC, salA, dfrG, aacA-aphD, and str, respectively. We identified the first methicillin- and multidrug-resistant camel S. epidermidis strain of sequence type (ST) 1136 in East Africa. The pool of virulence-encoding genes was largest in the S. aureus strains, as expected, although other rather commensal strains contained distinct virulence-encoding genes. We identified toxin-antitoxin (TA) systems such as the hicA/hicB and abiEii/abiEi families, reported here for the first time for certain species of Staphylococcaceae. All strains contained at least one intact prophage sequence, mainly belonging to the Siphoviridae family. We pinpointed potential horizontal gene transfers between camel and cattle strains and also across distinct Staphylococcaceae clades and species. IMPORTANCE Camels are a high value and crucial livestock species in arid and semiarid regions of Africa and gain importance giving the impact of climate change on traditional livestock species. Our current knowledge with respect to Staphylococcaceae infecting camels is very limited compared to that for other livestock species. Better knowledge will foster the development of specific diagnostic assays, guide promising antimicrobial treatment options, and inform about potential zoonotic risks. We characterized 84 Staphylococcaceae strains isolated from camels with respect to their antimicrobial resistance and virulence traits. We detected potentially novel Staphylococcus species, resistances to different classes of antimicrobials, and the first camel multidrug-resistant S. epidermidis strain of sequence type 1136.

Genomic analysis of plasmid content in food isolates of E. coli strongly supports its role as a reservoir for the horizontal transfer of virulence and antibiotic resistance genes

The link between E. coli strains contaminating foods and human disease is unclear, with some reports supporting a direct transmission of pathogenic strains via food and others highlighting their role as reservoirs for resistance and virulence genes. Here we take a genomics approach, analyzing a large set of fully-assembled genomic sequences from E. coli available in GenBank. Most of the strains isolated in food are more closely related to each other than to clinical strains, arguing against a frequent direct transmission of pathogenic strains from food to the clinic. We also provide strong evidence of genetic exchanges between food and clinical strains that are facilitated by plasmids. This is based on an overlapped representation of virulence and resistance genes in plasmids isolated from these two sources. We identify clusters of phylogenetically-related plasmids that are largely responsible for the observed overlap and see evidence of specialization, with some food plasmid clusters preferentially transferring virulence factors over resistance genes. Consistent with these observations, food plasmids have a high mobilization potential based on their plasmid taxonomic unit classification and on an analysis of mobilization gene content. We report antibiotic resistance genes of high clinical relevance and their specific incompatibility group associations. Finally, we also report a striking enrichment for adhesins in food plasmids and their association with specific IncF replicon subtypes. The identification of food plasmids with specific markers (Inc and PTU combinations) as mediators of horizontal transfer between food and clinical strains opens new research avenues and should assist with the design of surveillance strategies.

Characterization of toxin-antitoxin systems from public sequencing data: A case study in Pseudomonas aeruginosa

The toxin-antitoxin (TA) system is a widely distributed group of genetic modules that play important roles in the life of prokaryotes, with mobile genetic elements (MGEs) contributing to the dissemination of antibiotic resistance gene (ARG). The diversity and richness of TA systems in Pseudomonas aeruginosa, as one of the bacterial species with ARGs, have not yet been completely demonstrated. In this study, we explored the TA systems from the public genomic sequencing data and genome sequences. A small scale of genomic sequencing data in 281 isolates was selected from the NCBI SRA database, reassembling the genomes of these isolates led to the findings of abundant TA homologs. Furthermore, remapping these identified TA modules on 5,437 genome/draft genomes uncovers a great diversity of TA modules in P. aeruginosa. Moreover, manual inspection revealed several TA systems that were not yet reported in P. aeruginosa including the hok-sok, cptA-cptB, cbeA-cbtA, tomB-hha, and ryeA-sdsR. Additional annotation revealed that a large number of MGEs were closely distributed with TA. Also, 16% of ARGs are located relatively close to TA. Our work confirmed a wealth of TA genes in the unexplored P. aeruginosa pan-genomes, expanded the knowledge on P. aeruginosa, and provided methodological tips on large-scale data mining for future studies. The co-occurrence of MGE, ARG, and TA may indicate a potential interaction in their dissemination.

Evolution of plasmid mobility: origin and fate of conjugative and non-conjugative plasmids

Conjugation drives the horizontal transfer of adaptive traits across prokaryotes. One-fourth of the plasmids encode the functions necessary to conjugate autonomously, the others being eventually mobilizable by conjugation. To understand the evolution of plasmid mobility, we studied plasmid size, gene repertoires, and conjugation-related genes. Plasmid gene repertoires were found to vary rapidly in relation to the evolutionary rate of relaxases, for example, most pairs of plasmids with 95% identical relaxases have fewer than 50% of homologs. Among 249 recent transitions of mobility type, we observed a clear excess of plasmids losing the capacity to conjugate. These transitions are associated with even greater changes in gene repertoires, possibly mediated by transposable elements, including pseudogenization of the conjugation locus, exchange of replicases reducing the problem of incompatibility, and extensive loss of other genes. At the microevolutionary scale of plasmid taxonomy, transitions of mobility type sometimes result in the creation of novel taxonomic units. Interestingly, most transitions from conjugative to mobilizable plasmids seem to be lost in the long term. This suggests a source-sink dynamic, where conjugative plasmids generate nonconjugative plasmids that tend to be poorly adapted and are frequently lost. Still, in some cases, these relaxases seem to have evolved to become efficient at plasmid mobilization in trans, possibly by hijacking multiple conjugative systems. This resulted in specialized relaxases of mobilizable plasmids. In conclusion, the evolution of plasmid mobility is frequent, shapes the patterns of gene flow in bacteria, the dynamics of gene repertoires, and the ecology of plasmids.

Substrate recognition and cryo-EM structure of the ribosome-bound TAC toxin of Mycobacterium tuberculosis

Toxins of toxin-antitoxin systems use diverse mechanisms to control bacterial growth. Here, we focus on the deleterious toxin of the atypical tripartite toxin-antitoxin-chaperone (TAC) system of Mycobacterium tuberculosis, whose inhibition requires the concerted action of the antitoxin and its dedicated SecB-like chaperone. We show that the TAC toxin is a bona fide ribonuclease and identify exact cleavage sites in mRNA targets on a transcriptome-wide scale in vivo. mRNA cleavage by the toxin occurs after the second nucleotide of the ribosomal A-site codon during translation, with a strong preference for CCA codons in vivo. Finally, we report the cryo-EM structure of the ribosome-bound TAC toxin in the presence of native M. tuberculosis cspA mRNA, revealing the specific mechanism by which the TAC toxin interacts with the ribosome and the tRNA in the P-site to cleave its mRNA target.

Toxin-antitoxin systems are widespread in bacteria. Here the authors present structures of M. tuberculosis HigB^TAC alone and bound to the ribosome in the presence of native cspA mRNA, shedding light on its mechanism of translation inhibition.

Functional and Biochemical Characterization of the MazEF6 Toxin-Antitoxin System of Mycobacterium tuberculosis

The Mycobacterium tuberculosis genome harbors nine toxin-antitoxin (TA) systems that are members of the mazEF family, unlike other prokaryotes, which have only one or two. Although the overall tertiary folds of MazF toxins are predicted to be similar, it is unclear how they recognize structurally different RNAs and antitoxins with divergent sequence specificity. Here, we have expressed and purified the individual components and complex of the MazEF6 TA system from M. tuberculosis. Size exclusion chromatography-multiangle light scattering (SEC-MALS) was performed to determine the oligomerization status of the toxin, antitoxin, and the complex in different stoichiometric ratios. The relative stabilities of the proteins were determined by nano-differential scanning fluorimetry (nano-DSF). Microscale thermophoresis (MST) and yeast surface display (YSD) were performed to measure the relative affinities between the cognate toxin-antitoxin partners. The interaction between MazEF6 complexes and cognate promoter DNA was also studied using MST. Analysis of paired-end RNA sequencing data revealed that the overexpression of MazF6 resulted in differential expression of 323 transcripts in M. tuberculosis. Network analysis was performed to identify the nodes from the top-response network. The analysis of mRNA protection ratios resulted in identification of putative MazF6 cleavage site in its native host, M. tuberculosis. IMPORTANCE M. tuberculosis harbors a large number of type II toxin-antitoxin (TA) systems, the exact roles for most of which are unclear. Prior studies have reported that overexpression of several of these type II toxins inhibits bacterial growth and contributes to the formation of drug-tolerant populations in vitro. To obtain insights into M. tuberculosis MazEF6 type II TA system function, we determined stability, oligomeric states, and binding affinities of cognate partners with each other and with their promoter operator DNA. Using RNA-seq data obtained from M. tuberculosis overexpression strains, we have identified putative MazF6 cleavage sites and targets in its native, cellular context.

Linear plasmids in Klebsiella and other Enterobacteriaceae

Linear plasmids are extrachromosomal DNA elements that have been found in a small number of bacterial species. To date, the only linear plasmids described in the family Enterobacteriaceae belong to Salmonella, first found in Salmonella enterica Typhi. Here, we describe a collection of 12 isolates of the Klebsiella pneumoniae species complex in which we identified linear plasmids. Screening of assembly graphs assembled from public read sets identified linear plasmid structures in a further 13 K. pneumoniae species complex genomes. We used these 25 linear plasmid sequences to query all bacterial genome assemblies in the National Center for Biotechnology Information database, and discovered an additional 61 linear plasmid sequences in a variety of Enterobacteriaceae species. Gene content analysis divided these plasmids into five distinct phylogroups, with very few genes shared across more than two phylogroups. The majority of linear plasmid-encoded genes are of unknown function; however, each phylogroup carried its own unique toxin–antitoxin system and genes with homology to those encoding the ParAB plasmid stability system. Passage in vitro of the 12 linear plasmid-carrying Klebsiella isolates in our collection (which include representatives of all five phylogroups) indicated that these linear plasmids can be stably maintained, and our data suggest they can transmit between K. pneumoniae strains (including members of globally disseminated multidrug-resistant clones) and also between diverse Enterobacteriaceae species. The linear plasmid sequences, and representative isolates harbouring them, are made available as a resource to facilitate future studies on the evolution and function of these novel plasmids.

Inter-species geographic signatures for tracing horizontal gene transfer and long-term persistence of carbapenem resistance

BACKGROUND

Carbapenem-resistant Enterobacterales (CRE) are an urgent global health threat. Inferring the dynamics of local CRE dissemination is currently limited by our inability to confidently trace the spread of resistance determinants to unrelated bacterial hosts. Whole-genome sequence comparison is useful for identifying CRE clonal transmission and outbreaks, but high-frequency horizontal gene transfer (HGT) of carbapenem resistance genes and subsequent genome rearrangement complicate tracing the local persistence and mobilization of these genes across organisms.

METHODS

To overcome this limitation, we developed a new approach to identify recent HGT of large, near-identical plasmid segments across species boundaries, which also allowed us to overcome technical challenges with genome assembly. We applied this to complete and near-complete genome assemblies to examine the local spread of CRE in a systematic, prospective collection of all CRE, as well as time- and species-matched carbapenem-susceptible Enterobacterales, isolated from patients from four US hospitals over nearly 5 years.

RESULTS

Our CRE collection comprised a diverse range of species, lineages, and carbapenem resistance mechanisms, many of which were encoded on a variety of promiscuous plasmid types. We found and quantified rearrangement, persistence, and repeated transfer of plasmid segments, including those harboring carbapenemases, between organisms over multiple years. Some plasmid segments were found to be strongly associated with specific locales, thus representing geographic signatures that make it possible to trace recent and localized HGT events. Functional analysis of these signatures revealed genes commonly found in plasmids of nosocomial pathogens, such as functions required for plasmid retention and spread, as well survival against a variety of antibiotic and antiseptics common to the hospital environment.

CONCLUSIONS

Collectively, the framework we developed provides a clearer, high-resolution picture of the epidemiology of antibiotic resistance importation, spread, and persistence in patients and healthcare networks.

A hyperpromiscuous antitoxin protein domain for the neutralization of diverse toxin domains

Toxin–antitoxin systems are enigmatic and diverse elements of bacterial and bacteriophage genomes. We have uncovered remarkable versatility in an antitoxin protein domain that has evolved to neutralize dozens of different toxin domains. We find that antitoxins carrying this domain—Panacea—form complexes with their cognate toxins, indicating a direct neutralization mechanism, and that Panacea can be evolved to neutralize a noncognate and nonhomologous toxin with just two amino acid substitutions. This raises the possibility that this domain could be an adaptable universal or semi-universal protein neutralizer with significant biotechnological and medical potential.

Toxin–antitoxin (TA) gene pairs are ubiquitous in microbial chromosomal genomes and plasmids as well as temperate bacteriophages. They act as regulatory switches, with the toxin limiting the growth of bacteria and archaea by compromising diverse essential cellular targets and the antitoxin counteracting the toxic effect. To uncover previously uncharted TA diversity across microbes and bacteriophages, we analyzed the conservation of genomic neighborhoods using our computational tool FlaGs (for flanking genes), which allows high-throughput detection of TA-like operons. Focusing on the widespread but poorly experimentally characterized antitoxin domain DUF4065, our in silico analyses indicated that DUF4065-containing proteins serve as broadly distributed antitoxin components in putative TA-like operons with dozens of different toxic domains with multiple different folds. Given the versatility of DUF4065, we have named the domain Panacea (and proteins containing the domain, PanA) after the Greek goddess of universal remedy. We have experimentally validated nine PanA-neutralized TA pairs. While the majority of validated PanA-neutralized toxins act as translation inhibitors or membrane disruptors, a putative nucleotide cyclase toxin from a Burkholderia prophage compromises transcription and translation as well as inducing RelA-dependent accumulation of the nucleotide alarmone (p)ppGpp. We find that Panacea-containing antitoxins form a complex with their diverse cognate toxins, characteristic of the direct neutralization mechanisms employed by Type II TA systems. Finally, through directed evolution, we have selected PanA variants that can neutralize noncognate TA toxins, thus experimentally demonstrating the evolutionary plasticity of this hyperpromiscuous antitoxin domain.

Functional characterization of HigBA toxin-antitoxin system in an Arctic bacterium, Bosea sp. PAMC 26642

Toxin-antitoxin (TA) systems are growth-controlling genetic elements consisting of an intracellular toxin protein and its cognate antitoxin. TA systems have been spread among microbial genomes through horizontal gene transfer and are now prevalent in most bacterial and archaeal genomes. Under normal growth conditions, antitoxins tightly counteract the activity of the toxins. Upon stresses, antitoxins are inactivated, releasing activated toxins, which induce growth arrest or cell death. In this study, among nine functional TA modules in Bosea sp. PAMC 26642 living in Arctic lichen, we investigated the functionality of BoHigBA2. BohigBA2 is located close to a genomic island and adjacent to flagellar gene clusters. The expression of BohigB2 induced the inhibition of E. coli growth at 37°C, which was more manifest at 18°C, and this growth defect was reversed when BohigA2 was co-expressed, suggesting that this BoHigBA2 module might be an active TA module in Bosea sp. PAMC 26642. Live/dead staining and viable count analyses revealed that the BoHigB2 toxin had a bactericidal effect, causing cell death. Furthermore, we demonstrated that BoHigB2 possessed mRNA-specific ribonuclease activity on various mRNAs and cleaved only mRNAs being translated, which might impede overall translation and consequently lead to cell death. Our study provides the insight to understand the cold adaptation of Bosea sp. PAMC 26642 living in the Arctic.

High-resolution genomic surveillance elucidates a multilayered hierarchical transfer of resistance between WWTP- and human/animal-associated bacteria

BACKGROUND

Our interconnected world and the ability of bacteria to quickly swap antibiotic resistance genes (ARGs) make it particularly important to establish the epidemiological links of multidrug resistance (MDR) transfer between wastewater treatment plant (WWTP)- and human/animal-associated bacteria, under the One Health framework. However, evidence of ARGs exchange and potential factors that contribute to this transfer remain limited.

RESULTS

Here, by combining culture-based population genomics and genetic comparisons with publicly available datasets, we reconstructed the complete genomes of 82 multidrug-resistant isolates from WWTPs and found that most WWTP-associated isolates were genetically distinct from their closest human/animal-associated relatives currently available in the public database. Even in the minority of lineages that were closely related, WWTP-associated isolates were characterized by quite different plasmid compositions. We identified a high diversity of circular plasmids (264 in total, of which 141 were potentially novel), which served as the main source of resistance, and showed potential horizontal transfer of ARG-bearing plasmids between WWTP- and humans/animal-associated bacteria. Notably, the potentially transferred ARGs and virulence factors (VFs) with different genetic backgrounds were closely associated with flanking insertion sequences (ISs), suggesting the importance of synergy between plasmids and ISs in mediating a multilayered hierarchical transfer of MDR and potentiating the emergence of MDR-hypervirulent clones.

CONCLUSION

Our findings advance the current efforts to establish potential epidemiological links of MDR transmission between WWTP- and human/animal-associated bacteria. Plasmids play an important role in mediating the transfer of ARGs and the IS-associated ARGs that are carried by conjugative plasmids should be prioritized to tackle the spread of resistance. Video Abstract.