Quantification of Local Morphodynamics and Local GTPase Activity by Edge Evolution Tracking

Advances in time-lapse fluorescence microscopy have enabled us to directly observe dynamic cellular phenomena. Although the techniques themselves have promoted the understanding of dynamic cellular functions, the vast number of images acquired has generated a need for automated processing tools to extract statistical information. A problem underlying the analysis of time-lapse cell images is the lack of rigorous methods to extract morphodynamic properties. Here, we propose an algorithm called edge evolution tracking (EET) to quantify the relationship between local morphological changes and local fluorescence intensities around a cell edge using time-lapse microscopy images. This algorithm enables us to trace the local edge extension and contraction by defining subdivided edges and their corresponding positions in successive frames. Thus, this algorithm enables the investigation of cross-correlations between local morphological changes and local intensity of fluorescent signals by considering the time shifts. By applying EET to fluorescence resonance energy transfer images of the Rho-family GTPases Rac1, Cdc42, and RhoA, we examined the cross-correlation between the local area difference and GTPase activity. The calculated correlations changed with time-shifts as expected, but surprisingly, the peak of the correlation coefficients appeared with a 6–8 min time shift of morphological changes and preceded the Rac1 or Cdc42 activities. Our method enables the quantification of the dynamics of local morphological change and local protein activity and statistical investigation of the relationship between them by considering time shifts in the relationship. Thus, this algorithm extends the value of time-lapse imaging data to better understand dynamics of cellular function

Distinct predictive performance of Rac1 and Cdc42 in cell migration.

We propose a new computation-based approach for elucidating how signaling molecules are decoded in cell migration. In this approach, we performed FRET time-lapse imaging of Rac1 and Cdc42, members of Rho GTPases which are responsible for cell motility, and quantitatively identified the response functions that describe the conversion from the molecular activities to the morphological changes. Based on the identified response functions, we clarified the profiles of how the morphology spatiotemporally changes in response to local and transient activation of Rac1 and Cdc42, and found that Rac1 and Cdc42 activation triggers laterally propagating membrane protrusion. The response functions were also endowed with property of differentiator, which is beneficial for maintaining sensitivity under adaptation to the mean level of input. Using the response function, we could predict the morphological change from molecular activity, and its predictive performance provides a new quantitative measure of how much the Rho GTPases participate in the cell migration. Interestingly, we discovered distinct predictive performance of Rac1 and Cdc42 depending on the migration modes, indicating that Rac1 and Cdc42 contribute to persistent and random migration, respectively. Thus, our proposed predictive approach enabled us to uncover the hidden information processing rules of Rho GTPases in the cell migration

Computer vision profiling of neurite outgrowth mordphodynamics reveals spatio-temporal modularity of Rho GTPase signaling

