Semi-automated quantification of filopodial dynamics

3-D Quantification of Filopodia in Motile Cancer Cells

We present a 3D bioimage analysis workflow to quantitatively analyze single, actin-stained cells with filopodial protrusions of diverse structural and temporal attributes, such as number, length, thickness, level of branching, and lifetime, in time-lapse confocal microscopy image data. Our workflow makes use of convolutional neural networks trained using real as well as synthetic image data, to segment the cell volumes with highly heterogeneous fluorescence intensity levels and to detect individual filopodial protrusions, followed by a constrained nearest-neighbor tracking algorithm to obtain valuable information about the spatio-temporal evolution of individual filopodia. We validated the workflow using real and synthetic 3-D time-lapse sequences of lung adenocarcinoma cells of three morphologically distinct filopodial phenotypes and show that it achieves reliable segmentation and tracking performance, providing a robust, reproducible and less time-consuming alternative to manual analysis of the 3D+t image data.

Automated profiling of growth cone heterogeneity defines relations between morphology and motility

Growth cones are complex, motile structures at the tip of an outgrowing neurite. They often exhibit a high density of filopodia (thin actin bundles), which complicates the unbiased quantification of their morphologies by software. Contemporary image processing methods require extensive tuning of segmentation parameters, require significant manual curation, and are often not sufficiently adaptable to capture morphology changes associated with switches in regulatory signals. To overcome these limitations, we developed Growth Cone Analyzer (GCA). GCA is designed to quantify growth cone morphodynamics from time-lapse sequences imaged both in vitro and in vivo, but is sufficiently generic that it may be applied to nonneuronal cellular structures. We demonstrate the adaptability of GCA through the analysis of growth cone morphological variation and its relation to motility in both an unperturbed system and in the context of modified Rho GTPase signaling. We find that perturbations inducing similar changes in neurite length exhibit underappreciated phenotypic nuance at the scale of the growth cone.

The tumor promoter-activated protein kinase Cs are a system for regulating filopodia

Different protein kinase C (PKC) isoforms have distinct roles in regulating cell functions. The conventional (α, β, γ) and novel (δ, ɛ, η, θ) classes are targets of phorbol ester tumor promoters, which are surrogates of endogenous second messenger, diacylglycerol. The promoter-stimulated disappearance of filopodia was investigated by use of blocking peptides (BPs) that inhibit PKC maturation and/or docking. Filopodia were partially rescued by a peptide representing PKC ɛ hydrophobic sequence, but also by a myristoylated PKC α/β pseudosubstrate sequence, and an inhibitor of T-cell protein tyrosine phosphatase (TC-PTP). The ability to turn over filopodia was widely distributed among PKC isoforms. PKC α and η hydrophobic sequences enhanced filopodia in cells in the absence of tumor promoter treatment. With transcriptional knockdown of PKC α, the content of PKC ɛ predominated over other isoforms. PKC ɛ could decrease filopodia significantly in promoter-treated cells, and this was attributed to ruffling. The presence of PKC α counteracted the PKC ɛ-mediated enhancement of ruffling. The results showed that there were two mechanisms of filopodia downregulation. One operated in the steady-state and relied on PKC α and η. The other was stimulated by tumor promoters and relied on PKC ɛ. Cycles of protrusion and retraction are characteristic of filopodia and are essential for the cell to orient itself during chemotaxis and haptotaxis. By suppressing filopodia, PKC ɛ can create a long-term "memory" of an environmental signal that may act in nature as a mnemonic device to mark the direction of a repulsive signal.

A Graphical User Interface for Software-assisted Tracking of Protein Concentration in Dynamic Cellular Protrusions

