RalR (a DNase) and RalA (a small RNA) form a type I toxin–antitoxin system in Escherichia coli

For toxin/antitoxin (TA) systems, no toxin has been identified that functions by cleaving DNA. Here, we demonstrate that RalR and RalA of the cryptic prophage rac form a type I TA pair in which the antitoxin RNA is a trans-encoded small RNA with 16 nucleotides of complementarity to the toxin mRNA. We suggest the newly discovered antitoxin gene be named ralA for RalR antitoxin. Toxin RalR functions as a non-specific endonuclease that cleaves methylated and unmethylated DNA. The RNA chaperone Hfq is required for RalA antitoxin activity and appears to stabilize RalA. Also, RalR/RalA is beneficial to the Escherichia coli host for responding to the antibiotic fosfomycin. Hence, our results indicate that cryptic prophage genes can be functionally divergent from their active phage counterparts after integration into the host genome.


INTRODUCTION
Toxin/antitoxin (TA) systems are widespread among prokaryotes (1). Five different types of TA systems have been characterized, depending on the interaction of the TA and the nature of the antitoxin (2,3). For type I systems, an RNA antitoxin interacts with the toxin transcript and inhibits translation of the toxic protein (4). The toxins and antitoxins of type II systems interact through direct proteinprotein binding (4). Type III systems rely upon the direct interaction of an RNA antitoxin with the toxin protein (5). A type IV designation has been proposed for a TA system in which the protein antitoxin does not interact with the toxin directly but suppresses the toxicity of the toxin by stabilizing its target (6), and a type V designation has been proposed in which the proteic antitoxin cleaves specifically the mRNA of the toxin to prevent the translation of the toxin (7,8).
In Escherichia coli, most well-studied TA systems belong to type II TA systems, in which the labile proteic antitoxin binds to the more stable toxin and inhibits its activity. Type I and type III systems are less well studied probably due to the delay in the identification of small RNAs (sRNAs) (9). However, in recent years, bioinformatic searches for putative type I TA systems based on tandem copies of the full loci and the presence of transmembrane domains has led to the identification of multiple copies of potential type I loci within new hosts across 774 bacterial genomes (10).
The role of TA systems in cell physiology, specifically in biofilm formation (11,12), persister cell formation (13,14), the general stress response (15,16), phage inhibition (17,18), and differential mRNA decay (19,20) is becoming more clear. Type II toxins with mRNA endoribonuclease activity usually alter gene regulation by cleaving specific mRNAs (21). Type I toxin TisB and type V toxin GhoT cause membrane damage when overproduced and are stress-response elements that are actively involved in persistence (7,8,22). Type I toxin Hok utilizes post-segregational killing to stabilize the R1 plasmid (23), but most of the Hok-like toxins in the E. coli K-12 chromosome appear inactive (24). Another major group of plasmid-based toxins function as gyrase inhibitors, and homologues of the plasmid RK2 ParE/ParD TA system in the Vibrio cholerae genome help maintain the integrity of the two chromosomes (25). These results indicate that chromosomal toxins are functionally divergent from plasmid loci.
Prophages or prophage remnants carrying toxic genes have been found to harbor TA systems, and five pairs have been reported in E. coli K12: RelE/RelB (26), RnlA/RnlB (27), YpjF/YfjZ (28), YkfI/YafW (28) and CbtA/YeeU (28). Furthermore, two interesting protein toxins have been described in cryptic prophage rac in E. coli K 12. KilR is a toxic peptide that inhibits cell division by inhibiting FtsZ (29). It has been suggested that RalR alleviates restriction modification possibly to protect the bacterial chromosome when recombination generates unmodified products by the same mechanism as Ral of phage lambda (30,31). However, the amino acid identity between the Ral protein in phage lambda and the RalR protein in rac prophage is very low (24%) (30) and the function of RalR remains unclear. Lambdoid prophage rac has lost about 60% of its original DNA (32), and ralR is differentially regulated in the development of E. coli biofilms (33). Also, ydaC, downstream of ralR, has been named as rcbA for its ability to reduce the frequency of double-strand chromosome breaks (34). However, it remains unclear whether rcbA encodes RNA or protein (34), or whether ydaC is transcribed.
In this study, unexpectedly, we found that RalR does not inhibit restriction modification but instead is part of a toxin/antitoxin system and functions as a non-specific DNase. Furthermore, we found that the adjacent gene product, a trans-encoded sRNA RalA (for RalR antitoxin), functions as an antitoxin for RalR. The ralR and ralA genes are adjacent but in the opposite orientation, and RalA RNA has 16 nucleotides of complementarity to the coding region of RalR mRNA. We show that RalA RNA interacts with the mRNA of RalR via base-pairing, thus preventing the translation of RalR. The activity of antitoxin RalA requires RNA chaperone Hfq. Thus, RalR/RalA belongs to a type I TA system where the antitoxin sRNA interferes with translation of the toxin mRNA via the 16-nt base-pairing.

