Interplay of motility and polymer-driven depletion forces in the initial stages of bacterial aggregation †

Motile bacteria are often found in complex, polymer-rich environments in which microbes can aggregate via polymer-induced depletion forces. Bacterial aggregation has many biological implications; it can promote biofilm formation, upregulate virulence factors, and lead to quorum sensing. The steady state aggregation behavior of motile bacteria in polymer solutions has been well studied and shows that stronger depletion forces are required to aggregate motile bacteria as compared with their nonmotile analogs. However, no one has studied whether these same trends hold at the initial stages of aggregation. We use experiments and numerical calculations to investigate the polymer-induced depletion aggregation of motile Escherichia coli in polyethylene glycol solutions on short experimental timescales ( B 10 min). Our work reveals that in the semi-dilute polymer concentration regime and at short timescales, in contrast to what is found at steady state, bacterial motility actually enhances aggregate formation by increasing the collision rate in viscous environments. These unexpected findings have implications for developing models of active matter, and for understanding bacterial aggregation in dynamic, biological environments, where the system may never reach steady state.


Introduction
Bacteria thrive in a wide range of biological and ecological contexts and play important roles in the human gut, soil, wastewater sludge, and other complex environments. In these environments, bacterial motility has implications for biofilm formation and virulence. For example, in the gut, Salmonella typhimurium uses its flagella to burrow through the intestinal mucus layer and penetrate the host epithelium, causing infection. 2 Certain species of Pseudomonas aeruginosa require motility to form biofilms, such as those found on medical devices. 3 Microbial motility is further influenced by polymers, which are abundant in many environments. 4 For example, in the gut, polymers are secreted by the host [5][6][7][8] and dietary fibers are ingested regularly. 1 In wastewater treatment plants, sludge used to collect microbes and particulate waste also contains polymers. 9,10 Polymers are well known to aggregate bacteria, which is important because aggregation precedes biofilm formation, 11 correlates with altered gene expression, 3,12 including antibioticresistance genes, 13 and induces phenotypic changes such as quorum sensing. 14 Polymers can bind to bacteria via specific chemical interactions and cause aggregation through agglutination. These interactions are found in biological settings such as the gut, where mucins, 5 immunoglobulins, 6 and other host-secreted proteins 7 can aggregate bacteria via chemically mediated interactions. Polymers can also aggregate microbes via depletion forces. This mechanism of aggregation does not depend on microbes binding to specific chemical functional groups but is instead only a function of the physical parameters of the polymer (molecular weight (MW), hydrodynamic radius) and bacteria size. [15][16][17] In the presence of non-adsorbing polymers such as polyethylene glycol (PEG), bacteria aggregate through depletion interactions, which occur when two bacteria approach each other at a close enough distance that the polymer is excluded from the space between the bacteria, a region called the depletion zone. 13 The difference in polymer concentration between the depletion zone and the bulk solution results in an osmotic pressure difference that generates an attractive force between bacteria. Because depletion forces depend on the physical properties of polymers and bacteria, these forces can drive bacterial aggregation irrespective of bacterial surface chemistry. In the case of nonmotile bacteria (which are not auto-aggregating) in solutions of non-adsorbing polymers, the only driving forces for aggregation are polymer-induced depletion forces. 17 However, for motile bacteria in solutions of non-adsorbing polymers, there are forces due to the motility of the bacteria and due to depletion forces; the competition between these forces determines the steady-state aggregation behavior. 18 At sufficient polymer concentrations and long time scales, when the system reaches steady state, polymer-induced depletion attractions between bacteria can result in aggregation (i.e., phase separation 18 ). The addition of motility has been found to require significantly stronger depletion attraction to achieve the aggregation as compared with nonmotile bacteria. 18 Active matter is an area of intense research, and the field is currently working on a unified theoretical framework to understand these systems. 19,20 In particular, the aggregation behavior of motile bacteria in polymeric solutions at long time scales is widely studied; 18,21-24 however, to our knowledge there are no published studies on aggregation at short time scales.
The initial stages of aggregate formation are of particular interest because microbial responses to aggregation (e.g., upregulation of quorum sensing and virulence pathways) occur on short time scales (tens of minutes) and would be influenced by the initial stages of aggregation. Furthermore, the initial stages of aggregation are particularly relevant for biological systems where the system may never reach steady state; for example, in the gastrointestinal (GI) tract where food and ingested material are constantly in transit. 1,25 In these environments, aggregate formation is constantly disrupted because of shear, peristaltic contractions, and other forces, making early aggregate formation relevant to these biological systems.
In this study, we investigate how motility influences the polymer-induced depletion aggregation of bacteria at short time scales (t B 10 min). We quantify this experimentally by using confocal fluorescence microscopy to measure the size distribution of bacterial aggregates in PEG solutions with molecular weights and concentrations relevant to the murine small intestine. 1 Furthermore, we develop an understanding of which physical parameters influence the initial formation of these aggregates. We use a physical model for motile bacteria in PEG solutions that focuses on the balance of depletion and swim forces as well as the effective diffusivity of the bacteria.

