ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Software Tool Article
Revised

EpiLog: A software for the logical modelling of epithelial dynamics

[version 2; peer review: 3 approved]
PUBLISHED 11 Mar 2019
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

Cellular responses are governed by regulatory networks subject to external signals from surrounding cells and to other micro-environmental cues. The logical (Boolean or multi-valued)  framework proved well suited to study such processes at the cellular level, by specifying qualitative models of involved signalling pathways and gene regulatory networks. 

Here, we describe and illustrate the main features of EpiLog, a computational tool that implements an extension of the logical framework to the tissue level. EpiLog defines a collection of hexagonal cells over a 2D grid, which embodies a mono-layer epithelium. Basically, it defines a cellular automaton in which cell behaviours are driven by associated logical models subject to external signals. 

EpiLog is freely available on the web at http://epilog-tool.org. It is implemented in Java (version ≥1.7 required) and the source code is provided at https://github.com/epilog-tool/epilog under a GNU General Public License v3.0.

Keywords

Logical modelling, Multicellular regulatory networks, Cellular automaton, Hexagonal grid

Revised Amendments from Version 1

Following constructive comments from the reviewers, this revised version includes clarifications, which mostly concern:

  1. Related work on a compositional approach for logical models as well as related software tools supporting the simulation of populations of Boolean networks.
  2. Further discussion and illustration of the alpha-asynchronous updating, with an additional Supplementary File 2.
  3. Mention of the possibility to clone an epithelium model, a feature used for the model of the patterning of the fly eggshell.
'CNRS, Centrale Marseille, l'Institut de Mathématiques de Marseille, Aix-Marseille University, Marseille, France' has been added as an affiliation for Claudine Chaouiya.

See the authors' detailed response to the review by Alejandro F. Villaverde and David Henriques
See the authors' detailed response to the review by Stefano Nichele
See the authors' detailed response to the review by Arnau Montagud

Introduction

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 a natural choice of modelling framework, particularly well suited to define local rules governing the behaviours of cells over a lattice2. 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). However, 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)8,7.

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). In contrast to the previous compositional approach, EpiLog relies on a cellular automata framework, which allows to scale model simulations up to thousands of cells. Still, cell-cell interactions are specified using (logical) integration rules, in a way similar to that proposed by Mendes et al.9. These rules govern the evolution of cellular input components, depending on the values of internal components (representing e.g. secreted proteins) in neighbouring cells.

Multiple simulation platforms have been developed to support the modelling of multi-cellular systems. A detailed overview of existing approaches and application domains is beyond the scope of this paper. To our knowledge, only two tools exist that integrate Boolean models in a multi-cellular context. The first, named CoGNaC10, is a Chaste plugin. Chaste is an open-source simulation package that includes tissue modelling. As such, it supports distinct cell-based modelling approaches11,12. The second tool, named PhysiBoSS13, combines intra-cellular Boolean signalling networks as handled by MaBoSS14 and multi-cellular agent based modelling as defined in PhysiCell, an open source agent-based simulator15. While these tools involve explicit representations of intra-cellular networks in the form of Boolean models, and the consideration of biophysical processes governing multi-cellular behaviours, they do not allow explicit specification of cell-cell communication, a special feature of EpiLog.

Methods

Implementation

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/bioLQM16.

Operation

Model definition EpiLog defines projects that include multi-cellular models, called epithelium models, 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 (e.g. juxtacrine signals, or other secreted molecules influencing cells in the neighbourhood); 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 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:

73b19a73-f4d4-4f6b-a278-a7b758d9c466_figure1.gif

Figure 1.

(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 (normal arrows denote activations, whereas blunt arrows denote inhibitions); 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. (CF) Simulations under α-asynchrony (α = 0.25): (C) Stable state of the grid, referred to as “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).

  • 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 updates of internal component17;

  • 2. The epithelium update, with the definition of probabilistic updates of the cell states, inspired by N. Fatès’ alpha-asynchronism18,19: parameter α specifies the proportion of the cells randomly chosen for update (see Supplementary File 2 for an illustration).

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.

Model clones At any simulation step, a clone of the epithelium model can be defined, having the current state of the grid as initial condition.

