MqsR/MqsA Toxin/Antitoxin System Regulates Persistence and Biofilm Formation in Pseudomonas putida KT2440

Bacterial toxin/antitoxin (TA) systems have received increasing attention due to their prevalence, diverse structures, and important physiological functions. In this study, we identified and characterized a type II TA system in a soil bacterium Pseudomonas putida KT2440. This TA system belongs to the MqsR/MqsA family. We found that PP_4205 (MqsR) greatly inhibits cell growth in P. putida KT2440 and Escherichia coli, the antitoxin PP_4204 (MqsA) neutralizes the toxicity of the toxin MqsR, and the two genes encoding them are co-transcribed. MqsR and MqsA interact with each other directly in vivo and MqsA is a negative regulator of the TA operon through binding to the promoter. Consistent with the MqsR/MqsA pair in E. coli, the binding of the toxin MqsR to MqsA inhibits the DNA binding ability of MqsA in P. putida KT2440. Disruption of the mqsA gene which induces mqsR expression increases persister cell formation 53-fold, while overexpressing mqsA which represses mqsR expression reduces persister cell formation 220-fold, suggesting an important role of MqsR in persistence in P. putida KT2440. Furthermore, both MqsR and MqsA promote biofilm formation. As a DNA binding protein, MqsA can also negatively regulate an ECF sigma factor AlgU and a universal stress protein PP_3288. Thus, we revealed an important regulatory role of MqsR/MqsA in persistence and biofilm formation in P. putida KT2440.

Introduction

Toxin/antitoxin (TA) systems were initially discovered in 1983 as plasmid addiction systems on low-copy-number plasmids due to their ability to stabilize plasmids in the bacterial host population (Ogura and Hiraga, 1983). After that, thousands of TA loci have been identified on plasmids, prophages as well as on the chromosomes of bacteria and archaea (Pandey and Gerdes, 2005; Yamaguchi et al., 2011). TA systems have received increasing attention due to their prevalence, diverse structures, and important physiological functions. In the last decade, TA systems in bacterial chromosomes have been extensively studied. The roles of TA systems in cell physiology, specifically in biofilm formation, persister cell formation (Dorr et al., 2010; Kim and Wood, 2010), the general stress response and phage inhibition are becoming clearer (Pecota and Wood, 1996; Shah et al., 2006; Magnuson, 2007; Tsilibaris et al., 2007; Kim et al., 2009). Till now, TA systems have been classified in to six different groups (Types I–VI) according to the nature of the antitoxin and the mode of interaction between the toxin and antitoxin (Page and Peti, 2016). About 40 TA systems have been reported in Escherichia coli, including 18 type I, 19 type II, 1 type IV, 1 type V, and 1 type VI (Baba et al., 2006; Yamaguchi and Inouye, 2011; Masuda et al., 2012; Wang et al., 2012; Aakre et al., 2013; Guo et al., 2014). More than 88 type II TAs have been identified in human pathogen Mycobacterium tuberculosis (Winther et al., 2016).

The most-studied chromosomal encoded TA systems include MazE/MazF, RelB/RelE, MqsR/MqsA, HigB/HigA, and VapC/VapB which belong to the type II TA systems. Toxins of these TA systems function as endoribonucleases, and regulate protein production by selectively degrading RNAs (Wang and Wood, 2011). Physiological functions of some of these type II toxins are associated with persister cell formation and biofilm formation. It has been reported that persister cell formation was reduced after deleting 10 toxins of the type II TA loci in E. coli (Maisonneuve et al., 2011). An earlier study showed that mqsR was the most highly induced gene in persister cells compared to non-persisters (Shah et al., 2006). The MqsR/MqsA TA system in E. coli has been shown to regulate biofilm formation (Barrios et al., 2006; Kim et al., 2010; Wang et al., 2011; Hong et al., 2012). MqsR/MqsA pair of E. coli is the first TA linked to biofilm formation since mqsR was induced in a transcriptome study that identified genes that were differentially regulated in biofilm cells (Ren et al., 2004). MqsR encodes a ribosome-independent that cleaves mRNA primarily at GCU sites (Yamaguchi et al., 2009; Christensen-Dalsgaard et al., 2010). Overproduction of MqsR led to a global change in the transcriptome profile due to its mRNA interferase activity which has substrate specificity for cellular mRNAs (Kim et al., 2010). Specifically, 14 mRNAs that lack a GCU site (Yamaguchi et al., 2009) were enriched when MqsR was activated including a novel type I toxin RalR (Guo et al., 2014) and the first reported Type V toxin GhoT (Wang et al., 2012). Six out of these 14 mRNAs were differentially regulated in biofilm (Domka et al., 2007). Moreover, antitoxin MqsA represses the expression of stationary sigma factor rpoS which reduces the concentration of the internal messenger c-di-GMP and lead to increased motility and decreased biofilm formation in E. coli (Wang et al., 2011). Thus, both MqsR and MqsA play an important role in E. coli biofilm formation.

Many Pseudomonas are known for their adaptability to thrive in diverse and even hostile habitats (Silby et al., 2011). Pseudomonas putida is an archetypical environmental microbe and is commonly found in various soil and aquatic environments including those were heavily polluted (Ramos et al., 2009). P. putida strain KT2440 genome is fully sequenced, and it serves as a standard reference to compare with other P. putida strains and with pathogenic Pseudomonas aeruginosa strains (Dos Santos et al., 2004). P. putida and P. aeruginosa are closely related and share approximately 85% of the predicted coding regions, but P. putida lacks the key virulence genes such as the type III secretion systems (Nelson et al., 2002). Genomic analysis of P. putida KT2440 also revealed that the high diversity of mobile elements, extracytoplasmatic function (ECF) sigma factors and stress response regulators are indicative of its renowned tolerance of environmental stresses and adaptation to diverse habitats (Dos Santos et al., 2004). TA systems may also contribute to the stress tolerance of P. putida living in heavily polluted environments through regulating biofilm formation and persistence. The first TA system characterized in pseudomonads was the type II TA system GraT/GraA of P. putida KT2440 which belong to the HigB/HigA TA family (Tamman et al., 2014). GraT of the GraT/GraA system exhibited moderate toxin at optimal growth temperatures but caused a severe growth defect at lower temperatures (Tamman et al., 2014). Recent work demonstrated that GraT toxin can also increase resistance to different antibiotics (Tamman et al., 2014), and inhibit ribosome biogenesis at low temperatures by interacting with the DnaK chaperone (Ainelo et al., 2016). Recently study of the HigB/HigA TA system in P. aeruginosa showed that the activation of HigB reduced the production of virulence factors and biofilm formation (Wood and Wood, 2016), while another study showed the activation of HigB in P. aeruginosa increased type III secretion system-mediated cytotoxicity (Li et al., 2016).

Besides E. coli, MqsR/MqsA TA system is distributed in a large number of genera such as Pseudomonas, Yesinia, Burkholderia, and Xylella (Makarova et al., 2009). Recent work showed that in plant pathogen Xylella fastidiosa, overexpressing the MqsR homolog increased persister cell formation and biofilm formation while repressed cell motility (Merfa et al., 2016). PP_4205/PP_4204 in P. putida KT2440 share medium identity (∼50%) with the MqsR/MqsA in E. coli K-12, however, the function of MqsR/MqsA in Pseudomonas remains unknown. Our previous work showed that MqsR variants not only exhibit different toxicity but also exhibit different ability to induce persister formation in E. coli (Hong et al., 2012). Thus, to explore the role of MqsR and MqsA homologs in different bacterial hosts would help gain more insights into the prevalence of TA system and how they participate in bacterial stress response in different hosts.

In this study, we first characterized PP_4205 and PP_4204 and further explored their function in P. putida KT2440. The two genes encoding MqsR homolog and MqsA homolog are co-transcribed and both encode small proteins. MqsR and MqsA interact with each other directly in vivo and MqsA is a negative regulator of the TA operon through binding to the promoter. Consistent with the MqsR/MqsA pair in E. coli, the MqsR inhibits the binding of MqsA to its promoter in P. putida KT2440. Disruption of the mqsA gene which induces mqsR expression increases persister cell formation while overexpressing mqsA which represses mqsR expression reduces persister cell formation, suggesting an important role of MqsR in persistence in P. putida KT2440. However, different from E. coli, both MqsR and MqsA promote biofilm formation in P. putida KT2440. As a DNA binding protein, MqsA in E. coli also regulates the stationary phase sigma factor RpoS and CsgD. In contrast, we found that MqsA negatively regulates the ECF sigma factor AlgU and a universal stress protein PP_3288 in P. putida KT2440, suggesting a different role of MqsR/MqsA in P. putida KT2440.