Neurite outgrowth is essential to build the neuronal processes that produce axons and dendrites that connect the adult brain. In cultured cells, the neurite outgrowth process is highly dynamic, and consists of a series of repetitive morphogenetic sub-processes (MSPs), such as neurite initiation, elongation, branching, growth cone motility and collapse (da Silva and Dotti 2002). Neurons also actively migrate, which might in part reflect neuronal migration during brain development. Each of the different MSPs inherent to neurite outgrowth and cell migration is likely to be regulated by precise spatio-temporal signaling networks that control cytoskeletal dynamics, trafficking and adhesion events. These MSPs can occur on a range of time and length scales. For example, microtubule bundling in the neurite shaft can be maintained during hours, while growth cone filopodia dynamically explore their surrounding on time scales of seconds and length scales of single microns. This implies that a correct understanding of these processes will require analysis with an adequate spatio-temporal resolution. The Rho family of GTPases are signaling switches that regulate a wide variety of cellular processes, such as actin and adhesion dynamics, gene transcription, and neuronal differentiation (Boguski and McCormick 1993). Rho GTPases are activated by guanine nucleotide exchange factors (GEFs), and are switched off by GTPase activating proteins (GAPs). Upon activation, Rho GTPases can associate with effectors to initiate a downstream response. Current models propose that Rac1 and Cdc42 regulate neurite extension, while RhoA controls growth cone collapse and neurite retraction (da Silva and Dotti 2002). However, until now the effects of Rho GTPases on neurite outgrowth have mostly been assessed using protein mutants in steady-state experiments, most often at late differentiation stages, which do not provide any insight about the different MSPs during neurite outgrowth. However, our proteomic analysis of biochemically-purified neurites from N1E-115 neuronal-like cells (Pertz et al. 2008), has suggested the existence of an unexpectedly complex 220 proteins signaling network consisting of multiple GEFs, GAPs, Rho GTPases, effectors and additional interactors. This is inconsistent with the simplistic view that classical experiments have provided before. In order to gain insight into the complexity of this Rho GTPase signaling network, we performed a siRNA screen that targets each of these 220 proteins individually. We hypothesized that specific spatio-temporal Rho GTPase signaling networks control different MSPs occurring during neurite outgrowth, and therefore designed an integrated approach to capture the whole morphodynamic continuum of this process. Perturbations of candidates that lead to a similar phenotype might be part of a given spatio-temporal signaling network. This approach consisted of: 1) A high content microscopy platform that allowed us to produce 8000 timelapse movies of 660 siRNA perturbations; 2) A custom built, computer vision approach that allowed us to automatically segment and track neurite and soma morphodynamics in the timelapse movies (collaboration with the group of Pascal Fua, EPFL, Lausanne); 3) A sophisticated statistical analysis pipeline that allowed the extraction of morphological and morphodynamic signatures (MDSs) relevant to each siRNA perturbation (collaboration with the group of Francois Fleuret, IDIAP). Analysis of our dataset revealed that each siRNA perturbation led to a quantifiable phenotype, emphasizing the quality of our proteomic dataset. Hierarchical clustering of the MDSs revealed the existence of 24 phenoclusters that provide information about neurite length, branching, number of neurites, soma migration speed, and a panel of additional morphological and morphodynamic features that are more difficult to grasp using visual inspection. This complex phenotypic space can more easily be understood when classified according to the first 4 features. Our screen then suggests the existence of 4 major morphodynamic phenotypes that define distinct stages of the neurite outgrowth process. These consist of phenotypes with short neurites, multiple short neurites, long neurites, and long and branched neurites. Further subdivision using the other features provides more information, with cell migration features being very often affected. This implies a high overlap between the signaling machinery that regulates the neurite outgrowth and cell migration processes. The high phenotypical redundancy (24 clusters for 220 candidate genes) provides only limited information to deduce unambiguous signaling networks regulating distinct MSPs. Further knowledge acquired from other approaches we used to study Rho GTPase signaling (FRET biosensors, and other live cell imaging techniques), made us realize that some morphodynamic phenotypes can only be understood when growth cone dynamics are inspected at a much higher resolution. For this purpose, we decided to further investigate a defined subset of genes using high resolution live cell imaging and a custom built growth cone segmentation and tracking pipeline for accurate quantification (collaboration with the group of Gaudenz Danuser, Harvard Medical School, Boston). These results shed light into how distinct cytoskeletal networks enabling growth cone advance can globally impact the neurite outgrowth process. A clear understanding of spatio-temporal Rho GTPase signaling will therefore require multi-scale approaches. Our results provide the first insight into the complexity of spatio-temporal Rho GTPase signaling during neurite outgrowth. The technologies we devised and our initial results, pave the way for a systems biology understanding of these complex signaling systems

Lamellipodin tunes cell migration by stabilizing protrusions and promoting adhesion formation

