{'fraction_non_alphanumeric': 0.03409257039267858, 'fraction_numerical': 0.013040554369755271, 'mean_word_length': 4.9193492557978535, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 0, 'https://': 1, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Impact of changing cell-cell communication network in models of epithelial pattern formation

When modelling multi-cellular systems, one has to account for cell-cell signalling in addition to the molecular networks driving cell behaviours. Here, we aim at exploring how the topology of the cell-cell communication network impacts the behaviour of the whole multicellular system. More precisely, we focus on epithelial pattern formation, on which our question can be rephrased in terms of cell sizes and shapes. Relying on a logical modelling framework, and using a simple lateral inhibition model over a population of epithelial cells, we assess the model behaviours considering a variety of communication networks. This study suggests that reasonable deviations from a fixed grid (with regular hexagonal shaped cells) do not change much the resulting patterns. We further explore the impact of cell shapes and show that characteristics such as network regularity and number of shared neighbours of contacting cells are relevant to qualify such deviations.


INTRODUCTION
Logical modelling (Abou-Jaoudé et al., 2016) has shown to be useful in handling regulatory networks involved in pattern formation in the course of tissue development (Albert and Othmer, 2003;Sánchez et al., 2008;Azpeitia et al., 2010;Fauré et al., 2014). This qualitative framework allows to consider large regulatory networks for which detailed quantitative data are lacking. Logical models account for essential dynamical properties such as model attractors, which can be generally associated with different cell types displaying distinct patterns of gene expression. We have recently released EpiLog, a software tool implementing a hierarchical extension of the logical framework to account for cell-cell communication within a simple epithelium (Varela et al., 2018b). EpiLog defines a cellular automaton over a 2D hexagonal grid, i.e., a grid of hexagons), in which each cell carries its own logical regulatory model. Specific mechanisms govern the observed regular organisation of epithelia (e.g. Vetter et al. (2019)). In EpiLog, the underlying cell communication network is a fixed hexagonal grid, which convey the predominant hexagonal shape of the cells (Gibson et al., 2006). Here, we aim at exploring the implications of EpiLog communication network. To this end, we consider alternative networks diverging from that of the hexagonal grid. To uncover the impact of the topology of the communication network on the induced pattern, we selected a simple model of lateral inhibition, known to drive a precise salt-and-pepper patterning (Collier et al., 1996).
For convenience, and because of our case study, this introduction is restricted to Boolean models, in contrast to multi-valued models in which variables can take more discrete values. This paper is organised as follows. The background Section 2 provides a brief overview of the modelling framework for cellular and multi-cellular systems, as well as an introduction to the lateral inhibition model. Section 3 introduces the selected communication networks and associated characteristics considered in this study. Section 4 presents the methods used to generate the different communication networks and to simulate the lateral-inhibition model over those networks, getting the reached stable patterns. Section 5 presents the results, which are discussed together with some prospects in Section 6. Logical modelling (Abou-Jaoudé et al., 2016) has shown to be useful in handling regulatory networks involved in pattern formation in the course of tissue development (Albert and Othmer, 2003;Sánchez et al., 2008;Azpeitia et al., 2010;Fauré et al., 2014). This qualitative framework allows to consider large regulatory networks for which detailed quantitative data are lacking. Logical models account for essential dynamical properties such as model attractors, which can be generally associated with different cell types displaying distinct patterns of gene expression. We have recently released EpiLog, a software tool implementing a hierarchical extension of the logical framework to account for cell-cell communication within a simple epithelium (Varela et al., 2018b). EpiLog defines a cellular automaton over a 2D hexagonal grid, i.e., a grid of hexagons), in which each cell carries its own logical regulatory model.
Specific mechanisms govern the observed regular organisation of epithelia (e.g. Vetter et al. (2019)). In EpiLog, the underlying cell communication network is a fixed hexagonal grid, which convey the predominant hexagonal shape of the cells (Gibson et al., 2006 that of the hexagonal grid. To uncover the impact of the topology of the communication network on the induced pattern, we selected a simple model of lateral inhibition, known to drive a precise salt-and-pepper patterning (Collier et al., 1996).
For convenience, and because of our case study, this introduction is restricted to Boolean models, in contrast to multi-valued models in which variables can take more discrete values. This paper is organised as follows. The background Section 2 provides a brief overview of the modelling framework for cellular and multi-cellular systems, as well as an introduction to the lateral inhibition model. Section 3 introduces the selected communication networks and associated characteristics considered in this study. Section 4 presents the methods used to generate the different communication networks and to simulate the lateral-inhibition model over those networks, getting the reached stable patterns. Section 5 presents the results, which are discussed together with some prospects in Section 6.  -Jaoudé et al., 2016) has shown to be useful in handling regulatory networks involved in pattern formation in the course of tissue development (Albert and Othmer, 2003;Sánchez et al., 2008;Azpeitia et al., 2010;Fauré et al., 2014). This qualitative framework allows to consider large regulatory networks for which detailed quantitative data are lacking. Logical models account for essential dynamical properties such as model attractors, which can be generally associated with different cell types displaying distinct patterns of gene expression. We have recently released EpiLog, a software tool implementing a hierarchical extension of the logical framework to account for cell-cell communication within a simple epithelium (Varela et al., 2018b). EpiLog defines a cellular automaton over a 2D hexagonal grid, i.e., a grid of hexagons), in which each cell carries its own logical regulatory model.

