Keywords
Logical modelling, Multicellular regulatory networks, Cellular automaton, Hexagonal grid
Logical modelling, Multicellular regulatory networks, Cellular automaton, Hexagonal grid
Pattern formation emerges from the interplay of interaction networks, at the cellular and multi-cellular levels1. To uncover the complex mechanisms at stake, computational modelling is very needed. In this context, cellular automaton approaches are particularly well suited2.
On the other hand, the logical formalism has proved efficient to explore cellular regulatory networks driving developmental processes (for a recent review on the logical modelling approach, see 3). The consideration of ensembles of communicating cells is often required to recapitulate observed patterns (e.g. 4,5).
This motivated the development of EpiLog, which implements an extension of the logical framework to hexagonal grids embodying simple epithelia.
Briefly, a logical model assigns discrete (Boolean or multi-valued) variables to the regulatory components, and logical regulatory rules define the evolution of the variables depending on the variables associated with the component regulators. Input components embody cell receptors receiving external signals. Unlike internal components, inputs have no associated regulatory rules, and are generally maintained at constant values accounting for fixed environmental conditions. A model state is given by the values of the components, and the dynamics is generated according to specific updating schemes. In a synchronous update, a state has at most one successor state in which all the variables have been updated as prescribed by the regulatory rules. In contrast, the asynchronous update generates as many successors as the number of updated variables, yielding more complex and non-deterministic behaviours. Properties of the resulting dynamics relate to the model attractors (stable states or point attractors, and cyclical attractors), which can be associated with different phenotypes. Beyond identifying the attractors, a major question relates to reachability properties: given an initial condition, what is (are) the reachable attractor(s)?
The definition, analysis and simulation of logical models of cellular networks can be performed using one of the existing software tools6,7. Many of these support the exchange format SBML qual (Qualitative Models package for SBML)7,8.
To handle multi-cellular systems, Mendes et al. have established a compositional approach relying on a process algebra framework to assess stable states reachability in asynchronous dynamics9. This approach showed some limitations due to the huge size of the considered state spaces, and is thus limited to the analysis of models encompassing a reduced number of cells (up to a dozen cells). Nonetheless, in the cellular automaton framework considered here, we use the same approach to specify cell-cell communication through (logical) integration rules. These govern the evolution of cellular input components, depending on the values of internal components (representing e.g. secreted proteins) in neighbouring cells.
EpiLog is implemented in Java (requires Java 7 or higher), is freely available at http://epilog-tool.org, and launches a Graphical User Interface (GUI) based on JFC/Swing. It is provided as a single .jar file which can be launched from a graphical file manager by double clicking on it or through the command line as follows:
$ java -jar EpiLog-versionX.jar [--file /path/to/multicellular_project.peps]
where the flag --file can be used to optionally specify the file to be open. EpiLog .peps files represent projects containing multiple multicellular models, each with potentially distinct characteristics.
Additionally, an EpiLog model repository is provided at http://epilog-tool.org/models_repository where users can deposit their models by contacting support@epilog-tool.org. Each model page includes a title, a description, a taxon and process classification, model .peps file(s) and supporting paper(s).
Developers can freely access the code repository at http://github.com/epilog-tool/epilog, clone it and extend the code. For dependency management of external libraries, EpiLog relies on Apache Maven. A particular library is the bioLQM Java toolkit (Logical Qualitative Models of biological regulatory networks), for the representation and manipulation of logical, cellular models, available at https://github.com/colomoto/bioLQM10.
Model definition EpiLog defines projects that include models, called epithelia, each being specified by:
1. The size of the hexagonal grid, as well as its border conditions; vertical and horizontal borders can be connected, or not leading to rectangular, cylindrical or toroidal grids;
2. A (set of) logical model(s), each associated with cells of the grid; these cellular models can be defined in any tool supporting the SBML qual format8. EpiLog imports SBML files and identifies input components in these cellular models;
3. The cellular input components; these are defined as constant positional inputs (e.g. morphogens produced by external sources) or evolving integration inputs; for the latter, logical integration rules specify which signals are received from (internal) components of neighbouring cells and how these signals are combined;
4. An initial state defining the values of the internal components of all the cells.
Simulation settings For each epithelium, simulation settings define the updating schemes for both the cellular and epithelium models. The default scheme is synchronous, that is to say, at each iteration, all the cellular models are considered for a synchronous update of their internal components. To overcome spurious oscillations generated by these synchronous updates (see Figure 1A and Supplementary File 1), EpiLog allows to specify:
1. The cellular model update, with the definition of (synchronous) priority classes (increasing or decreasing) updates are gathered in ordered sets, allowing to account for different time scales for the mechanisms underlying component updates11;
2. The epithelium update, with the definition of probabilistic updates of the cell states, inspired by N. Fatès’ alpha-asynchronism12: parameter α specifies the proportion of the cells randomly chosen for update.
(A) An idealised Boolean model of lateral inhibition; the Green marker is induced by the Red marker secreted by the neighbouring cell, whereas Red is induced by low levels of Green in the same cell; the synchronous dynamics of this two cell model leads to a cyclic attractor (see Supplementary File 1). (B) Simulation in EpiLog, with the same cellular model over a square grid of 50 × 50 hexagonal cells under a synchronous update (Green induced if at least one contacting cell is Red), oscillations are due to the synchronous update of the cells, see Supplementary Movie 1. (C–F) Simulations under α-asynchrony (α = 0.25): (C) Stable pattern reached in 29 steps (Green induced if at least one contacting cell is Red), see Supplementary Movie 2; (D) Stable reverted pattern reached in 30 steps (Green induced if all contacting cells are Red); (E) Stable pattern reached in 63 steps (Green induced if at least 12 cells at distance up to 3 are Red); (F) Stable pattern reached in 70 steps, with the same setting as for panel E but considering a torus (i.e., no grid borders).
Model perturbations The logical formalism permits to easily specify perturbations (e.g. knock-out or ectopic expression of cellular model components) by blocking the values of perturbed components. EpiLog includes this feature with the definition of component perturbations that can be applied to all the cells or to restricted regions of the grid.
The reader is invited to consult the documentation available at EpiLog web site to get further information (http://epilog-tool.org/documentation).
As a first illustration, we consider a very simple model, in which the cellular model includes two markers. When considering a juxtacrine signal (i.e., only contacting cells can communicate), this model accounts for the well-known Delta-Notch signalling, which implements a lateral inhibition process (see e.g. 13). The model is provided in Supplementary Model File 1. In Figure 1A, a two-cell model is shown, to illustrate interactions within and between the cells. Remaining panels of Figure 1 show EpiLog simulation results, illustrating the different patterns obtained when varying the updating mode, the integration rules and the border conditions.
The development of EpiLog was originally motivated by our second illustration, which relates to the dorsal appendage-forming regions of the Drosophila Melanogaster eggshell14 (see Figure 2). The logical model defined by Fauré et al. reproduces the patterning of the anterior follicular epithelium of the oocyte, defining the floor (Rho expressing cells) and roof (Br expressing cells) of the future appendages5. The model is provided in Supplementary Model File 2. EpiLog window is shown in Figure 2A, with the simulation panel displaying the stable pattern obtained with this model. The remaining panels C–G show different states of the grid, with panels F–G recapitulating patterns resulting from model perturbations (see Figures 6–8 in 5).
(A) EpiLog main window with the stable pattern reached by a simulation starting with all the cells of the grid having their internal components at 0, and positional inputs defined as shown in panel C (phase 1). (B) Egg chamber with two dorsal respiratory appendages (DA); DA primordia are established as 2 regions on both sides of the oocyte midline, with follicle cells expressing Broad (Br in red), future roof of the DA, and cells expressing Rhomboid (Rho in blue), future floor of the DA. Establishment of these regions involve Gurken signalling (Grk, in orange) and Decapentaplegic signalling (Dpp, in violet). (C) Grk and Dpp gradient defined as positional inputs in EpiLog. (D) Stable pattern obtained with a simulation starting from the pattern displayed in panel A and without the Grk signal (phase 2), suggesting that Grk extinction is required to split the floor regions (see 5). (E) EpiLog simulation of a mild overexpression of Dpp. (F) EpiLog simulation of Pointed (Pnt, internal component) loss-of-function. (G) EpiLog simulation of Pnt gain-of-function clones. Note that in panels F–G, perturbed cellular models are indicated by bold borders in the grid.
To the best of our knowledge, EpiLog is the first software tool for the definition, simulation and visualisation of qualitative, logical (Boolean and multivalued) models over hexagonal grids. It provides a graphical user interface and tools to conveniently support the study of epithelial pattern formation, relying on a logical framework.
EpiLog is free, open source and implemented in Java for operating system independence. It relies on the SBML qual standard format to import cellular models.
EpiLog v1.1 is available from: http://epilog-tool.org/downloads
Source code available from: https://github.com/epilog-tool/epilog
Archived source code as at time of publication: https://doi.org/10.5281/zenodo.132050315
License: GNU General Public License v3.0
PV and CR developed the software. PTM acquired funding, developed the software and supervised the project. CC designed the project, acquired funding and supervised the project. PLV, CVR, PTM and CC wrote the article.
This work was supported by national funds through Fundação para Ciência e a Tecnologia (FCT) with reference PTDC/BEX-BCB/0772/2014, UID/CEC/50021/2013 and IF/01333/2013. PLV was supported by grant PTDC/BEX-BCB/0772/2014 and PTDC/EIA-CCO/099229/2008. CVR was supported by PTDC/BEX-BCB/0772/2014 and IF/01333/2013/CP1204/CT0001. Furthermore, CC and CVR acknowledge support from the Fundação Calouste Gulbenkian.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors would like to thank Adrien Fauré who implemented a first Python prototype that inspired the development of EpiLog, and who provided insightful feedback. We are also grateful to Céline Hernandez for her useful comments.
Supplementary File 1: Dynamics of the lateral inhibition model under synchronous and asynchronous updates.
Click here to access the data.
Supplementary Movie 1: Oscillatory behaviour of the lateral inhibition model under a synchronous update (see Figure 1 panel B).
Click here to access the data.
Supplementary Movie 2: Stable pattern of the lateral inhibition model under an α-asynchronous update (see Figure 1 panel C).
Click here to access the data.
Supplementary Model File 1: EpiLog file including variations of the lateral inhibition model as illustrated in Figure 1.
Click here to access the data.
Supplementary Model File 2: EpiLog file including variations of the DA model as illustrated in Figure 2.
Views | Downloads | |
---|---|---|
F1000Research | - | - |
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Is the rationale for developing the new software tool clearly explained?
Yes
Is the description of the software tool technically sound?
Yes
Are sufficient details of the code, methods and analysis (if applicable) provided to allow replication of the software development and its use by others?
Yes
Is sufficient information provided to allow interpretation of the expected output datasets and any results generated using the tool?
Partly
Are the conclusions about the tool and its performance adequately supported by the findings presented in the article?
Yes
References
1. Mendes ND, Lang F, Le Cornec YS, Mateescu R, et al.: Composition and abstraction of logical regulatory modules: application to multicellular systems.Bioinformatics. 2013; 29 (6): 749-57 PubMed Abstract | Publisher Full TextCompeting Interests: No competing interests were disclosed.
Reviewer Expertise: Systems Biology, Cancer, Boolean modelling, Multiscale modelling, Metabolic modelling, Data deconvolution
Is the rationale for developing the new software tool clearly explained?
Yes
Is the description of the software tool technically sound?
Yes
Are sufficient details of the code, methods and analysis (if applicable) provided to allow replication of the software development and its use by others?
Yes
Is sufficient information provided to allow interpretation of the expected output datasets and any results generated using the tool?
Partly
Are the conclusions about the tool and its performance adequately supported by the findings presented in the article?
Yes
References
1. Mendes ND, Lang F, Le Cornec YS, Mateescu R, et al.: Composition and abstraction of logical regulatory modules: application to multicellular systems.Bioinformatics. 2013; 29 (6): 749-57 PubMed Abstract | Publisher Full TextCompeting Interests: No competing interests were disclosed.
Reviewer Expertise: David Henriques: dynamic modelling, bioinformatics, logic-based modelling, systems biology. Alejandro F Villaverde: dynamic modelling, system identification, systems biology, control theory.
Is the rationale for developing the new software tool clearly explained?
Yes
Is the description of the software tool technically sound?
Yes
Are sufficient details of the code, methods and analysis (if applicable) provided to allow replication of the software development and its use by others?
Yes
Is sufficient information provided to allow interpretation of the expected output datasets and any results generated using the tool?
Yes
Are the conclusions about the tool and its performance adequately supported by the findings presented in the article?
Yes
Competing Interests: No competing interests were disclosed.
Reviewer Expertise: Artificial Life, Complex Systems, Evo-Devo, Cellular Automata
Alongside their report, reviewers assign a status to the article:
Invited Reviewers | |||
---|---|---|---|
1 | 2 | 3 | |
Version 2 (revision) 11 Mar 19 |
read | ||
Version 1 27 Jul 18 |
read | read | read |
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)