Modelling and analysis of bacterial tracks suggest an active reorientation mechanism in Rhodobacter sphaeroides

Most free-swimming bacteria move in approximately straight lines, interspersed with random reorientation phases. A key open question concerns varying mechanisms by which reorientation occurs. We combine mathematical modelling with analysis of a large tracking dataset to study the poorly understood reorientation mechanism in the monoflagellate species Rhodobacter sphaeroides. The flagellum on this species rotates counterclockwise to propel the bacterium, periodically ceasing rotation to enable reorientation. When rotation restarts the cell body usually points in a new direction. It has been assumed that the new direction is simply the result of Brownian rotation. We consider three variants of a self-propelled particle model of bacterial motility. The first considers rotational diffusion only, corresponding to a non-chemotactic mutant strain. Two further models incorporate stochastic reorientations, describing ‘run-and-tumble’ motility. We derive expressions for key summary statistics and simulate each model using a stochastic computational algorithm. We also discuss the effect of cell geometry on rotational diffusion. Working with a previously published tracking dataset, we compare predictions of the models with data on individual stopping events in R. sphaeroides. This provides strong evidence that this species undergoes some form of active reorientation rather than simple reorientation by Brownian rotation.     

Reconstructing the flight kinematics of swarming and mating in wild mosquitoes

We describe a novel tracking system for reconstructing three-dimensional tracks of individual mosquitoes in wild swarms and present the results of validating the system by filming swarms and mating events of the malaria mosquito Anopheles gambiae in Mali. The tracking system is designed to address noisy, low frame-rate (25 frames per second) video streams from a stereo camera system. Because flying A. gambiae move at 1–4 m s⁻¹, they appear as faded streaks in the images or sometimes do not appear at all. We provide an adaptive algorithm to search for missing streaks and a likelihood function that uses streak endpoints to extract velocity information. A modified multi-hypothesis tracker probabilistically addresses occlusions and a particle filter estimates the trajectories. The output of the tracking algorithm is a set of track segments with an average length of 0.6–1 s. The segments are verified and combined under human supervision to create individual tracks up to the duration of the video (90 s). We evaluate tracking performance using an established metric for multi-target tracking and validate the accuracy using independent stereo measurements of a single swarm. Three-dimensional reconstructions of A. gambiae swarming and mating events are presented.     

Dynamics and robustness of the cardiac progenitor cell induced pluripotent stem cell network during cell phenotypes transition

Robustness is a fundamental characteristic of biological systems since all living systems need to adapt to internal or external perturbations, unpredictable environments, stochastic events and unreliable components, and so on. A long‐term challenge in systems biology is to reveal the origin of robustness underlying molecular regulator network. In this study, a simple Boolean model is used to investigate the global dynamic properties and robustness of cardiac progenitor cell (CPC) induced pluripotent stem cell network that governs reprogramming and directed differentiation process. It is demonstrated that two major attractors correspond to source and target cell phenotypes, respectively, and two dominating attracting trajectories characterise the biological pathways between two major cell phenotypes. In particular, the experimentally observed transition between different cell phenotypes can be reproduced and explained theoretically. Furthermore, the robustness of major attractors and trajectories is largely maintained with respect to small perturbations to the network. Taken together, the CPC‐induced pluripotent stem cell network is extremely robustly designed for their functions.Inspec keywords: cellular biophysics, Boolean functions, perturbation theory, molecular biophysics, cardiologyOther keywords: cardiac progenitor cell induced pluripotent stem cell network, cell phenotypes transition, biological systems, living systems, internal perturbations, external perturbations, unpredictable environments, stochastic events, unreliable components, long‐term challenge, systems biology, molecular regulator network, Boolean model, global dynamic properties, directed differentiation process, CPC‐induced pluripotent stem cell network     

Three-dimensional traction forces of Schwann cells on compliant substrates

