C-di-GMP Regulates Motile to Sessile Transition by Modulating MshA Pili Biogenesis and Near-Surface Motility Behavior in Vibrio cholerae

In many bacteria, including Vibrio cholerae, cyclic dimeric guanosine monophosphate (c-di-GMP) controls the motile to biofilm life style switch. Yet, little is known about how this occurs. In this study, we report that changes in c-di-GMP concentration impact the biosynthesis of the MshA pili, resulting in altered motility and biofilm phenotypes in V. cholerae. Previously, we reported that cdgJ encodes a c-di-GMP phosphodiesterase and a ΔcdgJ mutant has reduced motility and enhanced biofilm formation. Here we show that loss of the genes required for the mannose-sensitive hemagglutinin (MshA) pilus biogenesis restores motility in the ΔcdgJ mutant. Mutations of the predicted ATPase proteins mshE or pilT, responsible for polymerizing and depolymerizing MshA pili, impair near surface motility behavior and initial surface attachment dynamics. A ΔcdgJ mutant has enhanced surface attachment, while the ΔcdgJmshA mutant phenocopies the high motility and low attachment phenotypes observed in a ΔmshA strain. Elevated concentrations of c-di-GMP enhance surface MshA pilus production. MshE, but not PilT binds c-di-GMP directly, establishing a mechanism for c-di-GMP signaling input in MshA pilus production. Collectively, our results suggest that the dynamic nature of the MshA pilus established by the assembly and disassembly of pilin subunits is essential for transition from the motile to sessile lifestyle and that c-di-GMP affects MshA pilus assembly and function through direct interactions with the MshE ATPase.


Introduction
Vibrio cholerae, the causative agent of the human intestinal disease cholera, is a natural inhabitant of aquatic ecosystems [1]. Cholera infection results from consumption of food and water contaminated with V. cholerae. Subsequently the bacteria turn on regulatory networks that facilitate bacterial growth and survival during the infection process [2]. They also activate production of virulence factors including the toxin-coregulated pilus (TCP), essential for intestinal colonization, and the cholera toxin (CT), responsible for production of massive watery diarrhea that results in dissemination of V. cholerae back to aquatic ecosystems [3]. V. cholerae's ability to cause epidemics is tied to its dissemination and survival in aquatic habitats and its transmission to the human host. One critical factor for dissemination, environmental survival and transmission of the pathogen is its ability to form matrix enclosed surface-associated communities termed biofilms [4][5][6]. V. cholerae forms biofilms on the surfaces of phytoplankton and zooplankton [7], and exists in the surface waters of cholera endemic areas in matrixenclosed aggregates thought to arise from biofilm-like populations of V. cholerae present in human stools [5]. Removal of particles >20 μm in diameter from water can reduce cholera incidence by 48% [8]. Additionally, growth in a biofilm induces a hyper-infectious phenotype [9]. Collectively, these studies highlight the importance of the biofilm growth mode in both the intestinal and aquatic phases of the V. cholerae life cycle.
Two cell-surface structures, a single polar flagellum and type IVa mannose-sensitive hemagglutanin pili (MshA), are critical for initial attachment and biofilm formation [11,20]. Type IVa pili have the ability to rapidly extend and retract, contributing to twitching and swarming motility in many bacteria [21]. Though V. cholerae produces the type IVa MshA pilus, twitching motility has not been reported. Genes required for biogenesis of MshA pilus are clustered into 16.7 kb region in V. cholerae chromosome-I and organized into two operons: the first operon harbors 9 genes from mshI-mshF predicted to encode proteins required for assembly and secretion and the second harboring 7 genes from mshB-mshQ encoding pilus structural components mshB-Q [22]. The pilus is comprised of repeats of the major pilin, MshA. Predicted function of the proteins encoded by MSHA gene locus are included in Table 1.
MshA pili and flagellum are also crucial for two distinct near-surface motility trajectories of V. cholerae: 'roaming' and 'orbiting', [10]. Low curvature 'roaming' trajectories meander over large distances and result from weak MshA-surface interactions. In contrast, orbiting trajectories repeatedly trace out tight circular tracks over the same region, and are the result of strong MshA-surface interactions. Cells that attach to the surface come from the orbiting subpopulation, while the roaming subpopulation pass over the surface without attaching. That orbiting motility is ablated when a mannose derivative is added to the medium to saturate MshA pilus binding further indicates that interactions between MshA pili with the surface are important [10].
The second messenger cyclic dimeric guanosine monophosphate (c-di-GMP) is an important promoter of the switch from motile planktonic growth mode to biofilm growth mode [23][24][25]. c-di-GMP synthesis is catalyzed by diguanylate cyclases (DGC) harboring a GGDEF domain while degradation is catalyzed by phosphodiesterases (PDE) harboring an EAL or HD-GYP domain. Subsequently, c-di-GMP is sensed by different classes of receptor proteins or RNAs and thereby converted to specific phenotypic outputs affecting motility, biofilm  [26,27]. In Pseudomonas aeruginosa and V. cholerae, c-di-GMP inhibits motility by repressing transcription of flagellar genes through the AAA+ ATPase enhancer binding class of c-di-GMP receptor and transcriptional regulators FleQ [28] and FlrA [29], respectively. In addition to the regulation of flagellum production and activity by c-di-GMP, there are reports of c-di-GMP regulating the assembly and activity of Type IV pili. In Klebsiella pneumoniae, c-di-GMP is bound by PilZ class of c-di-GMP receptor protein MrkH, which upregulates the transcription of the fimbrial subunit mrkA [30][31][32]. In P. aeruginosa, the degenerate GGDEF-EAL domain class of c-di-GMP receptor FimX modulates Type IV pili production in an intracellular c-di-GMP concentration-dependent manner [33,34]. The V. cholerae genome encodes 31 GGDEF domain, 12 EAL domain and 10 dual GGDEF/ EAL domain proteins [35]. Systematic analysis of in-frame deletion mutants of all V. cholerae genes encoding proteins with GGDEF and/or EAL domains for motility phenotypes revealed that four DGCs (CdgH, CdgK, CdgL, and CdgD) and two PDEs (CdgJ and RocS) affect motility in an LB soft agar motility assay [36]. Though deletion of the PDE cdgJ affected motility, no difference in intracellular c-di-GMP concentration was observed between the ΔcdgJ mutant and WT [36]. This measurement was conducted from a population, so there may be subcellular localized differences or population differences that affect motility. The molecular mechanism of c-di-GMP mediated motility repression and contribution of these DGCs and PDEs to switch from motile to surface-associated lifestyle remains elusive. In this study, we demonstrate that c-di-GMP inversely regulates motility and biofilm formation through direct regulation of the assembly and activity of the MshA pilus. Swimming motility is impaired in strains lacking the phosphodiesterase cdgJ, and disruption of the assembly or disassembly of the MshA pilus restores motility to WT levels by reducing the interactions with surfaces. Quantitative measurements indicate that c-di-GMP leads to increased production of MshA pili, which in turn bind surfaces and reduce motility. We demonstrate that the ATPase responsible for pilus polymerization, MshE, functions as a c-di-GMP receptor thereby providing an input for the c-di-GMP signal into the assembly of the MshA pilus. Collectively, this study elucidates how type IV pili and swimming motility are regulated by c-di-GMP in V. cholerae by presenting the first characterization of the complex involved in the assembly and disassembly of the MshA pilus and how c-di-GMP regulates the production and function of this complex.

