The biomechanical role of extra-axonemal structures in shaping the flagellar beat of Euglena gracilis
Figures

Inner structure of Euglena gracilis' flagellum.
(a) A specimen of freely swimming Euglena gracilis, and (b) a sketch of the cross-section of its flagellum, as seen from the distal end. The flagellar inner structure is composed by the paraflagellar rod (PFR, textured), and the axoneme (Ax). The PFR is connected via bonding links to the axonemal doublets 1, 2, and 3. The inner structure of the flagellum is enclosed by the flagellar membrane (dotted contour). By inhibiting MTs’ sliding, the PFR selects the spontaneous bending plane of the Ax (dashed line). As a key geometric feature, the solid line that joins the Ax center and the PFR center crosses at an angle the spontaneous bending plane. Doublets are numbered following the convention adopted in the electron microscopy studies Melkonian et al., 1982 and Bouck et al., 1990, to facilitate comparison. The opposite convention, in which microtubules are numbered in increasing order in the anti-clockwise direction when seen from the distal end of the Ax, is far more common in structural studies of cilia.

Flagellar beat kinematics of freely swimming Euglena gracilis.
(a) flagellar configurations in evenly spaced instants (phases) within a periodic beat. (b) The same configurations overlapped and color coded according to their phases. (c) Computed torsion as a function of the flagellar arc length . The plot is presented in terms of the normalized quantities and , where is the total length of the flagellum. Panels (a-b) are adapted from Figure 5.E of Rossi et al., 2017.
-
Figure 2—source code 1
Experimental flagellar waveforms and torsion calculator.
- https://cdn.elifesciences.org/articles/58610/elife-58610-fig2-code1-v2.zip

Flagellar beat of capillary trapped specimens.
(a) A specimen of Euglena gracilis trapped at the tip of a capillary (bottom). The typical outline of its beating flagellum is that of a looping curve, which is consistent with the outline of a simple curve with two concentrated torsional peaks of alternate sign along the length of the curve, that is, a torsion dipole (top). (b) Close-up images of the same specimen of capillary-trapped (CT) E. gracilis as seen from different viewpoints, upon successive ∼90° turns of the capillary tube. The body orientation with respect to the objective is estimated from the anatomy of the cell, and in particular from the position of the eyespot (a visible light-sensing organelle present on the cell surface). Microscopy images are decorated with the tracked outlines of the flagellum in different phases (same color coding as in Figure 2). The outlines (2d projections) of the 3d reconstructed flagellar beat of freely swimming (FS) specimens are shown for comparison.

Details of the mechanical model.
(a) Geometry of the Ax. MTs lie on a tubular surface parametrized by generalized polar coordinates and , where is the arc length of the axonemal centerline . The unit vectors and lie on the orthogonal cross sections of the Ax (light blue circles). The material sections of the Ax are given by the curves (red), which connect points of neighbouring axonemal MTs corresponding to the same arc length . Bend deformations of the axoneme are generated by the shear (collective sliding) of MTs. The shear is quantified by the angle between the orthogonal sections and the material sections of the Ax. (b) Geometry of the euglenid flagellum, detail of the Ax-PFR attachment. The unit vectors and generate the plane of the PFR’s cross sections. The vector is parallel to the outer unit normal to the axonemal surface , while is parallel to the tangent vector to the material section .

Flagellar non-planarity arising from structural incompatibility.
The Ax-PFR mechanical interplay is explained in a three-steps argument (left to right). Consider first the two separated structures in their spontaneous configurations (left). The Ax is bent into a planar arc while the PFR is straight. Then, the PFR is forced to match to the Ax, while the latter is kept in its spontaneous configuration (middle). The attachment constraint induces shear strains in the PFR, such that the composite system cannot be in mechanical equilibrium without external forcing. When the composite system is released (right), it reaches equilibrium by the relaxation of the PFR shear, which induces additional distortion of the Ax. At equilibrium, an optimal energy compromise is reached, which is characterized by an emergent non-planarity.

Geometry and mechanics of non-planar flagellar shapes.
(a) The bending vector traces a curve on the plane of the bending parameters and . The norm of the bending vector determines the curvature of the flagellum. The rate of change of the angle determines the torsion . (b–f) Bending vectors’ traces of flagellar equilibrium configurations under the same (steady) dynein actuation, but different values of the material parameter . Equilibria are minimizer of the energy . For small values of , the Ax component of the energy dominates. In this case, is close to the target bending vector where . For large values of the PFR component of the energy dominates, and equilibria are dragged closer to the line orthogonal to the vector (dashed green). The bending vector undergoes rotations which result in torsional peaks of alternating sign.
-
Figure 6—source code 1
Equilibrium equations solver.
- https://cdn.elifesciences.org/articles/58610/elife-58610-fig6-code1-v2.zip

Kinematics of the beating euglenid flagellum: comparison between theoretical model and experiments.
(a-b) Dyneins’ shear forces. (c) Resulting bending strains and torsion for an Ax actuated by the force pattern (a–b), beating in a viscous fluid, and free of extra-axonemal structures. The beat is planar (Chlamydomonas-like). (d) Resulting bending strains and torsion for an euglenid flagellum (composite structure Ax+PFR) actuated by (a–b) and beating in a viscous fluid. The Ax-PFR interaction generates torsional peaks with alternate sign traveling from the proximal to the distal end of the flagellum. (e) Resulting shapes for the euglenid flagellum at different instants within a beat, and comparison with experimental observations.
-
Figure 7—source code 1
Flagellar dynamics solver and visualization tool.
- https://cdn.elifesciences.org/articles/58610/elife-58610-fig7-code1-v2.zip

Swimming kinematics.
(a) Side view and (b) top view of swimming cell simulation resulting from the flagellar beat generated by our model. The dimension of the cell body is not to scale with displacements for visualization purposes.

Sketch of two MTs’ centerlines during deformation.
The sliding is defined as the difference between the arc lengths and . The latter is the arc length corresponding to the projection of on the curve . We have positive sliding when dyneins push the -th MT toward the distal end of the flagellum and the -th MT toward the proximal end.

Bending vector traces and torsion of equilibrium configurations of the euglenid flagellum under four different target bending strains, and for different values of the parameter .
The target bending strains are given by , with and , where (a) , (b) , (c) , and (d) .

Swimming kinematics: numerical results.
(a) Cell body and body frame unit vectors. (b) Body velocity and (c) rotational velocity projections on the body frame vectors during a flagellar beat.
-
Appendix 5—figure 1—source code 1
Swimming dynamics solver.
- https://cdn.elifesciences.org/articles/58610/elife-58610-app5-fig1-code1-v2.zip
Videos
Four views of a capillary-trapped specimen of E. gracilis recorded during periodic flagellar beating.
Comparison between observations of a beating euglenid flagellum and the numerical simulations of our mechanical model.
Simulations of swimming cell kinematics resulting from the flagellar beat generated by our model.
Additional files
-
Transparent reporting form
- https://cdn.elifesciences.org/articles/58610/elife-58610-transrepform-v2.docx
-
Appendix 5—figure 1—source code 1
Swimming dynamics solver.
- https://cdn.elifesciences.org/articles/58610/elife-58610-app5-fig1-code1-v2.zip