Comprehensive comparative-genomic analysis of type 2 toxin-antitoxin systems and related mobile stress response systems in prokaryotes

RES-Xre toxin-antitoxin locus knaAT maintains the stability of the virulence plasmid in Klebsiella pneumoniae

Hypervirulent Klebsiella pneumoniae isolates have been increasingly reported worldwide, especially hypervirulent drug-resistant variants owing to the acquisition of a mobilizable virulence plasmid by a carbapenem-resistant strain. This pLVPK-like mobilizable plasmid encodes various virulence factors; however, information about its genetic stability is lacking. This study aimed to investigate the type II toxin-antitoxin (TA) modules that facilitate the virulence plasmid to remain stable in K. pneumoniae. More than 3,000 TA loci in 2,000 K. pneumoniae plasmids were examined for their relationship with plasmid cargo genes. TA loci from the RES-Xre family were highly correlated with virulence plasmids of hypervirulent K. pneumoniae. Overexpression of the RES toxin KnaT, encoded by the virulence plasmid-carrying RES-Xre locus knaAT, halts the cell growth of K. pneumoniae and E. coli, whereas co-expression of the cognate Xre antitoxin KnaA neutralizes the toxicity of KnaT. knaA and knaT were co-transcribed, representing the characteristics of a type II TA module. The knaAT deletion mutation gradually lost its virulence plasmid in K. pneumoniae, whereas the stability of the plasmid in E. coli was enhanced by adding knaAT, which revealed that the knaAT operon maintained the genetic stability of the large virulence plasmid in K. pneumoniae. String tests and mouse lethality assays subsequently confirmed that a loss of the virulence plasmid resulted in reduced pathogenicity of K. pneumoniae. These findings provide important insights into the role of the RES-Xre TA pair in stabilizing virulence plasmids and disseminating virulence genes in K. pneumoniae.

Evaluating the antibacterial effect of meropenem-loaded chitosan/sodium tripolyphosphate (TPP) nanoparticles on Acinetobacter baumannii isolated from hospitalized patients

BACKGROUND

Acinetobacter baumannii is a health threat due to its antibiotic resistance. Herein, antibiotic susceptibility and its association with the Toxin-antitoxin (TA) system genes in A. baumannii clinical isolates from Iran were investigated. Next, we prepared meropenem-loaded chitosan nanoparticles (MP-CS) and investigated their antibacterial effects against meropenem-susceptible bacterial isolates.

METHODS

Out of 240 clinical specimens, 60 A. baumannii isolates were assessed. Antibiotic resistance of the isolates against conventional antibiotics was determined alongside investigating the presence of three TA system genes (mazEF, relBE, and higBA). Chitosan nanoparticles were characterized in terms of size, zeta potential, encapsulation efficiency, and meropenem release activity. Their antibacterial effects were assessed using the well diffusion method, minimum inhibitory concentration (MIC), and colony-forming unit (CFU) counting. Their cytotoxic effects and biocompatibility index were determined via the MTT, LDH, and ROS formation assays.

RESULTS

Ampicillin, ceftazidime, and colistin were the least effective, and amikacin and tobramycin were the most effective antibiotics. Out of the 60 isolates, 10 (16.7%), 5 (8.3%), and 45 (75%) were multidrug-resistant (MDR), extensively drug-resistant (XDR), and pandrug-resistant (PDR), respectively. TA system genes had no significant effect on antibiotic resistance. MP-CS nanoparticles demonstrated an average size of 191.5 and zeta potential of 27.3 mV alongside a maximum encapsulation efficiency of 88.32% and release rate of 69.57%. MP-CS nanoparticles mediated similar antibacterial effects, as compared with free meropenem, against the A. baumannii isolates with significantly lower levels of meropenem. MP-CS nanoparticles remarkably prevented A549 and NCI-H292 cell infection by the A. baumannii isolates alongside demonstrating a favorable biocompatibility index.

CONCLUSION

Antibiotic-loaded nanoparticles should be further designed and investigated to increase their antibacterial effect against A. baumannii and assess their safety and applicability in vivo settings.

The Chlamydia-related Waddlia chondrophila encodes functional type II toxin-antitoxin systems

Bacterial toxin-antitoxin (TA) systems are widespread in chromosomes and plasmids of free-living microorganisms, but only a few have been identified in obligate intracellular species. We found seven putative type II TA modules in Waddlia chondrophila, a Chlamydia-related species that is able to infect a very broad series of eukaryotic hosts, ranging from protists to mammalian cells. The RNA levels of Waddlia TA systems are significantly upregulated by iron starvation and novobiocin, but they are not affected by antibiotics such as β-lactams and glycopeptides, which suggests different mechanisms underlying stress responses. Five of the identified TA modules, including HigBA1 and MazEF1, encoded on the Waddlia cryptic plasmid, proved to be functional when expressed in a heterologous host. TA systems have been associated with the maintenance of mobile genetic elements, bacterial defense against bacteriophages, and persistence upon exposure to adverse conditions. As their RNA levels are upregulated upon exposure to adverse conditions, Waddlia TA modules may be involved in survival to stress. Moreover, as Waddlia can infect a wide range of hosts including free-living amoebae, TA modules could also represent an innate immunity system to fight against bacteriophages and other microorganisms with which Waddlia has to share its replicative niche.IMPORTANCEThe response to adverse conditions, such as exposure to antibiotics, nutrient starvation and competition with other microorganisms, is essential for the survival of a bacterial population. TA systems are modules composed of two elements, a toxic protein and an antitoxin (protein or RNA) that counteracts the toxin. Although many aspects of TA biological functions still await to be elucidated, TAs have often been implicated in bacterial response to stress, including the response to nutrient starvation, antibiotic treatment and bacteriophage infection. TAs are ubiquitous in free-living bacteria but rare in obligate intracellular species such as chlamydiae. We identified functional TA systems in Waddlia chondrophila, a chlamydial species with a strikingly broad host range compared to other chlamydiae. Our work contributes to understand how obligate intracellular bacteria react to adverse conditions that might arise from competition with other viruses/bacteria for the same replicative niche and would threaten their ability to replicate.

Toxinome-the bacterial protein toxin database

IMPORTANCE

Microbes use protein toxins as important tools to attack neighboring cells, microbial or eukaryotic, and for self-killing when attacked by viruses. These toxins work through different mechanisms to inhibit cell growth or kill cells. Microbes also use antitoxin proteins to neutralize the toxin activities. Here, we developed a comprehensive database called Toxinome of nearly two million toxins and antitoxins that are encoded in 59,475 bacterial genomes. We described the distribution of bacterial toxins and identified that they are depleted by bacteria that live in hot and cold temperatures. We found 5,161 cases in which toxins and antitoxins are densely clustered in bacterial genomes and termed these areas "Toxin Islands." The Toxinome database is a useful resource for anyone interested in toxin biology and evolution, and it can guide the discovery of new toxins.

Toxin:antitoxin ratio sensing autoregulation of the Vibrio cholerae parDE2 module

The parDE family of toxin-antitoxin (TA) operons is ubiquitous in bacterial genomes and, in Vibrio cholerae, is an essential component to maintain the presence of chromosome II. Here, we show that transcription of the V. cholerae parDE2 (VcparDE) operon is regulated in a toxin:antitoxin ratio-dependent manner using a molecular mechanism distinct from other type II TA systems. The repressor of the operon is identified as an assembly with a 6:2 stoichiometry with three interacting ParD2 dimers bridged by two ParE2 monomers. This assembly docks to a three-site operator containing 5'- GGTA-3' motifs. Saturation of this TA complex with ParE2 toxin results in disruption of the interface between ParD2 dimers and the formation of a TA complex of 2:2 stoichiometry. The latter is operator binding-incompetent as it is incompatible with the required spacing of the ParD2 dimers on the operator.