Bacterial strains, plasmids and growth conditions
The isogenic E. coli BW25113 K-12 strains and plasmids used in this study are listed in Table 1. For deleting and overexpressing single genes, the Keio collection (35) and the ASKA library (36) were used. The ralRA mutant strain was created using the Red method (37) using primers DralRA-f and DralRA-r (Supplementary Table S1). The kanamycin resistance cassette from the Keio and ralRA mutants was removed as previously described to ensure that only the impact of the deleted gene was studied (38). The ralR and ralRA mutations were verified by amplifying chromosomal DNA using the primers CralRA-f and CralRA-r, or CralRA-f2 and CralRA-r2, respectively, and sequencing the resulting fragments using CralRA-f (Supplementary Table S1) to also verify that genes upstream and downstream of the mutation (recT and ralR or ydaQ, respectively) were left intact. Removal of hfq was verified using the same method using primer pairs hfq-up and hfqr, hfq-up and hfq-down, hfq-f and hfq-r, or hfq-f and hfqdown as described previously (35). All experiments were conducted at 37 • C in Luria-Bertani medium (LB) (39) unless specified otherwise. Chloramphenicol (30 g/ml) was used for maintaining pCA24N-based plasmids, and ampicillin (50 g/ml) was used to maintain plasmid pBAD-ralA. nullnull Cloning of pCA24N-ralR-ralA, pCA24N-ralR-ralA-G289del and pBAD-ralA The full-length gene region ralR-ralA and partial gene region ralR-ralA-G289 were cloned as described previously (36). They were polymerase chain reaction (PCR) amplified from E. coli K12 BW25113 chromosomal DNA using one front primer (p-ralR-f) and two different rear primers (p-ralA-r and p-ralR-ralA-G289del-r) (Supplementary Table  S1). PCR products were phosphorylated using a kination kit (Takara, Dalian, China) and purified using a PCR product purification kit (Qiagen, Hilden, Germany). The purified products were then ligated into vector pCA24N that was digested with the StuI restriction enzyme (New England Biolabs, Beverly, MA, USA). The ligation mixture was transformed into BW25113 and the rac mutant. The constructs were verified by DNA sequencing using primer pair of p-f and p-r (Supplementary Table S1). Plasmid pBAD-ralA was constructed with the full-length gene ralA, and 124 bp upstream of ralA was controlled by the arabinose promoter using primer pair pBAD-ralA-f and -r (Supplementary Table S1) and cloning into NcoI-HindIII-digested pBAD/Myc-HisA (Invitrogen, Carlsbad, CA, USA).

Mutant library creation by epPCR
Error-prone PCR (epPCR) was conducted on plasmid pCA24N-ralR and pCA24N-ralR-ralA using primers epPCR-f and epPCR-r (Supplementary Table S1) as described earlier (40). The epPCR program was as follows: 94 • C for 5 min, 30 cycles of 1 min at 94 • C, 1 min at 55 • C, and 2 min at 72 • C, followed by 10 min at 72 • C final extension. The PCR products were gel-purified and digested using BseRI and HindIII before being ligated into pCA24N-ralR-ralA. The ligation mixture was transformed into BW25113 rac.

Single and multiple site-directed mutagenesis
Single site-directed mutagenesis (16) was used to mutate the putative start codon of ydaC from ATG to ACG using primer pair p-ralR-T189C-ralA-f and -r, and to introduce a stop codon TAA into the ninth putative coding region of ydaC in pCA24N-ralR-ralA using primer pair p-ralR-ralA-C161A-f and -r (Supplementary Table S1). For the 16-nt repeat in RalA, positions 63 (T to C) and 64 (G to C) were mutated using corresponding primer pairs. In addition, the AGC site in ralA was mutated to AAC using p-ralR-ralA-G83A-f and -r. Site-directed mutagenesis at multiple sites was applied to mutate the GAGC sequence into CAAC using p-ralR-ralA-G81CG83A-f and -r. Moreover, p-ralR-ralA-M1-f and -r were used to mutate the 16nt repeat sequence of ralA on pCA24N-ralR-ralA without changing the potential amino acid sequence. p-ralR-M1-f and -r were used to mutate the 16-nt repeat in ralR with no amino acid changes in pCA24N-ralR and pCA24N-ralR-ralA, and p-ralR-M2-f and -r were used to mutate the 16nt repeat in ralR in pCA24N-ralR and pCA24N-ralR-ralA with amino acid changes ( Table 1). The correct mutations were verified by DNA sequencing. This study pCA24N-ralR-ralA-G81CG83A G to C at position 81 in ralA; G to A at position 83 in ralA pCA24N-ralR-ralA-G83A G to A at position 83 in ralA pCA24N-ralR-ralA-G289del Sequences removed after G at position 289 in ralR-ralA pCA24N-ralR-T189C-ralA T to C at position 189 in ralR, changed putative ydaC start aa from Met to Thr, no change in RalR aa sequence pCA24N-ralR-ralA-C161A C to A at position 161 in ralA, introduced a stop codon TAA at the ninth aa of putative YdaC pCA24N-ralR-ralA-G64C G to C at position 64 in ralA pCA24N-ralR-ralA-T63C T to C at position 63, no change in aa sequence pCA24N-ralR-ralA-M1 Modified the 16-nt repeat in ralA to AGCCTCCTTCTTACCT, no change in aa of putative YdaC pCA24N-ralR-M2-ralA Modified the 16-nt repeat in ralR, change aa from Gly-Ser-Glu-Lys-Glu-Ala to Gly-Ser-Ala-Ile-Ala-Glu pCA24N-ralR-M1-ralA-M1 Modified the 16-nt repeat both in ralR and ralA to AAGTGAGAAGGAGGCT and AGCCTCCTTCTTACCT Mutated pCA24N-ralR plasmids This study pCA24N-ralR-A154G Encodes RalR-mutant RalR K52G pCA24N-ralR-M1 Modified the 16-nt repeat in ralR to AAGTGAGAAGGAGGCT, no change in aa pCA24N-ralR-M2 Modified the 16-nt repeat in ralR, change aa from Gly-Ser-Glu-Lys-Glu-Ala to Gly-Ser-Ala-Ile-Ala-Glu Note. Cm R and Amp R indicate chloramphenicol and ampicillin resistance, respectively, and aa indicates amino acid.

