Dynamic density shaping of photokinetic E. coli
Abstract
Many motile microorganisms react to environmental light cues with a variety of motility responses guiding cells towards better conditions for survival and growth. The use of spatial light modulators could help to elucidate the mechanisms of photomovements while, at the same time, providing an efficient strategy to achieve spatial and temporal control of cell concentration. Here we demonstrate that millions of bacteria, genetically modified to swim smoothly with a light controllable speed, can be arranged into complex and reconfigurable density patterns using a digital light projector. We show that a homogeneous sea of freely swimming bacteria can be made to morph between complex shapes. We model nonlocal effects arising from memory in light response and show how these can be mitigated by a feedback control strategy resulting in the detailed reproduction of grayscale density images.
https://doi.org/10.7554/eLife.36608.001eLife digest
Many bacteria can move in response to environmental signals. This helps guide them towards better conditions for growth and survival. Escherichia coli is a bacterium that can swim quickly through liquid, using tiny propellerlike structures that rotate many times per second. These ‘propellers’ are powered by a cellular motor, called the flagellar motor, which similar to an electric motor, requires an energy source to drive movement.
Proteorhodopsin, a protein originally isolated from freeswimming microorganisms in the ocean, is an alternative energy source that helps bacteria move. The protein is located close to the surface of the cell, where it acts like a solar panel and captures energy from light. In cells powered by proteorhodopsin, the intensity of light from their environment determines their swimming speed: brighter light means faster movement, and less light, slower movement. Proteorhodopsin is now also a useful tool in the laboratory. For example, genetically engineering bacteria to produce proteorhodopsin provides a way to control their movement remotely, using a light source.
Swimming bacteria, much like cars in city traffic, are known to accumulate in areas where their speed decreases. By controlling swimming speed with proteorhodopsin, researchers can manipulate the local density of bacteria simply by projecting different patterns of light.
To study the factors influencing this phenomenon, Frangipane et al. used genetically modified E. coli that could respond to light via proteorhodopsin to make layers of cells that could then have light patterns projected onto them. The results showed that the bacteria responded slowly to these stimuli, which was the main factor limiting the resolution of the final pattern they formed. A simple feedback mechanism, which compared the pattern formed by the cells to the desired image and updated the projected light accordingly, was enough to solve this problem. This way, the layers of E.coli could be turned into a nearperfect copy of the original image.
This work allows us to control the movement of large populations of bacteria more precisely than ever before. This could be extremely valuable for building the next generation of microscopic devices. For example, bacteria could be made to surround a larger object such as a machine part or a drug carrier, and then used as living propellers to transport it where it is needed.
https://doi.org/10.7554/eLife.36608.002Introduction
A local reduction of speed in a crowd of walking pedestrians or urban car traffic usually results in the formation of highdensity regions. Mathematically speaking, a selfpropelled particle, moving by an isotropic random walk with space dependent speed, explores the available space with a probability density that is inversely proportional to the local value of the speed (Schnitzer, 1993; Tailleur and Cates, 2009; Cates and Tailleur, 2015), spending a longer time in slow regions than in fast ones. As a consequence, any mechanism that allows to substantially modulate the local propulsion speed can be used as an efficient strategy for controlling the density (Stenhammar et al., 2016). In swimming bacteria, proteorhodopsin (PR) provides a particularly convenient way of achieving spatial speed control through the projection of light patterns. PR is a lightdriven proton pump that uses photon energy to pump protons out of the cytoplasm, thus modulating the corresponding electrochemical gradient across the inner membrane (Béjà et al., 2000). Since this proton motive force drives the rotation of the flagellar motor (Gabel and Berg, 2003), PR puts a ‘solar panel’ on every cell allowing to remotely control its swimming speed with light (Walter et al., 2007). This has been recently exploited to control the rotational speeds of biohybrid micromachines using bacteria as micropropellers (Vizsnyiczai et al., 2017). In light driven Janus colloids, spatially asymmetric speed profiles can also affect particle orientation giving rise to an artificial phototaxis mechanism (Lozano et al., 2016). Arlt et al. (2018) have recently shown that, by projecting a masked illumination pattern, ‘photokinetic’ (Wilde and Mullineaux, 2017) bacteria can be accumulated in dark regions and depleted from brighter ones thus forming binary patterns. Here we demonstrate that bacterial density can be controlled to form reconfigurable complex patterns. To do this we exploit a Digital Light Processing (DLP) projector to shape light with a megapixel resolution and with a dynamic range of 256 (8 bit) light intensity levels. We experimentally investigate the limits of validity of the predicted law of inverse proportionality between local density and speed. We demonstrate that observed deviations from this law can be quantitatively described by a theoretical model that explicitly takes into account memory effects in light response. Finally we show that a model independent feedback control loop allows density shaping with high spatial resolution and gray level accuracy.
Results
Density shaping in photokinetic bacteria or colloids relies on the connection between density and local speed. The density $\rho (\mathbf{\mathbf{r}})\propto 1/v(\mathbf{\mathbf{r}})$ is always a steadystate solution of the master equation of a selfpropelled particle with isotropic reorientation dynamics and spatially varying speed $v(\mathbf{\mathbf{r}})$ (Cates and Tailleur, 2015). In our case, $v(\mathbf{\mathbf{r}})=v(I(\mathbf{\mathbf{r}}))$ where $\mathbf{\mathbf{r}}=(x,y)$ is the position in the two dimensional image space and $I(\mathbf{\mathbf{r}})$ the light intensity at $\mathbf{\mathbf{r}}$. This result assumes that all bacteria have the same response to light and that they instantaneously adapt to intensity variations so that the speed only depends on the local value of $I$. However, suspensions of swimming bacteria, are known to be characterized by motility properties that occur in broad distributions and, in principle, a broad spectrum of light responses could be found. To quantify the effect of light on the speed distributions of our strain, we monitored bacterial dynamics while projecting a chessboard pattern consisting of 12 different levels of light intensity. To do this we couple a DLP projector to a custom videomicroscopy setup working in both brightfield and darkfield mode (see Figure 1 and Materials and methods).
Using differential dynamic microscopy (DDM) (Cerbino and Trappe, 2008; Wilson et al., 2011; Maggi et al., 2013) we extract the light dependent values of the mean $v$ and standard deviation $\sigma $ of the speed distributions (Figure 2(a)) (see Materials and methods). Although the absolute value of the swimming speed depends on several factors (strain, growth medium, etc.) we find a non linear speed response that is consistent with what previously reported for PR expressing bacteria. In particular the speed versus intensity curve is very well fitted by a hyperbola showing saturation for large light intensity values as already reported in (Walter et al., 2007; SchwarzLinek et al., 2016; Vizsnyiczai et al., 2017; Arlt et al., 2018).