A novel double-ribonuclease toxin-antitoxin system linked to the stress response and survival of Acidovorax citrulli

IMPORTANCE

Bacterial fruit blotch (BFB), which is caused by the seed-borne bacterium Acidovorax citrulli, is a devastating disease affecting cucurbit crops throughout the world. Although seed fermentation and treatment with disinfectants can provide effective management of BFB, they cannot completely guarantee pathogen-free seedstock, which suggests that A. citrulli is a highly stress-resistant pathogen. Toxin-antitoxin (TA) systems are common among a diverse range of bacteria and have been reported to play a role in bacterial stress response. However, there is currently much debate about the relationship between TA systems and stress response in bacteria. The current study characterized a novel TA system (Aave_1720-Aave_1719) from A. citrulli that affects both biofilm formation and survival in response to sodium hypochlorite stress. The mechanism of neutralization differed from typical TA systems as two separate mechanisms were associated with the antitoxin, which exhibited characteristics of both type II and type V TA systems. The Aave_1720-Aave_1719 system described here also constitutes the first known report of a double-ribonuclease TA system in bacteria, which expands our understanding of the range of regulatory mechanisms utilized by bacterial TA systems, providing new insight into the survival of A. citrulli in response to stress.

Exploring the Accessory Genome of Multidrug-Resistant Rhodococcus equi Clone 2287

Decades of antimicrobial overuse to treat respiratory disease in foals have promoted the emergence and spread of zoonotic multidrug-resistant (MDR) Rhodococcus equi worldwide. Three main R. equi MDR clonal populations-2287, G2106, and G2017-have been identified so far. However, only clones 2287 and G2016 have been isolated from sick animals, with clone 2287 being the main MDR R. equi recovered. The genetic mechanisms that make this MDR clone superior to the others at infecting foals are still unknown. Here, we performed a deep genetic characterization of the accessory genomes of 207 R. equi isolates, and we describe IME2287, a novel genetic element in the accessory genome of clone 2287, potentially involved in the maintenance and spread of this MDR population over time. IME2287 is a putative self-replicative integrative mobilizable element (IME) carrying a DNA replication and partitioning operon and genes encoding its excision and integration from the R. equi genome via a serine recombinase. Additionally, IME2287 encodes a protein containing a Toll/interleukin-1 receptor (TIR) domain that may inhibit TLR-mediated NF-kB signaling in the host and a toxin-antitoxin (TA) system, whose orthologs have been associated with antibiotic resistance/tolerance, virulence, pathogenicity islands, bacterial persistence, and pathogen trafficking. This new set of genes may explain the success of clone 2287 over the other MDR R. equi clones.

Functional diversification despite structural congruence in the HipBST toxin-antitoxin system of Legionella pneumophila

IMPORTANCE

Toxin-antitoxin (TA) systems are parasitic genetic elements found in almost all bacterial genomes. They are exchanged horizontally between cells and are typically poorly conserved across closely related strains and species. Here, we report the characterization of a tripartite TA system in the bacterial pathogen Legionella pneumophila that is highly conserved across Legionella species genomes. This system (denoted HipBSTLp) is a distant homolog of the recently discovered split-HipA system in Escherichia coli (HipBSTEc). We present bioinformatic, molecular, and structural analyses of the divergence between these two systems and the functionality of this newly described TA system family. Furthermore, we provide evidence to refute previous claims that the toxin in this system (HipTLp) possesses bifunctionality as an L. pneumophila virulence protein. Overall, this work expands our understanding of the split-HipA system architecture and illustrates the potential for undiscovered biology in these abundant genetic elements.

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.

Oligomerization-mediated activation of a short prokaryotic Argonaute

Although eukaryotic and long prokaryotic Argonaute proteins (pAgos) cleave nucleic acids, some short pAgos lack nuclease activity and hydrolyse NAD(P)+ to induce bacterial cell death1. Here we present a hierarchical activation pathway for SPARTA, a short pAgo consisting of an Argonaute (Ago) protein and TIR-APAZ, an associated protein2. SPARTA progresses through distinct oligomeric forms, including a monomeric apo state, a monomeric RNA-DNA-bound state, two dimeric RNA-DNA-bound states and a tetrameric RNA-DNA-bound active state. These snapshots together identify oligomerization as a mechanistic principle of SPARTA activation. The RNA-DNA-binding channel of apo inactive SPARTA is occupied by an auto-inhibitory motif in TIR-APAZ. After the binding of RNA-DNA, SPARTA transitions from a monomer to a symmetric dimer and then an asymmetric dimer, in which two TIR domains interact through charge and shape complementarity. Next, two dimers assemble into a tetramer with a central TIR cluster responsible for hydrolysing NAD(P)+. In addition, we observe unique features of interactions between SPARTA and RNA-DNA, including competition between the DNA 3' end and the auto-inhibitory motif, interactions between the RNA G2 nucleotide and Ago, and splaying of the RNA-DNA duplex by two loops exclusive to short pAgos. Together, our findings provide a mechanistic basis for the activation of short pAgos, a large section of the Ago superfamily.

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/.

Role of fourteen XRE-DUF397 pairs from Streptomyces coelicolor as regulators of antibiotic production and differentiation. New players in a complex regulatory network

Bacteria of the genus Streptomyces have a plethora of transcriptional regulators, among which the xenobiotic response element (XRE) plays an important role. In this organism, XRE regulators are often followed downstream by small proteins of unknown function containing a DUF397 domain. It has been proposed that XRE/DUF397 pairs constitute type II toxin-antitoxin (TA) systems. However, previous work carried out by our group has shown that one of these systems is a strong activator of antibiotic production in S. coelicolor and other Streptomyces species. In this work, we have studied the overexpression of fourteen XRE/DUF397 pairs present in the S. coelicolor genome and found that none behave as a type II TA system. Instead, they act as pleiotropic regulators affecting, in a dependent manner, antibiotic production and morphological differentiation on different culture media. After deleting, individually, six XRE/DUF397 pairs (those systems producing more notable phenotypic changes when overexpressed: SCO2246/45, SCO2253/52, SCO4176/77, SCO4678/79, SCO6236/35, and SCO7615/16), the pair SCO7615/16 was identified as producing the most dramatic differences as compared to the wild-type strain. The SCO7615/16 mutant had a different phenotype on each of the media tested (R2YE, LB, NMMP, YEPD, and MSA). In particular, on R2YE and YEPD media, a bald phenotype was observed even after 7 days, with little or no actinorhodin (ACT) production. Lower ACT production was also observed on LB medium, but the bacteria were able to produce aerial mycelium. On NMMP medium, the mutant produced a larger amount of ACT as compared with the wild-type strain.

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.

Treponema denticola Has the Potential to Cause Neurodegeneration in the Midbrain via the Periodontal Route of Infection-Narrative Review

