Physiological magnetic field strengths help magnetotactic bacteria navigate in simulated sediments

  1. Agnese Codutti
  2. Mohammad A Charsooghi
  3. Konrad Marx
  4. Elisa Cerdá-Doñate
  5. Omar Muñoz
  6. Paul Zaslansky
  7. Vitali Telezki
  8. Tom Robinson  Is a corresponding author
  9. Damien Faivre  Is a corresponding author
  10. Stefan Klumpp  Is a corresponding author
  1. Max Planck Institute of Colloids and Interfaces, Department Theory and Biosystems, Germany
  2. Max Planck Institute of Colloids and Interfaces, Department Biomaterials, Germany
  3. Physics Department, TU München, Germany
  4. University of Göttingen, Institute for the Dynamics of Complex Systems, Germany
  5. Charité - Universitätsmedizin Berlin, Department for Operative, Preventive and Pediatric Dentistry, Germany
  6. Institute for Bioengineering, School of Engineering, University of Edinburgh, United Kingdom
  7. Aix-Marseille Université, CEA, CNRS, BIAM, France
8 figures, 2 tables and 1 additional file

Figures

Figure 1 with 2 supplements
Construction of sediment-mimicking obstacle channels.

(a) Micro-computer tomography slices of sediment sampled (diameter 4mm) are cropped out to ignore border effect due to the reconstruction and binarized (b). The binarized images are then fitted with circles (c) for statistical analysis (grain size distribution in d) and for construction of arrays of pillar-shaped obstacles (e). The images of irregular grains (b) and fitted circles (e) are used as masks to design microfluidic channels with obstacle arrays (f), in which the bacteria are injected in the middle channel (3.5 mm × 1.2 mm × 10 μm) where they encounter the obstacles.

Figure 1—figure supplement 1
Two microfluidic channels with irregularly shaped (A, C) obstacles and the corresponding channels with rounded obstacles (B, D).

Note that in (C), a narrow channel has been cut out of the big obstacle in the center (indicated by the red square), as it would block the channel completely. The scale bar is 150 μm.

Figure 1—figure supplement 2
Masks for all obstacle arrays used in this study.

Left column: irregular shaped obstacles based on 2d μCT images of sand samples, right column: rounded obstacle.

Analysis of micro-computer tomography images of a sediment sample in water.

(a) Slice of the sediment sample, (b) cube from the center of the cuvette, (c) and (d) examples of water gaps rendered with the Amira software. (e) Distribution of the smallest grain dimension (i.e. trabecular thickness analyzed with CTAn) with Gaussian fit aexp((xμ)2/c2), μ = 3.83, c = 0.5903, corresponding to a mean sand grain diameter of 46 μm. (f) Distribution of the water gap size (i.e. trabecular separation) with Gaussian fit aexp((xμ)2/c2), μ = 42.94 μm, c = 28.66 μm. (g) Percentage of sand over water as function of the distance from the bottom of the cuvette as obtained from 2D slices of the μCT images (yellow filled line) and for the equivalent fit with circles (yellow dashed line). The red points show the mean radius of the fitted circles (± standard error of the mean) in μm (right axis).

Figure 3 with 3 supplements
Swimming of magnetotactic bacteria through obstacle channels.

(A) View of a channel, in which the entrance (IN) and exit regions (OUT) are indicated, in which the bacterial density was monitored. Bacteria enter the IN region through an inlet that connects to a syringe containing bacteria. The magnetic field points to the left. (B) Cumulative bacterial intensity (measuring the cumulative arrival of bacteria) in these two regions, for three different field strengths. (C) Bacterial throughput, quantified by the ratio of the cumulative intensities in the OUT and IN regions as function of the magnetic field. The reduction of throughput by B = 500 μT compared to B = 50 μT is significant (Mann-Whitney test, p<103), the increase between B = 0 and B = 50 μT is weakly significant (p=0.028), but agrees with the observation for a channel without obstacles in which the throughput increases monotonically with increasing field strength (Figure 3—figure supplement 3).

Figure 3—figure supplement 1
Sketch of the setup showing the microfluidic system (blue), the inlet connecting to a syringe with bacteria and growth medium, one of the three pairs of Helmholtz coils to control the magnetic field, and a detailed view of the microfluidic system.
Figure 3—figure supplement 2
Comparison of bacterial throughput as a function of the magnetic field (0, 50, 500 μT) in channels with irregular (sand-like) obstacles and in channels with rounded (pillar-like) obstacles.
Figure 3—figure supplement 3
Bacterial throughput in a channel without obstacles.

(A) Cumulative density of bacteria in the IN and OUT regions of the channel. (B) Ratio of the cumulative density in the IN and OUT regions, both as functions of time and for different magnetic field strengths. (C) Dependence of the latter ratio (taken after 30 min) on the magnetic field strength. In contrast to channels with obstacles (Figure 3), a monotonic dependence is observed.

Figure 4 with 2 supplements
Sliding of bacteria on obstacle surfaces: (a) Experimental trajectories of sliding bacteria (B=0).

