Next Article in Journal
A Theoretical Study of 8-Chloro-9-Hydroxy-Aflatoxin B1, the Conversion Product of Aflatoxin B1 by Neutral Electrolyzed Water
Next Article in Special Issue
Structure, Evolution, and Functions of Bacterial Type III Toxin-Antitoxin Systems
Previous Article in Journal / Special Issue
Toxin-Antitoxin Systems in Clinical Pathogens
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 8
Issue 7
10.3390/toxins8070228
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Type II Toxin–Antitoxin Systems in the Unicellular Cyanobacterium Synechocystis sp. PCC 6803

by
Stefan Kopfmann
,
Stefanie K. Roesch
and
Wolfgang R. Hess
*
Faculty of Biology, Institute of Biology 3, Genetics and Experimental Bioinformatics, University of Freiburg, Schänzlestr. 1, D-79104 Freiburg, Germany
*
Author to whom correspondence should be addressed.
Toxins 2016, 8(7), 228; https://doi.org/10.3390/toxins8070228
Submission received: 20 June 2016 / Accepted: 11 July 2016 / Published: 21 July 2016
(This article belongs to the Special Issue Toxin-Antitoxin System in Bacteria)
Browse Figures
Versions Notes

Abstract

:
Bacterial toxin–antitoxin (TA) systems are genetic elements, which are encoded by plasmid as well as chromosomal loci. They mediate plasmid and genomic island maintenance through post-segregational killing mechanisms but may also have milder effects, acting as mobile stress response systems that help certain cells of a population in persisting adverse growth conditions. Very few cyanobacterial TA system have been characterized thus far. In this work, we focus on the cyanobacterium Synechocystis 6803, a widely used model organism. We expand the number of putative Type II TA systems from 36 to 69 plus seven stand-alone components. Forty-seven TA pairs are located on the chromosome and 22 are plasmid-located. Different types of toxins are associated with various antitoxins in a mix and match principle. According to protein domains and experimental data, 81% of all toxins in Synechocystis 6803 likely exhibit RNase activity, suggesting extensive potential for toxicity-related RNA degradation and toxin-mediated transcriptome remodeling. Of particular interest is the Ssr8013–Slr8014 system encoded on plasmid pSYSG, which is part of a larger defense island or the pSYSX system Slr6056–Slr6057, which is linked to a bacterial ubiquitin-like system. Consequently, Synechocystis 6803 is one of the most prolific sources of new information about these genetic elements.

Graphical Abstract

1. Introduction

Toxin–antitoxin (TA) systems are small genetic elements composed of a stable toxic protein and its unstable cognate antitoxin. They are encoded on chromosomes as well as on episomal genetic elements and are widely distributed throughout the prokaryotic domain of life. Many functions have been assigned to TA modules, ranging from plasmid maintenance to persister cell formation and stress response (for reviews, see [1,2,3,4,5,6,7,8]). TA systems exist in surprisingly high numbers in all prokaryotes and a growing number of studies suggest TA systems with milder effects to act as mobile stress response systems which help certain cells of a population in persisting adverse growth conditions. In general, the toxin is counteracted by its cognate antitoxin unless certain circumstances, such as plasmid loss or environmental stress, disrupt the TA balance and consequently frees the toxin to target various cellular functions. Whereas the toxin is always a protein, both proteins as well as RNA molecules have been found as antitoxins.
TA systems currently are classified into six types, based on the molecular nature of the antitoxins and the mechanisms by which they counteract the cognate toxins. Type I and type III are characterized by a small non-coding RNA (sRNA) as the antitoxin [9,10,11], whereas types II, IV, V and VI feature antitoxic proteins [1,12,13,14,15]. From these, Type II TA systems are certainly the best characterized systems with thousands of identified loci (reviewed in references [2,4,6]).
In this work, we focus on the unicellular cyanobacterium Synechocystis sp. PCC 6803 (from here: Synechocystis 6803), a model in biotechnology and fundamental research. DNA replication and metabolism, translation and cell membrane organization are among the verified cellular targets of bacterial TA systems. In cyanobacteria, oxygenic photosynthesis and thylakoid membranes appear as attractive additional targets, but very few cyanobacterial TA system have been characterized thus far. The numerically dominating TA systems belong to the Type II TA class. Makarova et al. (2009) [3] identified a total of 29 Type II TA pairs in Synechocystis 6803, whereas the TA Database [16,17] listed a total of 36 TA pairs, which are distributed over the chromosome and the four larger plasmids (Chr: 19; pSYSM: 5; pSYSX: 2; pSYSA: 6; pSYSG: 4). From these, only four chromosomally located TA pairs [18,19,20,21], one pair on plasmid pSYSX [22] and five located on plasmid pSYSA [23] were targeted experimentally thus far.
Addressing the number of possible Type II TA systems, here we expand the list of such putative systems in Synechocystis 6803 from 36 to 69 TA pairs and seven stand-alone TA components. From these 69 Type II TA loci, 47 are located on the chromosome and 22 are plasmid-located. With these numbers, Synechocystis 6803 is in the top-10 list of prokaryotic organisms with regard to the number of TA systems. We find that toxins belonging to one family are not always associated with antitoxins from a single family, as proposed in the classical classification system, but rather that different types of toxins are associated with various different antitoxins in a mix and match principle as proposed by Leplae et al. [4] and Makarova et al. [3]. According to the identified protein domains, 81% of all toxins in Synechocystis 6803 are predicted to exhibit RNase activity, suggesting extensive potential for toxicity-related RNA degradation but potentially also for toxin-mediated transcriptome remodeling. Indeed, the molecular characterization of 16 of these loci revealed various types of RNA endonuclease activities.

2. Results

2.1. Global Characterization of Synechocystis 6803 Type II TA Systems

Our results obtained in the course of this work indicated that TA systems in Synechocystis 6803 are much more abundant and diverse than previously thought. To identify undetected type II TA loci in Synechocystis 6803, we examined its genome and plasmids using the RASTA system [24], leading to the identification of seven new TA loci, which were not addressed in the TA database [17]. To extend this automated search, we scanned for putative TA systems based on the presence of known toxin and antitoxin associated protein domains [3,4,16,25,26]. We retrieved the gene IDs of Synechocystis 6803 protein sequences containing these domains via the Pfam database [27] and searched in a case-by-case analysis the genomic neighborhood of the predicted toxins or antitoxins for potential cognate TA partners. This led to the correction of six open reading frames by adjusting the 5′ ends according to multiple sequence alignments and information on the transcriptional start sites (TSSs; information extracted from references [28,29]) and the identification of seven novel TA-related open reading frames that had previously not been annotated in Synechocystis 6803 at all (Table 1).
We therefore decided to classify the 69 TA loci in Synechocystis 6803 based on the putative toxin component. We observed that the vapC, relE and hicA families appeared most frequent, while mazF is under-represented in Synechocystis 6803 compared to other bacteria (Figure 1).
A special situation exists with the Cas1-Cas2 module, which encodes the only conserved proteins among all subtypes of the CRISPR–Cas system and is typically connected to spacer acquisition [30]. While Cas2 is a homolog of the VapD toxin and exhibits sequence-specific nuclease activity [26,31], Makarova et al. [26] proposed that cas1-cas2 could act as a TA module and is therefore included here (Table 2).
Toxins belonging to one family are not always associated with antitoxins from a single antitoxin family, as proposed in the classic classification system, but rather TA loci are formed in a mix and match principle [3,4]. We also observed this relationship (Figure 2).

2.2. Verification of Predicted TA Functionality

To test the validity of predictions presented in Table 2, we selected four toxins for their ability to cause growth inhibition in E. coli that becomes attenuated by the associated antitoxin. We chose three chromosomally encoded TA systems and one located on plasmid pSYSG. To this aim, we expressed the toxin and antitoxin from different inducible promoters in E. coli. The different E. coli strains containing both pQE70::toxin (expression controlled by the addition of IPTG) and pBAD::antitoxin (expression controlled by the addition of arabinose) were grown to the required optical density in LB broth plus the respective antibiotics. The cultures were diluted 100-fold in fresh LB (plus antibiotics and the corresponding inducer) and the effects of separate toxin, antitoxin or toxin–antitoxin co-expression were analyzed by measuring the optical density over time in intervals of 15 min. Each of the four tested strains showed an inhibition in bacterial growth when the toxin was expressed alone, compared to strains with no expression or lone expression of the antitoxin (Figure 3). Co-expression of toxin and antitoxin attenuated the growth inhibitory effect for all four systems (Figure 3). Therefore, we conclude that the four tested systems exhibited typical properties of a TA system in vivo.
Hence, growth inhibition by the respective toxin and corresponding attenuation by the antitoxin has been demonstrated for altogether 11 Synechocystis 6803 TA systems: five located on plasmid pSYSA were previously analyzed [23], one located on plasmid pSYSG in this work (Figure 3) and five chromosomal TA systems analyzed before [18,20] or in this work (Figure 3).

2.3. Synechocystis 6803 Type II TA Systems of Special Interest

These analyses demonstrated interesting features of several of these previously known or newly defined Type II TA systems alike.
The slr1767–ssr2962 TA system genes are located in an excludon-like arrangement in the genome, in which they are covered over their entire length by a long transcript in antisense orientation (TU520 [29]) that otherwise represents the mRNAs of genes sll1639–sll1641 encoding the urease accessory protein UreD, nitrilase and glutamate decarboxylase. The excludon concept describes an unusual long asRNA, which acts as an asRNA regulator for one ORF while simultaneously comprising the ORFs and UTRs of the neighboring divergently oriented genes [34]. Furthermore, genes or operons organized in this arrangement often possess mutually exclusive or related functions [34]. All genes encoded in TU520 are involved in nitrogen metabolism, whereas the biological function of the TA module is unknown, but possesses a significant TSS only under nitrogen depletion [29]. According to the excludon concept [34], enhanced expression of the slr1767–ssr2962 dicistron might lead to operon discoordination and lowered expression of the glutamate decarboxylase gene compared to the first two genes within TU520. The Slr1767–Ssr2962 TA system is functional, and has previously been characterized as a bona fide TA system [20]. However, while this is the first report of a TA system in such an arrangement, the comparison to the closely related strain Synechocystis 6714 revealed the homologs of the genes slr1639–slr1641 to be conserved and within one operon, while the TA locus was missing. Thus, the slr1767–ssr2962 dicistron has been inserted or deleted in one of the two strains. It has to be evaluated whether the slr1767–ssr2962 TA module has a true biological function related to its excludon-like structure.
Another interesting and previously unrecognized variation of type II TA systems is the toxin–antitoxin fusion into a multi-domain protein. Such a situation was described with the EzeT system in E. coli in which the C-terminal domain of this type II TA system harbors a zeta-toxin like UDP-N-acetylglucosamine kinase domain that interferes with peptidoglycan synthesis, whereas the N-terminal domain harbors the antitoxin [35]. In Synechocystis 6803, we identified with Slr1999 an unusual putative TA component, consisting of a toxin (HicA) as well as an antitoxin (HicB) domain within one protein (Table 2). This fused protein probably resulted from the gene fusion of a conventional bicistronic operon, especially because there is no fused homolog detectable in any other organism. Another interesting aspect is that such a covalent linkage between toxin and antitoxin into a single polypeptide should lead to an auto-regulated toxin.

2.4. Synechocystis 6803 Type II TA Systems in the Maintenance of Plasmids and Genomic Islands