Figure 2—source data 1
 https://doi.org/10.7554/eLife.36608.005
We also find that, throughout the intensity range, $v$ and $\sigma $ are directly proportional (see Figure 2(b)) suggesting a homogeneous growth law for individual cell speeds ${v}_{i}(I)={v}_{i}^{s}f(I)$, where ${v}_{i}^{s}$ is the saturation speed for the $i$th cell and $f(I)$ is a dimensionless function saturating at one for large intensities. Interestingly, in this scenario, the mean speed $v(I)$ will also be proportional to $f(I)$ so that the probability density ${\rho}_{i}$ for the generic $i$th cell will be proportional to the inverse of the local mean speed:
This implies that, despite broad speed distributions, the local density is expected to be only determined by the mean local speed, as set by the light intensity pattern and measured by DDM. To validate this hypothesis we use a DLP projector to display a complex light pattern onto a 400 µm thick layer of cells (Figure 1) that have been preliminarily exposed to a uniform bright illumination for 5 min. This time is much longer than the speed response time of bacteria (Figure 4) thus ensuring that cells are initialized to swim at maximal speed. The projecting system has an optical resolution of 2 µm approximately matching the size of a single cell which represents the physical ‘pixel’ of our density images. This value sets the limit for the minimum theoretical resolution of density configurations that would be achievable if bacteria could be precisely and statically arranged in space. In practice, as we will see, the real resolution will always be larger for two main reasons: (i) bacteria do not respond instantly to light temporal variations thus introducing a blur in the target speed map, (ii) the stationary state is an ensemble of noisy patterns that constantly fluctuate because of swimming and Brownian motions of bacteria. After projecting an inverted Mona Lisa image, bacteria start to concentrate in dark regions while moving out from the more illuminated areas. After about 4 min, dark field microscopy reveals a recognizable bacterial ‘replica’ of Leonardo’s painting, where brighter areas correspond to regions of accumulated cells. Once a stationary pattern is reached, we collect a time averaged image from which we extract the bacterial density ${\rho}^{*}$ shown in Figure 3 normalized to be one when uniform (see Materials and methods).

Figure 3—source data 1
 https://doi.org/10.7554/eLife.36608.007
We split the density image into components corresponding to the same illumination level $I$. We then average the density within each component and report in Figure 3(c) the obtained value as a function of the inverse mean speed $1/v(I)$ obtained from the speed vs intensity curve in Figure 2(a). The high speed side of the graph can be very well fitted by a straight line with a finite intercept $q=0.5$. This suggests that 50% of the total scattering objects in the field of view responds to light and can be spatially modulated while the remaining 50% can be attributed to cells which are nonmotile or insensitive to light and also to stray light generated away from the focal plane. The ratio between the maximum and minimum modulation of the light sensitive component can be obtained as $\text{max}[{\rho}^{*}q]/\text{min}[{\rho}^{*}q]=2$. A strong deviation from linearity is evident in the high density/low speed region. A violation of the $\rho \propto 1/v$ law can be attributed to many factors that are not included in the simple theory discussed above. In particular, the theory assumes that speed is a local function of space which would only be the case if bacteria instantly adapt to temporal changes in light intensity. However, previous studies have evidenced that speed response is not instantaneous but it displays a relaxation pattern characterized by multiple timescales. In addition to a fast relaxation time associated to the membrane discharge (Walter et al., 2007), ATP synthase and the dynamics of stator units of the flagellar motors introduce slower timescales in the speed relaxation dynamics (Tipping et al., 2013; Arlt et al., 2018). Using DDM we measured speed response to a uniform light pattern whose intensity switches between two values every 100 s. Results are reported in Figure 4 and clearly show the presence of two timescales. The short time scale is smaller than our time resolution (1 s) and appears as an instantaneous jump that accounts for about a fraction $\beta =0.44$ of the total relaxation. A slower relaxation follows and it is well fitted by an exponential function with a time constant ${\tau}_{\text{m}}=35$ s that is the same for both the rising and falling relaxations. As a result, bacteria will experience an effective speed map that is a blurred version of what we would expect for an instantaneous response $V(\mathbf{\mathbf{r}})=v(I(\mathbf{\mathbf{r}}))$. In the case of smooth swimming cells, with a twostep light response, we calculate that, for weak speed modulations around a baseline value ${V}_{0}$, the actual speed map can be obtained as a simple convolution (see Materials and methods):