Bacteria cell culture
Overnight cultures of naturally motile E. coli K12 MG1655 (ATCC 47076) were prepared in liquid lysogeny broth (LB) culture incubated at 35 1C to mid-exponential phase. These cultures were combined to reach the desired cell concentration for the experiment (10 9 cells per mL). Cells were first centrifuged at 4.8 kG for 10 min and then resuspended in motility buffer (MB; 10 mM potassium phosphate buffer, 0.1 mM EDTA, pH 7.0) to stain with SYTO 9 (live, ex/em 480/500 nm). Following staining, cells were centrifuged again to wash out excess stain. To obtain nonmotile E. coli, cells were treated with 0.5% (75 mM) sodium azide in MB after washing.
Confocal microscope imaging and bacterial aggregation in polyethylene glycol (PEG) All images and z-stacks were obtained using a Zeiss LSM 800 confocal fluorescence microscope (488 nm excitation, detection at 490-540 nm). Each stack was 200 Â 200 Â 45 mm in volume and contained about 135 slices. Nonmotile and motile E. coli K12 were prepared and stained using the method described above. PEG solutions were prepared by dissolving four times the overlap concentration into MB (for motile conditions) or by using MB with 0.5% sodium azide (for nonmotile condition). A range of 10 kDa, 100 kDa, and 1 MDa PEG solutions were achieved by serial dilution. A 5 mL aliquot of each respective PEG solution was combined with 0.5 mL E. coli for a final cell concentration of 10 9 cells per mL. We pipetted 2 mL of each combined suspension into an imaging chamber made from SecureSeal imaging spacer (Electron Microscopy Sciences; 0.12 mm depth and 9 mm diameter) and a glass slide, and the top of the chamber was immediately sealed with a #1.5 glass coverslip. A single z-stack of each PEG dilution sample was taken approximately 10 min after the imaging chamber was sealed. Each biological replicate was conducted with a new bacterial cell culture.

FIJI macro imaging and empirical bootstrapping
All imaging analysis was performed as previously described in ref. 1.
Measuring mean-squared displacement (MSD) of E. coli E. coli K12 were cultured and prepared as described above. The final cell concentration was diluted to 5 Â 10 8 cells per mL when added to each PEG solution in MB. A Leica DMI6000 with a Visitech Infinity3 confocal microscope was used to obtain 20 s videos of the cells at about 16 frames per s.
Videos were analyzed using an ImageJ plugin developed by the MOSAIC group for 2D/3D particle tracking using an algorithm developed in ref. 26. At least 1000 bacteria were analyzed per condition. Data output from the ImageJ plugin was further analyzed and the mean square displacement (MSD) was calculated using MATLAB code used in ref. 27 in conjunction with an in-house script. (Script will be provided by request.)

Estimating overlap concentration for PEG
The polymer overlap concentration c P * was estimated using the following relation: 28,29 where MW is the polymer molecular weight in kDa, N Avo is Avogadro's number, and R g (0) is the radius of gyration given in m.

Estimating solution viscosity
We estimated the solution viscosity via a virial expansion: 28 where Z s is the solvent viscosity in Pa s, c P is the polymer mass concentration in kg m À3 , [Z] is the intrinsic viscosity measured to be 452.8 mL g À1 (using literature measurements for PEG 1 MDa), 30 and k H is the Huggins parameter for PEG 1 MDa, approximated to be 0.4. Using this equation, Z at c P = 6.5 mg mL À1 is B7 mPa s compared with B1 mPa s at c P = 0.4 mg mL À1 . Literature measurements report similar values for the zero shear rate viscosity (Z 0 ) of PEG 1 MDa at high concentrations, where Z 0 = 10 mPa s at c P = 5 mg mL À1 (our experiments were conducted our quiescent conditions, i.e., no shear). 31