The reader is invited to consult the documentation available at EpiLog web site to get further information (http://epilog-tool.org/documentation).

Use cases

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. 20). 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 eggshell21 (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. EpiLog window is shown in Figure 2A, with the simulation panel displaying the stable pattern obtained with this model. Interestingly, the model shows that to get split regions (as displayed in Figure 2C), Gurken signal emitted by the nucleus, and defined as positional input in the model (Figure 2B), must be removed. This is simulated in a second phase, which is performed by cloning the model when reaching the stable pattern (Figure 2A), and by setting the Gurken signal to zero all over the grid. The remaining panels E-G show different states of the grid, recapitulating patterns resulting from model perturbations. For further detail on this modelling study, the reader is referred to 5.

73b19a73-f4d4-4f6b-a278-a7b758d9c466_figure2.gif

Figure 2. Model of the fly eggshell patterning of Fauré et al.5, implemented in EpiLog.

(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 (in graded orange, with 3 different levels) and Dpp (in purple) 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 FG, perturbed cellular models are indicated by bold borders in the grid.

Summary

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.

Software availability

EpiLog v1.1.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.132050322

License: GNU General Public License v3.0

Author information

PV and CVR 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.

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 27 Jul 2018
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Varela PL, Ramos CV, Monteiro PT and Chaouiya C. EpiLog: A software for the logical modelling of epithelial dynamics [version 2; peer review: 3 approved]. F1000Research 2019, 7:1145 (https://doi.org/10.12688/f1000research.15613.2)
NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 2
VERSION 2
PUBLISHED 11 Mar 2019
Revised
Views
11
Cite
Reviewer Report 18 Mar 2019
Arnau Montagud, Institut Curie, Inserm U900, Mines Paris Tech, PSL Research University, Paris, France;  Barcelona Supercomputing Center (BSC), Barcelona, Spain 
Approved
VIEWS 11
The authors have have done a good work in expanding and ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Montagud A. Reviewer Report For: EpiLog: A software for the logical modelling of epithelial dynamics [version 2; peer review: 3 approved]. F1000Research 2019, 7:1145 (https://doi.org/10.5256/f1000research.20225.r45550)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Version 1
VERSION 1
PUBLISHED 27 Jul 2018
Views
28
Cite
Reviewer Report 10 Oct 2018
Arnau Montagud, Institut Curie, Inserm U900, Mines Paris Tech, PSL Research University, Paris, France;  Barcelona Supercomputing Center (BSC), Barcelona, Spain 
Approved with Reservations
VIEWS 28
The paper describes a software used for qualitative simulations of logical models in a hexagonal grid, called EpiLog. The paper is well written, even though some parts could be slightly expanded for a better comprehension. The software is easy to ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Montagud A. Reviewer Report For: EpiLog: A software for the logical modelling of epithelial dynamics [version 2; peer review: 3 approved]. F1000Research 2019, 7:1145 (https://doi.org/10.5256/f1000research.17032.r37601)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 11 Mar 2019
    Claudine Chaouiya, Instituto Gulbenkian de Ciência, Oeiras, P-2780-156, Portugal
    11 Mar 2019
    Author Response
    We thank you for the detailed comments and suggestions. We have revised the manuscript accordingly, and we briefly answer your points below.

    We have revised the introduction, noting that ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 11 Mar 2019
    Claudine Chaouiya, Instituto Gulbenkian de Ciência, Oeiras, P-2780-156, Portugal
    11 Mar 2019
    Author Response
    We thank you for the detailed comments and suggestions. We have revised the manuscript accordingly, and we briefly answer your points below.

    We have revised the introduction, noting that ... Continue reading
Views
18
Cite
Reviewer Report 04 Oct 2018
Alejandro F. Villaverde, Bioprocess Engineering Group,  IIM-CSIC, Vigo, Spain 
David Henriques, Bioprocess Engineering Group, IIM-CSIC, Vigo, Spain 
Approved
VIEWS 18
This paper describes a new software tool for modelling epithelial dynamics in 2D. The software, called EpiLog, uses a logical framework that assigns discrete values to the model variables. The approach used for specifying the logical rules was previously presented ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Villaverde AF and Henriques D. Reviewer Report For: EpiLog: A software for the logical modelling of epithelial dynamics [version 2; peer review: 3 approved]. F1000Research 2019, 7:1145 (https://doi.org/10.5256/f1000research.17032.r38246)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 11 Mar 2019
    Claudine Chaouiya, Instituto Gulbenkian de Ciência, Oeiras, P-2780-156, Portugal
    11 Mar 2019
    Author Response
    We thank you for the comments and suggestions, and we briefly answer each point or the report below.

    1) We believe the misunderstanding was due to the confusion with ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 11 Mar 2019
    Claudine Chaouiya, Instituto Gulbenkian de Ciência, Oeiras, P-2780-156, Portugal
    11 Mar 2019
    Author Response
    We thank you for the comments and suggestions, and we briefly answer each point or the report below.

    1) We believe the misunderstanding was due to the confusion with ... Continue reading
Views
35
Cite
Reviewer Report 20 Aug 2018
Stefano Nichele, Department of Computer Science, Oslo Metropolitan University, Oslo, Norway 
Approved
VIEWS 35
The article describes EpiLog, a software tool to model artificial development of epitelial tissue through cellular automata with hexagonal cells. 

EpiLog provides a high degree of customization through the graphical user interface. This includes: (a) model definition, ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Nichele S. Reviewer Report For: EpiLog: A software for the logical modelling of epithelial dynamics [version 2; peer review: 3 approved]. F1000Research 2019, 7:1145 (https://doi.org/10.5256/f1000research.17032.r37297)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 11 Mar 2019
    Claudine Chaouiya, Instituto Gulbenkian de Ciência, Oeiras, P-2780-156, Portugal
    11 Mar 2019
    Author Response
    We than you for the suggestions.  We considered that a proper state of the art was beyond the scope of this application note. However, we now acknowledge the existence of ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 11 Mar 2019
    Claudine Chaouiya, Instituto Gulbenkian de Ciência, Oeiras, P-2780-156, Portugal
    11 Mar 2019
    Author Response
    We than you for the suggestions.  We considered that a proper state of the art was beyond the scope of this application note. However, we now acknowledge the existence of ... Continue reading

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 27 Jul 2018
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

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.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.