Figure 4—source data 1
 https://doi.org/10.7554/eLife.36608.009
with $\gamma (x,y)$ a convolution kernel that is analytic in Fourier space:
where $q=\sqrt{{q}_{x}^{2}+{q}_{y}^{2}}$ is the modulus of the wave vector and $k}^{1}={v}_{0}{\tau}_{m$. Assuming that the stationary distribution will remain isotropic even in the presence of memory, the relation $\rho \propto 1/w$ will then still be valid provided one uses the actual, blurred speed map $w$ and not the original map $V$ corresponding to instantaneous response. We can then anticipate the stationary density once we calculate the actual speed map $w$ with the convolution formula (1). Recalling that only a fraction $\alpha $ of cells is motile the estimated normalized stationary density will be
where $\overline{{w}^{1}}$ is the spatial average of the inverse actual speed $w(x,y)$. Plotting ${\rho}^{*}$ as a function of $1/V$ we find that this simple model quantitatively explains the density behavior in Figure 3 with the natural choice of free parameters $\alpha =0.5$, $\beta =0.44$, and ${k}^{1}=\overline{v}{\tau}_{m}=59$ µm where $\overline{v}$ is the spatial average of $V(x,y)$. This result is certainly encouraging and prompts for a deeper and more systematic investigation using simpler light patterns. Our aim here, however, is to investigate the resolution limits of density shaping with photokinetic bacteria. Although memory effects could be reduced by deleting ATPsynthase genes (Arlt et al., 2018), resulting in a suppression of the slow component in speed relaxation, response will never be truly instantaneous. In addition to memory effects, steric and/or hydrodynamic interactions will always affect density by hindering bacteria penetration and further accumulation in densely packed regions. The combined action of these effects on the stationary density map cannot be easily incorporated in a theoretical model that provides precise quantitative predictions with a priori known parameters.
An alternative and more practical strategy to improve the accuracy of density shaping, is to implement a feedback control loop. We set up an automated control loop that performs one iteration every 20 s, comparing the current density map to the desired target image and updating the illumination pattern accordingly (see Materials and methods). Light levels are increased in those regions with a density that exceeds the target value and decreased where density is lower. At each step the light intensity increment is simply proportional to the pixel by pixel difference of target and actual density images. We quantify the distance between the obtained density pattern and the target image by the sum of the squared pixelbypixel difference after a proper rescaling procedure (see Materials and methods).
As soon as the feedback is turned on at time ${t}_{0}=6$ min (Figure 5(a)) the distance from the target starts to decrease reaching a new stationary value after about 4 min. The final stationary density map shows a level of detail significantly higher than the image obtained before the feedback was turned on (Figure 5(b,c)). The image is stable and resolution does not deteriorate for hours. However a dark region surrounds the illuminated area so that swimming cells that exit the field of view do not come back and the number of motile cells in the illuminated area is reduced by about 50% in one hour. This issue is the main limiting factor for the lifetime of our patterns and could be solved in the future by using sample chambers that match the size of the illuminated area. As evidenced in Figure 5(d), the feedback precompensates this blurring effect due to memory by converging to a sharpened version of the initial target image. Results of comparable quality have been replicated four times in each of the two independent biological replicates.

Figure 5—source data 1
 https://doi.org/10.7554/eLife.36608.011
A remarkable feature of this optically controlled material lays in its intrinsic dynamic and reconfigurable nature. As a demonstration of this property we show the dynamic morphing of a bacterial layer from an Albert Einstein’s to a Charles Darwin’s portrait (see Figure 6 and Video 1). The morphing is triggered by instantly switching between the two target images on the DLP projector and exponentially converges to the second target in a characteristic time of 1.8 min.

Figure 6—source data 1
 https://doi.org/10.7554/eLife.36608.013
Discussion
Although light induced density modulations can be achieved in large samples of Brownian (passive) colloidal systems, active materials provide a largely superior performance in terms of both required power and response time. The stationary density distribution of Brownian colloids is governed by the Boltzmann law stating that the logarithmic ratio of maximum to minimum density is given by
where $\mathrm{\Delta}U$ is the optical energy difference between brightest and darkest regions. For a plastic or glass bead of radius $a=1$µm in water, this energy difference can be estimated in the Rayleigh regime as (Ashkin et al., 1986) $\mathrm{\Delta}U\approx {a}^{3}I/c$ where $I$ is the maximum power density and $c$ is the speed of light. For a contrast level comparable to the ones shown above, $\mathrm{log}({\rho}_{\mathrm{max}}/{\rho}_{\mathrm{min}})\approx 1$, the required power density would be $I\approx {k}_{B}Tc/{a}^{3}\approx 1$ W/mm${}^{2}$. This value is three orders of magnitude larger than the maximum power densities used in our experiments. Our active system is also much faster than its Brownian counterprt as evidenced by comparing the timescales of pattern formation. In both the active and passive case, dynamics is governed by the interplay of diffusion and drift with characteristic timescales ${\tau}_{\mathrm{diff}}$ and ${\tau}_{\mathrm{drift}}$. Calling $\mathrm{\ell}=1$ mm the largest length scale in the target pattern, we obtain for the active case ${\tau}_{\mathrm{drift}}\approx \mathrm{\ell}/v\approx 200$ s where we used $v=5$ µm/s as the typical swimming speed. A diffusion time scale can be obtained as ${\tau}_{\mathrm{diff}}\approx {\mathrm{\ell}}^{2}/2D$ where $D$ is the active translational diffusion coefficient. For wild type bacteria, moving with a run and tumble dynamics, $D={v}^{2}{\tau}_{\text{run}}$ where ${\tau}_{\text{run}}\approx 1$ s is the mean duration of a run. The corresponding characteristic diffusion time is ${\tau}_{\mathrm{diff}}\approx {\mathrm{\ell}}^{2}/2D\approx 2\times {10}^{4}$ s. This timescale can be considerably reduced if, as in our case, we suppress tumbling and use a smooth swimming strain with a much larger reorientation time ${\tau}_{\mathrm{rot}}\approx $20 s (Berg, 1993). In this case the diffusion timescale drops down to ${\tau}_{\mathrm{diff}}\approx {10}^{3}$ s leading to a faster relaxation towards the stationary state. In the Brownian case the drift timescale is given by ${\tau}_{\mathrm{drift}}\approx \mu {k}_{B}T/\mathrm{\ell}\approx $10${}^{6}$ s where $\mu =50$ µm/s pN is the mobility of a 1 µm radius microsphere. The same colloidal particle has a diffusivity $D=\mu {k}_{B}T=0.2$ µm${}^{2}$/s giving again ${\tau}_{\mathrm{diff}}\approx {\mathrm{\ell}}^{2}/2D\approx {\mathrm{\hspace{0.33em}10}}^{6}$ s. Summarizing, using optical forces to achieve a comparable density modulation for Brownian particles would require a thousand times larger powers and from ${10}^{3}$ to ${10}^{4}$ longer times.
In conclusion, we have shown that a suspension of swimming bacteria, with optically controllable speed, can provide a new class of light controllable active materials whose density can be accurately, reversibly and quickly shaped by employing a low power light projector. An alternative strategy to produce static patterns of bacteria is by ‘biofilm lithography’ where an optical template of blue light is used to induce the expression of membrane proteins that promote cellcell and cellsubstrate attachment (Jin and RiedelKruse, 2018). As opposed to our motility driven density shaping, these patterns are static and form over a time scale of several hours. Interestingly, the two techniques could be used in combination to achieve faster and more complex (e.g. three dimensional) biofilm lithography by fixing with blue light a template pattern obtained by speed modulation with green light. By further genetic engineering, bacteria could be eventually encapsulated in silicate shells (Müller et al., 2008) producing solid permanent structures for micromechanics or microoptics applications. Finally, the possibility of spatial and temporal control of motile bacteria density could also lead to novel strategies for the transport and manipulation of small cargoes inside microdevices Maggi et al. (2018).
Materials and methods
Production of photokinetic cells
Request a detailed protocolFor all the experiments we used the E. coli strain HCB437 (Wolfe et al., 1987) transformed with a plasmid encoding the PR under the control of the araCpBAD promoter (Biobricks, BBa_K1604010 inserted in pSB1C3 plasmid backbone). E. coli colonies from frozen stocks are grown overnight at 33$}^{\circ$C on LB agar plates supplemented with kanamycin (Kan 30 µg/mL) and chloramphenicol (Cam 20 µg/mL). A single colony is picked and statically cultivated overnight at 33$}^{\circ$C in 10 mL of M9 broth (M9 salts with 0.2% glucose, 0.2% casaminoacids) supplemented with antibiotics as before. The overnight culture is diluted 100fold into 5 mL of the previous medium, grown at 33$}^{\circ$C, 200 rpm. 5 mM arabinose and 20 µM retinal are added once OD_{590} ≈ 0.2, keeping the culture in the dark to avoid retinal degradation. Once OD_{590} ≈ 0.8, cells are collected by centrifugation (1500 rcf, 5’). The resulting pellet is washed twice by centrifugation (1500 rcf, 5’) with motility buffer (MB: 0.1 mM EDTA, 10 mM phosphate buffer and 0.2% Tween20). The bacterial suspension is eventually adjusted to a working OD${}_{590}\approx $ 2.0.
Sample preparation
Request a detailed protocol100 µL of the prepared bacterial suspension is injected into a glass capillary (CM Scientific  Rect. boro capillaries 0.40 $\times$ 4 mm), sealed at both sides with vacuum grease (SigmaAldrich). The capillary is then placed under the optical microscope and observed in bright and dark field illumination focusing on the plane in the middle of the capillary.
Optics
Request a detailed protocolBright field and dark field imaging were performed using a custom inverted optical microscope equipped with a 4 $\times $ magnification objective (Nikon; NA = 0.13) and a highsensitivity CMOS camera (Hamamatsu OrcaFlash 2.8) (see Figure 1). Light shaping was performed using a digital light processing (DLP) projector (Texas Instruments DLP Lightcrafter 4500) coupled to the same microscope objective used for imaging. The size of a (squared) DLP projector pixel imaged on the sample plane results to be 2 µm. All darkfield images have been divided by the flat field image obtained by applying a Gaussian filter on the respective homogeneous bacterial suspension (under uniform green illumination) as detailed in the following. Flat background subtraction (fifth percentile of the image histogram) was only applied to the frames appearing in Figure 6.
Differential Dynamic Microscopy (DDM)
Request a detailed protocolTo perform DDM we compute the quantity:
where $M(\mathbf{\mathbf{q}},t)$ is the spatial Fourier transform at the wavevector $\mathbf{\mathbf{q}}$ of the brightfield image captured at time $t$. Assuming timetranslational invariance and isotropy in the bacterial movement $g(\mathbf{\mathbf{q}},{t}^{\prime},t)$ depends only on the timelag $t$ in Equation (4) and on the modulus of the wave vector $q=\mathbf{\mathbf{q}}$, that is $g(\mathbf{\mathbf{q}},t,{t}^{\prime})=g(q,t)$. The $g(q,t)$ is connected to the intermediate scattering function (ISF) $F(q,t)$:
where $A(q)$ and $B(q)$ are timeindependent factors related, respectively, to the number and shape of bacteria, and to the background noise in the images. Following Wilson et al., 2011 we use the ISF model for independent smooth swimming cells:
where $\alpha $ is the fraction of motile cells, $D$ the Brownian diffusion coefficient and $P({v}^{\prime})$ the Schultz distribution:
Here $v$ is the mean speed, $\mathrm{\Gamma}$ is the Euler gamma function and $Z$ is related to the speed standard deviation $\sigma$ by the formula $Z=(v/\sigma {)}^{2}1$. The relevant parameters discussed in the main text $v$ and $\sigma $ are extracted by fitting the experimental $g(q,t)$ in the $q$range 0.45 μm^{1}<q<1.2 μm^{1}with Equations 5, 6 and 7.
Image analysis
Request a detailed protocolWe assume that dark field images are proportional to the bacterial density modulated by a slowly varying envelope due to inhomogeneities in illumination. This flat field correction ${\rho}_{0}(\mathbf{\mathbf{r}})$ is obtained by acquiring 100 frames at 25 fps at the beginning of the experiments when the bacterial density is homogeneous. These frames are averaged and then filtered with a Gaussian kernel with a large standard deviation ($\approx 100$µm). The normalized experimental density ${\rho}^{\ast}(\mathbf{\mathbf{r}})$ is obtained by first computing a raw density field $\rho (\mathbf{\mathbf{r}})$ (50 frames average at 25 fps) and then dividing by ${\rho}_{0}$, that is ${\rho}^{\ast}(\mathbf{\mathbf{r}})=\rho (\mathbf{\mathbf{r}})/{\rho}_{0}(\mathbf{\mathbf{r}})$.
For computing the distance between ${\rho}^{\ast}(\mathbf{\mathbf{r}})$ and the target density ${\rho}^{\mathrm{tar}}(\mathbf{\mathbf{r}})$ we first rescale ${\rho}^{\ast}(\mathbf{\mathbf{r}})$ so that the 10th and 90th percentile of the two image histograms coincide:
The distance $\mathrm{\U0001d5bd\U0001d5c2\U0001d5cc\U0001d5cd}$ is then computed as:
where the sum runs over the image pixels $\mathbf{\mathbf{r}}$.
Modeling the effects of memory in speed response
Request a detailed protocolA smooth swimming cell with a twostep response to light intensity and traveling in the $x$ direction over a light imposed speed profile $V(x)$, will have an actual speed at time $t$ given by:
where $x(t)$ is the position of the cell at time $t$. For a weakly varying speed profile $V(x)={V}_{0}+\delta V(x)$ we can transform the above time convolution into a space convolution by assuming $x(t)\approx x(0)+{V}_{0}t$:
where we have also assumed that an equal number of bacteria will be traveling in the opposite direction. Assuming isotropy and using the fact that, in the weak modulation limit, density is homogeneous at the zero order, we can easily generalize to the three dimensional case and write:
where $\mathrm{\Gamma}(\mathbf{\mathbf{r}})=(k/4\pi {r}^{2}){e}^{kr}$ is a 3D convolution kernel. Our projecting light system is such that light patterns do not vary significantly across the entire sample depth so that we can assume $V(x,y,z)=V(x,y)$ and express the effective speed distribution as a 2D convolution
with $\gamma (x,y)={\int}_{\mathrm{\infty}}^{\mathrm{\infty}}\mathrm{\Gamma}(x,y,z)\mathit{d}z$. Although the convolution kernel $\gamma $ is not analytic in real space, its Fourier transform is analytic making convolutions easy to calculate numerically:
where $q=\sqrt{{q}_{x}^{2}+{q}_{y}^{2}}$.
Feedback
Request a detailed protocolThe $(n+1)$th illumination pattern at the pixel $\mathbf{\mathbf{r}}$ is updated as follows:
where $P>0$ is a proportional control parameter and $\mathrm{\Delta}\rho $ is the difference between the scaled density and the target density:
If, for example, the scaled density is larger than the target at $\mathbf{\mathbf{r}}$ the projected power density at that pixel will be increased. That will increase the speed of bacteria and consequently reduce their local density.
References

Painting with lightpowered bacteriaNature Communications 9:768.https://doi.org/10.1038/s41467018031618

MotilityInduced phase separationAnnual Review of Condensed Matter Physics 6:219–244.https://doi.org/10.1146/annurevconmatphys031214014710

Differential dynamic microscopy: probing wave vector dependent dynamics with a microscopePhysical Review Letters 100:188102.https://doi.org/10.1103/PhysRevLett.100.188102

Phototaxis of synthetic microswimmers in optical landscapesNature Communications 7:12828.https://doi.org/10.1038/ncomms12828

Currents and fluxinversion in photokinetic active particlesSoft Matter 14:4958–4962.https://doi.org/10.1039/C8SM00788H

Motility fractionation of bacteria by centrifugationSoft Matter 9:10885–10890.https://doi.org/10.1039/c3sm51223a

Theory of continuum random walks and application to chemotaxisPhysical Review E 48:2553–2568.https://doi.org/10.1103/PhysRevE.48.2553

Escherichia coli as a model active colloid: A practical introductionColloids and Surfaces B: Biointerfaces 137:2–16.https://doi.org/10.1016/j.colsurfb.2015.07.048

Lightinduced selfassembly of active rectification devicesScience Advances 2:e1501850.https://doi.org/10.1126/sciadv.1501850

Quantification of flagellar motor stator dynamics through in vivo protonmotive force controlMolecular Microbiology 87:338–347.https://doi.org/10.1111/mmi.12098

Light controlled 3D micromotors powered by BacteriaNature Communications 8:15974.https://doi.org/10.1038/ncomms15974

Lightcontrolled motility in prokaryotes and the problem of directional light perceptionFEMS Microbiology Reviews 41:900–922.https://doi.org/10.1093/femsre/fux045

Differential dynamic microscopy of bacterial motilityPhysical Review Letters 106:018101.https://doi.org/10.1103/PhysRevLett.106.018101

Reconstitution of signaling in bacterial chemotaxisJournal of Bacteriology 169:1878–1885.https://doi.org/10.1128/jb.169.5.18781885.1987
Decision letter

Raymond E GoldsteinReviewing Editor; University of Cambridge, United Kingdom

Naama BarkaiSenior Editor; Weizmann Institute of Science, Israel
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Dynamic density shaping of photokinetic E. coli" for consideration by eLife. Your article has been reviewed by three peer reviewers, including Raymond E Goldstein as the Reviewing Editor and Reviewer #1, and the evaluation has been overseen by Naama Barkai as the Senior Editor.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Summary:
In this paper, the authors demonstrate the ability to dynamically arrange large collections of E. coli bacteria into complex twodimensional density patterns using light illumination using a digital projector. The method relies on E. coli mutants whose speed depends on the local level of illumination and is based on a previously proposed relationship between swimming speed and density in suspensions of selfpropelled particles. The authors further implement a feedback control strategy to improve resolution, where they adjust illumination based on current density. Their method is shown to be rapid (compared to methods based on passive colloids), to be able to capture fine spatial features, and to be reversible.
Essential revisions:
1) The role of runandtumble locomotion in this whole process should beexplained. Perhaps the authors could clarify the length/time scales associated with translational diffusion coming from that process in the directed pattern formation.
2) A more detailed discussion of resolution (and what controls/limits resolution) would be useful.
3) Presentation: Many concepts and methods warrant expanded explanation. Symbols are introduced without definition. There is inadequate description of experimental methods, a diagram of the experimental protocol would have been helpful, highlighting parameters such as size and geometry of chamber etc. The method of preparation of the bacterial suspension is also unclear. (In fact, supplementary figure 1 might be helpful as part of the main text.) Several algorithms for data extraction or image manipulation are glossed over and require careful definition.
4a) Interpretation: Discussion of the theoretical background can be improved. For instance, the assumed inverse proportionality law between local density and speed is itself a simplification, and the limitations of this simple theory could be discussed in the context of the experimental results. (The discussion of Figure 1B is ambiguous.)
4b) There is a growing body of work in which active matter is manipulated successfully with light, and these should be cited and discussed. Notably, recent papers from the RiedelKruse lab (Jin and RiedelKruse, 2018, PNAS), or the Poon group (Arlt et al., 2018, Nature Comms.), and several others besides. The authors should ensure previously published literature is clearly referenced, and importantly, how their work differs from or represents improvements upon these previous results.
5) Purpose: While this is a beautiful result, the wider significance of this work is unclear. We are unconvinced by the bioengineering argument that these active materials are somehow better/more controllable for micro optics applications; it is argued in terms of drift and diffusion timescales that active systems are much more responsive than their Brownian counterparts, this is undeniable – but why might it be desirable to control such samples in the first place? Furthermore, are these patterns stable? If continued to be subjected to the light, individual bacteria will adapt to the surrounding light conditions, how long will a typical Mona Lisa image remain recognisable before the pattern degrades?
We would have liked to have seen either (i) a clear theoretical interpretation in terms of the theory of motility induced phase separation, e.g. what are the physical reasons for the measured deviations from the simple inverse densityspeed relationship, or (ii) a clear biological interpretation of the timescales/transients involved in the pattern formation process in terms of some intracellular "memory". More interesting would be if the authors could measure, by combining precisely projected patterns with feedback modulation of population density, the actual spatiotemporal memory of the bacteria and how this might change with genetic manipulation.
https://doi.org/10.7554/eLife.36608.017Author response
Essential revisions:
1) The role of runandtumble locomotion in this whole process should beexplained. Perhaps the authors could clarify the length/time scales associated with translational diffusion coming from that process in the directed pattern formation.
We thank the reviewers for pointing that out. We now discuss translational
diffusion timescales for wild type (run and tumble) bacteria which
gives us the opportunity of highlighting the importance of using smooth
swimmers for faster density shaping. This is discussed in the revised
manuscript as follows:
“Calling ℓ = 1 mm the largest length scale in the target pattern, we obtain for the active case τ_{drift} ≈ ℓ/v ≈ 200 s where we used v = 5 μm/s as the typical swimming speed. A diffusion time scale can be obtained as τ_{diff} ≈ ℓ^{2}/2D where D is the active translational diffusion coefficient. For wild type bacteria, moving with a run and tumble dynamics, D = v^{2}τ_{run} where τ_{run} ≈ 1 s is the mean duration of a run. The corresponding characteristic diffusion time is τ_{diff} ≈ ℓ^{2}/2D ≈ 2 × 10^{4} s. This timescale can be considerably reduced if, as in our case, we suppress tumbling and use a smooth swimming strain with a much larger reorientation time τ_{rot} ≈ 20 s (Ashkin et al., 1986). In this case the diffusion timescale drops down to τ_{diff} ≈ 10^{3} s leading to a faster relaxation towards the stationary state.”
2) A more detailed discussion of resolution (and what controls/limits resolution) would be useful.
We now discuss resolution limits in greater detail:
“To validate this hypothesis we use a DLP projector to display a complex light pattern onto a 400 μm thick layer of cells (Figure 1) that have been preliminarily exposed to a uniform bright illumination for 5 min. This time is much longer than the speed response time of bacteria (Figure 4) thus ensuring that cells are initialized to swim at maximal speed. The projecting system has an optical resolution of 2 μm approximately matching the size of a single cell which represents the physical “pixel” of our density images. This value sets the limit for the minimum theoretical resolution of density configurations that would be achievable if bacteria could be precisely and statically arranged in space. In practice, as we will see, the real resolution will always be larger for two main reasons: (i) bacteria do not respond instantly to light temporal variations thus introducing a blur in the target speed map, (ii) the stationary state is an ensemble of noisy patterns that constantly fluctuate because of swimming and Brownian motions of bacteria.”
See also reply to point 4 for a deeper theoretical analysis of how memory affects resolution and how feedback can help alleviating this blurring effect.
3) Presentation: Many concepts and methods warrant expanded explanation. Symbols are introduced without definition. There is inadequate description of experimental methods, a diagram of the experimental protocol would have been helpful, highlighting parameters such as size and geometry of chamber etc. The method of preparation of the bacterial suspension is also unclear. (In fact, supplementary figure 1 might be helpful as part of the main text.) Several algorithms for data extraction or image manipulation are glossed over and require careful definition.
We thank the reviewers for addressing these issues. We fully agree with the above comments and have extensively expanded the Materials and methods section to describe the details of data extraction and image analysis.
The Materials and methods section now includes a detailed subsection about DDM and analysis.
We have also added a detailed explanation of the normalization procedure, the computation of the distances from the target and the feedback loop (subsection “Image analysis” and subsection “Feedback”).
We also moved the supplementary figure to the main text as Figure 1 after adding schematic and geometric details of sample preparation.
4a) Interpretation: Discussion of the theoretical background can be improved. For instance, the assumed inverse proportionality law between local density and speed is itself a simplification, and the limitations of this simple theory could be discussed in the context of the experimental results. (The discussion of Figure 1B is ambiguous.)
We thank the reviewers for pushing us to go further in this direction. Stimulated by their questions we arrived at an original theoretical result that provides an expression for the stationary density of smooth swimming bacteria exploring a nonhomogeneous light pattern and with a speed response to temporal light changes that is not instantaneous. The details of the theory are reported in the Materials and methods section. The revised manuscript now includes a discussion of this theoretical result in the context of experimental data:
“As a result, bacteria will experience an effective speed map that is a blurred version of what we would expect for an instantaneous response V (r) = v(I(r)). […]
Our aim here, however, is to investigate the resolution limits of density shaping with photokinetic bacteria.”
4b) There is a growing body of work in which active matter is manipulated successfully with light, and these should be cited and discussed. Notably, recent papers from the RiedelKruse (Jin and RiedelKruse, 2018, PNAS), or the Poon group (Arlt et al., 2018, Nature Comms.), and several others besides. The authors should ensure previously published literature is clearly referenced, and importantly, how their work differs from or represents improvements upon these previous results.
The paper by the Poon group was already referenced (Arlt et al., 2018). We now better highlight the main differences from our work as follows:
“[1] have recently shown that, by projecting a masked illumination pattern, “photokinetic” [11] bacteria can be accumulated in dark regions and depleted from brighter ones thus forming binary patterns. […] Finally we show that a model independent feedback control loop allows density shaping with high spatial resolution and gray level accuracy.”
We thank the reviewers for pointing us to the paper by the ReidelKruse lab which we have included in the bibliography and discussed in the concluding section:
“An alternative strategy to produce static patterns of bacteria is by “biofilm lithography” where an optical template of blue light is used to induce the expression of membrane proteins that promote cellcell and cellsubstrate attachment [4]. As opposed to our motility driven density shaping, these patterns are static and form over a time scale of several hours. Interestingly, the two techniques could be used in combination to achieve faster and more complex (e.g. three dimensional) biofilm lithography by fixing with blue light a template pattern obtained by speed modulation with green light.”
5) Purpose: While this is a beautiful result, the wider significance of this work is unclear. We are unconvinced by the bioengineering argument that these active materials are somehow better/more controllable for micro optics applications; it is argued in terms of drift and diffusion timescales that active systems are much more responsive than their Brownian counterparts, this is undeniable – but why might it be desirable to control such samples in the first place? Furthermore, are these patterns stable? If continued to be subjected to the light, individual bacteria will adapt to the surrounding light conditions, how long will a typical Mona Lisa image remain recognisable before the pattern degrades?
We have expanded the concluding section including novel application ideas:
“In conclusion, we have shown that a suspension of swimming bacteria, with optically controllable speed, can provide a new class of light controllable active materials whose density can be accurately, reversibly and quickly shaped by employing a low power light projector. […] Finally, the possibility of spatial and temporal control of motile bacteria density could also lead to novel strategies for the transport and manipulation of small cargoes inside microdevices [5].”
Regarding stability, as we now say in the main text:
“The image is stable and resolution does not deteriorate for hours. However a dark region surrounds the illuminated area so that swimming cells that exit the field of view do not come back and the number of motile cells in the illuminated area is reduced by about 50% in one hour. This issue is the main limiting factor for the lifetime of our patterns and could be solved in the future by using sample chambers that match the size of the illuminated area.”
We would have liked to have seen either (i) a clear theoretical interpretation in terms of the theory of motility induced phase separation, e.g. what are the physical reasons for the measured deviations from the simple inverse densityspeed relationship, or (ii) a clear biological interpretation of the timescales/transients involved in the pattern formation process in terms of some intracellular "memory". More interesting would be if the authors could measure, by combining precisely projected patterns with feedback modulation of population density, the actual spatiotemporal memory of the bacteria and how this might change with genetic manipulation.
We decided to move in the second direction and investigate the role of intracellular “memory” with new experiments and an original theoretical analysis. We have used DDM to accurately evaluate the time response of swimming speed to a square wave light intensity modulation. Obtained data are now reported in a new figure and discussed in the main text as follows:
“In particular, the theory assumes that speed is a local function of space which would only be the case if bacteria instantly adapt to temporal changes in light intensity. […] A slower relaxation follows and it is well fitted by an exponential function with a time constant τ_{m} = 35 s that is the same for both the rising and falling relaxations.”
Looking for a possible origin of the slow component we have performed a more detailed genetic screening of our strain. In particular, contrary to our original expectations, we confirmed with PCR the presence of the uncoperon coding for the ATPsynthase complex. The Materials and methods section has been modified accordingly.
https://doi.org/10.7554/eLife.36608.018Article and author information
Author details
Funding
European Research Council (307940)
 Roberto Di Leonardo
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP7/20072013)/ERC grant agreement no. 307940.
Senior Editor
 Naama Barkai, Weizmann Institute of Science, Israel
