Anette
Funfak
*a,
Cathy
Fisch
b,
Hatem T.
Abdel Motaal
a,
Julien
Diener
a,
Laurent
Combettes
b,
Charles N.
Baroud
a and
Pascale
Dupuis-Williams
bc
aDepartment of Mechanics, LadHyX, Ecole Polytechnique-CNRS, F-91128 Palaiseau, France. E-mail: anette.funfak@espci.fr
bUMR-S 757 Inserm, Université Paris-Sud, Bat 443, 91405, Orsay cedex, France
cESPCI, 10 rue Vauquelin, 75005 Paris, France
First published on 3rd November 2014
Paramecium cells swim and feed by beating their thousands of cilia in coordinated patterns. The organization of these patterns and its relationship with cell motility has been the subject of a large body of work, particularly as a model for ciliary beating in human organs where similar organization is seen. However the rapid motion of the cells makes quantitative measurements very challenging. Here we provide detailed measurements of the swimming of Paramecium cells from high-speed video at high magnification, as they move in microfluidic channels. An image analysis protocol allows us to decouple the cell movement from the motion of the cilia, thus allowing us to measure the ciliary beat frequency (CBF) and the spatio-temporal organization into metachronal waves along the cell periphery. Two distinct values of the CBF appear at different regions of the cell: most of the cilia beat in the range of 15 to 45 Hz, while the cilia in the peristomal region beat at almost double the frequency. The body and peristomal CBF display a nearly linear relation with the swimming velocity. Moreover the measurements do not display a measurable correlation between the swimming velocity and the metachronal wave velocity on the cell periphery. These measurements are repeated for four RNAi silenced mutants, where proteins specific to the cilia or to their connection to the cell base are depleted. We find that the mutants whose ciliary structure is affected display similar swimming to the control cells albeit with a reduced efficiency, while the mutations that affect the cilia's anchoring to the cell lead to strongly reduced ability to swim. This reduction in motility can be related to a loss of coordination between the ciliary beating in different parts of the cell.
Insight, innovation, integrationGenetic disorders in the cilia or the basal body are known to cause several human diseases, although the contribution of the ciliary beating and metachronal waves to these diseases are often poorly understood. This motivated us to study ciliary beating on Paramecium cells as they swim freely within microfluidic channels, in order to elucidate the relationship that exists between the beating frequency, organization, and swimming efficiency. Advanced image analysis was used to provide quantitative measurements of the ciliary movement, which we related to the swimming of the cells. These measurements were repeated for different RNAi silenced genes associated with the ciliary structure or their anchoring to the cell surface, showing a strong effect of some genes on the swimming efficiency. |
The behavior that is described above, which has long been known for protozoa such as Paramecium, is also observed in a wide variety of other situations. Indeed, it also describes many physiological processes in humans, such as mucous clearance from airways, egg transport in fallopian tubes, or the motion of the spinal fluid in the nervous system.4 As a result, a malfunction in the beating patterns in these organs, whether genetic or acquired, can lead to different diseases (ciliopathies) that result from the accumulation of fluid (hydrocephaly),5 a confusion of the left-right asymmetry during development6 or the loss of mechanical role of the airway epithelium.7 Therefore, the study of cilia beating and coordination, as well as the signalization which conditions them is an essential step to understand how cilia movement is affected by some genetic dysfunction.
The fact that the ciliated cells in the human organs are anchored to a fixed solid boundary while the fluid is pushed by the beating motion simply corresponds to a change of reference frame, in comparison with micro-swimmers, and does not change the underlying mechanics. The strong similarity in behavior between ciliates and human organs, in addition to the conserved structure of cilia and the genes that affect them across species,8 have motivated the study of micro-swimmers as simple models for the observation of the fluid-solid coupling due to ciliary beating.
These issues have motivated a surge of interest in recent years in the study of ciliary motion and fluid flow. In particular, a large body of theoretical and numerical work has looked to identify the relation between the swimming velocity, the emergence and organization of the metachronal waves, and the motion of individual cilia.9–11 Hydrodynamic interactions between individual cilia have also been studied in addition to the optimality of this organization in terms of energy use or feeding ability.12,13 In addition, experimental work on ciliated eukaryotes has recently allowed the measurement of ciliary beating dynamics and subsequent flow fields.14–17
They complete older works, such as those of Machemer;18 here they determined that Paramecium movement relies on the coordinated beating of some five thousand cilia covering the whole cell body, ensuring the cell to swim forward along a spiraling line under normal conditions. This movement may be converted to a backward swimming, due to cilia beating reversal, when the cell encounters an obstacle or changes direction.2 The forward straight swimming is also affected by mechanical or chemical signals which can change the swimming mode by stimulating ciliary ion channels.19 In addition to these global descriptions of Paramecium, two sub-populations of beating cilia can be distinguished in terms of cell swimming: body (somatic) cilia and the cilia of the peristomal region, which both contribute to the cell's swimming, and cilia of the gullet which drive the feeding bacteria into the oral cavity and a group of cilia at the rear of the cell which do not beat and whose role is both cell orientation and sensing. Peristomal cilia, involved both in cell locomotion and feeding, beat stronger backward and in a better coordination than the somatic cilia.2 However, the mechanism which allows neighboring cilia on the cell's surface to beat differently remains mysterious, as well as whether those two categories of cilia respond to the same stimuli.
A few recent experimental studies have nevertheless provided a new level of detail on the beating pattern and the swimming. Machemer20 and Jung21 have measured the cilia beat frequency (CBF) on swimming cells and found that the value is homogeneous along the cell body. Jung et al. report that the different cilia have distinct behavior in response to an increase in viscosity of the surrounding fluid.21
Here we perform detailed quantitative measurements of the cell swimming and its relation with the ciliary beating on five populations of cells: one control strain and four RNAi silenced strains, where the silenced genes are associated with cilia functions or ciliopathies that occur in humans. Two of the silenced genes encode proteins that participate in the cilia structure itself, while the other two affect the anchoring of the cilia to the cell base. We use high-speed video recordings (2000 frames per second) of the cells as they swim in the semi-confined environment of a microfluidic channel which limits their movement to a two-dimensional plane. Image analysis is then used to obtain the cell velocity as well as the movement of the cilia, in the moving reference frame, along the edge of the cell. This then allows us to reconstruct graphically the metachronal waves and to obtain the CBF at every location of the cell and to explore their relation in the different cell populations. The introduced method could further be used for other species of microswimmers and provides an interesting approach to study ciliary beating anomalies in ciliopathies.
Fig. 1 Mean velocity plot (mean ± SE) and swimming tracks of control (WT) and RNAi silenced cells swimming freely in depression slides. |
Fig. 2 Cell tracking and extraction of image sequences obtained from a Paramecium cell swimming in a microchannel. (A) Definition of the border (green) and cell area (blue) for the tracking of a swimming cell (1. indicates the rear of the cell; 2. indicates the front of the cell). (B) Extracted and unfolded image sequence of the ciliary motion (Movie S1 provided in ESI†). (C) Extraction of the space-time diagram, 1. Red line indicates the selected resliced plane for the image sequence. 2. Space-time diagram constructed for a resliced plane at position y = 0. |
The region corresponding to the cilia on each image is then unfolded into a rectangle and a stack is made of the individual images, yielding a video of the ciliary beating in Cartesian coordinates, in the reference frame of the swimming cell (Fig. 2B). A sample movie of the cilia image sequence is provided with the ESI† (see Movie S1). Finally, a space-time diagram (chronograph) of the ciliary movement is produced by selecting a horizontal line on the cilia video and plotting the gray scale along this line as a function of time, as shown in Fig. 2C. In all of the chronographs, the left hand side corresponds to rear of the cell (region labeled 1 on Fig. 2A). Increasing x corresponds to going around the cell in a clockwise fashion, with the front of the cell (region labeled 2 on Fig. 2A) appearing near the middle of the plot. In this particular example, the front of the cell corresponds to a position x = 140 μm. Moving further in x then covers the bottom part of the cell, before closing the loop at the right hand side of the graph. For all of the data presented here, the space-time diagram is obtained one pixel away from the cell border. This ensures that the number of pixels is preserved between the original image and the rectangular mapping, while also increasing the signal to noise ratio.
The cell velocity was determined by numerically differentiating the coordinates obtained from the tracking tool. The typical size of the cell was taken as its perimeter, measured as the length of the spatial vector on the space-time diagram (Fig. 2C).
Four genes have been selected of which two are involved in ciliopathies (DPCD and MKS1) and two (St6 and T16) are structural proteins associated with cilia or their anchoring structures, the basal bodies.
• DPCD (for “deleted in primary ciliary diskynesia”) deletion was shown to induce primary ciliary dyskinesia (PCD) in knockout mice, as a result of impaired organization of axonemal dyneins.24,25 In Paramecium two homologous genes encode DPCD proteins which share 43% identity with their human counterpart.
• Mutations in MKS1 are causing Meckel Gruber syndrome, a severe ciliopathy which affects the development and homeostasis of several organs such as the central nervous system, kidneys, liver, or the skeleton. MKS1 has been shown to be associated with sensory signaling function of cilia and is essential for the cilium stability.26
• St6 is one component of the striated rootlets which anchor the basal body to the intracellular side of the cortex and extend over the cell cortex along the basal body rows. A multigenic family encodes the several striated rootlet proteins that we have analyzed by mass spectrometry. We have ascertained by the synthesis of an antibody that St6 belongs to the striated rootlets and shown that St6 RNAi impairs their assembly.
• T16 is a non conventional Paramecium tubulin. Paramecium possesses a complex multigenic tubulin family with some highly divergent members which form a phyletic branch, distinct from the α, β, γ, δ and ε tubulins. Our analysis of the role of these tubulins has revealed a strong phenotypic effect on cell swimming whilst an anti-peptide antibody directed against the T16 protein labels the cilia tip, suggesting that this tubulin participates in the cilia microtubule shaft.
Eighteen St6, twenty T16, fifteen MKS1 and thirteen DPCD knock-down mutants were analyzed for their swimming behavior and ciliary beating in the microchannels. The phenotype of the control and all investigated mutants as well as the mean cell velocities are shown in Fig. 3.
The images in Fig. 3A–E are representatives of control and each knock-down cell type investigated during swimming in a microchannel. The mean swimming velocity for all cell types, extracted from cell tracking as described in materials and methods, are shown in Fig. 3G. The average length and width of the different cell types are listed in Fig. 3F. While the length of RNAi cells is comparable to control (less than 10% variation for St6, T16 and MKS1), their width displays larger variations, up to a 28% increase in the case of DPCD cells. However, these shape differences do not correlate with swimming velocities, as further highlighted by the examination of cell shape variation for each cell type (see S4 in ESI†).
The T16 phenotype which displays a similar shape as the control cell (Fig. 3A, C, F) has a swimming velocity reduced by 43% (see Movie S5, ESI†) with an average velocity of 413 ± 50 μm s−1 and a maximum velocity of 789 μm s−1.
The DPCD phenotype appeared to be much “fatter” compared to all other cell types (Fig. 3A, E, F). It showed no loss in cilia but a reduced swimming velocity of 471 ± 75 μm s−1 with variability including some cells swimming at a maximum velocity of 1060 μm s−1 (see Movie S6, ESI†). The most affected DPCD cells displayed impaired straight forward swimming patterns but conserved their avoiding reactions when meeting the channel wall.
The swimming velocity was most affected in the MKS1 and St6 knock-down cells with a decrease of more than 70% as compared to the average velocity of control cells. The phenotype of the St6 cells (Fig. 3B) led to a strong reduction of swimming velocity with an average value of 221 ± 43 μm s−1 and a maximum velocity of 644 μm s−1 (see Movie S7, ESI†). Some St6 mutants showed partly inactive cilia resulting in insufficient beating as well as an affinity to stick to the channel walls. This resulted in a staggered swimming behavior.
Similar to the St6, the MKS1 phenotype (Fig. 3D) resulted in abnormal ciliary beating patterns, with cilia being immotile or beating weakly. Furthermore, some MKS1 cells showed partial lack of cilia. This phenotype resulted in the loss of the straight forward swimming and an increased drifting of the cells (see Movie S8, ESI†). We determined a strongly reduced average swimming velocity of 196 ± 38 μm s−1 with a maximum velocity of 410 μm s−1. Interestingly, experiments realized in maze-like microchannels revealed that MKS1 cells have lost their mechano-sensation, evidenced by the inability of those cells to change direction when knocking into the channel walls (see Movies S9 MKS1 and S10 control in ESI†).
In conclusion, all the mutants examined display a significant reduction in swimming velocity which appears to be independent of their shape (see S11, ESI†). The most affected cells are St6 and MKS1 mutants.
The distribution of the CBF on this cell is plotted in Fig. 4Bb. Two frequencies can be distinguished on this cell: F1 = 39 Hz, across large regions of the cell, and F2 = 77 Hz corresponding to the mouth region, where the beating serves for cell feeding. We confirmed this in our experiments using tracer particles (1 μm polystyrene spheres) which regularly flow towards the oral groove (see movie S12, ESI†). This discrepancy in frequencies between the cell body and mouth region was observed for almost all control cells, yielding mean CBF values of 35 ± 4 Hz for the cell cortex and 66 ± 8 Hz for the mouth region. The minority of cells which escaped this rule were tracked while facing with their back to the video, thus masking the mouth waves and preventing the measurement of the mouth CBF.
The mouth CBF was generally observed to be nearly twice the cell body CBF. The fact that the frequency was not exactly twice however implies that it is not simply a harmonic of the CBF. Rather, it likely results from the specific organization of the peristomal cilia (see also ref. 27), which are paired in each cortical unit. In contrast with the mouth region where the CBF is very regular, the rear of the cell (x > 290 μm and x < 20 μm on Fig. 4Ba) is characterized by the presence of caudal cilia that are longer and that rotate without coordination. The CBF measurements in these areas appear as noise in the CBF plots.
The velocity trace of the cells was found to fluctuate even during a smooth swimming period, as shown in Fig. 4Bc where the instantaneous velocity of a control cell is plotted alongside the spatio-temporal diagram. We were not able to relate the velocity fluctuations to any variations in the ciliary activity. Instead, the typical time scale of the fluctuations, which can be estimated from Fig. 4Bc at around 10 ms, was consistent with the viscous diffusion time scale for Paramecium, as discussed by Hamel et al.28
In Fig. 4C the wave's propagation has been indicated by arrows in the space-time diagram. The oblique patterns correspond to the waves that travel from the rear of the cell towards the front, i.e. from the edges of the graph towards the central region (at 150 μm in abscissa). This observation of antiplectic metachronal waves agrees with previous observations. The slope of the oblique lines corresponds to the wave speed, projected in the plane of observation. This value is not constant along the cell body. Instead we observe that the wave speed vary in the range of 1600 to 460 μm s−1. In contrast, the peristomal wave displayed remarkable invariance, originating at the oral groove and travelling towards the tip of the cell, regardless of the cell's swimming direction or speed. The signature of this wave on the edge of the cell is observed in the range x = 100 to 140 μm on Fig. 4C. The features discussed here were observed on almost all control cells. However, although the wave velocities varied from cell to cell, it was not possible to find a correlation between the wave velocities and the cell swimming velocity.
The space-time diagrams and frequency plots of T16 and DPCD knock-down cells (Fig. 5A and B) are similar to the ones obtained with control cells: they display distinctive CBF values for body and mouth regions with respective values either equivalent to the control cells for DPCD (34 ± 5 Hz; 66 ± 9 Hz) or slightly reduced for T16 (32 ± 6 Hz; 60 ± 8 Hz) silenced cells.
On the opposite, the space-time diagrams of most St6 and MKS1 cells are characterized by fuzzy CBF signals with strong fluctuation, restricted to the cell body region for St6 knock-down cells but affecting both body and mouth regions for MKS1 silenced cells (Fig. 5C and D). In contrast with the other cell types, the St6 silenced cells displayed some phenotypical variance that we interpreted as a variability of RNAi efficiency. Indeed, some cells retained two well defined frequencies over most of the cell body in addition to visible wave propagation while others displayed large regions with immotile cilia or cilia beating very slowly (see S13, ESI†). These regions of uncoordinated ciliary beating generated the noisy appearance of the space-time diagrams over the cell body region while the mouth region was not concerned by these uncoordinated beating patterns.
Moreover, the most affected St6 cells displayed a distorted morphology with the presence of a bulge near the mouth. This impaired the visualization of waves which appeared to travel transverse to the cell's major axis. As a result, the mouth waves appeared as standing waves on the space-time diagrams (see image S14, ESI†). However, although those distortions had a strong effect on the mouth wave propagation, they did not affect the CBF which displays a homogeneous reduction of 20% compared to the control cells, both for body and mouth regions (29 ± 5 Hz and 61 ± 7 Hz, respectively).
MKS1 silencing caused the strongest reduction in the CBF of the body region of approximately 30% with a mean value of 22 ± 6 Hz. A striking characteristic of those cells is the absence of ciliary motion and wave propagation in the mouth region although we did not note cilia loss in this region. However, like in St6 RNAi cells, the most affected MKS1 silenced cells displayed partial loss of cilia in the body region. Cells lacking too many cilia were excluded from the analysis, the corresponding diagrams being too noisy to extract a mean CBF value. Overall, MKS1 RNAi cells were characterized by a loss of synchrony of the ciliary beating illustrated by the absence of regular waves on the cell body and in even more remarkably in the peristomal region, as well as the very noisy space-time diagrams (demonstrated in the space-time diagram S15 in the ESI†).
Interestingly, although all the mutants are characterized by a low swimming velocity, their cilia beating dynamics is not affected in the same way. Surprisingly, knock-down of T16 and DPCD, which are core proteins of the axoneme and involved in structural stability does not significantly impact beating frequency. This has to be compared to another Paramecium mutant Bug 2229 affected in cilia rigidity, shown to be strongly impacted in swimming velocity and cilia beating trajectory but beating frequency was not reduced. Thus, in these examples, it seems that beating frequency is more affected by mutations in proteins involved in cilia anchoring than in axonemal stability.
This also applies to the knock-down mutant DPCD where most cells swim with a CBF over 30 Hz (range of 22 to 42 Hz) with corresponding velocities between 500 and 1000 μm s−1. T16 knock-down cells display a wider range of swimming patterns with half of the cells being slower (with a cell velocity below 500 μm s−1) and corresponding CBF below 30 Hz. However, these two knock-down mutants display like the control cells a strong correlation between swimming velocity and body CBF, as attested by the curves in Fig. 6F. Accordingly, the curves for DPCD and T16 (Fig. 6F respectively violet and green lines), although shifted towards lower velocities, are almost parallel to the control cell's curve, showing that the same increase in CBF will produce the same progression in swimming velocity. This indicates a ciliary beating which is less efficient in terms of swimming in the DPCD and T16 than in control cells. However, the acceleration of the ciliary beating brings the same effect on swimming acceleration for the three cell types.
Contrasting with DPCD and T16, the MKS1 and St6 knock-down mutants are uniformly characterized by very low swimming velocities (below 500 μm s−1). However, whereas MKS1 knock-down cells all display very low CBF values (below 30 Hz), St6 knock-down cells have a wider distribution of their CBF (values ranging mostly from 20 to 35 Hz), similar to the T16 RNAi mutant. The correlation coefficients calculated from the best linear regression curves obtained from these two strains (Fig. 6F respectively blue and red lines) have smaller values, indicating that the correlation between the cell velocity and CBF is not so strong. Furthermore, the change in the slope for these two RNAi mutants indicate that in addition to a bad cilia beating efficiency, their ability to accelerate in response to the increase of CBF is impaired. In the case of the St6 phenotype, an increase of CBF about 10 Hz resulted in a halved velocity change in comparison to the control (Fig. 6F St6 plot, 180 μm s−1, and control plot, 360 μm s−1). However, the correlation curve obtained for MKS1 knock-down cells must be treated with some caution since it is underlying a statistical uncertainty because of the partial deciliation induced by MKS1-RNAi. Thus, while most cells swim slowly due to slow CBF, the slow swimming velocity of the rare cells which have a CBF over 30 Hz could result from a lack of cilia rather than a low cilia beating efficiency. However, since it is not possible to evaluate the ciliary loss of each MKS1 RNAi cell, and that we know that ciliary instability is a characteristic of the MKS1 silencing, we decided to consider the whole MKS1 RNAi cell population and treat them all together.
Strikingly, the St6 knock-down cells are the ones which display the least efficient cilia beating among the five strains examined in this experiment, with a ratio of cell velocity against CBF which is the lowest over the whole CBF scale. In conclusion, the St6 knock-down mutant is characterized both by an inefficient ciliary beating and an inability to increase its speed concurrent with an increased beating frequency.
Note that the linear relationship between the swimming velocity and the CBF is an indication that the ciliary beating pattern is independent of the frequency for a given cell population. This is expected in the low Reynolds number regime, if the fluid-structure coupling with the cilia is weak. The slope of the velocity/CBF curves defines a length scale which is in the range of 10 to 35 μm, suggesting the presence of a characteristic length scale for each cell population. We verified this by measuring the spatial correlation coefficient C on each of the space-time diagrams (see e.g. plot S16, ESI†). These curves display a decrease to zero over a length scale comparable with the slope of the velocity/CBF curves. Indeed, a typical correlation length LC is extracted from each measurement by looking for the distance at which C crosses the value of 0.1 (in order to reduce the sensitivity to the noise near C = 0). The values of LC are plotted on Fig. 7, with respect to the corresponding slopes of the velocity/CBF curves (see Fig. 6F). They show a nearly linear relationship and confirm the presence of such a typical length within the spatial organization of the ciliary beating. While we expect that this length scale LC is related to the observed wavelength of the metachronal waves, the current measurements do not allow us to clearly extract this information and more analysis is required to make a direct link between the lengths scales that interact.
Fig. 7 Slope–correlation length (LC) plot where each dot represent a cell population. LC was obtained using a spatial correlation function on each cell population (see plot S16, ESI†). |
Finally, we emphasize that it is possible for Paramecium to beat its cilia at a finite frequency while remaining nearly stationary. This confirms visual observations that the different ciliary strokes are possible and that the cell can perform different ciliary movements than the standard power-stroke/recovery stroke. This implies in turn that the CBF, while strongly affecting the cell velocity, is not the only factor determining the swimming pattern and speed. As regards the mouth CBF, a similar plot displayed correlation with the swimming velocity except for St6 and MKS1, for which no correlation could be found (see plot S17, ESI†).
Conversely, precise measurements were possible for the ciliary beat frequency (CBF) at all locations on the cell periphery. The cells were found to display two distinct beat frequencies: one frequency along the cell body, showing a strong correlation with the swimming velocity, and another higher frequency emanating from the peristomal region which showed a correlation with the swimming velocity for certain cell types. Since this value is in conflict with previous measurements,21 we attempted to check if this higher frequency was an artifact, due to the doubling of the cilia in this region, but this does not appear to be the case from our videos. Instead, the effect may be attributed to the higher effort that is exerted by the doubled cilia which would then beat at a frequency nearly double the frequency elsewhere. This observation of a distinctive behavior of the peristomal cilia, intended to create a flux towards the gullet, supports the pioneering observations of Kuznicki et al.,30 who noticed that those cilia beat independently of the other group of cilia and in a better coordination. This was recently confirmed by21 who demonstrated that the beating of those cilia, in contrast with the body cilia is insensitive to variations in viscosity, brushing up to date the questions of the control of the movement of neighboring groups of cilia in ciliated protozoa.31
Finally, the measurements were repeated on control cells and four different RNAi induced knock-down mutants. Those knock-down cells have been selected because their phenotype was characterized by a significant reduction in swimming velocity as revealed by cell tracking in microwells (Fig. 1). The reduced cell velocity was confirmed in the microfluidic channels for all knock down mutants although with distinctive characteristics. The mean velocity of control cells was lower in the channels (780 μm s−1) than in the wells (1100 μm s−1), suggesting that this could result from the confinement in the microchannels, as suggested by.3 However, although the slowdown in the microchannels of the St6 and MKS1 mutants, which are fatter than the control cells, agrees with this hypothesis, the increased velocity in microchannels of the DPCD mutant which has the largest width refutes it. We did not observe a clear correlation between the ratio cell width against channel diameter and cell velocity in our experiments. Nevertheless, it is likely that the spiraling parameter, which influences cell's velocity3 is affected differently in the various phenotypes, making it difficult to draw any conclusions about the drag effect of channel boundaries in our experiments. However, it highlights the advantages of microfluidic channels to compare mutants whose particular unconfined swimming behavior impedes precise comparative measurements of their swimming velocities.
This method also proved to be very valuable in the analysis of the impact of knock-down mutations on cilia beating and subsequent cell swimming behavior. Although all the knock-down mutants were characterized by a decreased cell velocity compared to control cells, two categories could be distinguished with T16 and DPCD cells in one part (with a reduction in mean velocity less than 50% compared to control cells), and St6 and MKS1 whose mean velocity is reduced by more than 70% compared to control cells. We also show that in our conditions, cell velocity variation is not attributable to differences in cell morphology.
Our analysis demonstrates that although the ciliary beating efficiency of all knock-down mutants is lower in terms of swimming velocity compared to control cells, as attested by the ratio cell velocity against CBF, their ability of increasing their swimming velocity according to the acceleration of ciliary beating is unchanged for DPCD and T16 cells and highly affected in MKS1 and furthermore for St6. This suggests that in addition of impairing individual cilia beating (such as beating shape, kinetics or strength), St6 and MKS1 knock-down affect the efficiency of collective cilia movement. This is further supported by our experiments which show distinctive space-time diagrams for the two knock-down mutants categories, from which consistent values of CBF were observed over the whole cell periphery in DPCD and T16 knock-down cells, whereas noisy graphs of MKS1 and St6 silenced cells attest of a lack of cilia beating coordination along the cell's surface. Interestingly, this analysis is pertinent when considering the nature of the knock-down mutants: DPCD and T16 are both structural proteins of the ciliary microtubule scaffold and therefore expected to influence a cilium's individual movement. On the other hand, MKS1 and St6 influence cilia anchoring on the cell surface: St6 being involved in kinetodesmal fibers which anchor the basal body in the Paramecium cortex, and MKS1, localized at the transition zone between basal body and cilium's core, is essential for cilium stability.32 Our work thus provides experimental evidence for the respective contribution of cilium structure and anchoring on cilia movement and synchronization, as previously analyzed in theoretical models.12
In conclusion, our results contribute to understanding the movements of micro-swimmers but can also contribute to the characterization of the defects involved in ciliary diseases, which are still often restricted to the characterization of the structural defects of cilia.33 Our analysis suggests that the coordination of the ciliary beating should be considered as a major factor in the etiology of ciliary dyskynesia and highlights the benefit of the use of ciliated protozoa like Paramecium for which sophisticated genetic tools enables to phenocopy mutation involved in ciliary diseases.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c4ib00181h |
This journal is © The Royal Society of Chemistry 2015 |