Filopodia are dynamic, finger-like cellular protrusions associated with migration and cell-cell communication. In order to better understand the complex signaling mechanisms underlying filopodial initiation, elongation and subsequent stabilization or retraction, it is crucial to determine the spatio-temporal protein activity in these dynamic structures. To analyze protein function in filopodia, we recently developed a semi-automated tracking algorithm that adapts to filopodial shape-changes, thus allowing parallel analysis of protrusion dynamics and relative protein concentration along the whole filopodial length. Here, we present a detailed step-by-step protocol for optimized cell handling, image acquisition and software analysis. We further provide instructions for the use of optional features during image analysis and data representation, as well as troubleshooting guidelines for all critical steps along the way. Finally, we also include a comparison of the described image analysis software with other programs available for filopodia quantification. Together, the presented protocol provides a framework for accurate analysis of protein dynamics in filopodial protrusions using image analysis software.

Quantifying Filopodia in Cultured Astrocytes by an Algorithm

Astrocytes in vivo extend thin processes termed peripheral astrocyte processes (PAPs), in particular around synapses where they can mediate glia-neuronal communication. The relation of PAPs to synapses is not based on coincidence, but it is not clear which stimuli and mechanisms lead to their formation and are active during process extension/ retraction in response to neuronal activity. Also, the molecular basis of the extremely fine PAP morphology (often 50 to 100 nm) is not understood. These open questions can be best investigated under in vitro conditions studying glial filopodia. We have previously analyzed filopodial mechanisms (Lavialle et al. PNAS 108:12915) applying an automated method for filopodia morphometry, which is now described in greater detail. The Filopodia Specific Shape Factor (FSSF) developed integrates number and length of filopodia. It quantifies filopodia independent of overall astrocytic shape or size, which can be intricate in itself. The algorithm supplied here permits automated image processing and measurements using ImageJ. Cells have to be sampled in higher numbers to obtain significant results. We validate the FSSF, and characterize the systematic influence of thresholding and camera pixel grid on measurements. We provide exemplary results of substance-induced filopodia dynamics (glutamate, mGluR agonists, EGF), and show that filopodia formation is highly sensitive to medium pH (CO₂) and duration of cell culture. Although the FSSF was developed to study astrocyte filopodia with focus on the perisynaptic glial sheath, we expect that this parameter can also be applied to neuronal growth cones, non-neural cell types, or cell lines.

Automated analysis of filopodial length and spatially resolved protein concentration via adaptive shape tracking

A novel approach based on the convex-hull algorithm is used for parallel analysis of growth dynamics and relative spatiotemporal protein concentration along flexible filopodial protrusions. Testing of filopodia formation in silico, in vitro, and in vivo validates the robustness and sensitivity of the proposed approach.

Filopodia are dynamic, actin-rich structures that transiently form on a variety of cell types. To understand the underlying control mechanisms requires precise monitoring of localization and concentration of individual regulatory and structural proteins as filopodia elongate and subsequently retract. Although several methods exist that analyze changes in filopodial shape, a software solution to reliably correlate growth dynamics with spatially resolved protein concentration along the filopodium independent of bending, lateral shift, or tilting is missing. Here we introduce a novel approach based on the convex-hull algorithm for parallel analysis of growth dynamics and relative spatiotemporal protein concentration along flexible filopodial protrusions. Detailed in silico tests using various geometries confirm that our technique accurately tracks growth dynamics and relative protein concentration along the filopodial length for a broad range of signal distributions. To validate our technique in living cells, we measure filopodial dynamics and quantify spatiotemporal localization of filopodia-associated proteins during the filopodial extension–retraction cycle in a variety of cell types in vitro and in vivo. Together these results show that the technique is suitable for simultaneous analysis of growth dynamics and spatiotemporal protein enrichment along filopodia. To allow readily application by other laboratories, we share source code and instructions for software handling.

The dynamics of growth cone morphology

BACKGROUND

Normal brain function depends on the development of appropriate patterns of neural connections. A critical role in guiding axons to their targets during neural development is played by neuronal growth cones. These have a complex and rapidly changing morphology; however, a quantitative understanding of this morphology, its dynamics and how these are related to growth cone movement, is lacking.

RESULTS

