Spatial modulation of individual behaviors enables an ordered structure of diverse phenotypes during bacterial group migration

  1. Yang Bai
  2. Caiyun He
  3. Pan Chu
  4. Junjiajia Long
  5. Xuefei Li
  6. Xiongfei Fu  Is a corresponding author
  1. CAS Key Laboratory for Quantitative Engineering Biology, Guangdong Provincial Key Laboratory of Synthetic Genomics, Shenzhen Institute of Synthetic Biology, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, China
  2. University of Chinese Academy of Sciences, China
  3. Yale University, Department of Physics, United States
15 figures, 1 video, 3 tables and 1 additional file

Figures

Statistics of single-cell behavior during collective group migration.

(A) Bacterial population of diverse phenotypes sorted according to their chemotactic abilities (increasing from light green to dark green) during collective migration following the self-generated attractant gradient (brown color). Meanwhile, individual cells perform run-and-tumble random motions biased toward the migration direction. Scale bar reflects 0.1 mm. (B) Bacterial density profile ρz is stable (black solid line) in the moving coordinate z=x  VGt. Where the width, represented by the black dashed-dotted line, is defined by two times the standard deviation of the bacterial relative position (2σ, black dashed line). The instantaneous velocity (VIz) (blue solid line) is uniform and is equal to the average group velocity VG (blue dashed line). This observation was independent of the sampling time interval (Appendix 1—figure 2). (C) Sample runs of bacteria aligned by their initial positions (black dots) from the three regions (defined by colored regions in B) showed that cells at the back of the group tended to run forward, compared with cells in other regions. b, m, and f stand for the back, middle, and front of the migration group, respectively. Scale bar reflects 5 µm. (D) The exponential distributions of forward runs (solid lines) and backward runs (dashed lines) suggested the increasing efficiency of runs from the back to the front. (E) The mean run length in different directions confirmed the increasing efficiency of runs. (F) The expected drift velocity VD(z)=lR(z)cosθR(z) τR(z)+τT(z) (black solid line) decreased from the back to the front, and is linearly fitted by VD(z)-rz+VD0 , with r=0.05min1 and VD0=0.17mmmin1 . VD (z) crossing with the average group velocity VG (black dashed line), implied that bacteria perform mean reversion motions. The VD(z) profile was cut to present the majority of cells (~90%) (±1.65σ). (G) The time evolution of the average position (z, solid lines) of cells starting from the back, middle, and front of the migration group (color code was defined in B) confirmed the reversion behavior. Shaded area represents the s.e.m. of more than 450 cells (see Materials and methods). Analytically, the OU-type model predicted that z=C0e-rt-VD0-VG/r (dashed lines), where C0 can be fitted by the starting position (see supplementary text). In panels B, F, and G, the shaded area represents the s.e.m. of three biological replicates. The spatial bin size was 240 µm .

Ordered location of phenotypes of active particles in a moving gradient as a result of a pushed wave.

(A) Illustration of the Langevin-type model where the migrating bacteria are modeled by active particles under a deterministic decreasing driving force and a stochastic force. These forces allow the particles to form collective migration. On the moving coordinate of the migration group, active particles are bounded in a potential well generated by the driving force and the motion relative to the moving group (-VG). (B) The perceived gradient g(z) (blue line) is deduced from S(z), which is further calculated from the experimentally measured density profile (Appendix 1—figure 4D). This gradient profile was then transformed to a moving gradient gx,t=gz+VGt, and applied in the simulations of different active particles following Langevin dynamics. A linear fit of g(z) around z=0 is plotted (dashed blue line) and is applied to the Ornstein-Uhlenbeck (OU)-type model. (C) Particles migrate collectively with a moving gradient, while they are located in ordered locations according to their phenotype χi (green lines). Colors from light to dark represent increasing χ. (D) The effective potential wells generated by the effective force (χg(z)) are spatially ordered. Circles present the minimum of the potential wells. (E) The peak positions (cross) and the mean positions (circles) increase with χ. These were predicted by the OU-type model (Equation S15) (black line). The width of the density profile is defined by two times the standard deviation of all cells (2σ), which decreases with χ, which was consistent with the prediction of Equation S18 (blue line). (F–H) Group migration of mixed phenotypes under decreasing, invariant, and increasing gradient. Simulations were performed in one dimension (1D) with the chemotaxis ability of the group following a Gaussian distribution with a mean of 0.11 mm2 min−1 and a variation of 0.02 mm2 min−1 (Appendix 1—figure 5C). The force field used in the simulation followed the linearized force field, as in Equation 1, while g1>0 in F, g1=0 in G, and g1<0 in H. All force fields were set as moving with velocity VG . The dashed line in G represents a group with a single phenotype of χ=0.11mm2min1. Insert plots represent the wave front in the moving coordinate after 10 min of simulation. The width of the moving group, defined by the standard deviation of all particles, converges under decreasing gradient (F) while it diverges in cases of G and H.