Materials and Methods

Bacterial Strains, Plasmids, and Growth Conditions

The bacterial strains and plasmids used in this study are listed in Table 1, and the sequences of the primers are listed in Supplementary Table S1. The E. coli and P. putida KT2440 strains were grown in Luria Bertani broth (LB) medium at 30°C except where indicated. Gentamicin (15 μg/ml) was used for pre-culturing isogenic knockout mutants and for maintaining pEX18Ap-based plasmids, carbenicillin (100 μg/ml) was used to maintain plasmid pMQ70, and kanamycin (50 μg/ml) was used to maintain plasmids pET28b, pHGE, and pHGR01.

TABLE 1

TABLE 1. Bacterial strains and plasmids used in this study.

Construction of Deletion Mutants

To create the in-frame deletion of mqsR, mqsA, and mqsRA in P. putida KT2440, 5′ and 3′ homologous arm segments of 0.9–1.5 kb were generated by PCR from genomic DNA, and each amplicon was then ligated with Gm^R-encoding cassette (from pPS856) into pEX18Ap to produce the replacement plasmids as described previously (Hoang et al., 1998). For example, to delete mqsA, 1334 bp upstream and 962 bp downstream fragments of mqsA ORF were generated by PCR (primers used were listed in Supplementary Table S1) using P. putida KT2440 genomic DNA. PCR product of each amplicon was then ligated with Gm^R-encoding cassette into pEX18Ap which contains the sacB-based counter-selection marker to create the replacement plasmid. Next, the replacement plasmid was transformed into KT2440 by electroporation and the resulting deletion mutants (mqsA gene was replaced by Gm via homologous recombination) were selected on LB plate supplemented with gentamicin (Hoang et al., 1998). Deletion of the mqsR gene and the mqsRA operon were performed following similar steps, and all three deletion mutants were confirmed by PCR and DNA sequencing with primers listed in Supplementary Table S1.

Construction of Plasmids

To overexpress mqsR, mqsA, and mqsRA in E. coli and P. putida KT2440 hosts, the full coding region of mqsR, mqsA, and mqsRA were amplified with primer pairs listed in Supplementary Table S1 using P. putida KT2440 genomic DNA as the template. PCR products were purified using a gel extraction kit (Qiagen, Hilden, Germany), digested with EcoRI (or NheI for pHGE-base plasmids) and HindIII, and were purified with a PCR product purification kit (Qiagen). The purified PCR products were ligated into the pMQ70 and pHGE expression plasmids and transferred into E. coli WM3064 and P. putida KT2440, respectively. The correct constructs were verified by DNA sequencing using primer pairs pMQ70-f/r and pHGE-f/r.

For purification of the MqsA protein and MqsR/MqsA complex, the coding region of mqsA and mqsRA were amplified with primer pairs pET28b-mqsA-f/-r to make pET28b-mqsA-His, pET28b-mqsRA-f/pET28b-mqsRA-His-r to make pET28b-mqsRA-CHis (Supplementary Table S1) using P. putida KT2440 genomic DNA as the template, PCR products were purified using a gel extraction kit (Qiagen), and digested with NcoI and HindIII, and were purified with a PCR product purification kit (Qiagen). The purified PCR products were ligated into the pET28b plasmid and transferred into E. coli BL21. The correct constructs were verified by DNA sequencing using primer T7.

Construction of Reporter Strains for Promoter Activity Assays

For the promoter activity assay, DNA fragments containing 409 nt upstream of the translational start site of mqsR, a 328 nt upstream of the translational start site of PP_3288, and a 357 nt upstream of the translational start site of algU were selected as the promoter regions. pHGR01-PmqsRA-lacZ contains the promoter region of mqsRA, pHGR01-PmqsRA-mqsA-lacZ contains the promoter region and the ORF of the mqsA. Both of the fragments were digested with EcoRI and HindIII and cloned into the promoter-less lacZ-fusion vector pHGR01 digested with the same two enzymes. Constructed pHGR01-PmqsRA-lacZ was verified by DNA sequencing using the primer pair pHGR01-f/-r listed in Supplementary Table S1. pHGR01-PPP_3288-lacZ and pHGR01-PalgU-lacZ were constructed following the same procedure. The resulting plasmids were verified by DNA sequencing and transformed into E. coli WM3064 strain. Mid-log phase (OD₆₀₀ ∼ 0.7) cells of the indicated strains carrying the reporter plasmids were collected by centrifugation and washed with phosphate buffered saline. The cell soluble protein and β-galactosidase activity were determined using previously described protocols (Wu et al., 2011).

Protein Expression and Purification

MqsA or the MqsR/MqsA complex with a hexahistidine tag at the C-terminus of MqsA and the MqsR/MqsA complex without any tag were purified using E. coli BL21 with pET28b-mqsA-His, pET28b-mqsRA-CHis, and pET28b-mqsRA. Strains were grown in LB with kanamycin (50 μg/ml) and were induced with 1 mM IPTG at a turbidity of 0.1 at 600 nm for 5 h. Cells were collected and resuspended in lysis buffer [50 mM potassium phosphate buffer (pH 8.0), 300 mM NaCl, and protease inhibitor cocktail (Sigma–Aldrich, USA)]. Samples were sonicated using a Sonic Dismembrator (Ningbo dongzhi, China) at level 2 for 5 min on 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 the Bi Yuntian BCA assay kit (Haimen, Jiangsu, China). Tricine–SDS–PAGE was performed as previously described (Baba et al., 2006). A total of 20 μg of protein from each sample was loaded for SDS–PAGE.

Quantitative Real-Time Reverse-Transcription PCR (qRT-PCR)

Total RNA was isolated as described previously (Ren et al., 2004). The cDNA synthesis was conducted using reverse transcription (Promega, Madison, WI, USA). Total cDNA (50 ng) was used for qRT-PCR using the Step One Real-Time PCR System (Applied Biosystems StepOne Real-Time PCR System, USA). Primers used for qRT-PCR are listed in Supplementary Table S1, and the level of 16S rRNA gene transcript was used to normalize the gene expression data, and fold change of each gene was calculated as described previously (Pfaffl, 2001).

Attached Biofilm Assay

Attached biofilm formation was assayed in 96-well polystyrene plates with 0.1% crystal violet staining as described previously (Pratt and Kolter, 1998). To remove growth effects, biofilm formation was normalized by dividing the maximal bacterial growth (measured by the turbidity at 620 nm) for each strain. The assay was repeated twice with at least three independent cultures.

Electrophoretic Mobility Shift Assay (EMSA)

Electrophoretic mobility shift assays (EMSAs) were performed as previously described (Kim et al., 2010). The DNA probes were amplified by PCR using primers listed in Supplementary Table S1. The PCR products were gel purified and labeled with biotin using the Biotin 3′-end DNA Labeling Kit (Pierce, Rockford, IL, USA). For the binding reactions, biotin-probes (0.05 pmol DNA) were incubated with purified MqsA-His protein for 1 h at room temperature. The DNA binding reactions were performed with the non-specific competitor DNA (poly dI-dC) and NP-40 in buffer containing 10 mM HEPES (pH 7.3), 20 mM KCl, 1 mM MgCl₂, and 5% glycerol in a total volume of 20 μl at 25°C for 1 h. Aliquots of binding reaction mixtures were then electrophoresed on 6% DNA acrylamide retardation gels (Invitrogen, Carlsbad, CA, USA) at 100 V in 0.5 × TBE buffer (10 mM Tris borate at pH 8.3 and 2 mM EDTA) for 90 min. DNA was transferred to a nylon membrane at 390 mA for 40 min, followed by UV crosslinking at 302 nm for 20 min. Chemiluminescence was performed with the LightShift Chemiluminescent EMSA Kit (Thermo Scientific, Rockford, IL, USA) according to the manufacturer’s protocol.

Persistence Assay

