Electric field-guided collective motility initiation of large epidermal cell groups.

Recent research has elucidated mechanochemical pathways of single cell polarization, but much less is known about collective motility initiation in adhesive cell groups. We used galvanotactic assays of zebrafish keratocyte cell groups, pharmacological perturbations, electric field switches, particle imaging velocimetry, and cell tracking to show that large cell groups initiate motility in minutes toward the cathode. Interestingly, while PI3K-inhibited single cells are biased toward the anode, inhibiting PI3K does not affect the cathode-directed cell group migration. We observed that control groups had the fastest cathode-migrating cell at the front, while the front cells in PI3K-inhibited groups were the slowest. Both control and PI3K-inhibited groups rapidly repolarized when the electric field direction was reversed, and the group migration continued after the electric field was switched off. Inhibiting myosin disrupted the cohesiveness of keratocyte groups and abolished the collective directionality and ability to switch direction when the electric field is reversed. Our data are consistent with a model according to which cells in the group sense the electric field individually and mechanical integration of the cells results in coherent group motility.

Supracellular organization confers directionality and mechanical potency to migrating pairs of cardiopharyngeal progenitor cells.

Physiological and pathological morphogenetic events involve a wide array of collective movements, suggesting that multicellular arrangements confer biochemical and biomechanical properties contributing to tissue-scale organization. The Ciona cardiopharyngeal progenitors provide the simplest model of collective cell migration, with cohesive bilateral cell pairs polarized along the leader-trailer migration path while moving between the ventral epidermis and trunk endoderm. We use the Cellular Potts Model to computationally probe the distributions of forces consistent with shapes and collective polarity of migrating cell pairs. Combining computational modeling, confocal microscopy, and molecular perturbations, we identify cardiopharyngeal progenitors as the simplest cell collective maintaining supracellular polarity with differential distributions of protrusive forces, cell-matrix adhesion, and myosin-based retraction forces along the leader-trailer axis. 4D simulations and experimental observations suggest that cell-cell communication helps establish a hierarchy to align collective polarity with the direction of migration, as observed with three or more cells in silico and in vivo. Our approach reveals emerging properties of the migrating collective: cell pairs are more persistent, migrating longer distances, and presumably with higher accuracy. Simulations suggest that cell pairs can overcome mechanical resistance of the trunk endoderm more effectively when they are polarized collectively. We propose that polarized supracellular organization of cardiopharyngeal progenitors confers emergent physical properties that determine mechanical interactions with their environment during morphogenesis.

Cellular memory in eukaryotic chemotaxis depends on the background chemoattractant concentration.

Cells of the social amoeba Dictyostelium discoideum migrate to a source of periodic traveling waves of chemoattractant as part of a self-organized aggregation process. An important part of this process is cellular memory, which enables cells to respond to the front of the wave and ignore the downward gradient in the back of the wave. During this aggregation, the background concentration of the chemoattractant gradually rises. In our microfluidic experiments, we exogenously applied periodic waves of chemoattractant with various background levels. We find that increasing background does not make detection of the wave more difficult, as would be naively expected. Instead, we see that the chemotactic efficiency significantly increases for intermediate values of the background concentration but decreases to almost zero for large values in a switch-like manner. These results are consistent with a computational model that contains a bistable memory module, along with a nonadaptive component. Within this model, an intermediate background level helps preserve directed migration by keeping the memory activated, but when the background level is higher, the directional stimulus from the wave is no longer sufficient to activate the bistable memory, suppressing directed migration. These results suggest that raising levels of chemoattractant background may facilitate the self-organized aggregation in Dictyostelium colonies.

Minimal Network Topologies for Signal Processing during Collective Cell Chemotaxis.

Cell-cell communication plays an important role in collective cell migration. However, it remains unclear how cells in a group cooperatively process external signals to determine the group's direction of motion. Although the topology of signaling pathways is vitally important in single-cell chemotaxis, the signaling topology for collective chemotaxis has not been systematically studied. Here, we combine mathematical analysis and simulations to find minimal network topologies for multicellular signal processing in collective chemotaxis. We focus on border cell cluster chemotaxis in the Drosophila egg chamber, in which responses to several experimental perturbations of the signaling network are known. Our minimal signaling network includes only four elements: a chemoattractant, the protein Rac (indicating cell activation), cell protrusion, and a hypothesized global factor responsible for cell-cell interaction. Experimental data on cell protrusion statistics allows us to systematically narrow the number of possible topologies from more than 40,000,000 to only six minimal topologies with six interactions between the four elements. This analysis does not require a specific functional form of the interactions, and only qualitative features are needed; it is thus robust to many modeling choices. Simulations of a stochastic biochemical model of border cell chemotaxis show that the qualitative selection procedure accurately determines which topologies are consistent with the experiment. We fit our model for all six proposed topologies; each produces results that are consistent with all experimentally available data. Finally, we suggest experiments to further discriminate possible pathway topologies.

The cAMP-induced G protein subunits dissociation monitored in live Dictyostelium cells by BRET reveals two activation rates, a positive effect of caffeine and potential role of microtubules.