Alzheimer's disease (AD) is a neurodegenerative disease and the most common example of dementia. The neuropathological features of AD are the abnormal deposition of extracellular amyloid-β (Aβ) and intraneuronal neurofibrillary tangles with hyperphosphorylated tau protein. It is recognized that AD starts in the frontal cerebral cortex, and then it progresses to the entorhinal cortex, the hippocampus, and the rest of the brain. However, some studies on animals suggest that AD could also progress in the reverse order starting from the midbrain and then spreading to the frontal cortex. Spirochetes are neurotrophic: From a peripheral route of infection, they can reach the brain via the midbrain. Their direct and indirect effect via the interaction of their virulence factors and the microglia potentially leads to the host peripheral nerve, the midbrain (especially the locus coeruleus), and cortical damage. On this basis, this review aims to discuss the hypothesis of the ability of Treponema denticola to damage the peripheral axons in the periodontal ligament, to evade the complemental pathway and microglial immune response, to determine the cytoskeletal impairment and therefore causing the axonal transport disruption, an altered mitochondrial migration and the consequent neuronal apoptosis. Further insights about the central neurodegeneration mechanism and Treponema denticola's resistance to the immune response when aggregated in biofilm and its quorum sensing are suggested as a pathogenetic model for the advanced stages of AD.

Structural insights of the toxin-antitoxin system VPA0770-VPA0769 in Vibrio parahaemolyticus

Toxin-antitoxin (TA) systems are involved in both normal bacterial physiology and pathogenicity, including gene regulation, antibiotic resistance, and bacteria persistence under stressful environments. In pathogenic Vibrio parahaemolyticus, however, TA interaction and assembly remain largely unknown. In this work, we identified a new RES-Xre type II TA module, encoded by gene cluster vpa0770-vpa0769 on chromosome II of V. parahaemolyticus. Ectopic expression of the VPA0770 toxin rapidly arrests the growth of E. coli cells, which can be neutralized by co-expression of the VPA0769 antitoxin. To decipher the action mechanism, we determined the crystal structure of the VPA0770-VPA0769 TA complex. VPA0770 and VPA0769 proteins can assemble into two types of large complexes, a W-shaped hetero-hexamer and a donut-like hetero-dodecamer, in a concentration-dependent manner in solution. Disruption of the TA interface results in a loss of the antitoxic phenotype. The toxicity of the VPA0770 toxin, which harbors a NAD⁺-binding pocket, may be largely ascribed to its highly effective capability to degrade intracellular NAD⁺. Our study provides a structural basis for a better understanding of diverse molecular mechanisms employed by human pathogens.

Isolation and characterization of a novel phage belonging to a new genus against Vibrio parahaemolyticus

BACKGROUND

Vibrio parahaemolyticus is a major foodborne pathogen that contaminates aquatic products and causes great economic losses to aquaculture. Because of the emergence of multidrug-resistant V. parahaemolyticus strains, bacteriophages are considered promising agents for their biocontrol as an alternative or supplement to antibiotics. In this study, a lytic vibriophage, vB_VpaM_R16F (R16F), infecting V. parahaemolyticus 1.1997^T was isolated, characterized and evaluated for its biocontrol potential.

METHODS

A vibriophage R16F was isolated from sewage from a seafood market with the double-layer agar method. R16F was studied by transmission electron microscopy, host range, sensitivity of phage particles to chloroform, one-step growth curve and lytic activity. The phage genome was sequenced and in-depth characterized, including phylogenetic and taxonomic analysis.

RESULTS

R16F belongs to the myovirus morphotype and infects V. parahaemolyticus, but not nine other Vibrio spp. As characterized by determining its host range, one-step growth curve, and lytic activity, phage R16F was found to highly effective in lysing host cells with a short latent period (< 10 min) and a small burst size (13 plaque-forming units). R16F has a linear double-stranded DNA with genome size 139,011 bp and a G + C content of 35.21%. Phylogenetic and intergenomic nucleotide sequence similarity analysis revealed that R16F is distinct from currently known vibriophages and belongs to a novel genus. Several genes (e.g., encoding ultraviolet damage endonuclease and endolysin) that may enhance environmental competitiveness were found in the genome of R16F, while no antibiotic resistance- or virulence factor-related gene was detected.

CONCLUSIONS

In consideration of its biological and genetic properties, this newly discovered phage R16F belongs to a novel genus and may be a potential alternate biocontrol agent.

Manesh MJ, Sivabalasarma S, Albers SV, Kelly RM. Stay or Go: Sulfolobales Biofilm Dispersal Is Dependent on a Bifunctional VapB Antitoxin

A type II VapB14 antitoxin regulates biofilm dispersal in the archaeal thermoacidophile Sulfolobus acidocaldarius through traditional toxin neutralization but also through noncanonical transcriptional regulation. Type II VapC toxins are ribonucleases that are neutralized by their proteinaceous cognate type II VapB antitoxin. VapB antitoxins have a flexible tail at their C terminus that covers the toxin's active site, neutralizing its activity. VapB antitoxins also have a DNA-binding domain at their N terminus that allows them to autorepress not only their own promoters but also distal targets. VapB14 antitoxin gene deletion in S. acidocaldarius stunted biofilm and planktonic growth and increased motility structures (archaella). Conversely, planktonic cells were devoid of archaella in the ΔvapC14 cognate toxin mutant. VapB14 is highly conserved at both the nucleotide and amino acid levels across the Sulfolobales, extremely unusual for type II antitoxins, which are typically acquired through horizontal gene transfer. Furthermore, homologs of VapB14 are found across the Crenarchaeota, in some Euryarchaeota, and even bacteria. S. acidocaldarius vapB14 and its homolog in the thermoacidophile Metallosphaera sedula (Msed_0871) were both upregulated in biofilm cells, supporting the role of the antitoxin in biofilm regulation. In several Sulfolobales species, including M. sedula, homologs of vapB14 and vapC14 are not colocalized. Strikingly, Sulfuracidifex tepidarius has an unpaired VapB14 homolog and lacks a cognate VapC14, illustrating the toxin-independent conservation of the VapB14 antitoxin. The findings here suggest that a stand-alone VapB-type antitoxin was the product of selective evolutionary pressure to influence biofilm formation in these archaea, a vital microbial community behavior. IMPORTANCE Biofilms allow microbes to resist a multitude of stresses and stay proximate to vital nutrients. The mechanisms of entering and leaving a biofilm are highly regulated to ensure microbial survival, but are not yet well described in archaea. Here, a VapBC type II toxin-antitoxin system in the thermoacidophilic archaeon Sulfolobus acidocaldarius was shown to control biofilm dispersal through a multifaceted regulation of the archaeal motility structure, the archaellum. The VapC14 toxin degrades an RNA that causes an increase in archaella and swimming. The VapB14 antitoxin decreases archaella and biofilm dispersal by binding the VapC14 toxin and neutralizing its activity, while also repressing the archaellum genes. VapB14-like antitoxins are highly conserved across the Sulfolobales and respond similarly to biofilm growth. In fact, VapB14-like antitoxins are also found in other archaea, and even in bacteria, indicating an evolutionary pressure to maintain this protein and its role in biofilm formation.

Identification and Characterization of HEPN-MNT Type II TA System from Methanothermobacter thermautotrophicus ΔH

