Under stressful conditions, bacterial RelA-SpoT Homolog (RSH) enzymes synthesize the alarmone (p)ppGpp, a nucleotide second messenger. (p)ppGpp rewires bacterial transcription and metabolism to cope with stress, and, at high concentrations, inhibits the process of protein synthesis and bacterial growth to save and redirect resources until conditions improve. Single-domain small alarmone synthetases (SASs) are RSH family members that contain the (p)ppGpp synthesis (SYNTH) domain, but lack the hydrolysis (HD) domain and regulatory C-terminal domains of the long RSHs such as Rel, RelA, and SpoT. We asked whether analysis of the genomic context of SASs can indicate possible functional roles. Indeed, multiple SAS subfamilies are encoded in widespread conserved bicistronic operon architectures that are reminiscent of those typically seen in toxin−antitoxin (TA) operons. We have validated five of these SASs as being toxic (toxSASs), with neutralization by the protein products of six neighboring antitoxin genes. The toxicity of Cellulomonas marina toxSAS FaRel is mediated by the accumulation of alarmones ppGpp and ppApp, and an associated depletion of cellular guanosine triphosphate and adenosine triphosphate pools, and is counteracted by its HD domain-containing antitoxin. Thus, the ToxSAS–antiToxSAS system with its multiple different antitoxins exemplifies how ancient nucleotide-based signaling mechanisms can be repurposed as TA modules during evolution, potentially multiple times independently.

Within the larger ABC superfamily of ATPases, ABCF family members eEF3 in Saccharomyces cerevisiae and EttA in Escherichia coli have been found to function as ribosomal translation factors. Several other ABCFs including biochemically characterized VgaA, LsaA and MsrE confer resistance to antibiotics that target the peptidyl transferase center and exit tunnel of the ribosome. However, the diversity of ABCF subfamilies, the relationships among subfamilies and the evolution of antibiotic resistance (ARE) factors from other ABCFs have not been explored. To address this, we analyzed the presence of ABCFs and their domain architectures in 4505 genomes across the tree of life. We find 45 distinct subfamilies of ABCFs that are widespread across bacterial and eukaryotic phyla, suggesting that they were present in the last common ancestor of both. Surprisingly, currently known ARE ABCFs are not confined to a distinct lineage of the ABCF family tree, suggesting that ARE can readily evolve from other ABCF functions. Our data suggest that there are a number of previously unidentified ARE ABCFs in antibiotic producers and important human pathogens. We also find that ATPase-deficient mutants of all four E. coli ABCFs (EttA, YbiT, YheS and Uup) inhibit protein synthesis, indicative of their ribosomal function, and demonstrate a genetic interaction of ABCFs Uup and YheS with translational GTPase BipA involved in assembly of the 50S ribosome subunit. Finally, we show that the ribosome-binding resistance factor VmlR from Bacillus subtilis is localized to the cytoplasm, ruling out a role in antibiotic efflux.

Summary Analysis of conservation of gene neighbourhoods over different evolutionary levels is important for understanding operon and gene cluster evolution, and predicting functional associations. Our tool FlaGs (Flanking Genes) takes a list of NCBI protein accessions as input, clusters neighbourhood-encoded proteins into homologous groups using sensitive sequence searching, and outputs a graphical visualization of the gene neighbourhood and its conservation, along with a phylogenetic tree annotated with flanking gene conservation. FlaGs has demonstrated utility for molecular evolutionary analysis, having uncovered a new toxin-antitoxin system in prokaryotes and bacteriophages. The web version of FlaGs (webFlaGs) can optionally include a BLASTP search against a reduced RefSeq database to generate an input accession list and analyse neighbourhood conservation within the same run. Availability and implementation FlaGs can be downloaded from https://github.com/GCA-VH-lab/FlaGs or run online at http://www.webflags.se/. Supplementary information Supplementary data are available at Bioinformatics online.

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

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.