Agent-based simulations recapture the ordered behavior of individuals.

(A) The expected drift velocity VDz of simulated bacteria decreased from the back to the front of the migration group, where the chemotactic ability χ ranged from 0.2 to 0.35 mm2min-1 , which was consistent with the experimental results shown in Figure 1F. The intersections between the VD curves with the preset group velocity VG (black dashed line) shifted toward the back of the migration group as χ decreased (circles). The different colors of the lines and circles correspond to different chemotactic abilities χ, as shown in the legend. The same color coding also applies to (B–D). (B) The reversion rate ri=|dVD,i(z)/dz| increased with chemotactic ability. (C) The effective potential well calculated by Uiz=z+Vd,i(z)dz. Positions of the potential minimum zmin are marked as circles. As illustrated, for a lower chemotaxis ability χ, the potential well is shallower and zmin shifts toward the back of the migration group. (D) The width of the density profile (measured by 2σ, see Figure 1B) decreases with the reversion rate ri as well as the chemotaxis ability χi . The mean square displacement (MSD) of bacteria (insert, solid lines) is bound to 2σi2 (insert, dashed lines) (see supplementary text). In panel (A, C), curves were cut to present the majority of cells (90%) (zmin±1.65σi). More details on the results of this simulation are presented in Appendix 1—figure 7.

Spatial ordered structure emerged from the behavioral modulation of cells with different chemoreceptors.

(A) Genetic circuit of the Tar-titratable strain. In the experiments, the expression level of Tar (a chemoattractant receptor protein) was titrated by the concentration of external inducer (aTc). The chemotactic ability χ of bacteria was then determined by the expression level of Tar (Adler, 1969). (B) The expected drift velocity VD,iz of the Tar-titratable strain JCY20 (colored solid line) was spatially modulated and decreased from the back to the front of the migration group and intersected with a group migration velocity VG0.15mm/min (black dashed line). The linear fit of VD,iz (colored dashed lines) intersected with VG at positions (circles) determined by the corresponding Tar expression level. Colors from dark to light green corresponded to inducer (aTc) concentrations of [1, 3, 6, 20]ng/mL. The black shaded area of VG represents the s.d. of four experiments, while the colored shaded area of the VD curves presents the s.e.m. of the counted runs. (C) In the experiment, the positions of the VD,i(z)VG intersections (circles, illustrated in B), together with the peaks (stars, illustrated in the insert figure), and the average positions (diamond) of the bacterial density profiles, all shifted toward the front of the migration group for strains with higher Tar expression levels, which had higher chemotactic abilities and migrated faster on agar plates (x-axis, see Materials and methods; Appendix 1—figure 9). The related density profiles (PDF) were shown in the insert plot and the color coding of lines/symbols in both panels C and D was the same as that in B. (D) The width of density profiles (2σ) of Tar-titrated bacteria decreased with the reversion rate r.

Appendix 1—figure 1
Experiment protocol, sample tracks, and fluorescent labeled cells.

