Emergence of chiral multi-armed spirals in an open system of migrating cells under continuous cell supply

Supporting data for Emergence of chiral multi-armed spirals in an open system of migrating cells under continuous cell supply Masayuki Hayakawa, Biplab Bhattacherjee, Hidekazu Kuwayama, Tatsuo Shibata Corresponding author: Tatsuo Shibata (tatsuo.shibata@riken.jp) Experimental Data The following datasets (EXP-01–EXP-21) correspond to the experimental results presented in the main text and supplementary figures. EXP-01_spiral — Multi-armed spiral formation without image processing (Fig. 1b,c) EXP-02_band_velocity_distribution.csv — Spatial distribution of velocity of traveling bands (Fig. 1d) EXP-03_band_speed_vs_distance.csv — Dependence of speed of the traveling bands on the distance from the center (Fig. 1e) EXP-04_density_profile_x.csv — Cell density profile along the x axis (Fig. S10b) EXP-05_density_profile_y.csv — Cell density profile along the y axis (Fig. S10c) EXP-06_radial_PIV_vectorfield.txt — Vector field obtained by PIV along the radial direction (Fig. 1f) EXP-07_core_PIV_vectorfield.txt — Vector field within rotating core obtained by PIV (Fig. S1a) EXP-09_polar_order_timeseries.csv — Time series of the polar order parameter (Fig. S1b) EXP-10_band_shape_theta_vs_ra.csv — Shape analysis of the traveling band (θ vs r/a) (Fig. 2b) EXP-11_band_shape_theta_vs_r.csv — Shape analysis of the traveling band (θ vs r) (Fig. S1c) EXP-12_density_intensity_map.csv — Intensity map of the cell density (Fig. S2b) EXP-13_a_and_VOmega.csv — Time evolution of the slope parameter a by fitting and V/Ω (Fig. 2c) EXP-14_band_speed_and_angular_speed.csv — Time evolution of the band speed V and the angular speed Ω (Fig. 2e) EXP-15_cell_distance_upon_collision.csv — Distance between two cells before and after collision (Fig. 3c) EXP-16_collision_angle_in_out.csv — Incoming and outgoing angle upon two-cell collision (Fig. 3d) EXP-17_cluster_size_distribution.csv — Distribution of cell cluster size in the cell chirality assay (Fig. S3a) EXP-18_angular_change_sample1.csv — Distribution of angular change (sample 1) (Fig. 3e) EXP-19_angular_change_sample2.csv — Distribution of angular change (sample 2) (Fig. S3b) EXP-20_angular_change_sample3.csv — Distribution of angular change (sample 3) (Fig. S3b) EXP-21_angular_velocity_3samples.csv — Angular velocity of the KI cells (3 samples) (Fig. 3f) --- Simulation Data The following datasets (SIM-01–SIM-27) correspond to the simulation results presented in the main text and supplementary figures. SIM-01_spiral_particle_coordinates.dat — Coordinates of each particle: Multi-armed spiral (Fig. 4a) SIM-02_local_phi_t_16.025.dat — Spatial distribution of local polar order parameter (t=16.025) (Fig. S4a) SIM-03_local_phi_t_26.71.dat — Spatial distribution of local polar order parameter (t=26.71) (Fig. S4a) SIM-04_local_phi_t_37.4.dat — Spatial distribution of local polar order parameter (t=37.4) (Fig. S4a) SIM-05_local_phi_t_48.075.dat — Spatial distribution of local polar order parameter (t=48.075) (Fig. S4a) SIM-06_concentric_particle_coordinates.dat — Coordinates of each particle: Concentric pattern (Fig. 4a) SIM-07_spot_particle_coordinates.dat — Coordinates of each particle: Spot (Fig. 4a) SIM-08_density_and_orientation.csv — Density profile and cell orientation along the radial direction (Fig. 4a) SIM-09_phase_diagram.csv — Phase diagram (Fig. 4b) SIM-10_band_shape_theta_vs_ra.csv — Temporal dynamics of the spiral shape (constant supply rate) (Fig. 4c) SIM-11_a_vs_time.csv — Time evolution of the slope parameter a (constant supply rate) (Fig. 4c inset) SIM-12_band_shape_constant_supply_theta_vs_ra.csv — Shape analysis of the traveling band (constant supply rate, θ vs r/a) (Fig. 4d) SIM-13_band_shape_constant_supply_theta_vs_r.csv — Shape analysis of the traveling band (constant supply rate, θ vs r) (Fig. S4d) SIM-14_a_vs_supply_rate.csv — Slope parameter as a function of supply rate (Fig. 4d inset) SIM-15_band_shape_decaying_supply_theta_vs_ra.csv — Temporal dynamics of the spiral shape (decaying supply rate, θ vs r/a) (Fig. 4e) SIM-16_band_shape_decaying_supply_theta_vs_r.csv — Temporal dynamics of the spiral shape (decaying supply rate, θ vs r) (Fig. S4e) SIM-17_density_intensity_map.npz — Intensity map of the cell density (Fig. S2e) SIM-18_a_and_VOmega.csv — Time evolution of the slope parameter a by fitting and V/Ω (decaying supply rate) (Fig. 2c) SIM-19_band_speed_and_angular_speed.csv — Time evolution of the band speed V and the angular speed Ω (decaying supply rate) (Fig. 2e) SIM-20_polar_order_density_gradient.csv — Dependence of polar order direction on the density gradient (Fig. S8) SIM-21_reflection_boundary.dat — Coordinates of each particle (reflection boundary) (Fig. S4b) SIM-22_KVmodel.dat — Coordinates of each particle (KVmodel) (Fig. S6a) SIM-23_band_shape_theta_vs_r.csv — Temporal dynamics of the spiral shape (constant supply rate) (Fig. S4c) SIM-24_self-alignment-based model.dat — Coordinates of each particle (Self-alignment-based model) (Fig. S6c) SIM-25_collision_angle_in_out.csv — Incoming and outgoing angle upon collision (Fig. S6b) SIM-26_chiv_on_omega.csv — negative chiral circulation (Fig. S5a) SIM-27_cq_on_omega.csv — spiral chirality measure (Fig. S5b) --- Simulation Codes The following codes (CODE-01–CODE-05) were used to perform the simulations presented in the manuscript. CODE-01_ConstSupply.f90 — Main model code for constant supply CODE-02_DecaySupply.f90 — Main model code for decaying supply CODE-03_VicsekKuramoto.f90 — Code for Vicsek Kuramoto model CODE-04_SelfAlignment.f90 — Code for STERIC repulsion and Self Alignment model CODE-05_DecayReflectBoundary.f90 — Code for reflecting boundary with the main model