Random walk and cell morphology dynamics in Naegleria gruberi

Amoeboid cell movement and migration are wide-spread across various cell types and species. Microscopy-based analysis of the model systems Dictyostelium and neutrophils over the years have uncovered generality in their overall cell movement pattern. Under no directional cues, the centroid movement can be quantitatively characterized by their persistence to move in a straight line and the frequency of re-orientation. Mathematically, the cells essentially behave as a persistent random walker with memory of two characteristic time-scale. Such quantitative characterization is important from a cellular-level ethology point of view as it has direct connotation to their exploratory and foraging strategies. Interestingly, outside the amoebozoa and metazoa, there are largely uncharacterized species in the excavate taxon Heterolobosea including amoeboflagellate Naegleria. While classical works have shown that these cells indeed show typical amoeboid locomotion on an attached surface, their quantitative features are so far unexplored. Here, we analyzed the cell movement of Naegleria gruberi by employing long-time phase contrast imaging that automatically tracks individual cells. We show that the cells move as a persistent random walker with two time-scales that are close to those known in Dictyostelium and neutrophils. Similarities were also found in the shape dynamics which are characterized by the appearance, splitting and annihilation of the curvature waves along the cell edge. Our analysis based on the Fourier descriptor and a neural network classifier point to importance of morphology features unique to Naegleria including complex protrusions and the transient bipolar dumbbell morphologies.


Introduction
Combinatorial use of persistent motion and reorientation is a common feature found in cell movement.Be it bacterial swimming or amoeboid crawling, persistent movement allows cells to gain most distance in one preferred direction so as to facilitate efficient escape from hazards or conversely attraction to nutrients.Reorientation on the other hand is not only required to adjust direction of persistent movement but also to facilitate cells to randomly explore and survey uncertain extracellular environments (Viswanathan et al., 1999;Bartumeus et al., 2002).In E. coli bacteria, the cell movement consists of a period of straight run interrupted by a stall or "tumble" where flagellar rotation reverses and cells reorient in random directions.The frequency of tumbling is regulated through a chemosensory system so as to provide orientation bias towards an attractant or away from a repellent.The exact nature of such motility pattern determines how well E. coli cells disperse (Taktikos et al., 2013).In the amoeboid movement, pseudopodal protrusions enriched in branched F-actin networks (Pollard, 2007) are formed transiently and can guide cells in different orientations.In addition, a confined region of the plasma membrane needs to retract in order to realize net displacement.In many cell types, cortical F-actin that is crosslinked with myosin II is enriched in such contractile membrane regions (Chi et al., 2014).Persistent movement arises when a cell has mono-polarity meaning that it has a single dominating leading edge and a retracting trailing end.The occurrence and location of these organizational events along the plasma membrane determine the sequential appearance of plasma membrane protrusions and rear retractions, ultimately influencing the direction, speed, and duration of cell movements.
Quantitative time-series analyses of cell displacement and cell shape change are important for explicit characterization of random cell motion.In many cases, cell displacement can be approximated as a particle obeying persistent random walk.Phenomenologically, the simplest form of differential equation that describes such stochastic dynamics is the Langevin equation (Dunn and Brown, 1987;Selmeczi et al., 2005;Selmeczi et al., 2008) where v is the velocity vector, β is the decay rate, σ is the noise strength, and ξ t is 2D white Gaussian noise.Random walk of E. coli can be approximated by Brownian motion having a short-term memory.In eukaryotic crawling, cell trajectories of fibroblast cells (Gail and Boone, 1970) and endothelial cells (Stokes et al., 1991) are also consistent with this simple persistent random walk model.In many other cell types, random walk includes memory that depends on the velocity and orientation which can be described by modifications to the above model (Takagi et al., 2008;Li et al., 2011).There are also random walk statistics called Lévy-walk with step lengths that follows a long tail (power-law) distribution (Viswanathan et al., 1999).There, the Mean Square Displacement (MSD) essentially diverges and the trajectories are characterized by self-similarity of the step lengths (Reynolds, 2010;Reynolds and Ouellette, 2016).Because Lévy-walk has very small probability of revisiting the same location, it is thought to arise in systems such as bird foraging that require an efficient search strategy.At the cellularlevel, effector T-cells (Harris et al., 2012) and swarming bacteria have been reported to exhibit Lévy-flight like statistics (Ariel et al., 2015;Huo et al., 2021).
To date, quantitative understanding of random walk behavior of amoeboid cells is limited to data from a handful of cell-types; these are mostly timelapse microscopy images of cultured metazoan cells and amoebozoa Dictyostelium.From microbial ethology and evolutionary biology perspectives, however, we should note that amoeboid movement is found not only in animals, fungi and amoebozoans (Prostak et al., 2021) but also in largely uncharacterized species in the excavate taxon Heterolobosea namely Naegleria spp.and the slime mold acrasids (Brown et al., 2012).The ancestors of the opisthokont lineage and Naegleria diverged more than a billion years ago (Parfrey et al., 2011).Knowing the details of motility characteristics in Naegleria should help us understand the common design of the motility trait that is either deeply conserved across taxa or acquired independently by strong selective advantages.
Among members of genus Naegleria, non-pathogenic Naegleria gruberi (hereafter refer to as N. gruberi) is the better characterized species whose genome has been sequenced (Fritz-Laylin et al., 2010).In its amoebic phase, N. gruberi grows and divides by feeding on bacteria through phagocytosis (Fulton, 1970).Under low electrolyte conditions, it quickly shifts to the non-feeding flagellated state by rapid de novo synthesis of microtubules (Walsh, 2007).In the amoebic state, the overall cell cortex is enriched in F-actin with marked accumulation around membrane ruffles (Velle and Fritz-Laylin, 2020).An early work using reflection interference microscopy has revealed that N. gruberi adhere and form discrete dot-like contacts to non-treated glass surfaces and migrate (Preston and King, 1978).These so-called "focal contacts" leave behind footprints of membrane residues on the glass substrate as the cells crawl away (Preston and King, 1978).With the advent of genomics and molecular cell biology, it has become clear that N. gruberi possess the essential side-branching nucleator of F-actin-the Arp2/3 complex and its activators WASP and SCAR (Fritz-Laylin et al., 2017;Velle and Fritz-Laylin, 2020;Prostak et al., 2021).Inhibition of Formin reduces directional persistence, and inhibition of the Arp2/3 complex reduces the cell speed (Velle and Fritz-Laylin, 2020).N. gruberi also has Myosin II (Sebé-Pedrós et al., 2014) and a potential orthologue of Integrin beta (Morales et al., 2022), although whether they exist in other groups in Excavata is unclear (Velle and Fritz-Laylin, 2019).
While the above works indicate likely similarity of actindependent processes involved in cell crawling in an evolutionary distant eukaryote, quantitative characterization of the cell-level motility pattern is so far lacking.Do N. gruberi cells exhibit persistent random walk behavior?What is the characteristic time scale of persistence and reorientation if any? How similar are their movements compared to the well-studied systems such as Dictyostelium and immune cells?In this report, we performed quantitative analysis on cell movements and shape change of N. gruberi.Our analysis demonstrates that N. gruberi cells exhibit persistent random walk driven by a large morphology change that involves appearance, splitting and annihilation of uniquely complex pseudopodium protrusions.