(A) The microfluidic chip used in this study is consisted of a long channel (20 mm×0.6 mm×0.02 mm) and a large chamber (3 mm× 3 mm). Prepared bacteria are first loaded in the chip through inlets with M9 motility buffer (see Materials and methods). Next, the chip is spun at 3000 rpm for 15 min so that the bacteria are concentrated to one end of the channel. Bacteria consume aspartate (Asp) as the single chemoattractant to form a moving gradient that attracts the group migration. Finally, the bacterial motions are tracked under microscope of 4× objective by frame rate of 9 fps. (B) Cells expressing fluorescent protein (strain JCY2) are premixed into strain RP437 by 1:400 and counted each frame. Cell counts in nine successive frames (1s) were plotted within bin size of Δx=100μm. Colors from light to dark red represent time from 1 to 10 min. (C) Sample tracks of bacteria (colored lines) and bacterial signal (black dots) on the last captured frame. (D) Fluorescent labeled cells (strain JCY1 and JCY2) and the wild-type strain (RP437) were tracked separately on growth medium. All three strains show the same distribution of tumble bias (left), mean speed (middle), and mean run speed (right) of each track. All these quantities are average of more than 9000 tracks, weighted by their trajectory time. Each curve show data of one experiment.

Appendix 1—figure 2
Determination of coherent migration group.

(A) Positions of first n-th cell (counted from the front, from blue to red n=1.24 2.48 3.72 4.96 6.2 7.44 8.68 9.92 1.16×104) increase linearly with time, of which the slope gives the traveling speeds. (B) The traveling speeds of the first n-th cells are close to the group velocity VG (black dashed line), representing a coherent migration group. (C). The instantaneous group velocity VGt was constant over time. VGt was averaged for all tracks over 30 s (blue dots). Black dashed line represents the average group velocity for all tracks on all times VG . (D). The constant instantaneous velocity VIz=Δx(z)/Δt for sampling time Δt=1s (blue solid line) and Δt=10s (red solid line), where shaded area was s.e.m. of three replicates. Both curves equal to the group velocity VG , indicating the migration group is coherent. Data in (A–C) was extracted from one experiment.

Appendix 1—figure 3
Modulation of bacterial behaviors.

(A) Distribution of forward run duration (solid line) and backward run duration (dash line) in three regions defined in Figure 1D and Appendix 1—figure 3E. (B) Mean run duration in different directions of bacteria in three regions. s.e.m. of four replicates is smaller than thickness of lines. Angular bin size is 15°. (C) Spatial structure of adapted tumbled bias (blue) and run speed (red). (D) The mean run length in different directions, with the angular bin size of 15°, also shows that cells in the back were better skewed to run forward. (E) The run length bias Bl(z)=lR(z)cosθR(z) lR(z) (black solid line) and the run time bias Bτ(z)=τR(z)cosθR (z)τR(z) (blue solid line) both decreased from the back to the front of the migration group. In panels C–E, shaded area represents s.e.m. of three biological replicates. The spatial bin size is 240 µm. (F) Representative examples of single-cell trajectories (three colors represent three different tracks) showed the reversion behavior of bacteria around their mean positions. (G) The mean reversion behavior results in limited diffusion of bacteria in the moving frame. So that the mean square displacement (MSD) in the direction of migration (black curve) was lower than the pure diffusion of bacteria measured on free liquid with M9 motility buffer (blue curve). The shaded area on black curve represents the s.e.m. of more than 2900 tracks, while the shaded area on blue curve represents the s.e.m. of more than 1900 tracks.

Appendix 1—figure 4
Behavior analysis and deduction of the chemoattractant profile.