Protein expression and purification
RalR, RalR-mutant, and Hfq each with six histidines and 10 additional cloning artifacts at the N-terminus were purified via BW25113 with pCA24N-ralR, pCA24N-ralR-A154G (encodes inactive RalR K52Q) and pCA24N-hfq (36). Strains were grown in LB with chloramphenicol and were induced with 1 mM IPTG at a turbidity of 0.1 for 5 h. Cells were collected and resuspended in 25 ml lysis buffer [50 mM potassium phosphate buffer (pH 8.0), 300 mM NaCl, and protease inhibitor cocktail (Sigma-Aldrich, USA)]. Then samples were sonicated using a Sonic Dismembrator (Ningbodongzhi, China) at level 2 for 5 min twice with 2 s sonication and 4 s break while being cooled in ice. Ni-NTA resin (Qiagen) was used according to the manufacturer's protocol. Purified proteins were desalted by passage on disposable Sephadex G-25 pre-packed PD-10 columns pre-equilibrated in 20 mM Tris-HCl buffer (pH 8.0), and the protein concentration was measured by using a Bi Yuntian BCA assay kit (Haimen, China).

DNA cleavage assay
The DNA cleavage assay was performed as described earlier (41,42) with some modifications. Purified RalR and RalR-mutant (150 pmol each) were incubated separately with different substrates at 37 • C in 250 mM NaCl, 100 mM Tris-HCl (pH 7.5), 10 mM CaCl 2 and 10 mM MgCl 2 for 30 min and 120 min. Plasmid pBR322 isolated from BW25113, genomic DNA of BW25113, as well as methylated and unmethylated lambda DNA (dam − , dcm − ) (1 g) were used as substrates, respectively. The reaction was stopped by the addition of a stop solution (25% glycerol, 0.5% SDS, 0.05% bromophenol blue and 50 mM EDTA) and was analyzed by electrophoresis on 1% agarose gels stained with SYBR safe (Invitrogen). The equal amount of deoxyribonuclease I (DNase I; New England Biolabs, Beverly, MA, USA) was used as a positive control, and excess EDTA (20 mM) was added in the reaction to inhibit the nuclease activity. To study whether Ca 2+ and Mg 2+ are cofactors for RalR, 250 mM NaCl and 100 mM Tris-HCl (pH 7.5) buffer with or without the addition of 10 mM Ca 2+ or 10 mM Mg 2+ were used to conduct the DNA cleavage assay. The RalR-mutant was used as a negative control.

In vitro transcription and RNA labeling
For in vitro synthesis of T7-RalA sRNA and T7-ompA mRNA, PCR products were obtained from genomic DNA using primer pairs T7-ompA-f and -r and T7-ralA-f and -r,  Table S1). The PCR-amplified products were used as templates for in vitro transcription with T7 RNA polymerase. The T7 RNA polymerase promoter sequence was included in the forward primers. About 1 g gel-purified PCR products were used as the templates for in vitro RNA reactions with the T7 High Yield RNA Synthesis Kit (New England Biolabs, Beverly, MA, USA). About 50 pmol of the in vitro transcribed RNA was labeled using the Pierce TM RNA 3 end biotinylation kit (Thermo Scientific, Hudson, NH, USA) according to the manufacturer's protocol.

RNA cleavage assay
The RNA cleavage assay was performed as described earlier (7). The reaction mixture for the endoribonuclease cleavage assay (10 l) contained 2 g in vitro synthesized T7-ompA mRNA, 50 mM Tris-HCl (pH 8.5), 100 mM KCl, 2.5 mM MgCl 2 and 30 g of purified RalR or RalR-mutant protein. The reaction mixture was incubated at 37 • C for 10, 20 and 30 min and quenched by the addition of an equal volume of 2 × TBE-urea sample buffer (Invitrogen). Inactive RalR-mutant was used as a negative control. The reaction products were resolved by 15% TBE-urea gels. RNA was stained with SYBR safe.

Tricine-SDS-PAGE and western blot analysis
To investigate RalR protein levels, Tricine-SDS-PAGE and a western blot were performed. BW25113 rac strains containing pCA24N, pCA24N-ralR, pCA24N-ralA and pCA24N-ralR-ralA were grown to a turbidity of 0.2 in LB with chloramphenicol, then 1 mM IPTG was added to produce RalR for 5 h, and cells were washed with TE buffer. Samples were sonicated and total protein was quantified as mentioned above. Protein was denatured at 95 • C for 5 min. The Tricine-SDS-PAGE was performed as described earlier (43). Total protein (25 g) of each sample was loaded for SDS-PAGE, and a western blot was performed using 2.5 g protein of each sample with primary antibodies raised against a His-tag (Cell Signaling Technology, Danvers, MA, USA) and horseradish-peroxidase-conjugated goat antimouse secondary antibodies (Bio-Rad, Richmond, CA, USA).