An overview of cell movement and cell morphology
To quantitate the movement of N. gruberi on a two-dimensional flat surface, we performed phase contrast time-lapse microscopy.A non-coated glass coverslip was employed as a cell substrate throughout this study.Figure 1A shows representative phase contrast images of N. gruberi in liquid growth media (Materials and Methods).The cells under our culture condition exhibited one or more hyaline protrusions that appeared dark in phase contrast images (Figure 1A arrows).In the example shown, protrusions extended along the glass surface for 15-50 s and the one that became dominant (i.e. the leading edge) extended in the direction of the overall cell movement (Figure 1A, 0 s).Marked cytoplasmic streaming from the center of the cell towards these extensions was observed (Supplementary Movie S1).A new protrusion appeared and extended first in the lateral direction (Figure 1A 20 s, 60 s arrow) and steered towards the front.It was then bent sideway before being retracted (Figure 1A 40 s, 80 s).Duration of the pseudopod extension/retraction cycle varied between 15-50 s (Supplementary Figure S1; Supplementary Movie S2).Concomitant reversal of cytoplasmic streaming was observed during retraction of pseudopods.A small bud-like bulge at the trailing end of a cell which we shall refer to as "uroid" appeared as a residue of a retracted pseudopod that was retained for an extended period of time (Figure 1A white circle).The uroid contained thin filopodia-like projections as described earlier (Preston and King, 1978).
Under our culture condition, the cells appeared to re-orient in random directions at irregular timing.We performed long time cell tracking by employing an automated stage that was programmed to track target cells (see Methods).Figure 1B shows representative cell trajectories obtained from the automated tracking.The trajectories consisted of a period of straight movement that lasted for about 30-200 s and a time period of relative low displacement and reorientation (Figure 1B).The movement is thus, at surface, akin to the run-and-tumble behavior of E. coli.There was a close link between the run/re-orientation dynamics with the cell shape.During a straight run, cells took a fan-like shape (Figure 2A; Supplementary Movie S3).The tail remained narrow while the front was occupied by a broad lamellipodia that expanded then split into branches of pseudopods (Figure 2A, 0-16 s).These bifurcating protrusions often fused to restore a large lamellar extension (Figure 2A, 20 s).On the other hand, cells re-oriented when the bifurcated protrusions remained separate.In most cases, the uroid persisted during front splitting and thus the cells took a Y-or trident shape (Figure 2B; Supplementary Movie S4).There were also cases where the uroid disappeared in Y-shaped cells (Figure 2C; Supplementary Movie S5).The two fronts expanded in the opposing directions and gave rise to a transient "dumbbell-like" bipolar morphology (Figure 2C, 70 s).After 10 s, one end shrunk and became the uroid while the other end became the next front (Figure 2C, 80 s).There was little centroid displacement during this period which lasted for about 40 s.