(A) All cell tracks are divided into three sectors: forward (|θ|45) (green line), sideward (45<|θ|<135) (blue line), backward (|θ|135) (red line). The run durations in all three sectors increase form back to front. (B) Spatial distribution of reorientation angle. (C) Reorientation angles defined in Saragosti et al., 2011 are plotted over the direction of previous runs. Such angles are biased to the back, indicating that bacteria running down the gradient reorient more in the next tumble. Shaded area shows s.e.m. of 14,892 tracks. The s.e.m. in panel C is smaller than the thickness of the curve. (D) The chemoattractant profile S(z) used for the simulation is interpolated by Δz = 0.001mm (blue line) from the experimental S(z) (black circles). While the experimental chemoattractant is deduced from the averaged density profile ρz (Figure 1B).

Appendix 1—figure 5
Particle-based simulation with chemoattractants consumption.

(A) Simulated density profiles of active particles of single phenotype of χ=0.11mm2/min . VG=3.0μm/s over time. The peak positions increase linearly with time (insert plot). (B) Simulated density profiles of active particles of heterogenous phenotypes distributed in (C). The peak positions follow a stable speed of VG=2.9μm/s (insert plot). (D) Time-shifted density profiles of different phenotypes exhibit spatially ordered in moving coordinate. Cells in the gray area drop the group during migration. Time-shifted chemoattractant profile increases from back to front (E) while the perceived gradient decreases with a quasi-linear region around 0. (F) Shaded area in D, E, F represents s.d. of 40 profiles recorded at stable region (from 21 to 60 min) with time step of 1 min.

Appendix 1—figure 6
Agent-based simulations with consumption of chemoattractant.

(A) Simulated density profiles of bacteria integrated with chemotactic pathway show stable migration after 20 min simulation. Tracks were obtained in from 30 to 32 min, and were calibrated into a moving coordinate through group migration speed. (B) Four subgroups of different receptor gain N were orderly distributed. (C) The expected drift velocity VDz is decreasing and is spatially ordered.

Appendix 1—figure 7
Agent-based simulations of cells with different receptor gains.

(A) Bacteria of different chemotatic ability (χ) form ordered structure in a moving chemoattractant field as in Appendix 1—figure 4D with fixed velocity VG=3.0μm/s. (B) The average positions of each phenotype (diamonds), the peak positions (crosses), and the cross positions between their VD curves and the group velocity VG (circles) increase with χ. These positions follow the same trend as predicted by Ornstein-Uhlenbeck (OU) model (black line). (C, D) Run duration bias Bτ and run length bias Bl of different subgroups decrease from back to front. Both curves were also spatially ordered by χ (color from light to dark green).

Appendix 1—figure 8
Agent-based simulations of cells with different CheY-P levels and adaptation times.

Simulations results for CheY-P concentration Yp0 ranging from 2.70 to 2.85 µM (A–D) and for adaptation time τM ranging from 1 to 10 s (E,H) are presented. In both cases the cell density ρz (A, E) and the expected drift velocity VDz (B, F) are spatially ordered according to the varying parameter. The mean position (crossed) together with the cross position between VD(z) and VG (circles) increase with Yp0 (C) while they show non-monotonic trend to τM (G). The width of the density profile (defined by s.d. of all cells) and the reversion rate decrease with Yp0 (D). And they showed a non-monotonic dependence on τM . The receptor gain N is set to 24 in both cases, while Yp0 is set to 2.77 µM for the τM varying case and τM is 1 s for the Yp0 varying case. Simulations are performed under a steeper moving chemoattractant gradient of g1=2.5mm2 .

Appendix 1—figure 9
Migration rates, growth rates, motilities and coherent migration of the Tar-titratable strain.

