Automated discovery and phylogenetic analysis of new toxin-antitoxin systems
© Guglielmini et al; licensee BioMed Central Ltd. 2008
Received: 07 September 2007
Accepted: 25 June 2008
Published: 25 June 2008
Although often viewed as elements "at the service of" bacteria, plasmids exhibit replication and maintenance mechanisms that make them purely "selfish DNA" candidates. Toxin-antitoxin (TA) systems are a spectacular example of such mechanisms: a gene coding for a cytotoxic stable protein is preceded by a gene coding for an unstable antitoxin. The toxin being more stable than the antitoxin, absence of the operon causes a reduction of the amount of the latter relative to the amount of the former. Thus, a cell exhibiting a TA system on a plasmid is 'condemned' either not to loose it or to die.
Different TA systems have been described and classified in several families, according to similarity and functional parameters. However, given the small size and large divergence among TA system sequences, it is likely that many TA systems are not annotated as such in the rapidly accumulating NCBI database. To detect these putative TA systems, we developed an algorithm that searches public databases on the basis of predefined similarity and TA-specific structural constraints. This approach, using a single starting query sequence for each of the ParE, Doc, and VapC families, and two starting sequences for the MazF/CcdB family, identified over 1,500 putative TA systems. These groups of sequences were analyzed phylogenetically for a better classification and understanding of TA systems evolution.
The phylogenetic distributions of the newly uncovered TA systems are very different within the investigated families. The resulting phylogenetic trees are available for browsing and searching through a java program available at http://ueg.ulb.ac.be/tiq/.
Plasmids are autonomously-replicating extra-chromosomal circular DNA molecules usually nonessential for cell survival under non-selective conditions and widely distributed in prokaryotic cells. Because plasmids sometimes bear genes that provide bacteria with functions (such as virulence, resistance to drugs, the ability to exploit a specific source of carbon) that can be adaptive in variable environments, they are often viewed as elements "at the service of" a (intra- or inter-specific) pool of bacteria, thus allowing the long-term survival of these lines or species. However, their ability to autonomously replicate makes plasmids possible purely "selfish DNA" candidates. Indeed, some plasmids exhibit features that seem to be strictly restricted to mechanisms related to their maintenance in cell lines (through replication and partitioning mechanisms) or dispersal across cell lines or species (through conjugation mechanisms).
Different TA systems have been described and classified in several families, according to the target of the toxin and/or the nature of the protease that degrades the antitoxin . Recently, about 150 toxin genes have been separated into 4 groups on the basis of sequence or structure similarities and gene neighborhood criteria : the "families" relE/parE, mazF/kid/ccdB, and Doc, as well as the family of proteins sharing a "PIN-domain". On the basis of phylogenetic analyses, these families have been suggested to be non-homologous , i.e., the TA systems would have appeared at least four times independently during evolution.
As known TA systems, identified on different plasmids, phages, and prokaryotic (including archaeal) genomes are all very small and potentially very divergent (TA systems originated from one or a few very old radiations), we hypothesize that many TA systems might not be annotated as such in the NCBI database. However, given that (i) TA-bearing plasmids with broad host range can be found in multiple bacterial species, and (ii) most systems exhibit the structural organization outlined above, we predicted that many more descendent systems than previously described should be detected across a wide range of prokaryotic genomes and plasmids. To detect these putative TA systems, we developed an algorithm, implemented into a computer program, TAQ V1.0 (for "TA Query"), that searches public databases on the basis of predefined similarity and TA-specific structural constraints. Our algorithm is complementary to that implemented in RASTA-Bacteria . The latter first identifies sequences exhibiting conserved putative TA domains and then uses structural constrains to further restrict and score the resulting set of putative TA systems. Our approach, using a single starting query sequence for each of the ParE, Doc, and PIN families and two starting sequences for the MazF/CcdB family identified over 1,500 putative TA systems, of which many were unknown. These five groups of sequences are analyzed phylogenetically for a better classification and understanding of TA systems evolution.
Results and Discussion
Taxonomic distribution of the in-silico inferred toxins identified for the 5 TAQ runs.
In silico inferred toxins
Total in silico inferred toxins
Without antinodes poisons
Artificial & Plasmids
Although induced expression of the bacterial RelE toxin in yeast and in human cell lines indicated the broad potential activity of TA systems [24, 25], none of the sequences that met all the sequence similarity and structural criteria defined in our algorithm are found in eukaryotic genomes. On the other hand, other categories include eukaryotic sequences: e.g., Tetrahymena thermophyla (Alveolata), Debaryomyces hansenii (Fungi), and Dictyostelium discoideum (Mycetozoa) in the ParE "Bad poisons" category, and Cryptosporidium hominis (Alveolata) in the ParE "without antidote poisons" category; Aspergillus fumigatus (Ascomycota), Coccidioides immitis (Ascomycota) Macaca mulatta (Mammalia), Drosophila melanogaster (Insecta), Mus musculus (Mammalia), Homo sapiens (Mammalia),Gallus gallus (Aves), Rattus norvegicus (Mammalia), Pan troglodytes (Mammalia), Canis lupus (Mammalia), Bos taurus (Mammalia) in the Doc "Bad poisons" category. Some of these sequences contain a domain (Fic domain in the Drosophila and Mus sequences; HYPE domain in other mammals and the chicken) that has been suggested to be homologous to the Doc domain .
Figure 3b shows the assignment of the 22 CcdB "in-silico inferred toxins" to functional categories. Only two sequences are "unknown", all the others are CcdB annotated.
Figure 3c shows the assignment of the 86 Doc "in-silico inferred toxins" to functional categories: about 10% of the sequences are annotated as unknown, whereas all remaining sequences are annotated as "Doc".
The run that was started with a single member of the PIN family generated 1418 putative toxins, 544 bad poisons, and 611 without-antidote poisons. Among the 1418 putative toxins, 381 are known VapC toxins, 163 are unknown, and 556 are simply annotated as containing a PIN domain. However, as the PIN domain is not specific to TA systems, it is highly likely that the first step of our TAQ algorithm (i.e. the BLASTP search) returns many false positives (i.e., sequences that are not TA systems) that are not efficiently filtered out by the structural-constrain criteria. One example is the FlbT protein clearly involved in the flagellum biosynthesis ; note that one cannot exclude the possibility that some proteins regulating gene expression (e.g., the FlbT gene product may act as a negative regulator of the flagellin fljK gene expression ) originate from TA systems. Surprisingly, TAQ even identified known Phd antitoxins. Hence, to avoid comparing non-homologous sequences, we decided to decrease the E-value to 10-3. Under that setting, TAQ did not recover any FlbT, Phd, or any other sequence that can readily be identified as false positives.
We compared the results of TAQ to those obtained by RASTA-Bacteria , a program that searches for putative TA systems in a specific organism but that restricts the BLAST search to known TA domains before applying additional structural constraints (ORF size, occurrence of an operon). RASTA-Bacteria finds many more putative TA elements (including isolated putative poisons or putative antidotes) than does TAQ because the latter uses more stringent size criteria (and constraint for the presence of both a toxin and an antitoxin genes) to minimize the risk of false positives that would seriously jeopardize the phylogenetic analyses. RASTA-Bacteria and TAQ have different objectives as the former attempts to score putative Toxin or Antitoxin elements whereas the latter attempts to generate a phylogeny among poison sequences that are very likely to belong to real TA systems (because they are all associated to a putative antidote and they all meet multiple size and localization criteria).
For a meaningful evolutionary analysis of the large sets of proteins recovered by TAQ, an additional criterion (e.g., the use of RASTA-Bacteria to test for the presence of an antitoxin domain in the putative antitoxin sequence) would be warranted. Unfortunately, RASTA-Bacteria does not allow, in its present form, to perform automated multiple searches using batch files.
Ancestral sequences tree
After initially grouping sequences by query, we manually separated or merged groups for generating low ambiguity alignments. We used 35 groups for the ParE analysis, 12 groups for the MazF analysis, 7 groups for the Doc analysis, and 54 groups for the VapC analysis. Sequences within each group were then aligned and used to produce ML phylogenies. Then, the 35 inferred ML root sequences (MLRS) for the ParE analysis (12, 7, and 54 MLRS for the MazF, Doc and VapC analyses, respectively) were themselves aligned and analyzed phylogenetically (see Methods for details). The final trees and their branch support values are available through TIQ v1.0, a Java program, available at  (using the login "tiq" and password "yqgWrj.81"). TIQ allows browsing the trees, select branches, perform searches of NCBI annotation fields (such as sequence, taxa names, Global Identifiers (GI), etc.) and filter sequences according to the host taxonomy. The 22 CcdB sequences were simply incorporated into a single MrBayes analysis.
For the ParE analysis, after grouping hits on the basis of the query sequence(s) that generated them, we obtained 710 groups of which many were highly redundant. Using the algorithm described in the Methods section, we reduced that number to 111 partially-overlapping groups that collectively contain all putative TA systems uncovered here. We then inferred the ML phylogenetic relationships among sequences within each group, generating 111 trees that are partially overlapping in terms of the included sequence. The overlapping trees were then used as input for "supertree" inference (see Methods section). The strict consensus among the 4 best supertrees (score 310.49) is available through TIQ v1.0 (see above). Similarly, for the MazF analysis, the initial 227 highly redundant groups were reduced to 76 partially-overlapping groups that collectively contain all uncovered putative TA systems. The single best supertree (score 157.57) is available through TIQ v1.0. For the Doc analysis, the initial 86 highly redundant groups were reduced to 9 partially-overlapping groups that collectively contain all uncovered putative TA systems. The strict consensus among the 17 best supertrees (score 11.72) is available through the java program TIQ v1.0. Finally, for the VapC analysis, the initial 860 highly redundant groups were reduced to 208 partially-overlapping groups that collectively contain all uncovered putative TA systems. The 50% majority-rule consensus among the 115 best supertrees (score 246.91) with a cut-off value of 50% is available through the java program TIQ v1.0. No supertree analysis was run for the 22 CcdB sequences given the small size of that dataset.
Consensus among methods
Origin of the TA systems
Very little objective data is available to shed light on the origin of the TA systems. Considering the structural relative similarity (i.e., the criteria implemented in our software TAQ V1.0, see Methods section) among known systems, Gerdes  suggested that they all share a common ancestor. However, the evolutionary relationships and functional (dis)similarities among TA families are unclear. For example, the similarities between, on one hand, toxin sequences from the RelE family and, on the other hand, those from the ParE family (see refs [8, 31] and analyses above) demonstrate that the two families are homologous (i.e., they share a common ancestor). However, RelE and ParE proteins are thought to exert their toxic activity on different targets: mRNA cleavage for RelE and DNA gyrase for ParE . Similarly, MazF and CcdB proteins are thought to be homologous because they share the same basic tertiary structure [31, 32]. Finally, Schmidt et al.  described an "hybrid" TA system, i.e., whose antidote sequence is similar to the MazE antitoxin, whereas its toxin sequence is similar to RelE, providing a putative evolutionary link between the RelE/ParE and MazF/CcdB superfamilies.
Conversely, other scholars suggested that TA systems evolved several times independently. For example, on the basis of protein domains and "gene-neighborhood analysis", Anantharaman and Aravind  proposed that RelE/ParE, MazF/CcdB, Doc, and PIN form different TA superfamilies that have been assembled more than once during the evolution, from a limited pool of protein domains.
To test for the possible homology among TA superfamilies, we run TAQ using PSI-BLAST criteria (e-value of 10.0 and PSI-BLAST threshold of 0.1, 10 iterations) more permissive than in a classical TAQ run (see Methods). Using ParE as input, we found 4ζ toxins from the ω -ε -ζ system, as well as a sequence containing the PIN domain (previously detected in some VapC toxins  of the VapBC system ). This result suggests possible homologies among ParE/RelE, ζ, and VapC families, and the rapid accumulation of new prokaryotic genome sequences might bridge additional putative phylogenetic gaps among TA families. Moreover, our first run of TAQ with a VapC toxin as initial query retrieved some Phd antitoxins (see above).
We also input 303 sequences of putative toxins (139 referred as "toxins" and 164 referred as "unknown") uncovered by TAQ V1.0 into RASTA-Bacteria . When RASTA-Bacteria gave more than one possible domain per sequence, we only considered the domain with the highest score. Strikingly, seven of our putative poisons belonging to the RelE superfamily (Figure 5) were assigned non-RelE domains (a Doc, PemK, CcdB, VapC, and PIN domain for 3, 1, 1, 1, and 1 sequences, respectively) by RASTA-Bacteria. Although these results exhibit low scores, they might prove to be the first objective link, in terms of homology, among TA superfamilies. In-vivo experiments, identifying the nature of the domains by which these proteins exert their putative toxicity, would probably shed light on this exciting hypothesis (see below).
Figure 5 indicates that among the "in-silico inferred ParE toxins", i.e., the new putative TA systems (red lineages) identified by our search algorithm (Fig. 10) are widely distributed into the phylogenetic tree of the ParE family. Conversely, the new putative TA systems uncovered here for the MazF, Doc, and VapC families are more restricted in their phylogenetic localization (Fig. 6, 7 and 9). However, for the MazF family, the new putative TA systems are particularly ancient (hence, divergent) in the phylogeny of the family, emphasizing the efficiency of our algorithm for identifying new TA systems.
All analyses presented here include exclusively sequences that met all similarity and structural constraints implemented in TAQ V1.0, i.e., the "in-silico inferred toxin" sequences. However, many additional hits meet the similarity criterion but fail to meet either the size criterion or the presence of a putative antidote ("bad toxins" or "without-antidote putative toxins" sequences, respectively). It would be particularly interesting to extend the phylogenetic analyses performed here to these additional sequences, especially for those found in eukaryotic genomes, to assess their likelihood to represent functional or degenerated poison genes.
Putative toxins particularly divergent in the phylogenies presented in figure 5, to 9, as well as those making the possible link between different TA superfamilies (see above), must be tested in vivo to (i) authenticate their toxic nature, (ii) identify their mode of action, and/or (iii) possibly confirm their hybrid character (i.e., with the toxic mode of action from family A, and the antitoxin structure of family B). Finally, it would be worth investigating which of the two proteases known to degrade antitoxins are active against the antitoxins of the tested putative TA systems.
Search for TA homologs
The next step of the algorithm consists into using structural parameters to check if BLAST hits are likely to belong to TA systems. To this end, TAQ V1.0 uses the annotated NCBI file (listing all hits) generated by BLAST to identify, for each amino-acid sequence hit, the corresponding genomic region and nucleotide sequence. Sequences are then filtered according to the following successive criteria (Figure 10a): (i) the length of the putative toxin (L pT ) must be 60aa <L pT < 150aa, a range encompassing the length distribution of all known toxin genes plus 40 aa; (ii) an additional ORF (putative antitoxin) of length 40aa <L pA < 150aa must be found (using the appropriate translation table and the three phases) in the 410 bp interval downstream of the putative toxin ORF; (iii) the separation or overlap between the putative antitoxin and putative toxin ORFs (i.e., between the stop codon of the former and the start codon of the latter) must be of maximum 30 bp. Only the sequences meeting all the above conditions are stored as "in-silico inferred toxins" (and the corresponding antitoxins are stored in a separate database) for further analyses. For each run, we assume that all sequences retained are (i) TA systems (active or inactive), and (ii) homologous (i.e., share a common ancestor). Every new "in-silico inferred toxin" sequence is automatically fed as a query into the TAQ V1.0 software until it converges (i.e., until no new sequence is found, Figure 10b). BLAST hits that do not meet the size criterion are stored in a "bad toxin" database, whereas BLAST hits that meet the size criterion but cannot be associated to a putative antitoxin are stored in a "without-antidote putative toxins" database. It is indeed conceivable that real toxin genes can lose their activity, hence, their associated antitoxin, during evolution.
In order to produce a single set of homologous sequences, we performed one TAQ V1.0 search for each family of TA systems starting with a single poison (the parE toxin from the E. Coli RK2 plasmid for the relE/parE family, the MazF toxin from the delivery vector pIEF16S for the mazF/kid/ccdB family, the Doc toxin from the Enterobacteria phage P1 for the Doc family, and the VapC toxin from Leptospira interrogans [serovar Lai str. 56601] for the PIN domain proteins). Since we did not find any CcdB sequence during the search using the MazF toxin as query (despite that CcdB and MazF are thought to belong to the same family), we also run one search using CcdB from the F plasmid as the starting sequence.
Alignment and phylogenetic analyses
ML phylogeny inference was carried out with the Bayesian approach  implemented in MRBAYES 3.1.2 [43, 44] with the "mixed protein model". The Markov chain Monte Carlo search was run with 6 chains for 106 millions generations and an initial temperature of 0.1, with trees sampled every 100 generations (the first 2,500 trees were discarded as "burnin"). Hence, Bayesian posterior probabilities were estimated as the majority-rule consensus tree among the 7,500 last sampled trees.
As an alternative to the above-mentioned approach (performing phylogenetic analysis among MLRS inferred from sub-groups of sequences), we also applied the supertree method [45, 46]. To this end, we requested our program TAQ V1.0 to output the query sequence(s) that generated each hit. Different poison sequences were grouped when they had been generated by the same query sequence. As different queries can generate the same hit, the groups are partially overlapping. As many of these groups are highly redundant (some group are even fully included into others), we reduced the number of groups using the following rules. First, groups of less than 3 sequences are discarded. Second, groups are sorted by decreasing size; groups 2 to N are compared each with group 1 (the largest group); any group sharing more than 2 sequences with group 1 are discarded and group 1 is set aside; the procedure is then iterated using the remaining groups. Once the two last groups have been compared, we checked that all toxin sequences in the initial list were included in the union of all groups set aside. Note that for the MazF family, the maximum allowed overlap between groups was raised from 2 to 3. We then inferred the phylogenetic relationships among sequences (using MRBAYES 3.1.2 with the "mixed protein model") within each of the N groups, generating N trees that are partially overlapping in terms of the included sequence. The overlapping trees were then used as input for "supertree" inference (Fig. 11b) using CLANN v 3.0.3d . This program implements 5 different methods of supertree reconstruction, of which only the "Most Similar Supertree" approach allowed analysis of the RelE and VapC datasets in practical computing time. This hill-climbing heuristics consists into generating a supertree topology T i (i.e., containing all taxa as leaves) at iteration i by performing a random branch swapping (following the SPR algorithm, ) on the topology Ti-1generated at iteration i-1. T i is accepted as a new starting topology if its score (evaluated under an optimality criterion described below) is better than that of tree Ti-1, otherwise it is rejected and a new branch swapping is performed on Ti-1. The score of a proposed supertree topology is evaluated as follows. First a taxa distance matrix is computed for each of the N source trees: the distance between a pair of taxa is the number of nodes that separate them in the source tree. Second, each source tree topology is compared to the proposed supertree topology: for each comparison, the supertree topology is pruned of all the taxa that are not present in the source tree and the supertree distance matrix is computed. Third, the score of the supertree topology relative to a given source tree is computed as the absolute difference between the distance matrix of the latter and that of the corresponding pruned supertree. Finally, the full score of the proposed supertree topology is computed as the normalized sum of all scores. For each major toxin family (ParE, MazF, and Doc) we performed 10 CLANN (v 3.0.3d) runs and kept all equally best trees. For the VapC family, we performed 5 runs (each using a starting trees generated using various heuristics since CLANN did not manage to build a starting tree from our set of trees) and kept all equally best trees.
This work was supported by grants from the 'Communauté Française de Belgique' (ARC 1164/20022770), the National Fund for Scientific Research Belgium (FNRS), and the Université Libre de Bruxelles. JG is PhD candidate at the Fonds pour la formation à la Recherche dans l'Industrie et dans l'Agriculture (FRIA), Belgium. We thank Raphaël Helaers for assistance to JG in java programming.
- Ogura T, Hiraga S: Mini-F plasmid genes that couple host cell division to plasmid proliferation. Proc Natl Acad Sci USA. 1983, 80 (15): 4784-4788.PubMed CentralView ArticlePubMedGoogle Scholar
- Aizenman E, Engelberg-Kulka H, Glaser G: An Escherichia coli chromosomal "addiction module" regulated by guanosine [corrected] 3',5'-bispyrophosphate: a model for programmed bacterial cell death. Proc Natl Acad Sci USA. 1996, 93 (12): 6059-6063.PubMed CentralView ArticlePubMedGoogle Scholar
- Brown JM, Shaw KJ: A novel family of Escherichia coli toxin-antitoxin gene pairs. J Bacteriol. 2003, 185 (22): 6600-6608.PubMed CentralView ArticlePubMedGoogle Scholar
- Gotfredsen M, Gerdes K: The Escherichia coli relBE genes belong to a new toxin-antitoxin gene family. Mol Microbiol. 1998, 29 (4): 1065-1076.View ArticlePubMedGoogle Scholar
- Gronlund H, Gerdes K: Toxin-antitoxin systems homologous with relBE of Escherichia coli plasmid P307 are ubiquitous in prokaryotes. J Mol Biol. 1999, 285 (4): 1401-1415.View ArticlePubMedGoogle Scholar
- Couturier M, Bahassi el M, Van Melderen L: Bacterial death by DNA gyrase poisoning. Trends Microbiol. 1998, 6 (7): 269-275.View ArticlePubMedGoogle Scholar
- Gerdes K, Christensen SK, Lobner-Olesen A: Prokaryotic toxin-antitoxin stress response loci. Nat Rev Microbiol. 2005, 3 (5): 371-382.View ArticlePubMedGoogle Scholar
- Pandey DP, Gerdes K: Toxin-antitoxin loci are highly abundant in free-living but lost from host-associated prokaryotes. Nucleic Acids Res. 2005, 33 (3): 966-976.PubMed CentralView ArticlePubMedGoogle Scholar
- Christensen SK, Mikkelsen M, Pedersen K, Gerdes K: RelE, a global inhibitor of translation, is activated during nutritional stress. Proc Natl Acad Sci USA. 2001, 98 (25): 14328-14333.PubMed CentralView ArticlePubMedGoogle Scholar
- Jensen RB, Gerdes K: Programmed cell death in bacteria: proteic plasmid stabilization systems. Mol Microbiol. 1995, 17 (2): 205-210.View ArticlePubMedGoogle Scholar
- Salmon MA, Van Melderen L, Bernard P, Couturier M: The antidote and autoregulatory functions of the F plasmid CcdA protein: a genetic and biochemical survey. Mol Gen Genet. 1994, 244 (5): 530-538.View ArticlePubMedGoogle Scholar
- Smith AS, Rawlings DE: Autoregulation of the pTF-FC2 proteic poison-antidote plasmid addiction system (pas) is essential for plasmid stabilization. J Bacteriol. 1998, 180 (20): 5463-5465.PubMed CentralPubMedGoogle Scholar
- Bahassi EM, O'Dea MH, Allali N, Messens J, Gellert M, Couturier M: Interactions of CcdB with DNA gyrase. Inactivation of Gyra, poisoning of the gyrase-DNA complex, and the antidote action of CcdA. J Biol Chem. 1999, 274 (16): 10936-10944.View ArticlePubMedGoogle Scholar
- Bernard P, Couturier M: Cell killing by the F plasmid CcdB protein involves poisoning of DNA-topoisomerase II complexes. J Mol Biol. 1992, 226 (3): 735-745.View ArticlePubMedGoogle Scholar
- Delphi Genetics. [http://www.delphigenetics.com]
- Szpirer CY, Milinkovitch MC: Separate-component-stabilization system for protein and DNA production without the use of antibiotics. Biotechniques. 2005, 38 (5): 775-781.View ArticlePubMedGoogle Scholar
- Gerdes K: Toxin-antitoxin modules may regulate synthesis of macromolecules during nutritional stress. J Bacteriol. 2000, 182 (3): 561-572.PubMed CentralView ArticlePubMedGoogle Scholar
- Anantharaman V, Aravind L: New connections in the prokaryotic toxin-antitoxin network: relationship with the eukaryotic nonsense-mediated RNA decay system. Genome Biol. 2003, 4 (12): R81-PubMed CentralView ArticlePubMedGoogle Scholar
- Sevin EW, Barloy-Hubler F: RASTA-Bacteria: a web-based tool for identifying toxin-antitoxin loci in prokaryotes. Genome Biol. 2007, 8 (8): R155-PubMed CentralView ArticlePubMedGoogle Scholar
- Roberts RC, Helinski DR: Definition of a minimal plasmid stabilization system from the broad-host-range plasmid RK2. J Bacteriol. 1992, 174 (24): 8119-8132.PubMed CentralPubMedGoogle Scholar
- Zhang XZ, Yan X, Cui ZL, Hong Q, Li SP: mazF, a novel counter-selectable marker for unmarked chromosomal manipulation in Bacillus subtilis. Nucleic Acids Res. 2006, 34 (9): e71-PubMed CentralView ArticlePubMedGoogle Scholar
- Pullinger GD, Lax AJ: A Salmonella dublin virulence plasmid locus that affects bacterial growth under nutrient-limited conditions. Mol Microbiol. 1992, 6 (12): 1631-1643.View ArticlePubMedGoogle Scholar
- Lehnherr H, Maguin E, Jafri S, Yarmolinsky MB: Plasmid addiction genes of bacteriophage P1: doc, which causes cell death on curing of prophage, and phd, which prevents host death when prophage is retained. J Mol Biol. 1993, 233 (3): 414-428.View ArticlePubMedGoogle Scholar
- Kristoffersen P, Jensen GB, Gerdes K, Piskur J: Bacterial toxin-antitoxin gene system as containment control in yeast cells. Appl Environ Microbiol. 2000, 66 (12): 5524-5526.PubMed CentralView ArticlePubMedGoogle Scholar
- Yamamoto TA, Gerdes K, Tunnacliffe A: Bacterial toxin RelE induces apoptosis in human cells. FEBS Lett. 2002, 519 (1–3): 191-194.View ArticlePubMedGoogle Scholar
- Brown NL, Stoyanov JV, Kidd SP, Hobman JL: The MerR family of transcriptional regulators. FEMS Microbiol Rev. 2003, 27 (2–3): 145-163.View ArticlePubMedGoogle Scholar
- Schoenlein PV, Ely B: Characterization of strains containing mutations in the contiguous flaF, flbT, or flbA-flaG transcription unit and identification of a novel fla phenotype in Caulobacter crescentus. J Bacteriol. 1989, 171 (3): 1554-1561.PubMed CentralPubMedGoogle Scholar
- Mangan EK, Malakooti J, Caballero A, Anderson P, Ely B, Gober JW: FlbT couples flagellum assembly to gene expression in Caulobacter crescentus. J Bacteriol. 1999, 181 (19): 6160-6170.PubMed CentralPubMedGoogle Scholar
- TiQ v1.0. [http://ueg.ulb.ac.be/tiq/]
- Adams EN: N-trees as nestings: complexity, similarity, and consensu. Journal of Classification. 1986, 3: 299-317.View ArticleGoogle Scholar
- Hargreaves D, Santos-Sierra S, Giraldo R, Sabariegos-Jareno R, de la Cueva-Mendez G, Boelens R, Diaz-Orejas R, Rafferty JB: Structural and functional analysis of the kid toxin protein from E. coli plasmid R1. Structure. 2002, 10 (10): 1425-1433.View ArticlePubMedGoogle Scholar
- Kamada K, Hanaoka F, Burley SK: Crystal structure of the MazE/MazF complex: molecular bases of antidote-toxin recognition. Mol Cell. 2003, 11 (4): 875-884.View ArticlePubMedGoogle Scholar
- Schmidt O, Schuenemann VJ, Hand NJ, Silhavy TJ, Martin J, Lupas AN, Djuranovic S: prlF and yhaV Encode a New Toxin-Antitoxin System in Escherichia coli. J Mol Biol. 2007, 372 (4): 894-905.PubMed CentralView ArticlePubMedGoogle Scholar
- National Center for Biotechnology Information (NCBI). [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/]
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.View ArticlePubMedGoogle Scholar
- Altschul SF, Koonin EV: Iterated profile searches with PSI-BLAST–a tool for discovery in protein databases. Trends Biochem Sci. 1998, 23 (11): 444-447.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402.PubMed CentralView ArticlePubMedGoogle Scholar
- Unit of Evolutionary Genetics. [http://www.ulb.ac.be/sciences/ueg]
- Loytynoja A, Milinkovitch MC: A hidden Markov model for progressive multiple alignment. Bioinformatics. 2003, 19 (12): 1505-1513.View ArticlePubMedGoogle Scholar
- Gardner PP, Wilm A, Washietl S: A benchmark of multiple sequence alignment programs upon structural RNAs. Nucleic Acids Res. 2005, 33 (8): 2433-2439.PubMed CentralView ArticlePubMedGoogle Scholar
- Edwards RJ, Shields DC: GASP: Gapped Ancestral Sequence Prediction for proteins. BMC Bioinformatics. 2004, 5: 123-PubMed CentralView ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F, Nielsen R, Bollback JP: Bayesian inference of phylogeny and its impact on evolutionary biology. Science. 2001, 294 (5550): 2310-2314.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17 (8): 754-755.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574.View ArticlePubMedGoogle Scholar
- Bininda-Emonds OR: The evolution of supertrees. Trends Ecol Evol. 2004, 19 (6): 315-322.View ArticlePubMedGoogle Scholar
- Bininda-Emonds OR: Supertree construction in the genomic age. Methods Enzymol. 2005, 395: 745-757.View ArticlePubMedGoogle Scholar
- Creevey CJ, McInerney JO: Clann: investigating phylogenetic information through supertree analyses. Bioinformatics. 2005, 21 (3): 390-392.View ArticlePubMedGoogle Scholar
- Swofford DL, Olsen GJ, Waddell PJ, Hillis DM: Phylogenetic inference. molecular systematics. Edited by: Mable BK. 1996, Sinauer & Associates, Sunderland, UK, 407-514.Google Scholar
- Hayes F: A family of stability determinants in pathogenic bacteria. J Bacteriol. 1998, 180 (23): 6415-6418.PubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.