MshA pilus impacts CdgJ mediated motility
CdgJ is a PDE and a cdgJ mutant displays a decrease in motility, enhanced VPS production, and increased biofilm formation compared to WT [36]. To begin investigating the mechanism by which CdgJ impacts motility, we performed transposon mutagenesis in a cdgJ deletion mutant (ΔcdgJ) and screened the resulting mutants for enhanced motility phenotype using LB soft agar motility assay. We screened 7054 transposon mutants and identified 42 extragenic suppressor mutants with increased motility and mapped the transposon insertions site to 22 different genes ( Table 1, Fig 1). As previously reported, we found that mutations in the DGC encoding genes cdgH and cdgK in a ΔcdgJ strain enhance motility [36]. We also found that insertion in the gene encoding a pilus retraction motor, PilT, and insertions into different genes predicted to be required in mannose-sensitive hemagglutinin type IV pilus (MshA) biogenesis restored swimming motility.
To further investigate the suppression of a ΔcdgJ motility phenotype, we generated in-frame deletions of several of the genes identified in the transposon screen in wild-type and ΔcdgJ strains and analyzed the mutants for motility phenotype using LB soft agar motility assay. We focused on mshA (encoding major pilin subunit), mshE (encoding putative polymerizing ATPase), and pilT (encoding putative depolymerizing ATPase) as they are crucial for production and function of MshA pili. As previously reported, ΔcdgJ has a significant motility defect compared to the parental WT strain [36] (Fig 2A). Deletion of mshA, mshE, and pilT in a ΔcdgJ background restored motility similar to the WT strain, confirming that those mutations mediate suppression of the flagellar motility defect in the ΔcdgJ mutant. Deletion of pilU (predicted to encode a second copy of putative depolymerizing ATPase) in ΔcdgJ had no effect on the motility compared to the parental ΔcdgJ mutant. Mutants of mshA, mshE, and pilT in a WT background were assayed for motility to determine if their enhanced motility in the ΔcdgJ background is dependent on this mutation, or if this could be a case of bypass suppression. The ΔmshA, ΔmshE, and ΔpilT mutants exhibited enhanced motility compared to the WT strain (Fig 2A). The ΔpilU strain had a similar motility phenotype to the WT strain, suggesting that the functions or expression profiles of pilT and pilU are different. These data demonstrate that deletion of mshA, mshE, and pilT enhances motility regardless of the presence of a wild-type copy of cdgJ. However, we could not rule out the possibility that cdgJ could directly or indirectly control the production or function of the MshA pilus.
To determine if mshA-mediated suppression of the flagellar motility phenotype is specific to the cdgJ mutation or if it occurs in other PDE deletion backgrounds, we mutated mshA in a ΔrocS strain. RocS has both GGDEF and EAL domains and is predicted to function mainly as a PDE as ΔrocS mutants have reduced motility along with enhanced VPS production and biofilm formation [36][37][38]. We generated a ΔrocSmshA double mutant and determined that this mutant exhibited a wild-type motility phenotype (S1 Fig). These findings suggest that MshA negatively impacts V. cholerae flagellar motility and is involved in general c-di-GMP mediated repression of motility.
To evaluate further the ability of MshA to repress motility, we analyzed the effect of the expression of a WT copy of mshA provided in trans in an expression plasmid with an IPTGinducible promoter. Motility assays confirmed that expression of mshA, upon induction with IPTG, significantly reduced motility in a ΔmshA strain ( Fig 2B). Additionally, expression of mshA in the WT strain reduced motility, suggesting that overproduction of MshA impairs motility. IPTG had no effect on motility in strains harboring an empty vector control.