Persistence of P. putida KT2440 was measured by time-dependent killing experiments (Dorr et al., 2010). To test the formation of persister cells in P. putida KT2440, overnight bacterial culture was re-inoculated into fresh LB medium (at a turbidity of 0.1 at 600 nm) and grown to an OD₆₀₀ of 1.0. For strains carrying pHGE-based plasmids, overnight cultures were diluted to a turbidity of 0.1 and grown to a turbidity of 0.5, then 1 mM IPTG was used to induce gene expression for 2 h, and the turbidity was adjusted to 1.0. Then the bacterial cultures were exposed to 1.25 μg/ml ciprofloxacin. Cell viability (CFU/ml) was determined at the time points indicated. Cell survival was determined by drop assay.

Phylogenetic Analysis

MqsR homologs were retrieved from ortholog clusters (K13651) in the KEGG ORTHOLOGY database and the data in a previous constructed phylogenetic tree of MqsR (Makarova et al., 2009). By excluding a few nearly identical sequences from different strains of the same bacterial species, the data set of total 133 genes were used for our phylogenetic analysis. The amino acid sequences were aligned with the MAFFT program (L-INS-i) (Katoh and Standley, 2013) and the poorly aligned sites were deleted by trimAl (with gt 0.8) (Capella-Gutierrez et al., 2009). A final data set of 98 amino acid sites was used for the maximum likelihood (ML) phylogenetic analysis by the W-IQ-TREE (Trifinopoulos et al., 2016). The best-fit substitution model was automatically determined and the reliability of internal branches was tested by 1000 ultrafast bootstrap replicates (Minh et al., 2013) in the W-IQ-TREE web interface. The ML tree was rooted to split the two groups of MqsR described previously (Makarova et al., 2009) and further annotated by the iTOL tool (Letunic and Bork, 2016).

Palindrome Search

The motifs 5′-ACCT (N)₃ AGGT-3′ were identified using the Biostrings and BSgenome libraries in the R statistical package (version 2.9.2) and Fuzznuc (EMBOSS) (Rice et al., 2000).

Statistical Analysis

Data are presented as means ± SE of three or more independent cultures. Statistical significance was assessed using two-tailed unpaired Student’s t-test.

Results

Phylogenetic Analyses of MqsR/MqsA Homologs

A total of 143 genes encoding for MqsR homologs were found in various bacterial genomes by searching the KEGG ORTHOLOGY database. As shown in Figure 1, MqsR toxin seems to be divergent in different genera. In consistent with what has been reported earlier (Makarova et al., 2009), MqsR toxin belongs to two phylogenetically distinct groups (Figure 1). Group 1 comprises mostly proteobacterial MqsR toxins, while Group 2 comprises MqsR toxins of Firmicutes, Bacteroidetes, and Chlorobi. MqsR homologs of Proteobacteria basically fall into two separate clades. The Proteobacteria clade 1 mainly comprises MqsRs of α-Proteobacteria and δ-Proteobacteria, and the Proteobacteria clade 2 comprises those of β-Proteobacteria and γ-Proteobacteria. MqsR homologs of P. putida seem more closely related with that of E. coli and Yersinia spp. than other Pseudomonas MqsR toxins.

FIGURE 1

FIGURE 1. Phylogenetic tree of MqsR homologs. The maximum likelihood tree was rooted to split into two groups of MqsR homologs in 133 bacterial genomes. The branch of MqsR of Pseudomonas putida KT2440 is highlighted in red. Branch lengths are proportional to the number of amino acid substitutions (see scale bar in left) and gray circle on branches show ultrafast bootstrap support values that are greater than 50%. Taxa are colored according to the taxonomy.

Identification of MqsR/MqsA as a TA Pair in P. putida KT2440

Putative type II TA systems in the P. putida KT2440 genome were predicted with the web-based tool RASTA-Bacteria (Sevin and Barloy-Hubler, 2007). Two neighboring genes, PP_4205 and PP_4204, encoding two proteins of 98 aa and 133 aa, were identified as a putative TA pair (Figure 2A and Supplementary Figure S1). PP_4205 shares 58 % identity with MqsR in E. coli, and the four amino acid residues (Lys⁵⁶, Gln⁶⁹, Tyr⁸¹, and Lys⁹⁶) which are important for MqsR toxicity of E. coli (Brown et al., 2009) are also conserved in P. putida KT2440 (Figure 2B). PP_4204 shares 45% identity with MqsA in E. coli at the amino acid sequence level, the two residues Asn⁹⁷ and Arg¹⁰¹ that are important for MqsA to bind DNA and for transcription regulation in E. coli (Brown et al., 2009) are also conserved in P. putida KT2440 (Figure 2C).

FIGURE 2

FIGURE 2. MqsR (PP_4205) is toxic and MqsA (PP_4204) neutralizes the toxicity of MqsR. (A) The mqsRA locus in P. putida KT2440. The number indicates the relative position to the start codon of mqsR. (B) Comparison of amino acid sequence of MqsR in Escherichia coli K-12 and in P. putida KT2440. Triangles below the consensus sequences (K59, Q69, V81, and K96) show the important residues for toxicity of E. coli MqsR and squares show the MqsA-interacting residues (Brown et al., 2009). (C) Comparison of amino acid sequences and predicted secondary structures of MqsA in E. coli K-12 and in P. putida KT2440. Triangles below the consensus sequences (N97 and R101) show the important residues for MqsA to bind DNA and for transcription regulation in E. coli (Brown et al., 2009). (D) Growth of the P. putida KT2440 strains harboring the pHGE-based plasmids were induced with 0.5 mM IPTG at OD₆₀₀∼0.1. (E) Cell viability (CFUs/ml) was determined at the time points indicated. Error bars indicate standard error of mean (n = 3) in D,E. (F) E. coli WM3064 hosts harboring pMQ70-based plasmids were streaked on LB plates supplemented with 100 μg/ml ampicillin with or without 10 mM L-arabinose.

To explore the function of this TA pair in P. putida KT2440, we first tested the toxicity of the two-gene cassette. The coding regions of the two genes were cloned into the pHGE plasmid to construct pHGE-PP_4204 (pHGE-mqsA) and pHGE-PP_4205 (pHGE-mqsR) using genomic DNA of P. putida KT2440 as template. When transformed into P. putida KT2440, cells expressing MqsR using pHGE-mqsR exhibited a notable decrease in cell growth and viability as shown by the reduction in turbidity at 600 nm (OD₆₀₀) and colony forming units (CFUs) (Figures 2D,E). In contrast, expression of MqsA using pHGE-mqsA did not affect cell growth in P. putida KT2440 (Figures 2D,E).

To further determine whether the neighboring protein MqsA can neutralize the toxicity of MqsR, we constructed the pHGE-mqsRA plasmid to co-express MqsA and MqsR in P. putida KT2440 cells. As expected, MqsA completely neutralized the toxicity of MqsR (Figures 2D,E). Moreover, co-expression of MqsA with MqsR effectively inhibited cell death caused by MqsR overproduction (Figures 2D,E). To further test the toxicity of these two proteins in a different host, we constructed and transferred pMQ70-based constructs into E. coli WM3064 host. As expected, MqsR was toxic while MqsA neutralized the toxicity of MqsR in E. coli host (Figure 2F). Therefore, PP_4205 and PP_4204 in P. putida KT2440 consist of a type II TA pair.

mqsR and mqsA Are Co-transcribed

To determine whether mqsR and mqsA are co-transcribed, we designed a forward primer in the first gene mqsR (mqsR-RT-f) and a reverse primer in the second gene mqsA (mqsA-RT-r) (Supplementary Figure S1) and conducted a reverse transcription polymerase chain reaction (RT-PCR) assay. As shown in Supplementary Figure S2, a single band of 707 bp was amplified using cDNA synthesized from total RNA as the template, and the PCR product was sequenced to be the region among the two primers, indicating mqsR and mqsA form a single operon (Supplementary Figure S2, lane 3). As controls, the same band was detected using genomic DNA (Supplementary Figure S2, lane 2) as the template but not for total RNA (Supplementary Figure S2, lane 4). This result was also supported by qRT-PCR results as the same levels of transcripts were determined for mqsR and mqsA gene in the normal growing P. putida KT2440 cells (Supplementary Table S2). Collectively, these results demonstrate that mqsR and mqsA are co-transcribed.

MqsR Interacts with MqsA In Vivo