(A) Tar-titratable strain (JCY20) allows us to change its chemotaxis ability (measured by its migration rate on agar plate) by adding inducer ([aTc]) of different concentrations. (B) This strain has constant growth rates measured in liquid batch culture. (C–D) Motility statistics were extracted from bacterial tracks in free liquid. The growth rate, tumble bias, and mean run speed were found constant with all tested aTc concentrations. Error bar represents s.d. of three biological replicates. (E) The density profile of the wild-type (WT) strain with red fluorescent protein (RFP) (strain JCY2) (red lines) and Tar-titratable strain with yellow fluorescent protein (YFP) (strain JCY20) (yellow lines) are plotted over space at time before (dashed lines) and after tracking (solid lines). Because the density profiles of both strains before and after tracking are considered to be identical, the migration group is stable during the tracking period (10 min ). The fluorescent strains (JCY2 and JCY20) are mixed with WT non-fluorescent strain (RP437) by ratio of 8:1:400. (F) During the 10 min tracking period, only the Tar-titratable strain (JCY20) is recorded. For each half a minute, we plotted the density profiles and find their peak positions (yellow circles). Such peak positions increase linearly with time (fitted by the black line). Data were extracted from one experiment. (G) The density profile of the Tar-only strain with YFP (strain UU1624 from Prof. Parkinson) are plotted over space at time before (dashed lines) and after 10 min (solid lines).

Appendix 1—figure 10
Particle based simulation of the signal secretion model.

(A) Simulated density profiles of heterogenous phenotypes distributed in C show (B) stable group migration with its peak positions following a constant speed of Vp=6.69μm/s. (C) The phenotypic diversity was set by a Gaussian distribution on χ, with mean value of 1 and s.d. of 0.3. (D) Time-shifted density profiles of different phenotypes exhibit spatially ordered in moving coordinate. Color code follows C, while cells in the gray area drop from the group. (E) Time-shifted signal profile Ss , which plays the role of ‘global driving force’, decreases from back to front.

Author response image 1

Videos

Video 1
Bacterial tracking in the migration group.

Each dot represents a bacterium captured under microscope. Lines represent typical tracks. The colors of tracks from light green to dark green represent the mean positions of the tracks from the back to the front of the migrating group.

Tables

Table 1
Summary of quantities.
QuantitiesDefinitionFormulations
VGGroup velocityVG=dxdt
zMoving coordinatez=x-VGt
VIzInstantaneous velocityVIz=dxzdt
VDzExpected drift velocityVD(z)=LR(z)cosθR(z) τR(z)+τT(z)
ρzCell densityρz=izaΔz
SzChemoattractant concentration-VGdSdz=Ds2Sz2-kρ
gzPerceived gradientg(z)=dln (1+S(z)/Koff1+S(z)/Kon) dz
Appendix 1—table 1
Strains and plasmids used in this study.
StrainsDescriptionSource
RP437Wild-type strain for chemotaxisFu et al., 2018
JCY1RP437+ a Cmr plasmid with YFPThis study
JCY2RP437+ a Kanr plasmid with mRFP1This study
JCY20Construction based on RP437, tar<> loxp, bla:Ptet-tetR-tar at attB site,+ a Cmr plasmid with YFPThis study
UU1624Δaer-1 Δtsr-7028 Δtap-3654 Δtrg-100Gosink et al., 2006
Appendix 1—table 2
Summary of parameters in agent-based simulation.
SymbolDefinitionValueSource
χChemotaxis abilityDefined in contextThis study
kConsumption rate of chemoattractantFitted parameter for migration speedThis study
KonUpper limit of chemoattractant sensing1 mMFu et al., 2018
KoffLower limit of chemoattractant sensing1 µMSi et al., 2012
DsDiffusion coefficient of chemoattractant1000 µm2/sFu et al., 2018
v0Run speed of bacteria30 µm/sKeller and Segel, 1971b
τMAdaptation time1 sWaite et al., 2016
NReceptor gainDeduced from chemotactic abilitySI notes in Section 7
MMotor number5Sneddon et al., 2012
Yp0Basic CheY-P level2.77 µMCorrespond toTB00.27

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Yang Bai
  2. Caiyun He
  3. Pan Chu
  4. Junjiajia Long
  5. Xuefei Li
  6. Xiongfei Fu
(2021)
Spatial modulation of individual behaviors enables an ordered structure of diverse phenotypes during bacterial group migration
eLife 10:e67316.
https://doi.org/10.7554/eLife.67316