Efficient migration on adhesive surfaces involves the protrusion of lamellipodial actin networks and their subsequent stabilization by nascent adhesions. The actin-binding protein lamellipodin (Lpd) is thought to play a critical role in lamellipodium protrusion, by delivering Ena/VASP proteins onto the growing plus ends of actin filaments and by interacting with the WAVE regulatory complex, an activator of the Arp2/3 complex, at the leading edge. Using B16-F1 melanoma cell lines, we demonstrate that genetic ablation of Lpd compromises protrusion efficiency and coincident cell migration without altering essential parameters of lamellipodia, including their maximal rate of forward advancement and actin polymerization. We also confirmed lamellipodia and migration phenotypes with CRISPR/Cas9-mediated Lpd knockout Rat2 fibroblasts, excluding cell type-specific effects. Moreover, computer-aided analysis of cell-edge morphodynamics on B16-F1 cell lamellipodia revealed that loss of Lpd correlates with reduced temporal protrusion maintenance as a prerequisite of nascent adhesion formation. We conclude that Lpd optimizes protrusion and nascent adhesion formation by counteracting frequent, chaotic retraction and membrane ruffling.This article has an associated First Person interview with the first author of the paper

The Dynamics of Dorsal Actin Waves

The recent years have shown that waves of actin polyermization are central to the morphodynamics of cells. This thesis is dedicated to deciphering of the propagation mechanism underlying actin waves known as Circular Dorsal Ruffles (CDRs). While these ring-shaped undulations on the dorsal cell side have been known to the biological community for several decades the mechanism underlying their formation and propagation has remained a puzzle. It is the hypothesis of this work that CDRs can be described as waves that form and propagate in an active medium that is constituted by the actin machinery of the cell. The identification of the corresponding functional elements is the aim of this work. For this, the structure, morphology and dynamics of CDRs are investigated in detail and with a view that is guided by the typical structure of models of active media. Throughout the whole thesis, the FitzHugh-Nagumo system serves as a prototype model for the explanation of the mechanisms underlying the phenomena observed for CDRs on an abstract level. Novel results are presented regarding the identification of the processes of actin dynamics within CDRs and their compartmentalization. The systematic analysis of the dynamics of CDR wavefronts reveals that they exhibit a number of previously unknown phenomena, among them breathing modes, spiral waves, and collision annihilation. All these features are well founded in the framework of active media. Since the dynamics of CDRs strongly depends on the cellular morphology, a novel method for their investigation is developed in which cells are forced into disc-shapes via microcontact printing for a quantitative analysis of data of identically shaped cells. This framework allows for direct comparability to numerical studies, which reveals that stochastic elements in protein dynamics are key for the understanding of CDRs

Automated characterization of cell shape changes during amoeboid motility by skeletonization

<p>Abstract</p> <p>Background</p> <p>The ability of a cell to change shape is crucial for the proper function of many cellular processes, including cell migration. One type of cell migration, referred to as amoeboid motility, involves alternating cycles of morphological expansion and retraction. Traditionally, this process has been characterized by a number of parameters providing global information about shape changes, which are insufficient to distinguish phenotypes based on local pseudopodial activities that typify amoeboid motility.</p> <p>Results</p> <p>We developed a method that automatically detects and characterizes pseudopodial behavior of cells. The method uses skeletonization, a technique from morphological image processing to reduce a shape into a series of connected lines. It involves a series of automatic algorithms including image segmentation, boundary smoothing, skeletonization and branch pruning, and takes into account the cell shape changes between successive frames to detect protrusion and retraction activities. In addition, the activities are clustered into different groups, each representing the protruding and retracting history of an individual pseudopod.</p> <p>Conclusions</p> <p>We illustrate the algorithms on movies of chemotaxing <it>Dictyostelium </it>cells and show that our method makes it possible to capture the spatial and temporal dynamics as well as the stochastic features of the pseudopodial behavior. Thus, the method provides a powerful tool for investigating amoeboid motility.</p

FRET binding antenna reports spatiotemporal dynamics of GDI-Cdc42 GTPase interactions