The mechanical interaction between Schwann cells (SCs) and their microenvironment is crucial for the development, maintenance and repair of the peripheral nervous system. In this paper, we present a detailed investigation on the mechanosensitivity of SCs across a physiologically relevant substrate stiffness range. Contrary to many other cell types, we find that the SC spreading area and cytoskeletal actin architecture were relatively insensitive to substrate stiffness with pronounced stress fibre formation across all moduli tested (0.24–4.80 kPa). Consistent with the presence of stress fibres, we found that SCs generated large surface tractions on stiff substrates and large, finite material deformations on soft substrates. When quantifying the three-dimensional characteristics of the SC traction profiles, we observed a significant contribution from the out-of-plane traction component, locally giving rise to rotational moments similar to those observed in mesenchymal embryonic fibroblasts. Taken together, these measurements provide the first set of quantitative biophysical metrics of how SCs interact with their physical microenvironment, which are anticipated to aid in the development of tissue engineering scaffolds designed to promote functional integration of SCs into post-injury in vivo environments.     

Rear actomyosin contractility-driven directional cell migration in three-dimensional matrices: a mechano-chemical coupling mechanism

Cell migration is of vital importance in many biological processes, including organismal development, immune response and development of vascular diseases. For instance, migration of vascular smooth muscle cells from the media to intima is an essential part of the development of atherosclerosis and restenosis after stent deployment. While it is well characterized that cells use actin polymerization at the leading edge to propel themselves to move on two-dimensional substrates, the migration modes of cells in three-dimensional matrices relevant to in vivo environments remain unclear. Intracellular tension, which is created by myosin II activity, fulfils a vital role in regulating cell migration. We note that there is compelling evidence from theoretical and experimental work that myosin II accumulates at the cell rear, either isoform-dependent or -independent, leading to three-dimensional migration modes driven by posterior myosin II tension. The scenario is not limited to amoeboid migration, and it is also seen in mesenchymal migration in which a two-dimensional-like migration mode based on front protrusions is often expected, suggesting that there may exist universal underlying mechanisms. In this review, we aim to shed some light on how anisotropic myosin II localization induces cell motility in three-dimensional environments from a biomechanical view. We demonstrate an interesting mechanism where an interplay between mechanical myosin II recruitment and biochemical myosin II activation triggers directional migration in three-dimensional matrices. In the case of amoeboid three-dimensional migration, myosin II first accumulates at the cell rear to induce a slight polarization displayed as a uropod-like structure under the action of a tension-dependent mechanism. Subsequent biochemical signalling pathways initiate actomyosin contractility, producing traction forces on the adhesion system or creating prominent motile forces through blebbing activity, to drive cells to move. In mesenchymal three-dimensional migration, cells can also take advantage of the elastic properties of three-dimensional matrices to move. A minor myosin isoform, myosin IIB, is retained by relatively stiff three-dimensional matrices at the posterior side, then activated by signalling cascades, facilitating prominent cell polarization by establishing front–back polarity and creating cell rear. Myosin IIB initiates cell polarization and coordinates with the major isoform myosin IIA-assembled stress fibres, to power the directional migration of cells in the three-dimensional matrix.     

Collective cell migration: leadership,invasion and segregation

A number of biological processes, such as embryo development, cancer metastasis or wound healing, rely on cells moving in concert. The mechanisms leading to the emergence of coordinated motion remain however largely unexplored. Although biomolecular signalling is known to be involved in most occurrences of collective migration, the role of physical and mechanical interactions has only been recently investigated. In this study, a versatile framework for cell motility is implemented in silico in order to study the minimal requirements for the coordination of a group of epithelial cells. We find that cell motility and cell–cell mechanical interactions are sufficient to generate a broad array of behaviours commonly observed in vitro and in vivo. Cell streaming, sheet migration and susceptibility to leader cells are examples of behaviours spontaneously emerging from these simple assumptions, which might explain why collective effects are so ubiquitous in nature. The size of the population and its confinement appear, in particular, to play an important role in the coordination process. In all cases, the complex response of the population can be predicted from the knowledge of the correlation length of the velocity field measured in the bulk of the epithelial layer. This analysis provides also new insights into cancer metastasis and cell sorting, suggesting, in particular, that collective invasion might result from an emerging coordination in a system where single cells are mechanically unable to invade.     

Turbulence triggers vigorous swimming but hinders motion strategy in planktonic copepods

