# The EAL-domain protein FcsR regulates flagella, chemotaxis and type III secretion system in Pseudomonas aeruginosa by a phosphodiesterase independent mechanism

Jessica Rossello, Analía Lima, Magdalena Gil, Jorge Rodríguez Duarte, Agustín Correa, Paulo C Carvalho, Arlinet Kierbel, Rosario Durán

<figure><img src="image_1.png" /></figure>

## To cite this version:

Jessica Rossello, Analía Lima, Magdalena Gil, Jorge Rodríguez Duarte, Agustín Correa, et al. The EAL-domain protein FcsR regulates flagella, chemotaxis and type III secretion system in Pseudomonas aeruginosa by a phosphodiesterase independent mechanism. Scientific Reports, 2017, 7 (1), pp.10281. 10.1038/s41598-017-09926-3. pasteur-05052343

HAL Id: pasteur-05052343

https://pasteur.hal.science/pasteur-05052343v1

Submitted on 30 Apr 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
SCIENTIFIC REPORTS

<figure><img src="image_3.png" /></figure>

OPEN

# The EAL-domain protein FcsR regulates flagella, chemotaxis and type III secretion system in *Pseudomonas aeruginosa* by a phosphodiesterase independent mechanism

Jessica Rossello<sup>1</sup>, Analía Lima<sup>1</sup>, Magdalena Gil<sup>1,5</sup>, Jorge Rodríguez Duarte<sup>1</sup>, Agustín Correa<sup>2</sup>, Paulo C. Carvalho<sup>3</sup>, Arlinet Kierbel<sup>4</sup> & Rosario Durán<sup>1</sup>

Received: 21 June 2016

Accepted: 2 August 2017

Published online: 31 August 2017

The second messenger c-di-GMP regulates the switch between motile and sessile bacterial lifestyles. A general feature of c-di-GMP metabolism is the presence of a surprisingly large number of genes coding for diguanylate cyclases and phosphodiesterases, the enzymes responsible for its synthesis and degradation respectively. However, the physiological relevance of this apparent redundancy is not clear, emphasizing the need for investigating the functions of each of these enzymes. Here we focused on the phosphodiesterase PA2133 from *Pseudomonas aeruginosa*, an important opportunistic pathogen. We phenotypically characterized *P. aeruginosa* strain K overexpressing PA2133 or its inactive mutant. We showed that biofilm formation and motility are severely impaired by overexpression of PA2133. Our quantitative proteomic approach applied to the membrane and exoprotein fractions revealed that proteins involved in three processes were mostly affected: flagellar motility, type III secretion system and chemotaxis. While inhibition of biofilm formation can be ascribed to the phosphodiesterase activity of PA2133, down-regulation of flagellar, chemotaxis, and type III secretion system proteins is independent of this enzymatic activity. Based on these unexpected effects of PA2133, we propose to rename this gene product FcsR, for Flagellar, chemotaxis and type III secretion system Regulator.

*Pseudomonas aeruginosa* is a ubiquitous Gram negative bacterium causative of a major public health problem. This opportunistic pathogen is able to infect immunocompromised individuals making it one of the most common causes of hospital-acquired infections<fnref n="1" />. *P. aeruginosa* is also responsible for chronic infections, and is the primary cause of mortality in cystic fibrosis patients². A crucial aspect for host colonization and persistence is the ability of *P. aeruginosa* to switch between motile and surface-attached lifestyles³, ⁴.

The cyclic nucleotide bis-(3'-5')-cyclic dimeric GMP (c-di-GMP) is a global second messenger that has a central role in the regulation of bacterial adhesiveness, controlling both cell-cell and cell-surface interactions. Cellular levels of c-di-GMP are determined by its synthesis and hydrolysis catalyzed by diguanylate cyclases and phosphodiesterases, respectively. The current model links low levels of this second messenger with motile planktonic bacteria, while high c-di-GMP intracellular concentrations are associated with increased biofilm formation

<footnote n="1">¹Analytical Biochemistry and Proteomics Unit, Institut Pasteur de Montevideo/Instituto de Investigaciones Biológicas Clemente Estable, Montevideo, Uruguay. ²Unidad de Proteínas Recombinantes, Institut Pasteur de Montevideo, Montevideo, Uruguay. ³Laboratory for Proteomics and Protein Engineering, Carlos Chagas Institute, Fiocruz-Paraná, Curitiba, Brazil. ⁴Instituto de Investigaciones Biotecnológicas Dr. Rodolfo A. Ugalde (IIB-INTECH), Universidad Nacional de San Martín, Consejo Nacional de Investigaciones Científicas y Técnicas (UNSAM-CONICET), San Martín, Buenos Aires, Argentina. ⁵Present address: Unit of Dynamics of Host-Pathogen Interactions, Institut Pasteur, Paris, France. Correspondence and requests for materials should be addressed to R.D. (email: duran@pasteur.edu.uy)</footnote>
and sessile lifestyle⁵⁻⁷. To exert its effects, c-di-GMP acts at different levels controlling gene expression through transcription factors and riboswitches, as well as protein activities or protein subcellular localization⁸⁻¹¹.

Genome sequencing revealed a surprisingly large number of genes encoding proteins with diguanylate cyclase (GGDEF) and phosphodiesterase (EAL and HD-GYP) characteristic domains in different bacteria¹². However, the physiological relevance of this apparent redundancy in the enzymes controlling c-di-GMP levels is still not clear. The analysis of *P. aeruginosa* PAO1 genome revealed the presence of 17 GGDEF domain containing proteins, 5 EAL domain phosphodiesterases, and 16 hybrid proteins carrying both domains¹³. In addition, many of these enzymes present accessory domains, including sensory input and two-component receiver domains, suggesting that c-di-GMP synthesis is regulated by different environmental signals¹⁴,¹⁵. Most of these enzymes are part of the core genome of *P. aeruginosa* and were found in several clinical and environmental isolates using a DNA microarray¹³. A growing body of evidence suggests that specific diguanylate cyclases and phosphodiesterases are linked to different processes, ranging from the control of flagella and pili function, to the regulation of expression of adhesins, exopolysaccharide operons and virulence genes⁶, ⁷, ¹⁶⁻²⁰. Analysis of the phenotypes related to the overexpression or deletion of each GGDEF and EAL containing gene in *P. aeruginosa* suggested that not all the diguanylate cyclases or phosphodiesterases have the same role on biofilm formation and cytotoxicity¹³. These observations points to a c-di-GMP mediated signalling network far more complex than previously thought; thus deserving further investigation of the molecular function(s) of each GGDEF and EAL containing protein.

In this work, we focused on the EAL domain containing phosphodiesterase PA2133. This 285 amino-acid predicted membrane protein (PSORTb V3.0)²¹ is an active enzyme containing a modified EAL domain (ETL)¹³. According to previous reports, overexpression of this enzyme impairs biofilm production and cytotoxicity, two processes relevant for disease development¹³. Here, we phenotypically characterized *P. aeruginosa* strain K (PAK) overexpressing PA2133 and showed that biofilm formation and motility are severely impaired. Further investigation of the membrane and exoprotein fractions by quantitative proteomics revealed that proteins involved in three processes are mainly affected, namely flagellar motility, type III secretion system (TTSS) and chemotaxis. While inhibition of biofilm formation can be ascribed to the phosphodiesterase activity of PA2133, other proteomic and phenotypic changes are also observed using an enzymatically inactive mutant. In particular, down-regulation of flagellar, TTSS, and chemotaxis proteins does not rely in phosphodiesterase activity. Thus, we suggest renaming PA2133 FcsR, for Flagellar, chemotaxis and type III secretion Regulator.

## Results

**Expression of FcsR increases c-di-GMP phosphodiesterase activity in whole cell extracts.** This work characterizes a *P. aeruginosa* PAK strain overexpressing the EAL phosphodiesterase FcsR or its inactive single point mutant FcsRE<sub>60</sub>A. It is well known that the substitution of the conserved glutamic acid residue of the EAL signature by alanine abolishes the activity of phosphodiesterases<sup>22</sup>. Both PAK/pJN-FcsR and PAK/pJN-FcsRE<sub>60</sub>A presented the same colony morphology and growth rate as wild-type PAK (Supplementary Figure 1).

To evaluate c-di-GMP hydrolyzing activity, cell extracts of PAK, PAK/pJN-FcsR or PAK/pJN-FcsRE<sub>60</sub>A were incubated with synthetic c-di-GMP at different time points and subjected to HPLC analysis<sup>13</sup>. Peaks with the same retention time as the synthetic standard were collected and analyzed by MALDI-TOF/TOF mass spectrometry to corroborate that corresponded to c-di-GMP (theoretical M<sup>+</sup> = 691.06 Da). The identity of the peaks was further certified by comparing their product ion spectra with the c-di-GMP standard and with previously published data<sup>7</sup> (Supplementary Figure 2a). Chromatographic peak area was used for c-di-GMP quantitation. As expected, PAK/pNJ-FcsR showed the highest phosphodiesterase activity. c-di-GMP levels drop approximately five-fold after 3 h of incubation in protein extracts from this strain, while under the same conditions the levels of this cyclic nucleotide remained practically unchanged in protein extracts from PAK and PAK/pJN-FcsRE<sub>60</sub>A (Supplementary Figure 2b).

**Expression of FcsR generates a biofilm-defective phenotype.** Previous work showed that overexpression of FcsR (PA2133) in *P. aeruginosa* PA14 strain significantly impairs biofilm formation¹³. To confirm this result in PAK background, we evaluated the ability of PAK/pJN-FcsR and PAK/pJN-FcsRE<sub>60</sub>A to form floating pellicle biofilms at the air/liquid interface, as well as solid-surface associated biofilms. Pellicle formation on glass tubes is shown in Fig. 1a. Wild-type PAK produced a well-developed pellicle; conversely, PAK/pJN-FcsR's ability to form these multicellular structures was significantly reduced. As expected, the strain overexpressing the inactive phosphodiesterase could form biofurils to a similar extent as wild-type PAK (Fig. 1a).