One of the proposed functions of TA systems is the maintenance of plasmids and genomic islands through the post-segregational killing mechanism. Furthermore, TA modules can be classified as anti-phage defense systems due to their ability to mediate dormancy and even cell death during phage attack [36,37,38]. In general, defense and TA systems can be found enriched in defense islands and on account of their addictive properties also in mobilome islands [39]. We inspected the localization and the genetic neighborhood of the putative TA loci to determine their potential association with plasmid, or genomic, defense and mobilome islands stabilization. All four large plasmids harbor at least three putative TA systems (Figure 1). However, the plasmid pSYSA encodes even ten different TA modules (one potential Type I, six Type II and three Cas1/Cas2) plus three stand-alone components.

2.4.1. The Slr6057–Slr6056 TA System May Stabilize a Genetic Cassette for a Bacterial Multiubiquitin System on Plasmid pSYSX

On plasmid pSYSX, Synechocystis 6803 possesses a cassette for a possible bacterial multiubiquitin system. The ubiquitin protein modification system is common to all eukaryotes, but is only rarely found in bacteria or archaea [40,41]. However, bacterial homologs to ubiquitin and the components of the ubiquitin modification system (E1, E2) have been found in divergent bacteria and archaea [42]. In Synechocystis 6803, these components are encoded on the plasmid pSYSX in a typical arrangement (multiubiquitin; E2; E1; metal binding protein) described by Iyer et al. [42], adjacent to the slr6057–slr6056 locus (Figure 4).
There are no homologs of this system nor divergent ubiquitin-like gene families in the closely related Synechocystis 6714, making this cassette likely to have originated from lateral gene transfer. Consistent with this idea, within cyanobacteria a homolog was detected only in distantly related species such as Pseudanabaena sp. PCC 7367, Oscillatoria nigro-viridis or Anabaena (Nostoc) sp. PCC 7120 (genes alr7502..alr7505) on plasmid beta [42], which were also located adjacent to a homologous TA locus (Figure 4). Therefore, we extended our screen to the genetic neighborhood of all genes in the Pfam database [43], coding for bacterial multiubiquitin domain (PF14452) proteins and their vicinity to TA modules. We observed that 70% of these proteins were adjacent to potential TA loci, in which all toxic components were characterized by the same COG2856/DUF955 domain we found in the Slr6057 toxin. Thus, a stabilization effect caused by the Slr6057–Slr6056 TA locus on this bacterial multiubiquitin cassette appears likely.

2.4.2. Is the hoxEFUYH Hydrogenase Operon in Synechocystis 6803 Stabilized by a TA System?

Our previous transcriptomic analyses [28,29] revealed that not all TA loci become transcribed from their own distinct TSS but some are part of a longer TU. Likewise, the Sll1225–Ssl2420 TA module and the TA associated DUF820-domain protein Sll1222 are encoded within the hox hydrogenase operon (hoxEFUYH). This operon is of special interest from a biotechnological perspective, since the bidirectional hydrogenase encoded by the hox operon could be used for biohydrogen production [44]. A very similar gene arrangement exists in Synechocystis 6714, where in extension of the conserved hoxEF_sll1222_hoxUY_sll1225_ssl2420_hoxH gene order, the TA system genes are interspersed with another hydrogenase-related gene, hypD (D082_26330), encoding the NiFe hydrogenase metallocenter assembly protein HypD, whereas the homologue in Synechocystis 6803, slr1498, is located 1.62 Mb away [45]. In contrast; other cyanobacteria seldom encode the hox operon in a coherent gene cluster according to the information obtained from the cyanobacterial genome database CyanoBase [46]. However, we identified in seven other strains TA modules or DUF820-domain proteins in the vicinity of some corresponding hox genes. In most cases, TA modules or DUF820-proteins clustered with the hoxEF operon, encoding part of the HoxEFU diaphorase subunit in cyanobacteria, which catalyzes the oxidation/reduction of NAD(P)H/NAD(P) [47]. HoxE can functionally be substituted by an unrelated subunit (HoxI) in other bacteria [47], but was shown to be an essential component for the connection between the NAD+/NADH and H+/H2 oxidoreduction reactions in Synechocystis [48]. In connection with the frequent clustering of TA systems and the hox operon, especially the hoxEF subunit, the question arises whether TA modules stabilize the whole or part of the hydrogenase hox operon. We consider the stabilization of the whole hox operon in Synechocystis as likely.

2.4.3. Two Type II TA Systems Are Part of the Large 40 kb Genomic Island in Synechocystis 6803

The slr1062–slr1062a (norf2) encoding a PemK-type TA system is central to TU355. This TU is with 22 genes transcribed into a single consecutive mRNA the longest multicistronic transcript in Synechocystis 6803 [29]. TU355 and the neighboring region comprise mainly genes usually not found in other cyanobacteria. Among the encoded functions are enzymes for the modification of cell surface structures and rfb-glycosyltransferases. The transcription into long consecutive mRNAs and the fact that these genes have no close homologs, even in the otherwise closely related Synechocystis sp. PCC 6714 [45] suggest that it is a genomic island. Moreover, there is with the gene pair ssr1765–ssr1766 a second Type II TA system in this region. Therefore, the location of slr1062–slr1062a and of ssr1765–ssr1766 is consistent with the idea that these stabilize this genomic island.

2.5. RNase Activity Is the Most Frequent Characteristics Associated with Type II TA Systems in Synechocystis 6803

According to the protein domains assigned to the here identified 69 Type II TA loci in Synechocystis 6803, 81% of all toxins are predicted to exhibit RNase activity, suggesting extensive potential for toxicity-related RNA degradation but potentially also for toxin-mediated transcriptome remodeling. To test this prediction experimentally, several of the predicted toxins were purified and their potential ribonuclease activity assayed. Initial results indicated a ribonuclease activity for the putative toxins Ssl2923 encoded on the chromosome, Slr5116 encoded on plasmid pSYSM, Ssl7007 and Ssl7039 encoded on plasmid pSYSA, and Sll8011 encoded on plasmid pSYSG. These proteins possess a PIN domain (Ssl2923, Slr5116, Ssl7007, Sll8011) or a RelE domain (Ssl7039). Furthermore, the tested toxins encoded on plasmid pSYSA Sll7003 and Ssl7039, exhibited RNase activity ([23] and Figure 5).

2.5.1. The RelE Protein Ssl7039 Is a Ribosome-Independent Ribonuclease in Vitro and Leads to the Degradation of E. coli lpp mRNA in Vivo

Co-expression of the respective toxin and antitoxin in vivo showed for multiple Synechocystis 6803 TA systems that the growth inhibition effect of the toxin was attenuated by the respective antitoxin ([18,20,23] and Figure 3). Therefore, we tested if the proposed toxic RNase activity of a purified toxin is abolished by the purified antitoxin. For this experiment, we chose the previously not studied Ssl7039–Ssl7038 TA system. RNase activity of Ssl7039 can be seen by a rapid increase of fluorescence over time after addition of the toxin, compared to an even level of fluorescence in the buffer and Ssl7038 antitoxin control. Simultaneous addition of both toxin and antitoxin resulted in a significant lower increase of fluorescence (Figure 5).
The RelE-family can be divided into five subfamilies, which show ribosome-dependent as well as ribosome-independent ribonuclease activity. The gp49-domain is assigned to the RelE-superfamily and has been suggested to represent a subclass of HigB [49]. The subfamily HigAB commonly features an inverted gene order and solely ribosome-dependent mRNA cleavage. In addition to the presence of a gp-49 domain in Ssl7039, the ssl7038–ssl7039 locus also exhibits an inverted gene order. To determine if the toxin Ssl7039 is exclusively ribosome-dependent we tested the ribonuclease activity of the purified toxin using a cell-free transcription/translation system and in vitro transcribed RNA, respectively. Addition of the purified toxin to the coupled transcription/translation approach, expressing the DHFR control template, caused a decrease in the amount of DHFR protein (Figure 6A). We assayed the template mRNA stability, in order to evaluate if protein translation arrest was mediated by the endoribonuclease activity of Ssl7039. Total RNA was extracted and separated by PAGE. Northern analysis confirmed degradation of dhfr mRNA, resulting in several specific cleavage products, exclusively after toxin addition (Figure 6B). To investigate if this endoribonucleolytic activity was ribosome dependent, we utilized in vitro synthesized dhfr mRNA as substrate. Substrate RNA was incubated with purified toxin and analyzed by PAGE separation and ethidium bromide staining. Free dhfr mRNA was also cleaved into several fragments by the toxin, resulting in a distinct cleavage pattern similar to the coupled cell-free transcription/translation assay (Figure 6C).
To extend these results further, we checked the stability of the lpp mRNA in vivo in the E. coli M15(pREP4)_pQE70::ssl7039 + pBAD::ssl7038 expression strains. The lpp transcript was chosen for examination because it is highly expressed and was used in several toxin cleavage assays before [50,51,52].
Total RNA was extracted and separated by PAGE. Northern analysis revealed a decreased transcript level of the lpp mRNA exclusively after Ssl7039 induction compared to no or Ssl7038 induction (Figure 6D). In addition, distinct degradation products were detected only after toxin expression. In contrast, co-expression of toxin and antitoxin resulted in a decreased transcript level with no distinct degradation products (Figure 6D). Therefore, the lpp expression level was lowered for other (non-toxin mediated) reasons, like, e.g., the metabolic burden from co-expressing both proteins. We conclude that Ssl7039 is a RelE-type ribosome-independent ribonuclease.

2.5.2. Heterologous Expression of Slr8014 in E. coli Causes Degradation of lpp mRNA in Vivo That Is Abolished by the Associated Antitoxin

We evaluated the potential RNase activity of the heterologously expressed toxin Slr8014 in vivo. For this aim, we performed northern analysis on lpp mRNA stability after toxin induction in the E. coli Bl21 (DE3) Rosetta pET28a(+)::slr8014 expression strain (Figure 7). Northern analysis showed time dependent lpp mRNA decay, which resulted in a prominent small and a faint large cleavage product compared to the non-induced strain (Figure 7A). Parallel growth inhibition assays showed that, compared to the strain with no expression, lone induction of Slr8014 caused significant growth inhibition (Figure 7B).
To extend these results further, we examined whether the cognate antitoxin Ssr8013 was able to neutralize the toxic Slr8014/VapC activity. The co-expression of toxin and antitoxin attenuated both the growth inhibitory effect (see Figure 3) and toxin RNase activity. Lone expression of antitoxin did not influence lpp mRNA integrity at any measured time. Co-expression of toxin and antitoxin led to increased lpp stability (Figure 7C). Compared to the solely expressed toxin, the prominent cleavage product was hardly seen 30 min past induction. Additionally, we observed a significant decrease of full length lpp mRNA after 180 min IPTG induction (Figure 7C).

3. Discussion

Here we present the comprehensive analysis of Type II TA systems in the cyanobacterial model organism Synechocystis 6803. Altogether, we define 76 likely Type II TA systems, roughly doubling the previously known number. From these, 69 occur in the classical form of a TA gene pair and seven as stand-alone components. We suggest three loci to be translated from leaderless mRNAs and seven to be associated with an asRNA. Altogether, 16 of the 69 two-component loci have been analyzed thus far by biochemical or molecular genetic approaches. The experimentally obtained results are consistent with the prediction that RNase activity is the major reason for the toxicity of these TA systems. In fact, based on the domain annotation it appears likely that 57 of the defined elements possess RNase activity.