Results and discussion
Measuring E. coli aggregation at short time scales To understand how motility affects the depletion-driven aggregation of bacteria at short time scales, we measured the initial formation of bacterial aggregates. As a model organism, we used E. coli K12 MG1655. This naturally motile strain of E. coli displays ''run and tumble'' dynamics, and can be rendered nonmotile by treating with 0.5 wt% sodium azide. 32,33 This method of rendering cells nonmotile has been used before in research focusing on the flagellar motility of E. coli. 32 Plating showed that the azide treatment killed nearly all of the cells 1 h after the initial treatment ( Fig. S1, ESI †). To ensure we selected biologically relevant physical parameters for our experiments, we used data from our previous gel permeation chromatography experiments on luminal fluid from the murine small intestine to determine the range of the polymer molecular weights and concentrations. 1 In these previous experiments, the polymers we found in the murine small intestine ranged in size from a few kDa to a few MDa. 1 Therefore, we chose to tune the depletion potential in this work with polymers within this size range. As our test polymer, we chose to use PEG solutions in a motility buffer (MB, 0.1 mM EDTA, 10 mM potassium phosphate, pH 7.0). PEG is a linear, chemically inert polymer 33,34 that is well characterized in inducing depletion forces in passive colloid solutions. [35][36][37][38][39][40][41][42] We used a range of 1 MDa PEG concentrations (0.05-6.5 mg mL À1 ) to adjust the depletion potential and rheology of the solution, which span both the dilute and semi-dilute polymer concentration regimes (the transition between these regimes is denoted by the overlap concentration, c P * = 1.6 mg mL À1 (see calculations in Experimental)). We also measured bacterial aggregation in 10 kDa and 100 kDa PEG near their respective overlap concentrations (c P * = 8.5 mg mL À1 and 85 mg mL À1 for 100 kDa and 10 kDa PEG, respectively). We previously detected polymers of a similar size and concentration in the murine small intestine. 1 A motility buffer control was implemented to confirm that the cells were not auto-aggregating and that the aggregation measured in each sample containing PEG was the result of the PEG in the solution.
To quantify the initial aggregation of bacteria, we measured the volume-weighted average aggregate sizes (N) using fluorescence confocal microscopy (Fig. 1a). After mixing the E. coli with the PEG, the bacterial suspension was placed into an imaging chamber, and sealed with a glass coverslip to eliminate drifting and evaporation effects. Z-stacks of cells in solutions of PEG at various concentrations were obtained after 10 min to focus on the behavior on short timescales ( Fig. 1b and c). Separate experiments were performed for motile and nonmotile cells. Imaging at short timescales also reduces the effects of sedimentation from gravity. 43 Imaging analysis was performed (using an ImageJ pipeline that we developed previously 1 ) to count each object, measure the volumes of each aggregate, and normalize by the singlet volume to obtain the volume-weighted average aggregate size.
In these experiments, we found that nonmotile bacterial aggregation in the presence of 1 MDa PEG (Fig. 2a) was qualitatively consistent with depletion-driven aggregation with similar trends observed at half (Fig. 2b) and double (Fig. 2c) the bacterial concentration. Additionally, we tested the effect of changing PEG MW and found trends that were qualitatively consistent with depletion-driven aggregation (Fig. S2, ESI †); the extent of aggregation generally decreased with MW. Nonmotile E. coli aggregated less in the presence of PEG 100 kDa (Fig. S2b, ESI †) as compared to in the presence of PEG 1 MDa (Fig. 2a) and no aggregation was observed in the presence of PEG 10 kDa. Somewhat counterintuitively, as the 1 MDa PEG concentration increased, aggregate size increased up to a limit and then started decreasing. We previously measured similar aggregation profiles with particles in PEG solutions; 1 the shape of those profiles was at least in part attributed to the increase in solution viscosity, which hindered the Brownian motion of the particles, limiting the inter-particle collisions that initiate aggregation. We suspect that a similar mechanism may be at play for the nonmotile bacteria in this study. The depletion attractions increase as a function of polymer concentration due to the contribution of osmotic pressure, which is mirrored by the increasing aggregate size of the nonmotile bacteria in the dilute PEG concentration regime (eqn (3) and (4)). We observed that as the PEG concentration increased and approached the semi-dilute concentration regime, the aggregate sizes became larger. In the semi-dilute regime, the depletion attractions continued to increase with the increase in PEG concentration, but the range of these attractions decreased with PEG concentration. At PEG concentrations far above the overlap concentration, we observed a decrease in aggregate sizes. The use of sodium azide to render E. coli nonmotile could alter cell surface properties in a manner that affects their depletion-driven aggregation. However, due to their lack of aggregation in the motility buffer control and the aggregation curve being qualitatively similar to bioinert particles in similar polymer solutions, 1 we assume that this effect is minimal. We hypothesized that at higher PEG concentrations, the enhanced viscosity decreased inter-bacterial collisions, thus hindering aggregate formation.
Motile bacteria demonstrated different aggregation trends compared with their nonmotile analogs. In the PEG 1 MDa solutions, we observed no aggregation in the dilute PEG concentration regime, but at c P = 0.8 mg mL À1 , motile E. coli abruptly began to aggregate, and continued to aggregate through the semidilute polymer concentration regime (Fig. 2a). Similar trends were observed when we halved (Fig. 2b) or doubled (Fig. 2c) the bacterial concentrations. In the lower MW PEG solutions, we found minimal aggregation for the motile bacteria in PEG 100 kDa (Fig. S2b, ESI †) and no aggregation for motile bacteria in PEG 10 kDa (Fig. 2a). Because the main physical difference between the motile and nonmotile bacteria is their motion, we hypothesized that the differences in depletion attractions required to aggregate motile versus nonmotile bacteria at short time scales are due to cell motility. Previously, researchers have modeled the steady-state depletion aggregation of motile bacteria by assuming that the swim force produced by bacterial motility directly counteracts the depletion force. 18 We hypothesized that a similar framework could be applied here; the swim force counteracts the depletion force at low PEG concentrations, negating any attractive force to aggregate the cells. Therefore, a stronger depletion force, or higher PEG concentration, is necessary to aggregate the motile bacteria to the same extent as the nonmotile bacteria. Support for this hypothesis was experimentally demonstrated at c P = 0.4 mg mL À1 whereby we observed the nonmotile bacteria aggregate but the motile do not.