Random walk statistics
To characterize the random-walk statistics, we quantified the mean square displacement (MSD) and the instantaneous speed defined by the centroid displacement in 1 s time interval from trajectories of N = 10 cells (Figure 3).Even with the help of automated stage tracking, fast movement of N. gruberi made it difficult to track cells for long time duration before they come close to the edge of the plate or collided with one another.Thus, to obtain MSD, single trajectories were each divided into sub-trajectories of 100-3,600 s time-window and treated as independent data samples (Figure 3A).The slope of the MSD from the individual trajectories was 1.5-2.0,where the mean and standard deviation are 1.77 and 0.08 (Figure 3B; Supplementary Figure S2A).The time-dependency of the MSD indicates that the random walk of N. gruberi falls somewhere between pure Brownian (exponent of 1) and ballistic (constant velocity) motion (exponent of 2). Figure 3C and Supplementary Figures S2B-F show the distribution of the instantaneous velocity.The distribution followed 2-dimensional Gaussian with zero-mean and standard deviation of 51 μm/min (0.86 μm/s) (Figure 3C, Supplementary Figures S2B-F).This feature is distinct from Dictyostelium random motility which is non-Gaussian (Takagi et al., 2008).The median of the absolute speed was 60 μm/min (1.0 μm/s) which is close to the average speed reported in earlier literatures (King et al., 1981;Thong and Ferrante, 1986).
The temporal auto-correlation of the centroid speed (velocity autocorrelation; VAC) (Supplementary Figure S2G) shows, on average, that there are two characteristic decay times that cross over at around 10 s (Figure 3D).By assuming that VAC follows the sum of two exponential function (Selmeczi et al., 2005) with velocity v: Protrusion dynamics and the cell shape change.(A) A fan-shaped cell with front splitting during a persistent run.(B) Front splitting followed by reorientation (curvature kymograph for the sequence is shown in Figure 5A).(C) Dumbbell shape arise after front splitting and disappearance of the uroid (curvature kymograph for the sequence is shown in Supplementary Figure S4E).Arrows: orientation of centroid movement.Inverted triangles: propagating curvature waves. (2) we obtained decay time T 1 and T 2 of approximately 6 s and 90 s, respectively, where the weight Φ 1 and Φ 2 are 0.36 μm 2 /s 2 and 0.87 μm 2 /s 2 (Figure 3D Red curve).Based on the Bayesian information criterion (BIC) (Schwarz, 1978), two exponential functions gave the lowest value compared to one or three (Supplementary Figure S2H).When the length of time sequence chosen for the analysis was doubled from 50 s to 100 s, deviation of the parameter values was within an order of magnitude (Supplementary Figure S2H magenta curve).Decay time T 2 of approximately 90 s was also evident from the time derivative of VAC (Supplementary Figure S2I).
To check for orientational preference in the memory, we plotted the relationship between velocity and acceleration (change in v in 1 s interval) (Supplementary Figures S2J, K).The mean acceleration orthogonal (a ⊥ ) to the velocity was near zero regardless of |v| (Figure 3E; green circle) with non-zero variance (Figure 3E green stars) which suggests that the orientation of N. gruberi has no apparent left-right asymmetry.On the other hand, the mean acceleration parallel (a ) to the velocity was near zero at small velocity then decreased towards the negative at large velocities (Figure 3E).The standard deviation (Figure 3E; blue stars) increased somewhat at high |v|, however rarity of these fast step events prevented us from obtaining reliable averages.These features of acceleration are similar to those reported for Dictyostelium (Takagi et al., 2008).The negative acceleration parallel to the migration direction at high |v| implies that the cells do not maintain high |v| during re-orientation.Accordingly, when we plot reorientation angle θ as a function of |v| (Figure 3F) we see that most of re-orientation occurs below |v| = 1 μm/ sec.Above 1 μm/s the cells are moving in a straight line; i.e. cosθ = 1.

