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

Toxin-antitoxin systems: why so many, what for?

Toxin-antitoxin (TA) systems are small genetic modules that are abundant in bacterial genomes. Three types have been described so far, depending on the nature and mode of action of the antitoxin component. While type II systems are surprisingly highly represented because of their capacity to move by horizontal gene transfer, type I systems appear to have evolved by gene duplication and are more constrained. Type III is represented by a unique example located on a plasmid. Type II systems promote stability of mobile genetic elements and might act at the selfish level. Conflicting hypotheses about chromosomally encoded systems, from programmed cell death and starvation-induced stasis to protection against invading DNA and stabilization of large genomic fragments have been proposed.

The PIN-domain ribonucleases and the prokaryotic VapBC toxin-antitoxin array

The PIN-domains are small proteins of ~130 amino acids that are found in bacteria, archaea and eukaryotes and are defined by a group of three strictly conserved acidic amino acids. The conserved three-dimensional structures of the PIN-domains cluster these acidic residues in an enzymatic active site. PIN-domains cleave single-stranded RNA in a sequence-specific, Mg²+- or Mn²+-dependent manner. These ribonucleases are toxic to the cells which express them and to offset this toxicity, they are co-expressed with tight binding protein inhibitors. The genes encoding these two proteins are adjacent in the genome of all prokaryotic organisms where they are found. This sequential arrangement of inhibitor-RNAse genes conforms to that of the so-called toxin-antitoxin (TA) modules and the PIN-domain TAs have been named VapBC TAs (virulence associated proteins, VapB is the inhibitor which contains a transcription factor domain and VapC is the PIN-domain ribonuclease). The presence of large numbers of vapBC loci in disparate prokaryotes has motivated many researchers to investigate their biochemical and biological functions. For example, the devastating human pathogen Mycobacterium tuberculosis has 45 vapBC loci encoded in its genome whereas its non-pathogenic relative, Mycobacterium smegmatis has just one vapBC operon. On another branch of the prokaryotic tree, the nitrogen-fixing symbiont of legumes, Sinorhizobium meliloti has 21 vapBC loci and at least one of these loci have been implicated in the regulation of growth in the plant nodule. A range of biological functions has been suggested for these operons and this review sets out to survey the PIN-domains and summarise the current knowledge about the vapBC TA systems and their roles in diverse bacteria.

Txe, an endoribonuclease of the enterococcal Axe-Txe toxin-antitoxin system, cleaves mRNA and inhibits protein synthesis

The axe-txe operon encodes a toxin-antitoxin (TA) pair, Axe-Txe, that was initially identified on the multidrug-resistance plasmid pRUM in Enterococcus faecium. In Escherichia coli, expression of the Txe toxin is known to inhibit cell growth, and co-expression of the antitoxin, Axe, counteracts the toxic effect of Txe. Here, we report the nucleotide sequence of pS177, a 39 kb multidrug-resistant plasmid isolated from vancomycin-resistant Ent. faecium, which harbours the axe-txe operon and the vanA gene cluster. RT-PCR analysis revealed that the axe-txe transcript is produced by strain S177 as well as by other vancomycin-resistant enteroccoci. Moreover, we determine the mechanism by which the Txe protein exerts its toxic activity. Txe inhibits protein synthesis in E. coli without affecting DNA or RNA synthesis, and inhibits protein synthesis in a cell-free system. Using in vivo primer extension analysis, we demonstrate that Txe preferentially cleaves single-stranded mRNA at the first base after an AUG start codon. We conclude that Txe is an endoribonuclease which cleaves mRNA and inhibits protein synthesis.

TubercuList--10 years after