RNA-protein electrophoretic mobility shift assay (EMSA)
Biotinylated RalA sRNA (4 × 10 −3 pmol) was mixed with different concentrations of Hfq protein in 20 l buffer containing 10 mM HEPES (pH 7.3), 20 mM KCl, 1 mM MgCl 2 , 1 mM DTT and 5% glycerol at 37 • C for 30 min. Samples from the reactions were loaded onto a pre-run 6% native polyacrylamide gel, and electrophoretic transfer was used to bind the proteins to nylon membranes (400 mA for 30 min) after the separation. After crosslinking for 5 min using UV-light, biotin-labeled RNA was detected using a LightShift R Chemiluminescent RNA EMSA kit (Thermo Scientific, Hudson, NH, USA) according to the manufacturer's protocol. In this assay, in vitro synthesized ompA mRNA was used as a negative RNA control and inactive Hfq protein was used as a negative protein control.

Quantitative real-time reverse-transcription PCR (qRT-PCR)
Total RNAs were isolated as described previously (11) and was used as the template for the qRT-PCR reaction using the SuperScript TM III Platinum SYBR R Green One-Step qRT-PCR Kit (Invitrogen). Primers for qRT-PCR of ralR-, ralA-and ralA-related fragments are listed in Supplementary Table S1. The level of rrsG transcript was used to normalize the gene expression data.

Primer extension
5 end FAM dye (6-carboxyfluorescein) labeled primers FAM-ralA-r2 and FAM-ralA-r3 (Supplementary Table S1) were ordered from Invitrogen. A total of 30 g RNA was added to 2 × 10 −4 pmol of 5 end labeled primer, and the mixture was added to 3 l of 10× first-strand and 37.5 U AMV reverse transcriptase (Promega, Madison, WI, USA). The RNA mix was annealed to the primers by incubating at 37 • C for 1 h. To enhance the concentrations of FAMlabeled reverse transcriptions, the incubated mixture was retranscribed by adding 6 l of 10× first-strand buffer and 75U AMV reverse transcriptase and re-incubated at 37 • C for another 1 h. The products were screened by an ABI3730 DNA Analyzer (Applied Biosystems, Foster City, California, USA) and the results were analyzed using Genemapper (Version 4.1). In this study, only RNA isolated from BW25113 rac/pCA24N-ralR-ralA induced with 1 mM IPTG for 20 min was used due to the low expression level of ralA in BW25113.

Prediction of the ralA Rho-independent terminator
To determine whether RalA sRNA contains a Rhoindependent structure, we folded its sequence with RNAS-TRUCTURE 4.6 (44). A Rho-independent terminator was predicted based on properties characterized previously (45).

Fosfomycin resistance assay
Metabolic activities in the presence of fosfomycin and glucose 6-phosphate (to enhance the activity of fosfomycin) were measured using reagents from Biolog, Inc. (Hayward, CA, USA). Cells were grown to OD 600 1.0, then diluted to a turbidity of 0.07 in IF-10a (Cat. No. 72264), and further diluted 200-fold into a solution containing IF-10a, BioLog Redox Dye D (Cat. No. 74224), and rich medium (2.0 g/l tryptone, 1.0 g/l yeast extract and 1.0 g/l NaCl) to a final OD 600 of 0.0003. Volumes of 100 l of this cell suspension (with 50 mg/ml glucose 6-phosphate and 0 or 0.2 g/ml fosfomycin) were transferred into 96-well microtiter plates and incubated at 37 • C. Metabolic activity was measured via the absorbance at 590 nm, which indicates the intracellular reduced state due to formazan (purple) formed from a tetrazolium dye. Growth in the presence of fosfomycin was determined in minimal medium; cells were washed once with M9 minimal medium with 0.4% glucose (M9-Glu) (46), then diluted to OD 600 0.01 in M9-Glu with or without 2 g/ml fosfomycin. Turbidities were measured at 600 nm.