Calanoid copepods represent a major component of the plankton community. These small animals reside in constantly flowing environments. Given the fundamental role of behaviour in their ecology, it is especially relevant to know how copepods perform in turbulent flows. By means of three-dimensional particle tracking velocimetry, we reconstructed the trajectories of hundreds of adult Eurytemora affinis swimming freely under realistic intensities of homogeneous turbulence. We demonstrate that swimming contributes substantially to the dynamics of copepods even when turbulence is significant. We show that the contribution of behaviour to the overall dynamics gradually reduces with turbulence intensity but regains significance at moderate intensity, allowing copepods to maintain a certain velocity relative to the flow. These results suggest that E. affinis has evolved an adaptive behavioural mechanism to retain swimming efficiency in turbulent flows. They suggest the ability of some copepods to respond to the hydrodynamic features of the surrounding flow. Such ability may improve survival and mating performance in complex and dynamic environments. However, moderate levels of turbulence cancelled gender-specific differences in the degree of space occupation and innate movement strategies. Our results suggest that the broadly accepted mate-searching strategies based on trajectory complexity and movement patterns are inefficient in energetic environments.     

Mechanics regulates ATP-stimulated collective calcium response in fibroblast cells

Cells constantly sense their chemical and mechanical environments. We study the effect of mechanics on the ATP-induced collective calcium response of fibroblast cells in experiments that mimic various tissue environments. We find that closely packed two-dimensional cell cultures on a soft polyacrylamide gel (Young''s modulus E = 690 Pa) contain more cells exhibiting calcium oscillations than cultures on a rigid substrate (E = 36 000 Pa). Calcium responses of cells on soft substrates show a slower decay of calcium level relative to those on rigid substrates. Actin enhancement and disruption experiments for the cell cultures allow us to conclude that actin filaments determine the collective Ca²⁺ oscillatory behaviour in the culture. Inhibition of gap junctions results in a decrease of the oscillation period and reduced correlation of calcium responses, which suggests additional complexity of signalling upon cell–cell contact. Moreover, the frequency of calcium oscillations is independent of the rigidity of the substrate but depends on ATP concentration. We compare our results with those from similar experiments on individual cells. Overall, our observations show that collective chemical signalling in cell cultures via calcium depends critically on the mechanical environment.     

Cost–benefit analysis of the mechanisms that enable migrating cells to sustain motility upon changes in matrix environments

Cells can move through extracellular environments with varying geometries and adhesive properties. Adaptation to these differences is achieved by switching between different modes of motility, including lamellipod-driven and blebbing motility. Further, cells can modulate their level of adhesion to the extracellular matrix (ECM) depending on both the level of force applied to the adhesions and cell intrinsic biochemical properties. We have constructed a computational model of cell motility to investigate how motile cells transition between extracellular environments with varying surface continuity, confinement and adhesion. Changes in migration strategy are an emergent property of cells as the ECM geometry and adhesion changes. The transition into confined environments with discontinuous ECM fibres is sufficient to induce shifts from lamellipod-based to blebbing motility, while changes in confinement alone within a continuous geometry are not. The geometry of the ECM facilitates plasticity, by inducing shifts where the cell has high marginal gain from a mode change, and conserving persistency where the cell can continue movement regardless of the motility mode. This regulation of cell motility is independent of global changes in cytoskeletal properties, but requires locally higher linkage between the actin network and the plasma membrane at the cell rear, and changes in internal cell pressure. In addition to matrix geometry, we consider how cells might transition between ECM of different adhesiveness. We find that this requires positive feedback between the forces cells apply on the adhesion points, and the strength of the cell–ECM adhesions on those sites. This positive feedback leads to the emergence of a small number of highly adhesive cores, similar to focal adhesions. While the range of ECM adhesion levels the cell can invade is expanded with this feedback mechanism; the velocities are lowered for conditions where the positive feedback is not vital. Thus, plasticity of cell motility sacrifices the benefits of specialization, for robustness.     

How much information can be obtained from tracking the position of the leading edge in a scratch assay?