Effective potentials describe E. coli aggregation in the dilute polymer concentration regime
To investigate the interplay of swim and depletion forces that give rise to the differences in the observed aggregation behavior between motile and nonmotile bacteria, we began by using effective potentials to describe aggregation in the dilute polymer concentration regime. We focused on explaining the discrepancy in aggregation at c P = 0.4 mg mL À1 , where we observed substantial aggregation of nonmotile bacteria but no aggregation of motile bacteria (Fig. 3a). In The dilute regime, nonmotile bacterial aggregation increased with the PEG polymer concentration from the osmotic pressure contribution in the depletion potential (Fig. 3b). The effective potential between nonmotile bacteria can be described by the Asakura-Oosawa depletion potential: 15,16,44 where U dep is the depletion potential (in J), P P is the osmotic pressure of the polymer solution (in Pa), a is the radius of the bacteria (approximated as a sphere, in m), R P is the characteristic polymer size (in m), and r is the separation distance between two bacteria surfaces (in m). The contribution of polymer concentration to the depletion potential is implicit in the osmotic pressure. PEG behaves as a polymer in good solvent in buffer, 45 and we can use the following crossover equation for the osmotic pressure of a polymer solution which spans the dilute and semidilute polymer concentration regimes: 46,47 where N Avo is Avogadro's number, k is the Boltzmann constant, T is the temperature (in Kelvin), MW is the molecular weight of the polymer (in kDa), c P is the polymer mass concentration (in kg m À3 ), and c P * is the overlap mass concentration (in kg m À3 ). We use the concentration-dependent radius of gyration to estimate the characteristic polymer size: 48,49 where R P (c P ) is the concentration-dependent radius of gyration or the characteristic polymer size (in m), and R g (0) is the radius of gyration (in m). R g (0) was estimated using literature values of the hydrodynamic radius of PEG 30 and the Kirkwood-Riseman relation. 50,51 We estimated R g (0) to be 62.6 nm for 1 MDa PEG, and using this value and the molecular weight of the polymer, we estimated c P * to be 1.6 mg mL À1 (see Experimental for calculation). Combining eqn (3)-(5) gives us an expression for the depletion potential that closely approximates the Asakura-Oosawa potential in the dilute polymer concentration regime and the potential derived by Joanny, Leibler, and De Gennes in the semi-dilute regime. 52 Similar crossover equations for the depletion potential have been previously used to quantitatively describe experimentally observed depletion-driven colloid aggregation in polymer solutions that span the dilute and semi-dilute concentration regimes. 1,53 As the polymer concentration increases, the depletion potential also increases. We observed that the aggregate size of the nonmotile E. coli increases with PEG concentration in the dilute polymer concentration regime, which suggests that depletion forces drive aggregation. We found that motile bacteria do not aggregate at low PEG concentrations until a certain PEG concentration threshold is reached (Fig. 3c). To estimate the effective potential of motile bacteria and the effect of the swim force on the aggregation of motile bacteria, we used a previously established theoretical framework. 18 We began by considering the forces that bacteria experience in solution. This model accounts for the swim force that arises from bacterial motility and the polymer-induced depletion force. The swim force can be described from the ellipsoid approximation to the Stokes-Einstein drag coefficient: 17,54 where Z is the solution viscosity in Pa s (see Experimental for details of estimate), a and b are the lengths of the semi-minor and semimajor axes for E. coli (in m), and V is the speed (in m s À1 ). For E. coli, a and b are approximated to be about 0.5 mm and 2 mm, 55 respectively, and V is assumed to be constant at about 10 mm s À1 . 4 The effective force is then calculated using a force balance on the bacteria, assuming that the swim force directly counteracts the depletion force: where F eff is the effective force (in Newtons) and the depletion force is given by the negative first derivative of U dep (r) with respect to r. To find the effective potential U eff , we integrate eqn (5) with respect to r: The integration constant, U 0 (J), is defined as described previously. 18 In a condition where two bacteria are swimming in the exact opposite direction, there is a separation distance r* where the effective force acting on each bacterium is zero. Beyond this range, the swim force overwhelms the depletion potential, and the effective potential on the bacteria is zero. U 0 is defined such that U eff is zero beyond r*.
In the dilute polymer concentration regime, we observed that the aggregation trends for both motile and nonmotile bacteria were qualitatively consistent with expectations based on the changes of the minima of their respective effective potentials (U eff (r = 0)) at each PEG concentration (Fig. 3d). For nonmotile bacteria, the effective potential consists of only the depletion potential, which increases with polymer concentration. Our observations of the aggregation of nonmotile bacteria were qualitatively consistent with what is predicted from the depletion potential. In contrast, for motile bacteria, theoretical calculations suggest that swim force will exceed the depletion force at low polymer concentrations, resulting in no effective potential. Our experimental observations were consistent with these calculations; we saw no aggregation in motile bacteria at PEG 1 MDa Although the PEG induces the same depletion potential (pink arrows) at a given concentration, the swim force (white arrows) from the bacterial motility decreases both the well depth and the range, reducing their effective potential (blue dotted line) and preventing aggregation in the dilute PEG concentration regime. (f) The effective potential at contact (U min /kT) is plotted for motile and nonmotile E. coli against PEG concentration (mg mL À1 ). The vertical black dotted lines at c P = 0.4 and c P = 1.6 mg mL À1 denote the potential minima taken from the complete potentials plotted in (g) and (h), respectively. (g) The full effective potential (U eff /kT) is plotted against distance from the bacterial surface for both motile and nonmotile E. coli at c P = 0.4 mg mL À1 PEG. (h) The full effective potential (U eff /kT) is plotted against distance from the bacterial surface for both motile and nonmotile E. coli at c P = 1.6 mg mL À1 PEG.
concentrations less than 0.8 mg mL À1 . This effect is further illustrated by looking at the shape of the effective potentials for motile and nonmotile bacteria at c P = 0.4 mg mL À1 (Fig. 3e). For nonmotile bacteria, the total depletion range spans about 120 nm from the bacteria surface, and at inter-bacterial contact, the depletion well is at its minimum at approximately À5kT. In contrast, theory suggests that the swim force of motile bacteria will completely dominate over the depletion potential, resulting in no net attractive potential. The predicted lack of net attractions between motile bacteria is supported experimentally as we did not observe aggregation of motile bacteria at c P r 0.4 mg mL À1 ; the depletion forces are simply not strong enough to compete with the swimming force of the motile bacteria to induce aggregation.
At c p 4 0.4 mg mL À1 and as the polymer concentration enters the semi-dilute regime, effective potentials are insufficient to explain the experimentally observed aggregation results. At c P = 1.6 mg mL À1 (overlap concentration), where both motile and nonmotile E. coli aggregate, the depletion potential is strong enough to induce attractions between motile bacteria (Fig. 3f). Beyond the overlap concentration, in which the solution enters the semi-dilute regime, the magnitude of the effective potential minima continues to increase (Fig. 3d). Based on these predictions, we expect to see greater aggregation of nonmotile bacteria in the semi-dilute regime. However, in our experiments, we observed that aggregation actually began to decrease for the nonmotile bacteria at PEG concentrations of 40.8 mg mL À1 (Fig. 2a). Because the theory did not match with our experimental observations at these short time scales, it suggests that effective potentials are insufficient to explain the observed aggregation behavior at PEG concentrations greater than 0.8 mg mL À1 . We therefore proposed an alternate explanation for bacterial aggregation in the semi-dilute polymer concentration regime that relies on effective diffusivity. We test this hypothesis next.
Reduction in bacterial diffusivity explains aggregation in the semi-dilute polymer concentration regime where effective potentials are insufficient We hypothesized that in the semi-dilute concentration regime, the higher solution viscosity limits Brownian motion more appreciably than it limits bacterial motility. If this were the case, one would expect a reduced collision frequency among nonmotile bacteria leading to reduced aggregation at short time scales in viscous environments. In contrast, because motile bacteria are able to swim in viscous environments, 4 their rate of collision would remain high and they would aggregate at short time scales even in highly viscous environments.
To test our hypothesis, we compared the displacement of both motile and nonmotile E. coli in the dilute concentration regime (c P = 0.4 mg mL À1 ) and in the semi-dilute concentration regime (6.5 mg mL À1 ) by measuring their MSD (Fig. 4a). The MSD (in mm 2 ) is defined as MSD(t) = h|r(t 0 + t) À r(t 0 )| 2 i, where r (in mm) is a vector describing the position of the bacteria at a given time scale t (in s). At c P = 0.4 mg mL À1 , where nonmotile bacteria aggregate but motile do not, displacement of both motile and nonmotile cells is high. At this PEG concentration, we had already established that the depletion force is strong enough to aggregate nonmotile bacteria but is not strong enough to overcome the swim force induced by the motile cells. Although motile cells are likely undergoing more collisions due to their motility, there are insufficient attractions to initiate aggregation. The displacement of the nonmotile bacteria at c P = 6.5 mg mL À1 was substantially lower than their displacement at c P = 0.4 mg mL À1 , supporting the assertion that their collision frequency, and therefore aggregation rate, was decreased. The displacement of motile bacteria at this PEG concentration was much larger than the nonmotile bacteria, confirming that the higher viscosity did not impede movement and therefore collisions (Fig. 4b). The higher collision frequency of the motile bacteria, combined with the high depletion potential at the c P = 6.5 mg mL À1 PEG concentration, resulted in aggregation of the motile bacteria (Fig. 2). These results suggest that the at short time scales, aggregation becomes kinetically limited and is controlled by collision frequency of bacteria and their aggregates.

