pith. sign in

arxiv: 2511.09534 · v3 · submitted 2025-11-12 · ⚛️ physics.flu-dyn

Local chemotactic response of Escherichia coli in fluid and near surfaces

Adam Gargasson , Julien Bouvard , Carine Douarche , Peter Mergaert , Harold Auradou This is my paper

Pith reviewed 2026-05-17 22:05 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords chemotaxisEscherichia colimicrofluidic deviceconcentration gradientbacterial motilitysurfaceschemotactic susceptibilitydrift velocity
0
0 comments X

The pith

E. coli exhibits chemotactic drift velocity proportional to the log-gradient of concentration in fluid but with inhibited flux near surfaces.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes the local chemotactic response of Escherichia coli by using an optimized three-channel microfluidic device to impose a stable linear concentration profile of attractants such as casamino acids. In bulk fluid the measured drift velocity follows vc = χ(c) ∇c with the susceptibility χ(c) = χ0 / [(1 + c/c−)(1 + c/c+)], which reduces to vc ∝ ∇(log c) when concentrations lie between the two characteristic scales. The identical device shows that this directed flux is suppressed once bacteria interact with nearby surfaces. A reader would care because the result supplies a concrete, locally measurable rule for how bacteria bias their motion in free fluid versus confined geometries, directly linking single-cell trajectories to population-level navigation.

Core claim

In the fluid, the chemotactic response is described by the equation vc=χ(c) ∇c, with χ(c) = χ0 /[(1 + c/c−)(1 + c/c+)] the chemotactic susceptibility. For c− ≪ c ≪ c+, the bacterial chemotactic velocity is proportional to the concentration gradient divided by the concentration and vc ∝ ∇c/c = ∇(log c). However, on surfaces, the chemotactic flux is inhibited.

What carries the argument

The chemotactic susceptibility χ(c) that scales the drift velocity according to local concentration and its gradient, measured via stable linear profiles in a three-channel microfluidic device.

If this is right

  • Population-level chemotactic response can be extracted rapidly from many individual trajectories without needing bulk population assays.
  • In the effective concentration window, navigation is governed by the spatial derivative of log concentration rather than the absolute gradient.
  • Chemotactic accumulation and transport are strongly reduced once bacteria are within a body length of a surface.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Surface inhibition may limit the ability of bacteria to colonize or escape along walls in microfluidic channels or porous media.
  • Models of bacterial search strategies in heterogeneous environments could be simplified by replacing absolute-gradient sensing with log-gradient sensing in the bulk.
  • Repeating the same gradient protocol in channels of varying height would test whether the surface effect scales with distance to the wall.

Load-bearing premise

The three-channel microfluidic device produces a truly stable and linear concentration profile across the observation region without significant advection or mixing artifacts that would distort the measured drift velocities.

What would settle it

A direct plot of measured drift velocity versus ∇c/c across a wide range of mean concentrations should collapse onto a single curve only in fluid; any systematic deviation from the predicted proportionality or any restoration of flux on surfaces would falsify the central claim.

Figures

Figures reproduced from arXiv: 2511.09534 by Adam Gargasson, Carine Douarche, Harold Auradou, Julien Bouvard, Peter Mergaert.