Moving cell fronts are an essential feature of wound healing, development and disease. The rate at which a cell front moves is driven, in part, by the cell motility, quantified in terms of the cell diffusivity D, and the cell proliferation rate λ. Scratch assays are a commonly reported procedure used to investigate the motion of cell fronts where an initial cell monolayer is scratched, and the motion of the front is monitored over a short period of time, often less than 24 h. The simplest way of quantifying a scratch assay is to monitor the progression of the leading edge. Use of leading edge data is very convenient because, unlike other methods, it is non-destructive and does not require labelling, tracking or counting individual cells among the population. In this work, we study short-time leading edge data in a scratch assay using a discrete mathematical model and automated image analysis with the aim of investigating whether such data allow us to reliably identify D and λ. Using a naive calibration approach where we simply scan the relevant region of the (D, λ) parameter space, we show that there are many choices of D and λ for which our model produces indistinguishable short-time leading edge data. Therefore, without due care, it is impossible to estimate D and λ from this kind of data. To address this, we present a modified approach accounting for the fact that cell motility occurs over a much shorter time scale than proliferation. Using this information, we divide the duration of the experiment into two periods, and we estimate D using data from the first period, whereas we estimate λ using data from the second period. We confirm the accuracy of our approach using in silico data and a new set of in vitro data, which shows that our method recovers estimates of D and λ that are consistent with previously reported values except that that our approach is fast, inexpensive, non-destructive and avoids the need for cell labelling and cell counting.     

Comparison of linear and classical velocity update rules in particle swarm optimization: notes on diversity

In this study, we investigate the significance of diversity in the particle swarm optimization (PSO) algorithm. To do so, we study two different implementations of the PSO, being the so‐called linear and classical PSO formulations. While the behaviour of these two implementations is markedly different, they only differ in the formulation of the velocity update rule. In fact, the differences are merely due to subtle differences in the introduction of randomness into the algorithm. In this paper, we show that in algorithms employing the linear PSO velocity update rule, particle trajectories collapse to line searches in n‐dimensional space. The classical formulation does not suffer this drawback. Instead, directional diverse stochastic search trajectories are retained, which in turn helps to alleviate premature convergence. The performance of the classical implementation is therefore superior for all test problems considered, due to the presence of adequate diversity. Copyright © 2006 John Wiley & Sons, Ltd.     

Relevance of the sterilization-induced effects on the properties of different hydroxyapatite nanoparticles and assessment of the osteoblastic cell response

Hydroxyapatite (Hap) is a calcium phosphate with a chemical formula that closely resembles that of the mineral constituents found in hard tissues, thereby explaining its natural biocompatibility and wide biomedical use. Nanostructured Hap materials appear to present a good performance in bone tissue applications because of their ability to mimic the dimensions of bone components. However, bone cell response to individual nanoparticles and/or nanoparticle aggregates lost from these materials is largely unknown and shows great variability. This work addresses the preparation and characterization of two different Hap nanoparticles and their interaction with osteoblastic cells. Hap particles were produced by a wet chemical synthesis (WCS) at 37°C and by hydrothermal synthesis (HS) at 180°C. As the ultimate in vivo applications require a sterilization step, the synthesized particles were characterized ‘as prepared’ and after sterilization (autoclaving, 120°C, 20 min). WCS and HS particles differ in their morphological (size and shape) and physicochemical properties. The sterilization modified markedly the shape, size and aggregation state of WCS nanoparticles. Both particles were readily internalized by osteoblastic cells by endocytosis, and showed a low intracellular dissolution rate. Concentrations of WCS and HS particles less than 500 μg ml⁻¹ did not affect cell proliferation, F-actin cytoskeleton organization and apoptosis rate and increased the gene expression of alkaline phosphatase and BMP-2. The two particles presented some differences in the elicited cell response. In conclusion, WCS and HS particles might exhibit an interesting profile for bone tissue applications. Results suggest the relevance of a proper particle characterization, and the interest of an individual nanoparticle targeted research.     

Role of cell polarity dynamics and motility in pattern formation due to contact-dependent signalling

A key challenge in biology is to understand how spatio-temporal patterns and structures arise during the development of an organism. An initial aggregate of spatially uniform cells develops and forms the differentiated structures of a fully developed organism. On the one hand, contact-dependent cell–cell signalling is responsible for generating a large number of complex, self-organized, spatial patterns in the distribution of the signalling molecules. On the other hand, the motility of cells coupled with their polarity can independently lead to collective motion patterns that depend on mechanical parameters influencing tissue deformation, such as cellular elasticity, cell–cell adhesion and active forces generated by actin and myosin dynamics. Although modelling efforts have, thus far, treated cell motility and cell–cell signalling separately, experiments in recent years suggest that these processes could be tightly coupled. Hence, in this paper, we study how the dynamics of cell polarity and migration influence the spatiotemporal patterning of signalling molecules. Such signalling interactions can occur only between cells that are in physical contact, either directly at the junctions of adjacent cells or through cellular protrusional contacts. We present a vertex model which accounts for contact-dependent signalling between adjacent cells and between non-adjacent neighbours through long protrusional contacts that occur along the orientation of cell polarization. We observe a rich variety of spatiotemporal patterns of signalling molecules that is influenced by polarity dynamics of the cells, relative strengths of adjacent and non-adjacent signalling interactions, range of polarized interaction, signalling activation threshold, relative time scales of signalling and polarity orientation, and cell motility. Though our results are developed in the context of Delta–Notch signalling, they are sufficiently general and can be extended to other contact dependent morpho-mechanical dynamics.     