Here we use eigenshape analysis (principal components analysis in shape space) to uncover the set of five to six basic shape modes that capture the most variance in growth cone form. By analysing how the projections of growth cones onto these principal modes evolve in time, we found that growth cone shape oscillates with a mean period of 30 min. The variability of oscillation periods and strengths between different growth cones was correlated with their forward movement, such that growth cones with strong, fast shape oscillations tended to extend faster. A simple computational model of growth cone shape dynamics based on dynamic microtubule instability was able to reproduce quantitatively both the mean and variance of oscillation periods seen experimentally, suggesting that the principal driver of growth cone shape oscillations may be intrinsic periodicity in cytoskeletal rearrangements.

CONCLUSIONS

Intrinsically driven shape oscillations are an important component of growth cone shape dynamics. More generally, eigenshape analysis has the potential to provide new quantitative information about differences in growth cone behaviour in different conditions.

Promoting filopodial elongation in neurons by membrane-bound magnetic nanoparticles

Filopodia are 5-10 μm long processes that elongate by actin polymerization, and promote axon growth and guidance by exerting mechanical tension and by molecular signaling. Although axons elongate in response to mechanical tension, the structural and functional effects of tension specifically applied to growth cone filopodia are unknown. Here we developed a strategy to apply tension specifically to retinal ganglion cell (RGC) growth cone filopodia through surface-functionalized, membrane-targeted superparamagnetic iron oxide nanoparticles (SPIONs). When magnetic fields were applied to surface-bound SPIONs, RGC filopodia elongated directionally, contained polymerized actin filaments, and generated retrograde forces, behaving as bona fide filopodia. Data presented here support the premise that mechanical tension induces filopodia growth but counter the hypothesis that filopodial tension directly promotes growth cone advance. Future applications of these approaches may be used to induce sustained forces on multiple filopodia or other subcellular microstructures to study axon growth or cell migration. From the clinical editor: Mechanical tension to the tip of filopodia is known to promote axonal growth. In this article, the authors used superparamagnetic iron oxide nanoparticles (SPIONs) targeted specifically to membrane molecules, then applied external magnetic field to elicit filopodial elongation, which provided a tool to study the role of mechanical forces in filopodia dynamics and function.

F-dynamics: automated quantification of dendrite filopodia dynamics in living neurons

BACKGROUND

Dendritic filopodia are highly motile and flexible protrusions that explore the surroundings in search for an appropriate presynaptic partner. Dendritic filopodia morphologically and functionally transform into postsynaptic dendritic spines, once the appropriate partner has been chosen. Therefore, proper formation of synapses depends on the dynamics of dendritic filopodia and spines. Thus, a rigorous assessment of dendrite filopodia behavior could be informative in providing a link between filopodia dynamics and synaptic development.

NEW METHOD

In this paper, a tool for automated tracking of filopodia dynamics, the Filopodia-dynamics program (F-dynamics), will be described, tested and applied. The aim of this study is to validate the accuracy and reliability of F-dynamics and to test the program in live neurons.

RESULTS

We demonstrate that filopodia dynamics can be reliably and accurately quantified using the F-dynamics program. In the present study, this program was used to successfully show that lithium treatment increases filopodia motility.

COMPARISON WITH EXCITING METHODS

F-dynamics is the first analysis program that is able to determine dendritic filopodia dynamics automatically across both the longitudinal and lateral dimensions.

CONCLUSION

Our data suggests that this analysis method can be used to differentiate between different experimental conditions and illustrates the potential of the program to measure pharmaceutical or genetic effects on filopodia dynamics.

CellGeo: a computational platform for the analysis of shape changes in cells with complex geometries

The open source MATLAB application CellGeo is a user-friendly computational platform that allows simultaneous, automated tracking and analysis of dynamic changes in cell shape, including protrusions ranging from filopodia to lamellipodia to growth cones.