Boolean modelling
Specific mechanisms govern the observed regular organisation of epithelia (e.g. Vetter et al. (2019)). In EpiLog, the underlying cell communication network is a fixed hexagonal grid, which convey the predominant hexagonal shape of the cells (Gibson et al., 2006). Here, we aim at exploring the implications of EpiLog communication network. To this end, we consider alternative networks diverging from that of the hexagonal grid. To uncover the impact of the topology of the communication network on the induced pattern, we selected a simple model of lateral inhibition, known to drive a precise salt-and-pepper patterning (Collier et al., 1996).
For convenience, and because of our case study, this introduction is restricted to Boolean models, in contrast to multi-valued models in which variables can take more discrete values. This paper is organised as follows. The background Section 2 provides a brief overview of the modelling framework for cellular and multi-cellular systems, as well as an introduction to the lateral inhibition model. Section 3 introduces the selected communication networks and associated characteristics considered in this study. Section 4 presents the methods used to generate the different communication networks and to simulate the lateral-inhibition model over those networks, getting the reached stable patterns. Section 5 presents the results, which are discussed together with some prospects in Section 6.

Impact of changing cell-cell communication network in models of epithelial pattern formation
Abstract: When modelling multi-cellular systems, one has to account for cell-cell signalling in addition to the molecular networks driving cell behaviours. Here, we aim at exploring how the topology of the cell-cell communication network impacts the behaviour of the whole multicellular system. More precisely, we focus on epithelial pattern formation, on which our question can be rephrased in terms of cell sizes and shapes. Relying on a logical modelling framework, and using a simple lateral inhibition model over a population of epithelial cells, we assess the model behaviours considering a variety of communication networks. This study suggests that reasonable deviations from a fixed grid (with regular hexagonal shaped cells) do not change much the resulting patterns. We further explore the impact of cell shapes and show that characteristics such as network regularity and number of shared neighbours of contacting cells are relevant to qualify such deviations.