3.1. Transcriptional Features of Synechocystis 6803 Type II TA Systems

Comparative analysis of the Synechocystis 6803 primary transcriptome established genome-wide TSS maps and consequently transcriptional units under different, environmentally relevant stimuli [28,29]. According to the criteria defined in these publications almost all of the putative TA pairs are transcribed, and distinct TSS could be identified for the majority of TA loci (Table 2).
We noticed that the TSS of the sll1651, sll1965 and ssl7039 toxin genes coincide with their respective start codons, pointing to their possible leader-less initiation of translation. Leader-less initiation of translation has been associated with the stress-dependent expression of certain genes [53,54]. Intriguingly, all TA loci affiliated to leader-less initiation of translation in Synechocystis, exhibit a reverse gene order.
This TSS data examination revealed furthermore the presence of internal TSSs within the 5′ segment of the coding region of the leading TA operon element in some cases. Based on these TSSs we suggest in this work possible alternative start codons (Table 1). In addition we noticed with 207 nt a quite long 5′ UTR for the sll7006–ssl7007 module. Moreover, the annotated start of ssl7007 is an alternative start codon (GTG), and there is not a single in frame stop codon within the 5′ UTR [23]. In general, blastp analysis of the possible translation products of the shortened or extended open reading frames revealed a very high degree of conserved residues with predicted protein sequences from several other cyanobacteria. Blast analysis of Ssr0756 also identified only an incomplete toxin domain, which could be completed by extending the ORF starting with an alternative start codon (Table 1). The prolonged open reading frame was further supported by the comparison of Synechocystis 6803 and 6714 homologs. We conclude that the genes listed in Table 1 should be re-annotated accordingly.
Type I TA systems are defined by their small non-coding RNA antitoxins, which are often located antisense to their cognate toxin. Genome-wide TSS mapping suggested massive antisense transcription in Synechocystis 6803 [28,29,55], consistent with the likely existence of multiple putative type I TA systems. Notwithstanding this possibility, we also detected seven Type II TA modules associated with antisense-located transcriptional units (Table 2). Usually, Type II TA systems are organized in a bicistronic operon and the transcriptional control is typically secured by autoregulation. Therefore, this organization could depict an additional regulatory element in Type II TA system control. Besides a possible antisense regulatory element we found other conspicuous features not commonly described for TA systems. We detected the putative TA systems sll0624a–sll0624, which exhibited mapped TSS for both the antitoxin and toxin component. On the other hand, the ssl0259–ssl2058 loci showed only a TSS for the subsequent antitoxin gene, whereas no TSS could be found for the leading toxin gene (Table 2). Three TA modules (ssl1300–sll0690; ssr0761–ssr0761a; and ssr2962–slr1767) held two alternative TSSs, indicating potential alternative transcripts (Table 2). For instance, the TA pair sll0690–ssl1300 displayed two TSSs which are differentially regulated under diverse environmental stimuli, showing the highest expression levels under heat stress and in the absence of light, respectively. Another eleven (ssl0350–ssl0350a; slr0770–slr0771; slr1062a–slr1062; slr1209–sl1210; ssl2420–sll1225; ssr1258–ssr1260; ssl2733–sll1400; ssr2067–ssr2066; ssl2920–ssl2921; ssl2922–ssl2923; and ssr3588–ssr3589) TA modules do not possess a distinct TSS, but are part of a larger TU (Table 2).

3.2. Similarities of Synechocystis 6803 Type II TA Systems to Those of Other Cyanobacteria

Based on their association with known mobile elements [56] or their frequently scattered distribution, TA systems are often considered as highly mobile genetic elements and, in many cases exist only in some isolates of a given species or genus [57,58,59]. Comparing the suite of Synechocystis 6803 Type II TA systems to those of Synechocystis 6714 (16S rDNA identity 99.4% [45,60]), we identified homologs for only 28% of toxin and 24% of antitoxin DNA sequences using blastn and a cut-off of 65% sequence identity. However, using blastp at a 30% cut-off of sequence identity, we detected 47% toxin and 37% antitoxin homologs. Furthermore, toxins and antitoxin are often small proteins defined by a single domain and we noticed that in some cases one and the same sequence had the same corresponding blastp hit. It cannot be excluded that these TA loci depict paralog systems, but we narrowed each case down to a single ortholog using reciprocal blastp analysis and synteny analysis of the genetic neighborhoods in strains 6803 and 6714, resulting finally in a number of 40% orthologous toxins and antitoxins, respectively (Supplementary Materials Tables S1 and S2). Kopf et al. [45] observed that approximately 77% of the Synechocystis 6803 protein-coding genes share homologs in Synechocystis 6714. Consequently, both strains share significantly less TA loci compared to the total number of shared protein coding genes, indicating a high mobility of the TA modules. This result is supported further by the observation that at least twelve TA loci insertion/deletion events happened for complete TA systems or components in one of the two strains. Due to their mobility, Type II TA systems are thought to move from one genome to another by horizontal gene transfer, resulting in a great variation of TA loci among bacterial genera and isolates [4]. As seen in the comparison of the closely related Synechocystis strains 6803 and 6714, several TA loci were not found in strain 6714, indicating insertion/deletion events or an alien origin. For instance, the blastp analysis of the slr1062a–slr1062 TA system located within the ~40 kb genomic island resulted in top-matching hits in non-cyanobacterial species and was most closely related to proteins in Thiocapsa marina (67% identity) and Halomonas sp. TG39a (62% identity), respectively. This diversity was observed for the complete island, pointing to an alien origin of the genomic region [45] stabilized by TA systems.
Furthermore, blastp analyses of Synechocystis 6803 TA loci against other cyanobacterial species, sharing the same clade based on 16S rRNA phylogeny [61], and the distantly related Synechocystis sp. PCC 7509 indicated a high variation of orthologous TA systems in these species. Based on the toxin sequence, seven orthologs were found in all closely related strains (Sll0205; Sll0286; Ssr0335; Ssl0350a; Slr1327; Sll0690; and Slr1906), whereas the TA systems sll0406–sll0405 and sll5003–sll5004 were not found in any other strain. However, in contrast to the situation for Slr1062–Slr1062a, the top-matching proteins for Ssl0406 and Sll5003 belonged to the cyanobacteria Synechococcus sp. PCC 7335 and Gloeocapsa sp. PCC 73106, respectively. TA systems occurring in all examined strains could be of special interest since the orthologs were restricted to cyanobacteria and absent from bacteria such as E. coli K-12 or Mycobacterium tuberculosis H37v.

3.3. Synechocystis 6803 Type II TA Systems with a Putative Function in the Maintenance of Plasmids and Genomic Islands

TA systems are frequently enriched on plasmids, in genomic islands and other mobile genetic elements [39] and this relation can also be seen in Synechocystis 6803. All four large plasmids harbor at least three putative TA systems and in the case of pSYSA even nine TA modules plus three stand-alone components (Figure 1). This extremely high number might in case of pSYSA be linked to its remaining functions: 75% of the pSYSA coding capacity encode three different CRISPR systems [33]. Besides the connection between the TA loci and the defense plasmid pSYSA, we observed the clustering of TA systems among each other, defense genes and mobile elements also on several other occasions. One potential defense island is located on plasmid pSYSG. This island encompasses 26 genes and is composed of three putative TA systems (ssl8008–sll8007; sll8012–sll8011; and ssr8013–slr8014), four genes associated with the type I restriction-modification system (ssl8049; sll8006; sll8009; and ssl8010); six mobile elements (ssl8041; slr8042–slr8046; and sll8017) and five genes related to TA modules (sll8001; sll8002; ssl8003; sll8004; ssl8005), which have not met our strict identification criteria. A possible mobilome island can be found on the chromosome, where three consecutive TA loci (ssl2922–ssl2923; ssl2920–ssl2921; and sll1504–sll1505) are interspersed between a partial transposase (ssr2699) and a mobile retron-type reverse transcriptase gene (sll1503).
Comparative genome analysis of the closely related Synechocystis strains 6803 and 714 revealed the occurrence of a large ~40 kb genomic island [45]. We examined whether the genomic island contained TA modules and noticed that two putative TA loci (slr1062a–slr1062; and ssr1765–ssr1766) were located within this region in Synechocystis 6803. It has to be mentioned that the here identified gene slr1062a has previously been annotated as a novel open reading frame (norf2) of unknown function on basis of transcriptomic evidence [28].
It is furthermore of high potential interest that functional gene cassettes, such as the bacterial multiubiquitin system on plasmid pSYSX or the chromosomally located hydrogenase operon appear to be stabilized by the presence of respective TA loci.

4. Conclusions

The comprehensive analysis of Type II TA systems in the cyanobacterial model organism Synechocystis 6803 revealed the presence of 76 likely Type II TA systems (including three CRISPR-Cas loci and seven stand-alone components). For comparison, analyzing the numbers and types of TA systems in the genomes of 10 different Xanthomonas strains, between five and 15 complete Type II TA loci were identified per strain [62], whereas the average total genome size in the genus Xanthomonas (~5 Mb [63]) is about one-third larger than Synechocystis 6803 (3.93 Mb, cf. Figure 1). The highest number of TA systems in TADB is predicted for the cyanobacterium, Microcystis aeruginosa NIES-843 with 113 loci [17]. Microcystis and Synechocystis are closely related genera. These numbers suggest that Synechocystis 6803 and similar cyanobacteria constitute one of the most prolific sources of new information about these genetic elements. From all Type II TA systems in Synechocystis 6803, 47 are encoded on the chromosome and altogether 57 were demonstrated or predicted to encode RNase activity. These two facts together suggest that in addition to the known functions, several of these systems could play a hitherto underestimated role in remodeling the transcriptome under certain conditions.

5. Materials and Methods

5.1. Comprehensive TA System Detection

Several sources were combined to identify novel Type II TA modules in Synechocystis 6803 in silico. First we employed the web-based identification tool RASTA-bacteria, an automated method allowing identification of TA loci in sequenced prokaryotic genomes, to predicted TA loci for the chromosome and all plasmids [24]. These loci were compared with Synechocystis 6803 TA modules listed in the TADB [16]. Additional predicted TA loci, with a probability score up to 50%, were evaluated case by case using blastp analysis and genetic neighborhood information. Furthermore, protein domains associated with TA systems were gathered by literature search [3,4,16,25,26,64,65,66,67]. Because of the differences between these programs, we also performed blastp searches (E < 10−4) in which all TA proteins were used as queries against the proteins from the studied Synechocystis genomes. Corresponding Synechocystis 6803 protein IDs, and eventually gene IDs, were retrieved from the Pfam database [43]. Again, potential TA components were evaluated by case to case analysis based on blastp and genomic neighborhood information. TA components were classified as part of a potential TA system, when a neighboring gene (up- or downstream) was located in close proximity and also featured a potential TA domain or was affiliated with TA systems. TA affiliation was evaluated based on the results of the blastp analysis and assigned when multiple homologous proteins were also adjacent to a TA component. Stand-alone TA components were identified by comparison with the genome sequence of the closely related strain Synechocystis 6714 [45,60] based on homology and genomic context. Information on TSSs was retrieved from their previous genome-wide mapping in Synechocystis 6803 [28,29] and Synechocystis 6714 [68].

5.2. Experimental Methods

5.2.1. Growth Inhibition Assay

The growth conditions for studying growth inhibition by toxins are given in [23].

5.2.2. Overexpression and Purification of Toxins and Antitoxins

