Cellular motility provides bacteria with the capacity to actively seek out favorable environments and avoid hazardous situations, thereby facilitating growth and survival in natural habitats[@b1]. Some bacterial species have evolved motility mechanisms that allow cells to move along the direction of their long axis on solid surfaces without the aid of flagella[@b2]. In *Myxococcus xanthus*, a Gram-negative soil bacterium, flagella-independent motility is regulated by two separate adventurous (A) and social (S) systems[@b3][@b4]. S motility was initially described as the coordinated movement of large cell groups as discovered through the observations of the cells on agar surfaces[@b3], and the force-generating mechanism involved in S motility has been shown to be the extension and retraction of type IV pili (TFP)[@b5][@b6], which is mechanistically equivalent to the twitching motility in *Pseudomonas aeruginosa* and *Neisseria gonorrhoeae*[@b7]. It has been further demonstrated that the retraction of TFP is dependent on the interaction of TFP with exopolysaccharides (EPS)[@b8], which is the reason S motility requires cell groups for movement and isolated A^−^S^+^ single-cells are not motile on agar surfaces. S motility is inherently a multi-cell phenomenon, so it is not clear how EPS impacts single-cell motility, or how single-cell motility is modified when more cells are involved.

It is known that *M. xanthus* S-motile (A^−^S^+^) cells are able to move as isolated cells on polystyrene surfaces when they are submerged in a highly viscous medium containing 1% methylcellulose[@b6]. Interestingly, the mutants defective in EPS production are found to perform TFP-dependent motility in this system[@b9], although EPS is absolutely required for *M. xanthus* S motility on agar[@b10][@b11]. It has been proposed that the interactions between TFP and polystyrene surfaces are favored by methylcellulose, which may eliminate the requirement for EPS and enable TFP-dependent single-cell motility[@b9]. Previous studies manually tracked a small number of isolated cells in 1% methylcellulose. EPS^−^ cells, *e.g.*, SW504 (Δ*difA*) cells, exhibited a faster motion than EPS^+^ cells, *e.g.*, wild-type DK1622 cells[@b9], which suggests that EPS may also affect the single-cell S motility in *M. xanthus*. Recent work has shown that the relation between EPS and TFP motility can be complex in *P. aeruginosa*[@b12][@b13][@b14][@b15]; therefore, it is likely that EPS may influence single-cell motility parameters other than cell speed. However, it is not clear whether *M. xanthus* can combine TFP activity with EPS production to produce different motility outcomes. Although tracking of single cells could be illuminating, the ability to track large populations should be useful for the time-resolved analysis of the underlying biological mechanisms of cell motility[@b16].

In this study, we leverage recent advances in the tracking of early *P. aeruginosa* biofilm communities to combine single-cell resolution with large sample populations in the motility analysis of *M. xanthus*[@b12]. Using community tracking algorithms[@b13][@b17][@b18], the motility histories of individual *M. xanthus* cells can be extracted by translating video microscopy movies into searchable databases of cell behavior, and motility patterns can be identified by tracking every cell in the database. Thus, we quantitatively characterized TFP-mediated single-cell motility of *M. xanthus* and correlated the differences in motility pattern to EPS production.

Results
=======

Horizontal *M. xanthus* cells with different amount of secreted EPS exhibit different characteristics in single-cell S motility
-------------------------------------------------------------------------------------------------------------------------------