INTRODUCTION
Logical modelling (Abou-Jaoudé et al., 2016) has shown to be useful in handling regulatory networks involved in pattern formation in the course of tissue development (Albert and Othmer, 2003;Sánchez et al., 2008;Azpeitia et al., 2010;Fauré et al., 2014). This qualitative framework allows to consider large regulatory networks for which detailed quantitative data are lacking. Logical models account for essential dynamical properties such as model attractors, which can be generally associated with different cell types displaying distinct patterns of gene expression. We have recently released EpiLog, a software tool implementing a hierarchical extension of the logical framework to account for cell-cell communication within a simple epithelium (Varela et al., 2018b). EpiLog defines a cellular automaton over a 2D hexagonal grid, i.e., a grid of hexagons), in which each cell carries its own logical regulatory model.
Specific mechanisms govern the observed regular organisation of epithelia (e.g. Vetter et al. (2019)). In EpiLog, the underlying cell communication network is a fixed hexagonal grid, which convey the predominant hexagonal shape of the cells (Gibson et al., 2006). Here, we aim at exploring the implications of EpiLog communication network. To this end, we consider alternative networks diverging from that of the hexagonal grid. To uncover the impact of the topology of the communication network on the induced pattern, we selected a simple model of lateral inhibition, known to drive a precise salt-and-pepper patterning (Collier et al., 1996).
For convenience, and because of our case study, this introduction is restricted to Boolean models, in contrast to multi-valued models in which variables can take more discrete values. This paper is organised as follows. The background Section 2 provides a brief overview of the modelling framework for cellular and multi-cellular systems, as well as an introduction to the lateral inhibition model. Section 3 introduces the selected communication networks and associated characteristics considered in this study. Section 4 presents the methods used to generate the different communication networks and to simulate the lateral-inhibition model over those networks, getting the reached stable patterns. Section 5 presents the results, which are discussed together with some prospects in Section 6. Abstract: When modelling multi-cellular systems, one has to account for cell-cell signalling in addition to the molecular networks driving cell behaviours. Here, we aim at exploring how the topology of the cell-cell communication network impacts the behaviour of the whole multicellular system. More precisely, we focus on epithelial pattern formation, on which our question can be rephrased in terms of cell sizes and shapes. Relying on a logical modelling framework, and using a simple lateral inhibition model over a population of epithelial cells, we assess the model behaviours considering a variety of communication networks. This study suggests that reasonable deviations from a fixed grid (with regular hexagonal shaped cells) do not change much the resulting patterns. We further explore the impact of cell shapes and show that characteristics such as network regularity and number of shared neighbours of contacting cells are relevant to qualify such deviations.

INTRODUCTION
Logical modelling (Abou-Jaoudé et al., 2016) has shown to be useful in handling regulatory networks involved in pattern formation in the course of tissue development (Albert and Othmer, 2003;Sánchez et al., 2008;Azpeitia et al., 2010;Fauré et al., 2014). This qualitative framework allows to consider large regulatory networks for which detailed quantitative data are lacking. Logical models account for essential dynamical properties such as model attractors, which can be generally associated with different cell types displaying distinct patterns of gene expression. We have recently released EpiLog, a software tool implementing a hierarchical extension of the logical framework to account for cell-cell communication within a simple epithelium (Varela et al., 2018b). EpiLog defines a cellular automaton over a 2D hexagonal grid, i.e., a grid of hexagons), in which each cell carries its own logical regulatory model.
Specific mechanisms govern the observed regular organisation of epithelia (e.g. Vetter et al. (2019)). In EpiLog, the underlying cell communication network is a fixed hexagonal grid, which convey the predominant hexagonal shape of the cells (Gibson et al., 2006). Here, we aim at exploring the implications of EpiLog communication network. To this end, we consider alternative networks diverging from that of the hexagonal grid. To uncover the impact of the topology of the communication network on the induced pattern, we selected a simple model of lateral inhibition, known to drive a precise salt-and-pepper patterning (Collier et al., 1996).
For convenience, and because of our case study, this introduction is restricted to Boolean models, in contrast to multi-valued models in which variables can take more discrete values. This paper is organised as follows. The background Section 2 provides a brief overview of the modelling framework for cellular and multi-cellular systems, as well as an introduction to the lateral inhibition model. Section 3 introduces the selected communication networks and associated characteristics considered in this study. Section 4 presents the methods used to generate the different communication networks and to simulate the lateral-inhibition model over those networks, getting the reached stable patterns. Section 5 presents the results, which are discussed together with some prospects in Section 6. ) has shown to be useful in handling regulatory networks involved in pattern formation in the course of tissue development (Albert and Othmer, 2003;Sánchez et al., 2008;Azpeitia et al., 2010;Fauré et al., 2014). This qualitative framework allows to consider large regulatory networks for which detailed quantitative data are lacking. Logical models account for essential dynamical properties such as model attractors, which can be generally associated with different cell types displaying distinct patterns of gene expression. We have recently released EpiLog, a software tool implementing a hierarchical extension of the logical framework to account for cell-cell communication within a simple epithelium (Varela et al., 2018b). EpiLog defines a cellular automaton over a 2D hexagonal grid, i.e., a grid of hexagons), in which each cell carries its own logical regulatory model.
Specific mechanisms govern the observed regular organisation of epithelia (e.g. Vetter et al. (2019)). In EpiLog, the underlying cell communication network is a fixed hexagonal grid, which convey the predominant hexagonal shape of the cells (Gibson et al., 2006). Here, we aim at exploring the implications of EpiLog communication network. To this end, we consider alternative networks diverging from that of the hexagonal grid. To uncover the impact of the topology of the communication network on the induced pattern, we selected a simple model of lateral inhibition, known to drive a precise salt-and-pepper patterning (Collier et al., 1996).
For convenience, and because of our case study, this introduction is restricted to Boolean models, in contrast to multi-valued models in which variables can take more discrete values. This paper is organised as follows. The background Section 2 provides a brief overview of the modelling framework for cellular and multi-cellular systems, as well as an introduction to the lateral inhibition model. Section 3 introduces the selected communication networks and associated characteristics considered in this study. Section 4 presents the methods used to generate the different communication networks and to simulate the lateral-inhibition model over those networks, getting the reached stable patterns. Section 5 presents the results, which are discussed together with some prospects in Section 6.