Conclusions
Bacterial aggregation often occurs in complex physiological environments such as the gastrointestinal tract or the lungs; in these settings it is important to understand the initial formation of aggregates before physiological factors such as peristalsis or swallowing play a role. Bacterial aggregation is thought to be the first step in colonization and biofilm formation; therefore, understanding the mechanism by which aggregates initially form could reveal important patterns in these processes. In this paper, we investigated how the interplay of swim forces that arise from motility and polymer-induced depletion forces influence bacterial aggregation. It has been reported previously that, at steady state, E. coli aggregate in the presence of non-adsorbing polymers, due to polymer-induced depletion forces. 18 At steady state, an effective potential incorporating contributions from swim forces and depletion forces adequately describes the aggregation behavior of the system, 18 wherein motility diminishes depletion-driven aggregation. Thus, polymers in greater quantities and/or sizes (i.e., a stronger depletion potential) are required for motile bacteria to aggregate compared with their nonmotile analogs. Here, we investigated the effect of bacterial motility on aggregation at short time scales (before the system reaches steady state).
We found that, in contrast to what has been demonstrated at steady state, at short time scales, motility can actually enhance aggregation, resulting in the formation of aggregates that are larger than those formed by nonmotile bacteria under the same conditions. For polymer-induced depletion forces to result in bacterial aggregation at short time scales, two conditions must be met: (i) sufficient inter-bacterial attractions, and (ii) sufficient inter-bacterial collisions. In the dilute polymer concentration regime, swimming competes with depletion to reduce the effective potential of a bacteria, resulting in less aggregation. Above the overlap concentration, motility allows the bacteria to overcome kinetic hindrances that may otherwise prevent aggregation at high viscosity environments, resulting in an increase in aggregation.
Future work should explore whether differences in motility impact aggregation of microbes and microbial physiology in more complex environments, e.g., in the gastrointestinal tract and the lungs. Beyond bacterial aggregation, this work may inform other studies of the behavior of active matter at early timescales.

Conflicts of interest
The technology described in this publication is the subject of a patent application filed by Caltech.