mshA, mshE, and pilT mutants are deficient in biofilm formation
Since the MshA pilus is critical for initial stages of surface attachment and subsequent biofilm formation [10,11,20], we hypothesized that the ΔmshE and ΔpilT mutants would phenocopy the reduced biofilm phenotype of a ΔmshA mutant. Biofilms of these mutants were grown using a flow cell system, imaged using confocal microscopy, and analyzed using the COMSTAT image analysis software package to evaluate biofilm structural properties. As expected, the ΔmshA mutant attached poorly to the substrate and formed biofilms with low biomass (Fig 3A and 3B). The ΔmshE and ΔpilT mutants grew biofilms with significantly less biomass, thickness, and substrate coverage than WT. COMSTAT analysis revealed that while biofilm biomass of ΔpilT and ΔmshA was similar, the biomass of ΔmshE was significantly less than the ΔpilT strain ( Fig 3B). Additionally, surface coverage of ΔpilT mutant was greater than that of ΔmshA and ΔmshE. These differences suggest that though the ΔmshA, ΔmshE, and ΔpilT mutants are deficient at forming biofilms, there are subtle differences in the phenotypes of the strains. The ΔpilU mutant produced biofilms that were indistinguishable from WT.
Since our initial interest in the MshA pilus was sparked by the discovery that mutations in pilus genes can suppress the motility defect in mutants of the PDE cdgJ, we determined whether Mutations in mshA, mshE, and pilT increase motility. A. The diameters of migration zones of the WT and mutants were measured after 16 h of incubation at 30°C on LB soft agar motility plates and normalized to the motility of the WT strain. Four biological replicates were performed in quadruplicate. Statistical significance determined with Oneway ANOVA followed by Dunnett's multiple comparison test comparing to the WT strain. (ns-not significant, *** p 0.001) B. Expression of mshA reduces motility in soft agar. The diameters of migration zones were measured after 16 hours of growth at 30°C with or without IPTG induction and normalized to the motility of the WT strain. Two biological replicates were performed in triplicate. Induced strains compared to uninduced strains with Student's t-test. (ns-not significant, *** p 0.001). the MshA pilus mutations were epistatic to the cdgJ mutation. As previously reported, the ΔcdgJ strain produced thicker biofilms than the WT strain ( Fig 3C and 3D) [36]. The ΔcdgJmshA, ΔcdgJmshE, and ΔcdgJpilT strains formed biofilms with significantly less biomass, thickness, and surface coverage than the ΔcdgJ strain. These findings are consistent with the hypothesis that these proteins are involved in formation of the MshA pilus and that a functional MshA pilus is required for biofilm formation, surface attachment, and the inhibition of flagellar motility. Biofilms formed by the ΔcdgJmshE mutant had a significant reduction in biomass, thickness, and surface coverage compared to the ΔcdgJmshA or ΔcdgJpilT strains. mshA, mshE, and pilT mutants are defective in surface attachment and near surface motility In addition to using bulk differences in biofilm formation as an index of the transition between motile and sessile behavior, we also directly monitored surface attachment of V. cholerae with single cell resolution using high-speed microscopy and cell tracking (at 5 ms frame rate) to elucidate how genes involved in MshA pilus production impact microscopic outcomes such as initial surface attachment. ΔmshA mutants do not attach to the surface in significant numbers, with no attached cells observed in the first 15min after inoculation. In contrast, well over 100 cells of the WT strain attach to the surface over the same time interval (Fig 4A). Using the same metrics, the ΔmshE mutant was also unable to attach to the surface, exhibiting binding profiles that were similar to the ΔmshA strain. The ΔpilT strain was able to attach to the surface more than the ΔmshA; however, it is unable to attach with the same efficiency as the WT strain. This suggests the following hierarchy of behavioral categories for comparison with non-WT backgrounds: the strong binding strain (WT), the intermediate binding strain (ΔpilT), and the weak binding strains (ΔmshA and ΔmshE). To dissect the origin of this aggregate statistical behavior of surface attachment, we needed more single-cell metrics for bacterial behavior near a surface. We examined the temporal aspects of single bacterium interactions with the surface in the form of single cell residence time (number of seconds that each stationary cell remained associated with the surface, Fig 4B). The WT strain formed prolonged associations with the surface, with a mean residence time of 2.6 seconds and a maximum of 78.9 seconds. The ΔmshA, ΔmshE, and ΔpilT mutant strains demonstrated more transient interactions with the surface, with reduced mean residence times compared to WT (0.69, 0.58, and 0.58 seconds, respectively). This is most evident in the inset of Fig 4B, where the entire adherent populations of these mutants have residence times of less than 6 sec, while the WT residence times extend out to nearly 80 seconds. mshA, mshE, and pilT mutants are defective in near surface motility We used high-speed microscopy and near-surface cell tracking to record the trajectories of cells within one micrometer of the coverslip surface in a microscopy chamber ( Fig 5). All strains with flagella are capable of swimming motility in 3D, and can exhibit trajectories that come in and out of focus. Consistent with previous reports, the WT strain exhibits orbiting and roaming behavior with regards to near surface motility [10]. All of the surface attached cells come from the orbiting subpopulation. The ΔmshA mutant does not show orbiting or roaming behavior, consistent with the model that MshA pili-surface interactions are responsible for these near-surface motility phenotypes. Moreover, also consistent with the model, the mutant shows greatly reduced surface attachment (Fig 4A) [10]. The ΔmshE strain tracks phenocopy ΔmshA, exhibiting predominantly a swimming phenotype with little observable attachment. The ΔpilT mutant has an intermediate phenotype; a few cells appear to exhibit behavior similar to WT orbiting, but has greatly reduced attachment compared to WT. These data are consistent with the biofilm and attachment data described in Figs 2 and 3.
mshA, mshE and pilT mutations are epistatic to a cdgJ mutation with regards to surface attachment To investigate the role of the PDE CdgJ on initial surface attachment, we observed mshA, mshE, pilT mutants in a ΔcdgJ background with high-speed microscopy using similar experiments. We determined that the ΔcdgJ mutant exhibits strong attachment to the surface, with a sharp increase in the number of attached cells as a function of time during the initial few minutes compared to WT. (Fig 4C, time = 1). This strain rapidly associated with the surface, with nearly all cells binding within the first two minutes of observation. As predicted from the biofilm and motility data (Figs 1 and 2), the ΔcdgJmshA, ΔcdgJmshE, and ΔcdgJpilT strains exhibited reduced surface attachment compared to the parental ΔcdgJ strain and WT ( Fig 4C). As in the WT background, the ΔcdgJpilT strain exhibited an intermediate level of attachment and the ΔcdgJmshA and ΔcdgJmshE strains were poor at attachment. The ΔcdgJmshA, ΔcdgJmshE, and ΔcdgJpilT strains exhibited short residence times with mean values of 0.43, 0.53, and 0.57 s, respectively ( Fig 4D). By contrast, the ΔcdgJ strain exhibits much longer mean surface residence times (3.59 s with and a maximum of 75.8 seconds), which surpass even WT (Fig 4D, red vs black bars). These longer residence times indicate a strong tendency of the ΔcdgJ strain to associate with surfaces, which is also evident in the rapid decrease in density of tracks of the ΔcdgJ over time (Fig 5): since only motile cells are displayed in each image, the density of tracks diminishes over time in the ΔcdgJ strain due to the increasing proportion of adherent, nonmotile cells. These results demonstrate the ΔcdgJ strain adheres to surfaces more rapidly than WT, and that mshA is epistatic to cdgJ.
Production of an extracellular MshA pilus is altered in ΔmshE, ΔpilT, and ΔcdgJ The motility and biofilm phenotypes of the ΔmshE and ΔpilT mutants, combined with homology to known Type IV pilus motor proteins suggests that these genes are involved in the production of a functional MshA pilus. MshE shares 75% amino acid similarity with the Type IV extension ATPase PilB of P. aeruginosa and 77% similarity to the Type II extension ATPase EpsE of V. cholerae, suggesting that MshE is the extension ATPase of the MshA pilus (S2 and S3 Figs). We investigated the role of these genes in the production of MshA pili using a surface MshA pilin ELISA (Fig 6A), which detects only assembled pili. ΔmshA, ΔmshE, and ΔpilT mutants produced significantly less surface MshA pili than the WT strain. This supports the hypothesis that MshE and PilT are involved in the production of a functional MshA pilus. The ΔpilT strain produced significantly more surface pili than ΔmshA or ΔmshE, correlating with the intermediate phenotype of this strain observed in biofilm and near surface motility assays (Fig 3A and 3A). Surface MshA pilus production was determined in deletion mutants of other genes in the secretory operon with multiple transposon insertions (S4 Fig). Deletion of mshL, which encodes the putative outer membrane pore protein, resulted in no surface pili. While deletion of mshM (predicted to encode an ATPase) or mshN (predicted to encode a tetratricopeptide repeat domain) resulted in similar pilus production to the WT strain, pilus production was increased in a ΔmshI mutant, The specific mechanisms by which lack of these genes results in suppression of a motility defect in a ΔcdgJ strain is yet to be determined. The ΔcdgJ mutant produced significantly more surface MshA pili than WT. Surface MshA pili were reduced in ΔcdgJmshA, ΔcdgJmshE, and ΔcdgJpilT mutants to the level of the ΔmshA strain, indicating that none of these strains could produce a functional pilus. Whole cell western analysis indicated that all of the strains except ΔmshA and ΔcdgJmshA produced similar amounts of MshA protein, suggesting that the altered surface MshA is not a result of altered MshA production ( Fig  6B and 6C). These data support the hypothesis that a functional MshA pilus inhibits motility and enhances surface attachment and biofilm formation. This effect is likely due to the interactions between the pili and surfaces, however, there may be additional mechanisms affecting motility and surface attachment, as discussed below. Additionally, the elevated production of MshA pili by the ΔcdgJ mutant could explain the motility and biofilm phenotypes observed in this mutant.
To confirm further that mshE was responsible for the lack of pili in a ΔmshE mutant, we generated chromosomal replacements at the native locus with either the WT mshE or a sequence encoding a mutation in the Walker A ATPase active site (K329A) (S5 Fig). Pilus production was restored to WT levels when the WT sequence was inserted, however the K329A produced no detectable surface pili. These data confirm that MshE, and specifically an intact ATPase domain, are required for MshA pilus production. We also utilized transmission electron microscopy (TEM) to assess presence of MshA on the cell surface (Fig 7). We determined that WT produces several pili along the cell body (range 2-5) that were about one half to one cell body length. These pili were not present in the mshA mutant, suggesting that the pili observed are in fact MshA pili. Similarly, the ΔmshE strain produced no visible pili, which is consistent with the prediction that MshE is the motor protein responsible for extension of the MshA pilus. These images reveal that there were no differences observed between pili produced by WT and the ΔpilT, ΔpilU and ΔcdgJ mutants. Due to the fragile nature of MshA pili, quantitative measurements were not possible with TEM though these images observing the presence or absence of MshA pili support the quantitative measurements observed by ELISA.