Direct protein-protein interaction between toxin and antitoxin in vivo is a typical feature of type II TA systems (Brown et al., 2011). To determine whether MqsR binds to MqsA, we performed a pull-down assay using pET28b-mqsRA-CHis to co-express a C-terminus hexahistidine tagged (His-tagged) MqsA along with an untagged toxin MqsR (Yao et al., 2015; Guo et al., 2016). Affinity purification using Ni-NTA agarose beads and subsequent Tricine–SDS–PAGE revealed that another small protein was pulled down in addition to His-tagged MqsA (15.51 kDa) (Figure 3A, lane 4), and this small protein was verified to be MqsR (11.13 kDa) toxin by mass spectrometry (Supplementary Table S3). We also constructed pET28b-mqsRA to express untagged MqsR and untagged MqsA, and neither of them could bind to Ni-NTA beads (Figure 3B, lane 4), suggesting MqsR was purified due to its interaction with MqsA-His. Hence, we showed that inhibition of MqsR toxicity by antitoxin MqsA was most likely due to the direct interaction between the two proteins.

FIGURE 3

FIGURE 3. MqsR and MqsA form a complex in vivo. (A) His-tagged MqsA and untagged MqsR were produced from pET28b-mqsRA-CHis in E. coli BL21. MqsA-CHis (15.5 kDa) and MqsR (11.1 kDa) were induced with IPTG (lane 3). MqsR was co-purified with MqsA-CHis (lane 4). Cells that were not induced were used as the negative control (lane 2). (B) Cells harboring pET28b-mqsRA was used as negative control. Neither untagged MqsA nor untagged MqsR bound to the Ni-NTA agarose beads (lane 4). Protein marker (M) was loaded in lane 1 in A,B.

MqsA Represses Transcription of the mqsRA Operon

It has been shown that the type II antitoxin alone or the TA complex binds to its promoter and negatively regulates the transcription of the TA operon (Yamaguchi et al., 2011). Analysis of the promoter region of mqsRA operon revealed that it contains a 17 bp palindrome 5′-TTAACCT GGA AGGTTAA-3′, and the palindrome is overlapped the ribosome-binding site (RBS) for mqsR (Figure 4A). To check whether MqsA binds to the promoter region, we first performed EMSA using purified MqsA protein (Supplementary Figure S3) and a PCR product of 130 bp that includes the promoter region of the TA operon (Figure 4A and Supplementary Figure S1). MqsA specifically bound to the mqsRA promoter region in a concentration-dependent manner (Figure 4B). To determine whether MqsA binds to the palindromic sequence, we performed EMSA using a mutant DNA probe with the palindrome disrupted, and results showed that MqsA did not bind to the mutant DNA probe (Figure 4C). When the MqsR/MqsA complex purified using pET28b-mqsRA-CHis was used for EMSA, the TA complex could not bind and shift the mqsRA promoter (Figure 4D). This result is consistent with the previous work in E. coli which demonstrated that the MqsR/MqsA complex was unable to bind DNA when MqsR and MqsA are present at 1:1 ratio (Brown et al., 2013). These in vitro results suggest that the autoregulation of the mqsRA operon relies on the ratio of the toxin protein to the antitoxin protein present in the cell. Thus, similar to MqsR in E. coli, MqsR in P. putida KT2440 does not function as a transcriptional co-repressor, but blocks the antitoxin MqsA to bind DNA when they are present at approximately 1:1 ratio in the cell.

FIGURE 4

FIGURE 4. MqsA binds the mqsRA operon. (A) The promoter region of the mqsRA operon contains a 17 bp palindrome 5′-TTAACCT GGA AGGTTAA-3′ and the ribosome-binding site (RBS) for mqsR which is marked with boxes. (B) Electrophoretic mobility shift assay (EMSA) results showed that MqsA-CHis bound and shifted the promoter region of mqsRA (PmqsRA). (C) The binding of MqsA to PmqsRA was abolished after random mutating the palindrome to 5′-TTAACCT GGA TCACAGC-3′. (D) MqsR/MqsA complex did not bind to PmqsRA. DNA probe PmqsRA used for EMSAs is the same as A. (E) Schematic diagram of the two constructed reporter plasmids for promoter activity assay. Mid-log phase WM3064 cells carrying two different reporter plasmids (lower left) and overexpressing mqsA or mqsRA in the presence of one of the reporter plasmids (lower right) were tested for β-galactosidase activity. Three independent cultures for each strain were used and the data are shown as means ± standard deviations. Asterisks represent a statistically significant difference in E (P < 0.01 was shown in ^∗ and P < 0.001 was shown in ^∗∗; n = 3).

Moreover, a promoter activity assay using a transcriptional fusion was employed to study the autoregulation of the MqsR/MqsA TA pair in vivo. We cloned the promoter region of mqsRA into the pHGR01 plasmid to create the pHGR01-PmqsRA-lacZ reporter plasmid. To examine the effect of MqsA, we also fused the promoter region of mqsRA and the coding region of mqsA to the same promoterless lacZ reporter gene to create pHGR01-PmqsRA-mqsA-lacZ (Figure 4E, up panel). The promoter activity was decreased from 63 ± 4.0 miller units (MU) in the E. coli WM3064 cells carrying the pHGR01-PmqsRA-lacZ plasmid to 8 ± 0.1 MU in cells harboring the pHGR01-PmqsRA-mqsA-lacZ (Figure 4E, left panel), suggesting that the presence of MqsA repressed the promoter activity of pHGR01-PmqsRA-lacZ. In addition, we also overexpressed mqsA using pMQ70-mqsA in the reporter strain WM3064/pHGR01-PmqsRA-lacZ. The promoter activity reduced from 66.0 ± 5.0 MU to 3.0 ± 0.2 MU after overexpressing MqsA in the E. coli WM3064 cells. However, when TA complex was overproduced by pMQ70-mqsRA, the promoter activity was not significantly reduced in the reporter strain (Figure 4E, right panel). Therefore, MqsA represses transcription of the mqsRA operon, and the presence of MqsR blocks the activity of MqsA to bind DNA when it is co-expressed.

MqsA and MqsR Affect Biofilm Formation and Persistence

To probe the physiological function of MqsR/MqsA in P. putida KT2440, we constructed in-frame deletion mutants of mqsR, mqsA, and mqsRA in P. putida KT2440. The deletion mutants were verified by PCR (Figure 5A) followed by DNA sequencing to make sure no other mutations occurred in the neighboring gene. Deletion of the toxin gene or the antitoxin gene did not significantly affect cell growth (Figure 5B). Since TA systems are related to biofilm formation (Wang and Wood, 2011), attached biofilm formation was tested for the deletion mutants using the 96-well polystyrene plate assay. Attached biofilm formation was reduced 8.0 ± 0.1-fold, 6.3 ± 0.1-fold, and 6.8 ± 0.1-fold after 2 h for the mqsR, mqsA, and mqsRA deletion mutants compared to the wild-type KT2440 strain, respectively (Figure 5C). Moreover, overexpressing mqsA via pHGE-mqsA in the mqsA deletion mutant strain increased biofilm formation and overexpressing mqsRA via pHGE-mqsRA also increased biofilm formation (Figure 5D). We did not perform complementation study for mqsR due to the high toxicity of MqsR when overexpressed in P. putida KT2440.

FIGURE 5

FIGURE 5. MqsA and MqsR affect biofilm formation. (A) Deletion of mqsR, mqsA, and mqsRA operon in P. putida KT2440 are checked by PCR. PCR products are indicated with the expected sizes using genomic DNA from the wild-type (WT) and the deletion mutants. M indicates DNA ladder (500, 750, 1000, 1500, 2000, 3000, and 5000 bp). (B) Growth of different strains in LB medium. (C) Attached biofilm formation of different strains quantified in the 96-well polystyrene plate. (D) Attached biofilm formation of different strains (with 0.5 mM IPTG) using the 96-well polystyrene plate assay in LB medium. Three independent cultures for each strain were used and the data are shown as means ± standard deviations. Asterisks represent a statistically significant difference in C and D. (P < 0.01 was shown in ^∗ and P < 0.001 was shown in ^∗∗; n = 3).

