Retinal stem cells modulate proliferative parameters to coordinate post-embryonic morphogenesis in the eye of fish.

Combining clonal analysis with a computational agent based model, we investigate how tissue-specific stem cells for neural retina (NR) and retinal pigmented epithelium (RPE) of the teleost medaka (Oryzias latipes) coordinate their growth rates. NR cell division timing is less variable, consistent with an upstream role as growth inducer. RPE cells divide with greater variability, consistent with a downstream role responding to inductive signals. Strikingly, the arrangement of the retinal ciliary marginal zone niche results in a spatially biased random lineage loss, where stem- and progenitor cell domains emerge spontaneously. Further, our data indicate that NR cells orient division axes to regulate organ shape and retinal topology. We highlight an unappreciated mechanism for growth coordination, where one tissue integrates cues to synchronize growth of nearby tissues. This strategy may enable evolution to modulate cell proliferation parameters in one tissue to adapt whole-organ morphogenesis in a complex vertebrate organ.

A 3D self-organizing multicellular epidermis model of barrier formation and hydration with realistic cell morphology based on EPISIM.

The epidermis and the stratum corneum (SC) as its outermost layer have evolved to protect the body from evaporative water loss to the environment. To morphologically represent the extremely flattened cells of the SC - and thereby the epidermal barrier - in a multicellular computational model, we developed a 3D biomechanical model (BM) based on ellipsoid cell shapes. We integrated the BM in the multicellular modelling and simulation platform EPISIM. We created a cell behavioural model (CBM) with EPISIM encompassing regulatory feedback loops between the epidermal barrier, water loss to the environment, and water and calcium flow within the tissue. This CBM allows a small number of stem cells to initiate self-organizing epidermal stratification, yielding the spontaneous emergence of water and calcium gradients comparable to experimental data. We find that the 3D in silico epidermis attains homeostasis most quickly at high ambient humidity, and once in homeostasis the epidermal barrier robustly buffers changes in humidity. Our model yields an in silico epidermis with a previously unattained realistic morphology, whose cell neighbour topology is validated with experimental data obtained from in vivo images. This work paves the way to computationally investigate how an impaired SC barrier precipitates disease.

Bridging the scales: semantic integration of quantitative SBML in graphical multi-cellular models and simulations with EPISIM and COPASI.

MOTIVATION: Biological reality can in silico only be comprehensively represented in multi-scaled models. To this end, cell behavioural models addressing the multi-cellular level have to be semantically linked with mechanistic molecular models. These requirements have to be met by flexible software workflows solving the issues of different time scales, inter-model variable referencing and flexible sub-model embedding. RESULTS: We developed a novel software workflow (EPISIM) for the semantic integration of Systems Biology Markup Language (SBML)-based quantitative models in multi-scaled tissue models and simulations. This workflow allows to import and access SBML-based models. SBML model species, reactions and parameters are semantically integrated in cell behavioural models (CBM) represented by graphical process diagrams. By this, cellular states like proliferation and differentiation can be flexibly linked to gene-regulatory or biochemical reaction networks. For a multi-scale agent-based tissue simulation executable code is automatically generated where different time scales of imported SBML models and CBM have been mapped. We demonstrate the capabilities of the novel software workflow by integrating Tyson's cell cycle model in our model of human epidermal tissue homeostasis. Finally, we show the semantic interplay of the different biological scales during tissue simulation. AVAILABILITY: The EPISIM platform is available as binary executables for Windows, Linux and Mac OS X at http://www.tiga.uni-hd.de. Supplementary data are available at http://www.tiga.uni-hd.de/supplements/SemSBMLIntegration.html. CONTACT: niels.grabe@bioquant.uni-heidelberg.de.

Wound healing revised: a novel reepithelialization mechanism revealed by in vitro and in silico models.

Wound healing is a complex process in which a tissue's individual cells have to be orchestrated in an efficient and robust way. We integrated multiplex protein analysis, immunohistochemical analysis, and whole-slide imaging into a novel medium-throughput platform for quantitatively capturing proliferation, differentiation, and migration in large numbers of organotypic skin cultures comprising epidermis and dermis. Using fluorescent time-lag staining, we were able to infer source and final destination of keratinocytes in the healing epidermis. This resulted in a novel extending shield reepithelialization mechanism, which we confirmed by computational multicellular modeling and perturbation of tongue extension. This work provides a consistent experimental and theoretical model for epidermal wound closure in 3D, negating the previously proposed concepts of epidermal tongue extension and highlighting the so far underestimated role of the surrounding tissue. Based on our findings, epidermal wound closure is a process in which cell behavior is orchestrated by a higher level of tissue control that 2D monolayer assays are not able to capture.

Modeling multi-cellular behavior in epidermal tissue homeostasis via finite state machines in multi-agent systems.

MOTIVATION: For the efficient application of multi-agent systems to spatial and functional modeling of tissues flexible and intuitive modeling tools are needed, which allow the graphical specification of cellular behavior in a tissue context without presuming specialized programming skills. RESULTS: We developed a graphical modeling system for multi-agent based simulation of tissue homeostasis. An editor allows the intuitive and hierarchically structured specification of cellular behavior. The models are then automatically compiled into highly efficient source code and dynamically linked to an interactive graphical simulation environment. The system allows the quantitative analysis of the morphological and functional tissue properties emerging from the cell behavioral model. We demonstrate the relevance of the approach using a recently published model of epidermal homeostasis as well as a series of cell-cycle models. AVAILABILITY: The complete software is available in binary executables for MS-Windows and Linux at tiga.uni-hd.de.