TubercuList (http://tuberculist.epfl.ch/), the relational database that presents genome-derived information about H37Rv, the paradigm strain of Mycobacterium tuberculosis, has been active for ten years and now presents its twentieth release. Here, we describe some of the recent changes that have resulted from manual annotation with information from the scientific literature. Through manual curation, TubercuList strives to provide current gene-based information and is thus distinguished from other online sources of genome sequence data for M. tuberculosis. New, mostly small, genes have been discovered and the coordinates of some existing coding sequences have been changed when bioinformatics or experimental data suggest that this is required. Nucleotides that are polymorphic between different sources of H37Rv are annotated and gene essentiality data have been updated. A host of functional information has been gleaned from the literature and many new activities of proteins and RNAs have been included. To facilitate basic and translational research, TubercuList also provides links to other specialized databases that present diverse datasets such as 3D-structures, expression profiles, drug development criteria and drug resistance information, in addition to direct access to PubMed articles pertinent to particular genes. TubercuList has been and remains a highly valuable tool for the tuberculosis research community with >75,000 visitors per month.

Virulence factors of the oral spirochete Treponema denticola

There is compelling evidence that treponemes are involved in the etiology of several chronic diseases, including chronic periodontitis as well as other forms of periodontal disease. There are interesting parallels with other chronic diseases caused by treponemes that may indicate similar virulence characteristics. Chronic periodontitis is a polymicrobial disease, and recent animal studies indicate that co-infection of Treponema denticola with other periodontal pathogens can enhance alveolar bone resorption. The bacterium has a suite of molecular determinants that could enable it to cause tissue damage and subvert the host immune response. In addition to this, it has several non-classic virulence determinants that enable it to interact with other pathogenic bacteria and the host in ways that are likely to promote disease progression. Recent advances, especially in molecular-based methodologies, have greatly improved our knowledge of this bacterium and its role in disease.

TADB: a web-based resource for Type 2 toxin-antitoxin loci in bacteria and archaea

TADB (http://bioinfo-mml.sjtu.edu.cn/TADB/) is an integrated database that provides comprehensive information about Type 2 toxin-antitoxin (TA) loci, genetic features that are richly distributed throughout bacterial and archaeal genomes. Two-gene and much less frequently three-gene Type 2 TA loci code for cognate partners that have been hypothesized or demonstrated to play key roles in stress response, bacterial physiology and stabilization of horizontally acquired genetic elements. TADB offers a unique compilation of both predicted and experimentally supported Type 2 TA loci-relevant data and currently contains 10,753 Type 2 TA gene pairs identified within 1240 prokaryotic genomes, and details of over 240 directly relevant scientific publications. A broad range of similarity search, sequence alignment, genome context browser and phylogenetic tools are readily accessible via TADB. We propose that TADB will facilitate efficient, multi-disciplinary and innovative exploration of the bacteria and archaea Type 2 TA space, better defining presently recognized TA-related phenomena and potentially even leading to yet-to-be envisaged frontiers. The TADB database, envisaged as a one-stop shop for Type 2 TA-related research, will be maintained, updated and improved regularly to ensure its ongoing maximum utility to the research community.

Analyzing the regulatory role of the HigA antitoxin within Mycobacterium tuberculosis

Bacterial chromosomally encoded type II toxin-antitoxin (TA) loci may be involved in survival upon exposure to stress and have been linked to persistence and dormancy. Therefore, understanding the role of the numerous predicted TA loci within the human pathogen Mycobacterium tuberculosis has become a topic of great interest. Antitoxin proteins are known to autoregulate TA expression under normal growth conditions, but it is unknown whether they have a more global role in transcriptional regulation. This study focuses on analyzing the regulatory role of the M. tuberculosis HigA antitoxin. We first show that the M. tuberculosis higBA locus is functional within its native organism, as higB, higA, and Rv1957 were successfully deleted from the genome together while the deletion of higA alone was not possible. The effects of higB-Rv1957 deletion on M. tuberculosis global gene expression were investigated, and a number of potential HigA-regulated genes were identified. Transcriptional fusion and protein-DNA-binding assays were utilized to confirm the direct role of HigA in Rv1954A-Rv1957 repression, and the M. tuberculosis HigA DNA-binding motif was defined as ATATAGG(N(6))CCTATAT. As HigA failed to bind to the next-most-closely related motif within the M. tuberculosis genome, HigA may not directly regulate any other genes in addition to its own operon.

Regulation of phenotypic variability by a threshold-based mechanism underlies bacterial persistence

In the face of antibiotics, bacterial populations avoid extinction by harboring a subpopulation of dormant cells that are largely drug insensitive. This phenomenon, termed "persistence," is a major obstacle for the treatment of a number of infectious diseases. The mechanism that generates both actively growing as well as dormant cells within a genetically identical population is unknown. We present a detailed study of the toxin-antitoxin module implicated in antibiotic persistence of Escherichia coli. We find that bacterial cells become dormant if the toxin level is higher than a threshold, and that the amount by which the threshold is exceeded determines the duration of dormancy. Fluctuations in toxin levels above and below the threshold result in coexistence of dormant and growing cells. We conclude that toxin-antitoxin modules in general represent a mixed network motif that can serve to produce a subpopulation of dormant cells and to supply a mechanism for regulating the frequency and duration of growth arrest. Toxin-antitoxin modules thus provide a natural molecular design for implementing a bet-hedging strategy.

Part II: defining and quantifying individual and co-cultured intracellular proteomes of two thermophilic microorganisms by GeLC-MS2 and spectral counting

Probing the intracellular proteome of Thermotoga maritima and Caldicellulosiruptor saccharolyticus in pure and co-culture affords a global investigation into the machinery and mechanisms enduring inside the bacterial thermophilic cell at the time of harvest. The second of a two part study, employing GeLC-MS(2) a variety of proteins were confidently identified with <1% false discovery rate, and spectral counts for label-free relative quantification afforded indication of the dynamic proteome as a function of environmental stimuli. Almost 25% of the T. maritima proteome and 10% of the C. saccharolyticus proteome were identified. Through comparison of growth temperatures for T. maritima, a protein associated with chemotaxis was uniquely present in the sample cultivated at the non-optimal growth temperature. It is suspected that movement was induced due to the non-optimal condition as the organism may need to migrate in the culture to locate more nutrients. The inventory of C. saccharolyticus proteins identified in these studies and attributed to spectral counting, demonstrated that two CRISPR-associated proteins had increased expression in the pure culture versus the co-culture. Further focusing on this relationship, a C. saccharolyticus phage-shock protein was identified in the co-culture expanding a scenario that the co-culture had decreased antiviral resistance and accordingly an infection-related protein was present. Alterations in growth conditions of these bacterial thermophilic microorganisms offer a glimpse into the intricacy of microbial behavior and interaction.

The relBE2Spn toxin-antitoxin system of Streptococcus pneumoniae: role in antibiotic tolerance and functional conservation in clinical isolates

Type II (proteic) chromosomal toxin-antitoxin systems (TAS) are widespread in Bacteria and Archaea but their precise function is known only for a limited number of them. Out of the many TAS described, the relBE family is one of the most abundant, being present in the three first sequenced strains of Streptococcus pneumoniae (D39, TIGR4 and R6). To address the function of the pneumococcal relBE2Spn TAS in the bacterial physiology, we have compared the response of the R6-relBE2Spn wild type strain with that of an isogenic derivative, Delta relB2Spn under different stress conditions such as carbon and amino acid starvation and antibiotic exposure. Differences on viability between the wild type and mutant strains were found only when treatment directly impaired protein synthesis. As a criterion for the permanence of this locus in a variety of clinical strains, we checked whether the relBE2Spn locus was conserved in around 100 pneumococcal strains, including clinical isolates and strains with known genomes. All strains, although having various types of polymorphisms at the vicinity of the TA region, contained a functional relBE2Spn locus and the type of its structure correlated with the multilocus sequence type. Functionality of this TAS was maintained even in cases where severe rearrangements around the relBE2Spn region were found. We conclude that even though the relBE2Spn TAS is not essential for pneumococcus, it may provide additional advantages to the bacteria for colonization and/or infection.

The Yersinia pestis chromosome encodes active addiction toxins

Toxin-antitoxin (TA) loci consist of two genes in an operon, encoding a stable toxin and an unstable antitoxin. The expression of toxin leads to cell growth arrest and sometimes bacterial death, while the antitoxin prevents the cytotoxic activity of the toxin. In this study, we show that the chromosome of Yersinia pestis, the causative agent of plague, carries 10 putative TA modules and two solitary antitoxins that belong to five different TA families (HigBA, HicAB, RelEB, Phd/Doc, and MqsRA). Two of these toxin genes (higB2 and hicA1) could not be cloned in Escherichia coli unless they were coexpressed with their cognate antitoxin gene, indicating that they are highly toxic for this species. One of these toxin genes (higB2) could, however, be cloned directly and expressed in Y. pestis, where it was highly toxic, while the other one (hicA1) could not, probably because of its extreme toxicity. All eight other toxin genes were successfully cloned into the expression vector pBAD-TOPO. For five of them (higB1, higB3, higB5, hicA2, and tox), no toxic activity was detected in either E. coli or Y. pestis despite their overexpression. The three remaining toxin genes (relE1, higB4, and doc) were toxic for E. coli, and this toxic activity was abolished when the cognate antitoxin was coexpressed, showing that these three TA modules are functional in E. coli. Curiously, only one of these three toxins (RelE1) was active in Y. pestis. Cross-interaction between modules of the same family was observed but occurred only when the antitoxins were almost identical. Therefore, our study demonstrates that of the 10 predicted TA modules encoded by the Y. pestis chromosome, at least 5 are functional in E. coli and/or in Y. pestis. This is the first demonstration of active addiction toxins produced by the plague agent.

Interaction specificity, toxicity and regulation of a paralogous set of ParE/RelE-family toxin-antitoxin systems

Toxin-antitoxin (TA) gene cassettes are widely distributed across bacteria, archaea and bacteriophage. The chromosome of the alpha-proteobacterium, Caulobacter crescentus, encodes eight ParE/RelE-superfamily toxins that are organized into operons with their cognate antitoxins. A systematic genetic analysis of these parDE and relBE TA operons demonstrates that seven encode functional toxins. The one exception highlights an example of a non-functional toxin pseudogene. Chromosomally encoded ParD and RelB proteins function as antitoxins, inhibiting their adjacently encoded ParE and RelE toxins. However, these antitoxins do not functionally complement each other, even when overexpressed. Transcription of these paralogous TA systems is differentially regulated under distinct environmental conditions. These data support a model in which multiple TA paralogs encoded by a single bacterial chromosome form independent functional units with insulated protein-protein interactions. Further characterization of the parDE(1) system at the single-cell level reveals that ParE(1) toxin functions to inhibit cell division but not cell growth; residues at the C-terminus of ParE(1) are critical for its stability and toxicity. While continuous ParE(1) overexpression results in a substantial loss in cell viability at the population level, a fraction of cells escape toxicity, providing evidence that ParE(1) toxicity is not uniform within clonal cell populations.

New toxins homologous to ParE belonging to three-component toxin-antitoxin systems in Escherichia coli O157:H7

Type II toxin-antitoxin (TA) systems are considered as protein pairs in which a specific toxin is associated with a specific antitoxin. We have identified a novel antitoxin family (paaA) that is associated with parE toxins. The paaA-parE gene pairs form an operon with a third component (paaR) encoding a transcriptional regulator. Two paralogous paaR-paaA-parE systems are found in E. coli O157:H7. Deletions of the paaA-parE pairs in O157:H7 allowed us to show that these systems are expressed in their natural host and that PaaA antitoxins specifically counteract toxicity of their associated ParE toxin. For the paaR2-paaA2-parE2 system, PaaR2 and Paa2-ParE2 complex are able to regulate the operon expression and both are necessary to ensure complete repression. The paaR2-paaA2-parE2 system mediates ClpXP-dependent post-segregational killing. The PaaR2 regulator appears to be essential for this function, most likely by maintaining an appropriate antitoxin : toxin ratio in steady-state conditions. Ectopic overexpression of ParE2 is bactericidal and is not resuscitated by PaaA2 expression. ParE2 colocalizes with the nucleoid, while it is diffusely distributed in the cytoplasm when PaaA2 is coexpressed. This indicates that ParE2 interacts with DNA-gyrase cycling on DNA and that coexpression of PaaA2 antitoxin sequesters ParE2 away from its target by protein-protein complex formation.

The Escherichia coli mqsR and ygiT genes encode a new toxin-antitoxin pair

Toxin-antitoxin (TA) systems are plasmid- or chromosome-encoded protein complexes composed of a stable toxin and a short-lived inhibitor of the toxin. In cultures of Escherichia coli, transcription of toxin-antitoxin genes was induced in a nondividing subpopulation of bacteria that was tolerant to bactericidal antibiotics. Along with transcription of known toxin-antitoxin operons, transcription of mqsR and ygiT, two adjacent genes with multiple TA-like features, was induced in this cell population. Here we show that mqsR and ygiT encode a toxin-antitoxin system belonging to a completely new family which is represented in several groups of bacteria. The mqsR gene encodes a toxin, and ectopic expression of this gene inhibits growth and induces rapid shutdown of protein synthesis in vivo. ygiT encodes an antitoxin, which protects cells from the effects of MqsR. These two genes constitute a single operon which is transcriptionally repressed by the product of ygiT. We confirmed that transcription of this operon is induced in the ampicillin-tolerant fraction of a growing population of E. coli and in response to activation of the HipA toxin. Expression of the MqsR toxin does not kill bacteria but causes reversible growth inhibition and elongation of cells.

Abundance of type I toxin-antitoxin systems in bacteria: searches for new candidates and discovery of novel families

Small, hydrophobic proteins whose synthesis is repressed by small RNAs (sRNAs), denoted type I toxin-antitoxin modules, were first discovered on plasmids where they regulate plasmid stability, but were subsequently found on a few bacterial chromosomes. We used exhaustive PSI-BLAST and TBLASTN searches across 774 bacterial genomes to identify homologs of known type I toxins. These searches substantially expanded the collection of predicted type I toxins, revealed homology of the Ldr and Fst toxins, and suggested that type I toxin-antitoxin loci are not spread by horizontal gene transfer. To discover novel type I toxin-antitoxin systems, we developed a set of search parameters based on characteristics of known loci including the presence of tandem repeats and clusters of charged and bulky amino acids at the C-termini of short proteins containing predicted transmembrane regions. We detected sRNAs for three predicted toxins from enterohemorrhagic Escherichia coli and Bacillus subtilis, and showed that two of the respective proteins indeed are toxic when overexpressed. We also demonstrated that the local free-energy minima of RNA folding can be used to detect the positions of the sRNA genes. Our results suggest that type I toxin-antitoxin modules are much more widely distributed among bacteria than previously appreciated.

Comprehensive functional analysis of Mycobacterium tuberculosis toxin-antitoxin systems: implications for pathogenesis, stress responses, and evolution

Toxin-antitoxin (TA) systems, stress-responsive genetic elements ubiquitous in microbial genomes, are unusually abundant in the major human pathogen Mycobacterium tuberculosis. Why M. tuberculosis has so many TA systems and what role they play in the unique biology of the pathogen is unknown. To address these questions, we have taken a comprehensive approach to identify and functionally characterize all the TA systems encoded in the M. tuberculosis genome. Here we show that 88 putative TA system candidates are present in M. tuberculosis, considerably more than previously thought. Comparative genomic analysis revealed that the vast majority of these systems are conserved in the M. tuberculosis complex (MTBC), but largely absent from other mycobacteria, including close relatives of M. tuberculosis. We found that many of the M. tuberculosis TA systems are located within discernable genomic islands and were thus likely acquired recently via horizontal gene transfer. We discovered a novel TA system located in the core genome that is conserved across the genus, suggesting that it may fulfill a role common to all mycobacteria. By expressing each of the putative TA systems in M. smegmatis, we demonstrate that 30 encode a functional toxin and its cognate antitoxin. We show that the toxins of the largest family of TA systems, VapBC, act by inhibiting translation via mRNA cleavage. Expression profiling demonstrated that four systems are specifically activated during stresses likely encountered in vivo, including hypoxia and phagocytosis by macrophages. The expansion and maintenance of TA genes in the MTBC, coupled with the finding that a subset is transcriptionally activated by stress, suggests that TA systems are important for M. tuberculosis pathogenesis.

Tuberculosis (TB) continues to be a major global health problem, causing 2 million deaths every year. A hallmark of TB pathogenesis is that the bacilli can enter into a slow or non-growing state in response to the host immune system. Because these persistent bacteria are resistant to antibiotic treatment, efforts to eliminate TB from the human population must include therapies to target dormant organisms as they can eventually resume replication to cause active disease. How Mycobacterium tuberculosis, the causative agent of TB, alters its replication dynamics in response to host cues is not understood. Toxin-antitoxin (TA) systems, which may control persistence in other bacteria, are massively expanded in M. tuberculosis, suggesting that they are important for TB pathogenesis. Surprisingly, the vast majority of these numerous TA systems are conserved only in pathogenic mycobacteria, suggesting their acquisition was important in M. tuberculosis evolution. Of the 88 putative TA systems identified, we show that 30 are functional in mycobacteria. A subset of these systems is activated upon exposure to stresses encountered during infection, indicating that specific TA systems are involved in adapting to environmental cues in the host. These genes are promising candidates for the development of novel therapies to target persistent bacteria.

Prokaryotic homologs of Argonaute proteins are predicted to function as key components of a novel system of defense against mobile genetic elements

Background

In eukaryotes, RNA interference (RNAi) is a major mechanism of defense against viruses and transposable elements as well of regulating translation of endogenous mRNAs. The RNAi systems recognize the target RNA molecules via small guide RNAs that are completely or partially complementary to a region of the target. Key components of the RNAi systems are proteins of the Argonaute-PIWI family some of which function as slicers, the nucleases that cleave the target RNA that is base-paired to a guide RNA. Numerous prokaryotes possess the CRISPR-associated system (CASS) of defense against phages and plasmids that is, in part, mechanistically analogous but not homologous to eukaryotic RNAi systems. Many prokaryotes also encode homologs of Argonaute-PIWI proteins but their functions remain unknown.

Results

We present a detailed analysis of Argonaute-PIWI protein sequences and the genomic neighborhoods of the respective genes in prokaryotes. Whereas eukaryotic Ago/PIWI proteins always contain PAZ (oligonucleotide binding) and PIWI (active or inactivated nuclease) domains, the prokaryotic Argonaute homologs (pAgos) fall into two major groups in which the PAZ domain is either present or absent. The monophyly of each group is supported by a phylogenetic analysis of the conserved PIWI-domains. Almost all pAgos that lack a PAZ domain appear to be inactivated, and the respective genes are associated with a variety of predicted nucleases in putative operons. An additional, uncharacterized domain that is fused to various nucleases appears to be a unique signature of operons encoding the short (lacking PAZ) pAgo form. By contrast, almost all PAZ-domain containing pAgos are predicted to be active nucleases. Some proteins of this group (e.g., that from Aquifex aeolicus) have been experimentally shown to possess nuclease activity, and are not typically associated with genes for other (putative) nucleases. Given these observations, the apparent extensive horizontal transfer of pAgo genes, and their common, statistically significant over-representation in genomic neighborhoods enriched in genes encoding proteins involved in the defense against phages and/or plasmids, we hypothesize that pAgos are key components of a novel class of defense systems. The PAZ-domain containing pAgos are predicted to directly destroy virus or plasmid nucleic acids via their nuclease activity, whereas the apparently inactivated, PAZ-lacking pAgos could be structural subunits of protein complexes that contain, as active moieties, the putative nucleases that we predict to be co-expressed with these pAgos. All these nucleases are predicted to be DNA endonucleases, so it seems most probable that the putative novel phage/plasmid-defense system targets phage DNA rather than mRNAs. Given that in eukaryotic RNAi systems, the PAZ domain binds a guide RNA and positions it on the complementary region of the target, we further speculate that pAgos function on a similar principle (the guide being either DNA or RNA), and that the uncharacterized domain found in putative operons with the short forms of pAgos is a functional substitute for the PAZ domain.

Conclusion

The hypothesis that pAgos are key components of a novel prokaryotic immune system that employs guide RNA or DNA molecules to degrade nucleic acids of invading mobile elements implies a functional analogy with the prokaryotic CASS and a direct evolutionary connection with eukaryotic RNAi. The predictions of the hypothesis including both the activities of pAgos and those of the associated endonucleases are readily amenable to experimental tests.

Reviewers

This article was reviewed by Daniel Haft, Martijn Huynen, and Chris Ponting.