Boolean modelling
A Boolean model of a regulatory network defines a set of discrete variables x taking their values in {0, 1}, where the value of x conveys the state of the corresponding regulatory component g x (active or not, expressed or not, etc.). A Boolean function f x defines the evolution of the variable x, i.e., the value of x depending on the state of the model (more precisely on the states of the components regulating g x ). Indeed, the function f x reflects the regulatory structure of g x (its regulators and their effects on g x ).
A transition function is thus defined: given any state X = (x 1 , . . . , x n ) of the model, F (X) = (f x1 , . . . , f xn ) defines its target state. In a deterministic, synchronous update, F (X) is the successor of X. However, as this updating scheme is often biologically unrealistic, it is preferable to consider the non-deterministic asynchronous update for which at most one variable is updated at each simulation step, each successor differing from X by exactly one component.
A cellular model relates to the regulatory network controlling a specific cellular process. Input components accounting for signals from the micro-environment have no regulation and are associated to constant variables (i.e., f x (X) = x).

Boolean modelling of multi-cellular systems & EpiLog
To handle multi-cellular systems in which cell-cell communication matters, we proposed to compose cellular models through their input components (see Chaouiya et al. (2013); Varela et al. (2018a) for further detail on this composition). Basically, cell-cell communication operates through signalling proteins that activate receptors of neighbouring cells. For each cell, the neighbouring cells are those connected to the cell through an edge in the underlying communication network.
A multi-cellular model is defined by a collection of as many cellular models as the number of cells, the communication network defining the neighbours of each cell (i.e., cells enabled to send signals to that cell) and, for each input component, an integration function that defines which signals and from how many neighbouring cells (expressed by a cardinality constraint of the type AT Lc, i.e., at least c) should send the signal to activate that input. Here, we focus on juxtacrine signals, i.e., signalling between contacting cells, i.e., that are connected in the communication network.
EpiLog (http://epilog-tool.org) is a software tool implementing this framework in the form of cellular automata restricted to 2D hexagonal grids (Varela et al., 2018b). Hexagonal grids have cylindrical or toroidal periodic conditions or not. Hence the sole communication network handled in EpiLog is a regular grid with degree 6 (except possibly at the grid borders). The purpose of this work is to assess to what extent this topological restriction is acceptable.
The default updating scheme in EpiLog is the synchronous one, where all the cellular components in all the cells are simultaneously updated. However, EpiLog allows to associate (synchronous) priority classes with the cellular components and to specify an α-asynchronous update over the cells (Fatès, 2013). This consists in overcoming the classical synchronous update of cellular automata by defining α, the common updating probability of the cells. The resulting update is synchronous when α = 1 and, by convention, when α = 0 the update is equivalent to the random asynchronous (a unique cell is updated at each simulation step, selecting a single trajectory among all the possible asynchronous behaviours) (Varela et al., 2018b). Finally within a simulation step, when updating a cell, its input nodes are updated first, according to the values of components in neighbouring cells and only then their targets are updated.

Lateral inhibition model
A simplified version of the Delta-Notch signalling system was chosen for the purpose of this study (Collier et al., 1996). This system is well-known for establishing distinct fates of contacting cells through a lateral inhibition mechanism: Notch signalling is initiated when Delta ligand binds to the Notch receptor expressed on a contacting cell, causing its cleavage and the release of the Notch intracellular domain. This then moves into the nucleus where it represses Delta production. The Boolean model of this system is illustrated in Figure 1.
This model is convenient to assess the influence of the communication network topology because it gives rise to archetypal patterns that can be characterised by a simple measure of the cell-type ratio, here assessed as the percentage of Notch expressing cells (see e.g. Matsuda et al. (2015)). Hereinafter, cells in which Notch is active (resp. inactive) and Delta is inactive (resp. active) are referred to as N (resp. D) positive cells.
Diverse mathematical models have been delineated to study patterns generated by the Delta-Notch signalling system (see e.g. Collier et al. (1996); Ghosh and Tomlin (2004); Palau-Ortin et al. (2015); Tonello and Siebert (2019)). In particular, in Tonello and Siebert (2019), prove that for their Boolean model, under the random asynchronous update, all the attractors are stable states (i.e., there are no stable oscillatory behaviours, and the model will eventually reach a stable pattern), and that all the attractors are reachable from any homogeneous initial configuration in which all the cells have the same state.

COMMUNICATION NETWORKS
This section presents a selection of network classes, from regular grids to variations of regular networks and non- regular ones. We then introduce a few network characteristics considered along this study.

Network classes
In what follows, we denote by n, the number of nodes (or cells) in the communication network.
Regular grids (RG) are defined by cells with a unique shape (hexagonal or square) placed on a 2D torus. They correspond to regular networks subject to a geometrical constraint (i.e., 2D lattices), with degree k = 4 (square cells with a von Neumann neighbourhood), k = 6 (hexagonal cells), or k = 8 (square cells with a Moore neighbourhood). See Figure 2 for an illustration. Note that, being a 2D torus and due to the hexagonal structure (k = 6), dimensions of the grid must be even.
Random regular (RR) networks are such that all their n nodes have the same specified degree k (Steger and Wormald, 1999). Note that n × k must be even and k < n.
Watts-Strogatz (WG) networks are obtained from a regular network with nodes distributed over a ring, each being connected to its k nearest neighbours and then by rewiring the edges given a probability p (Watts and Strogatz, 1998;Barrat and Weigt, 2000). Varying p, nodes in the obtained networks -so-called "small-world" -have a mean degreek, and at least k/2 neighbours. Networks obtained with a probability p are henceforth denoted WGp.
Erdös-Rényi (ER) networks are also associated with a parameter p specifying the probability for a pair of nodes to be connected, the mean degree beingk = (n − 1) × p.

Network measures and characteristics
Notation: • A network (V, E) is defined by a set of nodes V , and a set of edges Network regularity: among the considered network classes, we can distinguish between those which are regular, where all the nodes have the exact same degree, and those which are non-regular, with node degrees following a given distribution. Figure 3-A shows the degree distribution of the three non-regular networks considered in this work.
Network connectedness: all networks (RG, RR, and WG) but the Erdös-Rényi are connected.
Number of shared neighbours: Given two adjacent nodes v and w, S(v, w) denotes the number of neighbours shared by v and w: Czekanowski-Dice dissimilarity measure: This classical measure applies to pairs of nodes v and w and accounts for the numbers of shared neighbours and the numbers of nodes exclusively connected to v or to w (Angelelli et al., 2008). Here, we only consider the Czekanowski-Dice dissimilarity measure for pairs of connected nodes v and w ((v, w) ∈ E), given by: Hence, the higher the number of shared neighbours, the lower the dissimilarity. In particular, if two nodes have the same neighbours, their dissimilarity measure is zero. In the context of cell communication networks, this corresponds to a situation where two cells are subject to exactly the same signals from their (shared) neighbours, and will thus adopt the same fate.
Note that both shared neighbours and dissimilarity measures are relevant and can be compared for all connected networks, i.e., all except ER (see Table 1). For regular networks, given a degree k, the number of connected pairs of nodes is exactly k × n/2, whereas for non-regular networks this is the mean number of connected nodes.
All pairs of connected nodes in the RG network with k = 6 have exactly 2 shared neighbours. As shown in Figure 3   Because the WS network with p = 0.13 (denoted WS0.13) has 2 shared neighbours on average, we retained this specific network for comparisons in Section 5.

METHODS
The generation of the Random regular (RR), Watts-Strogatz (WS) and Erdös-Rényi (ER) networks was performed using the NetworkX 2.2 (https://networkx. github.io/) Python software package (Hagberg et al., 2008). The Regular grid (RG) was generated using our own Python script. Section 7 provides the location of all the scripts of network generation and analyses, as well as for simulations and pattern analyses.
For each network class and each parameter combination (see Table 2), we generated 100 communication networks. Note that for the RG, only one network was generated for each degree, which defines a unique network. Network statistics, such as degree distribution or number of shared neighbours, are retrieved immediately after each network generation using a Python script.
To carry out all the simulations, we used a modified version of EpiLog that permits to load a predefined communi-cation network. For each communication network (100 instances for all classes but the RG), 100 simulations were performed. For each instance of the RG, 10 000 simulations were done. All simulations were carried out under a synchronous update within cells, a random asynchronous update (α = 0) between cells, and started from the naive initial state, in which all the components are set to value 0. This is because all the stable patterns are reachable from this homogeneous initial state (Tonello and Siebert, 2019).
EpiLog creates a file for each stable pattern reached. Then a Python script retrieves the percentage of N positive cells from these files.

RESULTS
Considering the different classes of communication networks described above, we evaluated how the stable patterns of the lateral-inhibition model change depending on different parameters. Patterns are characterised by the percentages of N positive cells.

Network regularity and cardinality constraints
In this section, the (mean) degree of all the nodes in all the networks is fixed to 6. The cardinality constraint varies from AT L1 to AT L6.
Increased cardinality constraint means that more neighbouring D positive cells (i.e., connected in the communication network) are required to get a N positive cell (i.e.. As illustrated in the Figure 4-A for the regular hexagonal grid, the percentage of N positive cells (green) in the reached stable pattern decreases with the increase of the cardinality constraint. This decrease is exactly compensated by the emergence of D positive cells. Figure 4-B shows that for all the communications networks, varying the cardinality constraint (AT Lc for c ∈ {1, . . . , 6}) results in a decrease of the average percentages of N positive cells in the resulting stable patterns.
Moreover, as illustrated in Figure 5, for all the regular networks, the percentages of N positive cells for the cardinality constrains AT Lc and AT L6−c+1 (for c ∈ {1, 2, 3}) add  up to 100%. In other words, the percentage of N positive cells for AT Lc equals the percentage of D positive cells for AT L6 − c + 1. This property could be easily proved for any regular communication network. However, this property does not hold for the non-regular networks. Indeed, cells with less than c neighbours, eventually become D positive since the cardinality constraint AT Lc cannot be satisfied.

Network degree
Here, we consider a fixed cardinality constraint AT L1, but the result would be qualitatively similar for any cardinality constraint.
Within each network class, Figure 6 shows the impact of the degree of the nodes on the resulting patterns. Clearly, whatever the network class, the higher the (mean) node degree, the higher the percentage of N positive cells. This is easily explained by the fact that a higher degree allows a single D positive cell to activate N in more neighbouring cells.
Furthermore, for non-regular networks, the wider the degree distribution (see Figure 3- the percentages of N positive cells differ within the class of regular networks. We propose in the next section that the number of shared neighbours may explain these differences.

Number of shared neighbours
Recall that, given a pair of connected nodes, their number of shared neighbours relates to their dissimilarity measure.
Here, we consider the average number of shared neighbours, accounting for all pairs of connected nodes.
Although limited differences, a higher number of shared neighbours seems to lead to an increase of the percentages of N positive cells (see Figure 7). This can be informally explained as follows: two connected cells are either both N positive or of different types (Notch and Delta); in the first case, their shared neighbours can also be N positive (provided the presence of a D positive cell in the neighbourhood), whereas in the later case their shared neighbours are necessarily N positive cells.
From Figure 7, within the regular networks, it appears that when k = 4 (resp. k = 6) the percentage of N positive cells  grows with the number of shared neighbours, i.e., WS0.0 leads to more N positive cells than WS0.13, which leads to more N positive cells than RG (resp. WS0.0 leads to more N positive cells than RG, which leads to (slightly) more N positive cells than WS0.13). The same property holds for the (non-regular) Watts-Strogatz networks WS0.13 and WS1.0, where the former leads to a higher percentage of N positive cells and has a higher average number of shared neighbours (for bothk = 4 andk = 6).

DISCUSSION AND CONCLUSION
Our goal was to assess the impact of the communication network underlying the EpiLog grid. To do so, we considered various classes of regular and non-regular networks and, focusing on the Notch signalling system, we evaluated the obtained stable patterns in terms of percentages of N positive cells. Overall, these percentages do not change drastically, suggesting that a 2D hexagonal grid (RG, k = 6) is a reasonable approximation for the modelling of simple epithelial tissues. However, we showed that the network topology does have a measurable impact, which depends on regularity, degree and average number of shared neighbours.
From Figure 5, it appears that the closest percentage of N positive cells is obtained with the WS0.13 (k = 6) network. This network has a degree distribution quite similar to the one observed in nature (Gibson et al., 2006), and could be taken as a reference.
Network regularity matters in keeping a complementarity property between cardinality constraints, where percentages of N positive cells for complementary cardinality constraints, e.g. AT L1 and AT L6, add up to 100% (see Figure 5). This property is not maintained with nonregular communication networks, however the WS0.13 is the one which almost keep this complementarity.
Moreover, the numbers of shared neighbours seem to impact the resulting patterns (Figure 7). Here again, the WS0.13 network is the network whose number of shared neighbours is the closest to that of the RG.
The patterns reached with the Boolean Delta-Notch model used for this study were qualified by the percentage of N positive cells. This measure does not allow to assess the nature of the pattern in terms of local organisation, regularity and symmetry. In Palau-Ortin et al. (2015), the authors propose parameters that allow to do so. However, these parameters are defined for the hexagonal grid, and their extension to general communication networks might be difficult. Nonetheless, it would be interesting to have a measure that better qualify the organisation of the patterns resulting from the lateral-inhibition model for different communication networks.
Finally, we relied for this study on a simple lateralinhibition model. The observations made with this model could be generalised for other patterning systems.

AVAILABILITY
Python scripts used to generate the communication networks for each network class, to extract the network characteristics, to perform the simulations with EpiLog, and to extract the Notch percentage from the obtained patterns is available at http://epilog-tool.org/node/190.