RalR is a proteic toxin with endonuclease activity
RASTA-bacteria has predicted RalR of prophage rac to be a small toxin (64 aa) ( Figure 1A, Supplementary Figure S1) (47). Indeed, induction of ralR via pCA24N-ralR in E. coli results in growth inhibition ( Figure 1B), along with a 10 3fold reduction in colony-forming units (CFU) per milliliter in the wild-type host ( Figure 1C). For the strain with the complete rac prophage of 23 kb removed by natural excision, rac, which is devoid of all rac prophage genes (38), overproducing RalR was also very toxic ( Figure 1B and C). Hence, RalR is a potent toxin.
We investigated the cause of the RalR toxicity and found it was due to DNase activity. As shown by an in vitro DNA cleavage assay, purified RalR cleaves lambda DNA in a pattern similar to DNase I (Figure 2A). DNase I is an endonuclease that non-specifically cleaves DNA to release di-, tri-and oligonucleotides with 5 -phosphorylated and 3hydroxylated ends (48). After 30 min, ∼70% of the lambda DNA was cleaved by RalR, and all of the lambda DNA was cleaved by DNase I. After 120 min, all of the lambda DNA was cleaved by RalR (Figure 2A and B). Moreover, RalR also cleaved native genomic DNA of E. coli K-12 BW25113 and plasmid DNA pBR3222 (Supplementary Figure S2A).
To test whether the endonuclease activity of RalR depends on modification of substrates by methylation, unmethylated lambda DNA was also tested, and the results show that RalR cleaves unmethylated DNA equally well, indicating that methylation is not required ( Figure 2B). Moreover, Ca 2+ and/or Mg 2+ were essential for RalR endonuclease activity ( Figure 2C), and as expected, the addition of excess EDTA inhibited the activity of both RalR and DNase I (Figure 2A and B).
Some nucleases degrade all nucleic acids, i.e. both RNA and DNA (49). Thus, the endoribonuclease activity of RalR was also tested, but no activity was detected in the in vitro RNA cleavage assay (Supplementary Figure S2B). Therefore, RalR is a non-specific endonuclease which cleaves DNA but not RNA. Furthermore, we found that overproduction of RalR after 5 h leads to filamentous growth as a result of the SOS response induced by degrading DNA (Supplementary Figure S2C).
We further determined the key residues in RalR using ep-PCR using pCA24N-ralR as the template. After screening ∼400 colonies for toxicity, 39 of them with reduced toxicity were sequenced, and changes at amino acid position 14, 37, 39, 49, 51 and 52 reduced RalR toxicity (Supplementary Figure S3). One of the mutated proteins (at aa 52) was also purified (RalR-mutant, RalR K52Q) ( Supplementary Figure S4), and since it did not cleave DNA under the same conditions (Figure 2A and B), the DNase activity of RalR is not due to contamination.