In addition, biofilm formation on solid surface was evaluated. As shown in Fig. 1b, induction of FcsR expression with 0.02% of L-arabinose led to a substantial decrease in biofilm formation, with a relative biofilm mass of 41% compared to the control. Moreover, higher L-arabinose concentrations (0.05% and 0.2%) yielded a statistically significant reduction of biofilm mass (32% and 23% respectively) (ANOVA and Tukey's *post hoc* test, p < 0.05). In line with previous results, expression of the inactive form of FcsR restored biofilm formation to wild-type levels (Fig. 1b).

Overall, our results showed that the expression of the phosphodiesterase FcsR led to defective biofilm formation phenotype that was dependent on its c-di-GMP hydrolyzing activity. To further investigate the function of FcsR at the molecular level, the membrane proteome and exoproteome of PAK, PAK/pJN-FcsR and PAK/pJN-FcsRE<sub>60</sub>A were compared through quantitative approaches.

**Quantitative membrane proteomics reveals that flagella and chemotaxis are altered in PAK/pJN-FcsR and PAK/pJN-FcsRE<sub>60</sub>A.** Differentially abundant proteins in membrane enriched fractions of PAK and PAK/pJN-FcsR were shortlisted using a quantitative shotgun approach. In these experiments, the proteins present in biological replicates of each strain were identified using mass spectrometry. Next, the qualitative and quantitative comparisons of these datasets using bioinformatics tools allowed the identification of proteins
<figure><img src="image_4.png" /><figcaption>**Figure 1.** PAK/pJN-FcsR is impaired for biofilm formation. (a) Pellicle formation assays. The figure shows the pellicles associated to the glass tubes stained with crystal violet after removal of the liquid culture. This assay was performed for PAK, PAK/pJN-FcsR and PAK/pJN-FcsRE<sub>60</sub>A. The concentrations of L-arabinose used as inducer are indicated. (b) Biofilm formation assays. Biofilm formation on solid surface was evaluated using microtiter dish binding assay for PAK, PAK/pJN-FcsR and PAK/pJN-FcsRE<sub>60</sub>A. L-arabinose concentrations used are indicated. Crystal violet staining retained on the biofilms was measured spectrophotometrically at 570 nm for quantification. *Indicates statistically significant difference determined by ANOVA and Tukey's *post hoc* test (*p* < 0.05).</figcaption></figure>

exclusive to one condition and those whose levels are significantly different between the two conditions. LC-MS/MS analysis resulted in the identification of 651 and 780 proteins in PAK and PAK/pJN-FcsR membrane fractions respectively, with at least two peptides per protein when analysing all replicates per condition (Supplementary Table S1). We verified that 77.5% of the most abundant proteins found in PAK membrane fraction were indeed outer membrane or cytosolic membrane proteins according to *Pseudomonas* genome database<sup>23</sup> (Supplementary Figure 3).

Venn diagram module from PatternLab for Proteomics software was used to pinpoint proteins uniquely identified in each strain<sup>24</sup>. Twenty-five proteins were exclusively detected in wild type PAK (in at least 3 of 4 replicates but absent in all PAK/pJN-FcsR replicates); likewise, 9 were unique to PAK/pJN-FcsR (Table 1). Remarkably, the most abundant proteins among those exclusive to PAK were chemotaxis proteins. In fact, 11 out of the 25 PAK unique proteins were methyl-accepting chemotaxis proteins (MCPs). In addition, P- and M-ring proteins of the flagellar basal body were also exclusively detected in wild-type PAK. Among proteins unique to PAK/pJN-FcsR, we identified PelB, PelC and FimL. PelB and PelC are involved in extracellular polysaccharide biosynthesis, while FimL is required for proper twitching motility, biofilm development and TTSS function<sup>25-27</sup>. Next, we compared the proteins present in both strains to identify those that were differentially abundant according to their normalized spectral abundance factors (NSAF). Using the PatternLab's TFold module statistics and considering five or more replicates in all classes, 23 proteins were shortlisted as differentially abundant (q < 0.05)<sup>28</sup> (Fig. 2 and Table 2). In line with the above results, the most underrepresented proteins in PAK/pJN-FcsR were A-type flagellin (fold change 10.5) and two methyl-accepting chemotaxis proteins (PctA and PctB) with fold changes of 4.1 and
<table><caption>**Table 1.** Proteins uniquely identified in the membrane fraction of PAK or PAK/pJN-FcsR.</caption><thead><tr><th>Description</th><th># replicates</th><th>Spectrum count</th><th>Gene in PAK</th><th>Gene name</th><th>Gene in PAO1</th></tr></thead><tbody><tr><td colspan="6"><strong>Proteins exclusively detected in membrane fraction of PAK</strong></td></tr><tr><td>Chemotaxis protein</td><td>4</td><td>109</td><td>PAK_0510</td><td></td><td>PA4633</td></tr><tr><td>Chemotaxis protein</td><td>4</td><td>107</td><td>PAK_02646</td><td></td><td>PA2654</td></tr><tr><td>Chemotaxis transducer</td><td>4</td><td>84</td><td>PAK_02331</td><td></td><td>PA2867</td></tr><tr><td>Chemotaxis protein</td><td>4</td><td>84</td><td>PAK_03742</td><td></td><td>PA1608</td></tr><tr><td>Chemotaxis protein</td><td>4</td><td>73</td><td>PAK_02492</td><td></td><td>PA2788</td></tr><tr><td>Chemotaxis protein</td><td>4</td><td>62</td><td>Y880_0130565</td><td>PctC</td><td>PA4307</td></tr><tr><td>Chemotaxis protein</td><td>4</td><td>58</td><td>PAK_02648</td><td></td><td>PA2652</td></tr><tr><td>Chemotaxis protein</td><td>4</td><td>48</td><td>PAK_02745</td><td>CtpH</td><td>PA2651</td></tr><tr><td>Aerotaxis receptor Aer</td><td>4</td><td>45</td><td>PAK_03791</td><td>Aer</td><td>PA1561</td></tr><tr><td>Flagellar P-ring protein</td><td>4</td><td>41</td><td>PAK_04298</td><td>FlgI</td><td>PA1084</td></tr><tr><td>Membrane protein</td><td>4</td><td>29</td><td>PAK_01662</td><td></td><td>PA3526</td></tr><tr><td>Uncharacterized protein</td><td>3</td><td>26</td><td>PAK_04850</td><td></td><td>PA4390</td></tr><tr><td>Chemotaxis transducer</td><td>3</td><td>26</td><td>PAK_00397</td><td>CttP</td><td>PA0180</td></tr><tr><td>Uncharacterized protein</td><td>3</td><td>26</td><td>PAK_02120</td><td></td><td>Not found</td></tr><tr><td>Flagellar M-ring protein</td><td>3</td><td>24</td><td>PAK_04270</td><td>FliF</td><td>PA1101</td></tr><tr><td>Heme d1 biosynthesis protein NirF</td><td>4</td><td>22</td><td>PAK_00731</td><td>NirF</td><td>PA0516</td></tr><tr><td>Lipoprotein</td><td>3</td><td>16</td><td>PAK_04784</td><td></td><td>PA4326</td></tr><tr><td>Cell division ATP-binding protein FtsE</td><td>3</td><td>15</td><td>PAK_00589</td><td>FtsE</td><td>PA0374</td></tr><tr><td>Acyltransferase</td><td>3</td><td>14</td><td>PAK_00676</td><td></td><td>PA0461</td></tr><tr><td>Chemotaxis protein</td><td>3</td><td>14</td><td>PAK_03895</td><td></td><td>PA1465</td></tr><tr><td>Chemotaxis protein</td><td>3</td><td>13</td><td>PAK_05418</td><td></td><td>PA4915</td></tr><tr><td>Membrane protein</td><td>3</td><td>12</td><td>PAK_04001</td><td></td><td>PA1365</td></tr><tr><td>Histidine kinase</td><td>3</td><td>12</td><td>PAK_04858</td><td></td><td>PA4398</td></tr><tr><td>Uncharacterized protein</td><td>3</td><td>12</td><td>PAK_03420</td><td></td><td>PA1913</td></tr><tr><td>Amino acid ABC transporter substrate-binding protein</td><td>3</td><td>10</td><td>PAK_01924</td><td>FecA</td><td>PA3268</td></tr><tr><td colspan="6"><strong>Proteins exclusively identified in membrane fraction of PAK/pJN-FcsR</strong></td></tr><tr><td>30S ribosomal protein S12</td><td>4</td><td>26</td><td>PAK_00903</td><td>RpsL</td><td>PA4268</td></tr><tr><td>Uncharacterized protein</td><td>3</td><td>32</td><td>PAK_02131</td><td>PelB</td><td>PA3063</td></tr><tr><td>Glycosyltransferase family 1</td><td>3</td><td>14</td><td>PAK_02044</td><td></td><td>Not found</td></tr><tr><td>Ferrous iron transporter B</td><td>3</td><td>11</td><td>PAK_03517</td><td>FimL</td><td>PA1822</td></tr><tr><td>Transcriptional regulator</td><td>3</td><td>13</td><td>PAK_01156</td><td></td><td>PA4021</td></tr><tr><td>Heat-shock protein</td><td>3</td><td>27</td><td>PAK_02067</td><td>LbpA</td><td>PA3126</td></tr><tr><td>Aprp-1-p processing protein</td><td>3</td><td>15</td><td>PAK_03601</td><td></td><td>PA1746</td></tr><tr><td>Lipoprotein</td><td>3</td><td>34</td><td>PAK_02132</td><td>PelC</td><td>PA3062</td></tr><tr><td>30S ribosomal protein S17</td><td>3</td><td>13</td><td>PAK_00917</td><td>RpsQ</td><td>PA4254</td></tr></tbody></table>

4.7 respectively. As expected, the phosphodiesterase FcsR was among proteins overrepresented in PAK/pJN-FcsR (fold change 61.6 using 0.2% arabinose). In addition, the His-kinase chemotaxis protein ChpA was also more abundant in PAK/pJN-FcsR membrane enriched fractions (fold change 11.2). These findings were unexpected as it is well known that low c-di-GMP levels are associated with increased flagellar motility and flagellar protein synthesis⁷, ²⁹⁻³¹. Interestingly, the analysis of the membrane enriched fraction of PAK/pJN-FcsRE₆₀A indicated that flagellar and chemotaxis proteins were also altered when the inactive phosphodiesterase was expressed. We identified 22 proteins in at least 3 PAK replicates that were not detected either in PAK/pJN-FcsR or in PAK/pJN-FcsRE₆₀A proteomes, including the previously identified chemotaxis proteins as well as the flagellar P- and M-ring proteins (Supplementary Table S2). Moreover, both FcsR and ChpA were also enriched in membrane fractions of PAK/pJN-FcsRE₆₀A, and fold changes were very similar to those observed for PAK/pJN-FcsR (Supplementary Table S3).

Altogether, our results indicate that FcsR expression alters the membrane proteome and leads to underrepresentation of chemotaxis and flagellar proteins by a phosphodiesterase independent mechanism.

**Exoproteome analysis reveals that flagella and TTSS are altered in PAK/pJN-FcsR and PAK/pJN-FcsRE₆₀A.** The exoprotein fraction was analyzed with different workflows. Difference gel electrophoresis (DIGE) was used to compare four biological replicates from PAK and PAK/pJN-FcsR. A total of 823 spots were detected in the master gel, and 27 of them were differential between strains (95% significance and considering a fold change ≥ 1.25). 15 and 12 spots from PAK and PAK/pJN-FcsR, respectively, presented significantly increased volumes. Figure 3 shows a representative gel where the 10 spots identified by MS analysis are pointed
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Quantitative analysis of proteins identified in membrane fractions of PAK and PAK/pJN-FcsR using shotgun approach. The figure shows the Volcano plot generated using the PatternLab for Proteomics TFold module. Each dot in the plot represents a protein identified at least in 5 replicates of all conditions, plotted according its p-value (log2 (p)) and fold change (log2 (fold change)). Red dots represent proteins that do not satisfy neither fold change nor statistic criteria, and thus are considered unchanged between strains. Green dots satisfy fold change criterion but not statistical one. Orange dots correspond to low abundant proteins satisfying both fold change and q value criteria, but due to the low number of spectra they deserved further validation. Finally, blue dots correspond to proteins satisfying all statistical filters and represent the differentially expressed proteins between strains. The identity of each blue dot is shown in the figure (see also Table 2).</figcaption></figure>

(Supplementary Table S4). Remarkably, 7 out of 8 spots underrepresented in PAK/pJN-FcsR gels corresponded to proteoforms of flagellar structure proteins (Type A flagellin and A-type flagellar hook-associated protein). These spots presented large differences between strains, with fold changes that ranged from 4.5 to 15.2 (Supplementary Table S4). The analysis also identified two proteins overrepresented in PAK/pJN-FcsR, namely CdrA and OprG.

These results were validated by our shotgun proteomic analysis. Supplementary Table S5 shows the list of proteins identified in each replicate of PAK and PAK/pJN-FcsR. The analysis of proteins unique to each strain (detected in at least 2 of 3 replicates of one class and absent in the other class) revealed 7 proteins exclusive to PAK/pJN-FcsR and 80 proteins unique to PAK (Supplementary Table S6). The proteins unique to PAK were mainly related to flagella and TTSS. A-type flagellar hook-associated protein and flagellar hook protein FlgK, were detected with high number of spectral counts in PAK (total spectra 1152 and 341 respectively) but could not be detected in the exoproteome of PAK/pJN-FcsR. In line with this finding five additional flagellar structure proteins (FlgE, FlgL, FlgG, FliK and FlaG), plus the anti-sigma factor FlgM, were exclusively detected in PAK exoproteome. Our data also indicated that TTSS is altered in PAK/pJN-FcsR. In particular, two proteins of the translocation pore (PopB, PopD), the needle-tip protein PcrV, the TTSS regulator PopN, and the exotoxins ExoS and ExoT were detected in PAK but not in PAK/pJN-FcsR (Supplementary Table S6). Indeed, the translocator protein PcrV and the effector ExoS presented marked inter-strain differences, with high number of spectra in PAK replicates but absent in PAK/pJN-FcsR. Analysis of differentially abundant proteins, performed by spectral counting, further supported our previous shotgun and DIGE results (Supplementary Table S7). A-type flagellin was overrepresented in PAK (71.5 fold change) and CdrA was overrepresented in PAK/pJN-FcsR (79.8 fold change).

Next, exoproteome analysis of PAK/pJN-FcsRE<sub>60</sub>A clearly indicated that the proteomics changes described above were mostly independent of phosphodiesterase activity. As shown in Supplementary Tables S8, 7 flagella proteins (FlgK, FlgM, FlgE, FlgL, FlgG, FliK, and FlaG) as well as the group of TTSS related proteins (PopB, PopD, PopN, ExoS and ExoT), confidently identified in PAK, were not detected either in PAK/pJN-FcsR or in PAK/
<table><caption>**Table 2.** Proteins identified in membrane fractions of PAK and PAK/pJN-FcsR at significantly different levels.</caption><thead><tr><th>Description</th><th>Fold change</th><th>p value</th><th>Gene in PAK</th><th>Gene name</th><th>Gene in PAO1</th></tr></thead><tbody><tr><td colspan="6"><strong>Proteins overrepresented in PAK</strong></td></tr><tr><td>A-type flagellin</td><td>10.5</td><td>0.0003</td><td>PAK_4280</td><td>FliC</td><td>PA1092</td></tr><tr><td>Chemotaxis protein</td><td>4.7</td><td>0.0009</td><td>PAK_4767</td><td>PctB</td><td>PA4310</td></tr><tr><td>Fe(III) dicitrate transporter</td><td>4.1</td><td>0.0033</td><td>PAK_01276</td><td>FecA</td><td>PA3901</td></tr><tr><td>Chemotaxis protein</td><td>4.1</td><td>0.0017</td><td>Y880_0130550</td><td>PctA</td><td>PA4309</td></tr><tr><td>Cbb3-type cytochrome c oxidase subunit</td><td>2.7</td><td>0.0095</td><td>PAK_3801</td><td>Cytochrome oxidase</td><td>PA1552</td></tr><tr><td>Membrane protein</td><td>1.9</td><td>0.0141</td><td>PAK_1111</td><td>OprG</td><td>PA4067</td></tr><tr><td>Channel protein TolC</td><td>1.5</td><td>0.0168</td><td>PAK_5479</td><td>TolC</td><td>PA4974</td></tr><tr><td>Porin</td><td>1.5</td><td>0.0123</td><td>PAK_2526</td><td>OprQ</td><td>PA2760</td></tr><tr><td colspan="6"><strong>Proteins overrepresented in PAK/pJN-FcsR</strong></td></tr><tr><td>Cyclic-guanilate-specific phosphodiesterase</td><td>61.6</td><td>0.0007</td><td>PAK_03199</td><td>—</td><td>PA2133</td></tr><tr><td>ChpA protein</td><td>11.2</td><td>0.0135</td><td>PAK_00629</td><td>ChpA</td><td>PA0413</td></tr><tr><td>Chaperone protein DnaK</td><td>5.4</td><td>0.0056</td><td>PAK_5258</td><td>DnaK</td><td>PA4761</td></tr><tr><td>Chaperone protein HtpG</td><td>5.3</td><td>0.0143</td><td>PAK_03754</td><td>HtpG</td><td>PA1596</td></tr><tr><td>4-hydroxyphenylpyruvate dioxygenase</td><td>4.5</td><td>0.0051</td><td>PAK_04521</td><td>Hpd</td><td>PA0865</td></tr><tr><td>Aconitatehydratase</td><td>3.8</td><td>0.0114</td><td>PAK_03790</td><td>AcnA</td><td>PA1562</td></tr><tr><td>30S ribosomal protein S13</td><td>2.8</td><td>0.0004</td><td>PAK_00930</td><td>RpsM</td><td>PA4241</td></tr><tr><td>Polyhydroxyalkanoate synthesis protein PhaF</td><td>2.7</td><td>0.0172</td><td>PAK_05567</td><td>PhaF</td><td>PA5060</td></tr><tr><td>Large-conductance mechano sensitive channel</td><td>2.2</td><td>0.0029</td><td>PAK_05082</td><td>MscL</td><td>PA4614</td></tr><tr><td>50S ribosomal protein L21</td><td>2.1</td><td>0.0081</td><td>PAK_05034</td><td>RplU</td><td>PA4568</td></tr><tr><td>High-affinity branched-chain amino acid transport ATP-binding protein BraF</td><td>1.8</td><td>0.0076</td><td>PAK_04311</td><td>BraF</td><td>PA1071</td></tr><tr><td>Cytoplasmic membrane protein</td><td>1.7</td><td>0.0024</td><td>PAK_03572</td><td>CmpX</td><td>PA1775</td></tr><tr><td>Protein tolQ</td><td>1.6</td><td>0.0214</td><td>PAK_04411</td><td>TolQ</td><td>PA0969</td></tr><tr><td>30S ribosomal protein S2</td><td>1.6</td><td>0.0120</td><td>PAK_01531</td><td>RpsB</td><td>PA3656</td></tr><tr><td>Phosphomannomutase/phosphoglucomutase</td><td>1.5</td><td>0.0127</td><td>PAK_05836</td><td>AlgC</td><td>PA5322</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Comparative analysis of exoproteomes of PAK and PAK/pJN-FcsR by 2D DIGE. (a) Multiple images for one DIGE gel showing spot profiles for PAK (Cy3 labeled, green), PAK/pJN-FcsR (Cy5 labeled, red) and internal standard (Cy2 labeled, blue). (b) Overlay image of PAK (green) and PAK/pJN-FcsR (red) channels. Differential spots between strains (considering all four DIGE gels, p-values ≤ 0.05 and fold change ≥ 25%) are indicated and were further identified by MS analysis. Protein identities, fold changes and *p* values are listed in Supplementary Table S4. (c) 3D view of selected spots with significant normalized abundance fold changes: spot 1 CdrA (fold change 4.1) spot 4 flagellar hook-associated protein (fold change 11.2) and spot 7 flagellin (fold change 15.2).</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Flagellar motility inhibition is independent of FcsR phosphodiesterase activity. (a) Swimming motility assays. Assays were performed in triplicates for PAK, PAK/pJN-FcsR and PAK/pJN-FcsRE<sub>60</sub>A strains. The swimming areas were measured and the relative swimming area (% of wild type PAK) was plotted. The concentrations of L-arabinose used as inducer are indicated. *Indicates statistically significant difference determined by ANOVA and Tukey's *post hoc* test, *p* < 0.05. (b) Effect of FcsR and FcsRE<sub>60</sub>A on flagella assembly. Flagella isolated from PAK, PAK/pJN-FcsR (0.2% arabinose), PAK/pJN-FcsRE<sub>60</sub>A (0.2% arabinose) and PAK-FliC<sup>-</sup> strains were analyzed by SDS-PAGE followed by MALDI-TOF/TOF mass spectrometry. The major band of approximately 45 kDa detected in PAK (but absent in all other strains) was unambiguously identified as type A flagellin by mass spectrometry. (c) MALDI mass spectrum of tryptic peptides isolated from the 45 kDa gel band in (b). The *m/z* values observed can be assigned to FliC sequences (*m/z* 1053.4: sequence 126–135; *m/z* 1413.7: sequence 303–316; *m/z* 1615.7: sequence 2–16; *m/z* 2555.2: sequence 222–247; *m/z* 2613.3: sequence 67–91; *m/z* 2874.6: sequence 368–394).</figcaption></figure>

pJN-FcsRE<sub>60</sub>A. Quantitative analysis also revealed that CdrA is largely enriched in the exoproteome of PAK/pJN-FcsRE<sub>60</sub>A if compared with wild-type PAK (Supplementary Table S3).

Thus, the study of different protein fractions using various proteomic approaches showed, in a consistent manner, that overexpression of phosphodiesterase FcsR or its inactive mutant led to down-regulation of proteins involved in flagellar motility, chemotaxis and TTSS; and up-regulation of CdrA, FimL and ChpA. Some of these unexpected proteomic results were further validated at the functional level.

**Flagellar motility inhibition is independent of FcsR phosphodiesterase activity.** Flagellar motility was assessed with a standard swimming assay using increasing L-arabinose concentrations. As shown in Fig. 4a, flagellar motility was severely impaired in PAK/pJN-FcsR and the effect was dependent on L-arabinose concentration. Normalized areas of swimming zones for PAK/pJN-FcsR were 45%, 27% and 15% of wild type with 0.02%, 0.05% and 0.2% of L-arabinose respectively (ANOVA and Tukey's *post hoc* test, *p* < 0.05). PAK/pJN-FcsR swimming motility in medium without L-arabinose was undistinguishable from wild type (data not shown). Interestingly, PAK/pJN-FcsRE<sub>60</sub>A was also impaired in flagellar motility and this phenotype was more pronounced with increasing inducer concentration (Fig. 4a). Areas of swimming zones were 40%, 27% and 11% of wild type strain using 0.02%; 0.05% and 0.2% of L-arabinose respectively. No significant differences were found between PAK/pJN-FcsR and PAK/pJN-FcsRE<sub>60</sub>A flagellar motility for the same L-arabinose concentration, pointing to a phosphodiesterase activity independent swimming defect.

To further investigate the flagellar status of the different strains, flagella isolation was performed. A major gel band of around 45 kDa, identified as type A flagellin (Mascot ion score: 133, sequence coverage: 31%) was detected in PAK but not in PAK/pJN-FcsR, PAK/pJN-FcsRE<sub>60</sub>A or a flagellum negative controls strain (PAK-FliC<sup>-</sup>) (Fig. 4b and c). Moreover, the presence of type IV pili and flagella was clearly seen by transmission electron microscopy in PAK. However, only the pili were present in PAK/pJN-FcsR (Fig. 5 and Supplementary Figure 4). Finally, immunoblotting detected FliC in whole cell extracts of PAK but not in PAK/pJN-FcsR or PAK/pJN-FcsRE<sub>60</sub>A. This result points to the inhibition of biosynthesis of this flagellar structural protein, a process that takes place after assembly of the hook-basal body structure<sup>32</sup> (Supplementary Figure 5).

Altogether, our proteomic, motility, biochemical and microscopy studies revealed that expression of FcsR led to impairment of flagella motility, and that this effect is independent of FcsR-phosphodiesterase activity.

**TTSS regulation is independent of FcsR phosphodiesterase activity.** Proteomic analysis revealed that FcsR expression led to down-regulation of TTSS proteins. Since those experiments were not conducted under optimum conditions for induction of TTSS, the results were confirmed in Ca<sup>2+</sup> depleted medium<sup>33</sup>. The presence of PcrV and ExoS in exoproteomes obtained under these experimental conditions was tested by Western blot. In agreement with proteomic results, both proteins were detected in PAK but not in PAK/pJN-FcsR or PAK/pJN-FcsRE<sub>60</sub>A, indicating that the overexpression of FcsR impairs TTSS function and that this effect was independent of phosphodiesterase activity (Fig. 6).
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Transmission electron microscopy of PAK and PAK/pJN-FcsR. Electron microscope images showing the presence of flagellum and type IV pili in PAK strain. Only pili were observed in PAK/pJN-FcsR. Flagellum is indicated by a black arrow. Type IV pili are indicated by dashed arrows. Scale bar: 500 nm.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** TTSS regulation is independent of FcsR phosphodiesterase activity. Secreted extracts of PAK, PAK/pJN-FcsR and PAK/pJN-FcsRE<sub>60</sub>A grown in Ca<sup>2+</sup> depleted media were separated on SDS-PAGE and the presence of specific TTSS proteins was detected by Western blot. ExoS and PcrV were sequentially detected in the same blotting membrane. Original images are displayed in Supplementary Figure 7. (a) Cropped view of immunoblotting using anti-ExoS antibody. ExoS is only detected in PAK strain (b) Cropped view of immunoblotting using anti-PcrV antibody. PcrV is only detected in PAK strain.</figcaption></figure>

**Flagellar motility phenotype is partially reverted in a ΔChpA strain.** Strikingly, two proteins are highly overrepresented in PAK/pJN-FcsR and PAK/pJN-FcsRE<sub>60</sub>A membrane proteome besides FcsR: FimL and ChpA. Interestingly, FimL presents high degree of homology with the N-terminal domain of ChpA<sup>27</sup>. To evaluate if ChpA could mediate some of the phosphodiesterase-independent phenotypes, we transformed PAO1 and PAO1ΔChpA with the pJN-FcsR plasmid. Expression of FcsR in PAO1 background reproduced the motility phenotype observed in PAK, with an average relative swimming area of 22.7% using 0.2% L-arabinose (Supplementary Figure 6). Interestingly, PAO1ΔChpA/pJN-FcsR has an intermediate motility phenotype (swimming area 48.9% of control) (Supplementary Figure 6). Importantly, PAO1ΔChpA has no flagellar motility defect. While we are still far from understanding the mechanism behind this observation, our results suggest that this His-kinase may play a role in the FcsR-phosphodiesterase independent flagellar motility inhibition.
Discussion

This report examines the role of the EAL-phosphodiesterase PA2133 by overexpressing the active or inactive enzyme in *P. aeruginosa* PAK strain. Phenotypic and differential proteomic analyses revealed that this protein regulates four main processes; namely: biofilm formation, flagellar motility, chemotaxis, and TTSS. Our results clearly showed that while biofilm defective phenotype is mediated by FcsR enzymatic activity, the regulation of flagellar motility, chemotaxis, and TTSS does not rely on c-di-GMP hydrolysis.

We showed that PAK/pJN-FcsR is defective in biofilm formation, in agreement with previous reports for strains with low c-di-GMP levels⁶, ¹³, ³⁴⁻³⁶. In fact, PAK/pJN-FcsRE₆₀A can form biofilm to the same extent as wild type PAK.

On the other hand, two independent proteomic approaches plus phenotypic characterization showed, con-
clusively, that flagella protein biosynthesis and flagella mediated motility are severely affected in PAK/pJN-FcsR.
Overall, 10 flagella structural proteins (that form part of hook, basal body and filament) were underrepresented
in the PAK/pJN-FcsR strain when compared to PAK. Differences in abundance of flagella structural proteins
were systematically detected using two orthogonal experimental approaches and two different cellular fractions.
Particularly A-type flagellin was identified as a differential protein in membrane and exoprotein fractions, while
A-type flagellin and type A associated flagellar hook corresponded to 70% of the differential spots in DIGE anal-
ysis of exoproteome fraction. Next, motility assays validated proteomics data at the functional level and provided
conclusive evidence of swimming motility inhibition. Finally, both electron microscopy and flagella isolation
assays indicated that the overexpression of this EAL phosphodiesterase generates mainly an aflagellated strain.
The lower abundance of flagella structural proteins in PAK/pJN-FcsR was an unexpected result, as it was previ-
ously documented that decreased levels of c-di-GMP correlates with increased flagellin synthesis and flagellar
mediated motility⁷, ²⁹⁻³¹. Interestingly, we provide strong evidence that the regulation of flagellar motility is not
a consequence of low levels of c-di-GMP, as the expression of an inactive mutant produced the same proteomic
pattern and swimming defect.

Our results also showed major changes in chemotaxis proteins in PAK/pJN-FcsR strain. *P. aeruginosa* genome codifies for 26 methyl accepting chemotaxis receptors and only half of them have been characterized so far³⁷. These receptors sense changes in the concentration of attractants or repellents in the environment to generate a behavioral response that directs flagella or pili mediated motility towards or against the chemical gradient. Analysis of membrane proteins indicated that 13 MCP receptors were not detected, or significantly underrepresented, in PAK/pJN-FcsR (corresponding to genes in PAO1: Aer, CtpH, PctC, CttP, PA4633, PA2654, PA2867, PA1608, PA2788, PA2652, PA4915, PctA and PctB), and a similar proteomic profile was obtained for the expression of inactive FcsR. Flagellar motility and chemotaxis are closely related processes that required the ordered expression of a great number of genes. In *P. aeruginosa*, the flagella regulon includes 50 genes codifying structural and regulatory flagella proteins, plus chemotaxis system components, that can be classified in four classes according to its temporal expression and the transcriptional regulator involved³⁸. Our results show that proteins from Class II to IV were absent in the proteome analysis of PAK/pJN-FcsR and PAK/pJN-FcsRE₆₀A.

Finally, we showed that proteins that form part of the TTSS are also underrepresented, and that this secretion
system is not induced by calcium depletion in PAK/pJN-FcsR or PAK/pJN-FcsRE₆₀A. This multi-protein secre-
tion apparatus is evolutionarily related to flagella and plays key roles in cytotoxicity³⁹, ⁴⁰. These results concur with
previous reports showing that a strain overexpressing this phosphodiesterase is defective in TTSS-mediated cyto-
toxicity¹³. Along with the down-regulation of many flagellar, chemotaxis and TTSS proteins in PAK/pJN-FcsR, we
have observed that the levels of several proteins increased in the same strain, including the filamentous haemag-
glutinin adhesion CdrA. CdrA is a protein involved in adhesion and biofilm formation in *P. aeruginosa* whose
expression is up-regulated by c-di-GMP, while low levels of this second messenger leads to adhesin release into
the extracellular medium and biofilm disaggregation⁶, ¹⁰, ⁴¹⁻⁴³.

The coordinated regulation of flagella, CdrA and Pel proteins synthesis through FleQ has been previously reported<sup>10, 41, 44, 45</sup>. FleQ is a c-di-GMP responsive transcription factor that activates the expression of flagella genes and at the same time negatively regulates genes of Pel and Cdr operons. Thus, FleQ is central to the proteomic changes reported here. The expression of FcsR or FcsRE<sub>60</sub>A leaded to increased levels of proteins repressed by FleQ and concomitantly decreased levels of proteins whose expression is activated by FleQ. FleQ expression is regulated by another global transcriptional regulator, Vfr, that directly or indirectly coordinates the expression of more than 200 genes relying on c-AMP dependant and independent mechanisms<sup>46-48</sup>. Among them, Vfr activates TTSS gene expression and represses FleQ transcription<sup>48, 49</sup>. Our proteomics results indicate that the biosynthetic pathways regulated by FleQ and Vfr are altered in PAK/pJN-FcsR and PAK/pJN-FcsR-E<sub>60</sub>A. Indeed, the phenotypes and proteomic profiles described here are reminiscent of the ones observed for Δvfr and ΔfleQ<sup>46, 48, 50</sup>. However, both Vfr (Cyclic AMP receptor protein, protein identifier: S0JAB5) and FleQ (ATPasse AAA with identifier S0IAH2) were detected with low number of spectra in PAK and PAK/pJN-FcsR proteomic analysis, and no statistically significant differences were observed between strains (Supplementary Table S1). The presence of these proteins in the membrane enriched fractions mainly reflects the detection of low abundant cytosolic contaminants using a very sensitive method for protein identification. Under these experimental conditions it is not possible to conclude if the expression levels between strains are different. Thus, the direct interference of FcsR (and FcsRE<sub>60</sub>A) with the expression and/or function of these transcription factors is an interesting hypothesis that needs to be tested.

Together with FcsR, and FcsRE₆₀A, a component of the chemotactic signal transduction system, the protein
ChpA, was highly enriched in membrane fraction. This protein has a complex domain organization, with six
classical and two atypical phosphotransfer domains, in addition to a CheY-like receiver domain⁵¹. To evaluate
if ChpA enrichment in membrane fraction could be part of the molecular mechanism underlying FcsR effects,
we express the EAL phosphodiesterase in a strain lacking ChpA. Our results clearly demonstrate that flagellar
motility phenotype is at least partially reverted in the absence of ChpA, however, understanding the molecular mechanism behind this observation will require further experimental analysis.

The apparent redundancy of enzymes related to c-di-GMP metabolism in bacterial genomes, together with the observation that strains with similar global c-di-GMP levels can generate distinct phenotypes, have led to the hypothesis that specific interactions between enzymes and their effectors would be responsible of specific phenotypic outputs¹³, ⁵², ⁵³. Likewise, the analysis of individual mutants of all DGCs and PDEs codified by *P. aeruginosa* genome revealed that they display different phenotypes¹³, raising the question about the molecular mechanism underlying its action. The data presented here strongly suggest that FcsR has other roles besides its phosphodiesterase activity, regulating the expression of flagella, chemotaxis and TTSS proteins. Several pieces of evidence support the specific, phosphodiesterase-independent effect of FcsR. First, the swimming motility inhibition inversely correlates with the concentration of L-arabinose used, and we observed a marked motility phenotype even when using low L-arabinose concentrations. Second, the proteomics changes observed are related to specific processes that have been already reported to be jointly regulated. Finally, we showed that the absence of a single regulatory His-kinase, ChpA, partially reverts flagellar motility inhibition. Altogether, the results strongly suggest a specific mechanism of action of FcsR.

Recently phosphodiesterase independent roles were reported for two enzymatically active EAL containing proteins from *E. coli*, PdeR and PdeL. The principal role of these trigger phosphodiesterases is to control gene expression through specific interactions with proteins or DNA, while c-di-GMP hydrolysing activity has a secondary role modulating these interactions⁵⁴. PdeR interacts with the MerR-like transcription factor MlrA leading to inhibition of the expression of the biofilm regulator CsgD and its downstream regulated genes involved in curli fibres and biofilm matrix production. Interestingly, a point mutation in the EAL motif renders an inactive enzyme that exerts the same effects on gene expression, but loss its regulation by c-di-GMP⁵⁵. Another example of trigger phosphodiesterase is PdeL that controls transcription by binding directly to DNA⁵⁶. These two enzymes represent a novel type of c-di-GMP-sensing proteins that controls gene expression, and opens the possibility that FcsR may be acting by a similar molecular mechanism, directly or indirectly interfering with master regulators of flagella, chemotaxis and TTSS systems. In the case of FcsR, we could not recognize other protein domains that could be responsible of the phosphodiesterase independent actions. The discovery of the moonlighting phenomenon has largely confronted the one protein-one function dogma. It is now well demonstrated that many bacterial metabolic enzymes have important functions in several processes, including adhesion, virulence and cell signalling, that are independent of its catalytic activity¹⁴, ⁵⁷.

Overall our results consistently show that the EAL domain containing protein FcsR participates in the control of processes that are crucial for *P. aeruginosa* pathogenesis through a mechanism that does not involve phosphodiesterase activity. It is tempting to speculate that the concept of trigger phosphodiesterases introduced for *E. coli* could be a more general feature in c-di-GMP signalling. To unveil the molecular mechanism underlying FcsR's regulation of motility, chemotaxis and type III secretion system expression further experimentation will be required.

## Methods

**Bacterial strains, plasmids and growth conditions.** *Pseudomonas aeruginosa* strains PAK, PAO1, PAO1ΔChpA and its derivatives were grown until stationary phase in Luria-Bertani (LB) broth at 37°C with shaking, and maintained on LB agar plates. Gentamicin (200 µg/mL) was added to the medium when needed. To evaluate type III secretion system functionality strains were cultured overnight in the presence of 0.2% of L-arabinose, diluted 1/100 in fresh LB containing 0.2% of arabinose and 5 mM EGTA and incubated until late exponential phase. The inactive version of PA2133 used in this work, was generated by site-directed mutagenesis by GenScript® services. Codon GAA, coding for the catalytic glutamic acid residue at position 60 was replaced by GCG which codes for alanine. PAK/pJN-FcsR and PAK/pJN-FcsRE₆₀A strains were generated by transforming chemically competent PAK with pJN105 plasmid encoding arabinose-inducible phosphodiesterase PA2133⁶ or its inactive version. The same procedure was used for PAO1 and PAO1ΔChpA strains overexpressing FcsR. L-arabinose was supplemented into the medium at 0.2% (w/v) unless otherwise indicated. PAK-PilA⁻ and PAK-FliC⁻ mutants (kindly provided by J. Engel) were used as type IV pilus and flagellum negative controls respectively.

**Motility assays.** Swimming motility was monitored on 0.3% LB agar plates, as previously described⁵⁸. Briefly, bacteria were staked on the top of the agar plates, and grown for 24 h at 37°C. Plates were imaged and swimming areas were determined using Image J 1.48 s (Wayne Rasband, National Institutes of Health USA). Statistical analysis was performed using ANOVA and Tukey HSD (Honestly Significant Difference) test (p < 0.05).

**Biofilm formation assays.** Pellicle formation was assessed by bacteria inoculation into 2 mL of T-broth without salts and statically grown at room temperature. After 72 h cultures were gently removed, and pellicles were stained with crystal violet 0.1%⁵⁹.

Microtiter dish biofilm formation assay was performed as previously described⁶⁰. Biofilm was stained with crystal violet and quantification was performed by measuring absorbance at 570 nm after destaining with 30% acetic acid. Statistical analysis was performed using ANOVA test followed by Tukey test (p < 0.05).

**Electron Microscopy.** 2µL of stationary phase cultures were loaded on carbon/formvar grids and negative stained with 2% uranyl acetate. Images were acquired with a Jeol (JEM 10-10) microscope operating at 100 V (Transmission Electron Microscopy Facility, Faculty of Sciences, UdelaR).
**Sample preparation for proteomic analysis.** Bacterial cultures grown until stationary phase were centrifuged at 6,000 × g. Exoprotein fraction was obtained from supernatants after filtration (0.2 µm pore size) and concentration using centrifugal filter devices (10 kDa cutoff). Membrane enriched fractions were obtained from pellets resuspended in 0.2 M Tris, 1 M sucrose, 1 mM EDTA pH 8, and treated with 1 mg/mL lysozyme. After bacterial osmotic disruption, membranes were resuspended in 40 mM Tris, 4% CHAPS, 10 mM MgCl₂ pH 8.2, containing 20 µg/mL DNAse and 10 µg/mL RNase. Samples were centrifuged 30 min at 4 °C and 120,000 × g. Pellets containing membrane proteins were resuspended in 0.1 M ice cold sodium carbonate, incubated 1 h at 4 °C and centrifuged 1 h at 120,000 × g. Membrane pellets were resuspended in 7 M urea, 2 M thiourea, 4% CHAPS, 1% ASB and 40 mM dithiothreitol solution containing protease inhibitor (Roche).

**2D-Difference Gel Electrophoresis (DIGE).** Four independent biological replicates of exoproteins of each strain were concentrated using 2-D Clean-Up kit (GE Healthcare) and resuspended in labeling buffer (7 M urea, 2 M thiourea, 4% CHAPS, 30 mM tris pH 8.5). A pool of all samples of each experiment was used as an internal standard. For minimal labeling CyDye DIGE fluors (GE Healthcare) were used. Briefly, 50 µg of each sample were individually labeled with 400 pmol of Cy3 or Cy5. Dye swapping (Cy3 and Cy5) was used to avoid potential bias in labeling efficiencies between samples. The internal standard was labeled with Cy2. Pairs of samples were mixed with an equal amount of internal standard and separated by 2-D electrophoresis. Samples were loaded into 13-cm IPG-strips (pH 3–10) and then were separated by isoelectric focusing (IEF). Disulfide bonds were reduced with dithiothreitol (10 mg/mL) and subsequently alkylated with 25 mg/mL iodoacetamide. The second dimension was performed on hand-cast gels (12% SDS-PAGE). Gels were scanned using a Typhoon™ FLA 9500 scanner (GE Healthcare) using the excitation wavelength and filters settings indicated for each fluorophore. DeCyder 7.2 (GE Healthcare) was used for analysis. Spot co-detection, spot quantification by normalization and ratio calculation were performed using DeCyder Differential In-gel Analysis module (DIA). Gel matching and statistical analyses allowing quantitative comparisons of protein abundance across multiple gels were performed using DeCyder Biological Variation Analysis module (BVA). Student's t-test was used to assign statistical significance. Protein spots differentially expressed (p ≤ 0.05; normalized abundance ratios fold changes 25%) were further analysed by mass spectrometry (MS).

**Protein identification by MALDI-TOF/TOF.** In-gel digestion of selected protein bands or spots was performed overnight at 37 °C by incubation with trypsin (Sequencing grade, Promega). Peptide extraction was performed with 0.1% trifluoroacetic acid (TFA) in 60% acetonitrile as previously described⁶¹. Samples were vacuum-dried, resuspended in 0.1% TFA, and desalted using C18 OMIX tips (Agilent). Peptides were eluted with matrix solution (α-cyano-4-hydroxycinnamic acid in 60% acetonitrile, 0.1% TFA) directly into the MALDI sample plate. Spectra acquisition was performed on a 4800 MALDI TOF/TOF (Abi Sciex) operating in positive reflector mode. Spectra were externally calibrated using a mixture of peptide standards (Applied Biosystems). MS/MS analysis of selected precursor ions was performed. Database searching (NCBI nr 20130706) was performed with Mascot (http://www.matrixscience.com) using the following parameters: unrestricted taxonomy; one trypsin missed cleavage allowed; methionine oxidation as variable modification; carbamidomethylation of cysteine as fixed modification; peptide tolerance of 0.05 Da and a MS/MS tolerance of 0.3 Da. Significant protein scores (p < 0.05) and at least one peptide with significant ions score (p < 0.05) per protein were used as criteria for positive identification.

**LC-MS/MS.** Exoproteins and membrane enriched fractions were run on 1 cm long SDS gels (12% acrylamide). In-gel Cys alkylation was performed by incubation with 10 mM dithiothreitol for 1 h at 56 °C followed by incubation with 55 mM iodoacetamide at room temperature for 45 min prior to in-gel digestion. In gel-digestion and peptide extraction was performed as described above. Tryptic peptides were separated using a nano-HPLC (EASY-nLC 1000, Thermo Scientific) coupled to an LTQ Velos mass spectrometer (Thermo Scientific). Peptide mixtures were injected into an Acclaim® PepMap C18 nano-trap column (75 µm × 2 cm, Thermo Scientific) and separated on a 50 µm × 150 mm C18 Easy spray column (PepMap® RSLC, 2 µm, 100 Å) at a flow rate of 250 nL/min. Peptide elution was achieved with a 100 min gradient from 5% to 55% of mobile phase B (A: 0.1% formic acid; B: 0.1% formic acid in acetonitrile). Online MS analysis was carried out in a data dependent mode (full scan followed by MS/MS of the top 10 *m/z* in each segment) using a dynamic exclusion list (exclusion duration 45 s).

**LC-MS/MS data analysis.** LC-MS/MS data analysis was performed in accordance to the PatternLab for proteomics 4.0 software (http://www patternlabforproteomics.org) data analysis protocol²⁴. The proteome from *P. aeruginosa* strain K was downloaded from UniProt (May 2015) (http://www.uniprot.org). A target-reverse database including the 127 most common contaminants was generated using PatternLab's database generation tool. Thermo raw file were searched against the database using the integrated Comet⁶² search engine (v. 2015.2) with the following parameters: mass tolerance from the measured precursor *m/z* (ppm): 700; enzyme: trypsin, enzyme specificity: fully-specific, missed cleavages: 1; variable modifications: methionine oxidation; fixed modifications: carbamidomethylation of cysteine. Peptide spectrum matches were filtered using PatternLab's Search Engine Processor (SEPPro) module to achieve a list of identifications with less than 1% of false discovery rate (FDR) at the protein level⁶³. Results were post-processed to only accept peptides with six or more residues, proteins with at least two peptide spectrum matches. These last filters led to an FDR at the protein level, to be lower than 1% for all search results. PatternLab's Approximately-area proportional Venn Diagram module was used for pinpointing proteins exclusively identified in one biological condition. For the enriched membrane protein sample, the analysis only considered proteins present in at least three of four replicates of each biological condition but absent in all replicates of the other condition. Likewise, for exoproteins, only those found in at least two of the three replicates in one biological condition but in none of the other. PatternLab's TFold module was used for pinpointing proteins
found in both conditions but having a statistically differential abundance according their spectral counts²⁸. TFold module relies on the Benjamini-Hochberg theoretical FDR estimator to maximize the number of identifications that satisfy both a fold-change cutoff that varies with the t-test p-value as a power law and a stringency criterion that aims to filter out lowly abundant proteins that could produce false positives. Proteins present in five replicates in all conditions for membrane proteins or four replicates in all conditions for exoproteins were considered. Only proteins satisfying fold change and p-value criteria were considered as differentially expressed.

To allow comparison between PAK proteins and the reference strain PAO1, the protein blast tool of *Pseudomonas* genome database was used²³.

**Phosphodiesterase activity assays.** Phosphodiesterase activity assays were performed as previously described¹³. Briefly, 25 mL of *P. aeruginosa* late-exponential-phase cultures were centrifuged and pellets washed and resuspended in buffer containing 50 mM Tris pH 8, 10 mM MgCl₂, 250 mM NaCl, 5 mM mercaptoethanol, 1 µM PMSF and complete EDTA-free protease inhibitor (Roche). Cells were lysed by sonication and c-di-GMP standard (Kerafast) was added to the cell lysates of the different strains under investigation (PAK, PAK/pJN-FcsR and PAK/pJN-FcsRE₆₀A) to a final concentration of 6 µM. c-di-GMP hydrolysis was followed at different time points using reversed-phase HPLC coupled to MALDI TOF/TOF analysis. Samples were evaporated to dryness using a Speed Vac and resuspended in 200 µL of miliQ water. HPLC analysis was performed by injecting 50 µL of sample in an Agilent 1200 HPLC system fitted with a Hypersyl Gold aQ column (150 × 4.6 mm, Thermo) using a 2 min isocratic step at 1% B followed by a 20 min linear gradient up to 20% of B at a flow rate of 1 mL/min (A: 0.1% TFA; B: 0.1% TFA in methanol). Detection was performed at 280 nm. Collected samples were analyzed on a MALDI-TOF/TOF Instrument (4800, ABi Sciex) in reflector and MS/MS mode (matrix solution: α-cyano-4-hydroxycinnamic acid in 60% acetonitrile, 0.1% TFA).

**Flagella isolation.** Flagella isolation was performed according to previous reports⁶⁴. Briefly, *P. aeruginosa* cultures were centrifuged; pellets were resuspended in PBS and blended for 1 min at 4 °C. After centrifugation, the supernatants were recovered and centrifuged at 40,000 × g for 3 h. The pellet containing flagella was resuspended in 150 µL of PBS, separated in 12% SDS-PAGE gels and the main bands were processed and analysed for MALDI-TOF/TOF as described before.

**Western Blot.** For ExoS and PcrV detection, cultures grown in the presence of EGTA as described before were centrifuged and supernatants were concentrated using Amicon centrifuge-filters. Exoproteins of PAK, PAK/pJN-FcsR and PAK/pJN-FcsRE₆₀A strains were separated in SDS PAGE gels, and transferred into low fluorescence PVDF membranes before incubation with 1:2,000 dilution of anti ExoS rabbit polyclonal antibody⁶⁵ (kindly provided by J. Shouguang) or 1:500 dilution of anti PcrV rabbit polyclonal antibody⁶⁶ (kindly provided by A. Rietsch). For FliC detection, whole cell extracts were separated in SDS PAGE gels. Proteins were transferred to PVDF membranes, and incubated with 1:20,000 dilution of anti-flagellin polyclonal antibody⁶⁷ (kindly provided by D. Wozniak). Goat anti rabbit antibody coupled to Cy5 (GE Healthcare) was used as secondary antibody. Detection was performed using a Typhoon FLA 9500 (GE, Healthcare).

**Data availability.** Proteomics data generated in the present work have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository⁶⁸ with the dataset identifier PXD004365.

**References**

1. Bodey, G. P., Bolivar, R., Fainstein, V. & Jadeja, L. Infections caused by Pseudomonas aeruginosa. *Reviews of infectious diseases* **5**, 279–313 (1983).

2. Rudkjobing, V. B. et al. The microorganisms in chronically infected end-stage and non-end-stage cystic fibrosis patients. *FEMS immunology and medical microbiology* **65**, 236–244, doi:10.1111/j.1574-695X.2011.00925.x (2012).

3. Costerton, J. W. Cystic fibrosis pathogenesis and the role of biofilms in persistent infection. *Trends in microbiology* **9**, 50–52 (2001).

4. Costerton, J. W., Stewart, P. S. & Greenberg, E. P. Bacterial biofilms: a common cause of persistent infections. *Science* **284**, 1318–1322 (1999).

5. Boles, B. R. & McCarter, L. L. Vibrio parahaemolyticus scrABC, a novel operon affecting swarming and capsular polysaccharide regulation. *Journal of bacteriology* **184**, 5946–5954 (2002).

6. Hickman, J. W., Tifrea, D. F. & Harwood, C. S. A chemosensory system that regulates biofilm formation through modulation of cyclic diguanylate levels. *Proceedings of the National Academy of Sciences of the United States of America* **102**, 14422–14427, doi:10.1073/pnas.0507170102 (2005).

7. Simm, R., Morr, M., Kader, A., Nimtz, M. & Romling, U. GGDEF and EAL domains inversely regulate cyclic di-GMP levels and transition from sessility to motility. *Molecular microbiology* **53**, 1123–1134, doi:10.1111/j.1365-2958.2004.04206.x (2004).

8. Krasteva, P. V., Giglio, K. M. & Sondermann, H. Sensing the messenger: the diverse ways that bacteria signal through c-di-GMP. *Protein science: a publication of the Protein Society* **21**, 929–948, doi:10.1002/pro.2093 (2012).

9. Schirmer, T. & Jenal, U. Structural and mechanistic determinants of c-di-GMP signalling. *Nature reviews. Microbiology* **7**, 724–735, doi:10.1038/nrmicro2203 (2009).

10. Hickman, J. W. & Harwood, C. S. Identification of FleQ from Pseudomonas aeruginosa as a c-di-GMP-responsive transcription factor. *Molecular microbiology* **69**, 376–389, doi:10.1111/j.1365-2958.2008.06281.x (2008).

11. Sudarsan, N. et al. Riboswitches in eubacteria sense the second messenger cyclic di-GMP. *Science* **321**, 411–413, doi:10.1126/science.1159519 (2008).

12. Galperin, M. Y. A census of membrane-bound and intracellular signal transduction proteins in bacteria: bacterial IQ, extroverts and introverts. *BMC microbiology* **5**, 35, doi:10.1186/1471-2180-5-35 (2005).

13. Kulasakara, H. et al. Analysis of Pseudomonas aeruginosa diguanylate cyclases and phosphodiesterases reveals a role for bis-(3'-5')-cyclic-GMP in virulence. *Proceedings of the National Academy of Sciences of the United States of America* **103**, 2839–2844, doi:10.1073/pnas.0511090103 (2006).
14. Henderson, B. & Martin, A. Bacterial moonlighting proteins and bacterial virulence. *Current topics in microbiology and immunology* **358**, 155–213, doi:10.1007/82_2011_188 (2013).

15. Ryan, R. P., Fouhy, Y., Lucey, J. F. & Dow, J. M. Cyclic di-GMP signaling in bacteria: recent advances and new puzzles. *Journal of bacteriology* **188**, 8327–8334, doi:10.1128/JB.01079-06 (2006).

16. Ha, D. G. & O'Toole, G. A. c-di-GMP and its Effects on Biofilm Formation and Dispersion: a Pseudomonas Aeruginosa Review. *Microbiology spectrum* **3**, MB-0003-2014, doi:10.1128/microbiolspec.MB-0003-2014 (2015).

17. Kulasekara, H. D. et al. A novel two-component system controls the expression of Pseudomonas aeruginosa fimbrial cup genes. *Molecular microbiology* **55**, 368–380, doi:10.1111/j.1365-2958.2004.04402.x (2005).

18. Li, Y. et al. BdlA, DipA and induced dispersion contribute to acute virulence and chronic persistence of Pseudomonas aeruginosa. *PLoS pathogens* **10**, e1004168, doi:10.1371/journal.ppat.1004168 (2014).

19. Matsuyama, B. Y. et al. Mechanistic insights into c-di-GMP-dependent control of the biofilm regulator FleQ from Pseudomonas aeruginosa. *Proceedings of the National Academy of Sciences of the United States of America* **113**, E209–218, doi:10.1073/pnas.1523148113 (2016).

20. Valentini, M. & Filloux, A. Biofilms and c-di-GMP Signaling: Lessons from Pseudomonas aeruginosa and other Bacteria. *The Journal of biological chemistry*, doi:10.1074/jbc.R115.711507 (2016).

21. Yu, N. Y. et al. PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. *Bioinformatics* **26**, 1608–1615, doi:10.1093/bioinformatics/btq249 (2010).

22. Waters, C. M., Lu, W., Rabinowitz, J. D. & Bassler, B. L. Quorum sensing controls biofilm formation in Vibrio cholerae through modulation of cyclic di-GMP levels and repression of vpsT. *Journal of bacteriology* **190**, 2527–2536, doi:10.1128/JB.01756-07 (2008).

23. Winsor, G. L. et al. Pseudomonas Genome Database: improved comparative analysis and population genomics capability for Pseudomonas genomes. *Nucleic acids research* **39**, D596–600, doi:10.1093/nar/gkq869 (2011).

24. Carvalho, P. C. et al. Integrated analysis of shotgun proteomic data with PatternLab for proteomics 4.0. *Nat Protoc* **11**, 102–117, doi:10.1038/nprot.2015.133 (2016).

25. Colvin, K. M. et al. The pel polysaccharide can serve a structural and protective role in the biofilm matrix of Pseudomonas aeruginosa. *PLoS pathogens* **7**, e1001264, doi:10.1371/journal.ppat.1001264 (2011).

26. Shan, Z. et al. Identification of two new genes involved in twitching motility in Pseudomonas aeruginosa. *Microbiology* **150**, 2653–2661, doi:10.1099/mic.0.27131-0 (2004).

27. Whitchurch, C. B. et al. Pseudomonas aeruginosa fimL regulates multiple virulence functions by intersecting with Vfr-modulated pathways. *Molecular microbiology* **55**, 1357–1378, doi:10.1111/j.1365-2958.2005.04479.x (2005).

28. Carvalho, P. C., Yates, J. R. 3rd & Barbosa, V. C. Improving the TFold test for differential shotgun proteomics. *Bioinformatics* **28**, 1652–1654, doi:10.1093/bioinformatics/bts247 (2012).

29. Kuchma, S. L. et al. BifA, a cyclic-Di-GMP phosphodiesterase, inversely regulates biofilm formation and swarming motility by Pseudomonas aeruginosa PA14. *Journal of bacteriology* **189**, 8165–8178, doi:10.1128/JB.00586-07 (2007).

30. Merritt, J. H., Brothers, K. M., Kuchma, S. L. & O'Toole, G. A. SadC reciprocally influences biofilm formation and swarming motility via modulation of exopolysaccharide production and flagellar function. *Journal of bacteriology* **189**, 8154–8164, doi:10.1128/JB.00585-07 (2007).

31. Tan, H. et al. Comprehensive overexpression analysis of cyclic-di-GMP signalling proteins in the phytopathogen Pectobacterium atrosepticum reveals diverse effects on motility and virulence phenotypes. *Microbiology* **160**, 1427–1439, doi:10.1099/mic.0.076828-0 (2014).

32. Chilcott, G. S. & Hughes, K. T. Coupling of flagellar gene expression to flagellar assembly in Salmonella enterica serovar typhimurium and Escherichia coli. *Microbiology and molecular biology reviews: MMBR* **64**, 694–708 (2000).

33. McCaw, M. L., Lykken, G. L., Singh, P. K. & Yahr, T. L. ExsD is a negative regulator of the Pseudomonas aeruginosa type III secretion regulon. *Molecular microbiology* **46**, 1123–1133 (2002).

34. Barraud, N. et al. Nitric oxide signaling in Pseudomonas aeruginosa biofilms mediates phosphodiesterase activity, decreased cyclic di-GMP levels, and enhanced dispersal. *Journal of bacteriology* **191**, 7333–7342, doi:10.1128/JB.00975-09 (2009).

35. Jimenez-Fernandez, A., Lopez-Sanchez, A., Calero, P. & Govantes, F. The c-di-GMP phosphodiesterase BifA regulates biofilm development in Pseudomonas putida. *Environmental microbiology reports* **7**, 78–84 (2015).

36. Newell, P. D., Boyd, C. D., Sondermann, H. & O'Toole, G. A. A c-di-GMP effector system controls cell adhesion by inside-out signaling and surface protein cleavage. *PLoS biology* **9**, e1000587, doi:10.1371/journal.pbio.1000587 (2011).

37. Kato, J., Kim, H. E., Takiguchi, N., Kuroda, A. & Ohtake, H. Pseudomonas aeruginosa as a model microorganism for investigation of chemotactic behaviors in ecosystem. *Journal of bioscience and bioengineering* **106**, 1–7, doi:10.1263/jbb.106.1 (2008).

38. Dasgupta, N. et al. A four-tiered transcriptional regulatory circuit controls flagellar biogenesis in Pseudomonas aeruginosa. *Molecular microbiology* **50**, 809–824 (2003).

39. Galle, M., Carpentier, I. & Beyaert, R. Structure and function of the Type III secretion system of Pseudomonas aeruginosa. *Current protein & peptide science* **13**, 831–842 (2012).

40. Hauser, A. R. The type III secretion system of Pseudomonas aeruginosa: infection by injection. *Nature reviews. Microbiology* **7**, 654–665, doi:10.1038/nrmicro2199 (2009).

41. Borlee, B. R. et al. Pseudomonas aeruginosa uses a cyclic-di-GMP-regulated adhesin to reinforce the biofilm extracellular matrix. *Molecular microbiology* **75**, 827–842, doi:10.1111/j.1365-2958.2009.06991.x (2010).

42. Cooley, R. B. et al. Cyclic Di-GMP-Regulated Periplasmic Proteolysis of a Pseudomonas aeruginosa Type Vb Secretion System Substrate. *Journal of bacteriology* **198**, 66–76, doi:10.1128/JB.00369-15 (2015).

43. Starkey, M. et al. Pseudomonas aeruginosa rugose small-colony variants have adaptations that likely promote persistence in the cystic fibrosis lung. *Journal of bacteriology* **191**, 3492–3503, doi:10.1128/JB.00119-09 (2009).

44. Arora, S. K., Dasgupta, N., Lory, S. & Ramphal, R. Identification of two distinct types of flagellar cap proteins, FliD, in Pseudomonas aeruginosa. *Infection and immunity* **68**, 1474–1479 (2000).

45. Baraquet, C. & Harwood, C. S. Cyclic diguanosine monophosphate represses bacterial flagella synthesis by interacting with the Walker A motif of the enhancer-binding protein FleQ. *Proceedings of the National Academy of Sciences of the United States of America* **110**, 18478–18483, doi:10.1073/pnas.1318972110 (2013).

46. Dasgupta, N., Ferrell, E. P., Kanack, K. J., West, S. E. & Ramphal, R. fleQ, the gene encoding the major flagellar regulator of Pseudomonas aeruginosa, is sigma70 dependent and is downregulated by Vfr, a homolog of Escherichia coli cyclic AMP receptor protein. *Journal of bacteriology* **184**, 5240–5250 (2002).

47. Fuchs, E. L. et al. The Pseudomonas aeruginosa Vfr regulator controls global virulence factor expression through cyclic AMP-dependent and -independent mechanisms. *Journal of bacteriology* **192**, 3553–3564, doi:10.1128/JB.00363-10 (2010).

48. Wolfgang, M. C., Lee, V. T., Gilmore, M. E. & Lory, S. Coordinate regulation of bacterial virulence genes by a novel adenylate cyclase-dependent signaling pathway. *Developmental cell* **4**, 253–263 (2003).

49. Davinic, M., Carty, N. L., Colmer-Hamood, J. A., San Francisco, M. & Hamood, A. N. Role of Vfr in regulating exotoxin A production by Pseudomonas aeruginosa. *Microbiology* **155**, 2265–2273, doi:10.1099/mic.0.028373-0 (2009).

50. Jyot, J., Dasgupta, N. & Ramphal, R. FleQ, the major flagellar gene regulator in Pseudomonas aeruginosa, binds to enhancer sites located either upstream or atypically downstream of the RpoN binding site. *Journal of bacteriology* **184**, 5251–5260 (2002).

51. Whitchurch, C. B. et al. Characterization of a complex chemosensory signal transduction system which controls twitching motility in Pseudomonas aeruginosa. *Molecular microbiology* **52**, 873–893, doi:10.1111/j.1365-2958.2004.04026.x (2004).
52. Merritt, J. H. et al. Specific control of Pseudomonas aeruginosa surface-associated behaviors by two c-di-GMP diguanylate cyclases. *mBio* **1**, doi:10.1128/mBio.00183-10 (2010).

53. Dahlstrom, K. M., Giglio, K. M., Collins, A. J., Sondermann, H. & O'Toole, G. A. Contribution of Physical Interactions to Signaling Specificity between a Diguanylate Cyclase and Its Effector. *mBio* **6**, e01978-01915, doi:10.1128/mBio.01978-15 (2015).

54. Hengge, R. Trigger phosphodiesterases as a novel class of c-di-GMP effector proteins. *Philosophical transactions of the Royal Society of London. Series B, Biological sciences* **371**, doi:10.1098/rstb.2015.0498 (2016).

55. Lindenberg, S., Klauck, G., Pesavento, C., Klauck, E. & Hengge, R. The EAL domain protein YcIR acts as a trigger enzyme in a c-di-GMP signalling cascade in E. coli biofilm control. *The EMBO journal* **32**, 2001–2014, doi:10.1038/emboj.2013.120 (2013).

56. Reinders, A. et al. Expression and Genetic Activation of Cyclic Di-GMP-Specific Phosphodiesterases in Escherichia coli. *Journal of bacteriology* **198**, 448–462, doi:10.1128/JB.00604-15 (2015).

57. Henderson, B. & Martin, A. Bacterial virulence in the moonlight: multitasking bacterial moonlighting proteins are virulence determinants in infectious disease. *Infection and immunity* **79**, 3476–3491, doi:10.1128/IAI.00179-11 (2011).

58. Murray, T. S. & Kazmierczak, B. I. FlhF is required for swimming and swarming in Pseudomonas aeruginosa. *Journal of bacteriology* **188**, 6995–7004, doi:10.1128/JB.00790-06 (2006).

59. Friedman, L. & Kolter, R. Two genetic loci produce distinct carbohydrate-rich structural components of the Pseudomonas aeruginosa biofilm matrix. *Journal of bacteriology* **186**, 4457–4465, doi:10.1128/JB.186.14.4457-4465.2004 (2004).

60. O'Toole, G. A. Microtiter dish biofilm formation assay. *Journal of visualized experiments: JoVE*, doi:10.3791/2437 (2011).

61. Lima, A. et al. Serine/threonine protein kinase PrkA of the human pathogen Listeria monocytogenes: biochemical characterization and identification of interacting partners through proteomic approaches. *J Proteomics* **74**, 1720–1734 (2011).

62. Eng, J. K. et al. A deeper look into Comet-implementation and features. *Journal of the American Society for Mass Spectrometry* **26**, 1865–1874, doi:10.1007/s13361-015-1179-x (2015).

63. Carvalho, P. C. et al. Search engine processor: Filtering and organizing peptide spectrum matches. *Proteomics* **12**, 944–949, doi:10.1002/pmic.201100529 (2012).

64. Totten, P. A. & Lory, S. Characterization of the type a flagellin gene from Pseudomonas aeruginosa PAK. *Journal of bacteriology* **172**, 7188–7199 (1990).

65. Kim, J. et al. Factors triggering type III secretion in Pseudomonas aeruginosa. *Microbiology* **151**, 3575–3587, doi:10.1099/mic.0.28277-0 (2005).

66. Lee, P. C., Stopford, C. M., Svenson, A. G. & Rietsch, A. Control of effector export by the Pseudomonas aeruginosa type III secretion proteins PcrG and PcrV. *Molecular microbiology* **75**, 924–941, doi:10.1111/j.1365-2958.2009.07027.x (2010).

67. Garrett, E. S., Perlegas, D. & Wozniak, D. J. Negative control of flagellum synthesis in Pseudomonas aeruginosa is modulated by the alternative sigma factor AlgT (AlgU). *Journal of bacteriology* **181**, 7401–7404 (1999).

68. Vizcaino, J. A. et al. The PRoteomics iDentifications (PRIDE) database and associated tools: status in 2013. *Nucleic acids research* **41**, D1063–1069, doi:10.1093/nar/gks1262 (2013).

## Acknowledgements

This work was funded by grants from the Agencia Nacional de Investigación e Innovación, Uruguay (FCE_3_2013_1_100344), FOCEM (MERCOSUR Structural Convergence Fund, COF 03/11) and PICT 2014-3047, from Agencia de Promoción Científica y Tecnológica (ANPCyT), Argentina. JR was supported by a fellowship from ANII (POSNAC_2012_1_8950) and from the Comisión Académica de Posgrado, UdelaR, Uruguay. We thank M. Portela for the excellent technical support. We also thank Dr. C. S. Harwood for kindly providing us with the pJN2133 plasmid and Dr. Daniel Wozniak for providing the anti-FliC antibodies. We acknowledge Dr. J. Engel for the strains PAK-FliC⁻ and PAK-PilA⁻, and Dr. A. Rietsch for anti PcrV antibodies. We thank Dr. J. Shouguang for anti-ExoS antibodies and C. Witchurch for PAO1 and PAO1ΔChpA strains

## Author Contributions

J.R. performed experiments and data analysis. A.L. and M.G. supervised data analysis and execution of DIGE and shotgun respectively. J.R.D. assisted in HPLC analysis. A.C. supervised molecular biology experiments and designed the mutant. P.C. supervised bioinformatics data analysis and extensively revised the manuscript. A.K. participated in the experimental designs, discussions and manuscript writing. R.D. supervised the research and wrote the manuscript. All authors contributed to manuscript correction.

## Additional Information

Supplementary information accompanies this paper at doi:10.1038/s41598-017-09926-3

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017