Toxin-antitoxin (TA) systems are widespread in bacteria and archaea plasmids and genomes to regulate DNA replication, gene transcription, or protein translation. Higher eukaryotic and prokaryotic nucleotide-binding (HEPN) and minimal nucleotidyltransferase (MNT) domains are prevalent in prokaryotic genomes and constitute TA pairs. However, three gene pairs (MTH304/305, 408/409, and 463/464) of Methanothermobacter thermautotropicus ΔH HEPN-MNT family have not been studied as TA systems. Among these candidates, our study characterizes the MTH463/MTH464 TA system. MTH463 expression inhibited Escherichia coli growth, whereas MTH464 did not and blocked MTH463 instead. Using site-directed MTH463 mutagenesis, we determined that amino acids R99G, H104A, and Y106A from the R[ɸX]4-6H motif are involved with MTH463 cell toxicity. Furthermore, we established that purified MTH463 could degrade MS2 phage RNA, whereas purified MTH464 neutralized MTH463 activity in vitro. Our results indicate that the endonuclease toxin MTH463 (encoding a HEPN domain) and its cognate antitoxin MTH464 (encoding the MNT domain) may act as a type II TA system in M. thermautotropicus ΔH. This study provides initial and essential information studying TA system functions, primarily archaea HEPN-MNT family.

A broadly distributed predicted helicase/nuclease confers phage resistance via abortive infection

There is strong selection for the evolution of systems that protect bacterial populations from viral attack. We report a single phage defense protein, Hna, that provides protection against diverse phages in Sinorhizobium meliloti, a nitrogen-fixing alpha-proteobacterium. Homologs of Hna are distributed widely across bacterial lineages, and a homologous protein from Escherichia coli also confers phage defense. Hna contains superfamily II helicase motifs at its N terminus and a nuclease motif at its C terminus, with mutagenesis of these motifs inactivating viral defense. Hna variably impacts phage DNA replication but consistently triggers an abortive infection response in which infected cells carrying the system die but do not release phage progeny. A similar host cell response is triggered in cells containing Hna upon expression of a phage-encoded single-stranded DNA binding protein (SSB), independent of phage infection. Thus, we conclude that Hna limits phage spread by initiating abortive infection in response to a phage protein.

A functional selection reveals previously undetected anti-phage defence systems in the E. coli pangenome

The ancient, ongoing coevolutionary battle between bacteria and their viruses, bacteriophages, has given rise to sophisticated immune systems including restriction-modification and CRISPR-Cas. Many additional anti-phage systems have been identified using computational approaches based on genomic co-location within defence islands, but these screens may not be exhaustive. Here we developed an experimental selection scheme agnostic to genomic context to identify defence systems in 71 diverse E. coli strains. Our results unveil 21 conserved defence systems, none of which were previously detected as enriched in defence islands. Additionally, our work indicates that intact prophages and mobile genetic elements are primary reservoirs and distributors of defence systems in E. coli, with defence systems typically carried in specific locations or hotspots. These hotspots encode dozens of additional uncharacterized defence system candidates. Our findings reveal an extended landscape of antiviral immunity in E. coli and provide an approach for mapping defence systems in other species.

Analysis of lineage-specific protein family variability in prokaryotes combined with evolutionary reconstructions

Background

Evolutionary rate is a key characteristic of gene families that is linked to the functional importance of the respective genes as well as specific biological functions of the proteins they encode. Accurate estimation of evolutionary rates is a challenging task that requires precise phylogenetic analysis. Here we present an easy to estimate protein family level measure of sequence variability based on alignment column homogeneity in multiple alignments of protein sequences from Clade-Specific Clusters of Orthologous Genes (csCOGs).

Results

We report genome-wide estimates of variability for 8 diverse groups of bacteria and archaea and investigate the connection between variability and various genomic and biological features. The variability estimates are based on homogeneity distributions across amino acid sequence alignments and can be obtained for multiple groups of genomes at minimal computational expense. About half of the variance in variability values can be explained by the analyzed features, with the greatest contribution coming from the extent of gene paralogy in the given csCOG. The correlation between variability and paralogy appears to originate, primarily, not from gene duplication, but from acquisition of distant paralogs and xenologs, introducing sequence variants that are more divergent than those that could have evolved in situ during the lifetime of the given group of organisms. Both high-variability and low-variability csCOGs were identified in all functional categories, but as expected, proteins encoded by integrated mobile elements as well as proteins involved in defense functions and cell motility are, on average, more variable than proteins with housekeeping functions. Additionally, using linear discriminant analysis, we found that variability and fraction of genomes carrying a given gene are the two variables that provide the best prediction of gene essentiality as compared to the results of transposon mutagenesis in Sulfolobus islandicus.

Conclusions

Variability, a measure of sequence diversity within an alignment relative to the overall diversity within a group of organisms, offers a convenient proxy for evolutionary rate estimates and is informative with respect to prediction of functional properties of proteins. In particular, variability is a strong predictor of gene essentiality for the respective organisms and indicative of sub- or neofunctionalization of paralogs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13062-022-00337-7.

Borreliella burgdorferi Antimicrobial-Tolerant Persistence in Lyme Disease and Posttreatment Lyme Disease Syndromes

The annual incidence of Lyme disease, caused by tick-transmitted Borreliella burgdorferi, is estimated to be at least 476,000 cases in the United States and many more worldwide. Ten to 20% of antimicrobial-treated Lyme disease patients display posttreatment Lyme disease syndrome (PTLDS), a clinical complication whose etiology and pathogenesis remain uncertain. Autoimmunity, cross-reactivity, molecular mimicry, coinfections, and borrelial tolerance to antimicrobials/persistence have been hypothesized and studied as potential causes of PTLDS. Studies of borrelial tolerance/persistence in vitro in response to antimicrobials and experimental studies in mice and nonhuman primates, taken together with clinical reports, have revealed that B. burgdorferi becomes tolerant to antimicrobials and may sometimes persist in animals and humans after the currently recommended antimicrobial treatment. Moreover, B. burgdorferi is pleomorphic and can generate viable-but-nonculturable bacteria, states also involved in antimicrobial tolerance. The multiple regulatory pathways and structural genes involved in mediating this tolerance to antimicrobials and environmental stressors by persistence might include the stringent (rel and dksA) and host adaptation (rpoS) responses, sugar metabolism (glpD), and polypeptide transporters (opp). Application of this recently reported knowledge to clinical studies can be expected to clarify the potential role of bacterial antibacterial tolerance/persistence in Lyme disease and PTLDS.

Comparative Analysis of Diverse Acetyltransferase-Type Toxin-Antitoxin Loci in Klebsiella pneumoniae

Toxin-antitoxin (TA) modules containing a Gcn5-related N-acetyltransferase (GNAT) toxin domain regulate bacterial physiology under adverse environmental stresses. Multiple GNAT-ribbon-helix-helix domain (RHH) TA loci have been identified in single bacterial genomes. However, their diversity and interactions are still obscure. Our previous analysis showed that the GNAT toxin of Klebsiella pneumoniae, KacT, introduces antibiotic tolerance and the toxicity of GNAT is neutralized by KacA, an RHH antitoxin. We here present a phylogenetic analysis of GNAT toxins of more than 1,000 GNAT-RHH pairs detected in completely sequenced K. pneumoniae genomes, revealing that the GNAT toxins are diverse and grouped into four distinct clades. Overexpression of GNAT toxins representative of each of the four clades halts the cell growth of K. pneumoniae, while the coexpression of the cognate RHH antitoxin neutralizes GNAT toxicity. We also identify point mutations that inactivate the GNAT toxins. Moreover, we observe a cross-interaction between GNAT-RHH pairs encoded by different replicons, where a chromosomal toxin (KacT2) can be neutralized by its cognate RHH antitoxin (KacA2 on a chromosome) and another antitoxin (KacA3 on a plasmid). Finally, statistical analysis of the distribution of GNAT-RHH loci in K. pneumoniae strains shows pronounced deviation from random distribution within the same clades. Moreover, we also obtain statistically significant correlations between different clades, which we discuss in terms of the experimental results.