To study the dynamics and mechanisms controlling activation of the heterotrimeric G protein Gα2ßγ in Dictyostelium in response to stimulation by the chemoattractant cyclic AMP (cAMP), we monitored the G protein subunit interaction in live cells using bioluminescence resonance energy transfer (BRET). We found that cAMP induces the cAR1-mediated dissociation of the G protein subunits to a similar extent in both undifferentiated and differentiated cells, suggesting that only a small number of cAR1 (as expressed in undifferentiated cells) is necessary to induce the full activation of Gα2ßγ. In addition, we found that treating cells with caffeine increases the potency of cAMP-induced Gα2ßγ activation; and that disrupting the microtubule network but not F-actin inhibits the cAMP-induced dissociation of Gα2ßγ. Thus, microtubules are necessary for efficient cAR1-mediated activation of the heterotrimeric G protein. Finally, kinetics analyses of Gα2ßγ subunit dissociation induced by different cAMP concentrations indicate that there are two distinct rates at which the heterotrimeric G protein subunits dissociate when cells are stimulated with cAMP concentrations above 500 nM versus only one rate at lower cAMP concentrations. Quantitative modeling suggests that the kinetics profile of Gα2ßγ subunit dissociation results from the presence of both uncoupled and G protein pre-coupled cAR1 that have differential affinities for cAMP and, consequently, induce G protein subunit dissociation through different rates. We suggest that these different signaling kinetic profiles may play an important role in initial chemoattractant gradient sensing.

Inferring single-cell behaviour from large-scale epithelial sheet migration patterns.

Cell migration plays an important role in a wide variety of biological processes and can incorporate both individual cell motion and collective behaviour. The emergent properties of collective migration are receiving increasing attention as collective motion's role in diseases such as metastatic cancer becomes clear. Yet, how individual cell behaviour influences large-scale, multi-cell collective motion remains unclear. In this study, we provide insight into the mechanisms behind collective migration by studying cell migration in a spreading monolayer of epithelial MCF10A cells. We quantify migration using particle image velocimetry and find that cell groups have features of motion that span multiple length scales. Comparing our experimental results to a model of collective cell migration, we find that cell migration within the monolayer can be affected in qualitatively different ways by cell motion at the boundary, yet it is not necessary to introduce leader cells at the boundary or specify other large-scale features to recapitulate this large-scale phenotype in simulations. Instead, in our model, collective motion can be enhanced by increasing the overall activity of the cells or by giving the cells a stronger coupling between their motion and polarity. This suggests that investigating the activity and polarity persistence of individual cells will add insight into the collective migration phenotypes observed during development and disease.

Relationship between cancer mutations and parameter sensitivity in Rb pathway.

It has long been known that formation of all sorts of tumors is largely owing to the genomic variations. Oncogenic mutations are often found focused on one or more important pathways which indicate that it is meaningful to investigate oncogenic mutations and oncogenic mechanisms from the point of view of biological network. Recently, we found that in apoptosis pathway of mammalian cell, mutations that cause large variations on the bifurcation point are more probably oncogenic mutations. Here, we used the Rb-E2F pathway in mammalian cell in response to growth factor as another example to verify this correlation. To conduct this study, nonlinear dynamics equations that describe the behavior of the Rb-E2F pathway was first constructed. Then we identified sensitive parameters which have a great influence on the system's bifurcation point. And we found that the sensitive parameters are highly related to high-frequency oncogenic mutations after comparing the results of parameter sensitivity analysis with profile of known cancer mutations. Moreover, the position of bifurcation point rather than concentration of a certain protein is a better measurement to determine biological network's function. Our results further confirm that nonlinear dynamics analysis of biological networks is an important way to understand oncogenesis. And the analysis method can become a powerful tool to understand and analyze the function of biological network.

Cellular memory in eukaryotic chemotaxis.

Natural chemical gradients to which cells respond chemotactically are often dynamic, with both spatial and temporal components. A primary example is the social amoeba Dictyostelium, which migrates to the source of traveling waves of chemoattractant as part of a self-organized aggregation process. Despite its physiological importance, little is known about how cells migrate directionally in response to traveling waves. The classic back-of-the-wave problem is how cells chemotax toward the wave source, even though the spatial gradient reverses direction in the back of the wave. Here, we address this problem by using microfluidics to expose cells to traveling waves of chemoattractant with varying periods. We find that cells exhibit memory and maintain directed motion toward the wave source in the back of the wave for the natural period of 6 min, but increasingly reverse direction for longer wave periods. Further insights into cellular memory are provided by experiments quantifying cell motion and localization of a directional-sensing marker after rapid gradient switches. The results can be explained by a model that couples adaptive directional sensing to bistable cellular memory. Our study shows how spatiotemporal cues can guide cell migration over large distances.

Correlation between oncogenic mutations and parameter sensitivity of the apoptosis pathway model.

One of the major breakthroughs in oncogenesis research in recent years is the discovery that, in most patients, oncogenic mutations are concentrated in a few core biological functional pathways. This discovery indicates that oncogenic mechanisms are highly related to the dynamics of biologic regulatory networks, which govern the behaviour of functional pathways. Here, we propose that oncogenic mutations found in different biological functional pathways are closely related to parameter sensitivity of the corresponding networks. To test this hypothesis, we focus on the DNA damage-induced apoptotic pathway--the most important safeguard against oncogenesis. We first built the regulatory network that governs the apoptosis pathway, and then translated the network into dynamics equations. Using sensitivity analysis of the network parameters and comparing the results with cancer gene mutation spectra, we found that parameters that significantly affect the bifurcation point correspond to high-frequency oncogenic mutations. This result shows that the position of the bifurcation point is a better measure of the functionality of a biological network than gene expression levels of certain key proteins. It further demonstrates the suitability of applying systems-level analysis to biological networks as opposed to studying genes or proteins in isolation.