To test the role of MqsR/MqsA in persistence, persister cell assay was performed for the mqsA or mqsR mutant strains using a high dose of ciprofloxacin. As shown in Figure 6A, deletion of mqsA increased the formation of persister cells 52.7 ± 4.2-fold, while deletion of mqsR or mqsRA did not significantly affected persister cell formation. We further complemented MqsA in the ΔmqsA strain via pHGE-mqsA, and results demonstrated that overproduction of MqsA inhibited the formation of persister cells 219.7 ± 10.2-fold compared to the ΔmqsA/pHGE strain (Figure 6B). Since the deletion of mqsA would activate the expression of mqsR, we conducted qRT-PCR to check the expression level of mqsR in the mqsA deletion mutant strain. As shown in Supplementary Table S2, the transcription of the toxin gene mqsR was induced 2.0 ± 0.1-fold in the ΔmqsA strain compared to the wild-type strain. In addition, expression of mqsR was repressed 5.6 ± 0.2-fold in the ΔmqsA/pHGE-mqsA strain compared to the ΔmqsA/pHGE-mqsA strain. Since there was no difference in persister cell formation between the ΔmqsR and ΔmqsRA strain, and deletion of mqsR did not affect mqsA expression as shown in the ΔmqsR strain (Supplementary Table S2), MqsR should be mainly responsible for the increased persistence detected in the mqsA deletion mutant and also responsible for the reduced persistence in the mqsA complementation strain in P. putida KT2440.

FIGURE 6

FIGURE 6. Role of MqsR and MqsA in persister cell formation. (A) Survival for the P. putida KT2440 wild-type, ΔmqsR, ΔmqsA, and ΔmqsRA strains after challenging with 1.25 μg/ml ciprofloxacin (added at OD₆₀₀∼1.0 at T = 0). (B) Survival for the P. putida KT2440 wild-type, ΔmqsA, ΔmqsA/pHGE, and ΔmqsA/pHGE-mqsA strains after challenging with 1.25 μg/ml ciprofloxacin at OD₆₀₀∼1.0 (0.5 mM IPTG was added at T = 0). At indicated time points, survival rates were determined and error bars indicate standard error of mean (n = 3) in A,B.

MqsA Negatively Regulates Other Genes

We have previously demonstrated that MqsA in E. coli can bind to the promoter of rpoS and csgD via the MqsA-like palindrome 5′-ACCT (N)₃ AGGT-3′ where N is any nucleotide. We found that although MqsR and MqsA in E. coli K-12 and in P. putida KT2440 share medium level of identify at the amino acid sequence level, the palindrome where MqsA binds are the same in these two strains. However, different from E. coli K-12, there is no MqsA-like palindrome in the gene encoding the stationary sigma factor RpoS in P. putida KT2440, and moreover, the master regulator CsgD for curli and cellulose production is not found in P. putida KT2440. Therefore, we conducted a whole-genome search to identify intergenic regions with 5′-ACCT (N)₃ AGGT-3′ in P. putida KT2440. As shown in Table 2, the MqsA-like palindrome was found in the upstream region of another three genes in P. putida KT2440, including algU, PP_3288, and nadB.

TABLE 2

TABLE 2. List of all genes with MqsA-like palindrome 5′-ACCT (N)₃ AGGT-3′ in the promoter region in Pseudomonas putida KT2440 identified by the whole-genome search.

AlgU is one of the ECF sigma factors in Pseudomonas, which is required for full resistance of P. aeruginosa to oxidative stress and responsible for alginate biosynthesis (Martinezsalazar et al., 1996; Tart et al., 2005; Bazire et al., 2010). PP_3288 is predicted to be a universal stress family protein but its exact function has not been elucidated. NadB is L-aspartate amino acid oxidase and its activity is required for nicotinamide adenine dinucleotide (NAD) biosynthesis in both aerobic and anaerobic conditions in E. coli (Seifert et al., 1990). We first conducted qRT-PCR to check the expression of these three genes in the mqsA deletion mutant strain as compared to the wild-type P. putida KT2440 strain, and the transcription of the three genes were significantly induced (about twofold to threefold) as compared to the wild-type strain. In addition, overexpressing mqsA via pHGE-mqsA repressed the expression of these genes about threefold (Figure 7A). Next, we performed a promoter activity assay using a transcriptional fusion to study the regulation of the MqsA with algU in vivo. We cloned the promoter region of algU into the pHGR01 plasmid to create the reporter plasmid pHGR01-PalgU-lacZ. The promoter activity was decreased from 216.5 ± 2.3 MU to 75.2 ± 5.2 MU after overexpressing mqsA via pMQ70-mqsA in the E. coli WM3064 that harbored the reporter plasmid (Figure 7B). However, when the TA complex MqsR/MqsA was overproduced via pMQ70-mqsRA, the promoter activity was not significantly changed. Similar results were obtained for the PP_3288 gene. The promoter activity was decreased from 109.6 ± 5.1 MU to 8.8 ± 0.1 MU after overexpressing mqsA in the E. coli WM3064 that harbored the pHGR01-PPP_3288-lacZ plasmid, and when the MqsR/MqsA complex was overproduced by pMQ70-mqsRA, the promoter activity was not significantly reduced (Figure 7B). Collectively, these results suggest that MqsA regulates more than its own locus in P. putida KT2440, and this activity is also blocked by MqsR when it is co-expressed.

FIGURE 7

FIGURE 7. MqsA negatively regulates algU, PP_3288 and nadB. (A) Fold changes of expression of three genes in ΔmqsA vs. KT2440, and ΔmqsA/pHGE-mqsA vs. ΔmqsA/pHGE were quantified by qRT-PCR. (B) Mid-log phase cells of the reporter strains harboring pMQ70, pMQ70-mqsRA, and pMQ70-mqsA (with 10 mM L-arabinose) were tested for β-galactosidase activity, respectively. Three independent cultures for each strain were used and the data are shown as means ± standard deviations. Asterisks represent a statistically significant difference. (P < 0.01 was shown in ^∗ and P < 0.001 was shown in ^∗∗; n = 3).

Discussion

In this study, we present evidence to support that MqsR (PP_4205) and MqsA (PP_4204) in P. putida KT2440 form a type II TA system. These results are: (i) mqsR and mqsA form an operon (PP_4205-PP_4204), and they are co-transcribed; (ii) both of them encode small proteins; (iii) MqsR functions as a toxin that inhibits growth, and MqsA blocks the toxicity of MqsR by direct protein–protein interaction; (iv) the antitoxin MqsA regulates transcription of the mqsRA operon by binding to its promoter; and (v) the binding of MqsR to MqsA inhibits its DNA binding ability. These features fit well with the MqsR/MqsA pair in E. coli K-12. Although these homologous proteins share medium level of amino acid identity, the previously identified key residues for MqsR toxicity and for the DNA binding of MqsA in E. coli K-12 are also conserved in P. putida KT2440 (Figures 2B,C). MqsR homologs from diverse bacterial lineages belong to two groups (Figure 1). In our analysis, Group 1 basically falls into two separate clades comprising different classes of the Proteobacteria phylum. In particular, MqsR homologs of Pseudomonas spp. and different P. putida strains further nest into two different clades (Figure 1). Thus, further study is needed to explore the function of MqsR in different phylogenetic groups and in different clades.

The MqsR/MqsA system has been shown to regulate biofilm formation, curli/cellulose production, stress resistance, motility, and persistence in E. coli (Barrios et al., 2006; Kim et al., 2010; Wang et al., 2011; Hong et al., 2012). Deletion of mqsR reduced biofilm formation suggests that MqsR increases biofilm formation in P. putida KT2440. Furthermore, MqsR increases persister cell formation in P. putida KT2440. The MqsR toxin identified in E. coli belongs to the RelE superfamily and functions as ribosome-independent mRNase (Yamaguchi et al., 2009; Christensen-Dalsgaard et al., 2010). MqsR cleaves mRNAs in vivo and in vitro, the endoribonuclease activity of the MqsR was inhibited by MqsA (Yamaguchi et al., 2009). Primer extension analysis of the in vivo and in vitro cleaved RNA indicated that MqsR cleaves RNA before or after the G residue in GCU sequences, and cleavage of some GCA/G/C sites could also be detected in vivo (Yamaguchi et al., 2009; Christensen-Dalsgaard et al., 2010). Our attempt to purify the MqsR was unsuccessful due to its high toxicity, similar to previous studies (Robson et al., 2009). Co-purification of MqsR and MqsA resulted in strong interaction between MqsR and MqsA, which made it difficult to separate them without affecting the activity of the toxin. It has been reported that overproduction of MqsR caused elongated cells in E. coli and in X. fastidiosa (Yamaguchi et al., 2009; Merfa et al., 2016). However, overexpressing MqsR in P. putida KT2440 did not show any alteration in cell morphology (Supplementary Figure S4). Moreover, the type I toxin gene ralR and the type V toxin gene ghoT which are enriched when MqsR is produced in E. coli (Wang et al., 2013; Guo et al., 2014) are both absent in P. putida KT2440. These results suggest that MqsR should have different cellular targets and is involved in specific biological processes in different bacteria.