IMPORTANCE Elucidating the roles of multifaceted GNAT-RHH TA loci is essential for understanding how these TAs interact among themselves. Recently, the reaction mechanisms and structures of several GNAT-RHH pairs have been reported. While bacterial strains can carry multiple GNAT-RHH loci with diverse origins, studies on the possible cross-interactions of these TA pairs are still limited. Here, we find that 1,000 predicted GNAT toxins of K. pneumoniae can be grouped into four distinct clades. The distributions of TA loci within these clades in K. pneumoniae strains are highly nonrandom, with the presence of a single locus of each clade per strain being highly overrepresented. Moreover, the toxicity of a GNAT toxin encoded by a chromosome was alleviated by a noncognate RHH antitoxin on a plasmid. These results might yield a profound understanding of the widespread GNAT-RHH TA pairs and the cross-interactions between noncognate TA pairs located on different replicons.

GNAT toxins evolve toward narrow tRNA target specificities

Type II toxin–antitoxin (TA) systems are two-gene modules widely distributed among prokaryotes. GNAT toxins associated with the DUF1778 antitoxins represent a large family of type II TAs. GNAT toxins inhibit cell growth by disrupting translation via acetylation of aminoacyl-tRNAs. In this work, we explored the evolutionary trajectory of GNAT toxins. Using LC/MS detection of acetylated aminoacyl-tRNAs combined with ribosome profiling, we systematically investigated the in vivo substrate specificity of an array of diverse GNAT toxins. Our functional data show that the majority of GNAT toxins are specific to Gly-tRNA isoacceptors. However, the phylogenetic analysis shows that the ancestor of GNAT toxins was likely a relaxed specificity enzyme capable of acetylating multiple elongator tRNAs. Together, our data provide a remarkable snapshot of the evolution of substrate specificity.

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.

Toxin-Antitoxin Systems as Phage Defense Elements

Toxin-antitoxin (TA) systems are ubiquitous genetic elements in bacteria that consist of a growth-inhibiting toxin and its cognate antitoxin. These systems are prevalent in bacterial chromosomes, plasmids, and phage genomes, but individual systems are not highly conserved, even among closely related strains. The biological functions of TA systems have been controversial and enigmatic, although a handful of these systems have been shown to defend bacteria against their viral predators, bacteriophages. Additionally, their patterns of conservation-ubiquitous, but rapidly acquired and lost from genomes-as well as the co-occurrence of some TA systems with known phage defense elements are suggestive of a broader role in mediating phage defense. Here, we review the existing evidence for phage defense mediated by TA systems, highlighting how toxins are activated by phage infection and how toxins disrupt phage replication. We also discuss phage-encoded systems that counteract TA systems, underscoring the ongoing coevolutionary battle between bacteria and phage. We anticipate that TA systems will continue to emerge as central players in the innate immunity of bacteria against phage. 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.

Recruitment of Mobile Genetic Elements for Diverse Cellular Functions in Prokaryotes

Prokaryotic genomes are replete with mobile genetic elements (MGE) that span a continuum of replication autonomy. On numerous occasions during microbial evolution, diverse MGE lose their autonomy altogether but, rather than being quickly purged from the host genome, assume a new function that benefits the host, rendering the immobilized MGE subject to purifying selection, and resulting in its vertical inheritance. This mini-review highlights the diversity of the repurposed (exapted) MGE as well as the plethora of cellular functions that they perform. The principal contribution of the exaptation of MGE and their components is to the prokaryotic functional systems involved in biological conflicts, and in particular, defense against viruses and other MGE. This evolutionary entanglement between MGE and defense systems appears to stem both from mechanistic similarities and from similar evolutionary predicaments whereby both MGEs and defense systems tend to incur fitness costs to the hosts and thereby evolve mechanisms for survival including horizontal mobility, causing host addiction, and exaptation for functions beneficial to the host. The examples discussed demonstrate that the identity of an MGE, overall mobility and relationship with the host cell (mutualistic, symbiotic, commensal, or parasitic) are all factors that affect exaptation.

Comprehensive and functional analyses reveal the genomic diversity and potential toxicity of Microcystis

Microcystis is a cyanobacteria that is widely distributed across the world. It has attracted great attention because it produces the hepatotoxin microcystin (MC) that can inhibit eukaryotic protein phosphatases and pose a great risk to animal and human health. Due to the high diversity of morphospecies and genomes, it is still difficult to classify Microcystis species. In this study, we investigated the pangenome of 23 Microcystis strains to detect the genetic diversity and evolutionary dynamics. Microcystis revealed an open pangenome containing 22,009 gene families and exhibited different functional constraints. The core-genome phylogenetic analysis accurately differentiated the toxic and nontoxic strains and could be used as a taxonomic standard at the genetic level. We also investigated the functions of HGT events, of which were mostly conferred from cyanobacteria and closely related species. In order to detect the potential toxicity of Microcystis, we searched and characterized MC biosynthetic gene clusters and other secondary metabolite gene clusters. Our work provides insights into the genetic diversity, evolutionary dynamics, and potential toxicity of Microcystis, which could benefit the species classification and development of new methods for drinking water quality control and management of bloom formation in the future.

Genome binning of viral entities from bulk metagenomics data

Despite the accelerating number of uncultivated virus sequences discovered in metagenomics and their apparent importance for health and disease, the human gut virome and its interactions with bacteria in the gastrointestinal tract are not well understood. This is partly due to a paucity of whole-virome datasets and limitations in current approaches for identifying viral sequences in metagenomics data. Here, combining a deep-learning based metagenomics binning algorithm with paired metagenome and metavirome datasets, we develop Phages from Metagenomics Binning (PHAMB), an approach that allows the binning of thousands of viral genomes directly from bulk metagenomics data, while simultaneously enabling clustering of viral genomes into accurate taxonomic viral populations. When applied on the Human Microbiome Project 2 (HMP2) dataset, PHAMB recovered 6,077 high-quality genomes from 1,024 viral populations, and identified viral-microbial host interactions. PHAMB can be advantageously applied to existing and future metagenomes to illuminate viral ecological dynamics with other microbiome constituents.

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.

Streptococcal bacterial components in cancer therapy

The incidence rate of cancer is steadily increasing all around the world, and there is an urgent need to develop novel and more effective treatment strategies. Recently, bacterial therapy has been investigated as a new approach to target cancer, and is becoming a serious option. Streptococcus strains are among the most common and well-studied virulent bacteria that cause a variety of human infections. Everyone has experienced a sore throat during their lifetime, or has been asymptomatically colonized by streptococci. The ability of Streptococcus bacteria to fight cancer was discovered more than 100 years ago, and over the years has undergone clinical trials, but the mechanism is not yet completely understood. Recently, several animal models and human clinical trials have been reported. Streptococcal strains can have an intrinsic anti-tumor activity, or can activate the host immune system to fight the tumor. Bacteria can selectively accumulate and proliferate in the hypoxic regions of solid tumors. Moreover, the bacteria can be genetically engineered to secrete toxins or enzymes that can specifically attack the tumors.

Biology and evolution of bacterial toxin-antitoxin systems

