Supplementary Data for "A numerical algorithm for modeling cellular rearrangements in tissue morphogenesis" by R.Z. Mohammad, H. Murakawa, K. Svadlenka, and H. Togashi

<b>Supplementary Data 1. </b>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.<br><b>Supplementary Data 2. </b>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.<br><b>Supplementary Data 3.</b> 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.<br><br><b>Supplementary Data 4.</b> 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.<br><br><b>Supplementary Data 5.</b> 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.<br><b>Supplementary Data 6.</b> 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.<br><b>Supplementary Data 7.</b> 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.<br><br><b>Supplementary Data 8.</b> 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.<br><br><b>Supplementary Data 8.</b> 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.<br><br><br>