MshE binds c-di-GMP
The pilus motor proteins MshE and PilT contain Walker A ATPase domains, which are utilized to energize the assembly and disassembly of pili. Baraquet et al. demonstrated that the activity of the P. aeruginosa regulator FleQ is modulated by binding c-di-GMP at its Walker A site [39]. We hypothesized that one or both of the Msh pilus motor proteins could function as a c-di-GMP receptor. Isothermal calorimetry was utilized to investigate the interaction of these proteins with c-di-GMP. We determined that MshE binds c-di-GMP, while PilT and PilU were unable to bind c-di-GMP (Fig 8A). ATPase activity of purified MshE, PilT, and PilU confirmed that these preparations contain functional protein (S6 Fig). VpsT was purified and included as a positive control, as it has been demonstrated to bind c-di-GMP [40]. Fitting the data to a single binding site model indicates that MshE has a slightly lower affinity for c-di-GMP (K = 1.14x10 5 ± 3.24x10 4 M -1 ) than VpsT (K = 9.9x10 4 ± 2.06x10 4 M -1 ). These data suggest that c-di-GMP does affect MshA pilus production through interactions with MshE.
Since MshE, PilT, and PilU contain conserved Walker A ATPase domains, but only MshE bound c-di-GMP, we purified the N terminal domain of MshE (amino acids 1-180) for analysis of the interactions with c-di-GMP. This domain lacks homology with either PilT or PilU and therefore is a likely candidate for binding c-di-GMP. Fluorescence thermal shift assays were utilized to determine the interaction between MshE or MshE-N terminal domain and nucleotide. Purified proteins are stabilized by a bound ligand; therefore the temperature midpoint of unfolding (T m ) of the protein in the presence of ligand is higher than the T m of the protein in buffer [41,42]. We observed that the T m of full-length MshE in buffer was 38.43 ± 0.52°C (Fig 8B). When full-length MshE was incubated with ATP or c-di-GMP, the T m increased to 40.15 ±1.02°C and 42.77±0.70°C, respectively. This indicates that full-length MshE binds both ATP and c-di-GMP. The N-terminal domain of MshE had a T m of 66.36 ± 0.79°C in buffer. There was no increase in T m in the presence of ATP (T m = 62.72 ±0.46°C). In contrast, incubation of the MshE N-terminal domain with c-di-GMP increased the T m to 70.18 ±0.55°C, indicating that this domain binds c-di-GMP (Fig 8C). Additionally, these interactions are specific, as neither the full-length nor the N-terminus of MshE bind to cAMP, c-di-AMP, GTP, or cGMP in thermal shift assays (Fig 8B and 8C).

