Flagellar number governs bacterial spreading and transport eiciency Javad Najai, 1, M. Reza Shaebani, 2, Thomas John, 1 Florian ltegoer, 3 Gert Bange, 3 and Christian Wagner 1,4, 1 Department o Experimental Physics, Saarland University, 66123 Saarbrücken, Germany. 2 Department o Theoretical Physics, Saarland University, 66123 Saarbrücken, Germany. 3 Department o Chemistry and LOEWE Center or Synthetic Microbiology, Philipps University Marburg, 35043 Marburg, Germany. 4 Physics and Materials Science Research Unit, University o Luxembourg, 1511 Luxembourg, Luxembourg. arxiv:1901.05373v1 [physics.bio-ph] 16 Jan 2019 Peritrichous bacteria synchronize and bundle their lagella to actively swim while disruption o the bundle leads to tumbling. It is still not known whether the number o lagella represents an evolutionary adaptation towards optimizing bacterial navigation. Here, we study the swimming dynamics o dierentially lagellated Bacillus subtilis strains in a quasi-two-dimensional system. We ind that decreasing the number o lagella N reduces the average turning angle between two successive run phases and enhances the duration and directional persistence o the run phase. s a result, having less lagella is beneicial or long-distance transport and ast spreading, while having a lot o lagella is advantageous or the processes which require localization and slow dynamics, such as bioilm ormation. We develop a two-state random walk model that incorporates spontaneous switchings between the states and yields exact analytical expressions or transport properties, in remarkable agreement with experiments. The results o numerical simulations based on our two-state model suggest that the eiciency o searching and exploring the environment is optimized at intermediate values o N. The optimal choice o N, or which the search time is minimized, decreases with increasing the size o the environment in which the bacteria swim. INTRODUCTION Many bacterial species swim by the rotation o lagella [1]. Several lagellation patterns can be distinguished according to the lagellar arrangement on the cell body, ranging rom polar (one lagellum at one pole) to peritrichous (helical arrangement on the whole cell body) [2]. Each lagellum is anchored within the cell membrane to a reversible rotary motor that switches between clockwise and counter-clockwise rotation[3]. The ability o the motor to turn in either direction allows polarly-lagellated bacteria, such as Caulobacter crescentus, to change the swimming direction o the cell [4]. peritrichouslylagellated bacterium experiences an alternating sequence o active runs and tumbles [5], controlled by the rotational states o each lagellum. When all lagella rotate counter-clockwise, they orm a bundle and synchronize their rotation. s a result, the bacterium moves orward smoothly (run phase). However, the bundle is disrupted when some o the lagella switch their rotational direction, which acilitates changes in the direction o motion in relatively short time intervals (tumble phase). The run-and-tumble dynamics o bacteria allows them to change their direction o motion. They can also adjust the duration o stay in the run phase in response to environmental changes induced, e.g., by temperature, light, or chemical gradients [1, 5]. Such an ability is known to be highly advantageous to enhance search eiciency [6 8], allowing or an optimal navigation towards avorable or away rom harmul regions [9]. J.N. and M.R.S. contributed equally to this work. Corresponding author. Email: c.wagner@mx.uni-saarland.de Despite the advances in understanding o the underlying mechanisms o bundle ormation and disruption [10 17] and its relationship to the chemotaxis signaling network [1, 9, 18, 19], it remains poorly understood how the number o lagella inluences the undamental properties o bacterial propulsion and transport. ttempts to understand the role o lagellation have mainly ocused on Escherichia coli which has up to eight lagella, while the lagellar number o many types o bacteria is typically much higher. So ar, it has been observedthat the torque and thus the swimming speed [11] and the raction o time spent in run or tumble phase [12] o peritrichouslylagellated bacteria remain roughly independent o the lagellar number N. Recent numerical studies also predict that the swimming speed is only slightly aected by N, assuming that just a single bundle is ormed [20, 21]. Here, we provide evidences or the ormation o multiple bundles, revealing that the ormer assumption does not hold at least at high lagellar number. key question rom an evolutionary point o view is based on which criteria a bacterium chooses to have a certain number o lagella. In the present study, we investigate bacterial motility over a wide range o the lagellar number and clariy how the run-and-tumble dynamics is inluenced by the choice o N. We study Bacillus subtilis, a rod-shaped bacterium commonly ound, e.g., in soil and the gastrointestinal tract o humans. The number o lagella is regulated by the master regulator in Bacillus subtilis [22]. Deletion o the gene leads to a reduction rom 26 to 10 lagella, while cells overexpressing exhibit up to 40 lagella. Thereore, we employed a wild-type Bacillus subtilis (NCIB3610 strain) with 26 lagella, a strain
2 deicient or (9 lagella), and a one which carries under the control o an IPTG-inducible promoter (41 lagella) (Fig.1). While the speeds in both run and tumble phases and the mean duration o tumbling periods show no systematic dependence on N, the duration and directional persistence o the run phase and the average turning angle between two successive run phases vary monotonically with N. Smaller values o N are more avorable or long-distance transport and ast spreading because o more persistent trajectories and longer excursion times in the run phase. On the other hand, increasing N acilitates tumbling, leading to slow dynamics and localization which is beneicial e.g. or the ormation o bioilms. We develop a coarsegrained analytical approach to study the run-and-tumble dynamics o bacteria, which enables us to identiy the contributions o the inluential actors and map out the phase diagrams in the space o the relevant parameters o run-and-tumble statistics. Our simulation results suggest that by adopting intermediate values o N (as e.g. in the wild-type strain), the search eiciency and the ability to explore the environment can be optimized. 1.00 RESULTS We report and compare the dynamics o three strains o Bacillus subtilis with dierent lagellar number N obtained by genetic manipulation (see Materials and Methods section or details). The strains, wt (wild-type NCIB3610) and have 9±2, 26±6 and 41±6 lagella, respectively [23] (Fig. 1). In contrast to the studies on tumbling o trapped cells [12], here we track the motion o reely swimming cells. ew typical trajectories o each strain are shown in Fig. 1B, which consist o successive active run and tumble phases. Distinct motility patterns Strikingly, we ind that the curvatures o the trajectories strongly depend on the lagellar number. The paths are clearly less curved with decreasing N, which may originate either rom the change in their tumbling statistics or rom having a larger persistence length in the run phase. While the inluence o N on the run-and-tumble statistics will be discussed in the next section, here we ocus on the latter possibility, i.e. the variations o the persistency in the run phase with changing N. To characterize the running persistency, we note that each run trajectory comprises a set o recorded positions o the bacterium. Given these data, one can extract the local direction o motion α and, thus, the directional change θ along the trajectory, resulting in the distribution o directional changes P(θ) [24]. Then, the persistency p o the run phase can be obtained as the Fourier transorm o P(θ), i.e. p = π dθ eiθ P(θ)= cosθ [25, 26]. lternatively, the curvature o the trajectories can be characterized by calculating the rotational mean square displacement RMSD(τ)= ( α(t+τ) α(t) ) 2 =2Dr τ, rom which one can extract the asymptotic rotational-diusion coe- 0.98 0.96 0.94 0.92 FIG. 1. Bacterial motility patterns. () Illustration o three strains o Bacillus subtilis with dierent number o lagella. (B),(C) Planar projections o quasi two-dimensional experimental (B) and simulation (C) paths. ew trajectories are randomly chosen and translated so that the irst point is located at the origin. The simulation parameter values are extracted rom experimental data. (D) Directional persistency o the run phase p (circles), and asymptotic rotationaldiusion coeicient D r (squares) or dierent strains. icient D r. The results shown in Fig.1D reveal that the directional persistency in the run phase remarkably depends on the lagellar number; the run paths are less curved or smaller N. Run-and-tumble statistics Next, we turn to studying the statistics o run and tumble events. detailed description o the tumbling detection procedure can be ound in Materials and Methods section. ccording to the results shown in Fig. 2, the tumbling speed is smaller (but non-negligible) compared to the running speed [27]. Moreover, the speeds in both run and tumble phases show no systematic dependence on N. closer look at the probability distributions o the instantaneous speeds (see e.g. the running speed data in Fig.2B) reveals that even the distributions o a given strain slightly dier rom one culture to another due to biological heterogeneity. Thus, there is no signiicant di-
3 Running Tumbling C D Running Tumbling F (x ) 10 8 6 Speed 4 2 B (x ) 10 8 6 E G 0-180 -90 0 90 180 4 2 0 0 10 20 30 40 Running Speed 50 60 0 8 16 24 32 40 FIG. 2. Run-and-tumble statistics. () The mean run and tumble speeds or dierent strains (open symbols). The ull symbols indicate the results or dierent cultures o each strain. The horizontal lines show the mean values over all strains. (B) The probability distribution o the run speed. The thin solid curves (bright red) indicate the distributions or dierent cultures o wt strain. (C) Fluorescence images o the bundles o a wild-type strain. Scale bar, 3.5µm. (D) The average duration o run (circles) and tumble (diamonds) periods. The horizontal lines indicate the average values over all strains. (E) The probability distribution o observing a running time longer than a given duration t. (F) The turning-angle probability distribution R(φ). (G) The average turning angle φ (triangles) and R= cos φ (squares) or dierent strains. erence in the speeds o dierent strains. By approximating the bundle as a single helix [4, 20, 28], simple models based on resistive-orce theory predict a weak (logarithmic) growth o running speed with N, though experiments have shown that hydrodynamic dissipation along the cell body considerably weakens the eect [11]. dditionally, the assumption o a single bundle does not hold in general, according to our observations o luorescentlystained lagella o swimming bacteria as well as other reports in the literature[29, 30]; even multiple bundles may orm during the active running phase(several examples o a wt strain with multiple bundles are shown in Fig.2C). The diversity o the locations and interactions between the bundles makes the prediction o the propulsive orce and swimming speed highly complex. To understand whether and to what extent the lagellar number inluences the switching requencies between run and tumble phases, we plot the average excursion times in Fig.2D. It turns out that the mean tumbling time ( 0.2s)andthustheabilityorestoringtherunphaseare not inluenced by N. Given the mean tumbling time and speed, one inds that the displacement in the tumbling phase is o the order o the cell length (i.e. a ew micrometers). However, such a swimmer with the running speed o 30µm/s in an aqueous medium is expected to stop ater only 0.6µs and traveling approximately 0.1 i the pushing orce suddenly stops [31]. This suggests that the bundles are not ully disrupted in the tumbling phase, thus, the bacterium is partially propelled. The mean running times remarkably depend on the lagellar number. Increasing N enhances the probability o switching rom running to tumbling. This supports the veto model or bundle disruption [13, 19], which proposes that clockwise rotation o a single or a ew lagella (less than N 2 ) is suicient to disrupt the bundle. Figure 2E shows that the tail o the probability distribution o running times is nearlyexponential, with a N -dependent slope (data censoring or the irst and last incomplete periods sharpens the trend even more). We observe similar exponential decays or the tumbling time distributions (not shown), veriying that the switching events between the two states o motility happen spontaneously and can be described by Markovian processes. nother characteristic o the run-and-tumble dynamicsthat isaected bythe numberolagellais theturning angle between successive run phases. We deine the turning angle φ between two successive runs as the change in the direction o motion rom the end o one run to the beginning o the next run (the direction o motion was obtained by a linear it to our data points or each
4 run phase). The turning-angle distribution R(φ) obeys let-right symmetry and develops a peak or all strains. However, the peak and, more clearly, the mean o R(φ) shit towards larger angles with increasing N, as can be seen in Figs.2F and 2G. Similar to the curvature o run trajectories, we can use the Fourier transorm o R(φ), i.e. R = π dφ eiφ R(φ)= cosφ, as a measure o the directional change between consecutive run phases. R ranges between 1 and 1, with 1, 0, and 1 denoting 0, 90, and 180 turning, respectively. Indeed, R is correlated with the tumbling time and speed [32] as well as the strength o the torque exerted on the cell body during the re-ormation o the bundle [11]. Since the speed and excursion time in the tumbling phase are not considerably aected by N, we attribute the reduction o R with increasingn to experiencing a largertorque during bundle re-ormation at higher lagellar numbers [11]. We also checked that the direction o motion does not vary signiicantly when switching rom running to tumbling. This act and the relatively high speed in the tumbling phase once again indicate that the bundles are not ully disrupted in the tumbling phase. s a result, one expects that the tumbling is not a purely random motion and should exhibit directional persistency. To veriy this idea we assign a dimensionless number l tumble to each tumbling segment o the trajectory, obtained by dividing the end-to-end distance by the total length o the trajectory segment. The probability distribution P( l tumble ) shown in Fig.3develops a sharp peak at l tumble 0.9, indicating that most o the tumbling trajectories are nearly straight during the relatively short periods o tumbling. Finally, we measure the mean square displacement MSD(τ)= ( r(t+τ) r(t) ) 2 or dierent strains. ccording to our indings, the strains with lower lagellar number have less-curved run trajectories, switch less requently rom run to tumble phase, and experience a smaller turning angle between successive run phases. Consequently, the crossover rom a persistent motion to the asymptotic diusive dynamics is expected to occur at longer times or smaller values o N, which is conirmed by the experimental results shown in Fig.3B. The asymptotic diusion coeicient D is also a decreasing unction o N (see Fig.3C). These results indicate that or increasing the eiciency o long-distance transport, it is more beneicial to have less lagella. The strains with more lagella spend a larger raction o time in tumbling state and have a shorter persistence length when swimming. We note that the trends (as a unction o N ) discussed in this section are robust against variations o tumbling detection thresholds. Model Based on our indings, we develop a stochastic coarsegrained model or a random walk with spontaneous switchings between two states o motility. We irst validate the approach by comparing the analytical predictions with experimental results. Then the model enables us to better understand the inluence o the lagellar num- B C (ballistic) (di usive) FIG. 3. Spreading and transport properties. () Probability distribution o the dimensionless quantity l tumble, relecting the bending o tumbling trajectories. (B) Evolution o the mean square displacement or dierent strains, obtained rom experiments (thin dark-color lines), simulations (symbols), and the model via Eq. 13 (thick bright-color lines). (C) The asymptotic diusion coeicient D or dierent strains. The analytical predictions (circles) are compared to the experimental results (squares). ber on bacterial motility, by numerically exploring the phase space o the run-and-tumble statistical parameters and identiying the isolated role o each inluential actor. Moreover, our simulations based on the two-state model evidence or a nonmonotonic dependence o the mean irst-passage time on N. Stochastic two-state models have been widely employed to describe altering phases o motion o e.g. swimmers and cytoskeletal motor proteins, or the locomotive patterns in other systems [33 39]. iming at obtaining analytical insight by including only the most prominent characteristics o the bacterial motility pattern into the model, we introduce a random walk approach in which the walker experiences two states o motility: (i) running with the mean speed v R and persistency p, and (ii) tumbling with the mean speed v T. The switching probabilities rom run to tumble state and vice versaare supposed to be asymmetric and denoted by and, respectively. ssuming constant switching probabilities results in exponential distributions or the residence times in each state, which is in agreement with our experimental observations. In the run state, a persistent random walk in the continuous space is considered which is char-
5 C E 2.0 opt run tumble recorded positions MFPT 1.5 1.0 0.5 B D F run-phase persistency 1.0 eicient long-distance transport / ast spreading optimal global search and cover time slow dynamics / localization increasing lagellar number FIG. 4. Phase diagram o the inluential parameters. () typical sample trajectory o Bacillus subtilis with run-andtumble dynamics. (B)-(D) Three cross-sections corresponding to (B) R= cosφ =0.56, (C) =0.05, and (D) p= cosθ =0.97 o the three dimensional phase diagram in the (R, p, ) space. The color intensity relects the magnitude o the asymptotic diusion coeicient D. The marked regions indicate the accessible range o parameters in our experiments with Bacillus subtilis. (E) Mean irst-passage time MFPT, scaled by MFPT at p=0 versus the directional persistency p o the run phase. (inset) Optimal persistency p opt vs the eective system size L/l. (F) Schematic picture depicting how dierent aspects o transport eiciency vary with N. acterized by its speed v R and the distribution P(θ) o directional changes along the run trajectory. Thus, the directional persistence p= π dθ eiθ P(θ)= cosθ quantiies the curvature o the run trajectories. Based on the results shown in Fig. 3, the tumbling trajectories are approximated by straight lines along the last direction o motion in the run phase (see Fig.4). The last key property is the turning angle φ between two successive run phases. Using the turning-angle distribution R(φ), we calculate R= π dφ eiφ R(φ)= cosφ to quantiy the directional change between consecutive run phases. Denoting the time spacing between consecutive rames by t, we describe the process in discrete time by introducing the probability densities P R t (x,y γ) and P T t (x,y γ) to ind the walker at position (x,y) along the direction γ at time t in the run and tumble states, respectively. The dynamical evolution is deined by the ollowing set o master equations P R π t+ t (x,y γ) = (1 ) dβp(γ β)pt R (x v tcosγ,y v tsinγ β) R R π + dβr(γ β)pt T (x v tcosγ,y v tsinγ β), R R P T t+ t (x,y γ) = (1 )P T t (x v T tcosγ,y v T tsinγ γ) + P R t (x v T tcosγ,y v T tsinγ γ). (1) By the two terms on the right-hand side o each equation, we consider the possibilities o being at each o the two states in the previous step. On the right-hand side o the second equation, the possibilities o moving along a straight line or switching rom running to tumbling are considered. For simplicity we assume that the direction o motion does not change when switching rom run to tumble state. By means o a Fourier-z-transorm technique [40, 41], we obtain exact analytical expressions or temporal evolution o arbitrary moments o displacement (see Materials and Methods section or the details o the method and the explicit lengthy expression o MSD). The parameter values extracted rom experiments are given in Table1 (or the calculation o p we chose t= 1 6 s). Using these data as the input o the model, the time evolution o the mean square displacement can be predicted via Eq.13 (plotted in Fig.3B or dierent strains). We also perormed extensive Monte Carlo simulations o a stochastic process with the same parameter values. The results agree perectly with the analytical predictions. The agreement with experiments is also strikingly good; Despite all the simpliications, the model recovers the MSDs o all strains over the whole Table 1: Parameter values extracted rom experiments. strain v R ( µm ) v s T (µm) t (s) t (s) p R s R T 24.5 12.1 4.39 0.244 0.99 0.65 wt 29.8 14.0 2.27 0.224 0.98 0.59 23.3 10.3 1.18 0.240 0.95 0.43
6 range o time. By introducing 2 v 2 v2 and T R B 2 2 Rv 2+22 T v 2, we obtain the ollowing expression or the asymptotic diusion coeicient R [ D = 1 t + B+2 (1+R)v ] Rv ] T. 4 + [1 p (R p) (2) The measured diusion coeicients or dierent strains are in satisactory agreement with the analytical predictions, as shown in Fig.3C. Equation (2) reduces to the asymptotic diusion coeicient o a single-state persistent random walker or =0 and =1 [42]. While it is extremely diicult to systematically vary the lagellar number in experiments, we can numerically explore the ull phase space o the inluential parameters. mong the key elements in determining the transport properties o the Bacillus subtilis, the speeds v R and v T, and the switching probability do not vary signiicantly with N. Thus, we ix them at their average experimental values, which allows us to reduce the degrees o reedom to three sensitive parameters, i.e. R, the directional persistency p, and the switching probability. Figures 4B-4D show 2D proiles o D in the (R, p, ) phase space, revealing that D varies by several orders o magnitude by varying these key parameters. However, limiting the parameter values to the accessible range in experiments with Bacillus subtilis reveals that here D and the long-distance transport are mainly aected by the variation o the run-phase persistency p. For example, by ixing the two other degrees o reedom at their average experimental values and calculating the variations o the asymptotic diusion coeicient when the third parameter varies within the accessible range in our experiments, we get Dmax D min =1.2,1.8, and 3.6 or varying R,, and p, respectively. More generally, by changing the lagellar number o other types o bacteria one may deal with other regions o the 6-old phase space. Our analytical approach allows us to extract urther inormation about the transport properties, such as the crossovertimet c toasymptoticdiusiveregimeasaunction o the key parameters. The crossover time can be estimated by balancing the linear terms in time and the nonlinear contribution in the MSD equation. We ind that t c also varies by several orders o magnitude in the (R, p, ) phase space (not shown). Note that due to the interplay between the key parameters, multiple transitions between dierent types o anomalous diusive dynamics occur on varying time scales in general [35]. Finally, we study the mean irst-passage time (MFPT) o a two-state persistent random walk in coninement by means o extensive Monte Carlo simulations. The MFPT is deined as the mean time taken by the random walker to reach a particular position in the system or the irst time. For simplicity, random walk on a cubic lattice with periodic boundary conditions is considered, and all run-and-tumble parameters except p (i.e. the most sensitive actor upon varying N ) are ixed at their mean experimental values. The walker is supposed to search or hidden targets during its both states o motility. The results reveal that the MFPT is a nonmonotonic unction o p and admits a minimum (Fig.4E). From the monotonic dependence o p on N (Fig. 1D), one concludes that the search eiciency is optimized at intermediate values o N. Denoting the size o the coninement with L, and the reaction range o the searcher with l, the eicient persistency grows with increasing L/l (inset o Fig. 4E); thus, the optimal number o lagella N opt shits towards smaller values with increasing the eective system size. DISCUSSION Our investigation o the dynamics o Bacillus subtilis revealed that the strains with lower lagellar number N have higher persistency and longer excursion time in the run phase and, additionally, their abrupt directional change when switching back to the run phase is smaller; Consequently, having too ew lagella enhances the eiciency o long-distance transport (relected e.g. in their higher translational diusion coeicient), and causes ast spreading. Increasing the number o lagella considerably increases the probability o switching rom running to tumbling as well as the curvature o the trajectories in the run phase; Moreover, switching rom tumbling to running is accompanied by a larger directional change when the bacterium has more lagella. s a result, the overall orientation changes more requently i there are too many lagella. This leads to slow dynamics and can be beneicial or the processes which require localization, such as the ormation o bioilms. However, the lagellar number o the wild-type Bacillus subtilis alls in the middle o the range, where the search eiciency could be optimized according to our simulation results (see Fig. 4F). Intermittent random walks are known to be beneicial or search eiciency in conined geometries [43 45]. For example, when slow diusion periods are interrupted by ast relocating ballistic lights, it has been shown that there exists an optimal ratio between the excursion times in the two states which leads to a global minimum o the search time [43]. The optimal choice varies with the system size and the reaction range o the searcher. In addition, it has been recently shown or a single-state persistent random walk, that the mean irst-passage time to ind a target admits a minimum as a unction o the persistency [46]. In a general two-state random walk with tumbling and persistent running phases (such as the motion o Bacillus subtilis), one expects that urther complexity to determine the optimal search time arises due to the interplay between the directional persistency p o the run phase, the switching probabilities and between run and tumble (which determine the ratio between the excursion times in the two states), the size o the coninement, and the characteristics o the bacterial chemotaxis system (esp. the reaction range). However, rom our experimental indings we expect
7 that changing the lagellar number o Bacillus subtilis inluences the dynamics and, thus, the search time mainly via aecting the run-phase persistency p. I the number o lagella o wild-type Bacillus subtilis has been solely evolved towards optimizing the search eiciency, then our simulation results suggest that it should be most appropriate to live in environments with a typical size o the order o a ew hundreds larger than the reaction radius o the bacterium (e.g. an environment size o less than a millimeter or a reaction radius comparable to the typical size o Bacillus subtilis i.e. a ew microns). Our results thus provide a new insight into the nature o bacterial motility, which goes beyond the mechanochemical description o chemotaxis and mechanisms o lagellar bundle ormation and disruption. MTERILS ND METHODS Genetic manipulation o bacteria llelic replacement was perormed using the pmdsystem according to [43]. ll strains originated rom the recently described NCIB3610 Bacillus subtilis strain harbouring a point mutation in the comi gene [44]. Briely, the hag gene encoding lagellin was ampliied including up- and downstream lanking regions and cloned into the pmd plasmid. The T209 to C mutation was obtained by quick-change mutagenesis. NCIB3610 comi Q12I was transormed and positive clones selected according to [43]. Exchange o the native hag gene to the hag T209C was conirmed by sequencing and light microscopy. The strain is reerred to as wild-type NCIB3610 in terms o the lagellar number in the text. We are grateul to Daniel Kearns or strain DK1693 (, amye::p hyspank ; lac::p hag hag T209C mls). This strain was used as both and ++ when induced with IPTG. Sample preparation The temperature was set to 37 C or all steps o the experiment. First, 20 ml o rozen stock Bacillus subtilis was streaked onto a LB-agar plate. The plate was incubated or 16 hours and then a ew colonies o bacteria rom the plate were stirred in LB and grown over night. The cultures were diluted to OD 600 0.1 next morning and grown or two more hours to reach the early exponential phase within the optical density range o 0.5 and 0.8. ter dilution, strain was induced by 1 mm IPTG solution to synthesize more lagella. The optical density othe cells in the exponentialphase wasirst adjusted to 0.5 and then urther diluted ( 15 olds) in a resh LB that had been previously puriied by 0.4 µm syringe ilter. Because o strong sticking o strain to the surace, 0.005% PVP-40 (polyvinylpyrrolidones) was added beore the experiments [45]. For luorescent labeling o cells, the dye was prepared by solving 1 mg lexa Fluor R 568 C5-maleimide in 200 µl DMSO (Dimethyl suloxide). 1 ml o the cell culture in OD 600 1 was centriuged at 8000g or one minute and gently washed three times in PBS ph 7.4(1X). The pellet was resuspended in 200 µl PBS together with 5 µl o the dye solution. The suspension was mixed and incubated in the dark at room temperature or 20 min and washed three times again in PBS buer and re-energized by hal hour outgrowth in LB to observe motile cells. Fluorescence microscopy samples were prepared by adding 30 µl o labeled cell suspension into a FluoroDish FD35-100 and covering it by circular cover glass (VWR, diameter 22 mm, thickness NO. 1.5). Nikon Eclipse Ti microscope together with a Nikon N Plan po λ 60x, N.. 1.4 oil immersion objective were used or luorescence microscopy. The dye was excited by a mercury lamp and images were acquired using Hamamatsu ORC-Flash 4.0 camera with 25 ms exposure time and 2 2 binning. Tracking chamber was a superposition o a cover slip (VWR, 20 20 mm, thickness NO. 1) on a microscope slide (Carl Roth GmbH, Karlsruhe, 76 26 mm) separated by a thin layer o silicone grease (GE Bayer Silicones Baysilone, medium viscosity) as spacer. The chambers were quasi-two-dimensional with a height ranging between 30 and 50 µm. The lateral sides were sealed by silicone grease ater illing the cavity with bacterial suspension. Microscopy was perormed with a Nikon Eclipse TE2000-s microscope and a Nikon 4x, N.. 0.2 objective in the dark ield mode. For each sample, sequences o images or 2 min with 60 Hz rame rate were recorded using a Point Grey FL3-U3-88s2cc camera. The experiments were repeated with three dierent cultures or each strain. Trajectory selection and tumbling detection triangular smoothing ilter was applied to smooth the data. Circular trajectories in the vicinity o the surace, extremely short tracks, and abnormal ones belonging to the cells in late exponential or dividing phase were discarded. More than 2500 cell trajectories were eventually selected or urther analysis. The tumbling detection algorithm was based on identiying the dramatic changes in speed v(t) and magnitude o angular velocity ω(t) o the bacterium [33, 50]. First, all local minima (maxima) o v (ω) over time were deter- FIG. 5. Sample trajectory with detected tumbling events. () The trajectory is color coded with respect to speed. (B) The detected tumbling events are indicated with red color. The arrow shows the direction o swimming. B
8 mined. Each minimum (maximum) is surrounded by two maxia (minima) located at t 1 and t 2. The depth v o the minimum or the height ω o the maximum can be characterized as v=max[v(t 1 ) v(t min ),v(t 2 ) v(t min )] or ω = max[ω(t max ) ω(t 1 ),ω(t max ) ω(t 2 )], respectively. To identiy a tumbling event two criteria were imposed: (i) I v/v(t min ) 0.7, there can be a tumbling phase around t min. The possible tumbling period is limited to those times t around t min where v(t) v(t min ) 0.2 v. (ii) The local maximum o ω may indicate a tumbling event i the total directional change during the time interval t 2 t 1 exceeds 0.8(t 2 t 1 ). The corresponding tumbling period consists o those times t around t max where ω(t max ) ω(t) ω. Imposing both conditions to identiy a tumbling event ensures that the reduction o speed is accompanied by a sudden change in the direction o motion. See [33, 50] or more details. sample trajectory with detected tumbling eventsis shown in Fig.5 (see alsosuppl. Movie S1). nalytical approach To describe the persistent motion o bacteria interrupted by stochastic tumbling periods, we developed an analytical ramework or a random walker with two states o motility, as described by the set o master Eqs. 1. In the ollowing, we briely explain how arbitrary moments o displacement can be obtained by a Fourier-z-transorm technique. Here, a two-dimensional motion is considered but extension to three dimensions is straightorward. The Fourier transorm o P j t (x,y γ) is deined as π P j t (ω m) dγ e imγ dy dx e iω r P j t (x,y γ), (3) with j {R, T}. n arbitrary moment o displacement x k1 y k2 j (t) can be obtained as x k1 y k2 j (t) dγ dy dx x k1 y k2 P j t (x,y γ) (4) P j =( i) k1+k2 k1+k2 t (ω x,ω y m=0). ωx k1 ω k2 y (ωx,ω y)=(0,0) The master Eqs. 1 ater Fourier transormation read k= (5) P R (ω,α m) = t+ t [ i k e ikα J k (ωv R t) R(m+k)Pt T (ω,α m+k) +(1 )p(m+k)p R t (ω,α m+k) ], P T (ω,α m) = t+ t [ i k e ikα J k (ωv T t) Pt R (ω,α m+k) k= +(1 )P T t (ω,α m+k) ], (6) (7) where we used the kth order Bessel s unction J k (z) = 1 2πi k π dγ e izcosγ e ikγ, (8) the Fourier transorms o the turning-angle distribution R(m) = π dφ e imφ R(φ), (9) and the distribution o the directional change along the run trajectory p(m) = π dθ e imθ R(θ). (10) The Fourier transorm o the probability P j t (ω,α m) can be expanded as a Taylor series P j t (ω,α m) = Qj 0,t (α m)+iωv j t Qj 1,t (α m) 1 2 ω2 v 2 j ( t)2 Q j 2,t (α m)+, (11) and the moments o displacement can be read in terms o the Taylor expansion coeicients. For example, x 2 j (t) =v 2 ( t)2 Q j j 2,t (0 0). (12) Thus we expand both sides o the master equations 6 and 7 and collect all terms with the same power in ω. s a result, coupled recursion relations or the Taylor expansion coeicients o terms with the same power in ω can be obtained. Next, the time indices on both sides o these equations can be equalized by means o z-transorm, which enables us to obtain the moments o displacement in the z-space. For example, one obtains the ollowing exact expression or the mean square displacement x 2 (z)= z t x 2 (t)= t=0 ( t) 2[ ] v 2 QR R 2 (z,0 0)+v 2 QT T 2(z,0 0) = ( t) 2[ (1 )Q R 0 (z,0 0)+ QT 0 ] (z,0 0) [ [ ] ] z z (1 ) v 2 R (z 1)G(z) + z (z 1)G(z) v v 1 R T 2(z 1) v2 R +( t) 2[ ] Q R 0 (z,0 0)+(1 )Q T 0(z,0 0) [ [ ] ] z z (1 )p v 2 T (z 1)G(z) + z (z 1)G(z) Rv Rv T 1 2(z 1) v2, T (13) where [ ][ ] G(z) = z (1 ) z (1 )p R. (14)
9 By inverse z-transorming the moments o displacement in the z-space, the moments can be obtained as a unction o time. SUPPLEMENTRY MTERILS movie S1. sample bacterial trajectory. [1] H. C. Berg, E. coli in motion (Springer Verlag, New York), (2004). [2] J. Schuhmacher, K. M. Thormann, G. Bange, How bacteria maintain location and number o lagella. FEMS Rev. Microbiol. 39, 812-822 (2015). [3] H. C. Berg, R.. nderson, Bacteria swim by rotating their lagellar ilaments. Nature 245, 380-382 (1973). [4] E. Lauga, T. R. Powers, The hydrodynamics o swimming microorganisms. Rep. Prog. Phys. 72, 096601 (2009). [5] H. C. Berg, D.. Brown, Chemotaxis in escherichia coli analysed by three-dimensional tracking. Nature 239, 500-504 (1972). [6] O. Benichou, M. Coppey, M. Moreau, P. H. Suet, R. Voituriez, Optimal search strategies or hidden targets. Phys. Rev. Lett. 94, 198101 (2005). [7] F. Bartumeus, S.. Levin, Fractal reorientation clocks: Linking animal behavior to statistical patterns o search. Proc. Natl. cad. Sci. 105, 19072 (2008). [8] O. Benichou, C. Loverdo, M. Moreau, R. Voituriez, Intermittent search strategies. Rev. Mod. Phys. 83, 81-129 (2011). [9] G. H. Wadhams, J. P. rmitage, Making sense o it all: bacterial chemotaxis. Nat. Rev. Mol. Cell Biol. 5, 1024-1037 (2004). [10] S. Chattopahyay, R. Moldovan, C. Yeung, X. L. Wu, Swimming eiciency o bacterium escherichia coli. Proc. Natl. cad. Sci. 103, 13712-13717 (2006). [11] N. C. Darnton, L. Turner, S. Rojevsky, H. C. Berg, On torque and tumbling in swimming escherichia coli. J. Bacteriol. 89, 1756-1764 (2007). [12] P. J. Mears, K. Santosh, C. V. Rao, I. Golding, Y. R. Chemla, Escherichia coli swimming is robust against variations in lagellar number. Elie 3, e01916 (2014). [13] L. Turner, W. S. Ryu, H. C. Berg, Real-time imaging o luorescent lagellar ilaments. J. Bacteriol. 182, 2793-2801 (2000). [14] E.. Korobkova, T. Emonet, H. Park, P. Cluzel, Hidden stochastic nature o a single bacterial motor. Phys. Rev. Lett. 96, 058105 (2006). [15] S. Terasawa, H. Fukuoka, Y. Inoue, T. Sagawa, H. Takahashi,. Ishijima, Coordinated reversal o lagellar motors on a single escherichia coli cell. Biophys. J. 100, 2193-2200 (2011). [16] B. Hu, Y. Tu, Coordinated switching o bacterial lagellar motors: evidence or direct motor-motor coupling? Phys. Rev. Lett. 110, 158703 (2013). [17] J. Hu, M. Yang, G. Gompper, R. G. Winkler, Modelling the mechanics and hydrodynamics o swimming E. coli. Sot Matter 11, 7867-7876 (2015). [18] N. Vladimirov, D. Lebiedz, V. Sourjik, Predicted auxiliary navigation mechanism o peritrichously lagellated chemotactic bacteria. PLoS Comput. Biol. 6, e1000717 (2010). [19] N. W. Sneddon, W. Pontius, T. Emonet, Stochastic coordination o multiple actuators reduces latency and improves chemotactic response in bacteria. Proc. Natl. cad. Sci. 109, 805-810 (2012). [20] P. Kanehl, T. Ishikawa, Fluid mechanics o swimming bacteria with multiple lagella. Phys. Rev. E 89, 042704 (2014). [21] S. Y. Reigh, R. G. Winkler, G. Gompper, Synchronization and bundling o anchored bacterial lagella. Sot Matter 8, 4363-4372 (2012). [22] D. B. Kearns, R. Losick, Cell population heterogeneity during growth o Bacillus subtilis. Genes Dev. 19, 3083-3094 (2005). [23] S. B. Guttenplan, S. Shaw, D. B. Kearns, The cell biology o peritrichous lagella in Bacillus subtilis. Mol. Microbiol. 87, 211-229 (2013). [24] S. Burov et al., Distribution o directional change as a signature o complex dynamics. Proc. Natl. cad. Sci. 110, 19689-19694 (2013). [25] M. R. Shaebani, Z. Sadjadi, I. M. Sokolov, H. Rieger, L. Santen, nomalous diusion o sel-propelled particles in directed random environments. Phys. Rev. E 90, 030701(R) (2014). [26] Z. Sadjadi, M. R. Shaebani, H. Rieger, L. Santen, Persistent random walk approach to anomalous transport o selpropelled particles. Phys. Rev. E 91, 062715 (2015). [27] L. Turner, L. Ping, M. Neubauer, H. C. Berg, Visualizing lagella while tracking bacteria. Biophys. J. 111, 630-639 (2016). [28] B. Rodenborn, C.-H. Chen, H. L. Swinney, B. Liu, H. P. Zhang, Propulsion o microorganisms by a helical lagellum. Proc. Natl. cad. Sci. 110, E338-E347 (2013). [29] C. Valeriani, M. Li, J. Novosel, J. rlt, D. Marenduzzo, Colloids in a bacterial bath: simulations and experiments. Sot Matter 7, 5228-5238 (2011). [30] Y. Hyon, T. R. Powers, R. Stocker, H. C. Fu, The wiggling trajectories o bacteria. J. Fluid Mech. 705, 58-76 (2012). [31] E. M. Purcell, Lie at low Reynolds number. m. J. Phys. 45, 3-11 (1977). [32] J. Saragosti, P. Silberzan,. Buguin, Modeling E. coli tumbles by rotational diusion: Implications or chemotaxis. PLoS ONE 7, e35412 (2012). [33] M. Theves, J. Taktikos, V. Zaburdaev, H. Stark, C. Beta, bacterial swimmer with two alternating speeds o propagation. Biophys. J. 105, 1915-1924 (2013). [34] J. Taktikos, H. Stark, V. Zaburdaev, How the motility pattern o bacteria aects their dispersal and chemotaxis. PLoS ONE 8, e81936 (2013). [35]. E. Haner, L. Santen, H. Rieger, M. R. Shaebani, Runand-pause dynamics o cytoskeletal motor proteins. Sci. Rep. 6, 37162 (2016). [36] I. Pinkoviezky, N. S. Gov, Transport dynamics o molecular motors that switch between an active and inactive state. Phys. Rev. E 88, 022714 (2013). [37] M. R. Shaebani,. E. Haner, L. Santen, Geometrical considerations or anomalous transport in complex neuronal dendrites. submitted (2017). [38] M. Kong, Y. Wu, G. Liab, R. G. Larson, bead-spring model or running and tumbling o lagellated swimmers: detailed predictions compared to experimental data or E. coli. Sot Matter 11, 1572-1581 (2015). [39] N. Watari, R. G. Larson, The hydrodynamics o a run-andtumble bacterium propelled by polymorphic helical lagella. Biophys. J. 98, 12-17 (2010). [40] Z. Sadjadi, M. F. Miri, M. R. Shaebani, S. Nakhaee, Diusive transport o light in a two-dimensional disordered packing o disks: nalytical approach to transport mean ree path. Phys. Rev. E 78, 031121 (2008). [41] Z. Sadjadi, M. F. Miri, Diusive transport o light in twodimensional granular materials. Phys. Rev. E 84, 051305 (2011). [42] P. Tierno, M. R. Shaebani, Enhanced diusion and anomalous transport o magnetic colloids driven above a two-state lashing potential. Sot Matter 12, 3398-3405 (2016). [43] O. Benichou, C. Loverdo, M. Moreau, R. Voituriez, Twodimensional intermittent search processes: n alternative to Levy light strategies. Phys. Rev. E 74, 020102 (2006). [44] M. Chabaud et al., Cell migration and antigen capture are antagonistic processes coupled by myosin II in dendritic cells.
10 Nat. Commun. 6, 7526 (2015). [45] J. F. Rupprecht, O. Benichou, R. Voituriez, Optimal search strategies o run-and-tumble walks. Phys. Rev. E 94, 012117 (2016). [46] V. Tejedor, R. Voituriez, O. Benichou, Optimizing persistent random searches. Phys. Rev. Lett. 108, 088103 (2012). [47] M. rnaud,. Chastanet, M. Débarbouillé, New vector or eicient allelic replacement in naturally nontransormable, low- GC-content, gram-positive bacteria. ppl. Enviro. Microbiol. 70, 6887-6891 (2004). [48] M.. Konkol, K. M. Blair, D. B. Kearns, Plasmid-Encoded ComI Inhibits Competence in the ncestral 3610 Strain o Bacillus subtilis. J. Bacteriol. 195, 4085-4093 (2013). [49]. P. Berke, L. Turner, H. C. Berg, E. Lauga, Hydrodynamic attraction o swimming microorganisms by suraces. Phys. Rev. Lett. 101, 038102 (2008). [50] J. B. Masson, G. Voisinnea, J. Wong-Nga,. Celania, M. Vergassolaa, Noninvasive inerence o the molecular chemotactic response using bacterial trajectories. Proc. Natl. cad. Sci. 109, 1802-1807 (2012). cknowledgments: We thank Daniel B. Kearns or providing us with strain DK1693, and Z. Sadjadi or ruitul discussions. Funding: M.R.S. acknowledges supportby DFG within SFB 1027 (7). uthor contributions: C. W. and G. B. designed the research. F.. prepared the strains o the bacterium Bacillus subtilis. C. W., T.J., and J.N. designed and built the experimental setup. J.N. carried out the experiments and analyzed the results. M. R. S. developed the analytical model and perormed simulations. ll authors contributed to the interpretation o the results. M. R. S. drated the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: ll data needed to evaluate the conclusions in the paper are present in the paper. dditional data related to this paper may be requested rom the authors.