Cell biologists increasingly rely on computer-aided image analysis, allowing them to collect precise, unbiased quantitative results. However, despite great progress in image processing and computer vision, current computational approaches fail to address many key aspects of cell behavior, including the cell protrusions that guide cell migration and drive morphogenesis. We developed the open source MATLAB application CellGeo, a user-friendly computational platform to allow simultaneous, automated tracking and analysis of dynamic changes in cell shape, including protrusions ranging from filopodia to lamellipodia. Our method maps an arbitrary cell shape onto a tree graph that, unlike traditional skeletonization algorithms, preserves complex boundary features. CellGeo allows rigorous but flexible definition and accurate automated detection and tracking of geometric features of interest. We demonstrate CellGeo’s utility by deriving new insights into (a) the roles of Diaphanous, Enabled, and Capping protein in regulating filopodia and lamellipodia dynamics in Drosophila melanogaster cells and (b) the dynamic properties of growth cones in catecholaminergic a–differentiated neuroblastoma cells.

Improved method for the quantification of motility in glia and other morphologically complex cells

Cells such as astrocytes and radial glia with many densely ramified, fine processes pose particular challenges for the quantification of structural motility. Here we report the development of a method to calculate a motility index for individual cells with complex, dynamic morphologies. This motility index relies on boxcar averaging of the difference images generated by subtraction of images collected at consecutive time points. An image preprocessing step involving 2D projection, edge detection, and dilation of the raw images is first applied in order to binarize the images. The boxcar averaging of difference images diminishes the impact of artifactual pixel fluctuations while accentuating the group-wise changes in pixel values which are more likely to represent real biological movement. Importantly, this provides a value that correlates with mean process elongation and retraction rates without requiring detailed reconstructions of very complex cells. We also demonstrate that additional increases in the sensitivity of the method can be obtained by denoising images using the temporal frequency power spectra, based on the fact that rapid intensity fluctuations over time are mainly due to imaging artifact. The MATLAB programs implementing these motility analysis methods, complete with user-friendly graphical interfaces, have been made publicly available for download.

FiloDetect: automatic detection of filopodia from fluorescence microscopy images

Background

Filopodia are small cellular projections that help cells to move through and sense their environment. Filopodia play crucial roles in processes such as development and wound-healing. Also, increases in filopodia number or size are characteristic of many invasive cancers and are correlated with increased rates of metastasis in mouse experiments. Thus, one possible route to developing anti-metastatic therapies is to target factors that influence the filopodia system. Filopodia can be detected by eye using confocal fluorescence microscopy, and they can be manually annotated in images to quantify filopodia parameters. Although this approach is accurate, it is slow, tedious and not entirely objective. Manual detection is a significant barrier to the discovery and quantification of new factors that influence the filopodia system.

Results

Here, we present FiloDetect, an automated tool for detecting, counting and measuring the length of filopodia in fluorescence microscopy images. The method first segments the cell from the background, using a modified triangle threshold method, and then extracts the filopodia using a series of morphological operations. We verified the accuracy of FiloDetect on Rat2 and B16F1 cell images from three different labs, showing that per-cell filopodia counts and length estimates are highly correlated with the manual annotations. We then used FiloDetect to assess the role of a lipid kinase on filopodia production in breast cancer cells. Experimental results show that PI4KIII β expression leads to an increase in filopodia number and length, suggesting that PI4KIII β is involved in driving filopodia production.

Conclusion

FiloDetect provides accurate and objective quantification of filopodia in microscopy images, and will enable large scale comparative studies to assess the effects of different genetic and chemical perturbations on filopodia production in different cell types, including cancer cell lines.

Fluorescent saxitoxins for live cell imaging of single voltage-gated sodium ion channels beyond the optical diffraction limit

A desire to better understand the role of voltage-gated sodium channels (Na(V)s) in signal conduction and their dysregulation in specific disease states motivates the development of high precision tools for their study. Nature has evolved a collection of small molecule agents, including the shellfish poison (+)-saxitoxin, that bind to the extracellular pore of select Na(V) isoforms. As described in this report, de novo chemical synthesis has enabled the preparation of fluorescently labeled derivatives of (+)-saxitoxin, STX-Cy5, and STX-DCDHF, which display reversible binding to Na(V)s in live cells. Electrophysiology and confocal fluorescence microscopy studies confirm that these STX-based dyes function as potent and selective Na(V) labels. The utility of these probes is underscored in single-molecule and super-resolution imaging experiments, which reveal Na(V) distributions well beyond the optical diffraction limit in subcellular features such as neuritic spines and filopodia.

