Model ID: https://identifiers.org/morpheus/M4377 Software: Morpheus (open source). Download from: https://morpheus.gitlab.io Full title: Neuromast Regeneration in Zebrafish Authors: N. G. Lavalle, J. Miranda-Rodríguez, E. C. Costa, A. Borges, O. Viader-Llargués, H. López-Schier, O. Chara Submitter: N. G. Lavalle Curator: D. Jahn Date: 04.03.2026 Reference: This model is the original used in the publication, up to technical updates: N. G. Lavalle, J. Miranda-Rodríguez, E. C. Costa, A. Borges, O. Viader-Llargués, H. López-Schier, O. Chara: Local control of cellular proliferation underlies neuromast regeneration in zebrafish. J. Theor. Biol. 623: 112389, 2026. https://doi.org/10.1016/j.jtbi.2026.112389 Comment: This Morpheus XML model implements a 2D simulation based on the Cellular Potts Model (CPM) to study neuromast regeneration in Danio rerio (zebrafish) larvae following severe ablation. The simulation is initialized from experimentally acquired images, which define the initial spatial configuration and reduced number of surviving cells after ablation. Starting from this image-based initial condition, the model reproduces organ regeneration through divisions of sustentacular and mantle cells. The core mechanism of the model is a local control of proliferation based on cell neighborhood: each cell capable of dividing (sustentacular or mantle) halts its cell cycle when it reaches a critical number of neighboring cells of the same type. Specifically, sustentacular cells stop proliferating when they contact 3 neighboring sustentacular cells, while mantle cells do so when surrounded by 1 neighboring mantle cell. This mechanism reproduces both the regenerative dynamics and the spatial architecture of the organ observed experimentally.

3

sum_cell_center_y/total_cells sum_cell_center_x/total_cells celltype.mantle.size+celltype.hair.size+celltype.sustentacular.size tau*dia < time and rand_uni(0,1) < pm and Mantle_Neighbours < umbral 1

p == 1 0 p == 2 2 1 p == 3 and s == 1 2 1 (tau*dia < time) and (rand_uni(0,1) < ps )and (Sustentacular_Neighbours < umbral) and (r<=a) if(d<0.594,0,1) (tau*dia < time) and (rand_uni(0,1) < ps ) and (Sustentacular_Neighbours < umbral) and (r>a) if(daughter == 1, 1, 0 ) if(d < 0.9,0,if((0.9<=d) and (d< 0.95),2,3)) sqrt((cell.center.x-mean_center_x)^2+(cell.center.y-mean_center_y)^2)

rand_uni(0,1) 2 2