Reviewing Editor
 Raymond E Goldstein, University of Cambridge, United Kingdom
Publication history
 Received: March 13, 2018
 Accepted: July 12, 2018
 Version of Record published: August 14, 2018 (version 1)
Copyright
© 2018, Frangipane et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 12,431
 Page views

 667
 Downloads

 55
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Further reading

 Physics of Living Systems
Little is known about how muscle length affects residual force enhancement (rFE) in humans. We therefore investigated rFE at short, long, and very long muscle lengths within the human quadriceps and patellar tendon (PT) using conventional dynamometry with motion capture (rFE_{TQ}) and a new, noninvasive shearwave tensiometry technique (rFE_{WS}). Eleven healthy male participants performed submaximal (50% max.) EMGmatched fixedend reference and stretchhold contractions across these muscle lengths while muscle fascicle length changes of the vastus lateralis (VL) were captured using Bmode ultrasound. We found significant rFE_{TQ} at long (7±5%) and very long (12±8%), but not short (2±5%) muscle lengths, whereas rFE_{WS} was only significant at the very long (38±27%), but not short (8±12%) or long (6±10%) muscle lengths. We also found significant relationships between VL fascicle length and rFE_{TQ} (r=0.63, p=.001) and rFE_{WS }(r=0.52, p=.017), but relationships were not significant between VL fascicle stretch amplitude and rFE_{TQ} (r=0.33, p=.126) or rFE_{WS }(r=0.29, p=.201). PT shearwave speedangle relationships did not agree with estimated PT forceangle relationships, which indicates that estimating PT loads from shearwave tensiometry might be inaccurate. We conclude that increasing muscle length rather than stretch amplitude contributes more to rFE during submaximal voluntary contractions of the human quadriceps.

 Cell Biology
 Physics of Living Systems
In addition to diffusive signals, cells in tissue also communicate via long, thin cellular protrusions, such as airinemes in zebrafish. Before establishing communication, cellular protrusions must find their target cell. Here, we demonstrate that the shapes of airinemes in zebrafish are consistent with a finite persistent random walk model. The probability of contacting the target cell is maximized for a balance between ballistic search (straight) and diffusive search (highly curved, random). We find that the curvature of airinemes in zebrafish, extracted from livecell microscopy, is approximately the same value as the optimum in the simple persistent random walk model. We also explore the ability of the target cell to infer direction of the airineme’s source, finding that there is a theoretical tradeoff between search optimality and directional information. This provides a framework to characterize the shape, and performance objectives, of noncanonical cellular protrusions in general.