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.

Within the larger ABC superfamily of ATPases, ABCF family members eEF3 in Saccharomyces cerevisiae and EttA in Escherichia coli have been characterised as ribosomal translation factors.Several other ABCFs including VgaA and LsaA confer resistance to MLS-type ribosome-targeting antibiotics. However, the diversity of ABCF subfamilies, the relationships among subfamilies and the evolution of antibiotic resistance factors from other ABCFs have not been explored. To address this, we analysed the presence of ABCFs and their domain architectures in 4505 genomes across the tree of life. We find that there are 45 distinct subfamilies of ABCFs, which are widespread across bacterial and eukaryotic phyla, suggesting they were present in the last common ancestor of both. Surprisingly, currently known antibiotic resistance (ARE) ABCFs are not confined to a distinct lineage of the ABCF family tree. This suggests that either antibiotic resistance is a general feature of bacterial ABCFs, or it is relatively easy to evolve antibiotic resistance from other ABCF functions. While eEF3 was thought to be fungi-specific, we have found eEF3-like factors in a range of single celled eukaryotes, suggesting an ancient origin in this domain of life. Finally, we address

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 first member of the pleuromutilin (PLM) class suitable for systemic antibacterial chemotherapy in humans recently entered clinical use, underscoring the need to better understand mechanisms of PLM resistance in disease-causing bacterial genera. Of the proteins reported to mediate PLM resistance in staphylococci, the least-well studied to date is Sal(A), a putative ABC-F NTPase that—by analogy to other proteins of this type—may act to protect the ribosome from PLMs. Here, we establish the importance of Sal proteins as a common source of PLM resistance across multiple species of staphylococci. Sal(A) is revealed as but one member of a larger group of Sal-type ABC-F proteins that vary considerably in their ability to mediate resistance to PLMs and other antibiotics. We find that specific sal genes are intrinsic to particular staphylococcal species, and show that this gene family is likely ancestral to the genus Staphylococcus. Finally, we solve the cryo-EM structure of a representative Sal-type protein (Sal(B)) in complex with the staphylococcal 70S ribosome, revealing that Sal-type proteins bind into the E site to mediate target protection, likely by displacing PLMs and other antibiotics via an allosteric mechanism.

Data in our present study revealed the notion that the epitopes from G and M proteins might be the target for peptide-based subunit vaccine design against HeV and NiV. However, the biochemical analysis is necessary to experimentally validate the interaction of epitopes individually with the MHC molecules through elucidation of immunity induction.