Except for binding to its cognate toxin MqsR and to its own promoter, antitoxin MqsA also regulates other important loci. Specially, we have demonstrated that the decrease of biofilm in E. coli is due to MqsA binding to the promoter of rpoS and repressing expression of rpoS, decreased levels of RpoS reduces the concentration of c-di-GMP (Wang et al., 2011). The binding of MqsA to the promoter of cgsD is also reported which confirms the role of MqsA in regulating other genetic loci (Soo and Wood, 2013). The binding of MqsA to the promoter region of rpoS and csgD are through the MqsA-like palindrome 5′-ACCT (N)₃ AGGT-3′ present in the promoters of these two genes. Moreover, the sequences of the three nucleotides spacer do not seem to have any effect on MqsA binding affinity to DNA (Brown et al., 2011). In this study, we found that, although there is very low sequence similarity at the promoter region of the mqsRA operon in P. putida KT2440 and in E. coli K12, the palindromes where MqsA binds are the same in these two strains. In E. coli, the removal of mqsR increased biofilm formation (Kim et al., 2010), while the removal of the mqsRA operon reduced biofilm formation (Kasari et al., 2010). However, in P. putida KT2440, deletion of mqsR and deletion of mqsA both reduced biofilm formation, and the role of MqsR and MqsA in biofilm formation was further shown by the complementation study in which overexpressing the TA complex or MqsA alone increased biofilm formation. Thus, it seems that both MqsR and MqsA are able to regulate biofilm formation of P. putida KT2440. Moreover, different from E. coli K-12 in which the removal of mqsA gene was lethal, the deletion of the antitoxin gene mqsA was not lethal in P. putida KT2440. The non-lethal effect of mqsA deletion in KT2440 could be partially explained by a low expression of the mqsRA operon and/or a relatively low level of MqsR production in the absence of its cognate antitoxin in the planktonic growth condition (Supplementary Table S2). A relatively lower expression of the mqsRA operon was also found in the planktonic growth condition in E. coli (Domka et al., 2007). Considering that there are two identical palindromes for the binding of MqsA in the promoter region of the mqsRA operon in E. coli but only one in P. putida KT2440, it is possible that the repression of MqsA on the mqsR transcription in P. putida KT2440 is less efficient than in E. coli. Deletion of the antitoxin gene graA of the GraT/GraA type II TA pair was also successful in P. putida KT2440; however, the level of graT mRNA was much higher in the graA deletion mutant than that of the wild-type strain (Tamman et al., 2014).

As a DNA binding protein, MqsA in E. coli binds to the promoter of rpoS gene and the csgD genes. However, different from E. coli K-12, there is no MqsA-like palindrome in the promoter region of the stationary sigma factor RpoS in P. putida KT2440. In addition, the master regulator CsgD for curli and cellulose production is not found in P. putida KT2440. Interestingly, we found another gene encoding an alternative sigma factor AlgU containing the MqsA-like palindrome in the upstream region in P. putida KT2440. ECF sigma factors are a group of alternative sigma factors whose target gene products often function outside the cytoplasm and many of them are related to virulence and pathogenesis in Pseudomonas (Lonetto et al., 1994). AlgU is a positive regulator of alginate production and negative regulator of flagella-associated motility in P. aeruginosa (Tart et al., 2005). Both alginate and flagella are important for the biofilm formation of P. aeruginosa. Alginate promotes cells to adhere to biotic and abiotic surfaces and can also function as the scaffolding material of P. aeruginosa biofilms (Halverson, 2009). Motility is necessary for the bacteria to reach surfaces suitable for biofilm formation and spreading of the biofilms (Otoole and Kolter, 1998; Klausen et al., 2003). RpoS and AlgU are both global regulators that control expression of many proteins in conjunction with various regulatory factors (Hengge, 2008; Schulz et al., 2015). Thus, both the toxin MqsR and antitoxin MqsA may play important regulatory roles in P. putida including biofilm formation and virulence. It remains interesting to investigate the regulation and impact of TA systems on the modulation of a single cell or in a population or during host–pathogen interaction.

Author Contributions

XW conceived and designed the experiments, CS, YG, KT, ZW, BL, and ZZ performed the experiments and data analysis. XW and CS wrote the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by the National Science Foundation of China (Grant No. 31290233, 31625001, and 31500025), the Chinese Academy of Sciences (XDA11030402) and by National Science Foundation of Guangdong Province (2015A030310405). XW is a recipient of the 1000-Youth Elite Program (the Recruitment Program of Global Experts in China).

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2017.00840/full#supplementary-material

References

Aakre, C. D., Phung, T. N., Huang, D., and Laub, M. T. (2013). A bacterial toxin inhibits DNA replication elongation through a direct interaction with the beta sliding clamp. Mol. Cell. 52, 617–628. doi: 10.1016/j.molcel.2013.10.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Ainelo, A., Tamman, H., Leppik, M., Remme, J., and Horak, R. (2016). The toxin GraT inhibits ribosome biogenesis. Mol. Microbiol. 100, 719–734. doi: 10.1111/mmi.13344

PubMed Abstract | CrossRef Full Text | Google Scholar

Baba, T., Ara, T., Hasegawa, M., Takai, Y., Okumura, Y., Baba, M., et al. (2006). Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol. Syst. Biol. 2, 1–11. doi: 10.1038/msb4100050

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrios, A. F. G., Zuo, R. J., Hashimoto, Y., Yang, L., Bentley, W. E., and Wood, T. K. (2006). Autoinducer 2 controls biofilm formation in Escherichia coli through a novel motility quorum-sensing regulator (MqsR, B3022). J. Bacteriol. 188, 305–316. doi: 10.1128/Jb.188.1.305-316.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Bazire, A., Shioya, K., Soum-Soutera, E., Bouffartigues, E., Ryder, C., Guentas-Dombrowsky, L., et al. (2010). The sigma factor AlgU plays a key role in formation of robust biofilms by nonmucoid Pseudomonas aeruginosa. J. Bacteriol. 192, 3001–3010. doi: 10.1128/Jb.01633-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, B. L., Grigoriu, S., Kim, Y., Arruda, J. M., Davenport, A., Wood, T. K., et al. (2009). Three dimensional structure of the MqsR:MqsA complex: a novel TA pair comprised of a toxin homologous to RelE and an antitoxin with unique properties. PLoS Pathog. 5:e1000706. doi: 10.1371/journal.ppat.1000706

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, B. L., Lord, D. M., Grigoriu, S., Peti, W., and Page, R. (2013). The Escherichia coli toxin MqsR destabilizes the transcriptional repression complex formed between the antitoxin MqsA and the mqsRA operon promoter. J. Biol. Chem. 288, 1286–1294. doi: 10.1074/jbc.M112.421008

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, B. L., Wood, T. K., Peti, W., and Page, R. (2011). Structure of the Escherichia coli antitoxin MqsA (YgiT/b3021) bound to its gene promoter reveals extensive domain rearrangements and the specificity of transcriptional regulation. J. Biol. Chem. 286, 2285–2296. doi: 10.1074/jbc.M110.172643

PubMed Abstract | CrossRef Full Text | Google Scholar

Capella-Gutierrez, S., Silla-Martinez, J. M., and Gabaldon, T. (2009). trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25, 1972–1973. doi: 10.1093/bioinformatics/btp348

PubMed Abstract | CrossRef Full Text | Google Scholar

Christensen-Dalsgaard, M., Jorgensen, M. G., and Gerdes, K. (2010). Three new RelE-homologous mRNA interferases of Escherichia coli differentially induced by environmental stresses. Mol. Microbiol. 75, 333–348. doi: 10.1111/j.1365-2958.2009.06969.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Domka, J., Lee, J. T., Bansal, T., and Wood, T. K. (2007). Temporal gene-expression in Escherichia coli K-12 biofilms. Environ. Microbiol. 9, 332–346. doi: 10.1111/j.1462-2920.2006.01143.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dorr, T., Vulic, M., and Lewis, K. (2010). Ciprofloxacin causes persister formation by inducing the TisB toxin in Escherichia coli. PLoS Biol. 8:e1000317. doi: 10.1371/journal.pbio.1000317

PubMed Abstract | CrossRef Full Text | Google Scholar