Generalized Langevin equation
To gain further insights on the specifics of the random walk statistics, it is instructive to compare the data with the behavior of simple idealized equations.The velocity auto-correlation that follows the sum of two exponential indicates that random walk dynamics cannot be captured simply by the Ornstein-Uhlenbeck process (Eq. 1) which only has a single exponent (Dunn and Brown, 1987).A straight-forward and minimal extension to Eq. 1 is to include additional memory with the decay rate γ as an integral in the form of generalized Langevin-equation (Selmeczi et al., 2005) Here, α is the strength of memory effect, and ξ t is a normalized Gaussian white noise that satisfies 〉 is an ensemble average and δ(t) is the delta function, σ is the noise strength (Selmeczi et al., 2005).By introducing the equation of motion becomes Based on the values of T 1 , T 2 , Φ 1 , Φ 2 obtained above, we calculated the parameter values of the generalized Langevin Trajectories, the MSD and the VAC were obtained by numerically calculating Eqs 5a, b with the parameters obtained above (Table 1).The individual trajectories consist of combination of persistent movement and turns (Figure 4A).The slope of MSD had mean and standard deviation of 1.80 ± 0.06, which matched well with the experimental data (Figure 4B).The distribution of |v| showed a single peak that was slightly smaller compared to the experimental data (Figure 4C).The median was 56 μm/min (0.94 μm/s) in the simulation, which matched well with 60 μm/ min in the experiment.The velocity autocorrelation consists of two slopes that crossed over at around 10 s (Figure 4D red), which was similar to the crossover in the experimental data (Figure 4D black).Velocity dependence of acceleration also matched well with the experimental data (Figure 4E).On the other hand, the range of cell speed at which turning occurred in the simulations was somewhat broader (0-1.2 μm/s) compared to the real cell (0-0.8μm/s) (Figure 4F).While the MSD and the VAC characteristics were well captured by the memory effect described in Eq. 3, deviation from the model became evident when comparing autocorrelation separately for the centroid movement (absolute velocity | v|) and the orientation ( v/| v|) (Supplementary Figure S3).In the experimental data, it is only the autocorrelation of the orientation v/| v| not | v| that showed two decay times (Supplementary Figure S3A, B).In the generalized Langevin-equation, the velocity and the orientation share the same time scales, and thus the autocorrelation of both the orientation v/| v| and | v| decayed with the two exponents (Supplementary Figures S3C, D).