The strains used for heterologous protein expression and primers for construction of expression plasmids (pQE70::Toxin, pET28a(+):Toxin, and pBAD::Antitoxin) are given in Tables S3 and S4. Overexpression of toxin (His)6 and antitoxin (His)6 was induced by adding IPTG (0.1 or 1 mM) for 3 h at 37 °C with shaking. Lysis and wash buffer composition was constant at 50 mM NaH2PO4, 300 mM NaCl, whereas best imidazole concentrations was experimentally tested for concentrations between 10 and 30 mM (lysis buffer) and 20–60 mM (wash buffer).

5.2.3. In Vitro Synthesis of RNA

In vitro transcription of the DHFR mRNA was carried out using the control plasmid of the PURExpress® In Vitro Protein Synthesis Kit (NEB) or from the PCR product using dhfr specific primers (Tables S3 and S4.).

5.2.4. RNA Extraction

E. coli total RNA extraction was performed following Hein et al. [69] with certain modifications. Cultures were harvested through centrifugation (8000 rpm, 20 °C, 15 min), resuspended in 1 mL PGTX [70] and flash-frozen in liquid nitrogen. Samples were incubated for 15 min at 65 °C and cooled down on ice before 1 vol chloroform/isoamylalcohol (24:1) was added. Samples were incubated at room temperature for 10 min under repetitive agitation. Phases were separated by centrifugation (swing out rotor at 6000 rpm, 15 °C, 15 min) and the upper aqueous phase was transferred to a new vial. Again, 1 vol chloroform/isoamylalcohol (24:1) was added and phases separated as described above. The upper aqueous phase was transferred to a new tube and gently mixed with 1 vol isopropanol. RNA was precipitated overnight at −20 °C and pelleted through centrifugation (13,000 rpm, 4 °C, 30 min). The pellet was washed with 1 mL of 70% ethanol (13,000 rpm, 20 °C, 5 min), allowed to air dry and resuspended in an appropriate volume of water, depending on favored RNA concentration.

5.2.5. RNase Activity

Detection of RNase activity was measured by fluorescence RNaseAlert™ essentially as described in Kopfmann et al. [23]. RelE toxin Sll7039 reaction buffer was at a final concentration composed of 25 mM Tris-HCl (pH 7.5).
Specific in vitro cleavage assays were performed as described previously [23]. Purified toxin and antitoxin, 100–250 ng of each, were incubated with 50–400 ng of target in vitro transcripts at 30 or 37 °C for 30–120 min in 25 mM Tris-HCl (pH 7.5), 60 mM KCl, 100 mM NH4Cl, 5 mM MgCl2, and 0.1 mM DTT and 25 mM Tris-HCl (pH 7.5) respectively.
In vivo cleavage Assay—E. coli BL21 (DE3): Rosetta pET28(a)::toxin strains and E. coli BL21 (DE3) pET28(a)::toxin + pBAD::antitoxin strains were inoculated and grown at 37 °C with shaking in LB medium containing appropriate antibiotics until an A600 0.18–0.3 was reached. Strains were grown at 37 °C with shaking in LB medium for another 30–180 min upon IPTG and L-arabinose addition. Cultures expressing only toxin, only antitoxin, both proteins, or none of these protein were harvested by centrifugation (8000 rpm, 15 min), resuspended in 1 mL PGTX [70], immediately frozen in liquid nitrogen and stored at −80 °C. Total RNA was extracted as described above, followed by DNase treatment (2 units of Turbo DNase, Invitrogen) and phenol/chloroform purification. Five micrograms of treated total RNA per sample were separated on 10% polyacrylamide-urea gels, electroblotted onto Hybond-N+ membranes from Amersham and cross-linked via UV-light exposure. Northern analysis was performed following Stazic et al. [71]. Membranes were hybridized overnight at 57 °C in hybridization buffer (50% deionized formamide, 7% SDS, 250 mM NaCl, 120 mM NaPi buffer, pH 7.2). Membranes were washed at 52 °C in prewarmed 2× SSC/1% SDS for 10 min, followed by 1× SSC/0.5% SDS for 10 min, and 1–10 min, according to signal intensity, in 0.1× SSC/0.1% SDS. Membranes were exposed to imaging plates, and the resulting signals were visualized using a Personal Molecular Imager FX system with Quantity One software (Bio-Rad, Munich, Germany).
Coupled Transcription/Translation Assay: For the coupled transcription/translation assays, we employed the PURExpress® In Vitro Protein Synthesis (NEB) cell-free transcription/translation system. Protein Synthesis Reaction was executed using the DHFR control template and supplemental 2.5–25 ng purified toxin Sll7039 or elution buffer respectively. Synthesis Reaction mixture was incubated at 37 °C for 2 h and reactions were stopped by placing the tubes on ice. Ten percent of the synthesis reaction were separated on 10% SDS–PAA gel electrophoresis and visualized by Coomassie-staining. The RNA of the remaining reaction mixtures was extracted by phenol/chloroform treatment followed by ethanol precipitation. RNA was separated on 10% polyacrylamide-urea gels, electroblotted onto Hybond-N+ membranes (Amersham) and cross-linked via UV-light exposure. Detection of DHFR mRNA transcript was performed by Northern analysis as described above.

Supplementary Materials

The following is available online at www.mdpi.com/2072-6651/8/7/228/s1, Table S1: Comparison of toxin and antitoxin DNA sequences using blastn and a cut-off of 65% sequence identity, Table S2: Homology search of toxin and antitoxin protein sequences using blastp and a cut-off of 30% sequence identity, Table S3: Strains, Table S4: Primers.

Acknowledgments

This work was supported by the Deutsche Forschungsgemeinschaft, Bonn, grants HE 2544/6-1 and HE 2544/8-2 (to W.R.H.). We thank Nicola Huber for her work on additional TA systems.

Author Contributions

W.R.H. and S.K. conceived and designed the experiments; S.K. and S.K.R. performed the experiments; W.R.H. and S.K. analyzed bioinformatics data; and W.R.H. and S.K. wrote the paper. All authors have read and agreed to the final version.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

CRISPRclustered regularly interspaced short palindromic repeats
DUFdomain of unknown function
TAtoxin–antitoxin
TADBtoxin–antitoxin database
TSStranscriptional start site
TUtranscriptional unit

