Abstract
Cyanobacteria are key contributors to biogeochemical cycles through photosynthesis and carbon fixation. In filamentous, multicellular cyanobacteria these functions can be influenced through gliding motility, which enables filaments to localise in response to light and also form aggregates. Here, we use the aggregate forming species Fluctiforma draycotensis to study gliding motility dynamics in detail. We find that filaments move in curved and straight trajectories interspersed with re-orientation or reversal of direction. Most reversals take few seconds but some take substantially longer, resulting in a long-tailed distribution of stoppage times. Mean filament speeds range around a micron per second with a relatively uniform distribution against filament length, implying that all or fixed proportion of cells in a filament contribute to movement. We implement a biophysical model that can recapitulate these findings. Model simulations show that for filaments to reverse quickly, cells in a filament must achieve high coordination of the direction of the forces that they generate. To seek experimental support of this prediction, we track individual cells in a filament. This reveals that cells’ translational movement is fully coupled with their rotation along the long-axis of the filament, and that cellular movement remains coordinated throughout a reversal. For some filaments, especially longer ones, however, we also find that cellular coordination can be lost, and filaments can form buckles that can twist around themselves, resulting in plectonemes. The experimental findings and the biophysical model presented here will inform future studies of individual and collective filament movement.
Significance Statement
Cyanobacteria contribute to global oxygen production and carbon capture. Some cyanobacteria exist as multicellular filaments and display gliding motility that allows them to respond to light and to form aggregates, which influences their biological functions. Here, we study the dynamics of gliding motility. We find that filaments’ movement is interspersed with re-orientation or reversal of direction and that mean filament speed is mostly independent of filament length. We implement a biophysical model that predicts these features to relate to cells in a filament having high coordination of the direction of the forces that they generate. We find experimental support for this predicted cellular coordination, but also discover instances of longer filaments loosing coordination, resulting in buckling and entangling with other filaments.
eLife assessment
Using microscopy experiments and theoretical modelling, the authors present convincing evidence of cellular coordination in the gliding filamentous cyanobacterium Fluctiforma draycotensis. The results are important for the understanding of cyanobacterial motility and the underlying molecular and mechanical pathways of cellular coordination.
Significance of findings
important: Findings that have theoretical or practical implications beyond a single subfield
- landmark
- fundamental
- important
- valuable
- useful
Strength of evidence
convincing: Appropriate and validated methodology in line with current state-of-the-art
- exceptional
- compelling
- convincing
- solid
- incomplete
- inadequate
During the peer-review process the editor and reviewers write an eLife assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife assessments
Introduction
Cyanobacteria are key contributors to global primary production and oxygen generation [1, 2]. They display high diversity and are adapted to a wide range of habitats from the soil crust to freshwater, and the ocean [2, 3]. Within this diversity, some cyanobacterial species display multicellularity, existing as filaments composed of individual cells. Many of these filamenteous cyanobacteria display a specific type of motility, known as gliding motility [3, 4, 5, 6]. Gliding motility enables filaments to position themselves in response to light [7, 8, 9, 10] and also to form aggregates [11, 12, 13], thereby enhancing photosynthesis and additional functions such as nitrogen fixation and iron acquisition [14].
Gliding motility is a specific form of bacterial motility that does not involve flagella [15, 16, 17, 18, 6, 19]. In single celled bacteria, studies in the species Flavobacterium johnsoniae and Myxococcus xanthus have shown that gliding motility involves membrane-bound protein complexes moving in helical patterns across and around the cell body [20, 21, 22, 23]. These membrane proteins are suggested to reach so-called focal adhesion points, at cell-surface contact points, where they convert their motion into axial forces that translate the cell and rotate it around its long axis [24, 18]. Mutational studies in these species, as well as filamenteous cyanobacteria, have also identified proteins involved in the biosynthesis and secretion of a sugar polymer, so-called slime, to be essential for gliding motility [15, 16, 17, 18, 25, 26]. In filamentous cyanobacteria, helically organised surface fibril proteins [27, 28, 29, 30] and junctional protein complexes [31] were identified and suggested to link with motility and slime secretion [31, 27]. Type IV pili machinery is also implicated in filamenteous cyanobacteria gliding, both in Nostoc punctiforme, where only differentiated filaments called hormogonia display gliding motility[4], and in other species [32, 33, 34, 8].
Despite these ongoing efforts on identifying the molecular mechanisms of force generation, the dynamics of gliding motility in filamentous cyanobacteria remains poorly characterised. In particular, it is not clear how filamentous cyanobacteria achieve coordination of cellular propulsive forces during gliding [6]. A better understanding of motility dynamics can provide insights into cellular coordination and inform molecular studies of propulsion. Study of the motility dynamics at the filament level is also relevant to population-level observations, such as aggregate formation. Current theoretical explanations for these observations make use of specific assumptions about filament behaviour [35, 36], which can be better justified, or ruled out, upon better understanding of movement dynamics.
Here, we study gliding motility dynamics in a filamentous cyanobacterium Fluctiforma draycotensis, capable of gliding motility and extensive aggregate formation [13]. We focus specifically on characterisation of movement dynamics, rather than elucidating their underlying molecular mechanisms. To this end, we used extensive time lapse microscopy under phase, fluorescent, and total internal reflection (TIRF) modalities. We found that gliding motility involves movement on curved or straight trajectories interspersed with rapid re-orientation or reversal events and with mean speeds independent of filament length. We were able to recapitulate these findings with a physical model of filament movement, which predicts a high importance for cellular coordination to achieve rapid reversals. We were able to confirm this prediction by tracking individual cells during a reversal. The single-cell analysis also revealed a direct coupling of translational and rotational movement, which can generate torsional forces during reversals. In support of such forces, we found that longer filaments tend to readily buckle and twist upon themselves to form plectonemes. The formation of plectonemes is associated with loss of cell-to-cell coordination during reversals, leading to filament ends moving independently. The presented findings quantify individual filament movement and provide a physical model of gliding motility that is consistent with experimental observations. Together, they will inform future studies on the molecular and physical mechanisms of force generation and collective behaviour of filaments.
Results
Filament movement is interspersed with stoppage events, whose duration has a long-tailed distribution
Using time-lapse microscopy, we observed filament movement on glass slides and under agar pads (see Methods). We found that filaments move on straight or curved trajectories that are interspersed with re-orientation or reversal of direction (see Supplementary Movies S 1 and S 2). Under agar, it was possible to distinguish a region of different phase contrast, which aligns with the trajectories of the filaments, as forming a ‘track’ (Fig. 1A and Movie S 3) (we note that it is possible that the tracks are associated with the secreted slime, see Discussion section). Lone filaments’ trajectories are mostly confined to these tracks, involving a stoppage and reversal at each end of the track (Fig. 1B), while filaments entering a circular track did not readily reverse (Movie S 4). We also observed filaments sometimes reversing without reaching track ends, but when this happened, it was usually at a consistent location (Movie S 5). Thus, the observed tracks under phase microscopy constitute a confining space that bounds filament movement and dictates reversal frequency (see Supplementary Fig. S1). On glass, filament movement still displayed a go-stop-go pattern, but stoppages led to both reversal and re-orientation of direction, and as a result, trajectories were less well-defined (Movie S 2). Reversals and re-orientations were associated with a deceleration-acceleration of filaments (Fig. 1B and Fig. S2). This decrease in speed can either result from external counter forces or loss of propulsive force, for example through filaments reaching end of a track under agar or detaching from the surface on glass (see Movies S 1 and S 2). Mean filament speeds ranged around one micron per second and displayed a relatively uniform distribution against filament length (Fig. 1C and Fig. S3). Since filament speed results from a balance between propulsive forces and drag, the latter observation shows that all (or a fixed proportion of) cells in a filament contribute to propulsive force generation.
We can consider the reversal / re-orientation events of a filament akin to tumbling events seen in flagella-based bacterial motility. To this end, we were interested in the distribution of time spent during reversals, which we call the “dwell time”. Analysing over 1700 reversal events across 65 filaments, we found that dwell times display a long-tailed distribution (Fig. 1C). On glass, we found a similar long-tailed distribution from 256 reversal / reorientation events(Fig. S3). Most reversals happened within few seconds, while a minority involved significant time spent stationary. We found that there was no clear difference between dwell time of leading and trailing ends of the filament (Fig. S4). Taken together, these findings show that filament reversals are mostly rapid events and that they don’t involve delays across the filament.
A biophysical model captures filament movement dynamics and highlights a key role for cell-to-cell coordination for achieving fast reversals
As summarised in the introduction, molecular basis of propulsive forces are not fully understood. However, the presented quantification of motility dynamics highlights three key observations that can be used to develop a plausible biophysical model of filament movement. Firstly, and most obviously, cells in a filament remain physically linked during motion and therefore must exert mechanical forces onto each other. Secondly, all or most cells contribute to propulsive force generation, as seen from a uniform distribution of mean speed across different filament lengths. Thirdly, and finally, reversals are associated with reduction in speed (or stoppage) and usually happen within seconds.
We implement the first two observations by modelling cells in a filament as beads connected by springs, with each cell capable of generating propulsion and exerting pull and push forces onto others (see Fig. 2A and Methods). The third observation shows that cells have a mechanism to alter the direction of their propulsive force. To this end, we model reversals as resulting from a stochastic cellular signalling pathway that is linked to mechano-sensing of external and internal forces acting on a given cell (see Methods and Eqs. (3),(6)). This modelling choice effectively assumes that cells have an intrinsic reversal rate, but adjust this in such a way to reduce any compression or extension that they perceive. In addition, and to account for cellular sensing dynamics, we introduce a memory (i.e. a refractory period) that limits cells ability to switch force direction if they have recently done so (see Eqs. (3),(4) in Methods).
Given this model, we simulated movement of filaments in a 1-dimensional space and in the presence of an external force field. The latter allows us to implement a reduction of speed and analyse filament movement on a defined track, as observed experimentally for the movement under agar (Fig. 1). For appropriate choice of parameters, we found that this model can produce sustained back-and-forth movement (Fig. 2B and Movie S 6). The ability of the model to present reversals can be understood from the way it implements mechano-sensing at the cell level. When leading cells reach the track ends, the external force field causes an increase in their tendency to reverse the direction of their propulsive force (Fig. 2B). When this happens, trailing cells are compressed by both those in front and behind them, and, if cell-cell coupling is strong enough, are forced to reverse as well. Thus, the mechano-sensing creates a signal that travels across the filament and sustains a coordinated reversal event. The same dynamics can also be created from stochastic reversal of few individual cells, and as such, we observed also reversals in the model simulations without the filament reaching the track ends, for some parameter values (Fig. 2B).
To better understand how the key model parameters affect reversals, we run simulations across a wide range of parameter values (Fig. 2C and (2D). To quantify reversal dynamics, we considered two measures that are based on cell-to-cell coordination and extent of reversals. The reversal-related measure, M, gives the ratio of the number of actual reversals over the estimated number of reversals obtained by considering an idealised filament, that has a defined mean speed and reverses deterministically and instantaneously at track ends (Eq. (7). Thus, a value of M ≈ 1 can indicate regular and rapid reversals at track ends, however, it can also indicate simply a high reversal frequency. The coordination measure, S, gives the extent of homogeneity of cellular propulsive forces across a filament (Eq. (9)), with a high value indicating cells behaving in synchrony throughout their movement. Evaluating S and M together allows us to better quantify coordination in reversal dynamics along a track. We find that both measures are sensitive to the strength of the cell-to-cell coupling and the memory parameter (Fig. 2C). Indeed, when both parameters are reduced, reversals become more erratic, with filaments undergoing many reversals and individual cells having opposing propulsive force directions (Fig. 2B and Movie S 7). As expected from these results, these two parameters also affect the shape of the dwell time distribution resulting from model simulations (Fig. S5). For sufficiently high cell-to-cell coupling and memory parameters, simulated filaments reverse rapidly, and primarily at track ends, and cellular propulsive force directions remain homogeneous throughout their movement (Fig. 2B). To explore if real filaments’ movement conform to such expected dynamics arising from high cell-to-cell coupling, we quantified the number of expected reversals using individual filament and track lengths, and the observed mean filament speed from all observations. We found that for the high majority of the filaments, the number of observed reversals matches closely with the expected number of reversals, in qualitative agreement with the model simulations performed with sufficiently high cell-to-cell coupling parameter (Fig. S6).
Filament movement shows high cell-to-cell coordination and coupling of translation with rotation
Taken together, the experimental and model results presented so far show that most reversals are rapid events, and to achieve them the filaments should have high cell-to-cell coupling in their movement. To explore this model-highlighted idea of cellular coordination and to better understand what happens during a reversal at the single cell level, we undertook Total Internal Reflection Fluorescence (TIRF) microscopy, which illuminates only an ≈ 200 nm thick band of the sample and allows observation of surface features (Fig. 3A-B, top panels). This allowed us to clearly identify individual cells and their septa with neighboring cells, which we tracked to quantify individual cell movement. We found that cell movement across the filament mostly stays coordinated during a reversal, in terms of speed and direction of movement (Fig. 3A, bottom panels, see also associated Movie S 8).
In addition to identifying individual cells, TIRF microscopy revealed membrane-bound protein complexes localised predominantly near the poles of each cell (Fig. 3B, top panel). We have also tracked these complexes and found that their trajectories match 2-dimensional projections of a 3-dimensional helical trajectory, showing that filaments rotate as they translate (Fig. 3B, top panel and Movie S 8). The angle θ, that the trajectories make with the long axis of the filament shows a normal distribution with a mean of 21.2 degrees, while the angular speed of the protein complexes are well correlated with the linear velocity of the filament (Fig. 3B). These findings indicate either a helical force generation or a screw-like structure on the cell exterior that causes a strong coupling between rotation and translation. It is interesting to note in this context, that staining filaments with fluorescently labelled Concanavalin A (see Methods), which binds to slime, results in an helical pattern along the filament long axis (Fig. S7). In time-lapse TIRF movies, cellular protein complexes can be seen moving through this fluorescence pattern, suggesting that the slime surrounds the filaments as a tubular structure, fixed in position relative to filament movement (Movie S 9). The tubular nature of slime is further supported by scattering electron microscopy (SEM) images, which show features around the filaments that look like collapsed slime tubes Fig. S8.
Reversals in longer filaments can lead to buckling and plectoneme formation, causing loss of cellular coordination
The above findings show that filaments display high level of cell-to-cell coordination, supporting the model result that cellular coordination is one of the requirements for rapid reversals. In addition, we find a strong coupling between rotation and translation. Thus, it is possible that any differences in cell movement, for example the leading cells initiating reversal earlier, could lead to torsional forces developing along the filament. Supporting this, we note that cellular speeds display a steep acceleration immediately after a reversal (Fig. 3A), a charactheristic observed in all filaments analysed with TIRF. In a few filaments, we have also observed indication of leading cells initiating reversal earlier than trailing ones (Fig. S9 and Movie S10). To further explore this possibility of de-coordinated movement dynamics along the filament, and emergence of associated torsional forces, we sought to analyse longer filaments with the expectation that such dynamics are expected to arise more readily in longer filaments. We found that longer filaments, where observed, readily display buckling, and sometimes these buckled loop regions twist upon themselves and form plectonemes. Almost all of the 61 observed cases of buckling and plectoneme formation involved filaments longer than 400 microns (Fig. S10).
Given that plectonemes can form, we next wanted to understand cell-to-cell coordination during such events. While it was not possible to capture a filament with TIRF microscopy, during a plectoneme formation event, we were able to record whole-filament movies at low magnification. A representative case is shown on Figure 4 (and associated Movie S 11), where we were able to analyse movement dynamics for each end of the filament. We found that the filament initially moves in a coordinated manner, on a well-defined trajectory and the two ends of the filament being coordinated in speed (Fig. 4B). With the start of buckling and plectoneme formation, however, this coordination gets disrupted, with one end of the filament moving faster than the other, and finally reversing independently of the other end (Fig. 4B). This de-coordination seems to resolve towards the end of the recorded movement sequence, when the plectoneme resolves (Fig. 4 and Movie S 11). This representative case exemplifies the fact that buckling, and more specifically plectoneme formation, can result in de-coordination of cells within a filament and independent movement of filament ends. We have observed more extreme cases of such dynamics on glass, where plectonemes result in complete de-coordination in filament movement (Movie S 12). We also observed that plectonemes can get “dragged” and sometimes entangle other filaments (Movie S 13), a dynamic process that could be relevant for initiating aggregation of multiple filaments.
Discussion and Conclusions
Here, we analysed the motion dynamics of the filamentous, gliding cyanobacterium F. draycotensis. We show that gliding involves movement of filaments on curved or straight trajectories, interspersed with re-orientation or reversal events. Filament speed decreases and then increases around a reversal event, with mean speed around one micron per second, independently of filament length. We find that the time spent in reversal displays a long-tailed distribution, with most reversals taking few seconds, but some involving considerable time of filaments being stationary. We show that these experimental observations can be recapitulated by a physical model that assumes individual force generation by each cell, mechanical-coupling among cells, and mechano-sensory (coordinated) control of force direction. This model suggests that such coordination of propulsion direction among cells in a filament is important for being able to achieve fast reversals. In line with this insight, we find that cells in a filament remain highly coordinated during reversals and that their translational movement is coupled with rotation along the long axis of the filament. Following from the latter observation, we find that longer filaments can readily buckle, and buckled regions can twist upon themselves to form plectonemes. The formation of plectonemes results in de-coordination of cellular movement and independent dynamics for filament ends. We hypothesize such dynamics to be linked with the loss of sensory coupling among cells, and de-coordination of the direction of their propulsive forces. Taken together, these results provide a useful characterisation of gliding motility, which will inform future studies on both molecular mechanisms of force generation and collective behaviours of filaments.
In terms of molecular mechanisms of propulsive force generation, our findings that rotational and translational movement are closely coupled in F. draycotensis and that slime presents itself as a helical tubular structure around filaments could be highly relevant. Rotation, as well as its absence, has been commented on in several species of gliding, filamenteous cyanobacteria, but not quantified before [27, 30], while it has been highlighted as a key dynamical feature in the propulsive force generation in gliding, single-celled bacteria [24, 18]. Helical patterning of slime, as observed here, as well as surface fibrillar arrays, were described in several filamentous cyanobacteria [27, 30, 33]. Taken together, these observations point to the possibility that propulsive force generation involves helical tracks surrounding the cell and either motor complexes moving on them, as suggested for single-celled, gliding bacteria [24, 18] or pili being pushed across them. Recent studies from Nostoc species indicate propulsion to result from pili extension [6, 32], but whether that species presents filament rotation or fibrillar arrays is unclear [30].
Photoresponses in the motility behaviour is shown in several gliding filamentous cyanobacteria [7, 8, 9, 10], but a fully quantitative analysis of these at the single filament level is currently lacking. In this context, an interesting observation we made in this study is the presence of fluorescent, membrane-bound protein complexes. We hypothesize that these complexes are associated with light sensing, since our tracking of them in TIRF microscopy always coincided with reversal of filaments. Thus, it is possible that the wavelength at which they get excited (437nm laser used in TIRF microscopy) is also acting as a sensory signal for reversals. This possibility will be explored in a future study, along with further characterization of any phototaxis behavior in F. draycotensis.
Gliding motility in filamentous cyanobacteria is commonly associated with aggregate or macro-scale structure formation. These formations are relevant in the context of physiological functions of the cyanobacteria, and their associated microbial communities, as shown in the case of F. draycotensis [13]. Similarly, in the case of the ocean-dwelling species of the Trichodesmium genus, aggregate formation is linked with the biogeo-chemically relevant processes of iron acquisition [14]. The results presented here show that de-coordination of cellular forces within a filament can lead to formation of plectonemes, which can entangle multiple filaments together. Further studies of this process can elucidate how aggregate and macro-scale formations develop and inform their physical modelling.
Given that gliding motility is observed in phylogenetically diverse bacteria, it is possible for it to be under-pinned by different molecular mechanisms and present different motion dynamics. The inter-connectivity of propulsion and slime secretion is difficult to disentangle purely by molecular approaches, and development of physical models solely based on identified proteins can be limited. In particular, drastically different physical models can be constructed from observed molecular constituents. We advocate coupling of molecular studies with detailed analysis of motion dynamics, as quantification of these can provide an additional framework to constrain and test different models, ultimately providing a more complete understanding of gliding motility and associated collective behaviors.
Methods
Sampling and Culturing
Cultures of F. draycotensis were grown in medicinal flasks, in volumes of 30 mL BG11+ with added vitamin mix as described previously [13]. Cultures were kept under visible light irradiance of ≈ 30 µM·photons·m−2s−1, across a repeating 12-hour light/12-hour dark cycle. Samples were obtained from cultures aged between 4 and 8 weeks, by withdrawing 0.5 mL of material containing several pieces of biofilm, transferring the aliquot into a 1.5 mL microfuge tube and vortexing for 1 minute to break up the material. If the mixture was not homogeneous, additional mixing was performed using a 1 mL mechanical pipette.
Sample Preparation for Imaging
Prior to imaging, collected innoculate was diluted into fresh BG11+ media with vitamin mix to a desired dilution (typically 1:50). For microscopy performed with 1.5 − 2% agarose pads [37], 3µL of the homogenised diluted sample was added to the pad. The droplet was then left to dry for 10-15 minutes, depending on the age of the agarose pad and room humidity. Where stated, the agarose pad was flipped onto a 22×50 mm cover glass, sandwiching the sample between the agarose and the glass. Otherwise, the sample was left on top of the agarose and the pad placed in a cover glass bottomed dish (MatTek/CellVis) for microscopy.
As described, some samples were imaged inside a millifluidic chamber device (Ibidi µ-Slide III 3D Perfusion) or on glass slides, using a 1.5 × 1 cm Gene frame (Thermo Scientific) as a spacer. For glass slide samples, the Gene frame sticker was placed on the standard glass slide. To the frame 200 µL of diluted sample was added, and the chamber was sealed with a cover glass.
Regardless of sample preparation method, all prepared samples were left to settle for 30-60 minutes, being kept on a laboratory bench under room temperature and illumination.
Imaging Conditions
Light microscopy was performed on an Olympus IX83 microscope, using a CoolLED pE-300 illumination system. Imaging was performed using a CoolSnap-HQ2, with an imaging interval of 5-10 seconds or longer. The stage is enclosed in an environmental chamber and stabilised at 25°C. Phase imaging was performed with a low exposure rate (≈ 5 ms). Chlorophyll fluorescence was measured through a Chroma DSRed filter cube, with a recorded excitation light power of 17.4 mW and an exposure time of 100 ms. TIRF imaging was performed on a standard NanoImager system (ONI) using 473nm excitation laser and emission filters centered at 525nm and 625nm.
Image Analysis
Images were analysed with custom-written software in Python. Analyses involved tracking individual filaments, cells, or membrane complexes on each cell, across time-lapse images. In the case of filament tracking, movies of a larger field of view were cropped to focus on single filaments, and these were processed by fitting a contour on the filament and deriving a mid-point spline from this. By repeating this process across each image in a time-lapse, it was possible to track the position of middle, trailing, and leading ends of the filament and their speed. In the case of longer time-lapse movies, any drift in the images was stabilised using stationary filaments within the larger field of view. All code used in these analyses are made available on the Github repository; https://github.com/OSS-Lab/Cyano_Gliding_ImageAnalysis_and_Model.
Biophysical model
The model code is available at the Github repository;github.com/OSS-Lab/Cyano_Gliding_ImageAnalysis_and_Model. We model cyanobacteria as out-of-equilibrium filaments. For simplicity, we consider only a 1-dimensional model, but extension to higher dimensions would be straightforward. The filament is composed by Nf beads that represent cells in a cyanobacteria filament. The approach to model the filament as connected cells with springs between them is similar to several other physical models of active filaments and polymers[38, 39, 40]. The position of bead i is denoted by xi (i = 1..Nf). Beads are connected with their nearest neighbours, along the filament, by harmonic springs. Each bead follows an over-damped equation of motion, given by;
where γ is the friction coefficient and ξi(t) is a Gaussian white noise, following ⟨ξi(t)ξj(t′)⟩ = 2kBT/γδijδ(t − t′), T being the temperature. The potential energy V (x) reads
where β = 1/kBT is the inverse of the thermal energy (kB is the Boltzmann constant) and K is the spring constant, representing the elastic response of the cell upon compression or extension. The rest length l0 =1 represents the unperturbed (average) length of a cell. The remaining term in Eq. (1) represent the beads’ self-propulsion, i.e. there is a force , generated by the bead itself, that causes it to move in the positive or negative direction. The is defined as a constant in time, while si(x, t) sets the “active state” of bead i at position x and time t. We set si(x, t) to be a continuous variable, such that −1 < si(x, t) < 1. Thus, a cell can be propelled in the positive (si(x, t) ≈ +1) or negative (si(x, t) ≈ −1) x direction or it can be inert (si(x, t) ≈ 0).
To model the process of reversals, we envision that each cell can decide on the direction of its propulsive force. This idea is motivated by the fact that all studied uni-cellular gliding bacteria incorporate directional change, and individual cells in filamentous cyanobacteria are observed to change the location of their Type IV pili apparatus across the cell poles [32, 6]. To model cellular direction change, we introduce a function, ωi(x, t), describing the internal “state” of each cell and dictating its active state by the relation si(x, t) = tanh(ωi(x, t)). From a biological perspective, ωi(x, t) abstract the behavior of a cellular signalling network that integrates external and internal inputs for the cell and sets its propulsive direction. To this end, we define ωi(x, t) as integrating several inputs, as follows:
where γω is a friction coefficient, ηi(t) is a Gaussian white noise, Uc is a confinement potential, fext(x, t) is an (adimensional) external force field, and Vω is a mechano-sensing (adimensional) potential. The confinement potential is a mathematical abstraction that allows us to maintain the output of ωi(x, t) in a finite interval. This potential is given by:
with A = 1. The parameter ωmax determines the extreme values that ωi(x, t) can attain. As such, this parameter acts as a memory term; when given terms in Eq. (3) push it in one direction, ωmax will determine the time when their effects will start changing ωi(x, t).
The noise term captures any fluctuations, either internal to the cellular signalling networks or caused by external noise (thermal or otherwise). It ensures that the individual cells can also turn on or off or switch direction spontaneously. The noise term ηi(t) follows ⟨ηi(t)ηj(t′)⟩ = 2Dωδijδ(t − t′), with Dω = kBTω/γω a diffusion coefficient, and Tω being the temperature that characterises the magnitude of the fluctuations of ω in absence of external stimuli.
Following the experimental observations, we postulate the existence of an external, mechanical stimulus that keeps filaments on a well-defined trajectory on agar, i.e. a track. We assume for simplicity that the track is fixed and extends from x/l0 = 0 to x/l0 = L. The external stimulus field is given by the following functional form:
where fext is an adimensional constant. Given Eq. (3), the action of the external field drives the auxiliary variable to positive/negative values.
The potential Vω encodes the mechano-sensitivity of each cell to their neighbors and is given as;
where the quantity Kω = fext/l0 + δKω is an effective spring constant. The function Vω encodes the feed-back to a mechanical stress (compression/extension) applied on the cell, compelling the direction of the self-propulsion so that this stress is relieved.
Notice that, within the model, fext essentially encodes the basal strength of the response of the cells to stimuli, either extra-or inter-cellular. The mechano-sensing can be more (δKω > 0) or less (δKω < 0) relevant than the external stimulus; naturally, we impose δKω > −|fext|, i.e. there can not be a positive feedback to compression.
Simulation details
We integrate the equations of motion Eq. (1) using a stochastic predictor-corrector algorithm, while the equations Eq. (3) are integrated using the Euler-Maruyama algorithm. In both cases, we employ an integration time step Δt = 10−3τ, τ being the unit of time. As mentioned above, we set the average length of an unperturbed cell as the unit of length and set the unit of time as the diffusion time of a single inactive cell . For convenience, we set the unit of energy as the active “work” ϵ = f al0 = 1; as such, we fix f a =1 throughout the simulations. We further fix the friction coefficient to . We assume that the motion of the filament is dominated by its self-propulsion, i.e. f al0 ≫ kbT and we set kBT = 0.01ϵ.
Concerning the parameters that enters in the evolution of ω, the scale of the “energy” should be set here by kBTω. However, we set kBTω = 0.1 for convenience; we find that reasonably smooth reversals can be obtained for such a choice, maintaining the other parameters within the same order of magnitude. Indeed, we fix fext = 10.
Movement characterisation
The analysis of experimental and simulated trajectories focuses on the dwell time and the number of reversals. The dwell time is computed as follows: the trajectory is smoothed via a spline and the velocity is computed via numerical differentiation. A dwell time is defined as the interval within which the velocity remains, in absolute value, below 15% of its maximal value.
A reversal corresponds to a change of sign of the centre of mass velocity of the filament; we name nr the number of reversals counted along a trajectory. We define two quantities: the reversal efficiency M and the reversal rate ν. The former quantity is defined as
With and τm the measurement time (either in experiments or simulations). The quantity is the number of reversals one would expect from an ideal filament, that reverses smoothly and deter-ministically in a negligible time and then travels at constant speed. In simulations, v = f a/γ = 1l0/τ while in experiments we estimate v ≈ 1µm/s. The quantity τ (e) is, thus, the ideal time between two reversals, i.e. the time the centre of mass would take to travel from its initial condition (after a reversal) to the track boundary. The reversal rate is simply defined as the number of reversals over the measurement time
Finally, we define a syncronisation index S as
where the average is performed over time or over the different realisations. Essentially, for each conformation, we compute the (instantaneous) average propulsion direction; as we are interested in the overall synchronisation of the filament regardless of its direction, we take the absolute value.
Funding
This project is funded by the Gordon and Betty Moore Foundation (grant https://doi.org/10.37807/GBMF9200). E. L. acknowledges support from the MIUR grant Rita Levi Montalcini. C.V. acknowledges funding from MINECO grants IHRC22/00002 and PID2022-140407NB-C21. M.P. acknowledges funding from the Lever-hulme Trust (Grant number RPG-2018-345) and the fact that IMEDEA is an accredited “María de Maeztu Excellence Unit” (Grant CEX2021-001198, funded by MCIN/AEI/10.13039/ 501100011033).
Acknowledgements
We thank Douglas Risser, Jonasz Slomka, and Rebecca Poon for insightful comments on an earlier version of this manuscript. We acknowledge the help of imaging facility manager Ian Hands-Portman with microscopy setups and Nicole Robb for allowing access to her groups’ TIRF microscopy. We would also like to acknowledge the University of Warwick Electron Microscopy RTP for assistance in the research described in this paper.
Supplementary Information
10.1 Supplementary videos
1. Movie S1- Filament moving under 1.5 % agarose sandwich. Magnification of ×4 (1.61 µm/pixel). Time interval of 5 seconds. Fluorescence recorded in Red channel with 100 ms exposure. AVI produced with 20 fps compression. Link: https://youtu.be/RsANG2RBzTg
2. Movie S2- Filaments moving across glass slide. Magnification of ×10 (0.645 µm/pixel). Time interval of 10 seconds. Fluorescence recorded in Red channel with 200 ms exposure. AVI produced with 20 fps compression. Link: https://youtu.be/KOcwKboCFZg
3. Movie S3- Single filament in track on 2.5 % agarose. Magnification of ×10 (0.645 µm/pixel). Time interval of 5 seconds. Phase image with 5 ms exposure. AVI produced with 20 fps compression. Link: https://youtu.be/XVIv4FLYBas
4. Movie S4- Single filament on a circular trajectory on 1.5 % agarose. Magnification of ×4 (1.61 µm/pixel). Time interval of 5 seconds. Fluorescence recorded in Red channel with 100 ms exposure. AVI produced with 20 fps compression. https://youtu.be/PJdFEE6R6fk
5. Movie S5- Several filaments in track on 2.5 % agarose. Magnification of ×10 (0.645 µm/pixel). Time interval of 5 seconds. Phase image with 5 ms exposure. AVI produced with 20 fps compression. Link: https://youtu.be/hu-U810RMuw
6. Movie S6- Simulation movie of an active filament with a parameter set, resulting in coordinated reversals. The filament is composed of 8 beads (cells) in a track of length 50 units. AVI produced with 5 fps compression. Link: https://youtu.be/nyUmk8QXz1I
7. Movie S7- Simulation movie of an active filament with a parameter set, resulting in erratic reversals. The filament is composed of 8 beads (cells) in a track of length 50 units. AVI produced with 5 fps compression. Link: https://youtu.be/Llc5oOmdxm0
8. Movie S8- TIRF imaging of filament undergoing a reversal. This movie is associated with the analysis shown in (Fig. 3. Magnification of ×100 (0.116 µm/pixel). Time interval of 33 ms. Fluorescence in cyan is from the bundled complexes on the filament surface, observed between 500-560 nm, with 30 ms exposure from a 473 nm laser. Fluorescence in the red channel comes from chlorophyll fluorescence, with 30 ms exposure from a 473 nm laser and observed in a channel between 600-670 nm. AVI produced with 100 fps compression. Link: https://youtu.be/Doa-jOkl04I
9. Movie S9- TIRF imaging of stacked filament in a stained slime sheath. Magnification of ×100 (0.116 µm/pixel). Time interval of 33 ms. A 473 nm laser with an exposure time of 30 ms was used for excitation. The slime sheath is stained with concanavalin AlexaFluor-488 and coloured in yellow (observed 500-560 nm). Chlorophyll autofluorescence is coloured in cyano (observed 600-670 nm). AVI produced with 60 fps compression. Link: https://youtu.be/mYoqSqK9re0
10. Movie S10- TIRF imaging of filament undergoing a reversal, where fluorescence shows the protein complexes on the surface, This movie is associated with the analysis shown in (Fig. S9. Magnification of ×100 (0.116 µm/pixel). Time interval of 33 ms. Fluorescence observed between 500-560 nm, with 30 ms exposure from a 473 nm laser. AVI produced with 100 fps compression. Link: https://youtube.com/shorts/pa5iIcua8Ps?feature=share
11. Movie S11- Buckling filament under 1.5 % agarose. This movie is associated with the analysis in (Fig. 4 of the main text. Magnification ×4 (1.61 µm/pixel). Time interval of 5 seconds. Fluorescence recorded in Red channel with 100 ms exposure. AVI produced with 20 fps compression. Link: https://youtu.be/EmfIGJiiWUg
12. Movie S12- Filament twisted in middle, forming plectonemes on glass. Magnification of ×4 (1.61 µm/pixel). Time interval of 30 seconds. Fluorescence recorded in Red channel with 100 ms exposure. AVI produced with 20 fps compression. Link: https://youtu.be/VzHMagiGqfo
13. Movie S13- Filaments forming plectonemes on glass. Magnification of ×4 (1.61 µm/pixel). Time interval of 5 seconds. Fluorescence recorded in Red channel with 100 ms exposure. AVI produced with 20 fps compression. Link: https://youtu.be/OBm0IpUN6rI
10.2 Supplementary figures
References
- [1]Efficient CO2 fixation by surface Prochlorococcus in the Atlantic OceanISME J 8:2280–9https://doi.org/10.1038/ismej.2014.56
- [2]The role of biology in planetary evolution: cyanobacterial primary production in low-oxygen Proterozoic oceansEnviron Microbiol 18:325–40https://doi.org/10.1111/1462-2920.13118
- [3]Microbial mats: physiological ecology of benthic microbial communities
- [4]Hormogonium Development and Motility in Filamentous CyanobacteriaAppl Environ Microbiol 89https://doi.org/10.1128/aem.00392-23
- [5]Cyanobacteria in motionCurr Opin Plant Biol 37:109–115https://doi.org/10.1016/j.pbi.2017.03.018
- [6]Motility in cyanobacteria: polysaccharide tracks and Type IV pilus motorsMol Microbiol 98:998–1001https://doi.org/10.1111/mmi.13242
- [7]Collective self-caging of active filaments in virtual confinementarXiv https://doi.org/10.48550/arXiv.2403.03093
- [8]The involvement of type IV pili and the phytochrome CphA in gliding motility, lateral motility and photophobotaxis of the cyanobacterium Phormidium lacunaPLoS One 17https://doi.org/10.1371/journal.pone.0249509
- [9]Genetic analysis reveals the identity of the photoreceptor for phototaxis in hormogonium filaments of Nostoc punctiformeJ Bacteriol 197:782–91https://doi.org/10.1128/JB.02374-14
- [10]Photomovement in motile microorganisms–IIPhotochem Photobiol 47:837–69https://doi.org/10.1111/j.1751-1097.1988.tb01668.x
- [11]Undirected motility of filamentous cyanobacteria produces reticulate matsGeobiology 8:179–90https://doi.org/10.1111/j.1472-4669.2010.00235.x
- [12]Cellular Dynamics Drives the Emergence of Supracellular Structure in the Cyanobacterium, Phormidium sp. KSLife 4:819–36https://doi.org/10.3390/life4040819
- [13]Niche formation and metabolic interactions result in stable diversity in a spatially structured cyanobacterial communitybioRxiv https://doi.org/10.1101/2022.12.13.520286
- [14]Selective collection of iron-rich dust particles by natural Trichodesmium coloniesThe ISME journal 14:91–103https://doi.org/10.1038/s41396-019-0505-x
- [15]Gliding motility in cyanobacterial: observations and possible explanationsArch Microbiol 174:11–7https://doi.org/10.1007/s002030000187
- [16]Bacterial gliding motility: multiple mechanisms for cell movement over surfacesAnnu Rev Microbiol 55:49–75https://doi.org/10.1146/annurev.micro.55.1.49
- [17]Gliding motility revisited: how do the myxobacteria move without flagella?Microbiol Mol Biol Rev 74:229–49https://doi.org/10.1128/MMBR.00043-09
- [18]Bacterial Gliding Motility: Rolling Out a Consensus ModelCurr Biol 27:R154–R156https://doi.org/10.1016/j.cub.2016.12.035
- [19]Bacterial surface translocation: a survey and a classificationBacteriol Rev 36:478–503https://doi.org/10.1128/br.36.4.478-503.1972
- [20]Motor-driven intracellular transport powers bacterial gliding motilityProc Natl Acad Sci U S A 108:7559–64https://doi.org/10.1073/pnas.1101101108
- [21]Myxobacteria gliding motility requires cytoskeleton rotation powered by proton motive forceProc Natl Acad Sci U S A 108:2498–503https://doi.org/10.1073/pnas.1018556108
- [22]Helical flow of surface protein required for bacterial gliding motilityProc Natl Acad Sci U S A 110:11145–50https://doi.org/10.1073/pnas.1219753110
- [23]Flagella stator homologs function as motors for myxobacterial gliding motility by moving in helical trajectoriesProc Natl Acad Sci U S A 110:E1508–E1513https://doi.org/10.1073/pnas.1219982110
- [24]The mechanism of force transmission at bacterial focal adhesion complexesNature 539:530–535https://doi.org/10.1038/nature20121
- [25]Comparative transcriptomics with a motility-deficient mutant leads to identification of a novel polysaccharide secretion system in Nostoc punctiformeMol Microbiol 87:884–93https://doi.org/10.1111/mmi.12138
- [26]Identification of a hormogonium polysaccharide-specific gene set conserved in filamentous cyanobacteriaMol Microbiol 114:597–608https://doi.org/10.1111/mmi.14566
- [27]Gliding in a blue-green alga: a possible mechanismNature 225:1163–5https://doi.org/10.1038/2251163a0
- [28]Envelope Structure of Four Gliding Filamentous CyanobacteriaJournal of Bacteriology 177
- [29]Fibrillar Array in the Cell Wall of a Gliding Filamentous CyanobacteriumJournal of Bacteriology 181
- [30]Nanoscale visualization of a fibrillar array in the cell wall of filamentous cyanobacteria and its implications for gliding motilityJ Bacteriol 189:7361–6https://doi.org/10.1128/JB.00706-07
- [31]The junctional pore complex, a prokaryotic secretion organelle, is the molecular motor underlying gliding motility in cyanobacteriaCurr Biol 8
- [32]Evidence that a modified type IV pilus-like system powers gliding motility and polysaccharide secretion in filamentous cyanobacteriaMol Microbiol 98:1021–36https://doi.org/10.1111/mmi.13205
- [33]Type two secretion systems secretins are necessary for exopolymeric slime secretion in cyanobacteria and myxobacteriabioRxiv https://doi.org/10.1101/2022.03.08.483542
- [34]Dynamic localization of HmpF regulates type IV pilus activity and directional motility in the filamentous cyanobacterium Nostoc punctiformeMol Microbiol 106:252–265https://doi.org/10.1111/mmi.13761
- [35]Controlled motility in the cyanobacterium Trichodesmium regulates aggregate architectureScience 380:830–835https://doi.org/10.1126/science.adf2753
- [36]Active Spaghetti: Collective Organization in CyanobacteriaPhys Rev Lett 131https://doi.org/10.1103/PhysRevLett.131.158303
- [37]Preparation of agarose pads suitable for viewing filamentous cyanobacterial microbial communities using time-lapse imagingprotocols.io https://doi.org/10.17504/protocols.io.81wgb65p1lpk/v1
- [38]Self-propelled worm-like filaments: spontaneous spiral formation, structure, and dynamicsSoft matter 11:7181–7190
- [39]The physics of active polymers and filamentsThe journal of chemical physics 153
- [40]Globulelike conformation and enhanced diffusion of active polymersPhysical review letters 121
Article and author information
Author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
Copyright
© 2024, Rosko 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
- views
- 85
- downloads
- 0
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.