henhaobio
Sign in

bioRxivBiophysicsPreprint

Heterogeneity in cilia patterning enables multiple flow functions in Paramecium

Daphne M. LaaniD, Aikaterini Maria Kourkoulou, Guillermina R. Ramirez-San JuaniD

PDFDOICode

Abstract

Free-living unicellular organisms known as ciliates rely on fluid flows to perform essential functions. These flows emerge from the coordinated activity of thousands of cilia organized into arrays with highly diverse architectures. Within a species, the structure of the ciliary array is highly stereotyped and persists in time, raising the question of how a single cell achieves multiple distinct flow functions from one fixed array architecture. Here, we investigate how the ciliary array in the ciliate Paramecium tetraurelia enables this organism to feed and swim simultaneously. Using expansion microscopy and high-speed imaging, we measure ciliary organization and kinematics, from individual beat dynamics to collective metachronal wave patterns, together with the flows these cilia generate. We find that the cell's surface is partitioned into domains with distinct spatiotemporal cilia organization. In addition to the oral apparatus, there are two structurally different domains. One is a densely ciliated, high-frequency beating domain located anterior to the oral apparatus, and the other covers the remainder of the cell's surface, where cilia are more sparsely distributed and slower-beating. Selective removal of each results in impaired feeding or swimming, respectively, demonstrating that each domain is functionally specialized. Our work establishes heterogeneity in spatiotemporal ciliary patterning as a mechanism that allows a single cell to generate multiple flow functions. By measuring cilia organization and the flows they generate, we uncover general physical principles by which cilia arrays produce flow patterns that underlie diverse biological functions.

In-depth summary

Claude Opus 5.5 · xhigh effort · 8 Nov 2026 · from the full text of version 2 (16 Oct 2026) · CC BY 4.0

Paramecium tetraurelia partitions its ciliary array into a dense, fast-beating domain anterior to the oral apparatus that drives feeding flows and a sparser, slower domain that drives swimming.

A ciliate such as Paramecium feeds while it swims, using one fixed, stereotyped ciliary array and without the neural control that lets marine larvae retune ciliary kinematics to switch functions. Earlier Paramecium studies reported region-to-region variation in basal body arrangement and beat frequency, but no whole-cell quantitative map had linked cilia density, orientation and metachronal kinematics to the flows they produce and the functions those flows serve.

Key results

  1. ShownFeeding depends on domain D1 (dense cilia anterior to the oral apparatus, OA). Cells lacking D1 cilia ingest essentially no beads. Cells keeping only D1 and OA cilia feed like treated cells with minor, uniform cilia loss. OA cilia alone do not support uptake.

    Intracellular/background fluorescence ratio after 10-15 min with 1 µm beads: control 1.32 ± 0.10, Uniform 1.21 ± 0.11, D1+OA 1.18 ± 0.11 (ns vs Uniform), D2+OA 0.99 ± 0.03 (significantly below Uniform, like fully deciliated or OA-only cells). Untreated, unvortexed cells: 1.49 ± 0.14. Mann-Whitney U with Holm correction; 6 to 22 cells per condition. · DISTINCT CELLULAR REGIONS CARRY OUT SPECIFIC…; THE CILIARY ARRAY GENERATES STEREOTYPED FLOWS…; e. Feeding analysis; f. Statistical analysis

  2. ShownSwimming depends on domain D2 (the sparser cilia over the rest of the cell). Cells keeping only D1 and OA cilia barely move. Cells lacking D1 show no significant loss of speed (vs uniformly ciliated treated cells) and no change in gait time fractions or transitions (vs untreated cells).

    Speed: control 0.41 ± 0.04 mm/s, Uniform 0.43 ± 0.25, D2+OA 0.30 ± 0.06 (ns vs Uniform), D1+OA 0.09 ± 0.07 (significantly slower); OA-only cells immotile. D2+OA cells (n=11) spent 91 ± 6% of time in forward swimming vs 88 ± 8% for untreated cells (n=59). Transition probabilities from forward swimming to reverse/reorient: 0.7 ± 0.1/0.3 ± 0.1 vs 0.7 ± 0.2/0.3 ± 0.2, none significant. · DISTINCT CELLULAR REGIONS CARRY OUT SPECIFIC…; THE CILIARY ARRAY GENERATES STEREOTYPED FLOWS…; f. Statistical analysis; Tables

  3. ShownThe cortex divides into a dense domain (D1) anterior and right of the oral apparatus, where basal bodies sit in ciliated pairs, and a sparser domain (D2) over the rest of the ventral surface and all of the dorsal surface. Cilia density is about 1.4-fold higher in D1.

    Expansion microscopy: intra-row spacing 2.28 ± 0.22 µm between basal body pairs in D1 vs 1.89 ± 0.17 µm between basal bodies in D2 (n=10 cells); inter-row 2.28 ± 0.16 vs 1.93 ± 0.25 µm (n=27/28 cells). Because both basal bodies of each D1 pair are ciliated and fewer D2 basal bodies are, cilia density is 1.4-fold higher in D1. Areas: D1 16.1 ± 2.8%, D2 80.4 ± 3.0%, OA 3.5 ± 1.5% (n=12). · CILIA ARRAYS WITH DISTINCT DENSITY AND…; a. Measurement of cell shape and cilia array…; Figure legends

  4. ShownD1 cilia beat about twice as fast as D2 cilia but have the same metachronal wavelength, so waves travel roughly twice as fast in D1. In both domains waves head anteriorly at ∼30° to the anterior-posterior (A-P) axis, while stroke direction follows the local striated fibers.

    High-speed DIC (1000 fps) of immobilized cells, n=15: beat frequency 35.4 ± 3.0 Hz in D1 vs 15.7 ± 2.5 Hz in D2; wavelength 9.1 ± 1.0 µm vs 9.1 ± 0.6 µm; wave speed 320.6 ± 39.9 vs 143.9 ± 26.7 µm/s. Waves are dexioplectic in D1 and mixed dexioplectic/antiplectic in D2. Paroral membrane cilia: 30.2 ± 2.1 Hz, stroke at 80.6 ± 7.7° to the A-P axis. · INDIVIDUAL CILIA DYNAMICS AND COLLECTIVE…; b. Analysis of ciliary beat direction; c. Quantification of metachronal wave properties; Figure legends

  5. SuggestedParticle image velocimetry (PIV) shows two counter-rotating vortices and stronger ventral flow aimed at the oral apparatus. Removing D1 erases the mouth-directed anterior flow. Removing D2 keeps the ventral flow but leaves negligible dorsal flow.

    Untreated immobilized cells (n=4): ∼650 µm/s ventral vs ∼250 µm/s dorsal near the cell, decaying to 100 µm/s at 200 µm. D2+OA cells (n=5) kept posteriorly shifted vortices, and their normalized radial profiles matched uniformly ciliated treated cells (n=3). D1+OA cells (n=11) matched those cells (n=3) ventrally but had negligible dorsal velocity at every distance. · THE CILIARY ARRAY GENERATES STEREOTYPED FLOWS…; DISTINCT CELLULAR REGIONS CARRY OUT SPECIFIC…; g. Analysis of flow fields; Figure legends