References

  1. Schuster, C.F.; Bertram, R. Toxin-antitoxin systems are ubiquitous and versatile modulators of prokaryotic cell fate. FEMS Microbiol. Lett. 2013, 340, 73–85. [Google Scholar] [CrossRef] [PubMed]
  2. Chan, W.T.; Espinosa, M.; Yeo, C.C. Keeping the Wolves at Bay: Antitoxins of Prokaryotic Type II Toxin-Antitoxin Systems. Front. Mol. Biosci. 2016, 3. [Google Scholar] [CrossRef] [PubMed]
  3. Makarova, K.S.; Wolf, Y.I.; Koonin, E.V. Comprehensive comparative-genomic analysis of type 2 toxin-antitoxin systems and related mobile stress response systems in prokaryotes. Biol. Direct 2009, 4. [Google Scholar] [CrossRef] [PubMed]
  4. Leplae, R.; Geeraerts, D.; Hallez, R.; Guglielmini, J.; Drèze, P.; Van Melderen, L. Diversity of bacterial type II toxin-antitoxin systems: A comprehensive search and functional analysis of novel families. Nucleic Acids Res. 2011, 39, 5513–5525. [Google Scholar] [CrossRef] [PubMed]
  5. Page, R.; Peti, W. Toxin-antitoxin systems in bacterial growth arrest and persistence. Nat. Chem. Biol. 2016, 12, 208–214. [Google Scholar] [CrossRef] [PubMed]
  6. Rocker, A.; Meinhart, A. Type II toxin: Antitoxin systems. More than small selfish entities? Curr. Genet. 2016, 62, 287–290. [Google Scholar] [CrossRef] [PubMed]
  7. Schuster, C.F.; Bertram, R. Toxin-Antitoxin Systems of Staphylococcus aureus. Toxins 2016, 8, 140. [Google Scholar] [CrossRef] [PubMed]
  8. Yamaguchi, Y.; Park, J.-H.; Inouye, M. Toxin-antitoxin systems in bacteria and archaea. Annu. Rev. Genet. 2011, 45, 61–79. [Google Scholar] [CrossRef] [PubMed]
  9. Blower, T.R.; Short, F.L.; Rao, F.; Mizuguchi, K.; Pei, X.Y.; Fineran, P.C.; Luisi, B.F.; Salmond, G.P.C. Identification and classification of bacterial Type III toxin-antitoxin systems encoded in chromosomal and plasmid genomes. Nucleic Acids Res. 2012, 40, 6158–6173. [Google Scholar] [CrossRef] [PubMed]
  10. Fozo, E.M.; Makarova, K.S.; Shabalina, S.A.; Yutin, N.; Koonin, E.V.; Storz, G. Abundance of type I toxin-antitoxin systems in bacteria: Searches for new candidates and discovery of novel families. Nucleic Acids Res. 2010, 38, 3743–3759. [Google Scholar] [CrossRef] [PubMed]
  11. Brantl, S.; Jahn, N. sRNAs in bacterial type I and type III toxin-antitoxin systems. FEMS Microbiol. Rev. 2015, 39, 413–427. [Google Scholar] [CrossRef] [PubMed]
  12. Markovski, M.; Wickner, S. Preventing bacterial suicide: A novel toxin-antitoxin strategy. Mol. Cell 2013, 52, 611–612. [Google Scholar] [CrossRef] [PubMed]
  13. Aakre, C.D.; Phung, T.N.; Huang, D.; Laub, M.T. A Bacterial toxin inhibits DNA replication elongation through a direct interaction with the β sliding clamp. Mol. Cell 2013, 52, 617–628. [Google Scholar] [CrossRef] [PubMed]
  14. Brown, J.M.; Shaw, K.J. A novel family of Escherichia coli toxin-antitoxin gene pairs. J. Bacteriol. 2003, 185, 6600–6608. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, X.; Lord, D.M.; Cheng, H.-Y.; Osbourne, D.O.; Hong, S.H.; Sanchez-Torres, V.; Quiroga, C.; Zheng, K.; Herrmann, T.; Peti, W.; et al. A new type V toxin-antitoxin system where mRNA for toxin GhoT is cleaved by antitoxin GhoS. Nat. Chem. Biol. 2012, 8, 855–861. [Google Scholar] [CrossRef] [PubMed]
  16. Shao, Y.; Harrison, E.M.; Bi, D.; Tai, C.; He, X.; Ou, H.-Y.; Rajakumar, K.; Deng, Z. TADB: A web-based resource for Type 2 toxin-antitoxin loci in bacteria and archaea. Nucleic Acids Res. 2011, 39. [Google Scholar] [CrossRef] [PubMed]
  17. TADB. Available online: http://202.120.12.135/TADB2/index.php (accessed on 8 June 2016).
  18. Ning, D.; Ye, S.; Liu, B.; Chang, J. The proteolytic activation of the relNEs (ssr1114/slr0664) toxin-antitoxin system by both proteases Lons and ClpP2s/Xs of Synechocystis sp. PCC 6803. Curr. Microbiol. 2011, 63, 496–502. [Google Scholar] [CrossRef] [PubMed]
  19. Ning, D.; Zhao, W.; Qian, Y. A hypothetical gene pair, ssl2138 and sll11092, constitutes a functional TA system on the chromosome of Synechocystis PCC 6803. Wei Sheng Wu Xue Bao 2013, 53, 1043–1049. [Google Scholar] [PubMed]
  20. Ning, D.; Liu, S.; Xu, W.; Zhuang, Q.; Wen, C.; Tang, X. Transcriptional and proteolytic regulation of the toxin-antitoxin locus vapBC10 (ssr2962/slr1767) on the chromosome of Synechocystis sp. PCC 6803. PLoS ONE 2013, 8. [Google Scholar] [CrossRef] [PubMed]
  21. Krishna, P.S.; Rani, B.R.; Mohan, M.K.; Suzuki, I.; Shivaji, S.; Prakash, J.S.S. A novel transcriptional regulator, Sll1130, negatively regulates heat-responsive genes in Synechocystis sp. PCC6803. Biochem. J. 2013, 449, 751–760. [Google Scholar] [CrossRef] [PubMed]
  22. Čelešnik, H.; Tanšek, A.; Tahirović, A.; Vižintin, A.; Mustar, J.; Vidmar, V.; Dolinar, M. Biosafety of biotechnologically important microalgae: Intrinsic suicide switch implementation in cyanobacterium Synechocystis sp. PCC 6803. Biol. Open 2016, 5, 519–528. [Google Scholar] [CrossRef] [PubMed]
  23. Kopfmann, S.; Hess, W.R. Toxin-antitoxin systems on the large defense plasmid pSYSA of Synechocystis sp. PCC 6803. J. Biol. Chem. 2013, 288, 7399–7409. [Google Scholar] [CrossRef] [PubMed]
  24. Sevin, E.W.; Barloy-Hubler, F. RASTA-Bacteria: A web-based tool for identifying toxin-antitoxin loci in prokaryotes. Genome Biol. 2007, 8. [Google Scholar] [CrossRef] [PubMed]
  25. Blower, T.R.; Salmond, G.P.C.; Luisi, B.F. Balancing at survival’s edge: The structure and adaptive benefits of prokaryotic toxin-antitoxin partners. Curr. Opin. Struct. Biol. 2011, 21, 109–118. [Google Scholar] [CrossRef] [PubMed]
  26. Makarova, K.S.; Anantharaman, V.; Aravind, L.; Koonin, E.V. Live virus-free or die: Coupling of antivirus immunity and programmed suicide or dormancy in prokaryotes. Biol. Direct 2012, 7. [Google Scholar] [CrossRef] [PubMed]
  27. Finn, R.D.; Bateman, A.; Clements, J.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Heger, A.; Hetherington, K.; Holm, L.; Mistry, J.; et al. Pfam: The protein families database. Nucleic Acids Res. 2014, 42, D222–D230. [Google Scholar] [CrossRef] [PubMed]
  28. Mitschke, J.; Georg, J.; Scholz, I.; Sharma, C.M.; Dienst, D.; Bantscheff, J.; Voss, B.; Steglich, C.; Wilde, A.; Vogel, J.; Hess, W.R. An experimentally anchored map of transcriptional start sites in the model cyanobacterium Synechocystis sp. PCC6803. Proc. Natl. Acad. Sci. USA 2011, 108, 2124–2129. [Google Scholar] [CrossRef] [PubMed]
  29. Kopf, M.; Klähn, S.; Scholz, I.; Matthiessen, J.K.F.; Hess, W.R.; Voß, B. Comparative analysis of the primary transcriptome of Synechocystis sp. PCC 6803. DNA Res. 2014, 21, 527–539. [Google Scholar] [CrossRef] [PubMed]
  30. Nuñez, J.K.; Lee, A.S.Y.; Engelman, A.; Doudna, J.A. Integrase-mediated spacer acquisition during CRISPR-Cas adaptive immunity. Nature 2015, 519, 193–198. [Google Scholar] [CrossRef] [PubMed]
  31. Nuñez, J.K.; Kranzusch, P.J.; Noeske, J.; Wright, A.V.; Davies, C.W.; Doudna, J.A. Cas1-Cas2 complex formation mediates spacer acquisition during CRISPR-Cas adaptive immunity. Nat. Struct. Mol. Biol. 2014, 21, 528–534. [Google Scholar] [CrossRef] [PubMed]
  32. Trautmann, D.; Voß, B.; Wilde, A.; Al-Babili, S.; Hess, W.R. Microevolution in cyanobacteria: Re-sequencing a motile substrain of Synechocystis sp. PCC 6803. DNA Res. 2012, 19, 435–448. [Google Scholar] [CrossRef] [PubMed]
  33. Scholz, I.; Lange, S.J.; Hein, S.; Hess, W.R.; Backofen, R. CRISPR-Cas systems in the cyanobacterium Synechocystis sp. PCC6803 exhibit distinct processing pathways involving at least two Cas6 and a Cmr2 protein. PLoS ONE 2013, 8. [Google Scholar] [CrossRef] [PubMed]
  34. Sesto, N.; Wurtzel, O.; Archambaud, C.; Sorek, R.; Cossart, P. The excludon: A new concept in bacterial antisense RNA-mediated gene regulation. Nat. Rev. Microbiol. 2013, 11, 75–82. [Google Scholar] [CrossRef] [PubMed]
  35. Rocker, A.; Meinhart, A. A cis-acting antitoxin domain within the chromosomal toxin-antitoxin module EzeT of Escherichia coli quenches toxin activity. Mol. Microbiol. 2015, 97, 589–604. [Google Scholar] [CrossRef] [PubMed]
  36. Hazan, R.; Engelberg-Kulka, H. Escherichia coli mazEF-mediated cell death as a defense mechanism that inhibits the spread of phage P1. Mol. Genet. Genom. MGG 2004, 272, 227–234. [Google Scholar] [CrossRef] [PubMed]
  37. Fineran, P.C.; Blower, T.R.; Foulds, I.J.; Humphreys, D.P.; Lilley, K.S.; Salmond, G.P.C. The phage abortive infection system, ToxIN, functions as a protein-RNA toxin-antitoxin pair. Proc. Natl. Acad. Sci. USA 2009, 106, 894–899. [Google Scholar] [CrossRef] [PubMed]
  38. Dy, R.L.; Przybilski, R.; Semeijn, K.; Salmond, G.P.C.; Fineran, P.C. A widespread bacteriophage abortive infection system functions through a Type IV toxin-antitoxin mechanism. Nucleic Acids Res. 2014, 42, 4590–4605. [Google Scholar] [CrossRef] [PubMed]
  39. Makarova, K.S.; Wolf, Y.I.; Snir, S.; Koonin, E.V. Defense islands in bacterial and archaeal genomes and prediction of novel defense systems. J. Bacteriol. 2011, 193, 6039–6056. [Google Scholar] [CrossRef] [PubMed]
  40. Shigi, N. Posttranslational modification of cellular proteins by a ubiquitin-like protein in bacteria. J. Biol. Chem. 2012, 287, 17568–17577. [Google Scholar] [CrossRef] [PubMed]
  41. Humbard, M.A.; Miranda, H.V.; Lim, J.-M.; Krause, D.J.; Pritz, J.R.; Zhou, G.; Chen, S.; Wells, L.; Maupin-Furlow, J.A. Ubiquitin-like small archaeal modifier proteins (SAMPs) in Haloferax volcanii. Nature 2010, 463, 54–60. [Google Scholar] [CrossRef] [PubMed]
  42. Iyer, L.M.; Burroughs, A.M.; Aravind, L. The prokaryotic antecedents of the ubiquitin-signaling system and the early evolution of ubiquitin-like beta-grasp domains. Genome Biol. 2006, 7. [Google Scholar] [CrossRef] [PubMed]
  43. Finn, R.D.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Mistry, J.; Mitchell, A.L.; Potter, S.C.; Punta, M.; Qureshi, M.; Sangrador-Vegas, A.; et al. The Pfam protein families database: Towards a more sustainable future. Nucleic Acids Res. 2016, 44. [Google Scholar] [CrossRef] [PubMed]
  44. Appel, J. The Physiology and Functional Genomics of Cyanobacterial Hydrogenases and Approaches towards Biohydrogen Production. In Functional Genomics and Evolution of Photosynthetic Systems; Burnap, R., Vermaas, W., Eds.; Advances in Photosynthesis and Respiration; Springer Netherlands: Berlin, Germany, 2012; pp. 357–381. [Google Scholar]
  45. Kopf, M.; Klähn, S.; Pade, N.; Weingärtner, C.; Hagemann, M.; Voß, B.; Hess, W.R. Comparative genome analysis of the closely related Synechocystis strains PCC 6714 and PCC 6803. DNA Res. 2014, 21, 255–266. [Google Scholar] [CrossRef] [PubMed]
  46. Fujisawa, T.; Okamoto, S.; Katayama, T.; Nakao, M.; Yoshimura, H.; Kajiya-Kanegae, H.; Yamamoto, S.; Yano, C.; Yanaka, Y.; Maita, H.; et al. CyanoBase and RhizoBase: Databases of manually curated annotations for cyanobacterial and rhizobial genomes. Nucleic Acids Res. 2014, 42, D666–D670. [Google Scholar] [CrossRef] [PubMed]
  47. Eckert, C.; Boehm, M.; Carrieri, D.; Yu, J.; Dubini, A.; Nixon, P.J.; Maness, P.-C. Genetic analysis of the Hox hydrogenase in the cyanobacterium Synechocystis sp. PCC 6803 reveals subunit roles in association, assembly, maturation, and function. J. Biol. Chem. 2012, 287, 43502–43515. [Google Scholar] [CrossRef] [PubMed]
  48. Aubert-Jousset, E.; Cano, M.; Guedeney, G.; Richaud, P.; Cournac, L. Role of HoxE subunit in Synechocystis PCC6803 hydrogenase. FEBS J. 2011, 278, 4035–4043. [Google Scholar] [CrossRef] [PubMed]
  49. Jurėnaitė, M.; Markuckas, A.; Sužiedėlienė, E. Identification and characterization of type II toxin-antitoxin systems in the opportunistic pathogen Acinetobacter baumannii. J. Bacteriol. 2013, 195, 3165–3172. [Google Scholar] [CrossRef] [PubMed]
  50. Prysak, M.H.; Mozdzierz, C.J.; Cook, A.M.; Zhu, L.; Zhang, Y.; Inouye, M.; Woychik, N.A. Bacterial toxin YafQ is an endoribonuclease that associates with the ribosome and blocks translation elongation through sequence-specific and frame-dependent mRNA cleavage. Mol. Microbiol. 2009, 71, 1071–1087. [Google Scholar] [CrossRef] [PubMed]
  51. Rothenbacher, F.P.; Suzuki, M.; Hurley, J.M.; Montville, T.J.; Kirn, T.J.; Ouyang, M.; Woychik, N.A. Clostridium difficile MazF toxin exhibits selective, not global, mRNA cleavage. J. Bacteriol. 2012, 194, 3464–3474. [Google Scholar] [CrossRef] [PubMed]
  52. Sharp, J.D.; Cruz, J.W.; Raman, S.; Inouye, M.; Husson, R.N.; Woychik, N.A. Growth and translation inhibition through sequence-specific RNA binding by Mycobacterium tuberculosis VapC toxin. J. Biol. Chem. 2012, 287, 12835–12847. [Google Scholar] [CrossRef] [PubMed]
  53. Vesper, O.; Amitai, S.; Belitsky, M.; Byrgazov, K.; Kaberdina, A.C.; Engelberg-Kulka, H.; Moll, I. Selective translation of leaderless mRNAs by specialized ribosomes generated by MazF in Escherichia coli. Cell 2011, 147, 147–157. [Google Scholar] [CrossRef] [PubMed]
  54. Grill, S.; Moll, I.; Giuliodori, A.M.; Gualerzi, C.O.; Bläsi, U. Temperature-dependent translation of leaderless and canonical mRNAs in Escherichia coli. FEMS Microbiol. Lett. 2002, 211, 161–167. [Google Scholar] [CrossRef] [PubMed]
  55. Georg, J.; Voß, B.; Scholz, I.; Mitschke, J.; Wilde, A.; Hess, W.R. Evidence for a major role of antisense RNAs in cyanobacterial gene regulation. Mol. Syst. Biol. 2009, 5. [Google Scholar] [CrossRef] [PubMed]
  56. Iqbal, N.; Guérout, A.-M.; Krin, E.; Le Roux, F.; Mazel, D. Comprehensive functional analysis of the 18 Vibrio cholerae N16961 toxin-antitoxin systems substantiates their role in stabilizing the superintegron. J. Bacteriol. 2015, 197, 2150–2159. [Google Scholar] [CrossRef] [PubMed]
  57. Magnuson, R.D. Hypothetical functions of toxin-antitoxin systems. J. Bacteriol. 2007, 189, 6089–6092. [Google Scholar] [CrossRef] [PubMed]
  58. Pandey, D.P.; Gerdes, K. Toxin-antitoxin loci are highly abundant in free-living but lost from host-associated prokaryotes. Nucleic Acids Res. 2005, 33, 966–976. [Google Scholar] [CrossRef] [PubMed]
  59. Averina, O.V.; Alekseeva, M.G.; Abilev, S.K.; Il’in, V.K.; Danilenko, V.N. Distribution of genes of toxin-antitoxin systems of mazEF and relBE families in bifidobacteria from human intestinal microbiota. Genetika 2013, 49, 315–327. [Google Scholar] [CrossRef] [PubMed]
  60. Kopf, M.; Klähn, S.; Voss, B.; Stüber, K.; Huettel, B.; Reinhardt, R.; Hess, W.R. Finished genome sequence of the unicellular cyanobacterium Synechocystis sp. strain PCC 6714. Genome Announc. 2014, 2. [Google Scholar] [CrossRef] [PubMed]
  61. Shih, P.M.; Wu, D.; Latifi, A.; Axen, S.D.; Fewer, D.P.; Talla, E.; Calteau, A.; Cai, F.; de Marsac, N.T.; Rippka, R.; et al. Improving the coverage of the cyanobacterial phylum using diversity-driven genome sequencing. Proc. Natl. Acad. Sci. USA 2013, 110, 1053–1058. [Google Scholar] [CrossRef] [PubMed]
  62. Martins, P.M.M.; Machado, M.A.; Silva, N.V.; Takita, M.A.; de Souza, A.A. Type II Toxin-Antitoxin distribution and adaptive aspects on Xanthomonas genomes: Focus on Xanthomonas citri. Front. Microbiol. 2016, 7. [Google Scholar] [CrossRef] [PubMed]
  63. Sühnel, J. The Xanthomonas Genome Browser. Available online: http://xgb.fli-leibniz.de/cgi/index.pl (accessed on 14 June 2016).
  64. Guglielmini, J.; Van Melderen, L. Bacterial toxin-antitoxin systems: Translation inhibitors everywhere. Mob. Genet. Elem. 2011, 1, 283–306. [Google Scholar] [CrossRef] [PubMed]
  65. Holberger, L.E.; Garza-Sánchez, F.; Lamoureux, J.; Low, D.A.; Hayes, C.S. A novel family of toxin/antitoxin proteins in Bacillus species. FEBS Lett. 2012, 586, 132–136. [Google Scholar] [CrossRef] [PubMed]
  66. McNeil, M.B.; Iglesias-Cans, M.C.; Clulow, J.S.; Fineran, P.C. YgfX (CptA) is a multimeric membrane protein that interacts with the succinate dehydrogenase assembly factor SdhE (YgfY). Microbiology 2013, 159, 1352–1365. [Google Scholar] [CrossRef] [PubMed]
  67. Zhang, Y.; Inouye, M. RatA (YfjG), an Escherichia coli toxin, inhibits 70S ribosome association to block translation initiation. Mol. Microbiol. 2011, 79, 1418–1429. [Google Scholar] [CrossRef] [PubMed]
  68. Kopf, M.; Klähn, S.; Scholz, I.; Hess, W.R.; Voß, B. Variations in the non-coding transcriptome as a driver of inter-strain divergence and physiological adaptation in bacteria. Sci. Rep. 2015, 5. [Google Scholar] [CrossRef] [PubMed]
  69. Hein, S.; Scholz, I.; Voß, B.; Hess, W.R. Adaptation and modification of three CRISPR loci in two closely related cyanobacteria. RNA Biol. 2013, 10, 852–864. [Google Scholar] [CrossRef] [PubMed]
  70. Pinto, F.; Thapper, A.; Sontheim, W.; Lindblad, P. Analysis of current and alternative phenol based RNA extraction methodologies for cyanobacteria. BMC Mol. Biol. 2009, 10. [Google Scholar] [CrossRef] [PubMed]
  71. Stazic, D.; Lindell, D.; Steglich, C. Antisense RNA protects mRNA from RNase E degradation by RNA–RNA duplex formation during phage infection. Nucleic Acids Res. 2011, 39, 4890–4899. [Google Scholar] [CrossRef] [PubMed]