Dos Santos, V. A., Heim, S., Moore, E. R. B., Stratz, M., and Timmis, K. N. (2004). Insights into the genomic basis of niche specificity of Pseudomonas putida KT2440. Environ. Microbiol. 6, 1264–1286. doi: 10.1111/j.1462-2920.2004.00734.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, H., Jin, M., Ju, L., Mao, Y., and Gao, H. (2014). Evidence for function overlapping of CymA and the cytochrome bc(1) complex in the Shewanella oneidensis nitrate and nitrite respiration. Environ. Microbiol. 16, 3181–3195. doi: 10.1111/1462-2920.12457

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Y., Quiroga, C., Chen, Q., McAnulty, M. J., Benedik, M. J., Wood, T. K., et al. (2014). RalR (a DNase) and RalA (a small RNA) form a type I toxin-antitoxin system in Escherichia coli. Nucleic Acids Res. 42, 6448–6462. doi: 10.1093/Nar/Gku279

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Y., Yao, J., Sun, C., Wen, Z., and Wang, X. (2016). Characterization of the deep-sea Streptomyces sp. SCSIO 02999 derived VapC/VapB toxin-antitoxin system in Escherichia coli. Toxins 8:195. doi: 10.3390/toxins8070195

PubMed Abstract | CrossRef Full Text | Google Scholar

Halverson, L. J. (2009). Alginates: Biology and Applications. Berlin: Springer.

Google Scholar

Hengge, R. (2008). The two-component network and the general stress sigma factor RpoS (σS) in Escherichia coli. Adv. Exp. Med. Biol. 631, 40–53. doi: 10.1007/978-0-387-78885-2_4

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoang, T. T., Karkhoff-Schweizer, R. R., Kutchma, A. J., and Schweizer, H. P. (1998). A broad-host-range Flp-FRT recombination system for site-specific excision of chromosomally-located DNA sequences: application for isolation of unmarked Pseudomonas aeruginosa mutants. Gene 212, 77–86. doi: 10.1016/S0378-1119(98)00130-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Hong, S. H., Wang, X., O’Connor, H. F., Benedik, M. J., and Wood, T. K. (2012). Bacterial persistence increases as environmental fitness decreases. Microb. Biotechnol. 5, 509–522. doi: 10.1111/j.1751-7915.2011.00327.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasari, V., Kurg, K., Margus, T., Tenson, T., and Kaldalu, N. (2010). The Escherichia coli mqsR and ygiT genes encode a new toxin-antitoxin pair. J. Bacteriol. 192, 2908–2919. doi: 10.1128/Jb.01266-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, Y., Wang, X., Ma, Q., Zhang, X., and Wood, T. K. (2009). Toxin-antitoxin systems in Escherichia coli influence biofilm formation through YjgK (TabA) and fimbriae. J. Bacteriol. 191, 1258–1267. doi: 10.1128/Jb.01465-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, Y., Wang, X., Zhang, X., Grigoriu, S., Page, R., Peti, W., et al. (2010). Escherichia coli toxin/antitoxin pair MqsR/MqsA regulate toxin CspD. Environ. Microbiol. 12, 1105–1121. doi: 10.1111/j.1462-2920.2009.02147.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, Y., and Wood, T. K. (2010). Toxins Hha and CspD and small RNA regulator Hfq are involved in persister cell formation through MqsR in Escherichia coli. Biochem. Biophys. Res. Commun. 391, 209–213. doi: 10.1016/j.bbrc.2009.11.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Klausen, M., Heydorn, A., Ragas, P., Lambertsen, L., Aaesjorgensen, A., Molin, S., et al. (2003). Biofilm formation by Pseudomonas aeruginosa wild type, flagella and type IV pili mutants. Mol. Microbiol. 48, 1511–1524. doi: 10.1046/j.1365-2958.2003.03525.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Letunic, I., and Bork, P. (2016). Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 44, W242–W245. doi: 10.1093/nar/gkw290

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Long, Y., Liu, Y., Liu, Y., Chen, R., Shi, J., et al. (2016). HigB of Pseudomonas aeruginosa enhances killing of phagocytes by up-regulating the type III secretion system in ciprofloxacin induced persister cells. Front. Cell Infect. Microbiol. 6:125. doi: 10.3389/fcimb.2016.00125

PubMed Abstract | CrossRef Full Text | Google Scholar

Lonetto, M. A., Brown, K. L., Rudd, K. E., and Buttner, M. J. (1994). Analysis of the Streptomyces coelicolor sigE gene reveals the existence of a subfamily of eubacterial RNA polymerase sigma-factors involved in the regulation of extracytoplasmic functions. Proc. Natl. Acad. Sci. U.S.A. 91, 7573–7577. doi: 10.1073/pnas.91.16.7573

PubMed Abstract | CrossRef Full Text | Google Scholar

Magnuson, R. D. (2007). Hypothetical functions of toxin-antitoxin systems. J. Bacteriol. 189, 6089–6092. doi: 10.1128/JB.00958-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Maisonneuve, E., Shakespeare, L. J., Jorgensen, M. G., and Gerdes, K. (2011). Bacterial persistence by RNA endonucleases. Proc. Natl. Acad. Sci. U.S.A. 108, 13206–13211. doi: 10.1073/pnas.1100186108

PubMed Abstract | CrossRef Full Text | Google Scholar

Makarova, K. S., Wolf, Y. I., and Koonin, E. V. (2009). Comprehensive comparative-genomic analysis of type 2 toxin-antitoxin systems and related mobile stress response systems in prokaryotes. Biol. Direct 4:19. doi: 10.1186/1745-6150-4-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinezsalazar, J. M., Moreno, S., Najera, R., Boucher, J., Espin, G., Soberonchavez, G., et al. (1996). Characterization of the genes coding for the putative sigma factor AlgU and its regulators MucA, MucB, MucC, and MucD in Azotobacter vinelandii and evaluation of their roles in alginate biosynthesis. J. Bacteriol. 178, 1800–1808. doi: 10.1128/jb.178.7.1800-1808.1996

PubMed Abstract | CrossRef Full Text | Google Scholar

Masuda, H., Tan, Q., Awano, N., Wu, K., and Inouye, M. (2012). YeeU enhances the bundling of cytoskeletal polymers of MreB and FtsZ, antagonizing the CbtA (YeeV) toxicity in Escherichia coli. Mol. Microbiol. 84, 979–989. doi: 10.1111/j.1365-2958.2012.08068.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Merfa, M. V., Niza, B., Takita, M. A., and De Souza, A. A. (2016). The MqsRA toxin-antitoxin system from Xylella fastidiosa plays a key role in bacterial fitness, pathogenicity, and persister cell formation. Front. Microbiol. 7:904. doi: 10.3389/fmicb.2016.00904

PubMed Abstract | CrossRef Full Text | Google Scholar

Minh, B. Q., Nguyen, M. A. T., and von Haeseler, A. (2013). Ultrafast approximation for phylogenetic bootstrap. Mol. Biol. Evol. 30, 1188–1195. doi: 10.1093/molbev/mst024

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelson, K. E., Weinel, C., Paulsen, I. T., Dodson, R. J., Hilbert, H., Dos Santos, V. A., et al. (2002). Complete genome sequence and comparative analysis of the metabolically versatile Pseudomonas putida KT2440. Environ. Microbiol. 4, 799–808. doi: 10.1046/j.1462-2920.2003.00480.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogura, T., and Hiraga, S. (1983). Mini-F plasmid genes that couple host cell division to plasmid proliferation. Proc. Natl. Acad. Sci. U.S.A. 80, 4784–4788. doi: 10.1073/pnas.80.15.4784

PubMed Abstract | CrossRef Full Text | Google Scholar

Otoole, G. A., and Kolter, R. (1998). Flagellar and twitching motility are necessary for Pseudomonas aeruginosa biofilm development. Mol. Microbiol. 30, 295–304. doi: 10.1046/j.1365-2958.1998.01062.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Page, R., and Peti, W. (2016). Toxin-antitoxin systems in bacterial growth arrest and persistence. Nat. Chem. Biol. 12, 208–214. doi: 10.1038/Nchembio.2044

PubMed Abstract | CrossRef Full Text | Google Scholar

Pandey, D. P., and Gerdes, K. (2005). Toxin-antitoxin loci are highly abundant in free-living but lost from host-associated prokaryotes. Nucleic Acids Res. 33, 966–976. doi: 10.1093/Nar/Gki201

PubMed Abstract | CrossRef Full Text | Google Scholar

