Supplementary Data for "A numerical algorithm for modeling cellular rearrangements in tissue morphogenesis" by R.Z. Mohammad, H. Murakawa, K. Svadlenka, and H. Togashi
Supplementary Data 1. Dataset used to generate Figures 3b and 3f: Plot of relative intensity of β-catenin accumulations at OO, OS, and SS junctions during cellular rearrangements of WT and αN-catenin KO olfactory epithelia, obtained by interpolation of β-catenin intensity mean values at time instants E14, E16, E18 and P1.
Supplementary Data 2. Dataset used to generate Figures 3d, 3h, and 4c: Box plots of cell neighbor counts obtained from actual images WT and αN-catenin KO olfactory epithelia at time instants E14, E16, and P1; from their corresponding simulation images; and from WT simulation with hypothetical P1 tension values.
Supplementary Data 3. Dataset used to generate Figure 4a: Plot of evolution of interfacial tensions based on measurements up to embryonic stage E18 but is selected hypothetically for postnatal stage P1.
Supplementary Data 4. Dataset used to generate Figure 5e: Plot of evolution of interfacial tensions at SS, HS, and HH junctions used in WT and Nectin-3 KO simulations of an embryonic auditory epithelium.
Supplementary Data 5. Dataset used to generate Supplementary Figure 4b and 4c: Log-log plot of relative error of mean radius of growing circles and shrinking circles at t = 0.01875 for varying mesh-time configurations.
Supplementary Data 6. Dataset used to generate Supplementary Figure 5d and 5e: Log-log plot of relative error in radius of the best fitted circle to each interface and relative error in contact angles at triple junction for varying mesh-time configurations.
Supplementary Data 7. Dataset used to generate Supplementary Figure 6c: Plot of evolution of the magnitude of contact angle measured across topology change for varying mesh-time configurations.
Supplementary Data 8. Dataset used to generate Supplementary Figures 8d-f and 9: Plot of evolution of OO-junction length, total and surface energy, and error in cell volumes for vertex dynamics model with volume penalties ρ = 500 and ρ =100; and level-set based method.
Supplementary Data 8. Dataset used to generate Supplementary Figures 11: Log-log plots of CPU times for executing the vertex dynamics algorithm on test case A (with topology change) and test case B (without topology change), showing increased cost when T1-swap is present.
Funding
JSPS Grants-in-Aid for Scientific Research (KAKENHI) [Grant Numbers 19K03634, 19H04965, 18H04764, 18H05481, 18K06219, 18H01139, 20H01823]
Grant-in-aid for JSPS Fellows [Grant Number 18F18016]
PRESTO (Precursory Research for Embryonic Science and Technology) Program [Grant Number JPMJPR1946]