High-content neurite development study using optically patterned substrates

The study of neurite guidance in vitro relies on the ability to reproduce the distribution of attractive and repulsive guidance molecules normally expressed in vivo. The identification of subtle variations in the neurite response to changes in the spatial distribution of extracellular molecules can be achieved by monitoring the behavior of cells on protein gradients. To do this, automated high-content screening assays are needed to quantify the morphological changes resulting from growth on gradients of guidance molecules. Here, we present the use of laser-assisted protein adsorption by photobleaching (LAPAP) to allow the fabrication of large-scale substrate-bound laminin-1 gradients to study neurite extension. We produced thousands of gradients of different slopes and analyzed the variations in neurite attraction of neuron-like cells (RGC-5). An image analysis algorithm processed bright field microscopy images, detecting each cell and quantifying the soma centroid and the initiation, terminal and turning angles of the longest neurite.

Protein fluxes along the filopodium as a framework for understanding the growth-retraction dynamics: the interplay between diffusion and active transport

We present a picture of filopodial growth and retraction from physics perspective, where we emphasize the significance of the role played by protein fluxes due to spatially extended nature of the filopodium. We review a series of works, which used stochastic simulations and mean field analytical modeling to find the concentration profile of G-actin inside a filopodium, which, in turn, determines the stationary filopodial length. In addition to extensively reviewing the prior works, we also report some new results on the role of active transport in regulating the length of filopodia. We model a filopodium where delivery of actin monomers towards the tip can occur both through passive diffusion and active transport by myosin motors. We found that the concentration profile of G-actin along the filopodium is rather non-trivial, containing a narrow minimum near the base followed by a broad maximum. For efficient enough actin transport, this non-monotonous shape is expected to occur under a broad set of conditions. We also raise the issue of slow approach to the stationary length and the possibility of multiple steady state solutions.

Cofilin cooperates with fascin to disassemble filopodial actin filaments

Cells use a large repertoire of proteins to remodel the actin cytoskeleton. Depending on the proteins involved, F-actin is organized in specialized protrusions such as lamellipodia or filopodia, which serve diverse functions in cell migration and sensing. Although factors responsible for directed filament assembly in filopodia have been extensively characterized, the mechanisms of filament disassembly in these structures are mostly unknown. We investigated how the actin-depolymerizing factor cofilin-1 affects the dynamics of fascincrosslinked actin filaments in vitro and in live cells. By multicolor total internal reflection fluorescence microscopy and fluorimetric assays, we found that cofilin-mediated severing is enhanced in fascin-crosslinked bundles compared with isolated filaments, and that fascin and cofilin act synergistically in filament severing. Immunolabeling experiments demonstrated for the first time that besides its known localization in lamellipodia and membrane ruffles, endogenous cofilin can also accumulate in the tips and shafts of filopodia. Live-cell imaging of fluorescently tagged proteins revealed that cofilin is specifically targeted to filopodia upon stalling of protrusion and during their retraction. Subsequent electron tomography established filopodial actin filament and/or bundle fragmentation to precisely correlate with cofilin accumulation. These results identify a new mechanism of filopodium disassembly involving both fascin and cofilin.

Efficacy of synaptic inhibition depends on multiple, dynamically interacting mechanisms implicated in chloride homeostasis