(b) Trajectory of sliding particle in the simulation. The light blue region shows the spatial extent of the particle. Sliding is defined as motion in the sliding region (up to 2 μm from the surface, dotted line), provided that the particle also reaches the region where interactions with the obstacle take place (up to ≈ 0.6 μm from the surface, dashed line). The sliding distance is the distance covered tangentially to the surface. (c) Histogram of sliding distances as measured in the experiment (filled black circles) and the simulation (empty markers, for different values of the wall torque parameter α). We find good agreement for α = 0.2 μm.

Figure 4—figure supplement 1
Sliding distance as a function of the pillar radius as obtained from simulations for different wall torque parameters α.

The presented values show the median, error bars indicate 25th and 75th percentile. The sliding distances were extracted from obstacle channel simulation data, where pillar radii were binned into 20 μm bin sizes for statistical analysis.

Figure 4—figure supplement 2
Dependence of the sliding distance on the curvature of the pillar as obtained from the tracking bacteria in obstacle channels with cylindrical obstacles.
Simulated motion of bacteria in an obstacle channel: Heat maps of the time-averaged density of bacteria from simulations of 2000 bacteria for 30 min with a magnetic field B=50µT, pointing from left to right.

The blue map represents bacteria that did not arrive at the right end of the simulation box during the simulated time the red heat map represents those that did. The interplay between magnetic field and the channel’s geometry creates natural paths through the channel. These paths can include traps or end in traps that are formed by overlapping obstacles (examples are indicated by the green boxes).

Effect of magnetic field strength on simulated swimming through obstacle channels.

(a) Motion of bacteria in a channel for different values of the field strength B (heat maps of the densities of arriving and non-arriving bacteria, blue and red, respectively, as in Figure 5) (b) Fraction ϕ of simulated bacteria that have arrived up to time t for different B. Circles mark the first-arrival times tf; crosses with thin lines mark total arrival fractions ϕtot. The curve for B = 500 μT remains at 0. (c) First-arrival times tf (top) and total arrival fraction ϕtot (bottom) as a function of the field strength B. With stronger fields, tf converges to the first-arrival time of a persistent swimmer in an empty channel (dashed gray line). ϕtot shows a peak at an intermediate (optimal) field strength, which can be explained by the interplay of two opposing effects that arise with stronger fields: effectively faster motion in the direction of the field and higher susceptibility to trapping. Blue data points show the results for the channel shown in (a), gray lines show results in other channels.

Escape from traps - Simulation.

(a) Trajectories of 100 bacteria escaping from the apex of a symmetric trap (B = 50 μT). The trap consists of two obstacles of radius 50 μm, in a distance of Δy = 100 μm. This results in a trap depth d = 50 μm. (b) Distribution of escape times of bacteria from the trap in (a). An exponential curve exp(t/τ) (dashed line) is fitted to the data, with mean escape time τ = 22 s. (c) Dependence of the mean escape time τ on the field strength B for the trap in (a). An exponential curve, τ0;1exp(B/B¯), is fitted to the data (dashed line). (d) Trajectories of 100 bacteria escaping from the apex of a generic trap (B = 50 μT). (e) Dependence of τ on d for B = 50 μT. Blue points show symmetric traps with both obstacles of radius 50 μm, similar to the trap in (a), gray points show a representative sample of all possible traps. An exponential curve, τ0;2exp(d/d¯), is fitted to the blue data points (dashed line).

Escape from trap - Experiment.

The bacterial density in a trap (indicated by the yellow triangular area in the inset) is monitored over time. At the time points indicated by the arrows, the field strength is reduced (from 500 μT to 50 μT and from 50 μT to 0 μT, respectively), allowing the escape of trapped bacteria.

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Strain, strain background Magnetospirillum gryphiswaldenseMSR-1Schleifer et al., 1991MSR-1Obtained from Schüler Lab (Univ. Bayreuth, Germany)
Software, algorithmSimulation codeThis paperSee Methods, section Model and simulations
OtherMicrofluidic channelsThis paperSee Methods, section Microfluidic channel
Table 1
Model parameters and their values.
ParameterSymbolValue
Diameter of active particleσ1 μm (Klumpp and Faivre, 2016)
Magnetic momentμ4 × 10-16 J/T (Codutti et al., 2022; Nadkarni et al., 2013)
Self-propulsion speedv050 μm/s
Fluid temperatureT298 K
Fluid viscosityη8.9 × 10-4 Pa/s
Magnetic field strengthB[0 . . . 500] µT
Surface force parameterϵ4 kBT
Surface torque parameterα0.2 µm

Additional files

Download links

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

Downloads (link to download the article as PDF)

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

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

  1. Agnese Codutti
  2. Mohammad A Charsooghi
  3. Konrad Marx
  4. Elisa Cerdá-Doñate
  5. Omar Muñoz
  6. Paul Zaslansky
  7. Vitali Telezki
  8. Tom Robinson
  9. Damien Faivre
  10. Stefan Klumpp
(2025)
Physiological magnetic field strengths help magnetotactic bacteria navigate in simulated sediments
eLife 13:RP98001.
https://doi.org/10.7554/eLife.98001.3