Pecota, D. C., and Wood, T. K. (1996). Exclusion of T4 phage by the hok/sok killer locus from plasmid R1. J. Bacteriol. 178, 2044–2050. doi: 10.1128/jb.178.7.2044-2050.1996

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfaffl, M. W. (2001). A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29:e45. doi: 10.1093/nar/29.9.e45

CrossRef Full Text | Google Scholar

Pratt, L. A., and Kolter, R. (1998). Genetic analysis of Escherichia coli biofilm formation: roles of flagella, motility, chemotaxis and type I pili. Mol. Microbiol. 30, 285–293. doi: 10.1046/j.1365-2958.1998.01061.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramos, J. L., Krell, T., Daniels, C., Segura, A., and Duque, E. (2009). Responses of Pseudomonas to small toxic molecules by a mosaic of domains. Curr. Opin. Microbiol. 12, 215–220. doi: 10.1016/j.mib.2009.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, D., Bedzyk, L. A., Thomas, S. M., Ye, R. W., and Wood, T. K. (2004). Gene expression in Escherichia coli biofilms. Appl. Microbiol. Biotechnol. 64, 515–524. doi: 10.1007/s00253-003-1517-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Rice, P., Longden, I., and Bleasby, A. (2000). EMBOSS: the European molecular biology open software suite. Trends Genet. 16, 276–277. doi: 10.1016/S0168-9525(00)02024-2

CrossRef Full Text | Google Scholar

Robson, J., McKenzie, J. L., Cursons, R., Cook, G. M., and Arcus, V. L. (2009). The vapBC operon from Mycobacterium smegmatis is an autoregulated toxin-antitoxin module that controls growth via inhibition of translation. J. Mol. Biol. 390, 353–367. doi: 10.1016/j.jmb.2009.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Schulz, S., Eckweiler, D., Bielecka, A., Nicolai, T., Franke, R., Dotsch, A., et al. (2015). Elucidation of sigma factor-associated networks in Pseudomonas aeruginosa reveals a modular architecture with limited and function-specific crosstalk. PLoS Pathog. 11:e1004744. doi: 10.1371/journal.ppat.1004744

PubMed Abstract | CrossRef Full Text | Google Scholar

Seifert, J., Kunz, N., Flachmann, R., Läufer, A., Jany, K. D., and Gassen, H. G. (1990). Expression of the E. coli nadB gene and characterization of the gene product L-aspartate oxidase. Biol. Chem. Hoppe Seyler 371, 239–248. doi: 10.1515/bchm3.1990.371.1.239

PubMed Abstract | CrossRef Full Text | Google Scholar

Sevin, E. W., and Barloy-Hubler, F. (2007). RASTA-bacteria: a web-based tool for identifying toxin-antitoxin loci in prokaryotes. Genome Biol. 8:R155. doi: 10.1186/Gb-2007-8-8-R155

PubMed Abstract | CrossRef Full Text | Google Scholar

Shah, D., Zhang, Z., Khodursky, A. B., Kaldalu, N., Kurg, K., and Lewis, K. (2006). Persisters: a distinct physiological state of E. coli. BMC Microbiol. 6:53. doi: 10.1186/1471-2180-6-53

PubMed Abstract | CrossRef Full Text | Google Scholar

Shanks, R. M. Q., Caiazza, N. C., Hinsa, S. M., Toutain, C. M., and Otoole, G. A. (2006). Saccharomyces cerevisiae-based molecular tool kit for manipulation of genes from gram-negative bacteria. Appl. Environ. Microbiol. 72, 5027–5036. doi: 10.1128/AEM.00682-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Silby, M. W., Winstanley, C., Godfrey, S. A. C., Levy, S. B., and Jackson, R. W. (2011). Pseudomonas genomes: diverse and adaptable. FEMS Microbiol. Rev. 35, 652–680. doi: 10.1111/j.1574-6976.2011.00269.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Soo, V. W. C., and Wood, T. K. (2013). Antitoxin MqsA represses curli formation through the master biofilm regulator CsgD. Sci. Rep. 3:3186. doi: 10.1038/Srep03186

PubMed Abstract | CrossRef Full Text | Google Scholar

Stover, C. K., Pham, X. Q., Erwin, A. L., Mizoguchi, S. D., Warrener, P., Hickey, M. J., et al. (2000). Complete genome sequence of Pseudomonas aeruginosa PAO1, an opportunistic pathogen. Nature 406, 959–964. doi: 10.1038/35023079

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamman, H., Ainelo, A., Ainsaar, K., and Horak, R. (2014). A moderate toxin, GraT, modulates growth rate and stress tolerance of Pseudomonas putida. J. Bacteriol. 196, 157–169. doi: 10.1128/Jb.00851-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Tart, A. H., Wolfgang, M. C., and Wozniak, D. J. (2005). The alternative sigma factor AlgT represses Pseudomonas aeruginosa flagellum biosynthesis by inhibiting expression of fleQ. J. Bacteriol. 187, 7955–7962. doi: 10.1128/Jb.187.23.7955-7962.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Trifinopoulos, J., Nguyen, L. T., von Haeseler, A., and Minh, B. Q. (2016). W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. 44, W232–W235. doi: 10.1093/nar/gkw256

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsilibaris, V., Maenhaut-Michel, G., Mine, N., and Van Melderen, L. (2007). What is the benefit to Escherichia coli of having multiple toxin-antitoxin systems in its genome? J. Bacteriol. 189, 6101–6108. doi: 10.1128/Jb.00527-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Kim, Y., Hong, S. H., Ma, Q., Brown, B. L., Pu, M., et al. (2011). Antitoxin MqsA helps mediate the bacterial general stress response. Nat. Chem. Biol. 7, 359–366. doi: 10.1038/nchembio.560

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Lord, D. M., Cheng, H. Y., Osbourne, D. O., Hong, S. H., Sanchez-Torres, V., et al. (2012). A new type V toxin-antitoxin system where mRNA for toxin GhoT is cleaved by antitoxin GhoS. Nat. Chem. Biol. 8, 855–861. doi: 10.1038/Nchembio.1062

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Lord, D. M., Hong, S. H., Peti, W., Benedik, M. J., Page, R., et al. (2013). Type II toxin/antitoxin MqsR/MqsA controls type V toxin/antitoxin GhoT/GhoS. Environ. Microbiol. 15, 1734–1744. doi: 10.1111/1462-2920.12063

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., and Wood, T. K. (2011). Toxin-antitoxin systems influence biofilm and persister cell formation and the general stress response. Appl. Environ. Microbiol. 77, 5577–5583. doi: 10.1128/Aem.05068-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Winther, K., Tree, J. J., Tollervey, D., and Gerdes, K. (2016). VapCs of Mycobacterium tuberculosis cleave RNAs essential for translation. Nucleic Acids Res. 44, 9860–9871. doi: 10.1093/nar/gkw781

PubMed Abstract | CrossRef Full Text | Google Scholar

Wood, T. L., and Wood, T. K. (2016). The HigB/HigA toxin/antitoxin system of Pseudomonas aeruginosa influences the virulence factors pyochelin, pyocyanin, and biofilm formation. Microbiologyopen 5, 499–511. doi: 10.1002/mbo3.346

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, L., Wang, J., Tang, P., Chen, H., and Gao, H. (2011). Genetic and molecular characterization of flagellar assembly in Shewanella oneidensis. PLoS ONE 6:e21479. doi: 10.1371/journal.pone.0021479

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamaguchi, Y., and Inouye, M. (2011). Regulation of growth and death in Escherichia coli by toxin-antitoxin systems. Nat. Rev. Microbiol. 9, 779–790. doi: 10.1038/Nrmicro2651

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamaguchi, Y., Park, J., and Inouye, M. (2009). MqsR, a crucial regulator for quorum sensing and biofilm formation, is a GCU-specific mRNA interferase in Escherichia coli. J. Biol. Chem. 284, 28746–28753. doi: 10.1074/jbc.M109.032904

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamaguchi, Y., Park, J. H., and Inouye, M. (2011). Toxin-antitoxin systems in bacteria and archaea. Annu. Rev. Genet. 45, 61–79. doi: 10.1146/annurev-genet-110410-132412

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, J., Guo, Y., Zeng, Z., Liu, X., Shi, F., and Wang, X. (2015). Identification and characterization of a HEPN-MNT family type II toxin-antitoxin in Shewanella oneidensis. Microb. Biotechnol. 8, 961–973. doi: 10.1111/1751-7915.12294

PubMed Abstract | CrossRef Full Text | Google Scholar