Phenotyping single-cell motility in microfluidic confinement
Abstract
The movement trajectories of organisms serve as dynamic read-outs of their behaviour and physiology. For microorganisms this can be difficult to resolve due to their small size and fast movement. Here, we devise a novel droplet microfluidics assay to encapsulate single micron-sized algae inside closed arenas, enabling ultralong high-speed tracking of the same cell. Comparing two model species - Chlamydomonas reinhardtii (freshwater, 2 cilia), and Pyramimonas octopus (marine, 8 cilia), we detail their highly-stereotyped yet contrasting swimming behaviours and environmental interactions. By measuring the rates and probabilities with which cells transition between a trio of motility states (smooth-forward swimming, quiescence, tumbling or excitable backward swimming), we reconstruct the control network that underlies this gait switching dynamics. A simplified model of cell-roaming in circular confinement reproduces the observed long-term behaviours and spatial fluxes, including novel boundary circulation behaviour. Finally, we establish an assay in which pairs of droplets are fused on demand, one containing a trapped cell with another containing a chemical that perturbs cellular excitability, to reveal how aneural microorganisms adapt their locomotor patterns in real-time.
Editor's evaluation
This paper reports on the development of an impressive microfluidic platform for the study of motility, and motility transitions, exhibited by single algal cells in circular confinement. Building on previous work that showed a three-state motility repertoire for certain green algae, the present work uses extremely long time series and a variety of physical perturbations to show how those dynamics can be altered by environmental conditions. The work will be of interest to a wide range of scientists studying motility and non-equilibrium dynamics.
https://doi.org/10.7554/eLife.76519.sa0Introduction
All lifeforms are environmentally intelligent; even single cells sense and respond to their environment (Jennings, 1901). How living systems interact with the external world is a fundamental question that has fascinated researchers for centuries. The exploratory behaviours of living beings can be extraordinarily complex and context-dependent, and yet also stereotyped. Digitally-assisted and non-invasive tracking methods have revolutionised the study of behaving animals (Berman, 2018; Mathis et al., 2018). Algorithms for pose and gait-estimation based on computer vision and machine learning have bolstered the utility of animal models for neuroscience, providing a highly quantitative and objective framework for measuring behaviour. These nascent approaches are now being used in combination with genetic manipulation, to resolve how an organism’s sensory architecture (neural pathways, nervous systems) may be coupled to motor actuators (e.g. muscles) to drive motor actions (Brown et al., 2013).
Surprisingly, nervous systems are not necessary for the emergence of complex search patterns or behavioural motifs (Coyle et al., 2019; Wan, 2020). Micoorganisms also adopt highly stereotyped locomotion strategies, often employing motile appendages: bacterial flagella for swimming, pili for twitching and swarming, and in many eukaryotes, cilia for self-propulsion. Microbial motility is strongly influenced by boundaries, interfaces, and other obstacles (Ostapenko et al., 2018; Kantsler et al., 2013). Microorganisms naturally encounter heterogeneous environments, such as soil, sea ice, or porous materials that constitute complex networks of confining tubes, spaces, or interstices. By replicating such microenvironments in the laboratory, realistic situations can be created for tracking and understanding the mechanisms of movement and behavioural adaptation (Théry et al., 2021; Bhattacharjee and Datta, 2019).
Despite their small size, microswimmers self-propel at speeds of several tens (even hundreds) of body lengths per second, posing a formidable challenge for live-cell microscopy and tracking. Typically, the same individual can only be observed for relatively short periods, and motility statistics are then coarse-grained across entire populations. Across different species, even across different individuals of the same clonal population, microbes exhibit significantly divergent behaviours. This single-cell heterogeneity may be hidden if only population averages are considered. However, long-term tracking of the same individual was not possible until very recently (Hageman et al., 2019; Ostapenko et al., 2018; Codutti et al., 2022).
What are the behavioural signatures that distinguish one microswimmer from another? To address this, we phenotype two species of motile algae, and decode their responses to diverse environmental cues. The first, Chlamydomonas reinhardtii (hereafter CR) is a freshwater biflagellate and a model organism for motility studies and cilia biology (Sasso et al., 2018). The second, Pyramimonas octopus (hereafter PO), is a marine octoflagellate, which has a unique tripartite run-stop-shock motility pattern (Wan and Goldstein, 2018). Harnessing modern microfluidics technology to access long-term behaviour (Brückner et al., 2019; Tweedy et al., 2020), we developed a platform to stably encapsulate algae inside droplets contained in shallow microwells. Here, cells can explore and interact with their local environment over much longer timescales than previously possible. Single-cell responses were evaluated against physical confinement, light, and chemical stimulation. For the latter, this is achieved by fabricating a tailored device to directly perturb the chemical environment of a trapped cell, by droplet pairing and fusion controlled by surfactant-induced interface destabilization (Niu et al., 2009; Frenz et al., 2008; Fradet et al., 2011).
From these long-time recordings, we then explore universal metrics of microbial motility, including speed, gait switching probabilities, and emergent spatial probability fluxes in the time-averaged motion of single cells. We show how these behaviours are effectively captured by a model of a roaming cell in a circular corral. Taken together, our approach integrates experimental and computational techniques to provide an exciting opportunity to understand how individual living cells interact with local environmental cues at the microscopic scale.
Results
Microfluidic trapping of single motile algae
We encapsulated CR or PO cells in droplet microwells of four different diameters (Ø = 40, 60, 120, 200 μm). To study the cells’ baseline motility (in the absence of stimulation), we imaged in red light to prevent phototactic effects (long-pass IR filter). According to the absorption spectra of CR photoreceptors (ChR1 and ChR2), cells are insensitive to wavelengths of > (Harz and Hegemann, 1991). The identity of the photoreceptor in PO has not yet been confirmed, but is likely to be similar to CR, since green algal eyespots are highly conserved (Kreimer, 2009).
Only wells containing precisely one cell were imaged (see Materials and methods). These water-in-oil droplets were stably trapped and retained their shape and volume for more than 1 hr. The trapping procedure was sufficiently gentle and did not damage the cilia or cause spontaneous deciliation. Compared to the cell-scale, the arenas are shallow (∼ 30 μm deep in total) and quasi-2D, so that swimming motion is largely restricted to the plane (Figure 1A and B). We trapped each cell for a total of 1 hr, recording at 5-min intervals ( cells per trap size, 6 time points each) at 500 frames per second. High-speed imaging was necessary to resolve rapid motility changes. Cell-centroids were tracked with Trackmate in ImageJ and analysed with custom scripts (see Materials and methods).
The two species display very different swimming behaviours. Typical trajectories and speeds are displayed in Figure 2. Both species have front-mounted cilia (Figure 1C and D). CR executes a breaststroke gait where the two cilia alternate between long-periods of in-phase synchrony, and phase slips when one cilium transiently undergoes altered beating (Wan et al., 2014). This leads to noisy trajectories with straight runs whenever the cilia are synchronised, and turns whenever synchrony is lost. Asynchronous cilia dynamics are associated with reduced swimming speed (Polin et al., 2009). Meanwhile PO swims forward using a rotary breaststroke with frequent episodes of ultrafast backward swimming (called ’shocks’), where all eight cilia undulate synchronously in front of the cell for up to (Wan and Goldstein, 2018). During shocks cell speed can exceed .
Both species displayed stochastic transitions between different swimming modes (Figure 2). However, CR swims largely continuously albeit with noisy fluctuations in speed (fractal timeseries), PO exhibits excitable dynamics reminiscent of neuronal spiking, where bursts of high-speed movement are followed by periods of quiescence. This more episodic motility repertoire of PO is associated with more sparse trap coverage over time (Figure 2A and D) (Videos 1 and 2). The distinction persists in the long-time single-cell heatmaps (Figure 2C and F). In the smallest traps, PO cells showed bursting dynamics - with clusters of shocks in quick succession (Figure 2E). Over 1 hr, CR motility remained largely unchanged, in PO a slight reduction in activity was observed in the smallest traps (Figure 2—figure supplement 1).
Effect of physical confinement
We assayed four trap sizes (Ø = 40, 60, 120, 200 μm), varying from strong confinement where the cell is never more than 1.5 body lengths away from the boundary, to weak confinement where cells are free to roam along tight helices. Population-level trends and stereotyped behaviours were evident when we averaged trajectory statistics over all cells, across all timepoints.
Our findings are twofold. First, confinement affects cell swimming speed. As trap size is decreased, CR speeds are shifted towards lower values (Figure 3A, whereas in PO the characteristic speed and likelihood of runs are reduced (Figure 3B), compare location and width of the violin plot ’waist’). See also bivariate histograms of linear and angular speed (Figure 3—figure supplement 2). Second, increasing confinement changes the nature of how cells explore the circular traps. To see this, we partition the arena (scaled by trap size) and compute the relative probability density of occupancy in each box. Averaging over all trajectories, we find that cells preferentially move within an annulus of the solid boundary (Figure 3C, D - shaded blue regions). To account for the ‘arrow of time’ in the movement history of a single individual (Battle et al., 2016; Cammann et al., 2021), we also compute for each trap size a time-averaged probability flux at each position (see Appendix 1). These steady-state fluxes have a heading and strength according to the most probable trajectory direction starting from the current position. Unexpectedly, single-cell trajectories display non-equilibrium flux loops resulting from persistent circling behaviour inside the traps (Figure 3C, D). Circulation is predominantly counter-clockwise (CCW), viewed from above in the lab frame.
To explain why the cell not only follows the trap boundary, but also how a preferred circling direction can emerge from a perfectly symmetrical trap, we turn to a minimal computational model (Cammann et al., 2021). For brevity, we consider here the case of a CR-like microswimmer only. Briefly, the swimmer is modelled as an active dumbbell with unequal-sized spheres, whose dynamics obey overdamped Langevin equations, with a forcing term to account for steric interactions with the solid boundary (see Appendix 1). We further assume cells swim with a constant speed of 80 μms-1 in all traps. This results in boundary circulation behaviour without recourse to change the cell’s internal motility state (Ostapenko et al., 2018). In domains with constant curvature (e.g. circular), however, macroscopic flux loops should not arise, due to equal probabilities of circling clockwise or counter-clockwise. A further source of asymmetry must be present to account for the experimental observations. Closer inspection of high-magnification videos of single CR cells revealed that the two flagella often lose synchrony, one flagellum consistently beating faster than the other leading to bilateral asymmetry (Figure 3E). Motivated by this, we modify the dumbbell model so that the two spheres are misaligned with an offset angle (Figure 3F). For a small offset (where is the CCW direction), we simulated single-swimmer trajectories for the four trap sizes and evaluated the long-time behaviour. Radial probability densities, mean-squared displacement (MSD) profiles, and mean tangential fluxes of the simulated tracks agree rather well with measured data (Figure 3G-I). In the simulations, chiral fluxes became more prominent with increasing confinement, consistent with the data (Figure 3I).
Effect of white light stimulation
Next we explore the effect of an orthogonal parameter on motility - illumination wavelength. We imaged the samples with broad-spectrum white light (WL) instead of red light (RL), at an intensity expected to influence phototaxis in CR, but not sufficient for photoshock. For two trap sizes, 120 μm and 60 μm (again per trap size), we again recorded motility over 1 hour (Videos 3 and 4).
Some single-cell heterogeneity is evident (Figure 4). At early times CR cells swim faster in WL compared to RL in the same trap sizes at equivalent times, suggesting a photokinetic response (Figure 4A, Figure 4—figure supplement 1A). A decline in motility is observed across all sampled cells after min, especially in the smallest 60 μm traps (note wide ’waists’ at zero speed). All cells in the 60 μm traps were completely immotile after min. Similarly, PO became increasingly quiescent and immotile over time. By 1 hr, all cells stopped moving (Figure 4B), and in some cases deformation of cell shape and even deciliation was noted at the final time point. This was not seen in RL (Figure 2, Figure 2—figure supplement 1). Comparing PO speeds in WL and RL at early times, maximum speeds during ciliary reversals or shocks are unchanged, but their distributions are skewed toward higher speeds due to increased likelihood of forward runs ( μms-1) (Figure 4—figure supplement 1B).
As before, we compute the spatial occupancy of single-cell tracks and their associated steady-state (time-averaged) fluxes (Figure 4C). Again, directed flux loops emerged, due to circling behaviour, most prominently for CR. For both species, the flux strengths are larger and the circulation patterns more ordered in WL than in RL. The predominant sense of circulation also switched from CCW to clockwise (CW).
To explain this, we return to our minimal model, and again consider a CR-like swimmer for simplicity. We first increase the swimming speed from 80 μms-1 to 100 μms-1, to reflect the increase in speed in WL at early times compared to RL. As before, we hypothesize that biflagellar synchrony may be key. High-magnification imaging of CR swimming in WL revealed a dramatic improvement in flagellar synchrony, consistent with previous observations (Wan et al., 2014). Accordingly, we update the offset angle from in RL to in WL. These values are illustrative, as there is no clear mapping between flagellar synchrony and shape offset. With these two simple modifications, the model reproduces both the increase in flux strength and reversal in flux direction in WL (Figure 4F and G). The simulations can be extended to account for more complex motility patterns (e.g. episodic swimming in PO), or a variable swimming speed.
Behaviour is compressed into three motility macrostates
How must the dynamics of individual locomotor appendages change, to explain the above whole-cell behaviours? How does this depend on the species? To answer these questions, we use long-time statistics to provide a robust measure of behaviour that is independent of a cell’s physical environment. We first observe that the motility space is low-dimensional, comprising only a trio of movement states (Figure 5 and Appendix 1). Both species exhibit a quiescent or ‘stop’ state in which the cell body and cilia exhibits minimal movement, a ‘run’ state associated with constant speed and smooth forward swimming, and a ’transitional’ state for re-orientations. This is analogous to prokaryotic strategies such as run-and-tumble in E. coli (Berg, 2003; Perez Ipiña et al., 2019), or the run-reverse-flick in V. alginolyticus (Son et al., 2015).
We can decompose long-time tracks into discrete timeseries of these states (labelled ) according to a heuristic approach based on direct observations of ciliary beat patterns and associated trajectory parameters such as linear or angular velocity . At time , the cell is in state . For CR, we define ’runs’ () to be when the two cilia engage in synchronous breaststrokes (high , low ), ‘tumbles’ () when the cilia lose synchrony (low v, high ), and ’stops’ () where for a threshold speed (Figure 5A–D). Similarly, the behaviour of PO can also be mapped to a tripartite run-stop-shock repertoire (Wan and Goldstein, 2018), where PO ’shocks’ are analogous with CR ’tumbles’. ’Runs’ () result from coordinated breaststrokes involving 8-cilia (Wan, 2020), corresponding to moderate and low . ’Shocks’ () occur when all eight cilia switch to a symmetric beat and undulate in front of cell, producing rapid backward swimming followed by reorientation (high and high ). Finally, ’stops’ () are where for a threshold speed (Figure 5H–K). We note that for both species while coordinated strokes produce straight swimming across the trap interior (’runs’), curvature-guided interactions with the wall lead to circling behaviour (’runs’ along the boundary); both are considered ’runs’ here (Figure 5B and I).
This discrete state representation allows us to quantify how sub-cellular dynamics change over time or in response to environmental cues. We can estimate for each cell and for each assayed condition (see Appendix 1) the expected probability of being in state (pi), the list of sojourn times in state (), survival probabilities and expected residence times in state , as well as the pairwise state transition probability () and the transition rate (number of transitions per unit time) from state to . Together, these rate constants specify a unique network akin to a chemical reaction network. An example from the 120 μm traps is shown (Figure 5E and L). See Figure 5—figure supplement 1 for results in the 60 μm traps.
In terms of state probabilities for CR in the 120 μm traps (Figure 5G): (runs more likely in WL), (tumbles less likely in WL), but (both equally likely to stop swimming by 1 hr). WL also increased the tumble→run transition rate (q21) and decreased the run→tumble transition rate (q12) (Figure 5G), Empirically, this means that CR performs more frequent and longer runs in WL compared to RL. For PO, the difference between RL and WL is more obvious, see for example the run-state survival probability (Figure 5M). State probabilities remained largely constant over 1 hr, but decreased over time (runs and shocks became less likely) while increased (stop more likely) (Figure 5N). Initially , but after the 4th timepoint (about 40 min) this drops to below the corresponding RL value. We see that (run to shock transitions are much less likely in WL), and increased over time so that run→stop transitions became more likely. This means that in WL, PO runs are longer and more likely than in RL, but cells gradually lost motility over time.
Algal cell motility is light-switchable
Next we perform a dynamic assay to see how the same individual responds to photostimulation, and if such responses are reversible. Single CR or PO cells were trapped in 120 μm wells (again per trap size) and imaged continuously for 20 minutes in the following sequence: (1) 5 minutes in RL (RL1), (2) 5 min photostimulation in WL (WL1), (3) 10 min back in RL (RL2). For both species, motility patterns in RL1 and WL1 are similar to previous results under constant RL or WL illumination. As before, single-cell heterogeneity is observed (Figure 6A and B).
Both species display a reversible but asymmetric photokinetic response to WL1. Response to RL → WL (step-up) is fast and occurs within seconds, but recovery after WL → RL (step-down) is slow and persists for several minutes into the RL2 phase (Figure 6C and D). Most cells swim faster with longer runs in WL1 period than in the preceding RL1 period. In CR this is primarily by increasing the speed of prospective runs (Figure 6E), but in PO this is by reducing the likelihood of state transitions, particularly of shocks immediately in response to step-up (Figure 6D and F). The suppression of shocks in PO is consistent with previous results in constant WL. Cells gradually adapt by increasing the transition frequency over the remainder of WL1, and well into RL2.
Thus, algal swimming motility is light-switchable in the presence of short-lived stimuli, which induce perturbations to the cell’s internal motility state network. Cells retained a memory of the stimulation (WL1) for longer than the stimulus duration itself. Analysing ‘run’ phases exclusively to extract the most prominent frequencies from the helical trajectories (see Appendix 1), we further deduce that CR swims faster by increasing the frequency and synchronicity of ciliary beating (from to Hz), in response to step-up from RL1 to WL1 (Figure 6G). A similar analysis of PO trajectories yielded significant cell-to-cell variability but no clear frequency signatures that could be correlated with ciliary beating (Figure 6H). This is likely due to the unique coordination of the 8 cilia that suppresses the helicity of trajectories (Wan, 2020).
Effect of chemical stimulation
In many protists, ciliary motility is modulated by chemical sigalling, including Ca2+-influx through ion channels, suggesting a conserved Ca-dependent mechanism (Inaba, 2015; Wan and Jékely, 2021). Ionic fluxes elicit transitions in cilia beating, waveform and swimming behaviour (Beck and Uhl, 1994; Kung and Naito, 1973; Brette, 2021; Geyer et al., 2022). In Paramecium, Ca2+, and K+ ions are major regulators of motility, including transitions between forward and backward swimming (Kunita et al., 2014). We hypothesize that motile algae also possess a behavioural response to changes in external ion concentrations, and confirmed this in PO with a bulk motility assay (see Materials and methods). The introduction of a droplet pre-loaded with 50 mM KCl to a suspension of PO cells produces perturbed swimming (Figure 7A).
Traditional, bulk assays do not allow path-sampling of single-cell trajectories, nor do they have the ability to track in real-time how a single cell responds to chemical perturbations. We solve this by designing a microfluidic device to deliver controlled fusion of two droplets (Figure 7B). Using a cross-junction microfluidic chip, we generated droplets in alternation that contain either a 10 mM KCl solution or cell suspension. Pairs of droplets were then trapped in downstream doublet wells. We then identified droplet pairs comprising one with a trapped PO cell, and another containing the 10 mM KCl solution spiked with 1 µM of fluorescein which acted as a reporter. Fusion is triggered by starting the flow of 40% perfluoro-1-octanol (PFO) in HFE (a competing surfactant); see Materials and methods. Two such droplets fuse within 1 or 2 ms. Note that trap-occupancy is stochastic: not all droplet pairs have the target composition, since droplet generation does not alternate perfectly between the two channels. Doublet wells with both droplets containing KCl, or both droplets containing cell suspension, were not imaged (Figure 7B, inset).
We compare the motility of cells in these doublet traps in the presence or absence of the chemical perturbation. Control experiments were performed by substituting the KCl medium with normal culture medium (without cells). Prior to fusion, cell motility was similar in both cases (Figure 7C). After fusion, cells paired with a KCl droplet immediately showed erratic swimming, followed by gradual decay in motility over time (see Video 5), which is likely a combined effect of an increase in K+ and decrease in Ca2+ concentration (Table 1). Significant changes in the morphology and tortuosity of swimming trajectories are also observed after fusion (Figure 7D and E). Further examples of this response are presented in Figure 7—figure supplement 1.
Discussion
Stereotypy and the ‘arrow of time’ in long-term single cell behaviour
The emerging field of computational ethology seeks to track and annotate the complex behaviours of organisms (Han et al., 2018; Yaski et al., 2011; Yemini et al., 2013; Stephens et al., 2011). Due to their small size, variable morphology, and high-speed movement, detailed analysis of the behavioural patterns of microorganisms is technically challenging. Here, we leveraged droplet microfluidics to track the swimming trajectories of microbes for long periods of time. Comparing a freshwater biflagellate with a marine octoflagellate, we quantified how single-cells respond to a number of controlled environmental stimuli (including mechanical, light and chemical).
We showed that the behavioural space of roaming microbes comprise highly stereotyped movements, and proposed a new paradigm for phenotyping microbial motility in terms of a minimal set of three motility macrostates (runs, stops, tumbles/shocks). This inherent low-dimensionality has been observed in other organisms (Larson et al., 2022; Jordan et al., 2013; Tsang et al., 2018) and artificial microswimmers (Leoni et al., 2020; Hokmabad et al., 2021). Over long times, CR swimming is fairly continuous, yet PO transitions between bursts of fast swimming and lengthy periods of quiescence, reminiscent of ’wakeful’ versus ’sleep states’ observed in animals (Lawler et al., 2021).
Importantly, our long-term behavioural recordings retained the ’arrow of time’ (Gnesotto et al., 2018). This gives us access to key information that may be hidden by traditional, bulk-averaged measurements, such as chiral movements in confined spaces (see next section). From these ultralong timeseries, we quantified how behaviour changes at the sub-cellular level. Transition rates and probabilities between states were shown to not only change in time but also in response to environmental stimuli. Such timeseries can also be mined to infer non-equilibrium entropy production (Skinner and Dunkel, 2021; Wan and Goldstein, 2018). These species-specific movement signatures reflect different intracellular processes for controlling the algal motility apparatus (Guo et al., 2021), and could be key to the emergence of cognition and environmental responsiveness (Kunita et al., 2016; Wan and Jékely, 2021). Microbes inhabiting different ecological niches (e.g. freshwater or marine) likely evolved divergent strategies for sensorimotor control, despite all relying on the same organelle for motility. Future work should seek to map the space of ciliary waveforms, itself low-dimensional (Geyer et al., 2022), to whole-cell swimming behaviours.
Flux loops without curvature gradients
Microbes often interact with their physical environment, particularly with interfaces and boundaries. Interactions are commonplace in natural settings with exposed surface features or other heterogeneous structures, for example soil, foam, or particulate matter (Théry et al., 2021; Kantsler et al., 2013; Souzy et al., 2022). Here, we engineered PDMS chambers with precise shapes to explore how confinement affects cell motility. In red light, cells tolerated long-periods of confinement with little change to their overall motility. We also discovered a chiral circling behaviour, confirmed by flux representations of the time-averaged trajectories. While confinement has been shown previously to stabilise macroscopic chiral movement in swimmer suspensions (Wioland et al., 2013; Beppu et al., 2021), here chirality has emerged in the travel history of a single cell. This implicates a distinct route towards a chiral non-equilibrium steady-state that has not been reported previously.
That the sense of chirality also depends on light is strong indication that it arises from a mechanism under cellular control. We outlined a model for how this can occur, based on a minimal representation of a CR-like cell as a dumbbell with fore-aft asymmetry (Cammann et al., 2021). Such a swimmer will align with solid boundaries in a curvature-dependent manner, as shown previously (Ostapenko et al., 2018). In a circular trap (constant curvature), it should not produce flux loops on average (Cammann et al., 2021). We find that breaking of the dumbbell swimmer’s left-right symmetry was sufficient to produce trajectories with CW or CCW bias, depending on the sign of the offset . This is supported by data showing the bilateral flagellar dominance (asymmetry between cis and trans flagella) to be tunable by light. This is also important for helical swimming and klinotaxis (Cortese and Wan, 2021; Leptos et al., 2022). Thus, even minor cell-internal asymmetries (here represented by a very small ) are sufficient to break macroscopic symmetry. Further particulars of cell shape (Lushi et al., 2017; Zaferani et al., 2021), wall interactions (Lauga et al., 2006; Spagnolie et al., 2017), or cilia mechanosensitivity (Wan and Goldstein, 2018) may be needed to account for the more complex swimming repertoire of PO.
Light-dependent algal motility and phototoxicity
We also found that motility parameters in RL are more stable than in WL, for both species. Short periods of WL exposure increases the average run-speed of CR by increasing ciliary beat frequency, as found previously in micropipette-immobilised cells (Rüffer and Nultsch, 1990). In contrast, WL decreases the frequency of gait transitions in PO, producing longer stops, longer runs, and fewer shocks. Prolonged or excessive exposure (> 40 minutes) to WL eventually led to reduced or cessation of motility in both species. In CR this is likely due to WL-induced adhesion (Kreis et al., 2018). In PO, there is more evidence of irreversible photodamage. Increasing confinement (smaller traps) decreased the time taken to reach this state, suggesting that motility deterioration may be due to buildup of some metabolite or excretion from the trapped cell over time (Boedicker et al., 2009).
Shorter periods ( min) of WL-stimulation produced reversible changes in motility in the algae (Figure 6), again with species-dependent photokinetic responses occurring at the single cell level. While CR responded largely by modulating run speed, PO modulated the balance of motility macrostates (e.g. suppressed shocks). These distinct behaviours offer intriguing prospects for synthetic biology and bioengineering, such as light-guided patterning of microbes (Frangipane et al., 2018; Bittermann et al., 2021). It will be interesting to combine our assay with genetic perturbations (particularly in CR) to reveal the regulatory pathways responsible for photosynthesis, phototransduction, and phototaxis.
Future prospects of droplet microfluidics for assaying cell motility
Micro-encapsulation provides unique functionalities that enable controlled, on-demand creation and manipulation of a cell’s local environment (Son et al., 2015). A microfluidics-based pipeline allows us to stably trap motile cells, to establish their unique behavioural signatures. We use multilayer devices and static droplet trapping wells (Fradet et al., 2011) to keep droplets stable for hours, and isolated from possible fluctuations in the continuous phase (Han et al., 2020). This allows us to achieve long-time imaging, and reliable control over trap size, shape, and geometry. The trapping arrays we used maximise the probability of obtaining multiple singly-trapped cells in a single experiment (only limited by Poisson statistics), thus increasing analytical throughput. Our current design focuses on imaging a small number of traps at a time. Simultaneous, high-throughput imaging of multiple cells in multiple wells is also possible, but at the detriment of spatial resolution and thus may not be sufficient for identifying motility state or behavioural transitions.
Our work paves the way for new applications and opportunities for designing diagnostic tools in the absence of molecular tests. Microfluidics-based assays enable the detection of single-cell motility signatures and any heterogeneity in microbial behaviour. Phenotypic diversity is critical in microbial ecosystems, where they may give rise to distinct selection pressures and antimicrobial resistance (Dhar and McKinney, 2007). Our droplet fusion device represents a highly novel solution for assaying fast cell-environment interactions and responses to chemical perturbations (Figure 7). Contrary to traditional assays that use flow channels to deliver concentration gradients, our design applies real-time control and localised perturbations to single microbes. Further developments and extensions are on the way, including augmentation of on-demand single cell encapsulation with active cell sorting (Anagnostidis et al., 2020; Howell et al., 2022). The integration of lab-on-chip technologies, high-speed microscopy and computer vision has significant potential for reconstructing the species-specific sensorimotor pathways of microorganisms, and revealing their response thresholds to dynamic environmental perturbations.
Materials and methods
Cell culturing and maintenance
Request a detailed protocolC. reinhardtii cultures were prepared from axenic plates of the CC125 WT strain (Chlamydomonas Center). Individual colonies were transferred to volumes of Tris-minimal media. Liquid cultures were grown under a 14/10 light-dark cycle at 21°C and 40% humidity, with constant shaking at . Liquid cultures were sub-cultured when they approached the end of the exponential phase. For motility experiments, second- or third- generation cultures were harvested in late-exponential phase (6-9 days after inoculation). Cell density was measured as ∼1 × 106 cells per mL. Cells were centrifuged at 100 g for 10 min, and then concentrated 10-fold. The cells were then left in darkness for a minimum of 30 min to dark-adapt them before each experimental run.
We prepared our PO cultures from axenic liquid cultures of the WT of the species P. octopus (NIVA/NORCCA). 200 µL of axenic culture was transferred to of TL30 media. Cells were grown under continuous illumination at 21°C and 40% humidity, without shaking. For motility experiments, cultures were harvested during the latter half of the exponential phase (20-30 days after inoculation). Cell density was measured at ∼3 × 103 cells per mL. Cells were centrifuged at 100 g for 10 min, and then concentrated 10-fold. The cells were then left in darkness for a minimum of 30 min to dark-adapt them before each experimental run.
Microfluidic chip fabrication
Request a detailed protocolThe devices were designed with CAD software (DraftSight, Dassault Systems) and fabricated following classical soft-lithography procedures by using a high-resolution acetate mask (Microlithography Services Ltd.). Negative photoresist SU-8 3025 (MicroChem, Newton, MA) was spin-coated onto clean silicon wafers to a thickness of 10 μm, patterned by exposure to UV light through the photomask and hard bake at 95ºC for ∼7 minutes (Xia and Whitesides, 1998). Prior to development through immersion in propylene glycol monomethyl ether acetate (PGMEA, Sigma-Aldrich), a second layer of SU-8 3025 at 20 μm in height was spin-coated, UV exposed and hard baked (95ºC, ∼7 minutes) for the development of the trapping arrays. Uncured polydimethylsiloxane (PDMS) consisting of a 10:1 polymer to cross-linker mixture (Sylgard 184) was poured onto the master, degassed, and baked at 70ºC for 4 hours. The PDMS mould was then cut and peeled from the master, punched with a biopsy punch (Kai Medical) to create inlet ports for tubing insertion. A total of three holes were punched; two inlets for the continuous and aqueous phase and an outlet for waste collection. The PDMS mould was plasma bonded to thin cover slips ( mm, mm thick). Hydrophobic surface treatment was performed immediately after bonding by flushing with 1% (v/v) Trichloro (1 H, 1 H, 2 H, 2H-perfluorooctyl) silane (Aldrich) in HFE-7500, and placed in a 65ºC oven for 30 min.
Flow-focusing droplet generation
Request a detailed protocolMicrofluidic device fabrication was done using classical soft lithography techniques. A total of 4 devices were developed. All devices consisted of a flow-focusing junction for droplet generation (height, 10 μm), and a second layer with a trapping array made up of circular wells (height, 20 μm)(Figure 1A). The trap diameters were 40, 60, 120 and 200 μm. The dimensions of the flow-focusing junction varied depending on trap size. We designed a range of dimensions for the flow focusing junctions and matching trap sizes (i.e. circle diameter). Depending on the trapping size, the total number of traps was between 78 (∅ = 200 μm) and 840 (∅ = 40 μm).
Droplets were typically generated at rates approx ∼ 50 s-1. The flow rates were controlled using syringe pumps (Nemesys, Cetoni), 1 ML plastic syringes (BD PlastipakTM; sterile needles, 25G x 1” – NR. 18, 0.5 × 25 mm, BD MicrolanceTM 3), and portex tubing PE (Scientific Laboratory Supplies, 0.38 × 0.355 mm). The flow rates for oil and cell suspension were varied depending the size dimensions of each device. A 1:3 ratio was aimed for the continuous and aqueous phase, respectively.
The carrier oil phase was prepared using fluorinated oil HFE-7500 (Fluorochem Ltd) containing 0.25% (w/v) 008-Fluorosurfactant (RAN Biotechnologies, Inc). The aqueous phase consisted of liquid cell cultures (see Cell culturing and maintenance). The droplets were generated at the flow-focusing junction creating water-in-oil emulsions. A dilute suspension of algae was injected through the inlet. Following trapping, droplets were stably confined to the microwell during imaging acquisition. Once all the traps were filled, the aqueous phase flow was halted and the continuous phase flow ratewas reduced approximately 5-fold to flush away excess droplets.
Live-cell high-speed imaging
Request a detailed protocolBrightfield imaging was conducted with an inverted microscope (Leica Microsystems, DMi8), with a high-speed camera (Phantom Vision Research, V1212). We first scanned the array of trapped cells to locate traps matching our criteria (droplet fitting exactly into the trap, droplet containing only one cell). For the 40 μm, 60 μm and 120 μm trap sizes, as well as the cell fusion experiments, we used a 20 x long-working distance objective (HC PL/0.40). For the largest 200 μm traps, we lowered the magnification to 5 x (NPLAN/0.12) equipped with a 1.6 x tube lens, to reduce file size. All traps were imaged with the same intensity and aperture settings, and at 500 fps. For the 1 hour confinement experiments, cells were imaged continuously but 5-min recordings taken at 5 min intervals, to obtain a total of 6 timepoints per cell. Data from droplets that were disrupted at any point during imaging was discarded.
Light-modulation experiments
Request a detailed protocolFor brightfield imaging in WL, we used a standard broad-spectrum LED source to illuminate the specimen. Red light (RL) imaging was accomplished by insertion of an IR long-pass filter (610 nm, Chroma) to the light path. Spectra corresponding to the two possible illumination options are compared in Figure 8. For light-switching experiments, the red filter was removed or inserted manually.
Bulk cell motility assay
Request a detailed protocolWe assayed the effect of KCl on P. octopus behaviour via a simple open-air method. We first added 1 µL of a concentrated suspension of cells to a glass coverslip under red light in a dark room. We waited for flows from the placement of the droplet to subside, before imaging for (at 10 x, 100 fps). We then added either a 1 µL droplet of culture media (for the control) or a 1 µL drop of KCl (for the KCl test), and after waiting another for flows to subside we imaged for another .
Paired-droplet fusion assay
Request a detailed protocolDroplet pairs were generated using a cross junction microfluidic device in which cells and KCl solutions flowing in separate channels are encapsulated in alternation. We used a 25 degrees taper angle previously reported to produce the most stable alternation function (Saqib et al., 2018). When stable alternation was achieved, droplets were suddenly halted by removal of the inlet tubings of KCl and cell solutions followed by gentle removal of excess droplets by flowing oil at 3 µL/min for ∼1 min. We subsequently identified droplet pairs of expected volumes with one containing a single cell. To ensure fusion with KCl and the absence of mixing prior to fusion, the KCl solution was spiked with 1 µM fluorescein (Merck) which was imaged before triggering of fusion. A 15-min video of the cells was acquired displaying the entrapped droplets 7.5 min prior and 7.5 min post fusion. Fusion was induced by surfactant replacement with 1 H,1H,2H,2H-Perfluoro-1-octanol (Merck) (PFO). A solution of 40% PFO in HFE was flown at 5 µL/min and run until fusion was achieved (∼ 4 minutes). PFO competes with the other surfactant which destabilizes the droplet interface to induce rapid, reproducible fusion. The change in ion concentrations in the solution containing the cell pre- and post-fusion can be found in (Table 1). Control experiments where the KCl solution substituted with normal culture medium were performed to confirm the absence of confounding factors.
For the fusion experiments, the identification process consisted of two steps. Firstly we located a trap that had two equal droplets in place, one containing a cell. Secondly, we took a fluorescence image to verify that the droplet without the cell contained KCl. The fluorescence image was taken in the LASX software, using a broad-spectrum LED source (CoolLED-pE300) equipped with a triple-band filter set (including FITC, Ex: , Em: ). The fluorescence intensity was set to 60%, the exposure time was , and the gain was 2.0. The presence of fluorescence covering the whole of the empty droplet was sufficient to prove the presence of KCl.
Image processing and cell tracking
Request a detailed protocolRaw video data was exported to 8-bit grayscale and enhanced by subtracting an average image in MATLAB (Mathworks). Trap boundaries were identified manually to increase the fidelity of 2D cell tracking, which was performed automatically using the TrackMate plugin in ImageJ (Tinevez et al., 2017). A Laplace of Gaussian detector was used for spot identification, with slightly different blob diameters for CR and for PO (14 μm and 21 μm respectively). Single continuous tracks were obtained for each experimental run ( individuals per condition), and exported for further processing and extraction of detailed track features/other statistics (see Appendix 1). Video frames from the bulk motility assays were processed and analysed similarly.
Appendix 1
Analysis of trajectories
Trajectories obtained from cell tracking (see Materials and methods) consist of 2D coordinates at discrete times (t1, t2, t3…) with a constant time interval, i.e. for videos captured at 500 fps. The trajectories were analysed using custom MATLAB scripts. First, a smoothing filter was applied. For CR, a Savitzky–Golay filter with order 2 and frame length 201 was used to smooth out the helical trajectories such that a speed and angular velocity corresponding to the net forward motion of the cell could be obtained. For PO, a Savitzky–Golay filter with order 2 and frame length 21 was used. A smaller frame length was used due to the fast timescale of the shock behaviour and since PO has smoother trajectories than CR (Appendix 1—figure 1). The cell velocity at time ti was calculated as
where is the 2D coordinate of the centroid of the cell at time ti for the smoothed trajectories. The angular velocity of the cell was defined as
where is the normalised velocity vector at time ti. To reduce the noise of the angular velocity data, the results reported are a moving mean across 25 frames.
Violin plots, which combine box plot and histogram data representations into one diagram, were created using a MATLAB package (Bechtold et al., 2021). The width of the violin plots correspond to the probability density function and are scaled linearly for CR and logarithmically for PO.
Radial probability density
The radial probability density is defined as
where is the cell’s distance from the centre of the trap and is the count of the number of trajectory points that lie in a circular shell at distance with width .
Mean square displacement
The the mean square displacement (MSD) is defined as
where is the delay and the number of pairs of time points in the trajectory with . The MSD was calculated for the raw 2D coordinate positions using the ’@msdanalyzer’ MATLAB package (Tarantino et al., 2014). The MSD results reported here were calculated using 1/50th of the data points (i.e. every 50th frame) due to the data array size constraints of the package.
The actual trap diameters of the droplets were variable because the resolution of the film mask used to fabricate the PDMS devices was 4 µm. The relative error becomes more noticeable for smaller trap sizes and explains the discrepancy between the experimental and simulation results (Figure 3H).
Rotational diffusion
The effective rotational diffusion coefficient was obtained using the mean square angular displacement (MSAD) for CR in the largest trap size (i.e. ∅ = 200 μm). First an angle defining the orientation of the cell relative to the y-axis is computed as
where is a unit vector defining the orientation of the cell at time and, assuming that the cell always swims forwards, is given by , where is the velocity vector. The MSAD is then defined as
where is the delay and the number of pairs of time points in the trajectory with . The data is then fitted to for in order to obtain an estimate for the effective rotational diffusion coefficient (Appendix 1—figure 2).
Probability flux calculation and relative probability density
We calculated the probability fluxes using a method previously applied to trajectories of C. reinhardtii (Cammann et al., 2021) and first introduced by Battle et al., 2016. The 2D positional space is divided into a grid of equally sized square boxes with side length . For each time point the position is assigned to a box , where and are the box positions in the and direction respectively. From these coarse-grained trajectories, a time-series of transitions is obtained by constructing the following array
where and are the positions of consecutively visited boxes and is the length of time spent in the initial state before transitioning to the new state . In a small number of cases, the two success states and do not correspond to nearest neighbours. In such cases, the intermediate boxes are determined by linear interpolation and the corresponding extra transitions are inserted into the array (Equation 7) to ensure that all transitions in are between nearest neighbours.
The net transition rates between neighbouring boxes are calculated from the coarse-grained trajectories using
where is the number of transitions from box to box and is the total duration of the trajectory. The net transition rates are then used to calculate the probability flux
This method is summarised in Appendix 1—figure 3. Note that this calculation does not account for diagonal transitions, however since they account for only < 5% of all nearest neighbour transitions, we assume they have a minimal effect on the results.
To obtain the mean tangential flux plotted in Figures 3I and 4G, the tangential component of the flux vector was obtained for each box and the mean calculated for all boxes with the same radial distance from the centre of the trap.
The 2D trajectory data was also used to calculate the relative probability density
where is the number of trajectory points within box , is the box area and is the area of the trap (Ostapenko et al., 2018).
Conversion into three motility macrostates
We started from raw track data, and assigned states in CR and PO using a combination of linear and/or angular Ω speeds. In both cases, we took moving averages to reduce frame-to-frame noise (due to the high imaging frame rates).
States for Chlamydomonas reinhardtii: (run, stop, tumble)
We defined the ‘stop’ state as times when 5 μms-1. A smoothing filter was applied to the binary stop/move data to remove spurious state transitions. We then verified the stop states by visual inspection. To identify ’tumbles’, we use angular velocity to locate times when . We filtered the dataset to correct for false tumbles by removing tumbles that were under , and for false runs of very short duration found between successive tumbles. We again verified that tumble states corresponded with cell reorientation by visual inspection. Finally, frames that were neither a ’stop’ nor a ’tumble’, were designated ’runs’, with ’stops’ taking precedence over ’tumbles’. The workflow is summarised in Appendix 1—figure 4.
States for Pyramimonas octopus: (run, stop, shock)
Here, only the linear speed is sufficient to assign states to PO, since they are associated with very distinct speeds (Wan and Goldstein, 2018). We defined the ‘stop’ state as times when the smoothed speed 5 μms-1. A smoothing filter was applied to the binary stop/move data to remove spurious state transitions. We then verified the stop states by visual inspection.
To identify shocks, we first smoothed the data by computing a local median value and then a moving mean. We then identified each local minimum and the subsequent local maximum, and computed the increase between the two. Increases in median speed greater than 185 μms-1 were identified as shocks. The start point of shocks were chosen as the point in time between the local minimum and maximum where half of the total increase in displacement had occurred. Equivalently, the end of the shock was defined as the point in time between the current displacement maximum and post-shock minimum where half of the total decrease in displacement had occurred. Finally, frames that were neither a ’stop’ or a ’shock’, were designated ’runs’, with ’stops’ taking precedence over ’shocks’. Brief (<) ’run’ states between a ’stop’ and a ’shock’ were reclassified as ’stop’ states to remove spurious state transitions. The workflow is summarised in Appendix 1—figure 5.
State probability and transition rate analysis
We estimated transition probabilities between motility macrostates via a simple counting algorithm. State probability is given by
where ni is the number of frames in which the cell is classified to be in state .
Transition probability from state to is given by
where is the number of transitions from state to state .
Transition rate from state to is defined as
where is the mean duration of state .
The survival probabilities for state are defined by
where is the number of instances where the duration of state is longer than .
Beat frequency during a run
The cilia beat frequency for each run period longer than was estimated using a fast Fourier transform analysis of the raw speed (i.e. calculated using the raw centroid positions). A second order Savitzky–Golay filter with frame length 45 was used to reduce the noise of the Fourier transform and the beat frequency for each run was taken to be the highest peak within 30-90 Hz (or if the highest peak was twice the frequency of the second highest peak, then the latter was taken as the beat frequency).
Numerical simulation of a microswimmer in confinement
The CR cells are modelled as asymmetric dumbbells with a large sphere in front and a smaller sphere in the back, which capture in a minimal way the fore-aft asymmetry of body and appendages (Roberts and Deacon, 2002; Roberts, 2006, see Appendix 1—figure 6). The dumbbell’s spheres have different radii and are separated by a distance . Dumbbell models have been used previously to model such microswimmers (Ledesma-Aguilar and Yeomans, 2013; Wensink et al., 2014; Wysocki et al., 2015; Ostapenko et al., 2018; Cammann et al., 2021). We improve the model by taking into account the slightly different beating patterns of the cis and trans flagella (Cortese and Wan, 2021). The asymmetric beating pattern gives rise to an effective chirality in the direction of motion. We implement this by introducing an angle between the swimming direction and the axis of symmetry of the dumbbell. The equation of motion for the position of the active dumbbell is given by
Here, and is the self-propulsion speed of the cell in red and white light, respectively, and is a Gaussian white noise with correlator with translational diffusion coefficient (both values are taken from Ostapenko et al., 2018). The term accounts for steric wall interactions, where , labels the large and small sphere of the dumbbell, respectively. To compute the steric forces we use the following potential
if , and 0 otherwise, where is the distance of the sphere to the wall. The radii of the dumbbell’s spheres are , (see Ostapenko et al., 2018). The choice of a softer potential in Equation 16 is dictated by the rather broad peak in radial probability density compared to the results in Ostapenko et al., 2018. The direction of motion is defined with respect to the vector pointing from the small to the large sphere as where . We find that values of (RL) and (WL) model the behavior observed in the different lighting conditions well. The equation of motion for reads
Here is a Gaussian white noise with correlator and rotational diffusion coefficient (see Appendix 1—figure 2). The torque acting at the wall , where , , and . For the shear time at the wall we use (see Kantsler et al., 2013).
Data availability
New data and analysis codes generated as part of this study are available for download from Zenodo. The dataset includes all raw cell trajectories and motility states, as well as analysis and simulation codes.
-
ZenodoDataset for: Phenotyping single-cell motility in microfluidic confinement.https://doi.org/10.5281/zenodo.7226288
References
-
On the localization of voltage-sensitive calcium channels in the flagella of Chlamydomonas reinhardtiiThe Journal of Cell Biology 125:1119–1125.https://doi.org/10.1083/jcb.125.5.1119
-
The rotary motor of bacterial flagellaAnnual Review of Biochemistry 72:19–54.https://doi.org/10.1146/annurev.biochem.72.121801.161737
-
Bacterial hopping and trapping in porous mediaNature Communications 10:2075.https://doi.org/10.1038/s41467-019-10115-1
-
Integrative neuroscience of paramecium, a “swimming neuron.”ENeuro 8:ENEURO.0018-21.2021.https://doi.org/10.1523/ENEURO.0018-21.2021
-
Control of helical navigation by three-dimensional flagellar beatingPhysical Review Letters 126:088003.https://doi.org/10.1103/PhysRevLett.126.088003
-
Microbial phenotypic heterogeneity and antibiotic toleranceCurrent Opinion in Microbiology 10:30–38.https://doi.org/10.1016/j.mib.2006.12.007
-
Broken detailed balance and non-equilibrium dynamics in living systems: a reviewReports on Progress in Physics. Physical Society 81:066601.https://doi.org/10.1088/1361-6633/aab3ed
-
Intracellular coupling modulates biflagellar synchronyJournal of the Royal Society, Interface 18:20200660.https://doi.org/10.1098/rsif.2020.0660
-
Long-Term observation of Magnetospirillum gryphiswaldense in a microfluidic channelArchives of Microbiology 201:1427–1433.https://doi.org/10.1007/s00203-019-01713-0
-
Static array of droplets and on-demand recovery for biological assaysBiomicrofluidics 14:051302.https://doi.org/10.1063/5.0022383
-
Emergence of bimodal motility in active dropletsPhysical Review X 11:e43.https://doi.org/10.1103/PhysRevX.11.011043
-
Multi‐Object detector yolov4‐tiny enables high‐throughput combinatorial and spatially‐resolved sorting of cells in microdropletsAdvanced Materials Technologies 7:e2101053.https://doi.org/10.1002/admt.202101053
-
On the significance of the spiral swimming of organismsThe American Naturalist 35:369–378.https://doi.org/10.1086/277922
-
A ciliate memorizes the geometry of A swimming arenaJournal of the Royal Society, Interface 13:118.https://doi.org/10.1098/rsif.2016.0155
-
Swimming in circles: motion of bacteria near solid boundariesBiophysical Journal 90:400–412.https://doi.org/10.1529/biophysj.105.069401
-
Sleep analysis in adult C. elegans reveals state-dependent alteration of neural and behavioral responsesThe Journal of Neuroscience 41:1892–1907.https://doi.org/10.1523/JNEUROSCI.1701-20.2020
-
Enhanced motility of a microswimmer in rigid and elastic confinementPhysical Review Letters 111:13.https://doi.org/10.1103/PhysRevLett.111.138101
-
Scattering of biflagellate microswimmers from surfacesPhysical Review. E 96:023102.https://doi.org/10.1103/PhysRevE.96.023102
-
DeepLabCut: markerless pose estimation of user-defined body parts with deep learningNature Neuroscience 21:1281–1289.https://doi.org/10.1038/s41593-018-0209-y
-
Electro-coalescence of digitally controlled dropletsAnalytical Chemistry 81:7321–7325.https://doi.org/10.1021/ac901188n
-
Curvature-guided motility of microalgae in geometric confinementPhysical Review Letters 120:068002.https://doi.org/10.1103/PhysRevLett.120.068002
-
Bacteria display optimal transport near surfacesNature Physics 15:610–615.https://doi.org/10.1038/s41567-019-0460-5
-
Gravitaxis in motile micro-organisms: the role of fore–aft body asymmetryJournal of Fluid Mechanics 452:405–423.https://doi.org/10.1017/S0022112001006772
-
Mechanisms of gravitaxis in ChlamydomonasThe Biological Bulletin 210:78–80.https://doi.org/10.2307/4134597
-
Live from under the lens: exploring microbial motility with dynamic imaging and microfluidicsNature Reviews. Microbiology 13:761–775.https://doi.org/10.1038/nrmicro3567
-
Microbial narrow-escape is facilitated by wall interactionsPhysical Review Research 4:2029.https://doi.org/10.1103/PhysRevResearch.4.L022029
-
Tnf and IL-1 exhibit distinct ubiquitin requirements for inducing NEMO-IKK supramolecular structuresThe Journal of Cell Biology 204:231–245.https://doi.org/10.1083/jcb.201307172
-
Lag, lock, sync, slip: the many “ phases” of coupled flagellaJournal of the Royal Society, Interface 11:20131160.https://doi.org/10.1098/rsif.2013.1160
-
Time irreversibility and criticality in the motility of a flagellate microorganismPhysical Review Letters 121:058103.https://doi.org/10.1103/PhysRevLett.121.058103
-
Synchrony and symmetry-breaking in active flagellar coordinationPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 375:1792.https://doi.org/10.1098/rstb.2019.0393
-
Origins of eukaryotic excitabilityPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 376:20190758.https://doi.org/10.1098/rstb.2019.0758
-
Controlling active self-assembly through broken particle-shape symmetryPhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 89:010302.https://doi.org/10.1103/PhysRevE.89.010302
-
Confinement stabilizes a bacterial suspension into a spiral vortexPhysical Review Letters 110:268102.https://doi.org/10.1103/PhysRevLett.110.268102
-
Giant adsorption of microswimmers: duality of shape asymmetry and wall curvaturePhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 91:050302.https://doi.org/10.1103/PhysRevE.91.050302
-
Arena geometry and path shape: when rats travel in straight or in circuitous paths?Behavioural Brain Research 225:449–454.https://doi.org/10.1016/j.bbr.2011.07.055
-
A database of Caenorhabditis elegans behavioral phenotypesNature Methods 10:877–879.https://doi.org/10.1038/nmeth.2560
Article and author information
Author details
Funding
European Research Council (853560)
- Kirsty Y Wan
Academy of Medical Sciences (SBF003\1160)
- Kirsty Y Wan
Biotechnology and Biological Sciences Research Council (BB/T011777/1)
- Fabrice Gielen
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 853560 EvoMotion, to KYW), and a Springboard Award from the Academy of Medical Sciences and Global Challenges Research Fund (SBF003\1160 to KYW). The research was also funded by the Biological and Biotechnological Research Council (grant BB/T011777/1 to FG). Calculations were performed using the Sulis Tier 2 HPC Platform funded by EPSRC Grant EP/T022108/1 and the HPC Midlands +consortium. We gratefully acknowledge the use of the Lovelace HPC service at Loughborough University.
Copyright
© 2022, Bentley, Laeverenz-Schlogelhofer et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,044
- views
-
- 275
- downloads
-
- 21
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Physics of Living Systems
We propose the Self Returning Excluded Volume (SR-EV) model for the structure of chromatin based on stochastic rules and physical interactions. The SR-EV rules of return generate conformationally defined domains observed by single-cell imaging techniques. From nucleosome to chromosome scales, the model captures the overall chromatin organization as a corrugated system, with dense and dilute regions alternating in a manner that resembles the mixing of two disordered bi-continuous phases. This particular organizational topology is a consequence of the multiplicity of interactions and processes occurring in the nuclei, and mimicked by the proposed return rules. Single configuration properties and ensemble averages show a robust agreement between theoretical and experimental results including chromatin volume concentration, contact probability, packing domain identification and size characterization, and packing scaling behavior. Model and experimental results suggest that there is an inherent chromatin organization regardless of the cell character and resistant to an external forcing such as RAD21 degradation.
-
- Physics of Living Systems
A hallmark of biomolecular condensates formed via liquid-liquid phase separation is that they dynamically exchange material with their surroundings, and this process can be crucial to condensate function. Intuitively, the rate of exchange can be limited by the flux from the dilute phase or by the mixing speed in the dense phase. Surprisingly, a recent experiment suggests that exchange can also be limited by the dynamics at the droplet interface, implying the existence of an ‘interface resistance’. Here, we first derive an analytical expression for the timescale of condensate material exchange, which clearly conveys the physical factors controlling exchange dynamics. We then utilize sticker-spacer polymer models to show that interface resistance can arise when incident molecules transiently touch the interface without entering the dense phase, i.e., the molecules ‘bounce’ from the interface. Our work provides insight into condensate exchange dynamics, with implications for both natural and synthetic systems.