Environment-dependent swimming strategy of under magnetic field

Nicolas Waisbord,1 Christopher T. Lef`evre,2 Lyd´ericBocquet,3 Christophe Ybert,1 and C´ecileCottin-Bizonne1 1Institut Lumi`ere Mati`ere, CNRS UMR 5306, Universit´eClaude Bernard Lyon1, Universit´ede Lyon, France 2CNRS/CEA/Aix-Marseille University, Biosciences and Biotechnologies Institute, 13108 Saint Paul lez Durance, France 3Laboratoire de Physique Statistique, Ecole Normale Suprieure, UMR CNRS 8535, 75005 Paris, France Magnetotactic (MTB) are fascinating micro-organisms which possess embodied biomin- eralized nanomagnets providing them the ability to orient with the Earth’s magnetic field. This property is presumably related to an evolutionary advantage in finding the oxic-anoxic interface along the up and down direction in aquatic environments. So far the magnetic field response by MTB, called magnetotaxis, has been well described by a paramagnetic model where bacteria orient passively along the field lines according to a purely physical mechanism where magnetic torque and orientational Brownian noise compete. Here we demonstrate using Magnetococcus marinus strain MC-1 as MTB model that magnetotaxis shows more complex behaviors, which are affected by envi- ronmental conditions of different types. Indeed while MC-1 swimmers are found to essentially obey the paramagnetic paradigm when swimming in their growth medium, they exhibit a run-and-tumble dynamics in a medium devoid of energy source. Tumbling events are found to provide isotropic reori- entation capabilities causing the cells to escape from their prescribed field direction. This behavior has a major influence on the capabilities of the cells to explore their environment across field lines and represents an alternative search strategy to the back-and-forth motion along field-imposed tracks. Moreover, we show that aside chemical conditions, steric/geometrical constraints are also able to trigger tumbling events through obstacle encountering. Overall, physico-chemical environmental conditions appear to be important parameters involved in the swimming properties of MTB. De- pending on environmental conditions, the run-and-tumble mobility may provide advantages in the search for nutrient or ecological niche, in complement to classical magnetotaxis.

Many micro-organisms have the ability to move in the magnetic north (or south) for north-seekers in the their environment in order to respond to their needs, Northern Hemisphere (respectively south-seekers) [2, 7]. looking for instance for nutrients or for an optimum of Since the early studies on the motility of MTB, the re- oxygen concentration [1]. Among these micro-organisms sponse to magnetic field alone has been well described magneto-aerotactic bacteria (MTB) are aquatic prokary- by a paramagnetic model. In this framework, the bac- otes that show the specificity of a magnetically-assisted terium direction results from a competition between (i) aerotaxis: they passively orient with magnetic field lines a passive orientation of the bacterium along the mag- along which they actively swim, the so-called magneto- netic field lines associated to the torque exerted by the aerotaxis [2, 3]. Combined with the fact that –except field on the magnetic moment of the bacterium and (ii) at the equator– magnetic field lines are indicative of the thermal energy that tends to disorient the bacterium [8]. top to bottom location, this constitutes the basis of the Within this description, MTB have no active control on magneto-aerotaxis paradigm. Indeed, by switching their their swimming orientation apart from back and forth motions from 3-dimensional to 1-dimensional (1D), MTB switches. find more efficiently their ecological niche, that is anaero- bic or micro-aerobic conditions [4]. Thus, in the Northern Magneto-aerotaxis could be seen as a limitation in the Hemisphere, the local magnetic field is antiparallel to the ability of MTB for optimizing their chemical environment oxygen gradient and MTB are north-seeking while in the and thus one may wonder if MTB can override this para- Southern Hemisphere, the local magnetic field is parallel magnetic swimming behavior to let emerge other swim- to the oxygen gradient and MTB are south-seeking. The ming strategies [9]. This interrogation is legitimated by 1D motility of the cells is regulated back and forth to the fact that i) MTB are found around the equator [10] arXiv:1603.00490v1 [cond-mat.soft] 1 Mar 2016 remain in the vicinity of the oxic-anoxic interface. where the field lines no longer provide a relevant direction for exploring the oxygen gradient and ii) south-seeking Magneto-aerotaxis is a common feature shared by a MTB are sometimes found in the Northen Hemisphere large number of bacterial species disseminated all around [11] suggesting that alternate search mechanisms may be the globe in almost each and every freshwater and marine possible. In this contribution we explore the behavior environment [5]. MTB have in common the biomineral- of strain MC-1 under various environmental and mag- ization of magnetosomes, a prokaryotic organelle, gen- netic field conditions. Starting with a canonical exper- erally aligned in the cell and responsible for their mag- iment of chemical waves and bacterial aerotactic band, netic orientation [3]. Our study is focused on the model we first show that the field-mediated 1D motion strategy magnetotactic bacterium Magnetococcus marinus strain encompassed in the classical paramagnetic model cannot MC-1[6]. This bacterium has a polar magneto-aerotaxis: always account for the observed response of MC-1. Turn- in a homogenous chemical environment it swims towards ing to model situations of dilute bacteria in homogeneous 2

lar [7]). Unlike classical experiments however, where the magnetic field is pointing in the opposite direction of the A B 6 5 oxygen gradients, we explored different field configura-

B B 4 tions. In the present experiments, MC-1 cells are intro-

3 I I duced in a glass capillary filled with medium-rich solution 2 (see Materials and Methods) that is subsequently sealed. 0s 0.5 mm 30s 60s Speed ( μ m/s) Speed 1 Prior to introduction in the capillary, North seeking cells

0 0 0.2 0.4 0.6 0.8 1 were selected using a magnet. Thus placing the capillary 90° C B (mT) axis along the converging field direction obtained from Helmholtz coils leads to the formation of a high concen-

B V V tration of bacteria ”droplet” at the center of the capillary I I V V (see figure 1 A). Few minutes after, a band of bacte-

V V V ria that spontaneously propagates at a velocity around 3µm/s is formed. Due to the high concentration of bac- teria in the initial droplet and in the aerotactic band, chemical gradients develop and give rise to a propaga- 300 microns tion of the band (figure 1 A) as observed and described for other systems [2, 12, 13]. FIG. 1. A: Set up of the chemical wave experiment. Up: ini- Once the band is formed, it is possible to rotate the tial stage, the glass capillary aligned on the axis of a Helmoltz Coils set-up and MC-1 bacteria are concentrated in a central capillary axis and to use the Helmholtz coils to see how “droplet” using focusing magnetic fields from each Coil. After the band propagation responds to a uniform magnetic a few minute, a front grows from the droplet and a propagat- field. Most strikingly, we have measured the band veloc- ing chemical wave develops along the capillary. Down: 1D ity for magnetic fields perpendicular to the propagation propagation stage, the magnetic field is set homogeneous by direction (the capillary axis). As shown in figure 1 B, reversing one of the coils’ current and the capillary is either the remarkable observation is that the front propagation oriented perpendicular or along (not shown) the field direc- (i) persists and (ii) exhibits a velocity which is indepen- tion. Remarkably, bacteria band keep moving driven by the dent of the field magnitude. This observation is different chemical gradient even when the magnetic field is perpendic- from what was previously observed for Magnetospirillum ular to the propagation direction. B: Speed of the chemical gryphiswaldense strain MSR-1 [13] where the tilt of 90◦ wave band with the magnetic field intensity for: perpendicu- of the magnetic field drastically reduced the magneto- lar field (blue ◦), and parallel field towards propagation (green ♦) or against it (red ×). C: Example of tumbling tracks in the aerotaxis efficiency. This difference likely comes from the propagating front region, showing distinct tumbling events. different magnetotactic behaviors displayed by these two strains. Indeed, strain MSR-1 has an axial magnetotac- tic behavior where aerotaxis is prioritized over magne- chemical environments, we evidence the existence of a totaxis, while strain MC-1 is a dipolar MTB where the swimming strategy alternative to the paramagnetic re- cells give priority to following the direction indicated by sponse. This strategy takes the form of a run-and-tumble the magnetic field. swimming pattern by analogy to the well-known behav- Very clearly, such a motion associated to magneto- ior exemplified by Escherichia coli and more generally chemotaxis cannot fit into the existing framework of by a large variety of chemotactic bacteria. This strat- chemically-activated back and forth motion along the egy, which can override the magnetically-set orientation, magnetic field lines. Note that few experiments have also is found to be promoted by the chemical composition of been conducted where the field (i) was set homogeneous the surrounding environment but also by the geometry and parallel to the band propagation direction, or (ii) of the surrounding environment through collisions with was stopped. In all cases, band propagation persisted obstacles. with velocities not much affected. For instance figure Overall, our results provide a clear evidence that some 1 B shows the propagation speed along or against the MTB do not rely on the sole magnetically-assisted aero- magnetic field direction. While the field magnitude in- taxis while swimming toward the North or South but are deed either accelerates or decelerates the band, the field- also capable of supplementary strategies that deflect from assited contribution never exceeds the propagation veloc- their prescribed field direction. ity found for perpendicular field. Somehow, such an observation is coherent with the fact that magnetotactic cocci related to strain MC-1 can be I. CHEMICAL WAVES UNDER MAGNETIC found in most aquatic environments on Earth, includ- FIELD ing at the equator, where the magnetic field is horizontal and thus perpendicular to the oxygen concentration gra- We performed the benchmark aerotactic band ex- dient in water [10]. However, it requires extending the periment commonly used to characterize the magneto- magneto-aerotactic paradigm to incorporate alternative aerotactic behavior of MTB (axial, (di)polar, or unipo- swimming strategies which could overcome the constraint 3

pears closer to environmental conditions. Under these A B C 80 experimental configurations, single individuals dynamics could be followed under magnetic field, with no evolution 60 of the properties over 20 min.

40 Φ (%)

20 A. No field 250 microns 0 250 microns

D rich media E poor media 6 We start by examining the swimming behavior in the 1

) 5 absence of –added– magnetic field. Note that the earth )

0.8 t

4 (

f magnetic field at our latitude is around 17µT, for which 0.6 d 3

p our cultured strains showed no biased directional motion.

0.4 ( 2

g Figure 2 A shows a typical set of individuals tracks cap-

o 1 0.2 l tured for MC-1 in the rich medium. Each track has been θ (t) - (t+ τ )) >

with increasing directionality as the magnetic field in- A B creases. Under the classical description, we expect the B bacteria orientations to be passively set by the magnetic torque exerted on magnetosomes, with the finite width

PDF of the distribution arising from the orientation noise at- tributed to rotational Brownian motion. Indeed the fig-

200 microns ure 3 B shows that the PDF of smooth swimmers is per- fectly fitted by a Langevin paramagnetic prescription of θ (rad) the form PDF∝ exp([MB/kT ] cos θ). Moreover, the fit- C D 1 ted distribution yields the value of the reduced energy 0.9 MB/kT = 8 ± 1. Assuming that the only source of ori- 0.8 0.7 entational noise comes from thermal Brownian motion, 0.6

0.5 this provides a value of the magnetic momentum carried PDF os(θ)> 0.4 −16 < c by each bacterium of M ≈ 1.7 ± 0.2 10 S.I., consis- 0.3 0.2 tent with direct measurements of the momentum (see 0.1 Materials and Methods). Accordingly, the average axial 0 0 0.2 0.4 0.6 0.8 1 1.2 Speed (μm/s) B (mT) velocity (along the field direction) is perfectly fitted by the prediction for paramagnetic orientating swimmers as FIG. 3. A: Tracks under magnetic field magnetic field can be seen in figure 3 D. Note that the magnetic field (poor medium). B: PDF p(θ) of the swimmers orientation dependency was also checked and found to agree quanti- under field (B = 0.21 mT) for smooth (blue ◦) and tumbling tatively with the paramagnetic approach [3]. Therefore, (red ) swimmers. Inset: polar diagram for smooth swim- the smooth swimmers, which dominate the observed dy- mers fitted by paramagnetic Langevin model. The fit for namics in rich medium, are corresponding to the previ- the smooth swimmers by p(θ) ∝ exp(MBcos(θ)/kT ) gives ously reported response of MTB under field. It further M = 1.7 × 10−16A.m2 C: Same as in B, for the PDF of veloc- ensures the coherence of our observations with previous ity magnitude; dashed lines corresponds to the mean veloci- ties (respectively 91±25 µm/s and 104 µm/s). D: Normalized reports [20]. averaged axial velocity as a function of the magnetic field for Concerning the run-and-tumble dynamics, we first smooth (blue ◦) and tumbling (red ) swimmers. The blue compare the velocity PDFs of smooth and tumbling line is the fit of this orientation of the smooth swimmers by swimmers in figure 3 C. As can be seen, the swim- the Langevin function coth MB  − kT . kT MB ming activity appears unmodified in tumbling condi- tions, with a mean velocity essentially unchanged around 100 ± 25 µm/s and marginal changes in the distribution. tumbling event. As shown in figure 2 E, the measured Indeed, the tumbling PDF differs only by a little skew density probability function is found to be exponential, towards low velocities which is coherent with the fact with a characteristic decay time of 0.6±0.2 s in agreement that tumbling events are associated to locally high rota- with the orientation decorrelation time. This is reminis- tional velocities and low translational ones. An “active cent of the characteristic Poisson distribution of the run magnetotaxis” model was proposed for Magnetospiril- times as found for E.coli [16], and further legitimates the lum magneticum strain AMB-1, where the magnetotactic suggested analogy with run-and-tumble trajectories. In cells would sense the magnetic field thanks to a chemo- some aspects, it echoes recent work on another magne- taxis protein that would interact with the magnetosome totactic bacterium, Magnetospirillum magneticum strain chain and transfer the magnetic signal to the flagella [21]. AMB-1, which showed that the response to the magnetic We have measured the tumbling rate as a function of the field of MTB can incorporate active pathways [19]. Here magnetic field and have seen no dependency which leaves we do also identify an active response into the magneto- the question of the challenging model of “active magento- aerotactic bacteria dynamics; however not related to the taxie” not consistent with our observation in strain MC- magnetic field action but to the chemical environment. 1. The ability of bacteria to escape the direction set by the external magnetic field can be quantified through the B. Under homogeneous magnetic field diffusion coefficient D⊥ in the direction perpendicular to the imposed field. In practice, this is measured from the Under magnetic field, as shown in figure 3 A, cell’s correlation function according to the Green-Kubo rela- R ∞ trajectories continue to exhibit two distinct behaviors: tionship: D⊥ = 0 hv⊥(t)v⊥(0)i dt. Figure 4 shows the smooth and run-and-tumble. cross-diffusivity for smooth and tumbling swimmers as a For the smooth swimmers, in figure 3 B we see the po- function of the external field, the two appearing markedly lar orientation diagram (inset) together with the proba- different. Under low fields, smooth swimmers explore bility distribution function (PDF) of orientations. Start- unexpectedly more space around the magnetic field line ing from an isotropic distribution without field, the swim- than tumblers. Indeed the bacteria have an effective dif- 2 ming direction gets peaked around the field direction, fusion coefficient that is proportional to V τr where τr 5

A B C constraints can induce changes in their swimming behav- 5 10 1600 ior under the form of tumbling events. These constraints 2 7000

2 6000 are of steric or geometrical type as can be encountered

4

2 10 1200

D /D (B=0)

5000

B (mT) in natural environment such as sediments.

m /s) m (μ 4000 800 3

D 10 3000 Experimentally, these effects have been probed in a 400

2000 model geometry made up from a straight microfluidic

τ) m /s ) /s m (μ > 2

45

6 40 5

35 4

3 30

2

25 L.B ( μ m.mT) 1

20 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 B (mT)

L ( μ m) 15

10

5

0 B C 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 B (mT)

FIG. 5. A: Force-free rheological experiment: bacteria in structured micro-channels. B: Bacteria tumbles upon wall-collision, followed by trajectories during orientation relaxation at 140 µT. C: Same as B for a field of 490 µT. D: Velocity Profile averaged alongside the channel (from tooth extremities up to flat saturated velocity); Solid line correspond to an exponential profile A exp(−x/L), thus yielding the wall influence length L. E Velocity profiles for various magnetic fields. F: Variations of L under various magnetic field B. Inset: LB vs B.

regime and so is the influence length L ∝ V0τR. swim in one direction due to the magnetic forces they In the high field limit, the relaxation time towards the are still able to sense their environment and to swim field direction becomes shorter than the Brownian persis- in other directions if the chemotactic machinery of the tent time and thus dominates the extent to which wall in- cells detects the need to do so. In a rich medium bac- fluence propagates. The ratio between both times writes teria exhibit a smooth motility perfectly described by a Langevin paramagnetic model. In a poor medium, closer τ 8πηR3 8πηR3 −1 M to environmental conditions, we show that MC-1 cells can R = = × B. (1) τB kB.T M.B kB.T exhibit a run-and-tumble motion outstripping the sim- ple paramagnetic behavior. By changing their direction According to the orientation distribution (Fig. 3), this of motility in poor media bacteria increase the chance ratio is around 4 for B = 0.1 mT, for which we do ob- to meet their energy source while swimming around the serve the plateau in LB (figure 5 F). This confirms the oxic-anoxic interface. Furthermore we show that that domination of the magnetic alignment relaxation time tumbling can also occur under geometrical constrains as in this regime, for which we then expect the L ∝ 1/B encountered in crowded environments of sediments for dependency according to instance. Finally, in a wider perspective, let us note that when 3 8πηR under the conditions for a paramagnetic swimming re- L ∝ V0τB = V0 . (2) M.B sponse, magnetotactic bacteria form a remarkable ac- tive system of extremely efficient self-propeled individu- Overall, flow profiles close to sawtooth-shaped walls als whose orientation can be quantitatively controlled by suggest that the effect of an obstacle on bacteria dynam- a physical magnetic field [26]. In the blooming context ics can be summed up by a wall tumbling followed by of biological active matter or in the perspective of driv- an orientation relaxation over a length scale either domi- able in-body micro-swimmers for medical applications, nated by a persistent random walk exploration (low field) this system forms a very appealing model of controllable or a deterministic trajectory under magnetic field. In the active system. environment, obstacles are many and the probability to meet a particle is high particularly in sediment where MC-1 and other MTB are generally found. Such tum- bling behavior allows the cells to avoid obstacles and to V. MATERIALS continue swimming toward their preferred oxygen con- centration. Bacterial strain. Cells of Magnetococcus marinus strain MC-1 were grown in a slightly modified semi-solid medium described in [6].The medium consisted of an arti- IV. CONCLUSION. ficial seawater (ASW) containing (per litre): NaCl, 16.43 g; MgCl2.6H2O, 3.49 g; Na2SO4, 2.74 g; KCl, 0.465 Our results show that the motility of strain MC-1 is af- g; and CaCl2.2H2O, 0.386 g. To this was added (per fected by the environment. The magneto-aerotaxis model litre) the following prior to autoclaving: 5 mL modified is not just a passive alignment of the cells while they are Wolfe’s mineral elixir [2], 0.25 g NH4Cl, 2.4 g HEPES, swimming. Indeed, whereas the cells are constrained to 100 µL0.2% (w/v) aqueous resazurin and 2.0 g agar noble 7

(Difco). The medium was then adjusted to pH 7.0, boiled the screen on the acquisition system. to dissolve the agar and autoclaved. After the medium Preparation of the experiment. We used linear had cooled to about 45C, the following solutions were silicone polydimethylsiloxane (PDMS) micro-channels added (per liter) from stock solutions (except for the cys- (1200µm × 70µm × 5mm). The channels were filled teine, which was made fresh and filter-sterilized directly with the swimming medium by capillary suction, then into the medium): 2.8 mL 0.5 M potassium phosphate sealed on one side. A drop of a culture of strain MC-1 buffer, pH 6.9; L-cysteine, to give a final concentration was disposed at the extremity of the channel and steered of 0.4 g.L−1; 0.5 mL vitamin solution [2]; 3 mL 0.01 M inside the micro-channel with a simple stirring magnet. FeSO dissolved in 0.2 M HCl; 3 mL 40% sodium thiosul- 4 Once the bacteria inside the channel, this latter is sealed fate pentahydrate solution (i.e., 10mM of energy source); on the former entrance. Using a computer controlled and 2.7 mL 0.8 M NaHCO (autoclaved dry; sterile water 3 Helmholtz coils system we control the magnetic field in- added after autoclaving to make the fresh stock solution). side the channel. We concentrate the north seeking bac- The medium (10 mL) was dispensed into sterile, 15 x 125 teria inside the channel by using a configuration where mm screw-capped test tubes. All cultures were incubated coils generate fields that face one another and then we at 28◦C and, after approximately 1 week, a microaero- set the magnetic field to zero. bic band of bacteria formed at the oxic-anoxic interface (pink-colorless interface) of the tubes. This band con- Microscopy and data acquisition. To acquire im- tains a majority of north-seeking cells; the rare south- ages we use a Leica DMI 4000B Microscope in transmit- seeking cells were eliminated via a magnetic separation ted light with a Zeiss 5X long working distance objective, that consists in transferring the cells from their semi-solid equipped with a Baumer HXC40 4Mpx CCD camera con- medium to a liquid “swimming medium”. trolled by custom Labview codes, covering a 1200µm× 1200 µm area. We run both rheology and single par- Swimming media. The medium referred as rich ticle experiments at 35 fps. The long-working distance medium consisted of ASW to which was added (per litre) objective allows a depth of field of ≈100 µm, making the following: 5 mL modified Wolfe’s mineral elixir, 0.25 possible to follow bacteria and acquire the projection of g NH Cl, 2.4 g HEPES. The medium referred as poor 4 their trajectories. In single particle experiment, we as- medium consisted of the same medium as the rich one sume that the problem is axially symmetric, which allows devoid of thiosulfate. Both media have a pH adjusted to several calculus tricks for the relationships between pro- 7. jected observables and 3D ones. Microfabrication. Microchannels were fabricated by Tracking and trajectory analysis. Image treat- standard “soft lithography” technique: a mold made ment and tracking was performed using standard track- of negative photoresist (SU8 3100 Microchem) is ob- ing routine written with Matlab software (Mathworks). tained by conventional photolithography. Channels are The selection of the trajectories was based on duration then made of PDMS poured on these molds, at 70 ◦C (¿1.5 s) exploration radius (¿ 50µm) and average instant for 3 hours, with a typical thickness of 1mm. After speed (¿ 40µm.s−1). To avoid taking into account the He- cross-linking, the PDMS imprint is pealed of its mold. lical Motion of the bacteria, we smoothed the positions The channels are cut, the entrance and the outlet are by a gaussian filter on a width of 4 positions, so that the punched, and are finally bonded with a glass slide after projected orientation of motion we consider in this study an air plasma treatment. The obtained channels are then is the one of the main component of the motion. The pro- filled with swimming media by capillarity. The outlet is jected observed speed is assumed to be the real one. The then sealed with capillary wax. A bacteria drop is then Transverse projected MSD is assumed to be half of the taken from its growing medium and then guided inside real one. The correlation function of the projected mo- the channels with a magnet. Once bacteria are inside, tion is assumed to be the same as the one in 3D. Finally, the entrance is sealed with the wax. the projection of the Orientation Distribution function Coils control. Coils are controlled directly from the I1(ξ.cos(α))+L−1(ξ.cos(α)) are fitted by the 4ξ−1 sinh(ξ) , which is the computer, with a Labjack U12 device and an homemade statistics of the orientation expected of the Langevin Dis- current source. They can be set in different configura- tribution, projected in the focal plane. Motivated by the tions, generating Magnetic field parallel or anti-parallel. large persistence lengths of the trajectories (over 200µm We use the anti-parallel configuration to concentrate bac- for smooth swimmers), we computed D⊥(B), the diffu- teria, and at the same time sort the south seekers from sion coefficient perpendicular to the magnetic field, from the north seekers. The parallel configuration allows fields the Green-Kubo formula, which reads: ranging from 0 to 2.8mT, considered homogeneous in the observation zone of 1.2 mm. Alignment of channels with Z ∞ the coils is made as follows: coils are designed with a 2D⊥(B) = hV⊥(τ)V⊥(τ + t)iτ dt. (3) 0 rectangular hole where a Nikon Mir Glass Slide fits ex- actly: the field and the line of the mir are perpendicular. We obtain D⊥(B) from the correlation functions of the The camera observation is then aligned on the mir. Then speed in the transverse directions, which necessitates to the micro channel is aligned with a ROI crop selection on follow bacteria during a shorter period of time and can be 8

fitted with exponential functions both for the tumblers ried out an experiment where non swimming bacteria are and the smooth swimmers. submitted to a gradient of magnetic field: living bacteria are suspended in a liquid rich medium with 0.5% BSA Tumble detection. Inspired from the method of [15], and buffered at pH 10. This high pH kills the bacteria we define a function which combines the variation of the while preserving the cell integrity and the BSA prevents speed of the bacteria and their rotational acceleration, them from sticking to the glass capillary in which the both characteristic of tumbling event where a bacteria suspension is sucked. The capillary is then disposed in stops and turns: α the middle of two facing magnets (Supermagnets W-03-N f (t) = |ω(t)| α,β |V (t)|β Gottmadingen, Germany) creating a magnetic gradient For a smooth swimmer, this quantity cannot increase of ∇~ B~ = 80T.m−1. Bacteria are dragged across the a lot more than 4 times its median value computed along a trajectory, as the disorientation is well described by a channel under the magnetic force created by such a gra- White Gaussian Noise. For a tumbler, the median reflects dient of field and the magnetic momentum is estimated its smooth swimming behavior as, for a clear majority of by measuring the average velocities of the bacteria: 6πηR.V its time it is indeed smooth swimming, and the tumbling M = drag . k∇~ B~ k events are not visible in the computation of the median. We measure a magnetic momentum M = 1.3 × These aredetected, and we compare it to 10 times the 10−16A.m2 with a standard error 0.5 × 10−16A.m2. This median of fα,β(t) for each trajectory. This threshold and value is obtained by an overage over 500 trajectories of the exponents (α, β) are adjusted and for our data, a bacteria coming from three different batches and it is threshold of 10 and (α, β) = (2, 2) are satisfactory. within the range of magnetic momentum obtained for PIV analysis. PIV analysis was performed using the other magnetotactic bacteria [20]. PIVlab OpenSource Matlab toolbox. Images had their background subtracted before PIV was run. The average ACKNOWLEDGMENTS velocity was done over 400 images, along 15 periods of the waved wall. We have also shown, in a force free rhoelogy We thank Eleonora Secchi for her help with the data experiment, that tumbling at the wall can be induced analysis. We thank Veroniquie Utzinger, Guy Condemine by roughness which modifies the shape of the velocity and Nicole Cotte-Pattat for their help with the bacterial profile. cultures. We thank INL for access to their clean room Paramagnetic behavior in rich medium. An im- facilities and AXA research fund for its financial support. portant physical parameter to estimate for a magnetic C.T.L. was supported by the French National Research bacterium is its magnetic moment. For this we have car- Agency (GROMA: ANR-14-CE35-0018-01)

[1] Wadhams GH, Armitage JP (2004) Making sense of it [10] Frankel RB, Blakemore RP, De Araujo F (1981) Mag- all: bacterial chemotaxis. Nat Rev Mol Cell Biol 5:1024– netotactic bacteria at the geomagnetic equator. Science 1037. 212:1269–1270. [2] Frankel RB, Bazylinski DA, Johnson MS, Taylor BL [11] Simmons SL, Bazylinski DA, Edwards KJ (2006) South- (1997) Magneto-aerotaxis in marine coccoid bacteria. seeking magnetotactic bacteria in the Northern Hemi- Biophysj 73:994–1000. sphere. Science 311:371–374. [3] Blakemore R (1975) Magnetotactic bacteria. Science [12] Saragosti J, et al. (2011) Directional persistence of 190:377–379. chemotactic bacteria in a traveling concentration wave. [4] Frankel RB, Bazylinski DA (2009) Magnetosomes and P Natl Acad Sci Usa 108:16235–16240. magneto-aerotaxis. Contrib Microbiol 16:182–193. [13] Bennet M, et al. (2014) Influence of Magnetic Fields on [5] Lefevre CT, Bazylinski DA (2013) Ecology, Diversity, Magneto-Aerotaxis. PLoS ONE 9:e101150. and Evolution of Magnetotactic Bacteria. Microbiology [14] Zhu X, et al. (2014) Angle sensing in magnetotaxis of and Molecular Biology Reviews 77:497–526. Magnetospirillum magneticum AMB-1. Integrative Biol- [6] Bazylinski DA, et al. (2013) Magnetococcus marinus gen. ogy 6:706–18. nov., sp. nov., a marine, magnetotactic bacterium that [15] Masson JB, Voisinne G, Wong-Ng J, Celani A, Vergassola represents a novel lineage (Magnetococcaceae fam. nov., M (2012) Noninvasive inference of the molecular chemo- Magnetococcales ord. nov.) at the base of the Alphapro- tactic response using bacterial trajectories. P Natl Acad teobacteria. Int J Syst Evol Microbiol 63:801–808. Sci Usa 109:1802–1807. [7] Lef`evreCT, et al. (2014) Diversity of Magneto-Aerotactic [16] Berg HC, Brown DA (1972) Chemotaxis in Escherichia Behaviors and Oxygen Sensing Mechanisms in Cultured coli analysed by three-dimensional tracking. Nature Magnetotactic Bacteria. Biophysj 107:527–538. 239:500–504. [8] Frankel RB (1984) Magnetic guidance of organisms. [17] (1993) Random Walks in Biology (Princeton University Annu Rev Biophys Bioeng 13:85–103. Press). [9] Simmons SL (2006) South-Seeking Magnetotactic Bacte- [18] Meldrum FC, Mann S, Heywood BR, Frankel RB, ria in the Northern Hemisphere. Science 311:371–374. Bazylinski DA (1993) Electron Microscopy Study of Mag- 9