Figure 1
Figure 1. Figure 1: Sketch of the experimental setup. (a) The microfluidic chip, placed on an inverted microscope and observed under 20× magnification, consists of three parallel channels. The channels ⃝1 and ⃝3 are filled with a chemo-attractant at concentration c1 and c3 diluted in the motility buffer (MB). (b) The chemo-attractant diffuses through the agarose layer, creating a linear gradient between channels ⃝1 and ⃝3 , n… view at source ↗
Figure 2
Figure 2. Figure 2: Bacterial tracks and concentration profiles for two experiments. [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Evolution of the average net, diffusive and [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Measuring the bacterial velocity bias in the transient state allows for a much quicker quantification of [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Transient and local measurements allow for a quicker and more accurate quantification of the log-sensing [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Absence of chemotactic drift on surfaces. [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Bacterial swimming velocity ⟨v¯s⟩T as a func￾tion of chemoattractant concentration. Purple and or￾ange markers correspond respectively to α-methyl-DL-aspartic (MeAsp) and casamino acids. The average is done over all the image sequences recorded during one experiment. Horizontal dashed line: average velocity for three experiments without chemoattractant (c1 = c3 = 0). The grey zone indicates the range of va… view at source ↗
Figure 8
Figure 8. Figure 8: Average correlation time ⟨τ ⟩T estimated from the velocity correlation function. Purple and orange markers correspond respectively to MeAsp and casamino acids. The average is done over all the image sequences recorded during one experiment. Horizontal dashed line: correlation time averaged over three experiments performed without casamino acids in the reservoir channels (c1 = c3 = 0 mM). Grey zone: the ran… view at source ↗
Figure 9
Figure 9. Figure 9: Normalised chemotactic velocity ⟨v˜c(T)⟩T , averaged over all times T, as function of ∇c measured for an average MeAsp concentration of c¯ = 100 µM. The solid line is a linear fit of the data by v˜c = χ∇c with χ = (0.72±0.05) µm/µM. The vertical bars represent the second moment of the distribution of v˜c(T). C Verification of the steady-state of casamino acids’ gradient Unlike MeAsp, casamino acids can be … view at source ↗
read the original abstract

Bacteria can adjust their swimming behaviour in response to chemical variations, a phenomenon known as chemotaxis. This process is characterised by a drift velocity that depends non-linearly on the concentration of chemical species and its "local" gradient. To study this process more effectively, we optimised a 3-channel microfluidic device to generate a stable, linear concentration profile of chemoattractants. This setup allows us to monitor the response of $Escherichia$ $coli$ to casamino acids or $\alpha$-methyl-DL-aspartic acid at the individual level. By analysing the movement of a population of individuals both in fluid and on surfaces, we achieve faster, more accurate quantification of the population's chemotactic response. In the fluid, the chemotactic response is described by the equation $v_c=\chi(c) \nabla c$, with $\chi(c) = \chi_0 /[(1 + c/c_-)(1 + c/c_+)]$ the chemotactic susceptibility. For $c_- \ll c \ll c_+$, i.e. when bacteria perform chemotaxis, the bacterial chemotactic velocity is proportional to the concentration gradient divided by the concentration and $v_c \propto \nabla c/c = \nabla (\log c)$. However, on surfaces, the chemotactic flux is inhibited.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The paper reports the use of an optimized three-channel microfluidic device to establish a stable linear concentration gradient of chemoattractants, enabling single-cell tracking of Escherichia coli chemotaxis both in bulk fluid and near surfaces. In fluid, the authors fit the observed drift velocities to vc = χ(c) ∇c with the rational susceptibility χ(c) = χ0 / [(1 + c/c−)(1 + c/c+)], which yields the logarithmic-gradient regime vc ∝ ∇c/c for intermediate concentrations. They additionally report that chemotactic flux is inhibited when cells interact with surfaces.

Significance. If the central empirical description holds, the work supplies quantitative single-cell data supporting a specific functional form for chemotactic susceptibility across concentration ranges and documents a clear difference in response between bulk fluid and surface-proximal conditions. This is relevant for modeling bacterial navigation in confined or heterogeneous environments. The individual-cell tracking approach in a controlled gradient is a methodological strength that enables direct measurement of local responses rather than relying solely on population averages.

major comments (2)
  1. [Microfluidic device and concentration profile] Microfluidic device and concentration profile section: The claim of an 'optimised' device producing a 'stable, linear' profile is central to determining the imposed ∇c used to fit χ0, c−, and c+. However, the manuscript provides no quantitative validation (e.g., fluorescence intensity profiles, repeated measurements over time, or CFD simulations) demonstrating that advection, incomplete mixing at junctions, or temporal drift remain negligible across the observation window. This directly affects the reliability of the extracted functional form and the comparison to surface data.
  2. [Fluid-phase chemotaxis results] Results on fluid-phase chemotaxis (around the definition of χ(c) and the vc ∝ ∇c/c regime): The fitted parameters are presented without reported uncertainties, number of tracked trajectories, or explicit goodness-of-fit statistics. Because the functional form is obtained by fitting rather than derived from first principles, the absence of these metrics makes it difficult to evaluate whether the rational expression is robustly supported or could be equally well described by alternative models.
minor comments (3)
  1. [Abstract] Abstract: does not report error bars, sample sizes (number of cells or trajectories), or statistical details supporting the stated functional form and surface-inhibition observation.
  2. [Introduction or Methods] Notation: the definitions and physical interpretation of the fitted constants c− and c+ are introduced without explicit reference to prior literature on receptor occupancy or adaptation models that might motivate the chosen rational form.
  3. [Figures] Figure presentation: several trajectory or velocity plots would benefit from overlaid model curves with confidence bands and explicit indication of the concentration range used for each fit.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments, which have helped us strengthen the presentation of our results. We respond to each major comment below and indicate the revisions made to the manuscript.

read point-by-point responses

  1. Referee: [Microfluidic device and concentration profile] Microfluidic device and concentration profile section: The claim of an 'optimised' device producing a 'stable, linear' profile is central to determining the imposed ∇c used to fit χ0, c−, and c+. However, the manuscript provides no quantitative validation (e.g., fluorescence intensity profiles, repeated measurements over time, or CFD simulations) demonstrating that advection, incomplete mixing at junctions, or temporal drift remain negligible across the observation window. This directly affects the reliability of the extracted functional form and the comparison to surface data.

    Authors: We acknowledge that the original manuscript did not include explicit quantitative validation of the concentration profile stability and linearity. In the revised version we have added fluorescence intensity measurements using a tracer dye, acquired at multiple time points across the observation window, together with CFD simulations of the device geometry. These data confirm that advection and mixing effects remain negligible and that the profile remains linear and temporally stable to within experimental precision. The added material directly supports the ∇c values used for the susceptibility fits and the fluid-versus-surface comparison. revision: yes

  2. Referee: [Fluid-phase chemotaxis results] Results on fluid-phase chemotaxis (around the definition of χ(c) and the vc ∝ ∇c/c regime): The fitted parameters are presented without reported uncertainties, number of tracked trajectories, or explicit goodness-of-fit statistics. Because the functional form is obtained by fitting rather than derived from first principles, the absence of these metrics makes it difficult to evaluate whether the rational expression is robustly supported or could be equally well described by alternative models.

    Authors: We agree that uncertainties, trajectory counts, and goodness-of-fit statistics are necessary to assess the robustness of the fitted rational susceptibility. The revised manuscript now reports the total number of tracked trajectories, the standard errors on the fitted parameters χ₀, c₋ and c₊ obtained from nonlinear least-squares regression, and the associated R² and reduced-χ² values. We have also added a brief comparison showing that the rational form yields a lower residual sum of squares than a simple linear or Hill-type alternative over the measured concentration range, thereby supporting the choice of functional form. revision: yes

Circularity Check

0 steps flagged

Empirical description of chemotactic response with no circular derivation

full rationale

The paper presents experimental measurements of E. coli trajectories in a microfluidic device and describes the observed chemotactic velocity via the phenomenological equation vc = χ(c) ∇c, where the rational form of χ(c) is fitted to data. This is a data-driven characterization rather than a first-principles derivation whose output reduces to its inputs by construction. No self-citations, uniqueness theorems, ansatzes, or renamings of known results are invoked as load-bearing steps in the provided claims. The central result remains an observational model self-contained against the experimental trajectories and concentration profiles.

Axiom & Free-Parameter Ledger

3 free parameters · 1 axioms · 0 invented entities

The central claim rests on three fitted parameters in the susceptibility function and the assumption that the microfluidic gradient is linear and stable; no new physical entities are postulated.

free parameters (3)
  • χ0
    Amplitude of chemotactic susceptibility, fitted to measured drift velocities.
  • c-
    Lower concentration scale below which susceptibility drops, fitted to data.
  • c+
    Upper concentration scale above which susceptibility drops, fitted to data.
axioms (1)
  • domain assumption The microfluidic device produces a stable, linear concentration profile across the observation region.
    Invoked to justify that measured drift velocities reflect true chemotactic response to a known ∇c.

pith-pipeline@v0.9.0 · 5552 in / 1481 out tokens · 38595 ms · 2026-05-17T22:05:09.484233+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

55 extracted references · 55 canonical work pages

  1. [1]

    H. C. Berg,Annual Review of Biophysics and Bio- engineering, 1975,4, 119–136

    work page 1975

  2. [2]

    Adler,Annual Review of Biochemistry, 1975,44, 341–356

    J. Adler,Annual Review of Biochemistry, 1975,44, 341–356

    work page 1975

  3. [3]

    Karmakar,Journal of Basic Microbiology, 2021, 61, 366–379

    R. Karmakar,Journal of Basic Microbiology, 2021, 61, 366–379

    work page 2021

  4. [4]

    Raina, B

    J.-B. Raina, B. S. Lambert, D. H. Parks, C. Rinke, N. Siboni, A. Bramucci, M. Ostrowski, B. Signal, A. Lutz, H. Mendis, F. Rubino, V. I. Fernandez, R. Stocker, P. Hugenholtz, G. W. Tyson and J. R. Seymour,Nature, 2022,605, 132–138

    work page 2022

  5. [5]

    Raina, V

    J.-B. Raina, V. Fernandez, B. Lambert, R. Stocker and J. R. Seymour,Nature Reviews Microbiology, 2019,17, 284–294

    work page 2019

  6. [6]

    H.-H. Tsai, Y. Tang, L. Jiang, X. Xu, V. D. Tendon, J. Pang, Y. Jia, K. Wippel, J. Vacheron, C. Keel, T. G. Andersen, N. Geldner and F. Zhou,Science, 2025,390, eadu4235

    work page 2025

  7. [7]

    Lextrait, S

    G. Lextrait, S. Joardar, R. Cossard, Y. Kikuchi, T. Ohbayashi and P. Mergaert,ISME Journal, 2025, in press,. Preprint — November 13, 2025 11

    work page 2025

  8. [8]

    P. M. Merritt, T. Danhorn and C. Fuqua,Journal of Bacteriology, 2007,189, 8005–8014

    work page 2007

  9. [9]

    Armitano, V

    J. Armitano, V. Méjean and C. Jourlin-Castelli,En- vironmental Microbiology, 2013,15, 3108–3118

    work page 2013

  10. [11]

    J. M. Keegstra, F. Carrara and R. Stocker,Nature Reviews Microbiology, 2022,20, 491–504

    work page 2022

  11. [12]

    R. M. Ford and R. W. Harvey,Advances in Water Resources, 2007,30, 1608–1617

    work page 2007

  12. [13]

    Bhushan, S

    B. Bhushan, S. K. Samanta, A. Chauhan, A. K. Chakraborti and R. K. Jain,Biochemical and Biophysical Research Communications, 2000,275, 129–133

    work page 2000

  13. [14]

    D. Paul, R. Singh and R. K. Jain,Environmental Microbiology, 2006,8, 1797–1804

    work page 2006

  14. [15]

    Harms and L

    H. Harms and L. Wick,Engineering in Life Sciences, 2006,6, 252–260

    work page 2006

  15. [16]

    Yamamoto-Tamura, I

    K. Yamamoto-Tamura, I. Kawagishi, N. Ogawa and T. Fujii,Bioscience, Biotechnology, and Biochemistry, 2015,79, 926–936

    work page 2015

  16. [17]

    Adler,Journal of General Microbiology, 1973,74, 77–91

    J. Adler,Journal of General Microbiology, 1973,74, 77–91

    work page 1973

  17. [18]

    Dahlquist, P

    F. Dahlquist, P. Lovely and D. Koshland,Nature New Biology, 1972,236, 120–123

    work page 1972

  18. [19]

    Holz and S

    M. Holz and S. Chen,Biophysical Journal, 1978,23, 15–31

    work page 1978

  19. [20]

    R.M.Ford, B.R.Phillips, J.A.QuinnandD.A.Lauf- fenburger,Biotechnology and Bioengineering, 1991, 37, 647–660

    work page 1991

  20. [21]

    Lewus and R

    P. Lewus and R. M. Ford,Biotechnology and Bioengi- neering, 2001,75, 292–304

    work page 2001

  21. [22]

    Morales-Soto, M

    N. Morales-Soto, M. E. Anyan, A. E. Mattingly, C. S. Madukoma, C. W. Harvey, M. Alber, É. Déziel, D. B. Kearns and J. D. Shrout,Journal of visualized exper- iments: JoVE, 2015, e52338

    work page 2015

  22. [23]

    Cremer, T

    J. Cremer, T. Honda, Y. Tang, J. Wong-Ng, M. Ver- gassola and T. Hwa,Nature, 2019,575, 658–663

    work page 2019

  23. [24]

    S. Kim, H. J. Kim and N. L. Jeon,Integrative Biology, 2010,2, 584

    work page 2010

  24. [25]

    Ahmed, T

    T. Ahmed, T. S. Shimizu and R. Stocker,Integrative Biology, 2010,2, 604

    work page 2010

  25. [26]

    J. Wu, X. Wu and F. Lin,Lab on a Chip, 2013,13, 2484

    work page 2013

  26. [27]

    Grognot and K

    M. Grognot and K. M. Taute,Communications Biol- ogy, 2021,4, 669

    work page 2021

  27. [28]

    J. Diao, L. Young, S. Kim, E. A. Fogarty, S. M. Heil- man, P. Zhou, M. L. Shuler, M. Wu and M. P. DeLisa, Lab on a Chip, 2006,6, 381–388

    work page 2006

  28. [29]

    Cheng, S

    S.-Y. Cheng, S. Heilman, M. Wasserman, S. Archer, M. L. Shuler and M. Wu,Lab on a Chip, 2007,7, 763

    work page 2007

  29. [30]

    Y. V. Kalinin, L. Jiang, Y. Tu and M. Wu,Biophysical Journal, 2009,96, 2439–2448

    work page 2009

  30. [31]

    Ahmed, T

    T. Ahmed, T. S. Shimizu and R. Stocker,Nano Let- ters, 2010,10, 3379–3385

    work page 2010

  31. [32]

    K. Son, F. Menolascina and R. Stocker,Proceedings of the National Academy of Sciences, 2016,113, 8624– 8629

    work page 2016

  32. [33]

    Menolascina, R

    F. Menolascina, R. Rusconi, V. I. Fernandez, S. Sm- riga, Z. Aminzare, E. D. Sontag and R. Stocker,npj Systems Biology and Applications, 2017,3, 16036

    work page 2017

  33. [34]

    Morse, R

    M. Morse, R. Colin, L. G. Wilson and J. X. Tang, Biophysical journal, 2016,110, 2076–2084

    work page 2016

  34. [35]

    Gargasson, C

    A. Gargasson, C. Douarche, P. Mergaert and H. Au- radou,Bio-protocol, 2025,15, e55184

    work page 2025

  35. [36]

    E. F. Keller and L. A. Segel,Journal of theoretical biology, 1971,30, 225–234

    work page 1971

  36. [37]

    M. A. Rivero, R. T. Tranquillo, H. M. Buettner and D. A. Lauffenburger,Chemical Engineering Science, 1989,44, 2881–2897

    work page 1989

  37. [38]

    L. A. Segel,SIAM Journal on Applied Mathematics, 1977,32, 653–665

    work page 1977

  38. [39]

    E. F. Keller and L. A. Segel,Journal of Theoretical Biology, 1971,30, 235–248

    work page 1971

  39. [40]

    G. L. Hazelbauer, J. J. Falke and J. S. Parkinson, Trends in Biochemical Sciences, 2008,33, 9–19

    work page 2008

  40. [41]

    Laganenka, R

    L. Laganenka, R. Colin and V. Sourjik,Nature Com- munications, 2016,7, 12984

    work page 2016

  41. [42]

    Y. Tu, T. S. Shimizu and H. C. Berg,Proceedings of the National Academy of Sciences, 2008,105, 14855–14860. 12 Preprint — November 13, 2025

    work page 2008

  42. [43]

    G. Si, T. Wu, Q. Ouyang and Y. Tu,Physical Review Letters, 2012,109, 048101

    work page 2012

  43. [44]

    Colin, R

    R. Colin, R. Zhang and L. G. Wilson,Journal of The Royal Society Interface, 2014,11, 20140486

    work page 2014

  44. [45]

    S. I. Bibikov, A. C. Miller, K. K. Gosink and J. S. Parkinson,Journal of Bacteriology, 2004,186, 3730–3737

    work page 2004

  45. [46]

    Stricker, I

    L. Stricker, I. Guido, T. Breithaupt, M. G. Mazza and J. Vollmer,Journal of The Royal Society Interface, 2020,17, 20200559

    work page 2020

  46. [47]

    Bouvard, C

    J. Bouvard, C. Douarche, P. Mergaert, H. Auradou and F. Moisy,Physical Review E, 2022,106, 034404

    work page 2022

  47. [48]

    Tinevez, N

    J.-Y. Tinevez, N. Perry, J. Schindelin, G. M. Hoopes, G. D. Reynolds, E. Laplantine, S. Y. Bednarek, S. L. Shorte and K. W. Eliceiri,Methods, 2017,115, 80–90

    work page 2017

  48. [49]

    P. S. Lovely and F. Dahlquist,Journal of theoretical biology, 1975,50, 477–496

    work page 1975

  49. [50]

    Taktikos, H

    J. Taktikos, H. Stark and V. Zaburdaev,PloS one, 2013,8, e81936

    work page 2013

  50. [51]

    Lauga,The Fluid Dynamics of Cell Motility, Cam- bridge University Press, 2020

    E. Lauga,The Fluid Dynamics of Cell Motility, Cam- bridge University Press, 2020

    work page 2020

  51. [52]

    Alert, A

    R. Alert, A. Martínez-Calvo and S. S. Datta,Physical Review Letters, 2022,128, 148101

    work page 2022

  52. [53]

    Lauga and T

    E. Lauga and T. R. Powers,Reports on progress in physics, 2009,72, 096601

    work page 2009

  53. [54]

    Junot, T

    G. Junot, T. Darnige, A. Lindner, V. A. Martinez, J. Arlt, A. Dawson, W. C. Poon, H. Auradou and E. Clément,Physical Review Letters, 2022,128, 248101

    work page 2022

  54. [55]

    Ganesh,Theses, Université Paris-Saclay, 2023

    A. Ganesh,Theses, Université Paris-Saclay, 2023

    work page 2023

  55. [56]

    Gargasson,Theses, Université Paris-Saclay, 2025

    A. Gargasson,Theses, Université Paris-Saclay, 2025. Preprint — November 13, 2025 13 Supplementary Material A Effect of the chemoattractant concentra- tion and gradient on the swimming velocity and motility The influence of chemoattractants on bacterial motility was investigated. Each film was analysed to determine the swimming velocity of the bacteriavs a...

    work page 2025