Toxin-antitoxin systems are widespread in bacterial genomes. They are usually composed of two elements: a toxin that inhibits an essential cellular process and an antitoxin that counteracts its cognate toxin. In the past decade, a number of new toxin-antitoxin systems have been described, bringing new growth inhibition mechanisms to light as well as novel modes of antitoxicity. However, recent advances in the field profoundly questioned the role of these systems in bacterial physiology, stress response and antimicrobial persistence. This shifted the paradigm of the functions of toxin-antitoxin systems to roles related to interactions between hosts and their mobile genetic elements, such as viral defence or plasmid stability. In this Review, we summarize the recent progress in understanding the biology and evolution of these small genetic elements, and discuss how genomic conflicts could shape the diversification of toxin-antitoxin systems.

Cargo Genes of Tn7-Like Transposons Comprise an Enormous Diversity of Defense Systems, Mobile Genetic Elements, and Antibiotic Resistance Genes

Transposition is a major mechanism of horizontal gene mobility in prokaryotes. However, exploration of the genes mobilized by transposons (cargo) is hampered by the difficulty in delineating integrated transposons from their surrounding genetic context. Here, we present a computational approach that allowed us to identify the boundaries of 6,549 Tn7-like transposons. We found that 96% of these transposons carry at least one cargo gene. Delineation of distinct communities in a gene-sharing network demonstrates how transposons function as a conduit of genes between phylogenetically distant hosts. Comparative analysis of the cargo genes reveals significant enrichment of mobile genetic elements (MGEs) nested within Tn7-like transposons, such as insertion sequences and toxin-antitoxin modules, and of genes involved in recombination, anti-MGE defense, and antibiotic resistance. More unexpectedly, cargo also includes genes encoding central carbon metabolism enzymes. Twenty-two Tn7-like transposons carry both an anti-MGE defense system and antibiotic resistance genes, illustrating how bacteria can overcome these combined pressures upon acquisition of a single transposon. This work substantially expands the distribution of Tn7-like transposons, defines their evolutionary relationships, and provides a large-scale functional classification of prokaryotic genes mobilized by transposition.

Minigene as a Novel Regulatory Element in Toxin-Antitoxin Systems

The axe-txe type II toxin-antitoxin (TA) system is characterized by a complex and multilayered mode of gene expression regulation. Precise and tight control of this process is crucial to keep the toxin in an appropriate balance with the cognate antitoxin until its activation is needed for the cell. In this report, we provide evidence that a minigene encoded within the axe-txe operon influences translation of the Txe toxin. This is the first example to date of such a regulatory mechanism identified in the TA modules. Here, in a series of genetic studies, we employed translational reporter gene fusions to establish the molecular basis of this phenomenon. Our results show that translation of the two-codon mini-ORF displays an in cis mode of action, and positively affects the expression of txe, possibly by increasing its mRNA stability through protection from an endonuclease attack. Moreover, we established that the reading frame in which the two cistrons are encoded, as well as the distance between them, are critical parameters that affect the level of such regulation. In addition, by searching for two-codon ORFs we found sequences of several potential minigenes in the leader sequences of several other toxins belonging to the type II TA family. These findings suggest that this type of gene regulation may not only apply for the axe-txe cassette, but could be more widespread among other TA systems.

Conserved Amino Acid Moieties of Candidatus Desulforudis audaxviator MazF Determine Ribonuclease Activity and Specificity

The toxin-antitoxin (TA) system, inherent to various prokaryotes, plays a critical role in survival and adaptation to diverse environmental stresses. The toxin MazF, belonging to the type II TA system, functions as a sequence-specific ribonuclease that recognizes 3 to 7 bases. In recent studies, crystallographic analysis of MazFs from several species have suggested the presence of amino acid sites important for MazF substrate RNA binding and for its catalytic activity. Herein, we characterized MazF obtained from Candidatus Desulforudis audaxviator (MazF-Da) and identified the amino acid residues necessary for its catalytic function. MazF-Da, expressed using a cell-free protein synthesis system, is a six-base-recognition-specific ribonuclease that preferentially cleaves UACAAA sequences and weakly cleaves UACGAA and UACUAA sequences. We found that MazF-Da exhibited the highest activity at around 60°C. Analysis using mutants with a single mutation at an amino acid residue site that is well conserved across various MazF toxins showed that G18, E20, R25, and P26 were important for the ribonuclease activity of MazF-Da. The recognition sequence of the N36A mutant differed from that of the wild type. This mutant cleaved UACAAG sequences in addition to UACAAA sequences, but did not cleave UACGAA or UACUAA sequences, suggesting that Asn36 affects the loosening and narrowing of MazF-Da cleavage sequence recognition. Our study posits UACAAA as the recognition sequence of MazF-Da and provides insight into the amino acid sites that are key to its unique enzymatic properties.

A toxin-antitoxin system confers stability to the IncP-7 plasmid pCAR1

Toxin-antitoxin (TA) systems were initially discovered as plasmid addiction systems. Previously, our studies implied that the high stability of the IncP-7 plasmid pCAR1 in different Pseudomonas spp. hosts was due to the presence of a TA system on the plasmid. Bioinformatics approaches suggested that ORF174 and ORF175 could constitute a type II TA system, a member of the RES-Xre family, and that these two open reading frames (ORFs) constitute a single operon. As expected, the ORF175 product is a toxin, which decreases the viability of the host, P. resinovorans, while the ORF174 product functions as an antitoxin that counteracts the effect of ORF175 on cell growth. Based on these findings, we renamed ORF174 and ORF175 as prcA (antitoxin gene) and prcT (toxin gene), respectively. The prcA and prcT genes were cloned into the unstable plasmid vector pSEVA644. The recombinant vector was stably maintained in P. resinovorans and Escherichia coli cells under nonselective conditions following 6 days of daily subculturing. The empty vector (without the prcA and prcT genes) could not be maintained, which suggested that the PrcA/T system can be used as a tool to improve the stability of otherwise unstable plasmids in P. resinovorans and E. coli strains.

Genic Selection Within Prokaryotic Pangenomes

Understanding the evolutionary forces shaping prokaryotic pangenome structure is a major goal of microbial evolution research. Recent work has highlighted that a substantial proportion of accessory genes appear to confer niche-specific adaptations. This work has primarily focused on selection acting at the level of individual cells. Herein, we discuss a lower level of selection that also contributes to pangenome variation: genic selection. This refers to cases where genetic elements, rather than individual cells, are the entities under selection. The clearest examples of this form of selection are selfish mobile genetic elements, which are those that have either a neutral or a deleterious effect on host fitness. We review the major classes of these and other mobile elements and discuss the characteristic features of such elements that could be under genic selection. We also discuss how genetic elements that are beneficial to hosts can also be under genic selection, a scenario that may be more prevalent but not widely appreciated, because disentangling the effects of selection at different levels (i.e., organisms vs. genes) is challenging. Nonetheless, an appreciation for the potential action and implications of genic selection is important to better understand the evolution of prokaryotic pangenomes.

Identification of Type II Toxin-Antitoxin Loci in Levilactobacillus brevis

Levilactobacillus brevis are present in various environments, such as beer, fermented foods, silage, and animal host. Like other lactic acid bacteria, L. brevis might adopt the viable but nonculturable (VBNC) state under unfavorable conditions. The toxin-antitoxin (TA) system, known to regulate cell growth in response to environmental stresses, is found to control the dynamic of the VBNC state. Here, we investigate the type II TA locus prevalence and compare the TA diversity in L. brevis genomes. Using the TAfinder software, we identified a total of 273 putative type II TA loci in 110 replicons of 21 completely sequenced genomes. Genome size does not appear to correlate with the amount of putative type II TA in L. brevis. Besides, type II TA loci are distributed differently among the chromosomes and plasmids. The most prevalent toxin domain is MazF-like in the chromosomes, and RelE/RelE-like in the plasmids; while for antitoxin, Xre-like and Phd-like domains are the most common in the chromosomes and plasmids, respectively. We also observed a unique GNAT-like/ArsR-like TA pair that presents only in the L. brevis chromosome. Detection of 273 putative type II TA loci in 21 complete genomes of Levilactobacillus brevis.

