Phenotyping single-cell motility in microfluidic confinement

  1. Samuel A Bentley
  2. Hannah Laeverenz-Schlogelhofer
  3. Vasileios Anagnostidis
  4. Jan Cammann
  5. Marco G Mazza
  6. Fabrice Gielen  Is a corresponding author
  7. Kirsty Y Wan  Is a corresponding author
  1. Living Systems Institute, University of Exeter, United Kingdom
  2. Mathematics and Statistics, University of Exeter, United Kingdom
  3. Biosciences, University of Exeter, United Kingdom
  4. Physics and Astronomy, University of Exeter, United Kingdom
  5. Interdisciplinary Centre for Mathematical Modelling and Department of Mathematical Sciences, Loughborough University, United Kingdom
  6. Max Planck Institute for Dynamics and Self-Organization (MPIDS), Germany
14 figures, 5 videos, 1 table and 1 additional file

Figures

Schematic of the experimental set-up.

(A) A two-layer microfluidic device with embedded single-cell traps, and syringes used for perfusion of the carrier oil phase and an aqueous suspension containing live motile cells. (B) 3D rendering …

Figure 2 with 1 supplement
Behavioural ethograms of motile algae.

(A) High-speed trajectories of single C. reinhardtii cells in circular traps of varying sizes (illustrating coverage over 5 min), overlaid with a 5-second representative trajectory (colour-coded by …

Figure 2—figure supplement 1
Summary violin plots of the speed data for different trap sizes, over the 6 timepoints (respectively, 0-5, 10-15, 20-25, 30-35, 40-45, and 50-55 min).

Note that the scaling of the probability density function is linear for CR and logarithmic for PO.

Figure 3 with 2 supplements
Effect of trap size on single cell motility.

(A,B) Violin plots of pdfs of speed (v) across all cells, and all time points, for CR and PO. (Note the scaling of the pdfs is linear for CR, but logarithmic for PO.) (C,D) Time-averaged heatmaps …

Figure 3—figure supplement 1
The experimental mean and standard error (N = 5) of the radial probability densities P(r) (A) and mean-squared displacement curves (B) for each of the trap sizes, for PO.

Compared with CR, the P(r) curves are more irregular and the MSD curves show a more gradual increase toward the maximum.

Figure 3—figure supplement 2
Bivariate histograms of linear and angular speed (v,Ω) across all cells, and all time points for CR (A), and PO (B).

Bivariate histograms of speed and radial position (v,r) across all cells, and all time points, for CR (C), and PO (D).

Figure 4 with 1 supplement
Effect of white light on single cell motility.

(A,B) Heatmaps of cell swimming speed over time, for CR and PO. (C) Time-averaged heatmaps of the cells’ centroid position overlaid with arrows showing the direction and strength of the probability …

Figure 4—figure supplement 1
Violin plots of the pdfs of speed over time for all cells, for CR (A) and PO (B).

Note that the pdfs are scaled linearly for CR and logarithmically for PO.

Figure 5 with 1 supplement
The locomotor behaviour of single-cells is controlled by a trio of motility macrostates and their transition probabilities.

(A–G) For CR, distinct ciliary beat patterns and coordination states (A) produce different swimming trajectories (B), and a discrete sequence (C) of states (run, stop, tumble). For cells in 120 μm …

Figure 5—figure supplement 1
Survival probabilities and selected reaction network parameters comparing the 120 μm and 60 μm traps in RL and WL for CR (A–B) and PO (C–D).

Error bars correspond to the standard errors (N = 5).

Figure 6 with 1 supplement
Light-switchable motility of green algae.

(A) Heatmaps of CR swimming speed over the course of a 20-min photoresponse assay, insets show short, representative timeseries in each phase. (B) The same for PO. (C) Population-averaged swimming …

Figure 6—figure supplement 1
Light switchable states data over time for individual cells of CR (A) and PO (B).
Figure 7 with 1 supplement
A novel droplet-fusion assay for querying single-cell motility responses to rapid chemical perturbation.

(A) Example tracks from a bulk motility assay (see Materials and methods), immediately before (’normal motility phenotype’) and after addition of 50 mM KCl to the medium (‘chemically perturbed …

Figure 7—figure supplement 1
More sample trajectories for control (A) and KCl (B) droplet pairs.
Illumination spectra recorded with a spectrometer (OceanInsight OceanHDX, 200 μm-fibre).
Appendix 1—figure 1
Smoothing tracking data using a Savitzky-Golay filter (sgolayfilt).

Example raw and smoothed trajectories for CR (A) and PO (B).

Appendix 1—figure 2
The mean square angular displacement (MSAD) for CR in the 200 µm diameter droplets, fitted to MSAD=2DRτ for τ0.1s.
Appendix 1—figure 3
Probability flux analysis method.
Appendix 1—figure 4
Method For Identifying CR States (A) CR cell velocity for a 20 s segment, with a run/stop threshold of 5 μms-1.

(B) CR cumulative reorientation for the same 20 s segment. (C) CR angular velocity for the same 20 s segment, with a tumble threshold of 32.5rads1. (D) CR states for the same 20 s segment.

Appendix 1—figure 5
Method For Identifying PO States.

(A) PO cell velocity for a 20 s segment, with a run/stop threshold of 5 μms-1. (B) PO median movement over 9-frame windows for the same 20 s segment. (C) PO increases in cell median movement for the …

Appendix 1—figure 6
In simulations the cell is modelled as an asymmetric dumbbell of two spheres a distance l apart.

One sphere represents the cell body and the other models the stroke averaged shape of the flagella. To account for differences in the beating patterns of the two flagella the direction of motion is …

Videos

Video 1
Example videos showing single Chlamydomonas reinhardtii cells trapped by droplet microfluidics inside circular arenas of different diameters (∅ = 200, 120, 60, 40 μm), under red light illumination.
Video 2
Example videos showing single Pyramimonas octopus cells trapped by droplet microfluidics inside circular arenas of different diameters (∅ = 200, 120, 60, 40 μm), under red light illumination.
Video 3
A single Chlamydomonas reinhardtii cell is trapped inside a circular arena of diameter (∅ = 120 μm), under white light illumination.

Videos show different timepoints for the same individual.

Video 4
A single Pyramimonas octopus cell is trapped inside a circular arena of diameter (∅ =120 μm), under white light illumination.

Videos show different timepoints for the same individual.

Video 5
Example video showing a single Pyramimonas octopus cell before and after fusion with KCl solution, under red light illumination.

Tables

Table 1
A comparison of the concentrations of key ions (in mM), before and after paired droplet fusion.
IonPre-fusionPost-fusion (with 10 mM KCl solution)
Na+395197
K+8.339.17
Ca2+8.754.38
Mg2+44.522.3
Cl-459230
SO42-23.611.8
NO3-1.10.6
HPO42-0.1negligible

Additional files

Download links