Ecological memory of hydrodynamic cues shapes growth and migration of motile microorganisms

1. Problem Statement

Microorganisms in aquatic ecosystems inhabit inherently dynamic environments where fluid motion (shear, turbulence) acts as a primary environmental cue. While it is known that hydrodynamic forces can alter microbial physiology, the concept of "ecological memory"—the lasting imprint of prior environmental conditions on future behavior—remains poorly understood in the context of hydrodynamics.

Specifically, the study addresses three critical gaps:

• Timing Dependence: How does the timing of hydrodynamic exposure relative to the microbial growth stage (lag, exponential, stationary) influence fitness?
• History Dependence (Hysteresis): Do cells retain a "memory" of past flow conditions that alters their response even after the flow is removed?
• Mechanistic Link: How do mechanical cues translate into physiological changes (growth, motility, lipid storage) and how do these changes affect the balance of torques governing gravitactic (gravity-directed) swimming?

2. Methodology

The researchers utilized the motile marine raphidophyte Heterosigma akashiwo (strains HA452 and HA3107) as a model organism. They employed a combination of millifluidic flow control, high-resolution cell tracking, and quantitative growth assays.

Experimental Design

The study utilized two complementary perturbation scenarios to disentangle exposure history from duration:

1. Standard Scenarios (Static $\to$ Dynamic): Cultures were grown statically and then exposed to orbital shaking (110 rpm) at specific time points:
   • 0-hour: Immediate exposure post-inoculation.
   • 120-hour: Exposure during early exponential phase.
   • 168-hour: Exposure during mid-exponential phase.
   • 348-hour: Exposure during stationary phase.
2. Reverse Scenarios (Dynamic $\to$ Static): Cultures were exposed to shaking immediately after inoculation and then returned to static conditions at:
   • 135-hour: Early exponential phase.
   • 168-hour: Mid-exponential phase.

Hydrodynamic Characterization:

• The orbital shaker (110 rpm) generated a mean kinetic energy dissipation rate ($\epsilon$) of $\approx 9.54 \times 10^{-4} \, \text{m}^2\text{s}^{-3}$.
• This corresponds to a Kolmogorov length scale ($\eta$) of $\approx 0.18$ mm and a shear rate ($S_\eta$) of $\sim 30 \, \text{s}^{-1}$.
• Since cell diameters ($10\text{--}20 \, \mu\text{m}$) are much smaller than $\eta$, cells experienced a linear velocity gradient (local shear).

Quantitative Metrics

• Growth Kinetics: Cell concentration was tracked via image analysis and flow cytometry. Data were fitted to a logistic growth model to extract doubling time ($T_d$) and carrying capacity ($K$).
• Motility Analysis:
  • Vertical Velocity ($v_y$): Mean upward swimming speed.
  • Directional Statistics: Wind-rose plots of swimming vectors to assess anisotropy (upward bias).
  • Reorientation Time ($\tau_r$): Calculated by flipping the chamber $180^\circ$ and tracking the relaxation of swimming orientation. $\tau_r = 1/(2A)$, where $A$ is the amplitude of the sinusoidal fit of angular velocity vs. angle. Lower $\tau_r$ indicates stronger gravitactic stability.
• Lipid Accumulation: Nile Red staining was used to quantify intracellular lipid droplets (LDs). Metrics included normalized lipid area and total LD volume. Crucially, the spatial distribution of LDs relative to the cell centroid was mapped to determine the "effective lipid droplet" position, which dictates the gravitational torque.

3. Key Contributions

• Definition of Hydrodynamic Memory: The study establishes that prior exposure to flow creates a "legacy effect" that persists even after the flow ceases, fundamentally altering growth and motility trajectories.
• Strain-Specific Trade-offs: It reveals distinct evolutionary strategies between two strains of the same species:
  • HA452: Prioritizes metabolic robustness (maintains growth/lipid storage under stress) but sacrifices motility stability (loses gravitaxis).
  • HA3107: Prioritizes motility robustness (maintains gravitaxis) but suffers greater fitness costs (reduced growth/biomass) under specific perturbation histories.
• Mechanistic Framework: The paper links mechanical history to lipid droplet reorganization. It demonstrates that hydrodynamic cues alter the spatial distribution of intracellular lipids, shifting the cell's center of mass relative to its geometric center. This shift modulates the gravitational torque, thereby reprogramming the cell's ability to swim against gravity.

4. Key Results

Growth Kinetics

• Standard Scenarios:
  • Early Exposure (0-hour): Accelerated growth (shorter $T_d$) with no change in carrying capacity ($K$) for both strains.
  • Mid-Exposure (120/168-hour): Significant reduction in both growth rate and carrying capacity. HA3107 was more severely affected than HA452.
  • Late Exposure (348-hour): No significant change in growth parameters, though death phase occurred earlier.
• Reverse Scenarios (Memory Effect):
  • HA452: When returned to static conditions after early shaking, cells retained a "memory" of the stress. They showed a secondary lag phase and a reduced carrying capacity, despite the removal of shear.
  • HA3107: Showed a delayed collapse. Cells maintained normal growth while shaken but, upon returning to static conditions after 168 hours, exhibited a drastic reduction in growth rate and biomass yield, indicating a maladaptive regulatory state.

Motility and Gravitaxis

• HA452 (Sensitive to History):
  • 0-hour (Continuous): Cells lost upward bias, became isotropic (random swimming), and exhibited significantly increased reorientation times ($\tau_r$), indicating a loss of gravitactic stability.
  • Reverse Scenarios: Even after flow stopped, HA452 cells failed to recover their upward swimming bias. They remained in a "reprogrammed" state with slow, diffusive motion.
• HA3107 (Robust to History):
  • Maintained strong upward bias and short reorientation times even under continuous shaking (0-hour).
  • In reverse scenarios, HA3107 only showed motility collapse if the shaking duration was extended to 168 hours, suggesting a higher threshold for mechanical reprogramming.

Lipid Remodeling (The Mechanism)

• HA452: Early and continuous shear induced a massive accumulation of lipid droplets that relocated toward the geometric center of the cell.
  • Consequence: This reduced the "bottom-heaviness" (offset between center of mass and geometric center), weakening the gravitational torque. The cell became mechanically unstable, leading to the observed loss of gravitaxis.
• HA3107: Maintained lipid droplets in a configuration that preserved a significant offset from the centroid, sustaining strong gravitational torque and stable upward swimming, even under shear.

5. Significance

• Predictive Ecology: The findings challenge models that assume microbial responses depend solely on current environmental conditions. Instead, hydrodynamic history is a critical control parameter. This is vital for predicting phytoplankton blooms, vertical migration patterns, and carbon cycling in a changing climate where wind forcing and stratification are shifting.
• Mechanistic Insight: By linking fluid mechanics to intracellular lipid architecture, the study provides a physical explanation for "ecological memory." It shows that mechanical stress can trigger irreversible (on experimental timescales) metabolic re-allocation that alters the physical properties of the cell.
• Evolutionary Trade-offs: The work highlights a trade-off between metabolic performance (growth/lipid storage) and mechanical robustness (motility stability). Different strains occupy different regions of this trade-off space, suggesting that hydrodynamic regimes act as a selective pressure shaping phenotypic diversity.
• Engineering Applications: Understanding these memory effects is crucial for designing bioreactors and engineered systems where flow regimes fluctuate, as the "history" of the flow will dictate the productivity and behavior of the microbial culture.

In conclusion, the paper demonstrates that for motile microorganisms, the past is not just a precursor to the present; it is a structural component of the cell's future physiology and behavior, encoded through the mechanical reorganization of intracellular components.