Evaluation of Putative Toxin-antitoxins Systems in Clinical Brucella melitensis in Iran

BACKGROUND

Toxin-antitoxin systems (TAs) are two-component elements, which are extensive in the bacterial genome and have a regulatory role in many cellular activities including, growth arrest, survival, biofilm formation, and bacterial persistence.

OBJECTIVE

TAs have not well studied in Brucella spp.

METHODS

We evaluated the presence of different toxin-antitoxin systems, including relE- rhhlike, Fic- Phd, Cog- Rhh, and cogT- cogAT in 40 clinical Brucella melitensis isolates using PCRbased sequencing assay.

RESULTS

Our results showed the high presence of relE-rhh-like, Fic-Phd, Cog-rhh, and cogTcogAT s TAs genes in B. melitensis isolates that were 96.25%, 92.5%, 96.25%, and 95%, respectively.

CONCLUSION

A high presence of TAs genes in clinical B. melitensis isolates revealed that the TA system could be an antibacterial target in B. melitensis but more investigation is necessitated to elucidate the exact roles of these genes.

An Auto-Regulating Type II Toxin-Antitoxin System Modulates Drug Resistance and Virulence in Streptococcus suis

Toxin-antitoxin (TA) systems are ubiquitous genetic elements that play an essential role in multidrug tolerance and virulence of bacteria. So far, little is known about the TA systems in Streptococcus suis. In this study, the Xress-MNTss TA system, composed of the MNTss toxin in the periplasmic space and its interacting Xress antitoxin, was identified in S. suis. β-galactosidase activity and electrophoretic mobility shift assay (EMSA) revealed that Xress and the Xress-MNTss complex could bind directly to the Xress-MNTss promoter as well as downregulate streptomycin adenylyltransferase ZY05719_RS04610. Interestingly, the Xress deletion mutant was less pathogenic in vivo following a challenge in mice. Transmission electron microscopy and adhesion assays pointed to a significantly thinner capsule but greater biofilm-formation capacity in ΔXress than in the wild-type strain. These results indicate that Xress-MNTss, a new type II TA system, plays an important role in antibiotic resistance and pathogenicity in S. suis.

Comparative genome analysis of Salmonella enterica serovar Gallinarum biovars Pullorum and Gallinarum decodes strain specific genes

Salmonella enterica serovar Gallinarum biovar Pullorum (bvP) and biovar Gallinarum (bvG) are the etiological agents of pullorum disease (PD) and fowl typhoid (FT) respectively, which cause huge economic losses to poultry industry especially in developing countries including India. Vaccination and biosecurity measures are currently being employed to control and reduce the S. Gallinarum infections. High endemicity, poor implementation of hygiene and lack of effective vaccines pose challenges in prevention and control of disease in intensively maintained poultry flocks. Comparative genome analysis unravels similarities and dissimilarities thus facilitating identification of genomic features that aids in pathogenesis, niche adaptation and in tracing of evolutionary history. The present investigation was carried out to assess the genotypic differences amongst S.enterica serovar Gallinarum strains including Indian strain S. Gallinarum Sal40 VTCCBAA614. The comparative genome analysis revealed an open pan-genome consisting of 5091 coding sequence (CDS) with 3270 CDS belonging to core-genome, 1254 CDS to dispensable genome and strain specific genes i.e. singletons ranging from 3 to 102 amongst the analyzed strains. Moreover, the investigated strains exhibited diversity in genomic features such as virulence factors, genomic islands, prophage regions, toxin-antitoxin cassettes, and acquired antimicrobial resistance genes. Core genome identified in the study can give important leads in the direction of design of rapid and reliable diagnostics, and vaccine design for effective infection control as well as eradication. Additionally, the identified genetic differences among the S. enterica serovar Gallinarum strains could be used for bacterial typing, structure based inhibitor development by future experimental investigations on the data generated.

Alternative dimerization is required for activity and inhibition of the HEPN ribonuclease RnlA

The rnlAB toxin-antitoxin operon from Escherichia coli functions as an anti-phage defense system. RnlA was identified as a member of the HEPN (Higher Eukaryotes and Prokaryotes Nucleotide-binding domain) superfamily of ribonucleases. The activity of the toxin RnlA requires tight regulation by the antitoxin RnlB, the mechanism of which remains unknown. Here we show that RnlA exists in an equilibrium between two different homodimer states: an inactive resting state and an active canonical HEPN dimer. Mutants interfering with the transition between states show that canonical HEPN dimerization via the highly conserved RX4-6H motif is required for activity. The antitoxin RnlB binds the canonical HEPN dimer conformation, inhibiting RnlA by blocking access to its active site. Single-alanine substitutions mutants of the highly conserved R255, E258, R318 and H323 show that these residues are involved in catalysis and substrate binding and locate the catalytic site near the dimer interface of the canonical HEPN dimer rather than in a groove located between the HEPN domain and the preceding TBP-like domain. Overall, these findings elucidate the structural basis of the activity and inhibition of RnlA and highlight the crucial role of conformational heterogeneity in protein function.

The extracellular contractile injection system is enriched in environmental microbes and associates with numerous toxins

The extracellular Contractile Injection System (eCIS) is a toxin-delivery particle that evolved from a bacteriophage tail. Four eCISs have previously been shown to mediate interactions between bacteria and their invertebrate hosts. Here, we identify eCIS loci in 1,249 bacterial and archaeal genomes and reveal an enrichment of these loci in environmental microbes and their apparent absence from mammalian pathogens. We show that 13 eCIS-associated toxin genes from diverse microbes can inhibit the growth of bacteria and/or yeast. We identify immunity genes that protect bacteria from self-intoxication, further supporting an antibacterial role for some eCISs. We also identify previously undescribed eCIS core genes, including a conserved eCIS transcriptional regulator. Finally, we present our data through an extensive eCIS repository, termed eCIStem. Our findings support eCIS as a toxin-delivery system that is widespread among environmental prokaryotes and likely mediates antagonistic interactions with eukaryotes and other prokaryotes.

Evolutionary history of Caulobacter toxin-antitoxin systems

Toxin-antitoxin (TA) systems have been studied in many bacterial genera, but a clear understanding of the evolutionary trajectory of TA operons has not emerged. To address this issue, I identified 42 distinct TA operons in three genomes that represent the three branches of the Caulobacter phylogenetic tree. The location of each operon was then examined to determine if the operon was present in eight additional Caulobacter genomes. Most of the 42 TA operons were present at the same chromosomal location in genomes that represent at least two different branches of the Caulobacter phylogenetic tree. This result indicates that the chromosomal location of TA operons is conserved over evolutionary time scales. One the other hand, there were 177 instances where a TA operon was not present at an expected chromosomal location and four instances where only the antitoxin gene was present. Thus, the variable number of TA operons found in each genome appears to be due primarily to the loss of TA operons, and the addition of new TA operons to a genome was relatively rare. An additional feature of the TA operons was that they seemed to accumulate mutations faster than the adjacent genes.

