- 1College of Computer, Mathematical and Natural Sciences, Institute for Physical Science and Technology, University of Maryland, College Park, College Park, United States
- 2Department of Physics, College of Computer, Mathematical and Natural Sciences, University of Maryland, College Park, College Park, MD, United States
- 3Department of Chemistry and Biochemistry, College of Computer, Mathematical and Natural Sciences, University of Maryland, College Park, College Park, MD, United States
- 4National Institute of Neurological Disorders and Stroke (NIH), Bethesda, MD, United States
The development of axons and dendrites (neurites) in a neural circuit relies on the dynamic interplay of cytoskeletal components, especially actin, and the integration of diverse environmental cues. Building on prior findings that actin dynamics can serve as a primary sensor of physical guidance cues, this work investigates the role of nanotopography in modulating and guiding actin waves and neurite-tip dynamics during early neural circuit development. Although actin dynamics is well known to contribute to pathfinding in wide axonal tips, typically referred to as growth cones, we also observe dynamic actin remodeling throughout neurites and at other, narrower, neurite tips. We find that actin-wave speeds do not change significantly in the first 2 weeks of neurite development on flat substrates, but decrease over the same period in neurites on nanoridges. The ability of nanoridges to guide actin waves and the neurite-tip direction also decreases as neurites mature, both for narrow tips and wide growth cones. This change in responsiveness to physical guidance cues with neuronal maturation may impact the regenerative capacity of developing neural cells that are inserted into mature brains.
1 Introduction
The development of a functional nervous system is a marvel of Nature, relying on the precise navigation and connection of countless neurons. This intricate process of neural circuit formation relies on the dynamic behavior of developing neurons that extend neurites – projections that will ultimately become axons and dendrites (Ma et al., 2022; Yogev and Shen, 2017). This navigation depends critically on the active and tightly regulated remodeling of the actin cytoskeleton (Nakajima et al., 2024; Inagaki and Katsuno, 2017). Actin filaments (F-actin), in conjunction with microtubules (Zhao et al., 2017; Cammarata et al., 2016), underlie the extension and steering of axons, as well as the formation of synapses (Leite et al., 2021).
Actin waves in migrating eukaryotic cells have been studied and characterized extensively. Actin waves are domains of locally high actin density that translocate within the cell. These waves can, for example, be foci that move within a neurite (Winans et al., 2016; Flynn et al., 2009), or broad fronts within a cell body (Bretschneider et al., 2009; Beta et al., 2023), and are wave-like in that their motion reflects translocation of the density distribution and not of the actin monomers themselves. The motion of actin waves arises from the interaction of actin with its associated network of cytoskeletal regulators, which forms a system of biomechanochemical excitability that is deeply involved in cell migration. In the slime mold Dictyostelium discoideum, actin waves have been shown to be responsible for migration (Bretschneider et al., 2004) and for texture sensing (Jasnin et al., 2016; Jasnin et al., 2019). Other members of the signal-transduction network, particularly the actin-related protein 2/3 (Arp2/3) complex and PIP3, also form complex waves with specific localization patterns within a cell (Bretschneider et al., 2004; Jasnin et al., 2016; Jasnin et al., 2019; Gerhardt et al., 2014). Interestingly, actin dynamics tends not to be required for signal-transduction waves to occur (Banerjee et al., 2022; Banerjee et al., 2023; Arai et al., 2010; Iwamoto et al., 2025), yet both types of waves are necessary for important functions such as cell motility, and demonstrate an intricate interplay of signaling (Huang et al., 2013).
In developing neurons, actin waves are most commonly observed in association with neurite extension, where these waves precede and promote axon outgrowth. Although growth-cone-led axonal development has been studied widely, recent work has shown that cytoskeletal dynamics may also direct growth for narrow neurite tips (Marsh and Letourneau, 1984; Bentley and Toroian-Raymond, 1986). It was demonstrated recently, for example, that in the Drosophila wing disc, the tip of the pioneer neuron TSM1 is almost exclusively narrow and filopodial (Clarke et al., 2020a). In this case, genes for proteins such as Abelson tyrosine kinase regulate actin organization and, therefore, the advance of the axonal tip (Clarke et al., 2020a; Clarke et al., 2020b) through actin dynamics that occurs distant from the leading edge. Neurons also seem to possess mechanisms to regulate neurite elongation that act in regions other than the growth cones. For example, differential neurite growth, including the selection of a single axon during initial polarization and the alternating growth of axon branches, occurs even in the absence of growth-cone motility (Ruthel and Hollenbeck, 2000). This observation suggests that neurons possess mechanisms for local regulation of neurite elongation in regions other than the growth cones.
The morphology of neurons changes dramatically during development (Dotti et al., 1988). Dendrites begin as small protrusions from the cell body, and grow into a complex network of branches that ultimately determines how a neuron receives input (Miller and Jacobs, 1984). Neurons in an earlier stage of development – which we define as days in vitro (DIV) 2 through 5 in our study – are characterized by a simpler morphology, with primary and secondary dendrites and likely to have a specified axon. Neurons in the middle developmental stage, DIV 7 through 12 in our study, have more complex dendritic trees with tertiary dendrites and longer axons. Finally, neurons in a later developmental stage, DIV 14 through 23 in our study, have intricate dendritic structures, and are likely to be developing synaptic spines. Given how neurites (axons and dendrites) change significantly during development, one goal of our study is to determine whether the dynamics of the cytoskeleton also changes during development.
Further, although the actin waves that traverse the length of the axon are associated with growth and retraction, conflicting views have emerged regarding the precise roles of these waves. On the one hand, neurite elongation, which may reach lengths of hundreds of microns, can occur independently of actin waves (Mortal et al., 2017). Conversely, studies have shown that, during early neuronal development, actin waves drive fluctuations in neurite outgrowth by controlling microtubule transport (Winans et al., 2016). In some cases, competition among neurites, which is mediated by actin waves, is proposed to allow the neuron to explore its environment until external cues specify the axon (Hjorth et al., 2014). These findings suggest that actin waves act more as passive responders to stochastic internal forces rather than as active drivers of axon extension and growth-cone formation.
Classic studies of axon patterning, such as those described above, have largely focused on what happens inside the axon and growth cone. However, the physical properties of the microenvironment, such as substrate stiffness (Moore and Sheetz, 2011) and topography, also can exert a profound influence on neuronal growth behavior in different stages of development (Simitzi et al., 2017; Villard, 2023). In early stages (4–20 h post plating), specific topographic cues can direct neuronal polarization by attracting key markers, resulting in preferential axon growth along micropatterned features (Micholt et al., 2013). With further development (4–5 days post plating), growth cones on such substrates became more elongated and aligned with respect to the neurite shaft (Vecchi et al., 2024). Moreover, substrate stiffness is not only relevant early in development, but also later in the mature brain (Yin et al., 2018). Studies indicate that region-specific changes to stiffness occur with age, and particularly that the aging brain decreases in stiffness (Previtera et al., 2010a; Previtera et al., 2010b). In neurodegenerative diseases, such as Alzheimer’s (M et al., 2011; Murphy et al., 2016) or Lewy-body dementia (Pavuluri et al., 2022), brain stiffness is decreased even further in areas relevant to those conditions. These findings highlight the importance of considering the three-dimensional complexity of the in vivo environment when studying the role of actin dynamics in axonal guidance, growth-cone navigation, and longer-term neural circuit function.
The current work investigates the role of nanoridges in modulating actin waves in neurites and growth cones during neuronal development. The guidance of neurites by nanoridges has been examined previously (Marcus et al., 2017), and our goal is to characterize the behavior of actin waves under this well-studied mechanical cue (Gangatharan et al., 2018) and under control conditions. Although nanoridges are a simplified model for variable physical cues, the ridges may mimic pre-existing nerve structures along which follower axons grow, and that are mechanically distinct from the surrounding epithelium. Prior work suggests that axon-guidance mechanisms differ between follower and pioneer axons (Forghani et al., 2023), lending physiological relevance to our experimental setting. The nanoridges we used had a 3.2 μm spacing, allowing us to compare with our previous work on other cell types (Sun et al., 2015; Driscoll et al., 2014; Yang et al., 2023; Herr et al., 2022; Yang et al., 2021). We hypothesized that this spacing would enable the guidance of neurites (Fiala et al., 1999; Lamoureux et al., 2010). Although the nanoridges have a stiffness comparable to that of physiological substrates such as collagen fibers, the ridges are rigidly fixed to the underlying surface, and do not deform. This distinction highlights the fact that nanotopographic cues mimic some, but not all, of the mechanical aspects of the in vivo environment. We examine the speed and the guidance of actin waves, neurites, and growth cones by this nanotopography. Because actin waves have been identified as primary sensors of the physical microenvironment, including electric fields (Yang et al., 2022) and texture (Yang et al., 2023; Bull et al., 2025), our work also sheds light on whether neuronal development also alters the ability of cells to sense their physical environment.
2 Materials and methods
2.1 Primary neuronal dissections and cell culture
Primary neurons were obtained from Sprague–Dawley rats housed at the University of Maryland (with approval by the University of Maryland Institutional Animal Care and Use Committee; protocols R-JAN-18-05 and R-FEB-21-04). At embryonic day of gestation 18 (E18), primary cortical neurons were obtained via dissection, and were cultured as previously described (O'Neill et al., 2023). We did not distinguish between cortices from male and female embryos. Briefly, hippocampi and cortices were dissected from decapitated embryos. Cortices from all embryos were combined and triturated manually, and then were counted using a hemocytometer. A monolayer of cortical neurons was plated in 35-mm-diameter imaging dishes with glass centers (P35G-1.5-20-C, MATTEK) at a density of 5 × 105 cells per dish (520 cells/mm2). Prior to the plating, the imaging dishes were coated with 0.1 mg/mL poly-D-lysine (PDL; P0899, MilliporeSigma) for at least 1 h at 37 °C and then were washed three times with phosphate-buffered saline (PBS; SH30256.01, Cytiva). Monolayers were also plated on acrylic nanoridges, which received the following treatments prior to neuronal plating: soaking in ethanol (24 h), followed by soaking in medium (24 h), followed by PDL coating (24 h), and finally 10 μg/mL laminin coating (24 h; L2020, MilliporeSigma), all at 37 °C. Cultures were maintained in NbActiv4 medium (Nb4-500, BrainBits) supplemented with 1% penicillin/streptomycin (P/S; 15070063, Thermo Fisher Scientific) at 37 °C in 5% CO2. Half of the medium was changed twice per week. Note that we did not optimize the nanoridges for longer term neuronal survival (DIV14+), which is why our data for cells on nanoridges are only reported at the early and middle developmental stages.
2.2 Transduction with Actin-GFP
For visualization of actin dynamics within primary cortical neurons, we used CellLight Actin-GFP, BacMam 2.0 (C10506, ThermoFisher Scientific). This viral vector has low toxicity and labels cells sparsely. Given that neurons have complex morphological structures even in early stages of development, the latter feature is particularly helpful for analysis of single-cell actin dynamics. At 48–72 h prior to imaging, 5–10 μL of CellLight was mixed with the cell-culture medium (1–2 particles per cell). The medium was fully changed approximately 20 h after transduction, and imaging was performed 24–48 h later.
We chose to use actin-GFP transduction via a baculovirus. The low toxicity of the baculovirus promotes healthy primary cultures, and the low efficiency of transduction of this vector ensured that we could analyze actin dynamics of individual neurons. However, using transduction of actin-GFP rather than a fluorescent actin dye introduced certain limitations to our study. In particular, our nanotopopgrahic surfaces, although optimized for primary cell survival, are not ideal for longer-term (DIV 10+) neuronal survival. Our preliminary studies of transduced, late-stage neurons on surfaces did not yield healthy cultures. Therefore, we chose to limit our studies of neurons on ridges to the early and middle developmental stages, to ensure the integrity of, and consistency within, our results. The combined challenge of low transduction efficiency and primary culture health on non-standard surfaces is also the reason that the number of replicates for experiments on nanotopographic surfaces is lower than on control surfaces.
2.3 Confocal imaging
Live-cell imaging of primary cortical actin dynamics was performed at the University of Maryland Imaging Core using a PerkinElmer spinning-disk confocal microscope. For all experiments, we used an oil-immersion, ×40 objective (1.30 NA; 0.36 μm/pixel), and maintained live cultures under temperature, humidity, and CO2 control. We recorded 16-bit images using a Hamamatsu ImagEM X2 EM-CCD camera (C9100-23B) and acquired time-lapse images with PerkinElmer’s Volocity software (v. 6.4.0). We identified transduced neurons via a positive signal in the 488-nm excitation channel. We performed imaging in growth medium and acquired images every 2 s in the 488-nm channel for 5–10 min. We ensured a consistent z-plane with Perfect Focus (Nikon PFS).
2.4 Preparation of nanoridges
Nanoridged substrates for the neuronal experiments featured the same master pattern as in a recent study that focused on the photomodification and reshaping of nanotopographic structures (Abdelrahman et al., 2024). Briefly, a large-scale pattern with 750-nm high, 400-nm wide nanoridges at a 3.2-μm pitch (Abdelrahman et al., 2024), initially fabricated on a silicon wafer using interference lithography, was replicated using soft lithographic techniques, as described in more detail elsewhere (Sun et al., 2018), to produce a polydimethylsiloxane (PDMS) mold with negative relief of the master pattern. Molding of the patterned silicon wafer was facilitated by exposure to plasma treatment in an oxygen atmosphere and subsequent surface functionalization of the wafer with (tridecafluoro-1,1,2,2-tetrahydrooctyl) methyldichlorosilane (Gelest, CAS #: 73609-36-6) under vacuum inside a desiccator. These steps reduce the surface energy of the wafer to prevent adhesion with the PDMS during the molding process. A hard PDMS layer was spin-coated onto the pattern and Sylgard 184 (Dow) comprised the remainder of the mold, based on previous work. (Sun et al., 2018; Schmid and Michel, 2000; Kang et al., 2006). The mold was peeled from the master after several baking steps. The hard-PDMS was allowed to cure at room temperature for 2 h before baking at 60 °C for 1 h. After the Sylgard 184 elastomer base and curing agent were mixed (10:1 ratio) and poured atop the hard-PDMS layer, the composite was baked an additional 70 min at 60 °C. At this point, the mold was peeled from the master. Replicas of the original pattern were then produced using nanoimprint lithography. A photocurable resin was prepared using a 1:1 w/w mixture of ethoxylated (Cammarata et al., 2016) trimethylolpropane triacrylate (Sartomer, CAS #: 28961-43-5) and tris (2-hydroxy ethyl) isocyanurate triacrylate (Sartomer, CAS #: 40220-08-04), and 3 wt% of the photoinitiator Irgacure TPO-L (BASF, CAS #: 84434-11-7). Glass coverslips featuring pendant acrylate groups were produced by plasma treatment under oxygen for 3 min followed by stirring in a 2% solution of (3-acryloxypropyl) trimethoxysilane (Gelest, CAS #: 4369-14-6) for at least 12 h. A drop of the resin was sandwiched between the PDMS mold and an acrylated coverslip, and then was exposed to UV light for 5 min to transfer the pattern to the photoresist. The patterned coverslip was then detached from the mold.
2.5 Image analysis
2.5.1 Determination of actin-flow vectors using optical flow
In this paper, “actin activity” refers to dynamic actin-polymerization events, which are visible as transient increases in fluorescence. These intensity changes propagate in space and time, and are tracked as distinct movement trajectories to quantify single-cell actin dynamics. To prepare the videos for analysis, image-registration and intensity-adjustment functions from the scikit-image module (Python) were used to correct for stage drift (jitter) and photobleaching, respectively. To reduce random background noise, processed videos were smoothed temporally using a Gaussian-like filter applied across five successive frames (Simoncelli, 1994). The smoothing timescale was considerably shorter than the timescale of detectable actin-wave propagation, ensuring preservation of relevant dynamics. Subsequently, optical flow was performed on the smoothed image frames using the Lucas–Kanade method (Lucas and Kanade, 1981). A Gaussian weighting matrix (σ = 2 pixels = 0.72 µm) was applied around each pixel location to extract local flow vectors. A reliability threshold of 0.1 was used to remove noisy flow vectors. We verified that varying the size and shape of the Gaussian filter, as well as making modest changes in the reliability threshold, did not significantly affect motion capture or downstream analyses.
2.5.2 Neuronal masking
Scikit-image (Python) functions were used for the generation and refinement of binary masks. Mean frames were generated, and low-contrast frames were processed using contrast-limited adaptive histogram equalization. A white top-hat filter (with a disk-shaped structuring element of radius = 15 pixels = 5.4 µm) was applied to highlight darker neuronal processes. Adaptive local thresholding (Niblack, 1985) yielded an initial binary mask, which was refined by multiple iterative morphological opening steps, retaining the n largest connected components after manual review and background-noise exclusion. Subsequently, manual selection of a region corresponding to the cell body was employed due to the inherent lack of consistent alignment of the cell bodies. This approach enabled the isolation of neuronal processes from the larger cell mask, and ensured that subsequent analyses focused on the dynamic behavior of actin within the neurites. Neuronal process orientation was determined using a rotating Laplacian-of-Gaussian filter on the mean image-frame (Mennona et al., 2023).
2.5.3 Clustering of optical-flow vectors
Optical-flow vectors were spatially clustered by quantifying their local alignment with neighboring vectors. For each pixel, we computed the dot product between its flow vector and the average flow vector within an 11 × 11 Gaussian-weighted neighborhood (σ = 5 pixels ≈ 1.8 µm). This process enhanced locally coherent motion patterns by reinforcing flow directions that were aligned with their surroundings. (Lee et al., 2020). Cluster peaks were identified with the trackpy (Python) module using an approximate detection diameter of 5 pixels = 1.81 µm and a minimum separation of 15 pixels = 5.4 µm between two neighboring peaks. Detected particle positions for each frame were stored, and were subsequently linked across frames to form trajectories with a search range of 5 pixels = 1.8 µm between successive frames and a memory of one frame (2 s), enabling cluster-track analysis. To emphasize longer-scale actin flows, only tracks persisting for at least eight frames (14 s) were analyzed. We computed and stored both per-frame and per-track metrics. Per-frame metrics included the instantaneous displacement vector, the movement angle, and the instantaneous speed (magnitude of displacement per frame). Per-track scalar metrics included the total track length (the sum of frame-to-frame displacements), the duration (the number of frames), the net displacement (the straight-line distance between the starting and ending points), the average speed (the total track length divided by the duration), the net velocity (the net displacement divided by the duration, reflecting directional persistence), and the sinuosity (the ratio of the total track length to the net displacement). A summary of how these metrics vary with developmental stage and substrate type is presented in Supplementary Figure S1. In all our downstream analysis, we weighted the histograms of the actin-track metrics by the track duration, so that longer-lasting events contributed more to the distributions. To exclude static clusters, only tracks with an average instantaneous speed greater than 0.5 pixels/frame = 5.4 μm/min were retained.
2.5.4 Recurrence
To quantify how consistently new tracks emerged in similar spatial locations, we computed the mean spatial recurrence of track start positions across successive time windows. A track’s start position was defined as its earliest detected (x,y) coordinate. The full time-lapse video was divided into non-overlapping windows of a fixed frame duration and within each window, the spatial coordinates of newly appearing tracks were extracted. For each pair of adjacent windows (win1, win2), we identified recurrent tracks as those in win2 whose start positions were within a fixed spatial radius of any start position in win1. The recurrence fraction was computed as the number of recurrent tracks in win2 divided by the total number of new tracks in win2.
To ensure statistical robustness, windows with fewer than a prespecified minimum number of tracks (e.g., 5) were excluded. For each condition, the mean recurrence was calculated by averaging the recurrence fractions across all valid window pairs. The standard error of the mean (SEM) was computed as the standard deviation of these fractions divided by the square root of the number of valid pairs and is reported as error bars in figures.
2.5.5 Growth-cone orientation
Growth-cone boundaries were identified manually and were extracted using ImageJ. The alignment of the central shaft was then determined by identifying its major axis.
2.5.6 Wasserstein distance
The Wasserstein distance (W1) was used to quantify the dissimilarity between two probability distribution functions (PDFs) (Vaserstein, 1969). This metric represents the area between the corresponding cumulative distribution functions (CDFs). Mathematically, for two probability measures μ and ν on the real line with CDFs F and G, respectively, the Wasserstein distance (W1) is defined as W1(μ, ν) = ∫ |F(x) - G(x)| dx. The Wasserstein distance is a useful tool for quantifying modest changes in broad distributions.
2.5.7 Pooling of DIV videos and statistical testing
For dynamics analysis, on flat surfaces we categorized videos from DIV 2 to 5 as ‘early’ (N = 17), from DIV 7 to 12 as ‘middle’ (N = 21), and from DIV 14 to 23 as ‘late’ (N = 17) developmental stages. On nanoridges, DIV 2 to 4 were grouped as ‘early-stage’ (N = 9), and DIV 7 to 8 as ‘middle-stage’ (N = 6).
To compare multiple independent test groups, we first performed the Shapiro-Wilk test to assess normality (Shapiro and Wilk, 1965). If the data met the assumption of normality, we conducted a one-way ANOVA to test for overall group differences, followed by Tukey’s post hoc test to determine which group means differed significantly (Tukey, 1949). For data that did not meet normality assumptions, we used the Kruskal–Wallis test instead (Kruskal and Wallis, 1952), followed by Dunn’s post hoc test for pairwise comparisons (Dunn, 1964). For comparisons between two groups, we used the t-test (Student, 1908) when both groups passed the normality test (Shapiro and Wilk, 1965), and the Mann–Whitney U test when either or both groups did not pass the normality test (Mann and Whitney, 1947). Statistical test results have been presented in a Supplementary Table.
3 Results
In this work, we used primary embryonic cortical neurons transduced with actin-GFP to study how the dynamics of actin changes during development and when neurites interact with nanotopographic surfaces. In particular, we sought to examine actin dynamics on a timescale of seconds. Because it is known that the dendritic arbor of neurons changes significantly during development (Dotti et al., 1988), we observed primary neurons between DIV 2 and 23 to determine whether actin dynamics changes as significantly as does the dendritic arbor during development. For neurons cultured on PDL-coated glass (control), we grouped our findings into early (DIV 2-5), middle (DIV 7-12), and late (DIV 14-23) stages. For neurons cultured on nanotopographic surfaces with 3.2 μm ridges, we grouped our findings into early (DIV 2-4) and middle (DIV 7-8) stages.
3.1 Actin dynamics is persistent throughout different stages of development on both flat and nanoridged surfaces
To understand how actin-based motility varies during neuronal development and in response to substrate nanotopography, we first visualized actin activity under multiple conditions. Specifically, we sought to determine whether actin dynamics remains localized and transient, or shifts toward more directed, long-range behavior as neurons mature or encounter structural guidance cues. Figure 1 illustrates the cluster-tracking workflow’s ability to extract actin-wave dynamics from live-cell imaging videos for a range of developmental stages and substrate types. A representative mean-image from a neuron cultured on a flat surface (Figures 1a,b) highlights an actin-flow event, which appears as a moving bright spot (arrows). The corresponding cluster tracks extracted from this video reveal localized movements with no apparent long-range flow or strong directional bias.

Figure 1. Actin waves in primary rat cortical neurons in different DIV stages on flat and nanoridged surfaces. (a) Mean image of a representative DIV 7 neuron cultured on a flat surface. The green box indicates a region of interest, which is shown in a magnified version in (b). (b) The trajectory of an actin wave (arrows) moving along a neuronal process. Actin tracks overlaid from representative videos, color-coded by individual track, for each developmental stage on (c) flat surfaces and (d) 3.2-μm-spaced nanoridges. On flat surfaces, the imaging durations were 5 min, 4 min 10 s and 8 min, respectively, for the representative videos from early, middle and late developmental stages. On nanoridges, the imaging durations were 5 min, 5 min, 10 min and 10 min, respectively, for the representative videos from early, early, middle and middle developmental stages. The orange lines indicate the nanoridge alignment.
These spatially confined actin events are observed consistently at all developmental stages on flat substrates (Figure 1c), and in early- and mid-stage neurons on nanoridged substrates (Figure 1d), suggesting that the fundamental nature of actin activity remains transient (Supplementary Figure S2) and spatially restricted (Supplementary Figure S3) despite changes in growth conditions. Interestingly, in some cases, the cluster tracks appear to be aligned visually with underlying neuronal extensions, suggesting, qualitatively, that local cytoskeletal architecture influences the orientation of actin flow. Overall, the persistence of these short-range, dynamic actin events at different developmental stages and on different topographies (Supplementary Videos S1, 2) highlights the robust and adaptable nature of actin polymerization during neuronal growth and remodeling.
3.2 Localized actin activity is less recurrent and less frequent on nanoridged substrates
Building on our observation that actin-polymerization events are brief and spatially localized at different developmental stages and on different substrates (Figure 1), we next examined whether these events tend to recur in the same locations. Recurrence of activity in specific regions reflects spatial memory in the cytoskeleton, and suggests that certain sites undergo repeated remodeling. This spatial persistence is relevant for reinforcing cellular structures, guiding polarity, and directing neurite outgrowth. To test whether these properties are influenced by the cellular environment, we compared the recurrence and frequency of actin activity on flat and nanoridged substrates.
Figure 2 quantifies the recurrence and spatial distribution of localized actin tracks within neuronal processes. We first extracted representative optical-flow vectors from a single frame (Figure 2a), and the clustered these vectors in space and time to generate robust actin tracks (Figure 2b). We define mean recurrence as the fraction of actin tracks originating in the same spatial region in successive time windows, providing a measure of how consistently actin activity reappears in specific locations over time. Figures 2c,d explore how the mean recurrence depends on the spatial search radius and the temporal window size, respectively. We found that cells cultured on nanoridges consistently exhibited lower recurrence values compared to those on flat substrates (Figure 2e), suggesting that substrate nanotopography disrupts the spatial persistence of actin waves. This trend was statistically significant at both the early (p = 0.0022) and the middle (p = 2.7 × 10−5) developmental stages, as assessed by Mann–Whitney U tests. Developmental stage alone, in contrast, did not significantly affect recurrence within the same substrate type (flat: p = 0.70; nanoridges: p = 0.27).

Figure 2. Actin-track dynamics and recurrence depend on the microenvironment. (a) Representative image frame with optical-flow vectors for a DIV2 neuron cultured on a flat surface. (b) Actin tracks observed over a 5-min interval from the same video as in (a). The circle illustrates the 10-μm search radius used to compute mean recurrence values shown in panels (d) and (e). (c,d) Mean recurrence quantifies the fraction of new actin tracks originating within a circular region during each time window. Mean recurrence as a function of (c) radius (with a 1-min time window) and (d) time window (at a fixed search radius of 10 μm) for two representative early-stage videos on flat and nanoridged substrates. The error bars represent the SEM. (e) The mean recurrence for early (N = 17), middle (N = 21), and late (N = 17) videos, respectively, on a flat substrate and for early (N = 9) and middle (N = 6) videos, respectively, on nanoridged substrates. Search radius = 10 μm and time-window = 1 min (f) Tracks/min plotted as a function of the neuronal mask area (excluding the cell body) with linear fits (through origin) for flat (N = 55) and nanoridged substrates (N = 15). (g) Normalized track counts at different developmental stages. For (e–g), each data point corresponds to a unique video. Statistical tests were performed on pooled data. Significance is reported as p-value for 0.1 > p > 0.05, * for p ≤ 0.05, ** for p ≤ 0.01, *** for p ≤ 0.001, and **** for p ≤ 0.0001.
We next quantified the overall frequency of actin waves by measuring the number of tracks per minute and normalizing by the neuronal mask area. We observed a linear relationship between the total number of actin tracks per minute and the area of the neuronal mask for all conditions tested (Figure 2f). This relationship is expected, as larger neurons have a broader cytoskeletal landscape, increasing the likelihood of detecting dynamic actin events. The linearity suggests that actin activity scales proportionally with cell area, rather than being concentrated in a specific sub-region or regulated independently of cell area. Notably, the dependence of the frequency of actin waves on cell area was weaker in cells on nanoridges (slope = 0.005) compared to those on flat substrates (slope = 0.013). This trend is further supported by normalized track counts (Figure 2g), which were significantly lower on nanoridges at both the early (p = 6.7 × 10−5) and middle (p = 1.1 × 10−5) stages. In contrast, no significant differences (flat: p = 0.67; nanoridges: p = 0.088), were observed across developmental stages for either substrate, suggesting that the nanoridged surface exerts a stronger influence on actin dynamics than does age.
Taken together, these results demonstrate that neurons not only initiate actin waves less frequently on nanoridges than on flat surfaces, but also exhibit reduced spatial recurrence of actin activity. These findings highlight the influence of physical nanotopography on the dynamics and spatial memory of cytoskeletal remodeling.
3.3 Actin-wave speeds depend on developmental stage and substrate environment
Next, we examined whether substrate topography also affects the speed at which these actin waves propagate. Actin-wave speed provides insight into how quickly polymerization-driven movements propagate through the cell, and reflects the timescales of cytoskeletal reorganization, mechanical resistance, and traveling signaling cues. Figure 3 illustrates how actin-wave speed varies with developmental stage and substrate. A representative time-lapse image from a DIV8 neuron on nanoridges (Figure 3a) highlights a dynamic growth cone in which transient increases in fluorescence intensity mark actin activity over a 10-min period. Figure 3b shows alternating retraction and extension of this growth cone over time. This pulsatile behavior is typical of actin waves involved in guided protrusions in vivo. (Winans et al., 2016; Flynn et al., 2009; Clarke et al., 2020a;Clarke et al., 2020b).

Figure 3. Actin-wave speed as a function of developmental stage and substrate identity. (a) A representative actin time-lapse image of a DIV8 neuron on a nanoridged surface. The inset shows a highly active growth cone over a 10-min period. The orange line indicates the nanoridge alignment. (b) The growth-cone trajectory from the same time-lapse image, showing alternating retraction and forward extension along the ridge direction. We applied the AI-based model ZS-DeconvNet (Qiao et al., 2024) to denoise individual frames, exclusively for visualization purposes. The arrows show the locations of the leading edge of the growth cone. A fixed spatial reference point (asterisk) is shown in all three panels to aid visual comparison of the axon tip’s position over time. (c) CDF of the average wave speed for neurons cultured on flat surfaces for the three different developmental stages (N = 17, 21, 17 for the early, middle and late stages, respectively). (d) CDF of the average wave speed for neurons cultured on nanoridged surfaces for the early (N = 9) and middle (N = 6) developmental stages. In (c) and (d), the insets show the corresponding probability distribution function (PDF) near the peak. Number of bins = 20. The arrows in (d) point to the corresponding PDFs’ median values. The error bars represent the interquartile range of the CDF at each bin for all videos for a given condition, and the markers indicate the corresponding median values. For each pair of conditions, we calculated W1 between the median CDFs, as well as the maximum and minimum distances among all combinations of 25th and 75th percentile CDFs. These values are reported in the format: median (minimum to maximum). On flat surfaces, the W1 values were 0.53 (0.43–2.56) for early vs. middle, 0.54 (0.61–2.84) for middle vs. late, and 0.82 (0.66–2.91) for early vs. late stages. On nanoridges, the W1 value between early and middle-stage neurons was 1.59 (0.76–3.36).
To generalize these observations, we measured instantaneous actin-wave speeds under different conditions. On flat substrates (Figure 3c), the speed distributions were similar for the early, middle, and late developmental stages, with an average of approximately 21 μm/min. However, on nanoridged substrates (Figure 3d), a ∼10% decrease in wave speed was observed when going from the early to the middle developmental stage. We used the Wasserstein distance (W1) to compare the distribution of speeds between stages and conditions. We found that the distance between early and middle stages on nanoridges was greater than that observed for any developmental comparison on flat surfaces. These findings suggest that neuronal maturation slows actin-wave propagation on nanoridged surfaces, whereas wave speed remains relatively constant on flat substrates. This observation implies that the mechanical or structural constraints imposed by nanotopography interact with developmental changes in the cytoskeleton, modulating how rapidly actin-based protrusions can extend. Such modulation can both affect how neurons explore their environment and stabilize emerging processes during growth.
3.4 Cellular processes align actin waves more effectively than do nanoridges
Wave speed characterizes how rapidly actin structures move, but does not convey whether the waves propagate in a consistent direction or respond to local structural cues. To understand how actin dynamics is organized spatially, we next examined how actin tracks align with extracellular topography and with internal neuronal architecture. This analysis helps address whether external topography or internal cellular architecture plays a greater role in steering actin waves during neuronal development.
Figures 4a,c show 10 min of accumulated actin tracks for representative early- and middle-stage neurons cultured on nanoridges. Many of the tracks appear by eye to be aligned with the nanoridge direction. This visual impression is reinforced by zoomed-in regions (Figure 4b), which reveal that both the actin tracks and neuronal processes are generally aligned with the nanoridges. To assess whether all actin waves follow a clear path, Figure 4d highlights four representative middle-stage tracks that display prominent back-and-forth motion. These examples emphasize that many tracks are confined and lack a persistent direction, reinforcing the need to disentangle the roles of internal and external guidance cues. Figure 4e quantifies the alignment of actin tracks with nanoridges in the same representative early- and middle-stage videos shown in Figures 4a,c. Although most tracks follow the ridge direction, alignment is reduced in the middle-stage example.

Figure 4. The influence of nanoridges and developmental stage on actin-track alignment. (a,c) Actin-cluster tracks accumulated for 10 min for two DIV2 and DIV 8 neurons, respectively, on nanoridges. The orange line indicates the nanoridge alignment. (b) Zoomed-in view of selected regions with actin tracks from (a). (d) Four representative actin tracks from (c) with frequent back-and-forth motion shown within fixed 4 µm × 4 µm windows. The data points are color-coded by time progression within each video (purple to yellow). (e) Alignment of actin tracks with underlying nanoridges in the same early (DIV 2) and mid-stage (DIV 8) videos in panels (a) and (c). The dashed line represents a uniform distribution of track alignments for the same bin size. There are 20 bins. (f,g) Fraction of track segments aligned at various angles with respect to neuronal processes on flat and nanoridged surfaces, at the (f) early and (g) middle developmental stages. (h) Relative alignment of track segments on nanoridges with respect to neuronal processes and the nanoridge for the early and middle developmental stages. The error bars represent the interquartile range of the PDF (or CDF in insets) at each bin across all videos for each condition, and the markers indicate the corresponding median values. From (f–h), the insets show the corresponding CDF. For each pair of conditions, we calculated W1 between the median CDFs and the maximum and minimum distances among all combinations of 25th and 75th percentile CDFs. These values are reported in the format: median (minimum to maximum). For early (f) and middle (g) stage neurons, the W1 value between flat and ridge surface movies were 3.91 (1.42–7.07) and 1.39 (1.34–7.03) respectively. On nanoridges, for early and middle stage neurons, W1 value between ridge and process guidance CDFs were 6.80 (3.0–13.5) and 12.07 (6.12–20.45) respectively. For process and ridge guidance comparisons, in early vs. middle stage movies, the W1 values were 2.12 (1.16–5.96) and 7.23 (1.64–14.83), respectively.
We next explore the alignment of actin tracks with respect to the direction of the neural process in which they reside. At both the early (Figure 4f) and middle (Figure 4g) developmental stages, actin tracks remain strongly aligned with their respective processes, regardless of whether the neurons are cultured on flat or nanoridged substrates. A comparison of early and middle developmental stages on nanoridges (Figure 4h) shows a clear reduction in track alignment with the nanoridges over time, whereas alignment with processes remains relatively stable. The inset in Figure 4f further shows that track alignment with neuronal processes consistently exceeds alignment with nanoridges, suggesting that dynamic intracellular cues, such as membrane tension or cytoskeletal reorganization, may play a more dominant role in steering actin dynamics over time than does static external topography.
It is important to note that at early stages of development, the neuronal processes themselves are closely aligned with the nanoridges. This initial bias could give the impression that nanoridges directly guide actin dynamics, when in fact the processes may serve as the intermediate conduit. As development progresses, and the processes start diverging from the ridge direction, the decreasing alignment between actin tracks and nanoridges becomes more apparent. This observation raises a fundamental question: Is actin guidance influenced by the internal architecture of the cell versus the external physical environment, or does the physical environment guide the intracellular cues, and thus the actin dynamics of growing neurite tips, in a self-evolving feedback loop?
3.5 Impact of nanoridges on growth-cone alignment
To explore the contributions of internal vs. external physical cues to cytoskeletal guidance, we next consider the structural features of developing growth cones. Specifically, we analyzed how the growth-cone alignment is related to nanoridge orientation over the course of development. Figures 5a,b show representative examples of growth cones that are aligned either parallel (at DIV 2) or nearly perpendicular (at DIV 8) to the nanoridges. The growth cone in Figure 5b, marked by the green box, is especially notable for originating from a neighboring neuron located outside the field of view, and extending along a neurite that crosses multiple nanoridges. Thus, by later stages of development, neurites are capable of growing across the ridges, potentially enabling more diverse growth-cone orientations. In Figure 5c, we show the growth-cone alignment as a function of growth-cone width, revealing a developmental shift. At early time points, growth cones display a strong preference for alignment along the nanoridges. This preference decreases as development progresses.

Figure 5. Growth-cone orientations at different developmental stages for neurons on nanoridges. Representative growth cones oriented (a) parallel to (DIV 2) and (b) perpendicular to (DIV 8) the nanoridges. The nanoridge directions are indicated by orange lines, and the growth-cone orientations by arrows. (c) The distribution of growth-cone alignment angles with respect to the nanoridge as a function of the growth-cone width. Darker to brighter colors indicate neuronal development stages from DIV 2 to DIV 8. (d) The average number of growth cones per neuron as a function of the alignment angle with respect to the nanoridges for early- and middle-stage neurons. Growth cones at angles between 0° and 15° are considered to be parallel, whereas growth cones at angles between 75° and 90° are considered to be perpendicular. (e) The relationship between the cross-sectional area and width of the growth cones, color-coded for developmental stage. Videos with at least one growth cone were considered for analysis. N = 7, 4, 1, 4, 3 for DIV 2, 3, 4, 7, and 8, respectively. Due to low N, we did not run statistical tests on the average number of growth cones across DIV stages and alignment angles.
We also observe that growth cones rarely adopt orientations around 45°, instead favoring directions that are either approximately parallel to or nearly perpendicular to the nanoridges (Figure 5d). Although the average number of growth cones aligned parallel to the ridges remains relatively consistent over time, there is a substantial increase in the number of perpendicularly aligned growth cones during later stages.
This bimodal distribution suggests a shift in the balance of guidance mechanisms over time. Early in development, growth cones appear to align closely with the nanoridges, indicating that the nanoridges serve as a primary external cue. As neurons mature and more neurites populate the surface, spanning across multiple ridges, the physical landscape becomes increasingly complex. In this setting, new growth cones may interact not just with the substrate, but also with these pre-existing neurites, which create internal mechanical and geometric constraints.
Interestingly, we also observed perpendicular growth in areas lacking visible pre-existing neurites, implying that additional factors, such as developmental changes in adhesion or cytoskeletal dynamic, may be at play. Similar dual-mode guidance patterns have been reported previously (Herr et al., 2022; Ziebert and Aranson, 2016; Chen et al., 2019). These shifts likely lead to a transition from continuous to more discrete preferred growth orientations, contributing to the observed bimodality. One possible mechanism is that new growth cones follow earlier neurites that have already crossed the ridges. These existing neurites may exert mechanical tension while straightening during development, promoting a greater amount of perpendicular orientation. Alternatively, geometric constraints, such as filopodia reaching the base of adjacent ridges, may pull the growth cone downward, reinforcing perpendicular growth.
Finally, Figure 5e shows that although growth-cone cross-sectional area is positively correlated with width, neither the area nor the width changes significantly with developmental stage. This observation indicates that size alone does not determine alignment or guidance behavior. Instead, it is the evolving balance between external substrate features and internal structural constraints that shapes how growth cones navigate their surroundings over time.
4 Discussion
We have characterized actin dynamics in primary rat cortical neurons at developmental stages ranging from DIV 2 (early) to DIV 23 (late), and we have quantified these dynamics using an optical-flow based technique. Because primary cells have a finite lifetime in vitro, we focused on time points at which our primary cultures were healthy (assessed twice weekly during media changes) and only used for experiments at timepoints when successful transduction was possible. In particular, because our nanotopographic surfaces were not optimized for longer-term (DIV 10+) neuron survival, we only studied neurons on ridges during the early and middle developmental stages.
An additional consideration for the interpretation of our findings is the presence of glial cells in primary neuronal cultures. Recent studies have shown that glial cells can affect neurite guidance (Learte and Hidalgo, 2007), and likely also actin dynamics. We did not control for the chemical or growth factors released by either neurons or glial cells, both of which can affect neurite guidance (Short et al., 2021). Moreover, the presence of glial cells modulates the physical environment (Vernadakis et al., 1979), which may affect neurite guidance and actin dynamics. Future work will dissect the complex interplay between the mechanical and chemical cues that drive actin dynamics and the consequences of eliminating glial cells from the in vitro environment.
Because actin dynamics can serve as a primary sensor of physical signals (Yang et al., 2023; Bull et al., 2025; O'Neill et al., 2023), we assessed how a well-established physical guidance cue, nanoridges, affects actin dynamics and neurite tips. Most previous studies (Winans et al., 2016; Clarke et al., 2020a; Clarke et al., 2020b; San Miguel-Ruiz and Letourneau, 2014; Marsick et al., 2010; Suter and Forscher, 2000) used imaging frame rates on timescales of minutes, which is too slow to capture the actin-polymerization waves that act as sensors of the physical environment. By imaging the dynamics of actin at a rate of 0.5 frames/s, we were able to analyze how polymerization and depolymerization waves in primary neurons contribute to growth-cone development, neuronal maturation, and sensing of the physical microenvironment.
We demonstrated that, on flat surfaces, the actin-wave speeds remain constant throughout neuronal development. In contrast, on nanoridges, actin waves are suppressed, and the wave speed decreases slightly with increasing neuronal maturity. These findings suggest that actin is indeed functioning as a sensor of the environment, as we found in our previous work on astrocytes (O'Neill et al., 2023): when the environment is not as physiologically relevant (flat, PDL-coated glass), actin-wave speed does not change as neurons mature, but when the environment is more realistic (nanoridged substrates), the actin-wave speed decreases with neuronal maturity. It is possible that the differential effects on speed observed on flat vs. nanoridged substrates are related to the more physiologically-relevant environment of the nanoridged substrates. Moreover, neurite orientation appears to be the primary guidance cue for the direction of actin-wave propagation throughout development. As neurons mature, the guidance of actin by nanoridges is reduced, and both neurite-tip and growth-cone guidance by nanoridges diminish, with an increase in motion perpendicular to the nanoridges for actin waves, neurite tips, and growth cones.
It would be interesting in future work to explore, via pharmacological experiments, whether branched or linear actin contributes more to the guidance we observe. Previous studies have generally indicated that branched actin is responsible for swellings in dendrites, whereas linear actin is responsible for filopodial extension (Shi et al., 2021). A delicate balance of both branched and linear actin is necessary for axonal patterning (Kannan et al., 2017; Omotade et al., 2017). Indeed, a recent preprint studying contact guidance in epithelial cells suggests that formins (regulators of linear actin) increase contact guidance, whereas Arp2/3 (a regulator of branched actin) decreases contact guidance (Foroozandehfar et al., 2025). In our system, it is possible that in younger neurons there is more formin activity (and therefore more guidance by nanoridges) and that in older neurons there is more Arp2/3 activity (and therefore less guidance by nanoridges). Possibilities for future work include dissecting the roles of branched vs. linear in nanoridge guidance and exploring the spatiotemporal relationships between actin dynamics and dynamics of the signal-transduction network in primary neurons (Jasnin et al., 2016; Gerhardt et al., 2014; Banerjee et al., 2022; Arai et al., 2010; Iwamoto et al., 2025; Huang et al., 2013).
Our findings demonstrate that the ability of neurons to sense and follow physical guidance cues may change with developmental stage, and raise the question of how such a change in physical sensing capabilities could impact brain development and regeneration. Recent studies have shown that promoting actin dynamics as a sensor of the physical microenvironment can enhance axon regeneration in the adult nervous system, particularly through the regulation of specific actin-binding proteins such as ADF/Cofilin (Leite et al., 2021). Understanding how nanotopography influences actin dynamics during various stages of development may thus provide valuable insights into potential strategies for promoting axon regeneration after injury or disease. It is important to note that our conclusions are limited to those that can be derived from an embryonic system, because even the “late stage” neurons in our study are still young compared to the adult neurons in a living organism.
Finally, our findings may have implications for the development of novel therapeutic approaches for processes that are known to involve actin dynamics, in particular neuroregeneration (Blanquie and Bradke, 2018) and the treatment of neurodegenerative diseases (Bernstein et al., 2011). Many studies have shown that brain stiffness decreases with age (Blanquie and Bradke, 2018; Bernstein et al., 2011; Arani et al., 2015; Takamura et al., 2020), and that varying substrate stiffness affects neuronal morphology (Previtera et al., 2010a; Previtera et al., 2010b). Thus, our work provides a baseline for future studies investigating how to tune actin dynamics and thus neurite tip guidance to physical microenvironments that more closely mimic the brain at various developmental or disease stages.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was approved by The University of Maryland Institutional Animal Care and Use Committee. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
SP: Conceptualization, Data curation, Formal Analysis, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review and editing. KO: Conceptualization, Data curation, Investigation, Methodology, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing. ER: Data curation, Investigation, Writing – review and editing. MH: Methodology, Resources, Writing – review and editing. CH: Methodology, Software, Writing – review and editing. JF: Conceptualization, Funding acquisition, Supervision, Writing – original draft, Writing – review and editing. EG: Conceptualization, Writing – review and editing. WL: Conceptualization, Funding acquisition, Supervision, Writing – original draft, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the US Air Force Office of Scientific Research MURI grant FA9550-16-1-0052 and grant FA9550-25-1-0002 to WL and JF. It was also supported in part by funds from the Intramural Research Program of NINDS, NIH: ZIA NS003013 to, EG.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2025.1631520/full#supplementary-material
References
Abdelrahman, M. H., Shen, J., Fisher, N. C., Losert, W., and Fourkas, J. T. (2024). Assessing the stability of azopolymer nanotopography during live-cell fluorescence imaging. Front. Bioeng. Biotechnol. 12, 1409735. doi:10.3389/fbioe.2024.1409735
Arai, Y., Shibata, T., Matsuoka, S., Sato, M. J., Yanagida, T., and Ueda, M. (2010). Self-organization of the phosphatidylinositol lipids signaling system for random cell migration. Proc. Natl. Acad. Sci. U. S. A. 107 (27), 12399–12404. doi:10.1073/pnas.0908278107
Arani, A., Murphy, M. C., Glaser, K. J., Manduca, A., Lake, D. S., Kruse, S. A., et al. (2015). Measuring the effects of aging and sex on regional brain stiffness with MR elastography in healthy older adults. NeuroImage 111, 59–64. doi:10.1016/j.neuroimage.2015.02.016
Banerjee, T., Biswas, D., Pal, D. S., Miao, Y., Iglesias, P. A., and Devreotes, P. N. (2022). Spatiotemporal dynamics of membrane surface charge regulates cell polarity and migration. Nat. Cell Biol. 24 (10), 1499–1515. doi:10.1038/s41556-022-00997-7
Banerjee, T., Matsuoka, S., Biswas, D., Miao, Y., Pal, D. S., Kamimura, Y., et al. (2023). A dynamic partitioning mechanism polarizes membrane protein distribution. Nat. Commun. 14 (1), 7909. doi:10.1038/s41467-023-43615-2
Bentley, D., and Toroian-Raymond, A. (1986). Disoriented pathfinding by pioneer neurone growth cones deprived of filopodia by cytochalasin treatment. Nature 323 (6090), 712–715. doi:10.1038/323712a0
Bernstein, B. W., Maloney, M. T., and Bamburg, J. R. (2011). Actin and diseases of the nervous system. Adv. Neurobiol. 5, 201–234. doi:10.1007/978-1-4419-7368-9_11
Beta, C., Edelstein-Keshet, L., Gov, N., and Yochelis, A. (2023). From actin waves to mechanism and back: how theory aids biological understanding. Elife 12, e87181. doi:10.7554/eLife.87181
Blanquie, O., and Bradke, F. (2018). Cytoskeleton dynamics in axon regeneration. Curr. Opin. Neurobiol. 51, 60–69. doi:10.1016/j.conb.2018.02.024
Bretschneider, T., Diez, S., Anderson, K., Heuser, J., Clarke, M., Müller-Taubenberger, A., et al. (2004). Dynamic actin patterns and Arp2/3 assembly at the substrate-attached surface of motile cells. Curr. Biol. 14 (1), 1–10. doi:10.1016/j.cub.2003.12.005
Bretschneider, T., Anderson, K., Ecke, M., Müller-Taubenberger, A., Schroth-Diez, B., Ishikawa-Ankerhold, H. C., et al. (2009). The three-dimensional dynamics of actin waves, a model of cytoskeletal self-organization. Biophys. J. 96 (7), 2888–2900. doi:10.1016/j.bpj.2008.12.3942
Bull, A. L., Mosher, M., Rodriguez, P., Fox, S., Hourwitz, M. J., Fourkas, J. T., et al. (2025). Suppressing collective cell motion with bidirectional guidance cues. Phys. Rev. E 111 (2), 024409. doi:10.1103/PhysRevE.111.024409
Chen, S., Hourwitz, M. J., Campanello, L., Fourkas, J. T., Losert, W., and Parent, C. A. (2019). Actin cytoskeleton and focal adhesions regulate the biased migration of breast cancer cells on nanoscale asymmetric sawteeth. ACS Nano 13 (2), 1454–1468. doi:10.1021/acsnano.8b07140
Clarke, A., McQueen, P. G., Fang, H. Y., Kannan, R., Wang, V., McCreedy, E., et al. (2020a). Dynamic morphogenesis of a pioneer axon in Drosophila and its regulation by Abl tyrosine kinase. Mol. Biol. Cell 31 (6), 452–465. doi:10.1091/mbc.E19-10-0563
Clarke, A., McQueen, P. G., Fang, H. Y., Kannan, R., Wang, V., McCreedy, E., et al. (2020b). Abl signaling directs growth of a pioneer axon in Drosophila by shaping the intrinsic fluctuations of actin. Mol. Biol. Cell 31 (6), 466–477. doi:10.1091/mbc.E19-10-0564
Dotti, C. G., Sullivan, C. A., and Banker, G. A. (1988). The establishment of polarity by hippocampal neurons in culture. J. Neurosci. 8 (4), 1454–1468. doi:10.1523/JNEUROSCI.08-04-01454.1988
Driscoll, M. K., Sun, X., Guven, C., Fourkas, J. T., and Losert, W. (2014). Cellular contact guidance through dynamic sensing of nanotopography. ACS Nano 8 (4), 3546–3555. doi:10.1021/nn406637c
Dunn, O. J. (1964). Multiple comparisons using rank sums. Technometrics 6 (3), 241–252. doi:10.2307/1266041
Fiala, J. C., and Harris, K. M. (1999). Dendrite structure, in Dendrite. Editors G. Stuart, N. Spruston, and M. Häusser (Oxford; New York: Oxford University Press), 1–34.
Flynn, K. C., Pak, C. W., Shaw, A. E., Bradke, F., and Bamburg, J. R. (2009). Growth cone-like waves transport actin and promote axonogenesis and neurite branching. Dev. Neurobiol. 69 (12), 761–779. doi:10.1002/dneu.20734
Forghani, R., Chandrasekaran, A., Papoian, G., and Giniger, E. (2023). A new view of axon growth and guidance grounded in the stochastic dynamics of actin networks. Open Biol. 13 (6), 220359. doi:10.1098/rsob.220359
Foroozandehfar, A., Tohidi, S., Siddiqui, S. A., Namanda, F. R., Forrester, M., Cochran, E., et al. (2025). Formins and arp2/3 reciprocally regulate contact guidance on aligned collagen fibrils. bioRxiv. 2025 (01.03), 631234. doi:10.1101/2025.01.03.631234
Gangatharan, G., Schneider-Maunoury, S., and Breau, M. A. (2018). Role of mechanical cues in shaping neuronal morphology and connectivity. Biol. Cell 110 (6), 125–136. doi:10.1111/boc.201800003
Gerhardt, M., Ecke, M., Walz, M., Stengl, A., Beta, C., and Gerisch, G. (2014). Actin and PIP3 waves in giant cells reveal the inherent length scale of an excited state. J. Cell Sci. 127 (Pt 20), 4507–4517. doi:10.1242/jcs.156000
Cammarata, G. M., Bearce, E. A., and Lowery, L. A. (2016). Cytoskeletal social networking in the growth cone: how +TIPs mediate microtubule-actin cross-linking to drive axon outgrowth and guidance. Cytoskeleton 73 (9), 461–476. doi:10.1002/cm.21272
Herr, C., Winkler, B., Ziebert, F., Aranson, I. S., Fourkas, J. T., and Losert, W. (2022). Spontaneous polarization and cell guidance on asymmetric nanotopography. Commun. Phys. 5 (1), 114. doi:10.1038/s42005-022-00889-0
Hjorth, J. J., van Pelt, J., Mansvelder, H. D., and van Ooyen, A. (2014). Competitive dynamics during resource-driven neurite outgrowth. PLoS One 9 (2), e86741. doi:10.1371/journal.pone.0086741
Huang, C. H., Tang, M., Shi, C., Iglesias, P. A., and Devreotes, P. N. (2013). An excitable signal integrator couples to an idling cytoskeletal oscillator to drive cell migration. Nat. Cell Biol. 15 (11), 1307–1316. doi:10.1038/ncb2859
Inagaki, N., and Katsuno, H. (2017). Actin waves: origin of cell polarization and migration? Trends Cell Biol. 27 (7), 515–526. doi:10.1016/j.tcb.2017.02.003
Iwamoto, K., Matsuoka, S., and Ueda, M. (2025). Excitable Ras dynamics-based screens reveal RasGEFX is required for macropinocytosis and random cell migration. Nat. Commun. 16 (1), 117. doi:10.1038/s41467-024-55389-2
Jasnin, M., Ecke, M., Baumeister, W., and Gerisch, G. (2016). Actin organization in cells responding to a perforated surface, revealed by live imaging and cryo-electron tomography. Structure 24 (7), 1031–1043. doi:10.1016/j.str.2016.05.004
Jasnin, M., Beck, F., Ecke, M., Fukuda, Y., Martinez-Sanchez, A., Baumeister, W., et al. (2019). The architecture of traveling actin waves revealed by cryo-electron tomography. Structure 27 (8), 1211–1223. doi:10.1016/j.str.2019.05.009
Kang, H., Lee, J., Park, J., and Lee, H. H. (2006). An improved method of preparing composite poly(dimethylsiloxane) moulds. Nanotechnology 17 (1), 197–200. doi:10.1088/0957-4484/17/1/032
Kannan, R., Song, J. K., Karpova, T., Clarke, A., Shivalkar, M., Wang, B., et al. (2017). The Abl pathway bifurcates to balance Enabled and Rac signaling in axon patterning in Drosophila. Development 144 (3), 487–498. doi:10.1242/dev.143776
Kruskal, W. H., and Wallis, W. A. (1952). Use of ranks in one-criterion variance analysis. J. Am. Stat. Assoc. 47 (260), 583–621. doi:10.2307/2280779
Lamoureux, P., Heidemann, S. R., Martzke, N. R., and Miller, K. E. (2010). Growth and elongation within and along the axon. Dev. Neurobiol. 70 (3), 135–149. doi:10.1002/dneu.20764
Learte, A. R., and Hidalgo, A. (2007). The role of glial cells in axon guidance, fasciculation and targeting. Adv. Exp. Med. Biol. 621, 156–166. doi:10.1007/978-0-387-76715-4_12
Lee, R. M., Campanello, L., Hourwitz, M. J., Alvarez, P., Omidvar, A., Fourkas, J. T., et al. (2020). Quantifying topography-guided actin dynamics across scales using optical flow. Mol. Biol. Cell 31 (16), 1753–1764. doi:10.1091/mbc.E19-11-0614
Leite, S. C., Pinto-Costa, R., and Sousa, M. M. (2021). Actin dynamics in the growth cone: a key player in axon regeneration. Curr. Opin. Neurobiol. 69, 11–18. doi:10.1016/j.conb.2020.11.015
Lucas, B. D., and Kanade, T. (1981). “An iterative image registration technique with an application to stereo vision,” in Proceedings of the 7th International Joint Conference on Artificial Intelligence, 674–679.
Murphy, M. C., Huston, J., Jack, C. R., Glaser, K. J., Manduca, A., Felmlee, J. P., et al. (2011). Decreased brain stiffness in Alzheimer's disease determined by magnetic resonance elastography. J. Magn. Reson Imaging 34 (3), 494–498. doi:10.1002/jmri.22707
Ma, L. (2022). “Generation of axons and dendrites,” in Neuroscience in the 21st Century: from basic to Clinical. Editors D. W. Pfaff, N. D. Volkow, and J. L. Rubenstein (Cham: Springer International Publishing), 341–356.
Mann, H. B., and Whitney, D. R. (1947). On a test of whether one of two random variables is stochastically larger than the other. Ann. Math. Statistics 18, 50–60. doi:10.1214/aoms/1177730491
Marcus, M., Baranes, K., Park, M., Choi, I. S., Kang, K., and Shefi, O. (2017). Interactions of neurons with physical environments. Adv. Healthc. Mater. 6 (15), 1700267. doi:10.1002/adhm.201700267
Marsh, L., and Letourneau, P. C. (1984). Growth of neurites without filopodial or lamellipodial activity in the presence of cytochalasin B. J. Cell Biol. 99 (6), 2041–2047. doi:10.1083/jcb.99.6.2041
Marsick, B. M., Flynn, K. C., Santiago-Medina, M., Bamburg, J. R., and Letourneau, P. C. (2010). Activation of ADF/cofilin mediates attractive growth cone turning toward nerve growth factor and netrin-1. Dev. Neurobiol. 70 (8), 565–588. doi:10.1002/dneu.20800
Mennona, N. J., Sedelnikova, A., Echchgadda, I., and Losert, W. (2023). Filament displacement image analytics tool for use in investigating dynamics of dense microtubule networks. Phys. Rev. E 108 (3), 034411. doi:10.1103/PhysRevE.108.034411
Micholt, L., Gärtner, A., Prodanov, D., Braeken, D., Dotti, C. G., and Bartic, C. (2013). Substrate topography determines neuronal polarization and growth in vitro. PLoS One 8 (6), e66170. doi:10.1371/journal.pone.0066170
Miller, J. P., and Jacobs, G. A. (1984). Relationships between neuronal structure and function. J. Exp. Biol. 112, 129–145. doi:10.1242/jeb.112.1.129
Moore, S. W., and Sheetz, M. P. (2011). Biophysics of substrate interaction: influence on neural motility, differentiation, and repair. Dev. Neurobiol. 71 (11), 1090–1101. doi:10.1002/dneu.20947
Mortal, S., Iseppon, F., Perissinotto, A., D'Este, E., Cojoc, D., Napolitano, L. M. R., et al. (2017). Actin waves do not Boost neurite outgrowth in the early stages of neuron maturation. Front. Cell Neurosci. 11, 402. doi:10.3389/fncel.2017.00402
Murphy, M. C., Jones, D. T., Jack, C. R., Glaser, K. J., Senjem, M. L., Manduca, A., et al. (2016). Regional brain stiffness changes across the Alzheimer's disease spectrum. Neuroimage Clin. 10, 283–290. doi:10.1016/j.nicl.2015.12.007
Nakajima, C., Sawada, M., Umeda, E., Takagi, Y., Nakashima, N., Kuboyama, K., et al. (2024). Identification of the growth cone as a probe and driver of neuronal migration in the injured brain. Nat. Commun. 15 (1), 1877. doi:10.1038/s41467-024-45825-8
Niblack, W. (1985). An introduction to digital image processing. Birkeroed, Denmark: Strandberg Publishing Company.
O'Neill, K. M., Saracino, E., Barile, B., Mennona, N. J., Mola, M. G., Pathak, S., et al. (2023). Decoding natural astrocyte Rhythms: dynamic actin waves result from environmental sensing by primary Rodent astrocytes. Adv. Biol. (Weinh) 7 (6), e2200269. doi:10.1002/adbi.202200269
Omotade, O. F., Pollitt, S. L., and Zheng, J. Q. (2017). Actin-based growth cone motility and guidance. Mol. Cell Neurosci. 84, 4–10. doi:10.1016/j.mcn.2017.03.001
Pavuluri, K., Huston, J., Ehman, R. L., Manduca, A., Jack, C. R., Savica, R., et al. (2022). Regional brain stiffness analysis of dementia with Lewy Bodies. J. Magn. Reson Imaging 55 (6), 1907–1909. doi:10.1002/jmri.27970
Previtera, M. L., Langhammer, C. G., and Firestein, B. L. (2010a). Effects of substrate stiffness and cell density on primary hippocampal cultures. J. Biosci. Bioeng. 110 (4), 459–470. doi:10.1016/j.jbiosc.2010.04.004
Previtera, M. L., Langhammer, C. G., Langrana, N. A., and Firestein, B. L. (2010b). Regulation of dendrite arborization by substrate stiffness is mediated by glutamate receptors. Ann. Biomed. Eng. 38 (12), 3733–3743. doi:10.1007/s10439-010-0112-5
Qiao, C., Zeng, Y., Meng, Q., Chen, X., Chen, H., Jiang, T., et al. (2024). Zero-shot learning enables instant denoising and super-resolution in optical fluorescence microscopy. Nat. Commun. 15 (1), 4180. doi:10.1038/s41467-024-48575-9
Ruthel, G., and Hollenbeck, P. J. (2000). Growth cones are not required for initial establishment of polarity or differential axon branch growth in cultured hippocampal neurons. J. Neurosci. 20 (6), 2266–2274. doi:10.1523/JNEUROSCI.20-06-02266.2000
San Miguel-Ruiz, J. E., and Letourneau, P. C. (2014). The role of Arp2/3 in growth cone actin dynamics and guidance is substrate dependent. J. Neurosci. 34 (17), 5895–5908. doi:10.1523/JNEUROSCI.0672-14.2014
Schmid, H., and Michel, B. (2000). Siloxane Polymers for high-resolution, high-Accuracy soft lithography. Macromolecules 33 (8), 3042–3049. doi:10.1021/ma982034l
Shapiro, S. S., and Wilk, M. B. (1965). An analysis of variance test for normality (Complete Samples). Biometrika 52 (3/4), 591–611. doi:10.2307/2333709
Shi, R., Kramer, D. A., Chen, B., and Shen, K. (2021). A two-step actin polymerization mechanism drives dendrite branching. Neural Dev. 16 (1), 3. doi:10.1186/s13064-021-00154-0
Short, C. A., Onesto, M. M., Rempel, S. K., Catlett, T. S., and Gomez, T. M. (2021). Familiar growth factors have diverse roles in neural network assembly. Curr. Opin. Neurobiol. 66, 233–239. doi:10.1016/j.conb.2020.12.016
Simitzi, C., Ranella, A., and Stratakis, E. (2017). Controlling the morphology and outgrowth of nerve and neuroglial cells: the effect of surface topography. Acta Biomater. 51, 21–52. doi:10.1016/j.actbio.2017.01.023
Simoncelli, E. (1994). “Design of multi-dimensional derivative filters,” in Proceedings of 1st IEEE International Conference on Image Processing, Austin, Texas (IEEE), 790–793.
Sun, X., Driscoll, M. K., Guven, C., Das, S., Parent, C. A., Fourkas, J. T., et al. (2015). Asymmetric nanotopography biases cytoskeletal dynamics and promotes unidirectional cell guidance. Proc. Natl. Acad. Sci. 112 (41), 12557–12562. doi:10.1073/pnas.1502970112
Sun, X., Hourwitz, M. J., Baker, E. M., Schmidt, B. U. S., Losert, W., and Fourkas, J. T. (2018). Replication of biocompatible, nanotopographic surfaces. Sci. Rep. 8 (1), 564. doi:10.1038/s41598-017-19008-z
Suter, D. M., and Forscher, P. (2000). Substrate–cytoskeletal coupling as a mechanism for the regulation of growth cone motility and guidance. J. Neurobiol. 44 (2), 97–113. doi:10.1002/1097-4695(200008)44:2<97::aid-neu2>3.3.co;2-l
Takamura, T., Motosugi, U., Sasaki, Y., Kakegawa, T., Sato, K., Glaser, K. J., et al. (2020). Influence of age on Global and regional brain stiffness in young and middle-Aged adults. J. Magnetic Reson. Imaging 51 (3), 727–733. doi:10.1002/jmri.26881
Tukey, J. W. (1949). Comparing individual means in the analysis of variance. Biometrics 5 (2), 99–114. doi:10.2307/3001913
Vaserstein, L. N. (1969). Markov processes over denumerable products of spaces, describing large systems of automata. Probl. Peredachi Inf. 5 (3), 64–72.
Vecchi, J. T., Rhomberg, M., Guymon, C. A., and Hansen, M. R. (2024). The geometry of photopolymerized topography influences neurite pathfinding by directing growth cone morphology and migration. J. Neural Eng. 21 (2), 026027. doi:10.1088/1741-2552/ad38dc
Vernadakis, A., Nidess, R., Culver, B., and Bragg, A. E. (1979). Glial cells: modulators of neuronal environment. Mech. Ageing Dev. 9 (5), 553–566. doi:10.1016/0047-6374(79)90094-0
Villard, C. (2023). Spatial confinement: a spur for axonal growth. Semin. Cell Dev. Biol. 140, 54–62. doi:10.1016/j.semcdb.2022.07.006
Winans, A. M., Collins, S. R., and Meyer, T. (2016). Waves of actin and microtubule polymerization drive microtubule-based transport and neurite growth before single axon formation. Elife 5, e12387. doi:10.7554/eLife.12387
Yang, Q., Hourwitz, M., Campanello, L., Devreotes, P. N., Fourkas, J., and Losert, W. (2021). Physical control of intracellular waves with nanotopography and electric fields. Biophysical J. 120 (3), 113a. doi:10.1016/j.bpj.2020.11.906
Yang, Q., Miao, Y., Campanello, L. J., Hourwitz, M. J., Abubaker-Sharif, B., Bull, A. L., et al. (2022). Cortical waves mediate the cellular response to electric fields. eLife 11, e73198. doi:10.7554/eLife.73198
Yang, Q., Miao, Y., Banerjee, P., Hourwitz, M. J., Hu, M., Qing, Q., et al. (2023). Nanotopography modulates intracellular excitable systems through cytoskeleton actuation. Proc. Natl. Acad. Sci. 120 (19), e2218906120. doi:10.1073/pnas.2218906120
Yin, Z., Romano, A. J., Manduca, A., Ehman, R. L., and Huston, J. (2018). Stiffness and beyond: what MR elastography can tell us about brain structure and function under physiologic and pathologic conditions. Top. Magn. Reson Imaging 27 (5), 305–318. doi:10.1097/RMR.0000000000000178
Yogev, S., and Shen, K. (2017). Establishing neuronal polarity with environmental and intrinsic mechanisms. Neuron 96 (3), 638–650. doi:10.1016/j.neuron.2017.10.021
Zhao, B., Meka, D. P., Scharrenberg, R., König, T., Schwanke, B., Kobler, O., et al. (2017). Microtubules modulate F-actin dynamics during neuronal polarization. Sci. Rep. 7 (1), 9583. doi:10.1038/s41598-017-09832-8
Keywords: actin dynamics, growth cone navigation, neurite development, nanotopography, cytoskeletal remodeling, physical guidance cues, neural regeneration
Citation: Pathak S, O’Neill KM, Robinson EK, Hourwitz MJ, Herr C, Fourkas JT, Giniger E and Losert W (2025) Nanotopographic control of actin waves and growth cone navigation in developing neurons. Front. Cell Dev. Biol. 13:1631520. doi: 10.3389/fcell.2025.1631520
Received: 19 May 2025; Accepted: 13 August 2025;
Published: 10 September 2025.
Edited by:
Dhiman Sankar Pal, Johns Hopkins University, United StatesReviewed by:
Michael W. Gramlich, Auburn University, United StatesJennifer Landino, Dartmouth College, United States
Tatsat Banerjee, Johns Hopkins University, United States
Copyright © 2025 Pathak, O’Neill, Robinson, Hourwitz, Herr, Fourkas, Giniger and Losert. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wolfgang Losert, d2xvc2VydEB1bWQuZWR1