Methods, data and reproducibility

Expansion microscopy of basal bodies and striated fibers mapped cilia density and orientation across whole cells; 1000 fps DIC imaging of immobilized cells measured stroke direction, beat frequency and metachronal waves. Bead-tracer PIV, bead-ingestion assays and tracking of free swimmers quantified flows, feeding and swimming gaits, before and after ethanol- or sucrose-induced partial deciliation.

  • Fixed, expanded cells (U-ExM) imaged by confocal microscopy for basal body spacing, striated fiber orientation and domain areas · Paramecium tetraurelia (wild-type strain d4-2) · n = 10 intra-row; 27 D1/28 D2 inter-row; 12 areas; 10 ventral/9 dorsal fibers
  • Confocal immunofluorescence images (centrin) for cell length and width · Paramecium tetraurelia (wild-type strain d4-2) · n = 79 cells
  • High-speed DIC recordings (1000 fps, 13 s) of coverslip-immobilized cells for stroke direction, beat frequency and metachronal waves · Paramecium tetraurelia (wild-type strain d4-2) · n = 15 cells (beat frequency and waves)
  • Free-swimming trajectories (5 min at 30 fps) in custom acrylic chambers · Paramecium tetraurelia (wild-type strain d4-2) · n = 59 untreated; 11 D2+OA (gaits); 5 to 14 per deciliation condition (speed)
  • PIV flow fields from 0.5 µm tracer beads around Cell-Tak-immobilized (baseline) and deciliated or uniformly ciliated treated cells · Paramecium tetraurelia (wild-type strain d4-2) · n = 4 untreated; sucrose 3 uniform, 5 D2+OA; ethanol 3 uniform, 11 D1+OA
  • Bead-uptake feeding assays with 1 µm fluorescent beads · Paramecium tetraurelia (wild-type strain d4-2) · n = 9 untreated cells; 6 to 22 cells per deciliation condition
  • Screening of individual ethanol- or sucrose-treated cells for deciliation pattern · Paramecium tetraurelia (wild-type strain d4-2) · n = 114 ethanol-treated, 74 sucrose-treated cells
Methods & tools
Ultrastructure expansion microscopy (U-ExM), High-speed DIC microscopy (Phantom T2540 camera), Nikon Ti2 AX confocal with NSPARC detector, Particle image velocimetry (FluidImage), Alignment by Fourier Transform (AFT), micro-sam (Segment Anything for microscopy), TrackPy, scikit-image, SciPy, StackReg / pystackreg, ImageJ, Mann-Whitney U test with Holm correction
Reproducibility gaps
The data behind the plots are promised only once the manuscript is published; no repository or accession is given. Software versions and environment files are not stated.
Limitations (6)
  • Cilia kinematics and baseline flows were measured on cells stuck to a coverslip. The authors note that the boundary could change beating compared with free swimmers, and that resolving cilia on fast, rotating swimming cells is not currently feasible. · stated by the authors
  • Why D1 cilia beat about twice as fast as D2 cilia remains speculative (cilia density and coupling, or axonemal biochemistry), and any role of calcium signaling is untested. · stated by the authors
  • How D1 and D2 are specified and maintained within one continuous cortex is unknown. Genetic tests (for example RNAi of cortical scaffolds, or patterning mutants such as crochu) are left to future work. · stated by the authors
  • The dexioplectic pattern seen in D1 conflicts with an older antiplectic description of this region in P. multimicronucleatum. The authors attribute the difference to earlier fixed-sample methods or possibly to species differences. · stated by the authors
  • The domain-removal conditions come from two different chemicals with low, imperfect yields (7 of 74 sucrose-treated cells were D2+OA; sucrose often also stripped anterior D2, and ethanol left polar cilia). Cells were screened by hand under DIC. · noted in this summary
  • The flow-field conclusions rest on few cells (4 untreated; 3 to 11 per deciliation condition), and every experiment uses one wild-type strain (d4-2). · noted in this summary

Read this if you study ciliary coordination, metachronal waves or low-Reynolds-number flows and want whole-cell measurements tying cilia architecture and kinematics to swimming and feeding in one cell.

Citations
1

Reproducibility ledger

Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.

Sign in to file a reproducibility report.