In a liquid medium containing 1% methylcellulose, the TFP-driven S motility of *M. xanthus* dominates and A motility is not active[@b6][@b9]. For this reason, cells can be tracked in either A^+^ or A^−^ background with comparable results[@b19]. In order to investigate the effects of EPS production on single-cell S motility, isolated cells of wild-type strain DK1622 (EPS^+^), EPS deficient strain SW504 (EPS^−^, Δ*difA*) and EPS overproducing strain DK3088 (EPS^++^, *stk*)[@b20][@b21] were used. Representative behavior of horizontally attached cells in each data set was manually selected and tracked in 1% methylcellulose. The trajectory and velocity profiles for one representative cell of each strain are shown in [Fig. 1](#f1){ref-type="fig"}. During the observation period of 30 minutes, the selected cells exhibited different single-cell motility patterns. Horizontally oriented cells with less EPS moved with higher instantaneous velocity (lower panel of [Fig. 1](#f1){ref-type="fig"}) and lower directional persistence (upper panel of [Fig. 1](#f1){ref-type="fig"}). Moreover, more tethering events with cells in a vertical orientation were observed in the motion of SW504 cell (indicated with dark grey frames in upper panel of [Fig. 1](#f1){ref-type="fig"}). Consistent with the previous observation made from 20 tracked cells, SW504 cells exhibited a rapid motility phenotype compared with that of DK1622 cells[@b9]. To statistically analyze the behaviors of horizontally and vertically oriented cell populations, community tracking algorithms[@b12][@b13] were employed to further investigate large sample populations with single-cell resolution.

EPS production strongly impacts the occurrence of vertically oriented, tethered *M. xanthus* cells in methylcellulose medium
----------------------------------------------------------------------------------------------------------------------------

*M. xanthus* cells in methylcellulose medium exhibit tethering behavior, in which cells attach to a polystyrene surface by the tips of their pili and stand up from the surface[@b6][@b22], but exhibit no lateral movement. The percentage of tethered cells was calculated over every frame in the acquisition, and approximately 900 frames were randomly chosen from the videos (see Methods) and analyzed for each strain (N = 38888 WT cell images, 7113 DK3088 cell images and 11160 SW504 cell images, respectively). As shown in [Fig. 2A](#f2){ref-type="fig"}, the tethering ratio of SW504 (EPS^−^) cells was approximately 3 times higher than that of DK1622 (EPS^+^) cells, whereas DK3088 (EPS^++^) cells showed lower tethering ratio than DK1622 (EPS^+^) cells. Because EPS plays a key role in cell-substratum adhesion[@b23], we measured the adhesiveness of *M. xanthus* cells on polystyrene surfaces in 1% methylcellulose together with their EPS production. As shown in [Fig. 2B](#f2){ref-type="fig"}, cells producing more EPS exhibited stronger attachment on the polystyrene surfaces, which could be attributed to the additional adhesiveness provided by more EPS.

The influence of EPS production on the velocity distribution of single-cell S motility of horizontal *M. xanthus*
-----------------------------------------------------------------------------------------------------------------

It is known that EPS mediates *M. xanthus* cell attachment. Adhesion generally involves forces perpendicular rather than forces parallel to the surface, which are more relevant for surface motility. It is interesting to see how EPS impacts lateral motion along the polystyrene surface. We calculated the velocity distribution histograms of cells of three *M. xanthus* strains by tracking motile cells of each strain over 15 minutes (N = 4395 WT cells, 730 DK3088 cells and 829 SW504 cells, respectively). The observed distributions are complex, but a few comparisons can be made. If we compare the velocities at the peak of each distribution, EPS^−^ cells (SW504) showed slightly faster motion than the EPS^+^ WT (DK1622) cells, whereas EPS^++^ (DK3088) cells exhibited slightly slower motion than WT cells (upper panel, [Fig. 3](#f3){ref-type="fig"}). These velocity differences, although clearly present, are surprisingly small for the three strains, especially given the large differences in surface adhesion noted in [Fig. 2B](#f2){ref-type="fig"}. SW504 cells in particular displayed a broad velocity histogram, in contrast with the more sharply peaked WT and DK3088 histograms, which indicated a tendency to move at a peak velocity and are relatively close in magnitude. Moreover, the tethering percentages of the three strains were different ([Fig. 2A](#f2){ref-type="fig"}). To better determine the motility behavior of horizontal cells (which are more relevant to S-motility), we calculated the velocity histograms after excluding the segments of the tracks around times where the bacteria stopped and moved slower than 1 μm/min, which effectively filtered out the vertically tethered cells, and the fraction of which are significantly different for the EPS^+^ (DK1622) and EPS^++^ (DK3088) cells. One can see that the velocity distribution for DK1622 is peaked at a higher nominal velocity than the corresponding parameter for DK3088 and that the velocity distribution of DK1622 has high velocity tails that are absent in that of DK3088 ([Fig. 3](#f3){ref-type="fig"}). Meanwhile, the EPS^−^ Δ*difA* (SW504) strain has a very broad velocity histogram.

To assess whether the overexpression of EPS has a homogenizing effect on the velocity distribution of *M. xanthus*, we statistically tested the difference between DK1622 and DK3088 velocity distributions by performing an analysis of the Pearson moments of the velocity distributions and calculating their skewness and kurtosis, which are statistical measures of the asymmetry and the 'peakedness' of a given distribution respectively. As shown in [Table 1](#t1){ref-type="table"}, the skewness was positive for all distributions because they were all right-skewed and the velocities were all non-negative. The key difference between the DK1622 and the DK3088 distribution is the higher kurtosis value of the DK3088 distribution, which indicates that the velocity distribution for DK3088 is much more sharply peaked than that for WT (DK1622). To further assess whether underproduction of EPS indeed results in a high velocity subpopulation of *M. xanthus* cells, we statistically tested the difference between DK1622 and SW504 velocity distributions using the Larkin test, an algorithm for assessing bimodality vs. unimodality in a univariate distribution based on the F-test. In this context, a larger value of the bimodality index indicates greater evidence for rejecting the null hypothesis of unimodality. Because the Larkin test is based on the F-test, it allows for a statistical significance evaluation. The bimodality indexes are shown in [Table 1](#t1){ref-type="table"}, and we obtained a *p*-value of just over 0.01 for SW504 and just under 0.01 for moving SW504. This analysis indicates that the velocity distribution of SW504 is quantitatively different than that for DK1622, with the EPS underproducer having a significantly more bimodal distribution.

*M. xanthus* EPS affects the directional persistence of single-cell S motility
------------------------------------------------------------------------------

EPS^−^ cells also showed low preference for moving along their body axis in single-cell motility, whereas EPS^+^ cells preferentially and consistently moved lengthwise along their respective body axes ([Fig. 1](#f1){ref-type="fig"}). We quantified the extent of the sideways motion of cells (N = 4395 WT cells, 730 DK3088 cells and 829 SW504 cells, respectively) by measuring the parameter of angle deviation, defined as the angular difference between the long axis and the moving direction of a cell. Cells of all three strains were more likely to move along their body axis, at low angle deviation ([Fig. 4](#f4){ref-type="fig"}, upper panels). This effect was less pronounced in the SW504 (EPS^−^) strain, which had fewer cells moving at the lowest angle deviation but larger fractions of cells moving at a high angle deviation. The frequency distributions of motions taken at different angle deviations in SW504 was significant different with that of DK1622 and DK3088 (both P \< 0.05 in K-S test), whereas the difference between distributions of DK1622 and DK3088 was not large enough to reach significance at the P = 0.05 level in K-S test (0.05 \< P \< 0.10). Cells moved faster when moving along their body axis and at low angle deviation. Moreover, SW504 (EPS^−^) cells moved at nearly the same velocity while moving at a high angle deviation (perpendicular to their body axis), whereas WT (EPS^+^) and DK3088 (EPS^++^) cells moved faster when moving along their body axis and at a low angle deviation ([Fig. 4](#f4){ref-type="fig"}, lower panels). We further quantified this phenomenon by calculating the ratio between their velocity at low angle deviation and at high angle deviation, which was defined as the velocity falloff ([Fig. 5](#f5){ref-type="fig"}, N = 4395 WT cells, 730 DK3088 cells and 829 SW504 cells, respectively). DK1622 (EPS^+^) cells lost approximately 41% of their velocity when moving perpendicular to their body axis, and DK3088 (EPS^++^) cells were less able to do so. Conversely, the EPS^−^ cells (SW504) showed a very small velocity falloff (less than 18%), meaning that cells lost very little instantaneous velocity when moving sideways.

Discussion
==========

In this study, we investigate the single-cell TFP-driven motility that underlies collective, multi-cell S-motility of *M. xanthus*, without the complicating mutual influence of neighboring cells on a given cell's individual motility. To do so, we measured several motility parameters exhibited by single cells of *M. xanthus* in methylcellulose medium by tracking every cell in a library of microscopy videos with community tracking algorithms[@b13][@b17][@b18]. The average velocity is normally considered an important parameter of kinematics or motility. The data show that the specific influence of EPS on TFP-driven motility is more subtle than previously thought, and impacts a range of motility parameters. The nature of this influence can be seen in the calculated velocity distribution histograms of *M. xanthus* cells. Interestingly, the velocity distributions for DK1622 and DK3088 both peaked at velocity values that are not too different from one another, with DK3088 having a slightly lower peak velocity.

Increased EPS secretion can result in both more surface coverage by EPS and/or thicker deposits of EPS, depend on the wetting properties of EPS on the surface, which depend on both surface and interface energies, as well as van der Waals forces. Simply having more EPS coverage on the surface but at the same thickness of EPS would not lead to a significant slow-down of single-cell motility. However, more EPS secretion in the form of thicker EPS can lead to increased viscous drag. In the present context, data are consistent with more EPS secretion in the form of a slightly thicker wetting layer, but at significantly greater EPS surface coverage. The macromolecular organization of EPS may also contribute to the lack of a simple correlation between surface adhesion and resistance to surface motility. The biggest influence of EPS appears to be the shapes and widths of the velocity distributions. The distribution for DK1622 is asymmetric, with a slowly decaying 'tail' from having more cells at higher velocities; these high velocity tails are absent in the distribution for DK3088, which is more tightly and symmetrically centered around the peak velocity value ([Fig. 3](#f3){ref-type="fig"}). These observations were confirmed by the Pearson moments statistical analysis ([Table 1](#t1){ref-type="table"}), which suggests that EPS may play a role in coordinating and homogenizing multicellular motility. Interestingly, the narrowing of the distribution and decrease of the mean in the overexpressing DK3088 mutant implies that the drag force is velocity-dependent (impacting higher cell velocities more than it does lower cell velocities), consistent with fluid friction rather than solid friction on a surface, which in contrast, depends on the normal force and has no (or only weak) velocity dependence.

In the absence of EPS, which TFP are known to bind to, individual cells are able to move faster in the methylcellulose medium. Moreover, SW504 cells have a broad velocity histogram that appears to be bimodal, with a subpopulation of fast cells that are more than 2.5× faster than the slow population, in contrast with the behavior of WT and DK3088, both of which seemed to be peaked around a single, albeit different, velocity value. This observation is further supported by the analysis of histograms using the Larkin test ([Table 1](#t1){ref-type="table"}). Having no EPS causes the SW504 cell population to bifurcate into a high velocity subpopulation and a low velocity subpopulation, rather than to simply have a broader distribution about a higher mean value, which is the expectation if we had only a lower degree of velocity dependent drag from fluid friction in the EPS underproducer.

We hypothesize that not all TFP-EPS binding events lead to pilus retraction and force generation and that EPS has a homogenizing influence on the cell velocities. Consistent with this, SW504 cells generally move at disparate speeds, and may not be able to coordinate into groups of cells at homogeneous velocities. By comparing all three strains, the trend indicates that greater secretion of EPS progressively eliminates the high velocity subpopulation of motile cells rather than slowing all cells equally. There are a number of mutually compatible explanations for the existence of subpopulations with different average velocities. For the EPS-deficient mutant SW504, a significant subpopulation of cells may have TFP that can bind directly to the 'naked' solid surface directly rather than to the viscoelastic polymeric EPS, which results in more direct force transduction and thereby a faster population of cells. EPS secretion may also lead to intracellular signaling that inhibits high motility, and an optimum level of EPS secretion may homogenize velocity distributions in a manner that is amenable to S-motility.

To test the hypothesis that EPS can homogenize cell velocities, we investigated two additional cell motility parameters, the angle deviation and velocity falloff, which quantitatively characterize the directional persistence of *M. xanthus* cells in methylcellulose medium. All three strains were more likely to move along their body axis ([Fig. 4](#f4){ref-type="fig"}, upper panels), but the SW504 (EPS^−^) strain exhibited a systematically greater tendency to have more cells move at higher angles of deviation away from the cell axis. In addition, it was observed that velocity falls off faster for WT cells moving at finite angle deviations in WT strains than for EPS^−^, which is less sensitive, as seen in the percentage change of the falloff (18% vs. 41%). Together, these observations are consistent with a picture where EPS-TFP binding keeps *M. xanthus* 'on track', with less movement away from the cell axis and a smaller velocity falloff. In *P. aeruginosa*, it has been shown that motile cells deposit a trail of Psl exopolysaccharide on solid surfaces by twitching motility, and the subsequent cells that encounter these trails pull themselves towards the Psl-rich regions by their TFP with higher probability[@b15]. Moreover, in the photosynthetic cyanobacterium *Synechocystis* sp. PCC 6803, which has the ability to perform light-directed motility over solid surfaces, cells at the front of a moving finger-like projection secrete and leave behind an extracellular substance trail that is subsequently followed by other cells, which suggests that the extracellular substance modifies the physical properties of the substrate and leads to enhanced motility by guiding the dynamic spatial organization of multicellular communities[@b24].

The behavior of the vertically tethered *M. xanthus* cells was observed in the methylcellulose medium[@b6], and constitutes direct evidence for bacterial pilus retraction in the context of motility[@b22]. We find that the SW504 (EPS^−^) strain has systematically higher fractions of cells in the vertically tethered orientation. Similar effects are observed for *P. aeruginosa*, where it is found that mutants with decreased levels of exopolysaccharide Psl secretion have a higher fraction of standing cells[@b13]. We hypothesized that *M. xanthus* EPS promotes horizontal adhesion at the expense of the vertical population and that direct interactions between TFP and the bare surface may enhance the probability that cells are vertical tethered. Consistent with this, there is a 5× difference in the standing fraction when SW504 (EPS^−^) and DK3088 (EPS^++^) are compared.

It is interesting to note that *M. xanthus* cells generally stopped moving in the horizontal plane when they were vertically tethered ([Fig. 1](#f1){ref-type="fig"}), in contrast to the behavior of *P. aeruginosa*, where vertically oriented cells move diffusively on the surface and at high speed[@b12][@b13]. We hypothesize that this is due to the localization of active TFP and their associated directions of pull. Vertically tethered *M. xanthus* cells exhibit a simple up-down movement, which implies that the functioning TFP are pulling mostly in the vertical direction. This contrasts with vertical cells of *P. aeruginosa* which can 'walk' laterally via TFP and produces high speed diffusive movement[@b12][@b13], which implies that the TFP can be splayed and pulled with a force component along the surface. Taken together, these results are highly suggestive. TFP in *M. xanthus* pull rigidly in directions that are generally parallel to the cell body, whereas TFP in *P. aeruginosa* can pull in a multiplicity of directions relative to the cell body. This observation is consistent with the fact that TFP in *M. xanthus* drive S motility, where cells follow one another with velocity vectors parallel to their cell bodies, whereas TFP in *P. aeruginosa* drive twitching motility, where cells move in a jittery random fashion with complex relations between velocity vectors and cell body axes.

The *M. xanthus* extracellular matrix (ECM) is mainly composed of carbohydrate substances (EPS) with associated proteins and nucleic acids[@b25][@b26], and has important functions in motility, cell-cell cohesion and fruiting body formation[@b27][@b28][@b29]. It has been hypothesized that *M. xanthus* EPS may decrease the friction when S-motile cells move across a solid surface via the retraction of TFP[@b10]. The recent observation of aperiodic stick-slip movements in *M. xanthus* S motility due to TFP-generated force interacting with EPS indicates that EPS acts as much like a lubricant as it does like a glue[@b14] when multicellular movements are considered. The lubricating function of *M. xanthus* EPS explains the observation that the highly motile cells that actively drive fruiting body formation move in the presence of high concentrations of EPS[@b14][@b28][@b30].

The TFP-driven motility of *M. xanthus* is usually investigated on solid agar medium, in the form of S-motility. To elucidate the dynamic effects of EPS production on single-cell motility without the complexities of coordinated multicellular movement, we exploited the fact that A^−^S^+^ cells are able to move as isolated cells in the submerged liquid medium and tracked *M. xanthus* cells with large sample populations in the current study. Here, we demonstrated that EPS impacts the pattern of TFP-mediated single-cell motility in *M. xanthus*, influencing the velocity profile, tethering probability and directional persistence. Bacterial EPS normally forms viscous aqueous suspensions due to their large molecular mass and complicated interactions with other macromolecules, ions and cells, as well as with some low-molecular-mass solutes[@b31]. The rheological properties of isolated EPS from different *M. xanthus* cells suggest that EPS play fundamental roles in cellular agglutination and clumping in liquid cultures [@b32][@b33]. On solid surface, EPS act as a physical substratum to maintain *M. xanthus* biofilm structure and functions[@b28][@b34]. The phenotypes of the mutants with alterations in different genetic loci responsible for *M. xanthus* EPS production examined in this study suggested that the surface attachment of cells significantly depends on the amount of EPS. Consistent with this observation, greater average adhesion force of rupture events than forces obtained on most other microbial surfaces has been determined in *M. xanthus* by *in vivo* force spectroscopy using atomic force microscopy (AFM), and longer retraction length has been observed on the EPS-overproducing mutant DK3088 (*stk*) cells than for wild-type cells[@b35]. Further investigations will be aimed at obtaining better insight into the mechanism responsible for the coordination of grouped cell S motility by EPS through the balance of its lubricant and adhesive functions.

Materials and Methods
=====================

Bacterial growth conditions and strains
---------------------------------------

The *M. xanthus* strains used in this study were wild-type strain DK1622[@b36], EPS deficient strain SW504 (Δ*difA*)[@b20] and EPS overproducing strain DK3088 (*stk*)[@b21]. *M. xanthus* were grown in CYE medium[@b37] at 32 °C on a rotary shaker at 300 rpm for 24 hr and were harvested by 13,000 × g centrifugation for 5 min.

Methylcellulose assay for motility
----------------------------------

The TFP-mediated single-cell motility of *M. xanthus* was assayed as previously described[@b6][@b9]. The harvested *M. xanthus* cells were diluted in MOPS buffer (10 mM MOPS, 8 mM MgSO~4~, pH 7.6) to 5 × 10^6^ cell/ml, and 5 μl suspension was transferred into 12-well polystyrene plates (Costar, Fisher Scientific, USA) for 10 min. Then, the cells were overlaid with 500 μl of 1% methylcellulose in MOPS buffer, and placed at room temperature in the dark for 1 h before filming.

Microscopic imaging
-------------------

Cell movements were monitored with a Nikon Eclipse TE200 inverted microscope at ×60 or ×40, captured with a SPOT RT740 CCD camera (Diagnostic Instrument, USA) and recorded with SPOT advanced software (Version 4.6, Diagnostic Instrument). Recording was set at 10 s intervals and a TIFF image sequence file was outputted. The image scales were calculated by ImageJ software[@b38] with a calibration scale. Microsoft videos were exported at 50 frame/s with SPOT software, which was 500 times faster than real-time.

Bacteria tracking and motility analysis
---------------------------------------

In our previous work on *M. xanthus*, we used high speed tracking to investigate the spatiotemporal fine structure of motility exhibited by single cells or doublets of cells at short times[@b14]. Here, we use community tracking algorithms and population analysis over many hours, similar to those previously used for *P. aeruginosa*[@b12][@b13]. In total 13 videos of DK1622, 13 videos of DK3088 and 11 videos of SW504 were analyzed, for which the continuous images were taken at 10 s intervals for 30--60 min. We preprocessed the image by inverting it to make the bacteria brighter than the background and then applied a spatial filter to remove noise. Pixels that were local maxima in at least three directions, or were brighter than a predetermined threshold were designated as bacteria pixels[@b18], which were then computationally joined to reconstruct locations of all bacteria. A single position for each bacterium was calculated by averaging *x* and *y* positions of its pixels, and the orientation of each bacterium was calculated from the moments of its backbone distribution[@b17]. To determine the length and width of the bacterium, we first rotated each bacterium by its orientation angle and then calculated the maximum *x* and *y* distances between any pair of pixels. Bacteria positions were joined into tracks with the use of code written in IDL (ITT VIS, White Plains, NY)[@b18]. As a control analysis, representative motility trajectories and velocity profiles were generated by Manual Tracking[@b39] as previously described[@b9].

The five motility parameters listed below were calculated and further analyzed to characterize the TFP-mediated single-cell motility of different *M. xanthus* strains in 1% methylcellulose medium.

1\) Velocity: From the cell trajectories, the instantaneous velocities (*v*~n~) for each frame of each cell were calculated via *v*~n~ = \|X~n + 1~-X~n~\|/10, where X~n~ is the cell position in frame *n*.

2\) Velocity histogram: To calculate a velocity histogram, we first separated out cells that moved with an average velocity of over 1.0 μm/min, which ensured that we measured the velocity of an actively moving cell rather than noise due to low amplitude motions of surface-attached cells about their equilibrium positions with average velocities near zero[@b40]. The remaining velocities were binned in half-micron increments and plotted as a histogram for each strain. The significant difference between DK1622 and DK3088 velocity distributions was analyzed by the Pearson moments of the velocity distributions and by calculating the skewness and kurtosis of the distributions. The significant difference between the DK1622 and SW504 velocity distributions was analyzed using the Larkin test[@b41].

3\) Tethering ratio: To detect bacteria that were tethered to the surface by one end, we used a cutoff based on the length of *M. xanthus* cells. Tethered cells would occasionally orient themselves normal to the surface (in a vertical orientation), thus having an apparent length that was significantly shorter than that of a horizontally oriented motile cell. At the same time, tethered cells would move less than their body length over the course of the acquisition. The tethering ratio was then calculated by the percentage of tethered cells versus total cells as previously described[@b9][@b42].

4\) Angle deviation: We defined the angle deviation of a *M. xanthus* cell to be the angle between the instantaneous velocity and the long axis of the cell body. From the position of a cell in a frame *n*, the direction of the long axis of this cell was determined by the straight line between the two ends of the cell. The direction of motion was determined by the straight line between the leading end of the cell in frame *n* and the same leading end in frame *n* + *1* (10 s interval). At each time point when the cell is moving exactly along its body axis, a value of 0 is assigned; when it moves perpendicular to its body axis, a value of π/2 is assigned. The significant difference between the distributions of motions taken at different angle deviations was analyzed by the Kolmogorov-Smirnov (K-S) test.

5\) Velocity falloff: As *M. xanthus* cells prefer to move in the direction of their body axis[@b2], it is expected that they would move faster when moving at low angles of deviation between the velocity vector and the body axis. For each strain, we can measure the average velocity (*v*~A~) for the time points when the cell is moving along its body axis (to within π/6) and the average velocity (*v*~P~) for the time points when the bacterium is moving perpendicular to its body axis (to within π/6). The velocity falloff (*VF*) was calculated via *VF* = (*v*~A~--*v*~P~)/*v*~A~.

EPS quantification
------------------

EPS production of different strains in CYE liquid was measured using a trypan blue binding assay[@b43]. The relative amounts of EPS were calculated relative to trypan blue bound to the wild-type DK1622 cells. Triplicate experiments were performed.

Cell adhesion assay
-------------------

This assay was adapted from the previously published procedure[@b44]. First, 10 μl *M. xanthus* cells at 1 × 10^8^ cell/ml in MOPS buffer was added into 96-well polystyrene plates (Costar, Fisher Scientific) for 10 min, and overlaid cells with 100 μl of 1% methylcellulose in MOPs buffer. The cells were placed at room temperature in the dark for 1 h. Then, 100 μl 5% glutaraldehyde was added and methylcellulose solution were gently removed. The nonadherent and loosely attached cells were gently washed with two 100 μl aliquots of MOPS buffer. The attached cells were fixed by the addition of 100 μl 5% glutaraldehyde for 20 min at room temperature. Samples were washed with 500 μl MOPS buffer 3 times, stained with 100 μl 0.1% crystal violet for 60 min, and solubilized in 100 μl 10% acetic acid, and the absorbance was measured at 570 nm. Next, 50 μl MOPS buffer was used as blank, and 10 μl cells at 1 × 10^8^ cell/ml was directly stained with 100 μl 0.1% crystal violet for 60 min. The cell pellet collected by 5,000 × g centrifugation for 10 min was dissolved in 100 μl 10% acetic acid, and the absorbance at 570 nm was used to estimate total biomass. Triplicate experiments were performed.

Additional Information
======================

**How to cite this article**: Hu, W. *et al.* Interplay between type IV pili activity and exopolysaccharides secretion controls motility patterns in single cells of *Myxococcus xanthus. Sci. Rep.* **6**, 17790; doi: 10.1038/srep17790 (2016).

Supplementary Material {#S1}
======================

###### Supplementary Information

###### Supplementary Video S1

###### Supplementary Video S2

###### Supplementary Video S3

###### Supplementary Video S4

###### Supplementary Video S5

###### Supplementary Video S6

We thank Drs Larry Shimkets, Mitch Singer and Dale Kaiser for providing bacterial strains, Muhaiminu Hossain and Dr. Zhe Yang for technical help in material preparation, and Drs Xuesong He and Hongwei Pan for helpful discussion. This work was supported by the China National Natural Science Foundation Grant 31570065 and 31370110 (to WH), the United States National Institutes of Health Grant GM54666 (to WS), Shandong Natural Science Foundation Grant ZR2012CM003 (to WH), Shandong Innovation and Transformation of Achievements Grant 2014ZZCX02601 (to WH) and China Postdoctoral Science Foundation Funded Project 10000073311022 (to JW). MG and GCLW acknowledge by the Human Frontiers Science Program and NSF DMR1411329.

**Author Contributions** W.H., M.L.G., G.C.L.W. and W.S. designed the experiments and wrote the paper. W.H., M.L.G., J.W. and C.W. performed the experiments. W.H. and M.L.G. analyzed the data. R.L. and Y.L. contributed to the materials and reagents.

![Representative trajectory and velocity profile for an isolated *M. xanthus* cell in 1% methylcellulose.\
DK1622 (Wt, EPS^+^), DK3088 (*stk*, EPS^++^) or SW504 (Δ*difA*, EPS^−^) cells were placed on polystyrene surfaces submerged in 1% methylcellulose medium and cell movements were recorded by time lapse photography. Continuous images were taken at 10 s intervals for 30 min. The motility and tracks of two poles of one representative isolated cell were analyzed. Data are presented as tracking plots (upper channels, trajectory for 30 min) and as velocity profiles (lower panel). Red and grey lines represent the results of tracking different poles of the same cell, respectively. Dark grey frames in the upper panels indicate the identified tethering events (see [supplementary videos](#S1){ref-type="supplementary-material"}). The instantaneous velocity of cell pole (in lower panel, red and grey lines represent two poles of one cell, respectively) was calculated by measuring the distance from a starting to an ending point within two frames of TIFF image sequences.](srep17790-f1){#f1}

![Tethering ratio and cell adhesiveness of *M. xanthus* cells.\
(**A**) The percentage of tethered vs. total cells of DK1622 (Wt), DK3088 (*stk*) or SW504 (Δ*difA*) were calculated on polystyrene surfaces in 1% methylcellulose. The means ± SD are plotted, Student's *t*-test was employed for statistical analysis and \*\*P \< 0.01. (**B**) Quantitative analysis (grey columns) of EPS production of strain DK1622, DK3088 and SW504. Values for all strains were normalized to the wild-type DK1622. Attachment percentages (black columns) of DK1622, DK3088 and SW504 cells on a polystyrene surface in 1% methylcellulose were measured with the cell adhesion assay described in the *Materials and Methods*. The data represent triplicate experiments, and the means ± SD are plotted.](srep17790-f2){#f2}

![Velocity histograms of DK1622 (Wt), DK3088 (*stk*) and SW504 (Δ*difA*) cells on polystyrene in 1% methylcellulose.\
The upper panels show the histograms of the velocities of the full trajectories of all the motile cells; the lower panels show the trajectories when individual nonmotile (slower than 1 μm/min) segments are excluded.](srep17790-f3){#f3}

![Angle deviation of DK1622 (Wt), DK3088 (*stk*) and SW504 (Δ*difA*) cells on polystyrene in 1% methylcellulose.\
The angle deviation was defined as the angle between the instantaneous velocity and the long axis of the cell body and was measured as described in the *Materials and Methods*. At each time point, when the cell is moving exactly along its body axis, the value of angle deviation is 0; when it moves perpendicular to its body axis, the value of angle deviation is π/2. The upper panels show the fraction of motions taken at different angle deviations, and the lower panels show the velocity of motion at different angle deviation.](srep17790-f4){#f4}

![Velocity falloff of strain DK1622 (Wt), DK3088 (*stk*) and SW504 (Δ*difA*) on polystyrene in 1% methylcellulose.\
The velocity falloffs (mean ± SD) of different *M. xanthus* cells were calculated as described in the *Materials and Methods* and quantified the loss of cell velocity when moving perpendicular to the cell body axis. The WT and *stk* strains moved faster along their body axis, whereas the Δ*difA* strain showed very little change. Student's *t*-test was employed for statistical analysis and \*\*indicates P \< 0.01.](srep17790-f5){#f5}

###### Statistical analysis of velocity histograms shown in [Fig. 3](#f3){ref-type="fig"}.

  **Pearson Moments Analysis**    **Larkin Test**                                             
  ------------------------------ ----------------- ------- --------- ---------------- ------- -----
  DK1622(WT)                           Total        0.42     −0.47      DK1622(WT)     Total   1.3
  Nonstop                              0.44         −0.47   Nonstop        1.2                
  DK3088(*stk*)                        Total         3.4      16      SW504(Δ*difA*)   Total   2.1
  Nonstop                               3.6          15     Nonstop        2.2                

[^1]: These authors contributed equally to this work.