Guanine-nucleotide dissociation inhibitors (GDI) are negative regulators of Rho family GTPases that sequester the GTPases away from the membrane. Here we ask how GDI-Cdc42 interaction regulates localized Cdc42 activation for cell motility. The sensitivity of cells to overexpression of Rho family pathway components led us to a new biosensor design (GDI.Cdc42 FLARE), in which Cdc42 was modified with a FRET ‘binding antenna’ that selectively reported Cdc42 binding to endogenous GDI. Similar antennae could also report GDI-Rac1 and GDI-RhoA interaction. Through computational multiplexing and simultaneous imaging, we determined the spatiotemporal dynamics of GDI-Cdc42 interaction and Cdc42 activation during cell protrusion and retraction. This revealed a remarkably tight coordination of GTPase release and activation on a time scale of 10 seconds, suggesting that GDI-Cdc42 interactions are a critical component in the spatiotemporal regulation of Cdc42 activity, and not merely a mechanism for global sequestration of an inactivated pool of signaling molecules

Amoeboid Transition Occurs in Mammilian Tumor Cells in Response to Changes in Spacial Confinement and Adhesion

In this dissertation, we review how plasticity in the modes of cell migration can occur in response to changes in the extracellular matrix. Different modes of migration exhibit varying characteristics, such as cell adhesion, membrane protrusion, and proteolysis, allowing tumor cells to adapt to their current extracellular environment, thus enhancing invasive behavior in response to different antitumorigenic therapies. The combination of various cell migration characteristics results in a distinct mode of single-cell migration, which, for single-cell migration, falls under two general modes, mesenchymal-traction-force motility or amoeboid-propulsion-squeeze motility. Mesenchymal-to-amoeboid transition is known to occur when tumor cells are in a softer matrix with low adhesion and in a high-confinement environment. Under these conditions the tumor cells exhibit higher levels of cortical myosin activity. Amoeboid migration is characterized by high cortical contractility, mediated by increased myosin activity along the cell cortex. Understanding how contractility is increased under changes in adhesion and confinement is important to uncover the possible mechanisms a tumor cell uses to optimize motility. Here, I will introduce a model that explores how the differential regulation of RhoGTPases (Rac1 and RhoA) is modulated with changes in cell-matrix adhesion and cell-matrix confinement to induce mesenchymal-to-amoeboid-transition. In the model, tumor cells in soft matrix are exposed to fewer ligands to which they can bind their integrins and activate Rac1. Loss of Rac1 activation in soft matrix inhibits lamellipodia formation, and the double-negative relationship between Rac1 and RhoA will shift towards RhoA to promote increase cortical myosin activity for amoeboid migration (Figure 1). Myosin activity is further enhanced under confinement through retrograde flow and recycling of transmembrane protein like integrin and syndecan, which results in increase of RhoA-mediated contractility for amoeboid migration

A requirement for filopodia extension toward Slit during Robo-mediated axon repulsion

Axons navigate long distances through complex 3D environments to interconnect the nervous system during development. Although the precise spatiotemporal effects of most axon guidance cues remain poorly characterized, a prevailing model posits that attractive guidance cues stimulate actin polymerization in neuronal growth cones whereas repulsive cues induce actin disassembly. Contrary to this model, we find that the repulsive guidance cue Slit stimulates the formation and elongation of actin-based filopodia from mouse dorsal root ganglion growth cones. Surprisingly, filopodia form and elongate toward sources of Slit, a response that we find is required for subsequent axonal repulsion away from Slit. Mechanistically, Slit evokes changes in filopodium dynamics by increasing direct binding of its receptor, Robo, to members of the actin-regulatory Ena/VASP family. Perturbing filopodium dynamics pharmacologically or genetically disrupts Slit-mediated repulsion and produces severe axon guidance defects in vivo. Thus, Slit locally stimulates directional filopodial extension, a process that is required for subsequent axonal repulsion downstream of the Robo receptor.National Institutes of Health (U.S.) (Grant F32-CA165700)National Institutes of Health (U.S.) (Grant R01-GM068678)National Institutes of Health (U.S.) (Grant P30-CA014051