During morphogenesis, cells precisely generate forces that drive cell shape changes and cellular motion. These forces predominantly arise from contractility of the actomyosin cortex, allowing for cortical tension, protrusion formation, and cell division. Image-based force inference can derive such forces from microscopy images, without complicated and time-consuming experimental set-ups. However, current methods do not account for common effects, such as physical confinement and local force generation. Here we propose a force-inference method based on a biophysical model of cell shape, and assess relative cellular surface tension, adhesive tension between cells, as well as cytokinesis and protrusion formation. We applied our method on fluorescent microscopy images of the early C. elegans embryo. Predictions for cell surface tension at the 7-cell stage were validated by measurements using cortical laser ablation. Our non-invasive method facilitates the accurate tracking of force generation, and offers many new perspectives for studying morphogenesis.

MOTIVATION: Uncovering the cellular and mechanical processes that drive embryo formation requires an accurate read out of cell geometries over time. However, automated extraction of 3D cell shapes from time-lapse microscopy remains challenging, especially when only membranes are labeled. RESULTS: We present an image analysis framework for automated tracking and three-dimensional cell segmentation in confocal time lapses. A sphere clustering approach allows for local thresholding and application of logical rules to facilitate tracking and unseeded segmentation of variable cell shapes. Next, the segmentation is refined by a discrete element method simulation where cell shapes are constrained by a biomechanical cell shape model. We apply the framework on Caenorhabditis elegans embryos in various stages of early development and analyze the geometry of the 7- and 8-cell stage embryo, looking at volume, contact area and shape over time. AVAILABILITY AND IMPLEMENTATION: The Python code for the algorithm and for measuring performance, along with all data needed to recreate the results is freely available at 10.5281/zenodo.5108416 and 10.5281/zenodo.4540092. The most recent version of the software is maintained at https://bitbucket.org/pgmsembryogenesis/sdt-pics. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Actin protrusion dynamics plays an important role in the regulation of three-dimensional (3D) cell migration. Cells form protrusions that adhere to the surrounding extracellular matrix (ECM), mechanically probe the ECM and contract in order to displace the cell body. This results in cell migration that can be directed by the mechanical anisotropy of the ECM. However, the subcellular processes that regulate protrusion dynamics in 3D cell migration are difficult to investigate experimentally and therefore not well understood. Here, we present a computational model of cell migration through a degradable viscoelastic ECM. This model is a 2D representation of 3D cell migration. The cell is modeled as an active deformable object that captures the viscoelastic behavior of the actin cortex and the subcellular processes underlying 3D cell migration. The ECM is regarded as a viscoelastic material, with or without anisotropy due to fibrillar strain stiffening, and modeled by means of the meshless Lagrangian smoothed particle hydrodynamics (SPH) method. ECM degradation is captured by local fluidization of the material and permits cell migration through the ECM. We demonstrate that changes in ECM stiffness and cell strength affect cell migration and are accompanied by changes in number, lifetime and length of protrusions. Interestingly, directly changing the total protrusion number or the average lifetime or length of protrusions does not affect cell migration. A stochastic variability in protrusion lifetime proves to be enough to explain differences in cell migration velocity. Force-dependent adhesion disassembly does not result in faster migration, but can make migration more efficient. We also demonstrate that when a number of simultaneous protrusions is enforced, the optimal number of simultaneous protrusions is one or two, depending on ECM anisotropy. Together, the model provides non-trivial new insights in the role of protrusions in 3D cell migration and can be a valuable contribution to increase the understanding of 3D cell migration mechanics.

Perfusion bioreactors regulate flow conditions in order to provide cells with oxygen, nutrients and flow-associated mechanical stimuli. Locally, these flow conditions can vary depending on the scaffold geometry, cellular confluency and amount of extra cellular matrix deposition. In this study, a novel application of the immersed boundary method was introduced in order to represent a detailed deformable cell attached to a 3D scaffold inside a perfusion bioreactor and exposed to microscopic flow. The immersed boundary model permits the prediction of mechanical effects of the local flow conditions on the cell. Incorporating stiffness values measured with atomic force microscopy and micro-flow boundary conditions obtained from computational fluid dynamics simulations on the entire scaffold, we compared cell deformation, cortical tension, normal and shear pressure between different cell shapes and locations. We observed a large effect of the precise cell location on the local shear stress and we predicted flow-induced cortical tensions in the order of 5 pN/μm, at the lower end of the range reported in literature. The proposed method provides an interesting tool to study perfusion bioreactors processes down to the level of the individual cell's micro-environment, which can further aid in the achievement of robust bioprocess control for regenerative medicine applications.

The EU-funded eNanoMapper project proposes a computational infrastructure for toxicological data management of engineered nanomaterials (ENMs) based on open standards, ontologies and an interoperable design to enable a more effective, integrated approach to European research in nanotechnology. eNanoMapper's goal is to support the collaborative safety assessment for ENMs by creating a modular, extensible infrastructure for transparent data sharing, data analysis, and the creation of computational toxicology models for ENMs. The eNanoMapper database solution builds on previous experience of the consortium partners in supporting diverse data through flexible data storage, semantic web technologies, open source components and web services. A number of opportunities and challenges exist in nanomaterials representation and integration of ENM information, originating from diverse systems. A short summary, highlighting the pros and cons of the existing integration approaches and data models is presented. We demonstrate the approach of adopting an ontology-supported data model, describing the materials and measurements. The data sources supported include diverse formats (ISA-Tab, OECD Harmonized Templates, custom spreadsheet templates, various databases provided by consortia members). Besides retaining the data provenance, the focus on measurements provides insights into how to reuse the chemical structure database tools for nanomaterials characterization and safety.

Adhesion governs to a large extent the mechanical interaction between a cell and its microenvironment. As initial cell spreading is purely adhesion driven, understanding this phenomenon leads to profound insight in both cell adhesion and cell-substrate interaction. It has been found that across a wide variety of cell types, initial spreading behavior universally follows the same power laws. The simplest cell type providing this scaling of the radius of the spreading area with time are modified red blood cells (RBCs), whose elastic responses are well characterized. Using a mechanistic description of the contact interaction between a cell and its substrate in combination with a deformable RBC model, we are now able to investigate in detail the mechanisms behind this universal power law. The presented model suggests that the initial slope of the spreading curve with time results from a purely geometrical effect facilitated mainly by dissipation upon contact. Later on, the spreading rate decreases due to increasing tension and dissipation in the cell's cortex as the cell spreads more and more. To reproduce this observed initial spreading, no irreversible deformations are required. Since the model created in this effort is extensible to more complex cell types and can cope with arbitrarily shaped, smooth mechanical microenvironments of the cells, it can be useful for a wide range of investigations where forces at the cell boundary play a decisive role.