Click here to close
Hello! We notice that you are using Internet Explorer, which is not supported by Xenbase and may cause the site to display incorrectly.
We suggest using a current version of Chrome,
FireFox, or Safari.
Spatiotemporal development of coexisting wave domains of Rho activity in the cell cortex.
Hladyshau S
,
Kho M
,
Nie S
,
Tsygankov D
.
???displayArticle.abstract???
The Rho family GTPases are molecular switches that regulate cytoskeletal dynamics and cell movement through a complex spatiotemporal organization of their activity. In Patiria miniata (starfish) oocytes under in vitro experimental conditions (with overexpressed Ect2, induced expression of Ξ90 cyclin B, and roscovitine treatment), such activity generates multiple co-existing regions of coherent propagation of actin waves. Here we use computational modeling to investigate the development and properties of such wave domains. The model reveals that the formation of wave domains requires a balance between the activation and inhibition in the Rho signaling motif. Intriguingly, the development of the wave domains is preceded by a stage of low-activity quasi-static patterns, which may not be readily observed in experiments. Spatiotemporal patterns of this stage and the different paths of their destabilization define the behavior of the system in the later high-activity (observable) stage. Accounting for a strong intrinsic noise allowed us to achieve good quantitative agreement between simulated dynamics in different parameter regimes of the model and different wave dynamics in Patiria miniata and wild type Xenopus laevis (frog) data. For quantitative comparison of simulated and experimental results, we developed an automated method of wave domain detection, which revealed a sharp reversal in the process of pattern formation in starfish oocytes. Overall, our findings provide an insight into spatiotemporal regulation of complex and diverse but still computationally reproducible cell-level actin dynamics.
Figure 1. Cortical waves of Rho activity in Patiria miniata oocytes. Here oocytes have overexpressed Rho-GEF Ect2, induced expression of β90 cyclin B, and are treated with roscovitine. Both the experiment31 and our model show similar overall dynamics and the presence of wave domains. (A) Snapshots of wave propagation in the experiment. The images are generated from the data published in Bement et al.31 with the permission of the authors. The scale bar is 50 ΞΌm. (B) Snapshots of wave propagation in our simulation (0=0.2, 2=1, 1=3, 2=0.1). The cyan contours outline wave domains as defined in the following section, while the arrows indicate the direction of wave propagation in these domains. For the full dynamics see Supplementary Video 16 in Ref.31 and our Supplementary Video S1.
Figure 2. The model and an illustration of the wave domain concept. (A) Diagram of the signaling motif. π΄ and πΌ represent active and inactive forms of the signaling molecule, while πΉ represents a local inhibitor (see βMethodsβ for equations and Table 1 for parameter values). (B) The resulting distribution of component π΄ in a snapshot of the simulation with π0=0.15, π 2=0.8 using homogeneous initial conditions. (C) An image of two merged channels: green channel represents activator π΄ and red channel represents inhibitor πΉ. The concentrations were scaled as [πΆβπππ(πΆ)]/[max(πΆ)βπππ(πΆ)]. (D) The colormap of the direction angles, π, of the wave vectors at the same time point as B and C. Arrows represent the direction of wave vectors. (E) The colormap of the magnitude of a local change (gradient) in the directions of the wave vectors. Regions of high gradient are located at the edges of wave domains. (F) The colormap of our metric of the coherence distance. (G) Segmentation of the coherence distance with the modified watershed algorithm and after merging domains with similar values of the mean π. (H) The dynamics of wave domains. The number of regions detected by the automated segmentation pipeline is decreasing (left panel), while the mean area of regions is increasing over time (right panel).
Figure 3. Parameter scan and textural analysis of the emerging patterns. (A) The result of model simulations for a range of parameters π0 and π 2 (see also Supplementary Video S3). (B) The colormap of the excitation measure (see βMethodsβ) for the parameter space in (A). (C, D) The textural measures of pattern entropy and correlation for the parameter space in (A). Each measure was averaged over 100 simulation steps after the formation of waves from the initial homogeneous state.
Figure 4. Two stages of pattern development and characteristics of the regimes in the low-activity stage. (A,B) The amplitude of the activator pattern, max,βmin,, during the low- and high-activity stages for the quasi-static low-activity regime (0=0.2, 2=0.8) and the oscillatory low-activity regime (0=0.2, 2=1.1), respectively. (C,D) The same plots as in (A), (B), but in the log scale. (E,F) A series of snapshots of the system behavior from the same parameters as in (A,C) and (B,D), respectively (see also Supplementary Videos S4βS6). The concentrations were scaled as [βmin]/[maxβmin] both for the inhibitor (red channel) and the activator (green channel) before the channels were merged together. (G,H) (Left panels) Snapshots of the system in the quasi-static (0=0.2, 2=0.8) and oscillatory (0=0.2, 2=1.1) regimes of the low-activity stage; top and bottom panels respectively. (Central panels) Kymographs of the system dynamics along the vertical lines shown in the left panels. For comparison of patterns with different amplitude, the concentration values were scaled to zero mean and standard deviation equal to one (mean value was subtracted from each pixel and then divided by standard deviation). (Right panels) Temporal autocorrelation plots, showing the correlation coefficient as a function of the time lag (see βMethodsβ). (I) Classification of the parameter space with white and red colors indicating the quasi-static and oscillatory regimes, respectively. (J) Colormap showing the duration of the low-activity stage (in s).
Figure 5 The transition of the system between the low- and high-activity stages. (A) Snapshots of the patterns at low- and high- activity stages for different values of the inhibition strength: (1) weak (π0=0.18 , π 2=0.7 for quasi-static and π0=0.2, π 2=0.75 for oscillatory regimes), (2) balanced (π0=0.2, π 2=0.9 for quasi-static and π0=0.26, π 2=1.05 for oscillatory regimes), and (3) strong (π0=0.2, π 2=1.1 for quasi-static and π0=0.2, π 2=1.1 for oscillatory regimes). (B) Snapshots of the system at the beginning of the high-activity stage. Arrows indicate the direction of wave propagation. Dashed outlines indicate vector alignment on a larger scale. Images in (A), (B) show two merged channels: green for the activator and red for the inhibitor. (C) The activity amplitude and temporal autocorrelation for the oscillatory low-activity stage with weak negative feedback (left), balanced activation and inhibition (center), and strong negative feedback (right).
Figure 6. Classification and long-term dynamics of wave domains. (A) Different types of systemβs behaviors in the space of parameters representing the basal activation rate and the strength of negative feedback. (B) Textural entropy as a function of time during the growth of wave domains (π0=0.15, π 2=0.8). (C) Snapshots from the simulation in B with the simultaneous visualization of the inhibitor and activator in the system (merged channels). (D) Colormaps of the direction of wave propagation in (C) showing the increasing size and decreasing number of the wave domains over time. The values of wave vector angles were averaged with a running window of 83 s.
Figure 7. Wave domains analysis for experimental starfish data. (A) The pipeline of wave domain identification in starfish experimental images, including active Rho pattern, colormap of wave vector directions, colormap of the coherence distance, and wave domain segmentation. The scale bar is 50 ΞΌm. The image in the left panel is generated from the data published in Bement et al.31 with the permission of the authors. (B) The mean size of wave domains as a function of time for starfish oocyte. The grey area indicates the period of time after the sharp switch in the cell behavior.
Bement,
Activator-inhibitor coupling between Rho signalling and actin assembly makes the cell cortex an excitable medium.
2015, Pubmed,
Xenbase
Bement,
Activator-inhibitor coupling between Rho signalling and actin assembly makes the cell cortex an excitable medium.
2015,
Pubmed
,
Xenbase
Brauns,
Bulk-surface coupling identifies the mechanistic connection between Min-protein patterns in vivo and in vitro.
2021,
Pubmed
Chiou,
Principles that govern competition or co-existence in Rho-GTPase driven polarization.
2018,
Pubmed
Cusseddu,
A coupled bulk-surface model for cell polarisation.
2019,
Pubmed
Fukujin,
Self-organization of chemoattractant waves in Dictyostelium depends on F-actin and cell-substrate adhesion.
2016,
Pubmed
Gierer,
A theory of biological pattern formation.
1972,
Pubmed
Goryachev,
Compete or Coexist? Why the Same Mechanisms of Symmetry Breaking Can Yield Distinct Outcomes.
2020,
Pubmed
Goryachev,
Dynamics of Cdc42 network embodies a Turing-type mechanism of yeast cell polarity.
2008,
Pubmed
Goryachev,
How to make a static cytokinetic furrow out of traveling excitable waves.
2016,
Pubmed
Guan,
Rho GTPases and related signaling complexes in cell migration and invasion.
2020,
Pubmed
Hodge,
Regulating Rho GTPases and their regulators.
2016,
Pubmed
Holmes,
Regimes of wave type patterning driven by refractory actin feedback: transition from static polarization to dynamic wave behaviour.
2012,
Pubmed
Holmes,
Analysis of a minimal Rho-GTPase circuit regulating cell shape.
2016,
Pubmed
Kamps,
Optogenetic Tuning Reveals Rho Amplification-Dependent Dynamics of a Cell Contraction Signal Network.
2020,
Pubmed
Komatsu,
Development of an optimized backbone of FRET biosensors for kinases and GTPases.
2011,
Pubmed
Kondo,
Reaction-diffusion model as a framework for understanding biological pattern formation.
2010,
Pubmed
Kuo,
Inhibitory GEF phosphorylation provides negative feedback in the yeast polarity circuit.
2014,
Pubmed
Landino,
Rho and F-actin self-organize within an artificial cell cortex.
2021,
Pubmed
,
Xenbase
Liu,
Spots, stripes, and spiral waves in models for static and motile cells : GTPase patterns in cells.
2021,
Pubmed
Loose,
Spatial regulators for bacterial cell division self-organize into surface waves in vitro.
2008,
Pubmed
Machacek,
Coordination of Rho GTPase activities during cell protrusion.
2009,
Pubmed
Majumder,
Optogenetics enables real-time spatiotemporal control over spiral wave dynamics in an excitable cardiac system.
2018,
Pubmed
Marston,
Multiplexed GTPase and GEF biosensor imaging enables network connectivity analysis.
2020,
Pubmed
Meinhardt,
Orientation of chemotactic cells and growth cones: models and mechanisms.
1999,
Pubmed
Meinhardt,
Pattern formation by local self-activation and lateral inhibition.
2000,
Pubmed
Meinhardt,
Applications of a theory of biological pattern formation based on lateral inhibition.
1974,
Pubmed
,
Xenbase
Miao,
Wave patterns organize cellular protrusions and control cortical dynamics.
2019,
Pubmed
Mori,
Wave-pinning and cell polarity from a bistable reaction-diffusion system.
2008,
Pubmed
Otsuji,
A mass conserved reaction-diffusion system captures properties of cell polarity.
2007,
Pubmed
Ridley,
Rho GTPases and actin dynamics in membrane protrusions and vesicle trafficking.
2006,
Pubmed
Schneider,
NIH Image to ImageJ: 25 years of image analysis.
2012,
Pubmed
Shcherbakova,
Direct multiplex imaging and optogenetics of Rho GTPases enabled by near-infrared FRET.
2018,
Pubmed
Shoji,
Stripes, spots, or reversed spots in two-dimensional Turing systems.
2003,
Pubmed
Simon,
Pattern formation of Rho GTPases in single cell wound healing.
2013,
Pubmed
,
Xenbase
Taniguchi,
Phase geometries of two-dimensional excitable waves govern self-organized morphodynamics of amoeboid cells.
2013,
Pubmed
Turing,
The chemical basis of morphogenesis. 1953.
1990,
Pubmed
Tyson,
Cyclic AMP waves during aggregation of Dictyostelium amoebae.
1989,
Pubmed
de Beco,
Optogenetic dissection of Rac1 and Cdc42 gradient shaping.
2018,
Pubmed