netosomes in a Cultured Coccoid Magnetotactic Bac- Journal of Molecular Biology 400:309–322. terium. Proceedings of the Royal Society of London B: [22] Palacci J, Cottin-Bizonne C, Ybert C, Bocquet L (2010) Biological Sciences 251:231–236. Sedimentation and Effective Temperature of Active Col- [19] Zhu X, et al. (2014) Angle sensing in magnetotaxis loidal Suspensions. Phys Rev Lett 105:088304. of Magnetospirillum magneticum AMB-1. Integr Biol [23] Stocker R (2012) Marine Microbes See a Sea of Gradients. 6:706–8. Science 338:628–633. [20] Nadkarni R, Barkley S, Fradin C (2013) A Comparison [24] Reufer M, et al. (2014) Switching of Swimming Modes in of Methods to Measure the Magnetic Moment of Magne- Magnetospirillium gryphiswaldense. Biophysj 106:37–46. totactic Bacteria through Analysis of Their Trajectories [25] Zhang SD, et al. (2013) Swimming behaviour and mag- in External Magnetic Fields. PLoS ONE 8:e82064–12. netotaxis function of the marine bacterium strain MO-1. [21] Philippe N, Wu LF (2010) An MCP-Like Protein Inter- Environmental Microbiology Reports 6:14–20. acts with the MamK Cytoskeleton and Is Involved in [26] Bazylinski DA, Frankel RB (2004) Magnetosome forma- Magnetotaxis in Magnetospirillum magneticum AMB-1. tion in prokaryotes. Nat Rev Micro 2:217–230.