DNA associated with ydaC functions as an antitoxin
Next, we tested whether previously annotated YdaC works as the antitoxin for RalR. YdaC is annotated as a racencoded protein (69 aa), and the start codon of ydaC overlaps with the stop codon of ralR (http://ecogene.org) (Supplementary Figure S1). The chromosomal region from the translational start of ralR to the translational stop of ydaC was used to construct a new plasmid pCA24N-ralR-ralA ( Figure 1A). Results showed that when co-expressed with ralR, ralA, via the pCA24N-ralR-ralA construct, reduces the toxicity of RalR, demonstrating that the DNA associated with ydaC is necessary for repression of RalR toxicity ( Figure 1B and C). However, no ydaC-encoded protein was detectable via SDS-PAGE or western blot when ydaC was tagged with an N-terminal hexahistidine sequence ( Figure  3A).
To further test whether YdaC protein is the antitoxin, we introduced a stop codon by a single nucleotide change into ydaC at the putative amino acid position 9 (nucleotide position 212 from G to T) (pCA24N-ralR-ralA-C161A) and tested its impact on cell growth. However, the G212T substitution did not affect the ability of the ydaC-associated DNA to block the toxicity of RalR in cell growth ( Figure 3B). In addition, when the putative start codon ATG at position 189 relative to the first nucleotide of start codon of ralR (lies inside the coding region of RalR) was mutated into a non-start codon ACG (pCA24N-ralR-T189C-ralA), which would prohibit the translation of the hypothetical YdaC, the DNA associated with ydaC still functioned as an antitoxin ( Figure 3B). Thus, the DNA associated with ydaC inhibits the toxicity of RalR, but ydaC does not encode a protein.
When we aligned ralR with the region associated with ydaC, we found two 16-nt repeats inside ralR and inside ydaC, respectively ( Figure 1A). One of the 16-nt repeats lies in the coding region of the ralR, and the other, in the opposite orientation, lies 105-nt downstream of the end of the coding region of ralR on the complementary strand. In light of the recent discoveries that non-coding regulatory RNAs could act as antitoxins by base-pairing with the toxin mRNAs using limited complementarity (10-25 nt) (10,50-52), we hypothesized that the DNA strand complementary to ydaC encodes an RNA antitoxin for RalR and propose the name RalA for RalR antitoxin.

RalA is a trans-encoded sRNA
To further study the antitoxin activity of RalA, we constructed a dual plasmid system to allow expression of ralR and ralA independently. Plasmid pCA24N-ralR was induced by IPTG to express ralR and pBAD-ralA was induced by L-arabinose to express ralA using the complementary strand of the ydaC gene (from 99-nt downstream of the stop codon to the start codon of ydaC) ( Figure 1A). pBAD-ralA was constructed to have a stop codon (TAA) immediately after the ATG start codon so that only RNA was formed. Toxicity test results indicated that there was a 100fold reduction in viability when ralR is induced by IPTG for 4 h; however, when L-arabinose was also added to induce ralA, the toxicity of RalR was completely inhibited ( Figure  4A). These results show that the RalA sRNA provides the antitoxin activity for RalR.
To map the 5 end of RalA, we carried out primer extension experiments with oligonucleotide FAM-ralA-r2 with complementarity to RalA at regions 34-nt downstream from the 16-nt repeat region ( Figure 4B). The results revealed a major extension product that is 110 nt in size, which indicates that the start of the transcript is 62-nt upstream of Since non-coding sRNAs frequently terminate transcription via Rho-independent transcription terminators (53), we first searched for a potential Rho-independent termination signal in ralA. After folding RalA RNA by RNAS-TRUCTURE 4.6 (44), we identified a Rho-independent termination signal, with a stem-loop structure composed of 5 base pairs (GC rich) and a loop of 14 bases followed by four stretches of uridine residues ( Figure 4D). This stretch is likely a Rho-independent terminator since (i) it is predicted to form a stem-loop structure, (ii) the formed stem has 5-10 base pairs in length and is GC rich with at least 60% GC base pairs, (iii) the formed loop generally contains 3-8 bases, (iv) in most structures, the average free energy calculated for a Rho-independent terminator is around −7 kcal/mol and here it was found to be −6.8 kcal/mol, and (v) the stem structure usually is followed by a stretch of Us (the average number is 4) and here it was found to have 4 (45). The Rho-independent termination signal thus gives an RNA product of 179 nucleotides ( Figure 4D), which indicates that RalA sRNA is wholly contained within what was previously annotated as ydaC (on the opposite DNA strand). This explains why the ydaC DNA fragment was found to reduce RalR toxicity ( Figure 1B and C). Moreover, the 5 end of RalA also formed a stable stem-loop at the transcription start site determined by primer extension, which helps to stabilize the sRNA molecule ( Figure 4D). Furthermore, both transcripts were detected in the exponential phase (OD 600 0.8) and stationary phase ( OD 600 6.2) of the wild-type strain using qRT-PCR with primer pair ralR-f and ralR-r (for ralR) and primer pair ralA-f2 and ralA-r3 (for ralA) (Supplementary Table S1), which shows that both native promoters function in vivo.
We further confirmed the length of the RalA transcripts using qRT-PCR and using cDNA templates synthesized from strain-specific primers in exponentially growing cells of BW25113 and BW25113 rac carrying pCA24N-ralR-ralA ( Figure 4B). To test further whether RalA is an untranslated RNA, we searched for open reading frames in which a putative initiation codon (ATG, TTG or GTG) was preceded at an appropriate distance by a ribosome-binding  site. This transcript does not contain a potential initiation codon preceded by a potential ribosome-binding site, thus it is unlikely that RalA encodes a protein.

RalA functions through base-pairing with RalR
We next investigated whether RalA RNA functions through a base-pairing mechanism. epPCR was employed to screen for the mutations in ralR and ralA that result in reduced antitoxin activity using pCA24N-ralR-ralA as template. After screening ∼400 colonies for toxicity, 20 of them with toxicity were sequenced, and the results showed that half of them had mutations within the two 16-nt repeats ( Figure 5A) in the ralR gene and in the ralA gene. These results suggest that the 16-nt repeat is important for RalA RNA to inhibit RalR mRNA.
To further investigate the involvement of the 16 nt in the base-pairing of the toxin and antitoxin RNAs, we performed site-directed mutagenesis on the 16-nt repeats as well as other neighboring regions. Our results showed that single mutations in the 16-nt repeat of ralA reduced its antitoxin function (pCA24N-ralR-ralA-G64C and pCA24N-ralR-ralA-T63C), and four mutations in the 16-nt repeat of ralA completely abolished the antitoxin activity of RalA (pCA24N-ralR-ralA-M1). In contrast, a single point mutation (pCA24N-ralR-ralA-G83A) and two point mutations (pCA24N-ralR-ralA-G81CG83A) that changed one or two nucleotides in regions outside of the 16-nt repeat of ralA did not affect the antitoxin activity of the RalA RNA ( Figure  5B). In addition, to determine the minimal length of RalA RNA required for maintaining the function as an antitoxin, we removed all of the sequence after position 289, which eliminated the 16-nt repeat in RalA (pCA24N-ralR-ralA-G289del). As a result, the antitoxin function of RalA was completely abolished ( Figure 5B). These results confirm the importance of the 16-nt repeat in the RalR/RalA TA pair in regard to the interaction of the toxin and antitoxin.
Moreover, when the 16 nt in ralR was altered to cause a synonymous mutation and non-synonymous mutation using pCA24N-ralR, the synonymous mutation retained the RalR toxicity (pCA24N-ralR-M1) while the nonsynonymous mutation led to a RalR variant (pCA24N-ralR-M2, pCA24N-ralR-M2-ralA) that was no longer toxic ( Figure 5A). These results confirm that RalR functions as a protein and, moreover, that amino acid changes at the 16-nt repeat coding region affect RalR toxicity.
To further test whether changes in the ralA 16-nt repeat can neutralize the toxic effect of RalR when the corresponding mutation is encoded by ralR, we introduced two complementary mutations and constructed pCA24N-ralR-M1-ralA-M1, which contains two identically mutated 16-nt repeats. The results show that the mutated ralA retains its antitoxin function when the mutation matches that of ralR, in contrast to the result that mutated ralA could not attenuate the toxic effect of the unmutated ralR (pCA24N-ralR-ralA-M1) ( Figure 5A). Therefore, we provide several lines of evidence that RalA is a novel RNA antitoxin in which its activity to block the toxicity of RalR relies on a 16-nt repeat.

RalA likely inhibits RalR mRNA translation
To determine how RalA RNA attenuates the toxic effects of RalR, we checked the RalR protein levels with and with- out RalA. Through SDS-PAGE and western blot analysis, we found that the protein levels of RalR are reduced in the presence of RalA ( Figure 3A), which indicates that RalA interferes with the accumulation of RalR protein. Basepairing between sRNA and its target mRNA usually leads to repression of protein levels through translational inhibition and/or mRNA degradation (51). Our qRT-PCR results showed that in the presence of more RalA, the levels of RalR mRNA were not reduced substantially (Table 2), thus it is unlikely that RalA RNA promotes RalR mRNA degradation. Translational inhibition of sRNA with a target gene often occurs when sRNA base-pairs with the 5 -UTR of target mRNA; however, other locations for base-pairing and consequent mechanisms of regulation are possible (51).

Hfq is required for RalA activity
To facilitate limited base-pairing for the trans-encoded sR-NAs and the target mRNAs, in many cases the RNA chaperone Hfq is required for trans-encoded sRNA-mediated regulation (51). The antitoxin RalA was tested in an hfq deletion mutant strain and, as expected, the hfq mutation abolished antitoxin activity of RalA ( Figure 6A and B). Hence, Hfq is required for RalA RNA activity.
To provide additional evidence for the requirement of Hfq for RalA RNA activity, we investigated direct binding of Hfq to RalA RNA in vitro. Hfq interacts strongly with small non-coding RNAs at single-stranded AU-rich regions (54), and RalA is AU rich (63.1%) and has four stretches of more than seven continuous U/As. RNA-protein binding was performed, and Hfq was found to bind RalA RNA ( Figure 6C). In contrast, Hfq did not bind to the negative control (3 end of the coding region of OmpA mRNA, Figure 6C). Moreover, it seems that Hfq helps to stabilize the RalA transcript as RalA RNA was degraded faster in the absence of Hfq as shown by qRT-PCR, while both RalR and RalA RNAs appeared stable in the wild-type cells with half-times of more than 20 min (Table 3).

RalR increases resistance to fosfomycin
Large-scale phenotypic screening of all Keio mutants from PortEco showed that the ralR knockout strain has a negative fitness score when exposed to fosfomycin (55). Fosfomycin is a broad-spectrum antibacterial agent that targets bacterial mucopeptide synthesis by inhibiting phosphoenolpyruvate transferase, the first enzyme involved in the synthesis of peptidoglycan (56). Our results show that both the ralR and the ralRA strain are more sensitive to fosfomycin (0.2 g/ml) compared to the wild-type strain based on their reduced metabolic activity ( Figure 7A and B). Corroborating these results, the ralRA strain is more sensitive to fosfomycin (2 g/ml) compared to the wild-type in terms of growth as well ( Figure 7C and D). Hence, the RalR/RalA TA locus has an impact on cell physiology and is beneficial for resistance to an antibiotic.

DISCUSSION
Collectively, our results support that RalR/RalA forms a type I TA pair. These results are: (i) both genes are small and are adjacent, (ii) RalR functions as a toxin that inhibits growth, (iii) RalA functions as an antitoxin and is a non-coding RNA, (iv) RalA blocks RalR-mediated toxicity by the base-pairing of 16 key nucleotides and prevents Exponentially growing cells (OD 600 0.8) were induced with 1% L-arabinose (+ ara) or without L-arabinose (− ara) for 20 min and 8 h. Lower Ct values indicate higher expression levels, and rrsG was used to normalize total RNA levels. Two independent cultures were used for the assay. Fold changes in the transcription of various targets with or without L-arabinose were calculated as 2 ∧ −(C ttarget(+ ara) −C t rrsG(+ ara) )/2 ∧ −(C ttarget(− ara) −C trrsG(− ara) ). Means and standard errors are indicated. the translation of RalR mRNA, and (v) Hfq is required for RalA antitoxin activity. As a novel type I TA system, RalR is the first toxin that functions as an endonuclease in E. coli. A recent study showed that when compared with a wild-type strain, a ydaC deletion mutant has an increased level of double-strand breaks that is remedied by expression of ydaC from a plasmid, suggesting that ydaC helps to maintain the integrity of the bacterial chromosome (34). These results fit well with our results here in that deletion of ralA (on the opposite strand as ydaC) leads to activation of toxin RalR that should result in double-stranded DNA breaks. Although overproduction of RalR leads to filamentous growth due to RalR-mediated DNA damage that activates the SOS response, it is not clear whether RalR can trigger cell death under physiological conditions. Another novel feature of this TA system is that RalA is a trans-encoded sRNA antitoxin that neutralizes the protein toxin by base-pairing within the toxin mRNA coding region. Traditional type I toxin and antitoxin genes are usually located at the same locus, and sRNA antitoxins can have greater than 60-nt base-pairing with the 5 UTR (e.g. Ldr of the Ldr/Rdl family in E. coli (57)) or the 3 UTR Overnight cultures were diluted to OD 600 0.1, and re-grown in LB till OD 600 1. Rifampicin (50 g/ml) was added and samples were taken after 5 min, 10 min and 20 min. Fold changes in the concentrations of the targets at different times points were calculated as 2 ∧ −(C ttarget T = 5 or 10 or 20 −C trrsG T = 5 or 10 or 20 )/2 ∧ −(C ttarget T = 0 −C trrsG T = 0 ). Values less than one are indicated as negative fold changes (i.e. the amount of RNA that is reduced). Means and standard errors are indicated. Figure 7. The RalR/RalA TA system increases resistance to fosfomycin. (A) Metabolic activity of the wild-type, ralR and the ralRA strains subjected to 0.2 g/ml fosfomycin with 50 g/ml glucose-6-phosphate to enhance fosfomycin activity and (B) metabolic activity with 50 g/ml glucose-6-phosphate alone (negative control). This assay was performed twice with two independent cultures, and one standard deviation is shown. Growth of the wild-type and ralRA strains (C) with 2 g/ml fosfomycin and (D) without fosfomycin. Two independent cultures were evaluated, and one standard deviation is shown.
of the toxin mRNA (e.g. TxpA/RatA in Bacillus subtilis (58)), and the two transcripts can also overlap (e.g. Ibs/Sib in E. coli) (50,52). Recently, four unconventional type I TA  systems have been reported where the antitoxins are transcribed divergently from the toxin gene and have limited  complementation regions to the toxin RNA (often 18-21 nucleotides), including TisB/IstR-1 (59), ShoB/OhsC (50) and ZorO/OrzO loci (50). Here with RalR/RalA, we provided evidence that RalA reduces RalR toxicity via the complementarity of 16 nt. For the trans-encoded base-pairing sRNAs, the potential base-pairing with target mRNA typically is 10-25 nt. Studies on ZorO/OrzO showed that for the 18-nt complementarity region, not all of the base-pairing interaction is necessary, and it is dependent on whether the complementary region is contiguous or not (52). Our results confirmed this and showed that 16 nt of complementarity (with 11 nt of contiguous perfect complement) is sufficient for base-pairing interaction between toxin and antitoxin. Cryptic prophage rac is a DNA fossil since it was acquired over 4.5 million years ago (60), yet it has become a stable resident in the E. coli chromosome. Ral in active phage lambda and RalR in prophage rac share a low sequence identity of 24%, indicating that the two proteins are not related. In contrast, a high identity of 77% exists for P22 Ral and lambda Ral with highly conserved Nand C-termini (Supplementary Figure S5). Moreover, an even higher identity of 94% is found for recombinant Enterobacteria phage phi21 Ral and lambda Ral with highly conserved N-and less conserved C-termini (Supplementary Figure S5). Therefore, Ral is more conserved in phages while RalR is more conserved in E. coli strains. Ral of lambda phage acts by modulating the restriction and modification (RM) activities of the type I restriction systems in E. coli (61), and RM systems and TA systems are related (3). In general, they originate from phages and plasmids and are used to invade hosts or to keep them in the host by post-segregational killing. However, after integration into the chromosome, the function of these two systems has diverged. We found that RalR in E. coli functions differently, acting as a toxin by cleaving DNA, and that it belongs to a type I toxin-antitoxin system. We show here that removal of ralRA reduces the resistance of the cell to the antibiotic fosfomycin; hence, there is a clear benefit to the cell for harboring the RalR/RalA TA system. Therefore, our results provide further evidence that the genes of cryptic prophage impact cell physiology and that these genes may be used to increase resistance to stress (38).
Besides RalR/RalA, five other TA systems located in the E. coli K12 cryptic prophages have been reported, including the well-studied RelE/RelB (in Qin prophage) (26) and four recently identified TA pairs YpjF/YfjZ (in CP4-57) (28), RnlA/RnlB (in CP4-57) (27), YkfI/YafW (in CP4-6) (28) and CbtA/YeeU (in CP4-44) (6). These TA systems are located in prophages that are not inducible by SOS responses (38), suggesting a relatively stable residence in the host genome. Type I TA pair TxpA/RatA in B. subtilis is also in a phage-like element, and the phage-like element is retained in the host chromosome for spore formation (58). TA systems have been linked to persistence, biofilm formation and the general stress responses, and ralR was found to be highly induced in the later stage of E. coli biofilm de-velopment (24 h biofilms compared to 4 h biofilms) (33). Indubitably, other uncharacterized TA systems in cryptic prophages might also contribute to the stress response of the host, and there are probably additional intricate regulation mechanisms that need to be investigated.
Unconventional type I TA systems are more limited in distribution than typical type I TA systems, and are found mainly in Escherichia, Shigella and Salmonella species. By comparing the phylogenetic tree of those identified type I loci with the host taxonomy, earlier studies suggest that type I loci have not been freely disseminated by horizontal gene transfer but instead may have a common ancient ancestor (10). This agrees with what we found for TA pair RalR/RalA in rac prophage. The discovery of a high occurrence of TA systems in E. coli prophages will help to search for new TA systems in other prokaryotes. Since TA genes from bacterial prophages are actively expressed in the human gastrointestinal tract, TA systems may play a dynamic role in more ecosystems than expected (62). The mechanism of harboring TA systems in prophages, the origins of these TA systems, and the regulation of these loci in response to various stress conditions need to be investigated further.

SUPPLEMENTARY DATA
Supplementary Data are available at NAR Online.