Coordination of gene expression noise with cell size: analytical results for agent-based models of growing cell populations

The chemical master equation and the Gillespie algorithm are widely used to model the reaction kinetics inside living cells. It is thereby assumed that cell growth and division can be modelled through effective dilution reactions and extrinsic noise sources. We here re-examine these paradigms through developing an analytical agent-based framework of growing and dividing cells accompanied by an exact simulation algorithm, which allows us to quantify the dynamics of virtually any intracellular reaction network affected by stochastic cell size control and division noise. We find that the solution of the chemical master equation—including static extrinsic noise—exactly agrees with the agent-based formulation when the network under study exhibits stochastic concentration homeostasis, a novel condition that generalizes concentration homeostasis in deterministic systems to higher order moments and distributions. We illustrate stochastic concentration homeostasis for a range of common gene expression networks. When this condition is not met, we demonstrate by extending the linear noise approximation to agent-based models that the dependence of gene expression noise on cell size can qualitatively deviate from the chemical master equation. Surprisingly, the total noise of the agent-based approach can still be well approximated by extrinsic noise models.     

The influence of surface morphology and rigidity of the substrata on cell motility

The purpose of this study is to measure the cell motility under different stiffness and pattern of the substrata. Human melanoma cells were used for this study. Three surface patterns including flat, 6 µm-cone, and 6 µm-groove on polydimethylsiloxane (PDMS) were created by the micro-electro-mechanical system (MEMS). Glass dish was used as control substrate. After cells were seeded for 30 min, the cell images were taken every minute for 2 h. Each group had 4-5 dishes and 27-38 cells were calculated. The cell motility was 1.13 ± 1.30, 2.12 ± 1.21, 2.39 ± 1.11 and 3.08 ± 1.49 µm/min on glass, flat, cone, and groove PDMS, respectively. Cells in PDMS groups moved significantly faster than the control group (glass) due to smaller stiffness of the former substrates. More than 80% of cells on grooved PDMS moved along the grooves, indicating the grooved surface morphology could control the direction of cell movement. Our results display that substrate modulus and pattern can influence cell motility.     

Novel temperature-responsive polymer brushes with carbohydrate residues facilitate selective adhesion and collection of hepatocytes

Temperature-responsive glycopolymer brushes were designed to investigate the effects of grafting architectures of the copolymers on the selective adhesion and collection of hypatocytes. Homo, random and block sequences of N-isopropylacrylamide and 2-lactobionamidoethyl methacrylate were grafted on glass substrates via surface-initiated atom transfer radical polymerization. The galactose/lactose-specific lectin RCA₁₂₀ and HepG2 cells were used to test for specific recognition of the polymer brushes containing galactose residues over the lower critical solution temperatures (LCSTs). RCA₁₂₀ showed a specific binding to the brush surfaces at 37 °C. These brush surfaces also facilitated the adhesion of HepG2 cells at 37 °C under nonserum conditions, whereas no adhesion was observed for NIH-3T3 fibroblasts. When the temperature was decreased to 25 °C, almost all the HepG2 cells detached from the block copolymer brush, whereas the random copolymer brush did not release the cells. The difference in releasing kinetics of cells from the surfaces with different grafting architectures can be explained by the correlated effects of significant changes in LCST, mobility, hydrophilicity and mechanical properties of the grafted polymer chains. These findings are important for designing ‘on–off’ cell capture/release substrates for various biomedical applications such as selective cell separation.     

Enhanced antimicrobial activity of silver nanoparticles‐Lonicera Japonica Thunb combo