Elucidating the Regulatory Elements for Transcription Termination and Posttranscriptional Processing in the Streptomyces clavuligerus Genome

Identification of transcriptional regulatory elements in the GC-rich Streptomyces genome is essential for the production of novel biochemicals from secondary metabolite biosynthetic gene clusters (smBGCs). Despite many efforts to understand the regulation of transcription initiation in smBGCs, information on the regulation of transcription termination and posttranscriptional processing remains scarce. In this study, we identified the transcriptional regulatory elements in β-lactam antibiotic-producing Streptomyces clavuligerus ATCC 27064 by determining a total of 1,427 transcript 3'-end positions (TEPs) using the term-seq method. Termination of transcription was governed by three classes of TEPs, of which each displayed unique sequence features. The data integration with transcription start sites and transcriptome data generated 1,648 transcription units (TUs) and 610 transcription unit clusters (TUCs). TU architecture showed that the transcript abundance in TU isoforms of a TUC was potentially affected by the sequence context of their TEPs, suggesting that the regulatory elements of TEPs could control the transcription level in additional layers. We also identified TU features of a xenobiotic response element (XRE) family regulator and DUF397 domain-containing protein, particularly showing the abundance of bidirectional TEPs. Finally, we found that 189 noncoding TUs contained potential cis- and trans-regulatory elements that played a major role in regulating the 5' and 3' UTR. These findings highlight the role of transcriptional regulatory elements in transcription termination and posttranscriptional processing in Streptomyces sp.IMPORTANCE Streptomyces sp. is a great source of bioactive secondary metabolites, including antibiotics, antifungal agents, antiparasitic agents, immunosuppressant compounds, and other drugs. Secondary metabolites are synthesized via multistep conversions of the precursor molecules from primary metabolism, governed by multicomplex enzymes from secondary metabolite biosynthetic gene clusters. As their production is closely related with the growth phase and dynamic cellular status in response to various intra- and extracellular signals, complex regulatory systems tightly control the gene expressions related to secondary metabolism. In this study, we determined genome-wide transcript 3'-end positions and transcription units in the β-lactam antibiotic producer Streptomyces clavuligerus ATCC 27064 to elucidate the transcriptional regulatory elements in transcription termination and posttranscriptional processing by integration of multiomics data. These unique features, such as transcript 3'-end sequence, potential riboregulators, and potential 3'-untranslated region (UTR) cis-regulatory elements, can be potentially used to design engineering tools that can regulate the transcript abundance of genes for enhancing secondary metabolite production.

Functional characterization of Kid-Kis and MazF-MazE in Sf9 cells and Mythimna separata embryos

Toxin-antitoxin (TA) systems are comprised of a toxin and its antidote antitoxin and are widely present in bacterial and in eukaryotic systems. However, no work regarding TA systems has been reported in insects. We characterized the Kid-Kis and MazF-MazE TA systems in Spodoptera frugiperda cells and Mythimna separata embryos and observed that the Kid and MazF toxins were highly toxic. In Sf9 cells transfected with Kid plasmid and MazF alone, the apoptosis rate was 24.37% and 29.47%, respectively. Whereas the toxicity of their cognate antitoxins were limited. Both apoptosis and necrosis were induced by the two toxins. Both the Kis and MazE antitoxins partly neutralized toxicity in a dose-dependent manner, with MazE accomplishing almost full neutralization at a 1:4 toxin:antitoxin ratio, the cell survival rate was 81% and 97%, respectively. Our results indicate that MazF-MazE is a good candidate module for application in insects, such as in developing new sterile insect technique (SIT).

Coexistence of SOS-Dependent and SOS-Independent Regulation of DNA Repair Genes in Radiation-Resistant Deinococcus Bacteria

Deinococcus bacteria are extremely resistant to radiation and able to repair a shattered genome in an essentially error-free manner after exposure to high doses of radiation or prolonged desiccation. An efficient, SOS-independent response mechanism to induce various DNA repair genes such as recA is essential for radiation resistance. This pathway, called radiation/desiccation response, is controlled by metallopeptidase IrrE and repressor DdrO that are highly conserved in Deinococcus. Among various Deinococcus species, Deinococcus radiodurans has been studied most extensively. Its genome encodes classical DNA repair proteins for error-free repair but no error-prone translesion DNA polymerases, which may suggest that absence of mutagenic lesion bypass is crucial for error-free repair of massive DNA damage. However, many other radiation-resistant Deinococcus species do possess translesion polymerases, and radiation-induced mutagenesis has been demonstrated. At least dozens of Deinococcus species contain a mutagenesis cassette, and some even two cassettes, encoding error-prone translesion polymerase DnaE2 and two other proteins, ImuY and ImuB-C, that are probable accessory factors required for DnaE2 activity. Expression of this mutagenesis cassette is under control of the SOS regulators RecA and LexA. In this paper, we review both the RecA/LexA-controlled mutagenesis and the IrrE/DdrO-controlled radiation/desiccation response in Deinococcus.

Redox signaling through zinc activates the radiation response in Deinococcus bacteria

Deinococcus bacteria are extremely resistant to radiation and other DNA damage- and oxidative stress-generating conditions. An efficient SOS-independent response mechanism inducing expression of several DNA repair genes is essential for this resistance, and is controlled by metalloprotease IrrE that cleaves and inactivates transcriptional repressor DdrO. Here, we identify the molecular signaling mechanism that triggers DdrO cleavage. We show that reactive oxygen species (ROS) stimulate the zinc-dependent metalloprotease activity of IrrE in Deinococcus. Sudden exposure of Deinococcus to zinc excess also rapidly induces DdrO cleavage, but is not accompanied by ROS production and DNA damage. Further, oxidative treatment leads to an increase of intracellular free zinc, indicating that IrrE activity is very likely stimulated directly by elevated levels of available zinc ions. We conclude that radiation and oxidative stress induce changes in redox homeostasis that result in IrrE activation by zinc in Deinococcus. We propose that a part of the zinc pool coordinated with cysteine thiolates is released due to their oxidation. Predicted regulation systems involving IrrE- and DdrO-like proteins are present in many bacteria, including pathogens, suggesting that such a redox signaling pathway including zinc as a second messenger is widespread and participates in various stress responses.

The Disordered C-Terminus of the Chaperone DnaK Increases the Competitive Fitness of Pseudomonas putida and Facilitates the Toxicity of GraT

Chaperone proteins are crucial for proper protein folding and quality control, especially when cells encounter stress caused by non-optimal temperatures. DnaK is one of such essential chaperones in bacteria. Although DnaK has been well characterized, the function of its intrinsically disordered C-terminus has remained enigmatic as the deletion of this region has been shown to either enhance or reduce its protein folding ability. We have shown previously that DnaK interacts with toxin GraT of the GraTA toxin-antitoxin system in Pseudomonas putida. Interestingly, the C-terminal truncation of DnaK was shown to alleviate GraT-caused growth defects. Here, we aim to clarify the importance of DnaK in GraT activity. We show that DnaK increases GraT toxicity, and particularly important is the negatively charged motif in the DnaK C-terminus. Given that GraT has an intrinsically disordered N-terminus, the assistance of DnaK is probably needed for re-modelling the toxin structure. We also demonstrate that the DnaK C-terminal negatively charged motif contributes to the competitive fitness of P. putida at both high and optimal growth temperatures. Thus, our data suggest that the disordered C-terminal end of DnaK enhances the chaperone functionality.