Cell shape dynamics
Rather than pursuing extensions of the particle-based formalism such as those that treat the two timescales separately (Li et al., 2008;Takagi et al., 2008), we sought to more directly characterize cell reorientation by analyzing the cell shape dynamics.Based on binarized cell mask images and a boundary tracking algorithm (Nakajima et al., 2016; see also Supplementary Figure S4A), 500 points along the edge of cell masks were tracked in the laboratory frame for the local curvature and the normal velocity (Figure 5A, B; see also Supplementary Figures S4B, D, F for another sample).A protruding edge can be seen as a positive local-maximum in the curvature (Figure 5A yellow regions).The advancing front of a cell can be discerned by its positive velocity (Figure 5B, yellow regions), and the trailing uroid by the negative velocity (Figure 5B, blue regions).At the cell front, a new protrusion frequently appeared to split off from a pre-existing protrusion (Figures 5A, B white  The curvature wave dynamics are surprisingly similar to those obtained for Dictyostelium and neutrophil-like HL60 cells (Driscoll et al., 2012;Imoto et al., 2021) with a noticeable difference that splitting was more frequent and thus numerous.The other difference compared to Dictyostelium and HL60 cells is the occasional and transient appearance of dumbbell-like shape (Supplementary Figures S4C, E, G; Supplementary Movie S5).When it appears, the centroid velocity orientation showed discontinuous change (Supplementary Figure S4C, black arrow).In the kymograph representation, a dumbbell-like cell shape appears as two or three stable curvature waves (Supplementary Figure S4E, black arrow).Most positions had zero velocity (Supplementary Figure S4G, black arrow), indicating stalling of cell shape change.
These observations indicate that as the dumbbell shape appeared, the cell stopped and randomized its orientation.There were also rare cases where the cell maintained mono-polarity for an extended period of time (Figures 5E-G; Supplementary Movie S6; see Supplementary Figure S5 for additional samples).There, new curvature waves emerged frequently and traveled fast before disappearing at the tail (Figure 5E).The position where curvature waves appeared always showed positive velocity, while the positions where curvature waves disappeared showed negative velocity (Figures 5E, F).These patterns in the kymograph correspond well with the observation of fast curvature waves that propagate from the advancing cell front and disappear at the uroid (Figure 5G).
A further analysis showed a close relationship between the curvature wave and the centroid movements.The protruding and the retracting membrane regions can be identified as positive curvature regions with positive (Figure 6A, white dots) or negative (Figure 6A black dots) velocity respectively.The orientation of the normal vector at the protruding region showed high correlation with the direction of centroid velocity (Figure 6B blue).The retracting regions oriented in the opposite directions which appeared somewhat broader in distribution (Figure 6B orange).To further analyze the dynamics of the curvature wave, high curvature regions (Figure 6C white) at each time frame were assigned as individual protrusions (Figure 6C green).While there were multiple protrusions in the protruding region, a dominant leading edge can be detected from identifying a single protrusion whose normal vector angle was the closest to that of the centroid velocity (Figure 6C magenta).Once a curvature wave became the leading edge, it remained so for about 2.8 s as measured from its average lifetime (Figure 6D).Another interesting feature of the membrane extensions is that they gave birth to secondary pseudopods or were steered to other directions.The typical angular velocity associated with this dynamic was 0.1 rad/s (Figure 6E).Together with the two timescales of decay (Supplementary Figure S3B), these behaviors indicate that the centroid velocity angle by itself follows 1D persistent random walk.From experimentally obtained parameters of the leading edge lifetime (2.8 s) and the angular velocity 0.1 rad/s, the 1D model (see Materials and Methods, Cell Boundary Analysis section) yielded decay time of 142 s on average which matched well with the experimental data (Figure 6F).

Fourier-based morphology space analysis
To obtain a quantitative morphometry, we chose by eye 21 representative mask images each for the 3 shapes; fan-shape, split and dumbbell (Supplementary Figure S6A) and calculated the Fourier power spectrum of the cell edge coordinates and their principal components were calculated (see Materials and Methods).We found that the first two principal components were sufficient to obtain well separated clusters that represented the shape class (Figure 7A).All cell masks analyzed were distributed within a confined domain in the PC1 fourier -PC2 fourier space (Supplementary Figure S6B).The fan-shaped data were located at a low PC1 fourier and high PC2 fourier region (Figure 7A circles).The split-shape were found in the low PC1 fourier -low PC2 fourier region (Figure 7A asterisks).The dumbbell-shape was located at high PC1 fourier and high PC2 fourier (Figure 7A triangles).To see what shape features the principal components represented, we reverse calculated an artificial form by obtaining Fourier spectrum as a product of synthetic principal component vector to the eigen vector matrix composed of the basis of Fourier spectrum (see Methods).In brief, PC1 fourier indicated the aspect ratio i.e., elongation, PC2 fourier the head width, PC3 fourier the rear steepness (Supplementary Figure
How the cell shape changed during turning can be analyzed by tracking the time sequence in the PC1 fourier -PC2 fourier space.Figure 7B shows three independent samples of re-orienting cells.Here, cells were mainly located in the negative PC1 fourier region with occasional visits to the positive PC1 fourier region.This is consistent with the above observation that cells took fan-or branched-shape (negative PC1 fourier ) in addition to rare occurrence of dumbbellshape (positive PC1 fourier ). Figure 7C shows three independent samples of the dumbbell-shape forming cells.The counterclockwise circular trajectories in the PC1 fourier -PC2 fourier space signify a transition from the fan-shape to splitting then to the dumbbell-shape.On the other hand, Figure 7D shows three independent trajectories that remained in the negative PC1 region for extended period of time.These cells at least during the time window of observation fluctuated between the fan-shape and the bifurcating fingers.There was no clear relationship between the morphometry state (PC1 fourier , PC2 fourier ) and the cell speed (Supplementary Figures S7A, B).
There was, however, negative correlation between the centroid speed and the rate of state transition d PC1 fourier { } /dt but not with d PC2 fourier { } /dt (Supplementary Figures S7C, D).As the former relation was seen in the negative direction d PC1 fourier { } /dt <0, it signifies that cells accelerate when recovering from dumbbell-shape.
Besides the rate of state transition in the principal components, there should be a direct relationship between the Fourier components C n themselves and the centroid movement.Autocorrelation analysis showed that the decay rates for C -3 , C 3 , and C 4 were 7.6, 8.9, and 10.6 s, respectively (Figure 7E) and thus matched most closely to the short decay time of VAC.As for the centroid velocity itself, according to the deformation tensor-based theory of cell movement (Ohta et al., 2016), it should be proportional to where −n + m = 1.More specifically, C nm is a complex number whose absolute value |C nm | and the angle arg (C nm ) are expected to be proportional to the speed and the velocity angle of the centroid respectively.In NIH3T3 cells, it has been shown that velocity is proportional to the elongation C -2 and triangular C 3 modes of deformation multiplied by their time derivatives; i.e. (Ebata et al., 2018).However, in N. gruberi, we found little correlation between C 23 and the centroid velocity (Supplementary Figures S7E-G).Instead, we found that it was C 34 _ C −3 C 4 − C −3 _ C 4 that correlated highly with the centroid velocity angle (Figure 7F) and x-and the y-component of the centroid velocity (Figures 7G, H).The difference between Naegleria and NIH3T3 may be attributed to the fact that Naegleria has many pseudopods that are complex in shape as analyzed below.

Deep learning-based morphology analysis
To further investigate the cell shape characteristics, we employed a convolutional neural network that was previously trained to classify cell shapes based on similarity to Dictyostelium-like, HL60-like, or fish keratocyte-like shapes (Imoto et al., 2021).While the method is not suited to track shape change over time due to discrete change in the morphometry space that is sometimes introduced by uncertainty in the cell orientation during mask alignment, it has an advantage of providing an objective morphometry that is independent of known feature basis.On average, Naegleria was classified as Dictyostelium-like (high PC1 cnn , low PC2 cnn ) (Figure 8A).This was natural as it has been shown to pick up branching shapes that are elongated overall in the migrating direction (Imoto et al., 2021).We noticed substantial variability, however, in the individual cell shapes (Figure 8B; black) that exceeded those normally observed in Dictyostelium (Figure 8B; green).Shapes that deviated in the PC1 cnn direction were mapped to dumbbell-like domain in the Fourier descriptor-based morphometry (Figure 8C; orange).Those that deviated towards low PC1 cnn were mapped to the domain that showed numerous pseudopods (Figure 8D; orange).Datapoints that fell at or near the HL60-like domain (low PC1 cnn , low PC2 cnn ) were mostly fan-like (Figure 8C; blue; Figure 8D; blue) and their occurrence per timeseries showed positive correlation with the MSD exponent (Figure 8E).This is consistent with the notion that more monopolarized the cells are, the more ballistic the cell trajectories become.

Discussion
In this report, we analyzed movements of N. gruberi cells by quantifying their speed, directionality, and shape change.The locomotive speed of N. gruberi cells was around 60 μm/min, which is similar to that reported in early literatures (King et al., 1981;Thong and Ferrante, 1986).It is substantially larger in magnitude compared to that of fibroblast ~0.4-1.0 um/min (Welf et al., 2012;Passucci et al., 2019), and even larger compared to fast migrating cells such as vegetative Dictyostelium 5 μm/min (Li et al., 2008), and neutrophils 17 μm/min (Hartman et al., 1994).Despite the large speed difference, we found common features between N. gruberi and other cell types whose random motility have previously been characterized.The exponent of MSD was approximately 1.8 meaning that the random walk is non-Fickian and nonballistic at least at surface.Stronger deceleration at higher velocity implies non-ballistic movement, where the non-memory term, i.e., fluctuating components plays a dominant role in determining the next move.Similar exponent is known in MDCK cells (Dieterich et al., 2008), A549 cancer cells (Kwon et al., 2019) hematopoietic progenitor cells (Partridge et al., 2022), and T cells (Jerison and Quake, 2020).Of particular note is that the timescale where such exponent was observed for N. gruberi was about 10-100 s which is within the order of magnitude required for a cell to move one cell-body length.This seems also to be the case for MDCK cells where the exponent of 1.8 was observed at much longer time-scale of 4-20 min with corresponding length scale 4 μm-20 μm.All in all, our data combined with the observations above from earlier literatures suggest that the time scale at which cells move in a straight line is the major determinant of cells' displacement.
The other common feature found in this study was the presence of two characteristic decay time in the VAC (Selmeczi et al., 2008).For N. gruberi, these were T 1 = 6 and T 2 = 90 s, which are in the same order of magnitude as that of Dictyostelium in the vegetative (T 1 = 5.2 and T 2 = 228 s) and the starved (T 1 = 11 and T 2 = 108 s) states (Li et al., 2008).Although equivalent measurements have not been documented for neutrophils, their cell shape changes had typical time scale of 8 s (Hartman et al., 1994) and the persistence time during chemotaxis was 103 s (Itakura et al., 2013;Haastert, 2021).From the MSD measurement, the persistent time of Dictyostelium and neutrophil-like HL60 were 151 s and 278 s, respectively (Imoto et al., 2021).Interestingly, VAC of Human keratinocyte-like cells (HeCaT) whose speed was much slower (0.18 μm/min) could also be fit with the sum of two exponential functions [T 1 = 76 s and T 2 = 860 s; (Selmeczi et al., 2005)].The characteristic time scale of around 10 s was attributed to the time scale of actin polymerization in the protruding pseudopodia (Haastert, 2021).However, the pseudopod lifetime in N. gruberi was rather long; about 15-50 s.The discrepancy may be attributed to the sister pseudopods that formed from the main pseudopods which were not analyzed in our manual tracking.In some cases, the pseudopod itself also appeared to bend in one direction.In support of this notion, the autocorrelation of C -3 and C 4 had decay time of 7-10 s which matched well with the first decay time of VAC.
On the other hand, the second decay time of VAC (90 s; Figure 3D) was close to the timescale of directional persistence i.e. "run" phase estimated from the curvature wave dynamics (142 s; Figure 6F).As for the average cell speed, we found strong correlation between the centroid velocity with the coupling of deformation modes C -3 and C 4 , instead of C -2 and C 3 .This suggests that the orientation of N. gruberi cells depends not on the primary membrane protrusions but on their sister sub-structures.A further pseudopod-level analysis at finer time-scale is required to clarify the relation between the deformation modes and the branching pseudopods.The rare cells with high persistency did not take high PC1 fourier value (Figure 7D) which was opposite of Dictyostelium (Tweedy et al., 2013).This likely stems from the fact that, in N. gruberi, the elongated form was usually dumbbell-shaped which occurred when the cells stalled and reoriented.
The splitting pseudopod may entail a mechanism similar to those found in amoebozoan and metazoan cells where dendritic actin meshworks are regulated by excitable and oscillatory dynamics (Huang et al., 2013).The presence of local inhibitor of pseudopod formation in neutrophils and Dictyostelium (Xu et al., 2017) and potential lack of such in Naegleria may underlie the difference in the number of pseudopods.Alternatively, there may be local reduction in the actin cortex that are stochastic in nature.Although protrusions observed under our culture conditions did not appear as blebs, marked flow of cytosol towards the membrane observed during extension of a protrusion suggests local pressure release.Protruding form triggered by the pressure difference at a fluid-fluid interface is known as viscous fingering.The movement speed of N. gruberi was 5 times as large as that of neutrophils and Dictysotelium, but closer to that known for fragments of Physarum which also exhibit marked cytoplasmic streaming (Rieu et al., 2015) and persistent random walk (Rodiek and Hauser, 2015).Such high velocity and potential interface instability may underlie the observed branching of pseudopods.Another unique shape feature was the dumbbell-like cell shape.According to a recent theoretical model of lamellipodia-based dynamics, a similar "two-arc shape" appeared when the protrusive force was high (Sadhu et al., 2023).The dumbbell-shape may thus be a prevalent shape feature that was heretofore overlooked due to peculiarity of the model cells.Indeed, a similar dumbbell-shape has been reported in fragmented Physarum polycephalum (Rieu et al., 2015).
In the E. coli run-and-tumble, the underlying biochemical network has been proposed to be optimally designed to extract binary information in a noisy environment (Nakamura and Kobayashi, 2021).Some bacterial species make use of multiple run modes that differ in how they are modulated in the presence of chemoattractants (Alirezaeizanjani et al., 2020) suggesting diversity and depth at which random walk strategies are likely employed in prokaryotes.In Dictyostelium amoebae, the run length increases under starvation (Haastert and Bosgraaf, 2009) which may be related to their foraging strategy.In immune cells, high correlation between cell speed and persistence is thought to underlie their search efficiency in vivo (Shaebani et al., 2020;Shaebani et al., 2022).Cancer cells show persistent random walk in the metastatic state while weakly persistent in non-metastatic state (Huda et al., 2018).Although chemoattractants for N. gruberi are so far unknown, in Naegleria fowleri, formylated peptides are known to act as chemokine (Marciano-Cabral and Cline, 1987) meaning that it enhances cell polarity and movement in the absence of gradient.Cell-cell variability in such response may explain how a minority of N. gruberi cells under our experimental condition showed persistent monopolarity.Naegleria fowleri are one of several known "brain-eating" amoebae that cause fatal central nervous system infection called amebic meningoencephalitis.Their pathogenicity is thought to be related to their capacity to enter brain by penetrating nasopharygeal mucosa and migrate along olfactory nerves (Thong and Ferrante, 1986).In future works, it should be informative to study how the properties quantified in this work are modulated by chemotactic and chemokinetic factors and how they are related to exploratory and invasive strategies.

Cell culture
Naegleria gruberi strain NEG-M was obtained from American Type Culture Collection (ATCC 30224).For routine cell propagation, small bits of frozen stock were scraped off using a sharp needle onto a fresh lawn of Klebsilella aerogenes on a NM agar plate (Peptone, Dextrose, K 2 HPO 4 , KH 2 PO4, 2% bactoagar) (Fulton, 1970).The two-member culture plate was incubated at 30 °C for a few days until cleared plaques appeared.To start axenic culture, growing cells were picked from the edge of a plaque and suspended in Milli-Q water.10 μL of the cell suspension was added to 25 mL modified HL5 media (Fulton, 1970) supplemented with 40 ng/mL vitamin B12 and 80 ng/mL folic acid, 10% fetal bovine serum (FBS, Sigma 172,012) and 1% Penicillin-Streptomycin (Gibco) in a 75 cm 2 canted-neck plastic flask (Corning 431464U).Cells were allowed to attach to the bottom of the flask and incubated at 30 °C for 3 days before harvesting for imaging.

Time-lapse imaging
Axenic growing cells were dislodged from the flask bottom by gentle agitation.Cells were pelleted by centrifugation at 7 × 10 2 G for 3 min and resuspended in fresh HL5 media.The medium contains 5 mM KH 2 PO 4 /Na 2 HPO buffer and thus provides required electrolytes (King et al., 1979) for optimal migration.The cell density was adjusted to 3.3 × 10 2 cells/mL for the observations.3 mL of the cell suspension was plated on a 35 mm glass bottom dish (No. 0 20 mm hole diameter, MatTek).The plate was set to the stage of an inverted microscope (IX81, Olympus) equipped with either a thermal plate or a closed stage-top incubator set to 30 °C and left still for 30 min before starting time-lapse image acquisition.All image acquisition was performed at 30 °C.
Phase contrast images were obtained by ×40 (LUCPLFLN) objective lens and a sCMOS camera (Prime 95B, Photometrics).To track target cells at multiple non-overlapping fields of view, Micromanager software with a custom written plugin was employed.Timelapse images were obtained from 2 or 3 positions at an interval of 1 s for up to 1 h.Each position was chosen so that initially only a single cell at the center existed in the entire field of view.In between each image acquisition, the cell centroid was calculated from a mask obtained by applying the "Make Binary" function in ImageJ to the most recent image.The automated stage was then recentered to cancel out the centroid displacement.
i 0 on a relative position on Γ: which is rescaled according to the total number of pixels A in the mask, and the coordinate was set so that the origin is at the cell centroid.
The elliptic Fourier descriptor was calculated by taking the discrete Fourier transformation of r i with wave number k: where R(•) is a rotational matrix.Its power spectrum S k | q k | 2 was calculated.C n and C -n are complex number equivalents of q n−1 and q 161−n .

FIGURE 3
FIGURE 3 Statistics of the centroid movement.(A) The trajectories used for detailed analysis (N = 35).The origin is set at the initial position.Scale bar: 500 μm.(B) The ensemble averaged MSD.Red solid curve line: exponent 1.8, which is the result of fitting the ensemble averaged MSD.Shaded region: standard deviation.(C) Probability distribution of | v|.Red solid curve: a fitting curve from a 2-dimensional isotropic Gaussian distribution with a standard deviation of 0.86 μm/s.(D) The ensemble averaged VAC.Red solid curve: the sum of two exponential functions (Eq.2).(E) Acceleration parallel (blue) and orthogonal (green) to the velocity.The binned average (circle) and the standard deviation (star).(F) Persistence of displacement orientation in a unit time step.Cosine of the angle change θ in the velocity in 1 s interval.The binned average (circles) and the standard deviation (error bars).
FIGURE 4 Statistics of the persistent random walker trajectories.(A) Simulated trajectories.(B) MSD; the ensemble average (circle) and the standard deviation (shade) (red).Experimental data (black) are duplicated from Figure 3B for comparison.(C,D) Probability distribution of | v| (C) and VAC (D) (red).Experimental data from Figures 3C, D (black) are duplicated for comparison.(E) Acceleration, parallel (red) and orthogonal (black) to the velocity.The average (cross mark) and the standard deviation (triangle).Blue and green markers show the experimental data in Figure 3E.(F) The average of cos θ (red circle) and the standard deviation (red error bar).Experimental data from Figure 3F (black circle and error bar) are duplicated for comparison.

FIGURE 5
FIGURE 5 Cell boundary analysis.(A,B) The curvature (A) and the normal velocity (B) of the cell boundary taken from a representative cell exhibiting random walk.White arrows: splitting.Black arrows: pair annihilation.(C,D) Magnified view of a subsection in (A,B).Orange arrows indicate protrusions that split (C) or annihilated (D).(E,F) The curvature (E) and the normal velocity (F) of the boundary taken from a cell with high persistence.(G) Snapshots of the cell analyzed in (E,F).The white arrow: the direction of the centroid movement.The inverted triangles mark the uroid.

FIGURE 6
FIGURE 6 Relation between the membrane protrusions and the centroid velocity angle.(A) The protruding front (white) and the retracting rear (black) detected from the velocity kymograph are overlaid on top of the curvature kymograph (see Methods).(B) The angular histogram of the protruding front (blue) and the retracting rear (orange) relative to the cell orientation as determined by the centroid velocity.(C) The position of protrusive regions ("curvature waves;" green).The region that co-extended most closely in the direction of the cell centroid motion ("leading protrusion;" magenta).The binarized mask of the protrusion region (white) obtained from the curvature kymograph is shown in the background.(D) Duration time histogram of the leading protrusion (magenta in (C)).(E) Histogram of the angular velocity in the protrusion orientation [the vector normal to the cell contour at positions indicated in green in (C)].Solid lines are exponential fit to the data (D,E).(F) Estimated VAC decay time T est for the representative data.

FIGURE 7
FIGURE 7 Fourier analysis of the cell contour.(A) Principal component space (PC1 fourier , PC2 fourier ) obtained from 63 manually selected binarized snapshots (left panel).Representative cell masks (right panels).(B) Time series in PC1 fourier -PC2 fourier space from 3 representative timelapse sequences (colors).The time spent in the negative PC1 fourier region per total trajectory time was 676 s/814 s (blue), 287 s/382 s (orange), and 459 s/578 s (green).(C) Time evolution of PC1 fourier and PC2 fourier of 10 s around a large turn that involves transition to the dumbbell shape (3 representative events; colors).Black arrow indicates the direction of time evolution.(D) Time evolution of PC1 fourier and PC2 fourier during persistent migration.Colors indicate different time series [duration: 269 s (blue), 1,039 s (orange), or 3,600 s (green)].(E) Autocorrelation of C n (n = ± 2, ± 3, ± 4).The decay rate: 12.2 (C 2 ), 18.0 (C -2 ), 8.9 (C 3 ), 7.6 (C -3 ), 10.6 (C 4 ), and 4.4 (C -4 ) seconds.(F) Distribution of angles of centroid velocity and C 34 .(G,H) The x-(G) and y-components (H) of centroid velocity plotted against real (G) and imaginary (H) parts of C 34 .Red circles and blue bars indicate average and standard deviation of centroid velocity binned with the value of C 34 .Orange lines indicate the result of fitting with linear proportionality.

FIGURE 8
FIGURE 8 Shape analysis by a CNN-based classifier.(A,B) Time-average (A) or snapshot (B) of PC1 CNN and PC2 CNN from N. gruberi images (black) were superimposed on the PC1 CNN -PC2 CNN phase space of D. discoideum (green), HL60 (red), and fish keratocyte (yellow).(C) Snapshot data in (B; black) was classified into split (black), fan (blue) or dumbbell (orange) based on PC1 fourier and PC2 fourier .(D) PC1 fourier and PC2 fourier of HL60-like (blue), cells with PC1 CNN lower than −1.6×10 6 (orange), or the other cells (black) classified with CNN.(E) Ratio of frames whose shape was classified as HL60 in deep leaning-based classification.

TABLE 1
Parameters for the generalized Langevin equation.The experimental data were fitted with the analytical VAC (Eq.29).