Silver metals have long been known to possess antimicrobial properties. Recently, even the nanoparticle version of silver (AgNPs) has also been established as antimicrobials. In this study AgNPs were combined with extracts of the medicinal plant Chinese honeysuckle, Lonicera japonica Thunb. The antimicrobial activity of the AgNPs‐herb was tested against pathogenic Escherichia coli CMCC44113. Using different AgNPs or herb (honeysuckle water extract or HWE) ratios in the presence of a fixed concentration of E. coli CMCC44113, potencies were found to be proportional with concentrations. The antimicrobial activities of AgNPs‐HWE combo were significant enhanced, when compared with solely AgNPs or HWE. Thus, atomic force microscopic and propidium monoazide‐PCR were used to probe the damages caused by AgNPs‐HWE combo on the cell morphology and cell membrane integrity of E. coli. The mechanism of AgNPs‐HWE combo against E. coli may attribute to AgNPs leads to cell wall lysis and damages cell membrane integrity, and thus increases the penetration of HWE into the bacterium, which results in more serious damage to bacterial cells. These findings indicated that AgNPs‐herb was more potent than the AgNPs alone and holds promise for the development of nanoparticle enhanced herbal pharmaceuticals.Inspec keywords: microorganisms, cellular biophysics, silver, nanoparticles, nanomedicine, antibacterial activity, vegetation, biomembranes, biomedical materialsOther keywords: enhanced antimicrobial activity, silver nanoparticle‐Lonicera japonica Thunb combo, silver metals, antimicrobial properties, medicinal plant Chinese honeysuckle, AgNP‐herb, pathogenic Escherichia coli CMCC44113, E. coli CMCC44113 concentration, antimicrobial activity, AgNP‐HWE combo, atomic force microscopy, propidium monoazide‐PCR, cell morphology, E. coli, cell wall lysis, cell membrane integrity, HWE penetration, bacterial cells, nanoparticle enhanced herbal pharmaceuticals, Ag     

Sources of variability in cytosolic calcium transients triggered by stimulation of homogeneous uro-epithelial cell monolayers

Epithelial tissue structure is the emergent outcome of the interactions between large numbers of individual cells. Experimental cell biology offers an important tool to unravel these complex interactions, but current methods of analysis tend to be limited to mean field approaches or representation by selected subsets of cells. This may result in bias towards cells that respond in a particular way and/or neglect local, context-specific cell responses. Here, an automated algorithm was applied to examine in detail the individual calcium transients evoked in genetically homogeneous, but asynchronous populations of cultured non-immortalized normal human urothelial cells when subjected to either the global application of an external agonist or a localized scratch wound. The recorded calcium transients were classified automatically according to a set of defined metrics and distinct sub-populations of cells that responded in qualitatively different ways were observed. The nature of this variability in the homogeneous cell population was apportioned to two sources: intrinsic variation in individual cell responses and extrinsic variability due to context-specific factors of the environment, such as spatial heterogeneity. Statistically significant variation in the features of the calcium transients evoked by scratch wounding according to proximity to the wound edge was identified. The manifestation of distinct sub-populations of cells is considered central to the coordination of population-level response resulting in wound closure.     

Shear-induced orientational dynamics and spatial heterogeneity in suspensions of motile phytoplankton

Fluid flow, ubiquitous in natural and man-made environments, has the potential to profoundly impact the transport of microorganisms, including phytoplankton in aquatic habitats and bioreactors. Yet, the effect of ambient flow on the swimming behaviour of phytoplankton has remained poorly understood, largely owing to the difficulty of observing cell–flow interactions at the microscale. Here, we present microfluidic experiments where we tracked individual cells for four species of motile phytoplankton exposed to a spatially non-uniform fluid shear rate, characteristic of many flows in natural and artificial environments. We observed that medium-to-high mean shear rates (1–25 s⁻¹) produce heterogeneous cell concentrations in the form of regions of accumulation and regions of depletion. The location of these regions relative to the flow depends on the cells'' propulsion mechanism, body shape and flagellar arrangement, as captured by an effective aspect ratio. Species having a large effective aspect ratio accumulated in the high-shear regions, owing to shear-induced alignment of the swimming orientation with the fluid streamlines. Species having an effective aspect ratio close to unity exhibited little preferential accumulation at low-to-moderate flow rates, but strongly accumulated in the low-shear regions under high flow conditions, potentially owing to an active, behavioural response of cells to shear. These observations demonstrate that ambient fluid flow can strongly affect the motility and spatial distribution of phytoplankton and highlight the rich dynamics emerging from the interaction between motility, morphology and flow.