Chloride homeostasis is a critical determinant of the strength and robustness of inhibition mediated by GABA(A) receptors (GABA(A)Rs). The impact of changes in steady state Cl(-) gradient is relatively straightforward to understand, but how dynamic interplay between Cl(-) influx, diffusion, extrusion and interaction with other ion species affects synaptic signaling remains uncertain. Here we used electrodiffusion modeling to investigate the nonlinear interactions between these processes. Results demonstrate that diffusion is crucial for redistributing intracellular Cl(-) load on a fast time scale, whereas Cl(-)extrusion controls steady state levels. Interaction between diffusion and extrusion can result in a somato-dendritic Cl(-) gradient even when KCC2 is distributed uniformly across the cell. Reducing KCC2 activity led to decreased efficacy of GABA(A)R-mediated inhibition, but increasing GABA(A)R input failed to fully compensate for this form of disinhibition because of activity-dependent accumulation of Cl(-). Furthermore, if spiking persisted despite the presence of GABA(A)R input, Cl(-) accumulation became accelerated because of the large Cl(-) driving force that occurs during spikes. The resulting positive feedback loop caused catastrophic failure of inhibition. Simulations also revealed other feedback loops, such as competition between Cl(-) and pH regulation. Several model predictions were tested and confirmed by [Cl(-)](i) imaging experiments. Our study has thus uncovered how Cl(-) regulation depends on a multiplicity of dynamically interacting mechanisms. Furthermore, the model revealed that enhancing KCC2 activity beyond normal levels did not negatively impact firing frequency or cause overt extracellular K(-) accumulation, demonstrating that enhancing KCC2 activity is a valid strategy for therapeutic intervention.

Engineered in vitro/in silico models to examine neurite target preference

Research on spinal cord injury (SCI) repair focuses on developing mechanisms to allow neurites to grow past an injury site. In this article, we observe that numerous divergent paths (i.e., spinal roots) are present along the spinal column, and hence guidance strategies must be devised to ensure that regrowing neurites reach viable targets. Therefore, we have engineered an in vitro micropatterned model in which cultured E7 dorsal root ganglia (DRG) explants may enter alternate pathways (?roots?) along a branching micropattern. Alongside this in vitro model, we have developed an in silico simulation that we validate by comparison with independent experiments. We find in both in silico and in vitro models that the probability of a neurite entering a given root decreases exponentially with respect to the number of roots away from the DRG; consequently, the likelihood of neurites reaching a distant root can be vanishingly small. This result represents a starting point for future strategies to optimize the likelihood that neurites will reach appropriate targets in the regenerating nervous system, and provides a new computational tool to evaluate the feasibility and expected success of neurite guidance in complex geometries.

Fasciculation and defasciculation of neurite bundles on micropatterned substrates

We describe experiments of fasciculation, i.e., bundling, of chick sensory neurites on 2D striped substrates. By Fourier decomposition, we separate left-going and right-going neurite components from in vitro images, and we find first that neurite bundles orient toward preferred angles with respect to the stripe direction, and second that in vitro bundles travel in leftward and rightward directions nearly uninterrupted by crossings of bundles traveling in the opposing direction. We explore mechanisms that lead to these behaviors, and summarize implications for future models for neurite outgrowth and guidance.

Quimp3, an automated pseudopod-tracking algorithm

To understand movement of amoeboid cells we have developed an information tool that automatically detects protrusions of moving cells. The algorithm uses digitized cell recordings at a speed of approximately 1 image per second that are analyzed in three steps. In the first part, the outline of a cell is defined as a polygon of approximately 150 nodes, using the previously published Quimp2 program. By comparing the position of the nodes in place and time, each node contains information on position, local curvature and speed of movement. The second part uses rules for curvature and movement to define the position and time of start and end of a growing pseudopod. This part of the algorithm produces quantitative data on size, surface area, lifetime, frequency and direction of pseudopod extension. The third part of the algorithm assigns qualitative properties to each pseudopod. It decides on the origin of a pseudopod as splitting of an existing pseudopod or as extension de novo. It also decides on the fate of each pseudopod as merged with the cell body or retracted. Here we describe the pseudopod tool and present the first data based on the analysis of approximately 1,000 pseudopodia extended by Dictyostelium cells in the absence of external cues.