Toxins 08 00228 g001 1024
Figure 1. Frequency and distribution of Type II TA families in Synechocystis 6803. The classification of the 69 TA loci is based on the respective toxin component and given in %. The location of these loci over the chromosome and the four large plasmids is shown on the genomic maps. The color code of the TA loci matches the color code in the tabular overview. The length of plasmid pSYSA is according to reference [32]. Here, the three cas1-cas2 gene pairs on plasmid pSYSA were not counted.
Figure 1. Frequency and distribution of Type II TA families in Synechocystis 6803. The classification of the 69 TA loci is based on the respective toxin component and given in %. The location of these loci over the chromosome and the four large plasmids is shown on the genomic maps. The color code of the TA loci matches the color code in the tabular overview. The length of plasmid pSYSA is according to reference [32]. Here, the three cas1-cas2 gene pairs on plasmid pSYSA were not counted.
Toxins 08 00228 g001
Toxins 08 00228 g002 1024
Figure 2. Relationships between different toxin and antitoxin families in Synechocystis 6803, illustrated in a network graph inspired by Makarova et al. [3]. Toxin and antitoxin families are depicted in red and green, respectively. The edges connect toxin and antitoxin families that were identified as cognate toxic and antitoxic components organized in a putative TA operon. The thickness of an edge is proportional to the abundance of the toxin–antitoxin combination. The scale represents the correlation between operon abundance and line thickness.
Figure 2. Relationships between different toxin and antitoxin families in Synechocystis 6803, illustrated in a network graph inspired by Makarova et al. [3]. Toxin and antitoxin families are depicted in red and green, respectively. The edges connect toxin and antitoxin families that were identified as cognate toxic and antitoxic components organized in a putative TA operon. The thickness of an edge is proportional to the abundance of the toxin–antitoxin combination. The scale represents the correlation between operon abundance and line thickness.
Toxins 08 00228 g002
Toxins 08 00228 g003 1024
Figure 3. Effect of potential cyanobacterial toxin and antitoxin overexpression on growth of E. coli M15(pREP4) and Bl21(DE3) Rosetta™ harboring the plasmids pQE70::toxin or pET28a(+)::toxin and pBAD::antitoxin. Multiple TA systems were tested: genes sll0525–ssl1004; sll1092–ssl2138; and ssl2923–ssl2922; slr8014–slr8013. Expression of the potential toxins causes growth inhibition in the expression strains (red triangle), which can be prevented by co-expression of the respective antitoxin (purple crossed). Strains without expression are indicated by the blue diamond lines, and lone expression of antitoxin is indicated by the green square lines. The experiments were performed in duplicates or triplicates, and standard deviations are indicated. T-AT, toxin–antitoxin; O. D., optical density; T & AT co-expression, toxin and antitoxin co-expression. For studying slr8014–slr8013, the culture was grown in LB medium to an optical density (O. D.) of ~0.3 and then split (indicated by the black arrow). Subsequently, the effect of lone toxin expression (red triangle), antitoxin expression (green square lines), co-expression (purple crossed) and the untreated control (blue diamond) was measured as described above.
Figure 3. Effect of potential cyanobacterial toxin and antitoxin overexpression on growth of E. coli M15(pREP4) and Bl21(DE3) Rosetta™ harboring the plasmids pQE70::toxin or pET28a(+)::toxin and pBAD::antitoxin. Multiple TA systems were tested: genes sll0525–ssl1004; sll1092–ssl2138; and ssl2923–ssl2922; slr8014–slr8013. Expression of the potential toxins causes growth inhibition in the expression strains (red triangle), which can be prevented by co-expression of the respective antitoxin (purple crossed). Strains without expression are indicated by the blue diamond lines, and lone expression of antitoxin is indicated by the green square lines. The experiments were performed in duplicates or triplicates, and standard deviations are indicated. T-AT, toxin–antitoxin; O. D., optical density; T & AT co-expression, toxin and antitoxin co-expression. For studying slr8014–slr8013, the culture was grown in LB medium to an optical density (O. D.) of ~0.3 and then split (indicated by the black arrow). Subsequently, the effect of lone toxin expression (red triangle), antitoxin expression (green square lines), co-expression (purple crossed) and the untreated control (blue diamond) was measured as described above.
Toxins 08 00228 g003
Toxins 08 00228 g004 1024
Figure 4. Genetic organization of a likely bacterial multiubiquitin and functionally associated genes in connection to a TA system on plasmid pSYSX of Synechocystis 6803 (Sy.6803), compared to a similar system in Anabaena (Nostoc) sp. PCC 7120 (Ana.7120). mbD, metal binding domain; E1, functionally related to ubiquitin activating enzyme E1; E2, functionally related to ubiquitin transferring enzyme E2; m-Ubi , multiubiquitin; AT, antitoxin; T, toxin.
Figure 4. Genetic organization of a likely bacterial multiubiquitin and functionally associated genes in connection to a TA system on plasmid pSYSX of Synechocystis 6803 (Sy.6803), compared to a similar system in Anabaena (Nostoc) sp. PCC 7120 (Ana.7120). mbD, metal binding domain; E1, functionally related to ubiquitin activating enzyme E1; E2, functionally related to ubiquitin transferring enzyme E2; m-Ubi , multiubiquitin; AT, antitoxin; T, toxin.
Toxins 08 00228 g004
Toxins 08 00228 g005 1024
Figure 5. The RelE domain protein Ssl7039 exhibits ribonuclease activity, which is significantly hampered by antitoxin Ssl7038. The RNase Alert™ nuclease detection system was used. Fluorescence is indicated as function of time. The addition of the toxin Ssl7039 (arrow) shows a significant increase in fluorescence (red triangle) as compared with the buffer control (blue diamond). The simultaneous addition of toxin and antitoxin (arrow) results in a significantly lower increase in fluorescence (purple crossed). RFU, relative fluorescent units. The experiment was repeated independently (see also Figure 6).
Figure 5. The RelE domain protein Ssl7039 exhibits ribonuclease activity, which is significantly hampered by antitoxin Ssl7038. The RNase Alert™ nuclease detection system was used. Fluorescence is indicated as function of time. The addition of the toxin Ssl7039 (arrow) shows a significant increase in fluorescence (red triangle) as compared with the buffer control (blue diamond). The simultaneous addition of toxin and antitoxin (arrow) results in a significantly lower increase in fluorescence (purple crossed). RFU, relative fluorescent units. The experiment was repeated independently (see also Figure 6).
Toxins 08 00228 g005
Toxins 08 00228 g006 1024
Figure 6. Purified toxin Ssl7039 exhibits ribosome-independent RNase activity. (A) Addition of Ssl7039 to the coupled transcription-translation approach, expressing the DHFR control template, decreases DHFR protein abundance. The DHFR band is labeled by an arrow. The synthesis reaction were separated on 10% SDS–PAA gel electrophoresis and visualized by Coomassie staining. M, marker (PageRuler, Thermo Fisher Scientific); -T, no addition of purified Ssl7039; +T, addition of purified Ssl7039. (B) In vitro synthesized dhfr mRNA was incubated for 30 min with purified toxin Ssl7039, resulting in a clear cleavage pattern. Incubation with elution buffer (EB) served as control. The reaction were separated on 6%–10% 7 M urea polyacrylamide gels (6 mA, 1.5–2.5 h) and visualized by ethidium bromide staining. FL, full length dhfr mRNA (~500 nt). (C) Total RNA from the coupled transcription/translation approach was extracted 2 h after toxin addition, followed by northern analysis targeting the dhfr mRNA. EB, elution buffer control; Ssl7039, addition of purified Ssl7039. (D) Expression of Ssl7039 decreases lpp mRNA stability in vivo. The expression strain E. coli Bl21 (DE3) Rosetta, harboring the plasmid pET28a(+)::ssl7039 and pBAD::ssl7038, was grown in LB medium to an optical density (O. D.) of 0.3–0.4 and then split before addition of the appropriate inducer. Total RNA was extracted 30 min after inductor addition, followed by northern analysis aimed at the lpp mRNA. FL, full length lpp transcript; C, no inducer; T, toxin expression; AT, antitoxin expression; T & AT, both toxin and antitoxin expression.
Figure 6. Purified toxin Ssl7039 exhibits ribosome-independent RNase activity. (A) Addition of Ssl7039 to the coupled transcription-translation approach, expressing the DHFR control template, decreases DHFR protein abundance. The DHFR band is labeled by an arrow. The synthesis reaction were separated on 10% SDS–PAA gel electrophoresis and visualized by Coomassie staining. M, marker (PageRuler, Thermo Fisher Scientific); -T, no addition of purified Ssl7039; +T, addition of purified Ssl7039. (B) In vitro synthesized dhfr mRNA was incubated for 30 min with purified toxin Ssl7039, resulting in a clear cleavage pattern. Incubation with elution buffer (EB) served as control. The reaction were separated on 6%–10% 7 M urea polyacrylamide gels (6 mA, 1.5–2.5 h) and visualized by ethidium bromide staining. FL, full length dhfr mRNA (~500 nt). (C) Total RNA from the coupled transcription/translation approach was extracted 2 h after toxin addition, followed by northern analysis targeting the dhfr mRNA. EB, elution buffer control; Ssl7039, addition of purified Ssl7039. (D) Expression of Ssl7039 decreases lpp mRNA stability in vivo. The expression strain E. coli Bl21 (DE3) Rosetta, harboring the plasmid pET28a(+)::ssl7039 and pBAD::ssl7038, was grown in LB medium to an optical density (O. D.) of 0.3–0.4 and then split before addition of the appropriate inducer. Total RNA was extracted 30 min after inductor addition, followed by northern analysis aimed at the lpp mRNA. FL, full length lpp transcript; C, no inducer; T, toxin expression; AT, antitoxin expression; T & AT, both toxin and antitoxin expression.
Toxins 08 00228 g006
Toxins 08 00228 g007 1024
Figure 7. Expression of Slr8014 causes growth inhibition of E. coli and decreases lpp mRNA stability in vivo that is attenuated by the associated antitoxin. (A) Overexpression of the Synechocystis 6803 toxin Slr8014 in expression strain E. coli Bl21 (DE3) Rosetta harboring the plasmid pET28a(+)::slr8014. Total RNA was extracted at the indicated time points after IPTG addition, followed by northern analysis aimed at the lpp mRNA. Cleavage products are labeled by red arrows. FL, full length lpp transcript; ni, no IPTG induction; i, IPTG induction; nic + ic, negative controls (Bl21 (DE3) Rosetta strain harboring no plasmid). (B) Expression of the Synechocystis 6803 toxin Slr8014 causes strong growth inhibition in the expression strain E. coli Bl21 (DE3) Rosetta harboring plasmid pET28a(+)::slr8014 when the inducer IPTG was added (red triangles), but only a mild growth retardation when no IPTG is added (blue diamonds). Cultures continued to grow with and without IPTG addition in the no plasmid strain (yellow and blue dots). The experiments were performed in triplicates, and standard deviations are indicated. T-AT, toxin–antitoxin; O. D., optical density; T & AT co-expression, toxin and antitoxin co-expression. The culture was grown in LB medium to an optical density (O. D.) of 0.3-0.4 and then split (indicated by the black arrow). IPTG (final concentration 1 mM) was added to one-half (red triangle), and the other half was left untreated (blue diamond). Control treatment was analogously performed with the E. coli Bl21 (DE3) Rosetta expression strain harboring no plasmid. The experiments were performed in triplicates, and standard deviations are indicated. (C) Inhibition of the toxic RNase activity by co-expression in E. coli. Total RNA was extracted at the indicated time points after IPTG addition, followed by northern analysis aiming at lpp mRNA. Cleavage products are labeled by arrows. FL, full length lpp transcript; ni, no IPTG induction; i, IPTG induction; T, only toxin; AT, only antitoxin; T&AT, toxin and antitoxin co-expression.
Figure 7. Expression of Slr8014 causes growth inhibition of E. coli and decreases lpp mRNA stability in vivo that is attenuated by the associated antitoxin. (A) Overexpression of the Synechocystis 6803 toxin Slr8014 in expression strain E. coli Bl21 (DE3) Rosetta harboring the plasmid pET28a(+)::slr8014. Total RNA was extracted at the indicated time points after IPTG addition, followed by northern analysis aimed at the lpp mRNA. Cleavage products are labeled by red arrows. FL, full length lpp transcript; ni, no IPTG induction; i, IPTG induction; nic + ic, negative controls (Bl21 (DE3) Rosetta strain harboring no plasmid). (B) Expression of the Synechocystis 6803 toxin Slr8014 causes strong growth inhibition in the expression strain E. coli Bl21 (DE3) Rosetta harboring plasmid pET28a(+)::slr8014 when the inducer IPTG was added (red triangles), but only a mild growth retardation when no IPTG is added (blue diamonds). Cultures continued to grow with and without IPTG addition in the no plasmid strain (yellow and blue dots). The experiments were performed in triplicates, and standard deviations are indicated. T-AT, toxin–antitoxin; O. D., optical density; T & AT co-expression, toxin and antitoxin co-expression. The culture was grown in LB medium to an optical density (O. D.) of 0.3-0.4 and then split (indicated by the black arrow). IPTG (final concentration 1 mM) was added to one-half (red triangle), and the other half was left untreated (blue diamond). Control treatment was analogously performed with the E. coli Bl21 (DE3) Rosetta expression strain harboring no plasmid. The experiments were performed in triplicates, and standard deviations are indicated. (C) Inhibition of the toxic RNase activity by co-expression in E. coli. Total RNA was extracted at the indicated time points after IPTG addition, followed by northern analysis aiming at lpp mRNA. Cleavage products are labeled by arrows. FL, full length lpp transcript; ni, no IPTG induction; i, IPTG induction; T, only toxin; AT, only antitoxin; T&AT, toxin and antitoxin co-expression.
Toxins 08 00228 g007
Table
Table 1. Genomic location of TA module components re-annotated or newly defined in the course of this work. Initial and terminal positions are given according to the Synechocystis 6803 reference sequences for the chromosome (chr) and plasmid pSYSA (GenBank accession numbers NC_005230 and CP003267). Newly defined gene received the ID of the neighboring TA gene with the suffix ”a”. Gene slr1062a was previously identified as protein-coding gene norf2 (“novel ORF 2”) based on its expression and the conserved reading frame [28] but not recognized as TA-associated. Note that all nt positions refer to the forward strand, irrespectively of the location of the respective gene. All mentioned genes were found expressed in the comparative transcriptome analysis of Synechocystis 6803 [29].
Table 1. Genomic location of TA module components re-annotated or newly defined in the course of this work. Initial and terminal positions are given according to the Synechocystis 6803 reference sequences for the chromosome (chr) and plasmid pSYSA (GenBank accession numbers NC_005230 and CP003267). Newly defined gene received the ID of the neighboring TA gene with the suffix ”a”. Gene slr1062a was previously identified as protein-coding gene norf2 (“novel ORF 2”) based on its expression and the conserved reading frame [28] but not recognized as TA-associated. Note that all nt positions refer to the forward strand, irrespectively of the location of the respective gene. All mentioned genes were found expressed in the comparative transcriptome analysis of Synechocystis 6803 [29].
Gene IDAnnotationRe-Annotation
slr0725Chr:110808..111224Chr:110727..111224
ssr0756Chr:2084575..2084760Chr:2084455..2084760
ssr0761Chr:2085826..2086110Chr:2085859..2086110
ssl1004Chr:3262471..3262749Chr:3262471..3262707
sll1714Chr:968893..969294Chr:968893..969132
sll7006pSYSA:3367..3591pSYSA:3367..3768
ssl0350aChr:2738918..2739133novel ORF
sll0624aChr:3341880..3342128novel ORF
ssr0757aChr:2085211..2085417novel ORF
ssr0761aChr:2086107..2086211novel ORF
slr1062aChr:365077..365283novel ORF
ssl1493aChr:529033..529416novel ORF
ssr2784aChr:3474807..3474998novel ORF
Table
Table 2. Type II TA systems of Synechocystis 6803 according to the TADB [16], Makarova et al. 2009 [3], more specific studies, or newly defined according to this work (TW). Systems investigated in the course of this work biochemically are labeled “BC”. The respective gene loci, protein domains, the distance between the start codon and the start site of transcription (TSS) are given. Moreover, the replicon (chr, chromosome or plasmid name pSYS A, G, M or X) and the respective arrangement of toxin (T) and antitoxin (AT) component (Ord.) are indicated. Symbols and domain names: N/K, not known; COG, Cluster of Orthologous Genes; DUF, Domain of Unknown Function; TSS, Transcriptional Start Site (R, reverse strand; F, forward strand); TU, Transcriptional Unit; asRNA, antisense RNA; PIN, PIN domain (homologues of the ”pilT N“-terminal domain); Hic, Hic domain (hif-contiguous); YcfA, YcfA like domain; GNAT, GCN5-related .-acetyltransferases; RelE, RelE domain (relaxed response); YhaV, YhaV domain; YoeB, YoeB domain; PemK, PemK domain (plasmid emergency maintenance); MNT, MNT domain (minimal nucleotidyl transferase); HigB, HigB domain (Host Inhibition of Growth); cas, CRISPR-associated; HEPN; HEPN domain (higher eukaryotes and prokaryotes nucleotide-binding); PhdYeFM-antitox, prevent-host-death protein Phd/YefM antitoxin; HTH, helix-turn-helix; VagC, VagC domain (virulence associated gene); MazE, MazE domain; AbrB, AidB regulator domain; ChpI, chromosomal homolog of pem; Xre, xenobiotic response element. Information on TSSs was taken from published studies targeting the primary transcriptome of Synechocystis 6803 [28,29].
Table 2. Type II TA systems of Synechocystis 6803 according to the TADB [16], Makarova et al. 2009 [3], more specific studies, or newly defined according to this work (TW). Systems investigated in the course of this work biochemically are labeled “BC”. The respective gene loci, protein domains, the distance between the start codon and the start site of transcription (TSS) are given. Moreover, the replicon (chr, chromosome or plasmid name pSYS A, G, M or X) and the respective arrangement of toxin (T) and antitoxin (AT) component (Ord.) are indicated. Symbols and domain names: N/K, not known; COG, Cluster of Orthologous Genes; DUF, Domain of Unknown Function; TSS, Transcriptional Start Site (R, reverse strand; F, forward strand); TU, Transcriptional Unit; asRNA, antisense RNA; PIN, PIN domain (homologues of the ”pilT N“-terminal domain); Hic, Hic domain (hif-contiguous); YcfA, YcfA like domain; GNAT, GCN5-related .-acetyltransferases; RelE, RelE domain (relaxed response); YhaV, YhaV domain; YoeB, YoeB domain; PemK, PemK domain (plasmid emergency maintenance); MNT, MNT domain (minimal nucleotidyl transferase); HigB, HigB domain (Host Inhibition of Growth); cas, CRISPR-associated; HEPN; HEPN domain (higher eukaryotes and prokaryotes nucleotide-binding); PhdYeFM-antitox, prevent-host-death protein Phd/YefM antitoxin; HTH, helix-turn-helix; VagC, VagC domain (virulence associated gene); MazE, MazE domain; AbrB, AidB regulator domain; ChpI, chromosomal homolog of pem; Xre, xenobiotic response element. Information on TSSs was taken from published studies targeting the primary transcriptome of Synechocystis 6803 [28,29].
ToxinTypeA.toxinTypeTSSDis.Repl.Ord.ReferenceComment
sll0205PINssl0385COG5559/DUF22812516776 R27ChrAT-TTW-
ssl0258HicA (YcfA)ssl0259HicB (COG1598)2193884 R72ChrT-ATTWonly second component exhibits TSS
sll0286GNATsll0284HEPN (DUF86)3024345 R74ChrAT-TTW-
ssr0335HicA (YcfA)ssr0336HicB (COG1598)2754240 F401ChrT-ATTW-
ssl0350aHicA (YcfA)ssl0350HicB (COG1598)2740742 R1609ChrAT-TTWTSS of TU
sll0406PINsll0405hypothetical2552996 R34ChrAT-TTWpotential asRNA regulation
sll0525PINssl1004PhdYeFM_antitox326273730ChrAT-T[3,16]see Figure 3
sll0624PIN3sll0624aHicB (COG1598)3342723 R281ChrT-ATTWboth components exhibit TSS
sll0658PINssl1255HTH (COG2886)3121956 R3ChrAT-T[3,16]-
slr0664RelE (DUF1044)ssr1114hypothetical3347861 F80ChrAT-T[16,18]-
sll0690PINssl1300VagC2632007 R139ChrAT-T[3,16]alternative TSS (2631867 R)
slr0725YhaVslr0724MazE (AbrB)110366 F26ChrAT-T[3,16]-
ssl0739HicA (YcfA)ssl0738HicB (COG1598)N/KN/KChrAT-TTW-
ssr0756PINssr0755hypothetical2084115 F28ChrAT-TTW-
ssr0757DUF497/COG2929ssr0757aDUF4415/COG35142084877 F45ChrT-ATTW-
ssr0761aYoeBssr0761PhdYeFM_antitox2085833 F27ChrAT-TTW-
slr0771PINslr0770COG2442/Duf4332393421 F811ChrAT-T[3,16]joint tricistronic TU2489 with slr0769
slr1062PemKslr1062ahypothetical3577247353ChrAT-TTWTSS of TU
sll1092PINssl2138COG5559/DUF2281N/KN/KChrAT-T[16,19]see Figure 3
sll1130PemKssl2245hypothetical1049139 R123ChrAT-T[16,21]potential asRNA regulation; involved in the heat shock response
slr1210PINslr1209COG2442/Duf433891086 F2175ChrAT-T[3,16]TSS of TU; potential asRNA regulation
sll1225PINssl2420hypothetical1678733 R4435ChrAT-T[16]-
slr1241MNT slr1240hypothetical1053924 F1650ChrAT-TTWTSS of TU
ssr1260HicA (YcfA)ssr1258HicB (COG1598)3428302 F660ChrAT-TTWTSS of TU
slr1327PIN3ssr2201hypothetical1660041 F28ChrAT-T[16]-
ssl1377HicA (YcfA)ssl1376HicB (COG1598)3427318 R115ChrAT-TTW-
slr1383PINssr2317hypothetical679231 F459ChrT-ATTW-
sll1400PINssl2733COG2442/Duf43350262 R1627ChrAT-T[3,16]TSS of TU
slr1421DUF4258ssr2377HicB (COG1598)3273485 F24ChrAT-TTW-
ssl1493aPINssl1493PhdYeFM_antitox3044721 R66ChrAT-TTW-
sll1504MNT sll1505HEPN (DUF86)477658 R1915ChrT-AT[3,16]TSS of TU
sll1543PINssl2999hypothetical3475237 R53ChrT-ATTWboth components exhibit TSS
sll1651DUF497/COG2929sll1652DUF3680/COG5304282821 R0ChrT-ATTWleader-less
sll1715PINsll1714hypothetical969538 R244ChrAT-T[16]-
ssr1766HicA (YcfA)ssr1765HicB (COG1598)382338 F54ChrAT-T[3,16]-
slr1767PINssr2962COG2442/Duf433544995 F34ChrAT-T[3,16,20]alternative TSS (2085834 F); excludon-like arrangement
slr1906GNATslr1907hypothetical609442 F27ChrT-ATTW-
sll1912PINssl3615hypothetical776369 R27ChrAT-T[16]-
sll1965DUF497/COG2929ssl3719DUF4415/COG3514915052 R0ChrT-ATTWleader-less
slr1999HicA (YcfA)slr1999HicB (COG1598)1444698727ChrAT/TTWboth domains in one protein
ssr2066HicA (YcfA)ssr2067HicB (COG1598)1053924 F2563ChrT-ATTWTSS of TU
ssl2749MNT sll1411HEPN (DUF86)1619636 R1540ChrT-AT[3,16]TSS of TU
ssr2755YoeBssr2754PhdYeFM_antitox2074673 F21ChrAT-TTW-
N/KN/Kssr2784hypothetical (ChpI homolog)3474436 F26ChrATTWstand-alone
N/KN/Kssr2784aCOG2442/Duf433N/KN/KChrATTWstand-alone; potential asRNA regulation
ssl2921HigBssl2920hypotheticalN/KN/KChrAT-TTWTSS of TU
ssl2923PINssl2922VagCN/KN/KChrAT-T[3,16]see Figure 3
N/KN/Kssr3154HTH (COG2886)1211302 F99ChrATTWstand-alone; potential asRNA regulation
ssr3572YoeBssr3571PhdYeFM_antitox1620247 F29ChrAT-TTW-
ssr3589PINssr3588hypothetical1644850 F3314ChrAT-TTWTSS of TU
sll5003PINsll5004COG2442/Duf4331846 R20pSYSMAT-T[3,16]-
slr5017MNT slr5018HEPN (DUF86)23222 F33pSYSMT-AT[3,16]-
ssr5020HigBslr5021HTH-Xre24162 F19pSYSMT-AT[3,16]-
N/KN/Kssl5025COG2442/Duf43328079 R26pSYSMATTWstand-alone
sll5094PINssl5095PhdYeFM_antitox93540 F20pSYSMAT-T[3,16]-
slr5102DUF1814slr5101hypothetical96110 F14pSYSMAT-TTW-
slr5116PIN3ssr5117HicB106534 F571pSYSMT-AT[3,16]-
slr6049PINssr6048hypothetical47161 F21pSYSXAT-TTW-
slr6057COG2856/Duf955slr6056HTH-Xre54326 F17pSYSXAT-T[3,16]potential asRNA regulation; linked to possible bacterial multiubiquitin system; see Figure 4
slr6100RelE (DUF2136)slr6101HTH-Xre93754 F66pSYSXT-AT[3,16,22]-
sll7003PINssl7004plasmid stability2074 R42pSYSAAT-T[3,16,23]-
ssl7007PINsll7006HicB3798 R207pSYSAT-AT[3,16,23]-
ssr7017cas2slr7016cas1N/KN/KpSYSAT-AT[33]potential asRNA regulation
sll7030GNAT sll7031COG4453/DUF177829663 R (27264 R)24pSYSAAT-T[3,16]-
N/KN/Kslr7032HEPN (DUF86)30317 F (27918 F)35pSYSAATTWstand-alone
sll7033PINsll7034COG2442/Duf43331387 R (28988 R)18pSYSAAT-T[3,16,23]-
ssl7039RelEssl7038HTH-Xre36479 R (34080 R)0pSYSAT-AT[3,16,23]leader-less; see Figure 5 and Figure 6
slr7041MazFssr7040MazE36567 F (34168 F)28pSYSAAT-T[3,16,23]-
ssl7046PIN3N/KN/K40539 R (38140 R)25pSYSATTWstand-alone
N/KN/Kssl7048MazE (AbrB)N/KN/KpSYSAATTWstand-alone
ssr7072cas2slr7071cas164765 F (62366 F)57pSYSAAT-T[33]-
ssr7093cas2slr7092cas188315 F312pSYSAAT-T[33]-
sll8007PINssl8008COG2442/Duf4336626 R25pSYSGAT-T[3,16]-
sll8011PINsll8012hypothetical9650 R29pSYSGAT-T[3,16]-
slr8014PIN3ssr8013DUF22818658 F1055pSYSGAT-T[3,16]see Figure 3 and Figure 7
sll8027PIN3ssl8028MazE (AbrB)24116 R39pSYSGAT-T[3,16]-

Share and Cite

MDPI and ACS Style

Kopfmann, S.; Roesch, S.K.; Hess, W.R. Type II Toxin–Antitoxin Systems in the Unicellular Cyanobacterium Synechocystis sp. PCC 6803. Toxins 2016, 8, 228. https://doi.org/10.3390/toxins8070228

AMA Style

Kopfmann S, Roesch SK, Hess WR. Type II Toxin–Antitoxin Systems in the Unicellular Cyanobacterium Synechocystis sp. PCC 6803. Toxins. 2016; 8(7):228. https://doi.org/10.3390/toxins8070228

Chicago/Turabian Style

Kopfmann, Stefan, Stefanie K. Roesch, and Wolfgang R. Hess. 2016. "Type II Toxin–Antitoxin Systems in the Unicellular Cyanobacterium Synechocystis sp. PCC 6803" Toxins 8, no. 7: 228. https://doi.org/10.3390/toxins8070228

APA Style

Kopfmann, S., Roesch, S. K., & Hess, W. R. (2016). Type II Toxin–Antitoxin Systems in the Unicellular Cyanobacterium Synechocystis sp. PCC 6803. Toxins, 8(7), 228. https://doi.org/10.3390/toxins8070228

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop