BackgroundBacterial toxin-antitoxin (TA) systems are formed by potent regulatory or suicide factors (toxins) and their short-lived inhibitors (antitoxins). Antitoxins are DNA-binding proteins and auto-repress transcription of TA operons. Transcription of multiple TA operons is activated in temporarily non-growing persister cells that can resist killing by antibiotics. Consequently, the antitoxin levels of persisters must have been dropped and toxins are released of inhibition.ResultsHere, we describe transcriptional cross-activation between different TA systems of Escherichia coli. We find that the chromosomal relBEF operon is activated in response to production of the toxins MazF, MqsR, HicA, and HipA. Expression of the RelE toxin in turn induces transcription of several TA operons. We show that induction of mazEF during amino acid starvation depends on relBE and does not occur in a relBEF deletion mutant. Induction of TA operons has been previously shown to depend on Lon protease which is activated by polyphospate accumulation. We show that transcriptional cross-activation occurs also in strains deficient for Lon, ClpP, and HslV proteases and polyphosphate kinase. Furthermore, we find that toxins cleave the TA mRNA in vivo, which is followed by degradation of the antitoxin-encoding fragments and selective accumulation of the toxin-encoding regions. We show that these accumulating fragments can be translated to produce more toxin.ConclusionTranscriptional activation followed by cleavage of the mRNA and disproportionate production of the toxin constitutes a possible positive feedback loop, which can fire other TA systems and cause bistable growth heterogeneity. Cross-interacting TA systems have a potential to form a complex network of mutually activating regulators in bacteria.

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.

The endoribonuclease toxins of the E. coli toxin-antitoxin systems arrest bacterial growth and protein synthesis by targeting cellular mRNAs. As an exception, E. coli MazF was reported to cleave also 16S rRNA at a single site and separate an anti-Shine-Dalgarno sequence-containing RNA fragment from the ribosome. We noticed extensive rRNA fragmentation in response to induction of the toxins MazF and MqsR, which suggested that these toxins can cleave rRNA at multiple sites. We adapted differential RNA-sequencing to map the toxin-cleaved 5'- and 3'-ends. Our results show that the MazF and MqsR cleavage sites are located within structured rRNA regions and, therefore, are not accessible in assembled ribosomes. Most of the rRNA fragments are located in the aberrant ribosomal subunits that accumulate in response to toxin induction and contain unprocessed rRNA precursors. We did not detect MazF- or MqsR-cleaved rRNA in stationary phase bacteria and in assembled ribosomes. Thus, we conclude that MazF and MqsR cleave rRNA precursors before the ribosomes are assembled and potentially facilitate the decay of surplus rRNA transcripts during stress.

No abstract

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