c-di-GMP promotes MshA pilus production
To further evaluate the link between c-di-GMP production and pilus assembly, we performed a surface pilin ELISA to determine the production of MshA pili over a range of c-di-GMP concentrations (Fig 9). To generate a range of c-di-GMP concentrations, we introduced a construct harboring an IPTG-inducible copy of the DGC VCA0956 (Ptac0956) into the chromosome of the V. cholerae O1 El Tor strain A1552. This strain was grown in the presence of varying concentrations of IPTG, followed by detection of extracellular MshA with the surface pilin ELISA.  We observed that upon induction of the expression of VCA0956 with IPTG, there was an increase in intracellular c-di-GMP (Fig 9, Grey bars). Additionally, the surface pilin ELISA detected increased amounts of extracellular MshA that coincided with the increased c-di-GMP (Fig 9, Black bars

Discussion
High c-di-GMP levels inhibit motility of bacteria and studies have highlighted some of the mechanisms involved. c-di-GMP can repress transcription of flagellar genes, or can act posttranscriptionally to regulate flagellar reversals by interactions with particular flagellar motor proteins or by altering the chemotaxis signal transduction system [26][27][28][29][43][44][45]. Our work revealed a role for c-di-GMP in the regulation of MshA pilus production with effects on near surface motility, motile to sessile transition, and biofilm formation via a post-translational mechanism. c-di-GMP also promotes the assembly and activity of Type IV pili in P. aeruginosa [34,46]. In this example, ectopic expression of a DGC results in elevated amounts of c-di-GMP and increased surface pili. We demonstrate that production of the MshA pilus in V. cholerae is increased in response to high concentrations of c-di-GMP. These similar results in two organisms suggest that c-di-GMP may promote type IV pili assembly and activity as a more general mechanism of pilus regulation than previously identified. Many bacteria rely on type IV pili for motility and attachment, so integration of c-di-GMP in post-translational control of these structures could be a conserved mechanism across species.
The V. cholerae genome encodes 31 GGDEF domain proteins but we have found that only a subset of these impact motility, biofilm formation, or both. We previously reported that four DGCs CdgH, CdgK, CdgL, and CdgD and two PDEs CdgJ and RocS affect motility. Furthermore, a strain lacking all four DGC-encoding genes (ΔcdgDΔcdgHΔcdgKΔcdgL) has a markedly high motility phenotype, suggesting the effect of these proteins on motility is additive [36]. Screening for suppressor mutants of ΔcdgJ that restored swimming motility identified CdgH and CdgK, suggesting that c-di-GMP produced by these DGCs could be a substrate for CdgJ. To test if CdgJ and DGCs that control motility physically interact, we analyzed interactions of CdgJ with CdgH, CdgK, CdgL, and CdgD using the commercially available Bacterial Adenylate Cyclase Two-Hybrid System (BACTH) (Euromedex Strasbourg, France). We did not observe any interaction between DGCs and CdgJ or between CdgJ and MshE or PilT, suggesting that these proteins do not need to be in physical contact, or the interaction is too weak or brief to be detected with this method.
Pili are dynamic structures that are generated by the assembly and disassembly of pilin subunits. The motor proteins responsible for this process have been characterized in many other systems. The two putative motor proteins of interest in this report, MshE and PilT, were identified based on homology to the PilT motor protein in P. aeruginosa. These proteins belong to AAA+ ATPase family proteins and are responsible for energizing the addition and disassembly of pilin subunits. Previous studies revealed presence of a bacterial AAA+ ATPase enhancer binding class of c-di-GMP receptors [28]. Here, we demonstrate that MshE binds to c-di-GMP. This is an important finding, as it establishes that ATPases beyond enhancer binding proteins are also capable of binding c-di-GMP. We note that PilT and PilU, which both harbor AAA+ ATPase domain, are unable to bind to c-di-GMP under the conditions tested. We also showed that N-terminal domain of MshE, which is not present in PilT and PilU is capable of binding to c-di-GMP. Future studies will elucidate the effect of c-di-GMP binding by MshE, as well as the specific mechanisms of this interaction.
We have also characterized several proteins necessary for the production of a functional MshA pilus. Although the function of members of the MshA operons were predicted based on homology to proteins of known function in other bacteria, function of these genes in MshA biogenesis were not characterized [22]. We have demonstrated that MshE is responsible for the assembly of MshA pilin subunits into a functional pilus. These data indicate that the decreased attachment and biofilm phenotypes, as well as the enhanced motility of the mshA, mshE and pilT mutants relies on the dynamic nature of the MshA pilus. If the presence of a pilus could enhance attachment and biofilm formation, a ΔpilT strain would in principle produce biofilms that had similar, or even greater biomass than the WT strain. The observation that the ΔpilT strain phenocopies ΔmshA and ΔmshE, which lack extracellular pili, suggests that both extension and retraction of the pili are critical for normal substrate attachment. The ELISA indicates that the ΔpilT strain produces fewer pili than WT, though the TEM images confirm that there are pili on the surface. Future studies will further investigate the role of PilT in the retraction of the MshA pilus. It is important to note that PilT has already been characterized as the retraction pilus of the ChiRP pilus (chitin-induced competence), suggesting that PilT can function in more than one type IV pilus system [47]. This could be a mechanism for genomic conservation, where one promiscuous retraction ATPase is encoded, while several extension ATPases allow for specificity of the system in regulation. Future studies will address the mechanism of co-regulation of these systems to determine whether pilT is expressed constitutively while the specific extension ATPases are regulated, or if there is overlap in the regulation between the extension and retraction ATPases.
Several studies have investigated how bacteria sense and respond to surfaces, often times by rapidly upregulating production of adhesins and polysaccharides [48][49][50]. This regulation is typically mediated by c-di-GMP [23,28,40,50,51]. Many bacteria also regulate motility in response to surfaces. E. coli uses the resistance to flagellar rotation as a mechanosensor and adapt by adding force-generating motor subunits to the stator complex [43]. This allows the bacterium to adjust the force of flagellar rotation to match the viscosity of the environment. An additional example of "stator swapping" to modulate flagellar force was recently published for P. aeruginosa [44]. In this example, c-di-GMP represses motility by excluding the swarming proficient MotC/D proteins from the stator complex in favor of the swarming deficient MotA/B proteins. P. aeruginosa also utilizes the altered chemotaxis protein WspA as a surface sensor, which results in production of c-di-GMP by the DGC WspR upon interaction with a surface [52,53]. In B. subtilis, production of flagella, and therefore swarming motility, is regulated by Lon-dependent proteolysis of the master regulator of flagellar biosynthesis SwrA upon surface contact [54]; and flagellar function is modulated by EpsE which synergizes exopolysaccharide biosynthesis with flagellar motility by acting as a clutch through interaction with the flagellar protein FliG to limit rotation and therefore motility [55].
This study presents a possible mechanism for how c-di-GMP production affects motility and biofilm formation through modulating MshA pilus production. Both pili and flagella contribute to near surface motility and initial attachment [10]. Besides the generation of near-surface motility modes conducive to surface attachment, it is known that van der Waals forces depend crucially on the extent of surface contact [56]. That V. cholerae has a commalike helicoid shape with smaller surface contact areas implies that adhesive forces between the surface and the cell body will be decreased relative to more cylindrically symmetrical species such as P. aeruginosa for most cell orientations, so adhesive contributions from appendages like MSHA to 'anchor' the cell on surface will be comparatively more important. In fact, recent measurements of TFP adhesive forces show that they can be quite strong, in the hundreds of pico-Newton (pN) range, and are surface chemistry dependent, in agreement with our results [57]. That V. cholera select for surfaces that interact with MSHA strongly (and thereby generate 'orbiting' motility) implies that the existence of more functional MSHA induced by c-di-GMP can better anchor a cell mechanically and mitigate against flagellum driven motion. This work further strengthens the notion that there is a mechanistic link between c-di-GMP and initial attachment through modulation of flagellar motility and pilus activity.

Recombinant DNA techniques
DNA manipulations were carried out by standard molecular techniques according to manufacturer's instructions. Restriction and DNA modification enzymes were purchased from New England Biolabs (Ipswitch, MA). Polymerase chain reactions (PCR) were carried out using primers purchased from Bioneer Corporation (Alameda, CA) and the Phusion High-Fidelity PCR kit (New England Biolabs, Ipswitch, MA), unless otherwise noted. Sequences of the primers used in the present study are available upon request. Sequences of constructs were verified by DNA sequencing (UC Berkeley DNA Sequencing Facility, Berkeley, CA).

Generation of a c-di-GMP overproducing strain
A region encompassing the P lacIq-lacI and the P tac promoter elements was amplified from the pMAL-c5x plasmid (New England Biolabs, Ipswitch, MA) by PCR. The amplified product was joined by overlapping PCR to amplicons of~500 bp that correspond to sequences upstream and downstream of the VCA0956 translational start site. The resulting amplicon was cloned into the suicide plasmid pGP704sacB28 and mobilized into Vibrio cholerae A1552 by biparental mating. The selection of double recombinants with the desired insertion of the P lacIq-lacI and P tac promoter elements was performed as described in [14]. Sequences of constructs were verified by DNA sequencing (UC Berkeley DNA Sequencing Facility, Berkeley, CA).

Transposon mutagenesis
To generate a library of transposon mutants, V. cholerae O1 El Tor strain A1552 ΔcdgJ was conjugated with the donor E. coli S-17-l λpir containing the Mariner transposon on the pSC189 backbone [58]. Transconjugants were selected on LB agar containing kanamycin 50μg/ ml and rifampicin 100μg/ml. A total of 7054 mutants were isolated and screened for motility phenotypes on LB soft agar (0.3%) motility plates.

Plate motility assay
Motility plates consist of LB containing 0.3% agar supplemented with 100μM IPTG where appropriate. Plates were poured and allowed to dry at room temperature for 4 h prior to inoculation. Colonies from overnight LB agar plates grown at 30°C were transferred to motility plates and incubated for 16 h at 30°C. Motility diameter was measured and normalized to the average of WT on each plate. Experiments were performed with three biological replicates in triplicate and data were analyzed with a Oneway ANOVA followed by Dunnett's multiple comparison test.

Confocal laser scanning microscopy (CLSM) and flow cell biofilm studies
Inoculation of flow cells was done by diluting overnight-grown cultures to an OD600 of 0.04 and injecting into a μ-Slide VI0.4 (Ibidi, Martinsried, Germany). To inoculate the flow cell surface, bacteria were allowed to adhere at room temperature for 1 h. Flow of 2% v/v LB (0.02% tryptone, 0.01% yeast extract, 1% NaCl; pH 7.5) was initiated at a rate of 7.5 ml/h and continued for 24 h. Confocal images were obtained on a Zeiss LSM 5 PASCAL Laser Scanning Confocal microscope. Images were obtained with a 40X dry objective and were processed using Imaris (Bitplane, Zurich, Switzerland). Quantitative analyses were performed using the COM-STAT software package [59]. Statistical significance was determined using Oneway ANOVA with Dunnett's Multiple Comparison test. Three biological replicates were performed in triplicate. Images presented are from one representative experiment.

High-speed imaging and cell tracking
Bacteria were cultured in full strength Luria-Bertani (LB) broth overnight under shaking at 30°C. Immediately prior to inoculation, cultures were diluted into 2% LB (containing 171 mM NaCl) to an OD600 0.01-0.03. The V. cholerae cells were then injected into a sterile flow-cell containing the same media and imaged immediately.
Imaging was performed with a Phantom V12.1 high-speed camera (Vision Research) col-lecting~20,000 bright-field images at 5 ms resolution with a 100x oil objective on an IX71 Olympus microscope. All movies were recorded at the same frame rate, for the same duration, after the first, third and 15th minute post inoculation.
For cell-tracking algorithms and analysis protocol, every frame of a movie was preprocessed in Matlab (Mathworks) by subtracting the background, scaling, smoothing and thresholding. Image processing this way causes the bacteria appear as bright regions. Tracking is done by locating all bright objects that overlap objects in the next frame by combining the two frames into a three-dimensional (3D) matrix and then by locating 3D connected components.
Results are stored in a tree-like data structure with multiple roots; every newly detected bacterium that appears is recorded as a 'root' of the tree. When bacteria interact, they are recorded as a 'node' of the tree; when they depart, they are recorded as a 'leaf'. Each root or node stores the sequence of pixel lists that comprise the bacterium in all frames until the next interaction or detachment event. We measure the instantaneous shape properties of the bacteria using the Matlab regionprops function [10].

Surface Pilin ELISA
Surface pili composed of MshA were quantified using an ELISA based on a previously published protocol [46]. Briefly, overnight culture was diluted 1:100 in fresh LB medium and grown to OD 600 0.5 at 30°C. Cells (125μL) were added to a 96-well plate (Greiner Bio-One, Monroe, NC) and incubated at 30°C for one hour. Cells were fixed with 100μL of methanol for 10 minutes at room temperature, then washed twice with PBS. Samples were blocked in 5% nonfat dry milk and immunoblotted with polyclonal rabbit anti-MshA (1:1000 dilution, gift of J. Zhu) and horseradish peroxidase (HRP)-conjugated secondary antibody (Santa Cruz Biotechnology, Santa Cruz, CA). After three washes in PBS, 100μL of TMB (eBioscience, San Diego, CA) was added and incubated for 30 minutes at room temperature followed by the addition of 100μL of 2N H 2 SO 4 . Absorbance was recorded at 490nm and the samples were normalized to the change in WT. Three biological replicates were assayed in triplicate and statistical significance was determined with a Oneway ANOVA followed by a Dunnett's Multiple Comparison test.

Western blot
Samples were grown to mid-exponential phase (OD 600 0.5) in LB or LB with IPTG. Cells were collected via centrifugation and cell pellets were resuspended in 2% SDS and boiled for 5 minutes. Lysates were cleared via centrifugation and total protein was quantified via BCA assay (Pierce, Rockford, IL). Two hundred μg of protein was separated on a 12% SDS PAGE gel and transferred to a PVDF membrane using a semi-dry transfer apparatus (Bio-Rad, Hercules, CA). Blots were blocked in 5% nonfat dry milk and immunoblotted with polyclonal rabbit anti-MshA (1:2000 dilution, gift of J. Zhu) and horseradish peroxidase (HRP)-conjugated secondary antibody (Santa Cruz Biotechnology, Santa Cruz, CA). Chemiluminescence was detected with the SuperSignal West Pico reagents (Pierce, Rockford, IL) on the ChemiDoc MP Imager (Bio-Rad Hercules, CA). Densitometry was performed using the Image Lab software v4.0.1 (Bio-Rad, Hercules, CA) using the band in the WT lane as a reference. Blots were performed in triplicate for densitometry analysis and a representative image is shown.

Transmission Electron Microscopy
Bacteria were prepared for electron microscopy by inoculating a single colony of each V. cholerae strain in LB broth grown overnight at 30°C with shaking at 200 rpm after which each culture was diluted 1:100 in LB broth and allowed to grow similarly to OD 0.4. An aliquot of each culture was diluted to yield an optical density of 0.1-0.2 and then applied to a 300 mesh carbon-coated Formvar grid (Electron Microscopy Sciences, Hatfield, PA). After 2 minutes, each grid was washed five times with deionized water, and negatively stained with 2% (w/v) aqueous uranyl acetate solution for 90 seconds. Imaging was performed with a JEOL JEM-1400 transmission microscope.
Protein purification E. coli BL21 harboring plasmids for gene expression were grown to an OD 600 of 0.4 at 30°C in LB containing 100μg/mL ampicillin. Cultures were shifted to 18°C and IPTG was added to a final concentration of 100μM. 16h post induction, cells were harvested by centrifugation at 10,000 x g for 15 minutes and stored at -80°C.
Cell pellets were resuspended in GST Lysis Buffer (50mM Tris (pH 8.0), 1M NaCl, 0.5% Tween-20 containing PI cocktail tablets (Roche Life Science, Indianapolis, IN). Cells were lysed by sonication and cell lysate was cleared via centrifugation. Cleared lysate was loaded onto GST FPLC column as follows.

Isothermal calorimetry
MshE (18.9μM), PilT (20.2μM), PilU (21.1μM), VpsT (19.5μM) and c-di-GMP (250μM) were prepared in 25mM TrisHCl pH 7.5, 150mM NaCl, and 200μM DTT and degassed prior to analysis. ITC was performed in with a VP-ITC (MicroCal, Northampton, MA) with the following parameters: 3 initial injections of 2μL followed by 40 injections of 10μL spaced at 180 seconds. The data were normalized to a run injecting c-di-GMP into buffer to account for the heat of dilution. Data were processed in Origin v7.0 software (OriginLab, Northampton, MA) and fit to a single site model.

Fluorescence thermal shift
Thermal shift assays were performed as previously described with modifications [41,42]. Briefly, purified protein was added to the reaction to a final concentration of 3μM in the presence or absence of 2mM concentration of the indicated nucleotide in buffer (25Mm TrisHCl pH 7.5, 100mM NaCl, 1:1000 dilution of SYPRO Orange Dye (Invitrogen), and 0.2mM MgCl 2 . A melt curve protocol was run on an Applied Biosystems ViiA7 qPCR instrument. The fluorescence was measured using the ROX reporter with a temperature gradient of 20-95°C in 0.5°C increments at 30 second intervals. Melt curve data were trimmed to three data points after maximum and the data were plotted with Boltzmann model to obtain the temperature midpoint of unfolding (T m ) of the protein in each condition using Prism 5.0 software (GraphPad). The fluorescence baseline of each sample was normalized to the buffer control for visualization purposes. Three biological replicates were assayed in triplicate and statistical significance was determined with a Oneway ANOVA followed by a Dunnett's Multiple Comparison test.

ATPase assay
ATPase activity of purified proteins was determined by measuring the production of inorganic phosphate from ATP using the Enzchek Phosphate Assay Kit (Invitrogen). The standard reaction mixture was prepared with the addition of 2mM MgCl 2 , 10mM KCl, and 1mM DTT. Purified protein in buffer (25mM TrisHCl pH 7.5, 100mM NaCl) was added to the standard reaction mixture to a final concentration of 5μM. After a 10 minute incubation at room temperature, ATP was added to a final concentration of 10mM and reactions were incubated at 22°C for 30 minutes. Production of inorganic phosphate was monitored by reading OD 360 and compared to a standard curve of solutions of KH 2 PO 4 . The data are reported as specific activity (nmol Pi/min/mg of protein). BSA was included as a negative control. Three independent experiments were run in triplicate.
c-di-GMP measurement c-di-GMP extraction was performed as described previously [36]. Briefly, 40 ml of culture grown to OD 600~0 .4 was centrifuged at 3220 x g for 30 min. Cell pellets were allowed to dry briefly then re-suspended in 1 ml extraction solution (40% acetonitrile, 40% methanol, 0.1% formic acid, 19.9% water), and incubated on ice for 5 min. Samples were then centrifuged at 16,100 g for 5 min and 800 μl of supernatant was dried under vacuum and lyophilized. Samples were re-suspended in 50 μl of 184 mM NaCl and analyzed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) on a Thermo-Electron Finnigan LTQ mass spectrometer coupled to a surveyor HPLC (Thermo, Waltham, MA). The Synergin Hydro 4u Fusion-RP 80A column (150 mm x 2.00 mm diameter; 4-μm particle size) (Phenomenex, Torrance, CA) was used for reverse-phase liquid chromatography. Solvent A was 0.1% acetic acid in 10 mM ammonium acetate, solvent B was 0.1% formic acid in methanol. The gradient used was as follows: time (t) = 0-4 min, 98% solvent A, 2% solvent B; t = 10-15 minutes, 5% solvent A, 95% solvent B. The injection volume was 20 μl and the flow rate for chromatography was 200 μl/ minutes.
The amount of c-di-GMP in samples was calculated with a standard curve generated from pure c-di-GMP suspended in 184 mM NaCl (Biolog Life Science Institute, Bremen, Germany). Concentrations used for standard curve generation were 50 nM, 100 nM, 500 nM, 2 μM, 3.5 μM, 5 μM, 7.5 μM, and 10 μM. The assay is linear from 50 nM to 10 μM with an R 2 of 0.999. c-di-GMP levels were normalized to total protein per ml of culture.
To determine protein concentration, 4 ml from each culture was pelleted, the supernatant was removed, and cells were lysed in 1 ml of 2% sodium dodecyl sulfate. Total protein in the samples was estimated with BCA assay (Pierce, Rockford, IL) using bovine serum albumin (BSA) as standards. Each c-di-GMP quantification experiment was performed with four biological replicates. Levels of c-di-GMP were compared to WT with Oneway ANOVA followed by a Dunnett's Multiple Comparison test.
Supporting Information S1 Fig. Motility of ΔrocS is restored to WT by deletion of mshA. The diameters of migration zones of the WT and mutants were measured after 16 h of incubation at 30°C on LB soft agar motility plates and normalized to the motility of the WT strain. Three biological replicates were performed in quadruplicate. Statistical significance determined with Oneway ANOVA followed by Dunnett's multiple comparison test comparing to the WT strain. agar motility plates containing a range of IPTG and normalized to the motility of the WT strain. Three biological replicates were performed in quadruplicate. B. Western blots detect MshA production over a range of IPTG concentrations (0, 6.25, 25, 100, 400μM) from cell pellets using α-MshA antibody. Blots were performed in triplicate, with one representative image included. C. Three-dimensional biofilm structures of the V. cholerae strains formed 24 h post inoculation in a flow cell system. Scale bar represents 40μm. Comstat analysis of two independent experiments in triplicate indicated that Ptac0956 with 400μM IPTG had significantly increased maximum thickness compared to WT (WT 13.42 ±1.50μm, Ptac0956 17.60±1.25μm p 0.01). (TIFF) S8 Fig. MshA pilus production is not affected by increased c-di-GMP in a ΔmshE strain. Expression of the DGC VCA0956 was induced with varying amounts of IPTG. Surface MshA pilin was determined by ELISA. Three biological replicates were tested in triplicate. Results were normalized to MshA production in the WT strain. Surface MshA pili production in the ΔmshE strain was significantly decreased in all conditions compared to WTPtac0956 (Oneway ANOVA, Dunnett's Multiple Comparison Test.) (TIFF) S1 Table. Strains and plasmids used in this study. (DOCX)