3D cell neighbour dynamics in growing pseudostratified epithelia

Abstract
Introduction
Results
Discussion
Materials and methods
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

During morphogenesis, epithelial sheets remodel into complex geometries. How cells dynamically organise their contact with neighbouring cells in these tightly packed tissues is poorly understood. We have used light-sheet microscopy of growing mouse embryonic lung explants, three-dimensional cell segmentation, and physical theory to unravel the principles behind 3D cell organisation in growing pseudostratified epithelia. We find that cells have highly irregular 3D shapes and exhibit numerous neighbour intercalations along the apical-basal axis as well as over time. Despite the fluidic nature, the cell packing configurations follow fundamental relationships previously described for apical epithelial layers, that is, Euler's polyhedron formula, Lewis’ law, and Aboav-Weaire's law, at all times and across the entire tissue thickness. This arrangement minimises the lateral cell-cell surface energy for a given cross-sectional area variability, generated primarily by the distribution and movement of nuclei. We conclude that the complex 3D cell organisation in growing epithelia emerges from simple physical principles.

Introduction

Common to all animals and plants, epithelia are a fundamental tissue type whose expansion, budding, branching, and folding is key to the morphogenesis of organs and body cavities. Characterised by apical-basal polarity (Figure 1a), epithelial cells adhere tightly to their apical neighbours in a virtually impermeable adhesion belt, form lateral cell-cell junction complexes along the apico-basal axis to provide mechanical stabilisation, and bind tightly to the basal lamina and extracellular matrix (ECM) on the basal side (Drubin and Nelson, 1996; Rodriguez-Boulan and Macara, 2014; Shin and Margolis, 2006). How cell neighbour relationships are organised in these tightly adherent layers, and how these change during tissue and concomitant cell shape changes is poorly understood, despite their importance for cell-cell signalling and the fluidity of the tissue.

Figure 1

Download asset Open asset

Principles of epithelial organisation.

(a) Schematic representation of an epithelial tissue layer. The cells are polarised between an apical and a basal side. Near the apical side, cells adhere tightly via adhesion junctions (green). Nuclei are depicted in blue. (b) Apical surface projection of an embryonic lung bud at E12.5 imaged using light-sheet microscopy. Cell contour segmentations (left) coloured according to neighbour relationships (middle) and area quantifications (right). (c) Current shape representations of 3D epithelial cells: prism, frustum, prismatoid, and scutoid. (d) Planar cell neighbour exchange (T1 transition). (e) Tissues differ widely in the frequency of neighbour numbers. The legend provides the measured average number of cell neighbours for each tissue and the references to the primary data (Classen et al., 2005; Escudero et al., 2011; Etournay et al., 2015; Farhadifar et al., 2007; Gibson et al., 2006; Heller et al., 2016; Sánchez-Gutiérrez et al., 2016). Data points for n < 3 were removed as they must present segmentation artefacts. (f) The measured average number of cell neighbours is close to the topological requirement ( $\bar{n} = 6$ ) in all tissues; see panel e for the colour code. (g) Epithelial tissues follow the AW law (black line). The AW law formulates a relationship between the average number of neighbours, n, that a cell has and that its direct neighbours have, $m_{n}$ . The product $m_{n} ∙ n$ can be determined by summing over all *n_i*. (h) The relative average apical cell area, $\bar{A_{n}} / \bar{A}$ , increases with the number of neighbours, n, and mostly follows the linear Lewis’ law (Equation 2, black line), or the quadratic relationship (Equation 3, yellow line) in case of higher apical area variability. (i) The average internal angle by polygon type is close to that of a regular polygon, $θ_{n} = (n - 2) / n ∙ 180 °$ (yellow lines). To form a contiguous lattice, the angles at each tricellular junction must add to 360°, and the resulting observed deviation in the angles follows the prediction (red line). (j) The average normalised side length by polygon type. (k) Observed fraction of hexagons versus area coefficient of variation (CV). The curves mark theoretical predictions when polygonal cell layers follow either the linear Lewis’ law (Equation 2, black line) or the quadratic law (Equation 3, yellow line). The colour code in panels g-k is as in panel e, but data is available only for a subset of tissues. The abbreviations in panel (e) are as follows: cNT refers to the Chick neural tube epithelium, EYE to the *Drosophila* eye disc, dPE to the *Drosophila* peripodal membrane from the larval eye disc, dWL to the *Drosophila* larval wing disc, dWP to the *Drosophila* pre-pupal wing disc, dMWP to the *Drosophila* mutant wing pre-pupa with reduced expression of myosin II, dPW to the *Drosophila* pupal wing disc, and dMWL to the wing disc with gigas RNAi clones. TPx indicates subsequent but not further specified pupal time points.

© 2019, Kokic et al. Panels e, f, g and k are reproduced with modifications from Figures 1,5 from Kokic et al., 2019, which are published under the Creative Commons Attribution Non-Commercial 4.0 International License

© 2019, Vetter et al. Panels g and i are reproduced with modifications from Figures 1,5 from Vetter et al., 2019, which are published under the Creative Commons Attribution Non-Commercial 4.0 International License

Cell neighbour relationships can be most easily studied on epithelial surfaces, and the polygonal arrangements of apical surfaces (Figure 1b) have been meticulously analysed (Classen et al., 2005; Escudero et al., 2011; Etournay et al., 2015; Farhadifar et al., 2007; Gibson et al., 2006; Gómez-Gálvez et al., 2018; Heller et al., 2016; Kokic et al., 2019; Ramanathan et al., 2019; Sánchez-Gutiérrez et al., 2016). Widely considered to be a reliable proxy for three-dimensional (3D) cell shape, 3D epithelial cell shapes are often depicted as prisms with polygonal faces that retain the same neighbour relationships along the entire apico-basal axis (Figure 1c). Cells in curved epithelia are pictured as frustra, which have the same number of sides, but different apical and basal areas. If the curvature differs substantially along the principal axes, as is the case in epithelial tubes, neighbour relationships must change along the apical-basal axis. Prismatoids accommodate the neighbour change at the surface, while scutoids undergo the neighbour change somewhere along the apical-basal axis (Figure 1c,d; Gómez-Gálvez et al., 2018). However, even though the curvature is the same in both principal directions of spherically shaped epithelia, the neighbour relationships still differ between the apical and basal sides (Gómez-Gálvez et al., 2018), suggesting that effects other than curvature must determine the 3D neighbour arrangements of cells in epithelia.

Given the challenges in visualising 3D neighbour arrangements, most studies to date have focused on apical cell arrangements, and have revealed striking regularities. First, even though the frequencies of neighbour numbers differ widely between epithelial tissues (Figure 1e), cells have on average (close to) six neighbours (Classen et al., 2005; Escudero et al., 2011; Etournay et al., 2015; Farhadifar et al., 2007; Gibson et al., 2006; Heller et al., 2016; Kokic et al., 2019; Ramanathan et al., 2019; Sánchez-Gutiérrez et al., 2016; Figure 1f). This can be explained with topological constraints in contiguous polygonal lattices, as expressed by Euler's polyhedron formula (Gibson et al., 2006; Rivier and Lissowski, 1982). Thus, if three cells meet at each vertex, the average number of neighbours in infinitely large contiguous polygonal lattices is exactly

\bar{n} = 6.

While the average number of neighbours in the entire lattice is (close to) six, the local averages deviate from six, and instead rather closely follow a phenomenological relationship, termed Aboav-Weaire's law (Aboav, 1970). According to Aboav-Weaire's law (Figure 1g), the average number of neighbours of all n cells that border a cell with n neighbours follows as

m (n) = 5 + \frac{8}{n} .

Finally, the average apical area, $\bar{A_{n}}$ , of cells with n neighbours is linearly related to the number of cell neighbours, n (Figure 1h, black line), a relation termed Lewis’ law (Lewis, 1928),

\frac{\bar{A_{n}}}{\bar{A}} = \frac{(n - 2)}{4} .

Here, $\bar{A}$ refers to the average apical cell area in the tissue.

We have recently shown that Aboav-Weaire's law and Lewis’ law are a direct consequence of a minimisation of the lateral cell-cell contact surface energy (Kokic et al., 2019; Vetter et al., 2019). The lowest lateral cell-cell contact surface energy is obtained in a regular polygonal lattice because regular polygons have the smallest perimeter per polygonal area. The distribution of apical cell sizes that emerges from cell growth and division is, however, such that epithelial tissues cannot organise into perfectly regular polygonal lattices. By adhering to Aboav-Weaire's law and Lewis’ law, cells assume the most regular lattice. In particular, by following Aboav-Weaire's law, the internal angles are closest to that of a regular polygon, while adding up to 360° at each tricellular junction (Figure 1i; Vetter et al., 2019). And by following the relationship between polygon area and polygon type as stipulated by Lewis' law (Equation 2), the difference in side lengths, ${\bar{S}}_{n} / \bar{S}$ , is minimised between cells (Figure 1j; Kokic et al., 2019). The side lengths would be equal (Figure 1j, yellow line), if cells followed a quadratic relation of the form

\frac{A_{n}}{\bar{A}} = \frac{n}{6} \cdot \frac{\tan (\frac{π}{6})}{\tan (\frac{π}{n})} \sim {(\frac{n}{6})}^{2} .

This quadratic relation (Figure 1h, yellow line), however, requires a larger area variability than is observed in most epithelia imaged to date. Accordingly, the predicted quadratic relation had not been previously reported, but could be confirmed experimentally by us by increasing the apical area variability (Kokic et al., 2019).

Given the relationship between apical area and neighbour numbers as stipulated by Equations 1,3,4, the apical area variability emerges as the key determinant of apical epithelial organisation, and the theory correctly predicts how the fraction of hexagons in the tissue depends on the apical area variability, as can be quantified by the coefficient of variation (CV = std/mean) (Figure 1k; Kokic et al., 2019). As such, active processes such as growth, cell division, cell death and extrusion, cell intercalation and apical constriction determine the variability of the apical areas and thus determine apical organisation indirectly. Taken together, the apical organisation of epithelia can be understood based on the principles of lateral cell-cell contact surface energy minimisation.

In this work, we leverage these theoretical insights along with light-sheet fluorescence microscopy to study 3D epithelial organisation, both in cleared and growing pseudostratified epithelia. We find that cells have complex 3D shapes with numerous neighbour transitions along their apical-basal axis as well as over time. We show that much as on the apical side, the variation of the cross-sectional areas along the apical-basal axis defines the epithelial organisation at all times and across the entire tissue thickness. The observed neighbour arrangement minimises the lateral cell-cell surface energy for a given cross-sectional area variability. The cross-sectional areas vary as a result of active cell processes, most prominently including interkinetic nuclear migration (IKNM). We conclude that the complex 3D cell organisation in growing epithelia emerges from simple physical principles.

Results

Apical and basal epithelial organisation

We started by exploring the apical and basal cellular organisation in epithelial tubes and buds (Figure 2a). To this end, we imaged CUBIC-cleared mouse embryonic (E12.5) lung rudiments from a Shh^GC/+; ROSA^mT/mG background using light-sheet microscopy, and segmented the fluorescent membrane boundaries of over 400 cells per dataset in 2.5D (Figure 2b, Figure 2—figure supplement 2). The apical and basal surfaces are both curved and thus differ in their total areas, that is, the total segmented apical area is about 5-fold smaller than the basal area (Figure 2b). We detected less than half as many cells on the apical side, and the mean cross-sectional cell area of apical cells is therefore on average only twofold smaller than that of basal cells, while the area variability, measured as area CV, is higher (Figure 2c). Notably, the frequencies of the different neighbour numbers are not identical on the apical and basal side (Figure 2d), suggesting that the neighbour relationships change along the apical-basal axis, both in the tube and tip datasets. This observation is consistent with previous reports (Gómez-Gálvez et al., 2018; Ramanathan et al., 2019; Rupprecht et al., 2017). The change in neighbour relationships has previously been attributed to a curvature effect in tubes, but the neighbour changes in spherical geometries cannot be explained with curvature alone, since prismatic cells fully accommodate the cell dilation from equal curvature changes in both directions (Gómez-Gálvez et al., 2018; Rupprecht et al., 2017), suggesting that mainly other effects determine epithelial organisation.

Figure 2 with 4 supplements see all

Download asset Open asset

Apical and basal epithelial organisation.

(a) Epithelium of E12.5 Shh^GC/+; ROSA^mT/mG mouse embryonic lung imaged using light-sheet microscopy. Scale bar 200 µm. Corresponding 2D sections are shown in Figure 2—figure supplement 1, and Figure 2—video 1. (b) Apical and basal 2.5D cell segmentation overlays on imaged tip and tube sections (dotted boxes in panel a). An illustration of the 2.5D segmentation workflow is presented in Figure 2—figure supplement 2. Number of cells (N) and segmented surface areas (SA) are given. Cells are coloured using random labels. Scale bars 20 µm. (c) Normalised apical and basal cell area distributions in the tip (broken lines) and tube (solid lines) datasets. The colour code in panel c is reused in panels d-f. (d) Frequencies of neighbour numbers on the apical and basal sides in the tip and tube datasets. (e) The apical and basal layers follow the AW law (black line). SEM is smaller than symbols. (f) The normalised average cell area, $\frac{\bar{A_{n}}}{\bar{A}}$ , increases with the number of neighbours, n. The basal cells (green) follow Lewis’ law (Equation 2, black line), while the apical cells (purple) follow the quadratic relationship (Equation 3, yellow line). SEM as error bars.

Figure 2—source data 1 2.5D neighbour number and area data.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig2-data1-v2.xls
Download elife-68135-fig2-data1-v2.xls
Figure 2—source data 2 2.5D AW data.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig2-data2-v2.xls
Download elife-68135-fig2-data2-v2.xls

So how can we explain the difference in apical and basal epithelial organisation in both datasets? We have previously shown that the apical organisation emerges from the minimisation of the overall lateral cell-cell contact surface energy (Kokic et al., 2019; Vetter et al., 2019). Aboav-Weaire's law (Equation 2, Figure 1g) and the linear Lewis’ law or the quadratic relationship (Equations 3,4, Figure 1h) emerge as global organisation laws from this physical constraint, and ensure that the angles are closest to that of a regular polygon (Figure 1i, yellow lines), and that the side lengths are the most equal (Figure 1j, yellow line). We now find that these hold not only for the apical, but also for the basal datasets (Figure 2e,f). Consistent with our theory, the apical layers, which have a larger area variability than the basal layers (Figure 2c), follow the quadratic law (yellow line) rather than the linear Lewis’ law (black line).

We conclude that basal layers follow the same organisational principles as apical layers, such that their organisation can also be explained with a minimisation of the lateral cell-cell contact surface energy. Accordingly, the observed difference in overall neighbour relationships (Figure 2e,f) is a consequence of the difference in the cross-sectional area distributions (Figure 2c). So, why do the normalised area distributions differ between the apical and basal sides in both the tube segment and the bud, and how do they change along the apical-basal axis?

3D organisation of epithelia

To explore the physical principles behind 3D epithelial cell organisation, we 3D segmented 140 cells from a tube segment and 59 cells from a bud segment in CUBIC-cleared, light-sheet imaged embryonic lung explants (Figure 3a, Figure 2—figure supplement 1, Figure 2—video 1, Figure 3—videos 1–3). By interpolating between equally spaced sequential contour surfaces (every 1.66 µm in the tube and every 1.72 µm in the bud dataset) along the apical-basal axis, accurate volumetric reconstructions of cell morphology were obtained that allowed for the extraction of morphometric quantifications along the apical-basal axis. In both datasets, the 3D organisation of epithelial cells is highly complex, and cell neighbour relationships change continuously along the apical-basal axis (Figure 3b). As a result, cells are in direct physical contact not only with the cells that are neighbours on the apical side, but also with cells that appear two or even three cell diameters apart (Figure 3c).

Figure 3 with 4 supplements see all

Download asset Open asset

3D epithelial organisation.

(a) 3D iso-surfaces of segmented epithelial tip (N=59) and tube (N=140) cells from a Shh^GC/+; ROSA^mT/mG E12.5 mouse lung rudiment imaged using light-sheet microscopy. Morphometric quantifications of cell boundary segmentations along the apical-basal axis were used to study spatial T1 transitions. The 3D segmentation workflow is introduced in Figure 3—figure supplement 1 and Figure 3—videos 1–3 illustrate the rendered epithelial tip, tube and all segmented volumes. Scale bars 10 µm. (b) Frequency of neighbour numbers as quantified along the apical basal axis in the tip (solid lines) and tube (broken lines) datasets. (c) Extent of neighbour contacts (center cell in yellow and neighbours in blue) in 3D as viewed from the apical side. Scale bars 5 µm. (d) Probability distributions of the lateral T1 transitions for tip (total=169, mean=2.86, N=59) and trunk (total=746, mean=5.41, N=140) datasets are consistent with Poisson distributions. (e) Normalised apical-basal distribution of T1 transitions for all cells shows no apical-basal bias, except for fewer transitions close to the basal surface. (f). Schematic of a tubular epithelium with elliptic cross section. The analysed cells are located in the cusp (brown) of the tube, where the local tissue curvature is close to that of the minor ellipse vertex. Subpanel illustrates ellipse fitting to apical (purple) and basal (green) domain boundaries. Lighter outlines correspond to the most proximal segment of the tube in Figure 2b, while darker ones to the most distal section. (g) The predicted impact of a curvature effect on T1 transitions decreases with increasing cell neighbour numbers (blue line). The measured T1 transitions for different neighbour numbers do not support a curvature effect (dots, boxplots, and red line). Boxplots indicate the median, 25% and 75% percentiles of the data.

Figure 3—source data 1 3D neighbour number and area data.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig3-data1-v2.xls
Download elife-68135-fig3-data1-v2.xls
Figure 3—source data 2 3D curvature change data.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig3-data2-v2.xls
Download elife-68135-fig3-data2-v2.xls

Remarkably, we record up to 14 cell neighbour changes per cell in the tube and up to eight in the tip, between adjacent cross-sections along the apical-basal cell axis (Figure 3d). We will refer to these neighbour changes as lateral T1 transitions, or T1L. The mean relative apical-basal position for the lateral T1 transitions is 0.489±0.020 (95% CI), and there is no clear apical or basal tendency, though fewer transitions are observed close to the basal surface (Figure 3e). The dispersion index, that is, the ratio of the variance σ² and the mean number µ of transitions per cell, which equals unity for a Poisson distribution, is close to unity for both samples (Figure 3d). The chi-squared test also confirms that the number of apical-basal T1 transitions per cell is Poisson-distributed (Figure 3d). A Poisson distribution models the probability of a number of independent random events occurring in a given interval at a constant average rate. The consistency with a Poisson distribution, therefore, suggests a stochastic basis to the 3D organisation of epithelial cells.

The large number of observed T1L transitions and their distribution along the apical-basal axis challenges the recently popularised notion of curvature-driven scutoids as cell building blocks for epithelia (Gómez-Gálvez et al., 2018). To further examine the potential influence of tissue curvature on T1L transitions, we measured the apical-basal distance between two consecutive neighbour number changes for each cell in the tube dataset and recorded at which local tissue curvature they occur. For this analysis, we excluded ce ll portions from the apical end to the first transition and from the last transition to the basal end to reduce boundary effects, that is, only interior segments between transitions were considered. The mean apical-basal distance between two transitions is 17.89 ± 0.66 µm (95% CI). Local tissue curvature was approximated by fitting ellipses to the apical and basal surface boundaries of the tubular epithelium in 624 equidistant sections perpendicular to the main tube axis. (Figure 3f). Epithelial tubes in the developing lung are often collapsed (Conrad et al., 2021), making their apical and basal surfaces nearly elliptic in shape. The semi-axes of the fitted ellipses were then averaged over all sections to obtain the semi-axes a_apical, a_basal, b_apical, b_basal. Since our sample of 140 cells was segmented from a region close to the cusp of the nearly elliptical tube, a reasonably close estimate of the local tissue curvature where a T1L transition occurs is given by a linear interpolation between the curvature at the minor vertices of the apical and basal ellipses, according to the relative apical-basal position of the transition. The minor curvature of an ellipse with major and minor semi-axes a and b is given by b/a². Therefore, we estimate the local radius of curvature R by

R (x) = \frac{{{(a}_{a p i c a l} + x (a_{b a s a l} - a_{a p i c a l}))}^{2}}{b_{a p i c a l} + x (b_{b a s a l} - b_{a p i c a l})}

where x∈[0,1] is the relative apical-basal location of the T1L transition. The examined tissue exhibits an average curvature fold change of R(1)/R(0)=2.21 from the basal to the apical side. Denoting by R₁ and R₂ the radii of curvature between two adjacent T1 transitions along the apical-basal axis of a cell, we find that the distribution of curvature fold change R₂/R₁ shows no significant dependency on the number of neighbours n the cell has along that portion of the cell (Figure 3g). The mean curvature fold change per apical-basal T1L transition per cell is <R₂/R₁> = 1.142 ± 0.010 (95% CI). By extending the theory of scutoids (Gómez-Gálvez et al., 2018) to multiple T1L transitions per cell, we have derived a quantitative estimate of how tissue curvature would translate into the number of neighbour exchanges within that framework (Supplementary Material). If curvature changes were a main driver of T1L transitions, cells with smaller neighbour numbers n would be expected to change n over a much larger curvature fold change than cells with many neighbours (Figure 3g, blue line). However, we observe no systematic dependency of the curvature fold change on the number of neighbours the cell has along that portion of the cell in the developing mouse lung epithelium (Figure 3g). From this, we conclude that tissue curvature affects cell neighbourhood rearrangements through the tissue thickness at most mildly.

Neighbour changes along the apical-basal axis are driven by changes in cross-sectional area variation

Other than curvature effects, what else could drive the observed changes in neighbour relationships along the apical-basal axis? We notice that much as the apical and basal layers, each layer along the apical-basal axis behaves according to the three relationships previously described for the apical side, that is, Euler's polyhedron formula (Equation 1, Figure 4a), Aboav-Weaire's law (Equation 2, Figure 4b), and Lewis’ law (Equations 3, 4, Figure 4c). As predicted by the theory based on the minimisation of the lateral cell-cell energy (Kokic et al., 2019), the layers with a large area variability (Figure 4a) follow the quadratic law (yellow line) and those with a lower area variability the linear Lewis’ law (black line). The fraction of hexagons also follows the predicted relationship with the cross-sectional area variability (Figure 4d). We conclude that the neighbour relationships of epithelial cells along the entire apical-basal axis can be explained with a minimisation of the lateral cell-cell contact surface energy, as previously revealed for the apical layer.

Figure 4

Download asset Open asset

Neighbour changes along the apical-basal axis are driven by changes in cross-sectional area.

(a) Average number of neighbours (black) and area CV (blue) along the apical basal axis in the tip and tube datasets. (b) All epithelial layers follow the AW law (black line). The colour code in panels c, d follow that in panel b. (c) All epithelial layers follow Lewis’ law (Equation 2, black line) in case of low, and the quadratic relationship (Equation 3, yellow line) in case of high cell area variability. (d) Observed fraction of hexagons versus area CV for segmented cell layers along the apical-basal axis. The lines mark the theoretical prediction if polygonal cell layers follow either the linear Lewis’ law (black line) or the quadratic law (yellow line). (e) 3D iso-surfaces of four segmented epithelial cells in a CUBIC-cleared Shh^GC/+; ROSA^mT/mG E12.5 distal lung tube, with 140 3D segmented epithelial cells (Figure 3—figure supplement 1c). (f) Cross-sectional area and cell neighbour number along the apical-basal axis for marked cells in panel e. Dotted lines indicate contact with one another. (g) Lateral cross-sections illustrating a T1 transition along the apical-basal axis (0.664 µm in-between frames). Scale bar 6 µm. (h) An increase in the cell cross-sectional area increases the frequency of a neighbour number increasing spatial T1 transitions, and vice versa. (i) Apical-basal distribution of spatial T1 transitions according to neighbour increase or decrease and cross-sectional area variation.

Figure 4—source data 1 3D area CV and avg. neighbour number data.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig4-data1-v2.xls
Download elife-68135-fig4-data1-v2.xls
Figure 4—source data 2 3D AW data.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig4-data2-v2.xls
Download elife-68135-fig4-data2-v2.xls
Figure 4—source data 3 Neighbour exchange/type data.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig4-data3-v2.xls
Download elife-68135-fig4-data3-v2.xls

If epithelial cell neighbour relationships are indeed driven by a minimisation of the total lateral cell-cell contact surface energy then the T1L transitions along the apical-basal axes should be driven by changes in the cross-sectional area along the apical-basal axis. If we analyse four 3D segmented cells (Figure 4e) in detail, we indeed see how an increase in the cross-sectional area results in an increase in the neighbour number, and vice versa (Figure 4f) via lateral T1 transitions (Figure 4g). As the cell neighbour arrangements represent global minima, the local analysis does, of course, not provide a perfect correlation. When we consider all 140 segmented cells in the tube segment with their 746 cell neighbour exchanges between adjacent cross-sections (Figure 3e), then we find that the frequency of T1L transitions along the apical-basal axis is indeed higher, the larger the increase in cross-sectional area, and vice versa (Figure 4h,i).

Changes in cross-sectional area as a result of interkinetic nuclear migration (IKNM)

So, what determines the cross-sectional cell areas in each layer? In pseudostratified epithelia, mitosis is restricted to the apical surface (Gundersen and Worman, 2013). Depending on the average diameter of nuclei and the average apical cross-sectional area, there is insufficient space for all nuclei to be accommodated apically. Therefore, as a cell exits mitosis, the nucleus moves from the apical towards the basal side (G1 and S phase) and then back to the apical side (G2 phase) to undergo another round of mitosis, a process referred to as interkinetic nuclear migration (IKNM) (Meyer et al., 2011). Consequently, nuclei are distributed along the entire apical-basal axis, giving the tissue a pseudostratified configuration (Norden, 2017). We wondered to what extent the nuclear distribution, and its effect on the 3D cell shape, explains the observed area distributions and lateral T1 transitions.

To this end, we stained the nuclear envelope with fluorescently tagged antibodies against lamin B1 (Figure 5a, Figure 5—figure supplement 1), and 3D segmented all nuclei within epithelial cells in a tube segment (Figure 5b, Figure 3—figure supplement 1, Figure 3—video 2). The nuclei were distributed along the entire apical-basal axis (Figure 5c), and consistent with the pseudostratified appearance of the epithelium, nuclei in neighbouring cells had different positions along the apical-basal axis (Figure 5b). The nuclear shapes, volumes, and cross-sectional areas (Figure 5d–f) all varied along the apical-basal axis. As expected, nuclei are largest and most spherical at the apical side, where they undergo mitosis (Figure 5d). Thus, a one-sided, two-sample Welch t-test revealed a significantly reduced ellipticity of nuclei located in the first 25% of the apical-basal axis compared to those in the middle 50% (p=0.0002). The nuclear volumes are on average about 50% smaller than the cell volumes and largely correlate (r = 0.79, Figure 5f), likely reflecting parallel expansion during the cell cycle. Where present, the nuclear cross-sectional areas are only slightly smaller than those of the entire cell, and the cross-sectional areas of the cell and the nucleus are strongly correlated (r = 0.94, Figure 5g). The strong correlation can be accounted for by the opposing actions of cells and nuclei in the columnar epithelium. The nuclear volumes are too large to allow for a spherical nucleus to fit into a cylindrical cell of the measured height (Figure 5h). Accordingly, to fit into the cell, the nucleus necessarily has to deform. Nuclei respond to external forces with anisotropic shape changes (Haase et al., 2016; Neelam et al., 2016), which is consistent with the elliptical nuclear shapes that we observe (Figure 5d). However, there is a limit to how much the stiff nucleus can deform (Lammerding, 2011; Shah et al., 2021), resulting in a local widening of the cell where the nucleus is present. Cell sections without nucleus typically have smaller cross-sectional areas, thereby leading to a higher frequency of small cross-sections in cells compared to nuclei. Accordingly, as previously seen for the cell cross-sectional areas, the observed changes in cell neighbour numbers correlates with the observed changes in nuclear cross-sectional areas (Figure 5i) such that most T1L transitions occur where the nucleus starts and ends (Figure 5j).

Figure 5 with 2 supplements see all

Download asset Open asset

Changes in cross-sectional area as a result of interkinetic nuclear migration (IKNM).

(a) Light-sheet microscopy longitudinal sections of an E12.5 CUBIC-cleared lung tube carrying the Shh^GC/+; ROSA^mT/mG reporter allele (green epithelium) and immunostained for lamin B1 (blue nuclear envelopes). Morphometric quantifications of 3D iso-surfaces (N=140) and cell segmentations along the apical-basal axis were used to study the nature of cross-sectional area variation and the effect of IKNM (Figure 3—figure supplement 1). Scale bar 20 µm. (b) Sequential cell membrane contour surfaces and nuclear iso-surfaces for six epithelial cells. By interpolating between contours and creating iso-surfaces, 3D shapes can be accurately extracted (Figure 3—video 2 and Figure 5—video 1). Scale bar 7 µm. (c) Distribution of nuclei center of mass along the apical-basal axis. (d) Nuclear ellipticity and volume distributions along the apical-basal axis. (e) Average cross-sectional area distribution along the apical-basal axis for all cells (black) and nuclei (blue). (f). Nuclear and cellular volumes of 140 segmented cells are correlated (r=0.79). Line of best fit in red. (g). The cell and corresponding nuclear cross-sectional areas along the apical-basal axis are highly correlated (r=0.94). (h). Diameters of nuclei and cells based on the measured nuclear and cell volumes if nuclei were perfect spheres, and cells were perfect cylinders of the measured height. Given the larger nuclear diameter, nuclei must deform in order to fit within cells. (i) An increase in the nuclear cross-sectional area increases the frequency of a neighbour-number-increasing spatial T1 transitions, and vice versa. (j) The largest number of changes in neighbour relationships occur at the apical and basal limits of the nucleus for all cells, where cross-sectional areas change sharply.

Figure 5—source data 1 3D neighbour number, cell area, and nuclear area data.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig5-data1-v2.xls
Download elife-68135-fig5-data1-v2.xls
Figure 5—source data 2 Volume and ellipticity data.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig5-data2-v2.xls
Download elife-68135-fig5-data2-v2.xls
Figure 5—source data 3 3D cell and nuclear diameter data.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig5-data3-v2.xls
Download elife-68135-fig5-data3-v2.xls

We conclude that the positions of nuclei can explain much of the observed variability in the cross-sectional cell areas. During the cell cycle, nuclei migrate, and the cell volumes first increase, and subsequently halve due to cell division. As all these processes affect the cross-sectional areas of the cells along the apical-basal axis, one would expect continuous spatial-temporal T1L transitions in growing pseudostratified epithelia.

3D cell organisation in growing epithelia

To follow 3D cellular dynamics during epithelial growth and deformation, we cultured embryonic lungs from a Shh^GC/+; ROSA^mT/mG background and imaged every 20 min for a total of 10 hr using light-sheet microscopy (Figure 6—video 1). We used a subset of this dataset (11 time points, >3 hr) (Figure 6—video 2) to 2.5D segment the apical and basal surfaces, and to explore 3D cell shape dynamics and neighbour relationships in a growing lung bud (Figure 6a, Figure 6—video 3).

Figure 6 with 4 supplements see all

Download asset Open asset

Dynamics of apical and basal epithelial organisation.

(a) Timelapse light-sheet microscopy series of a cultured mouse E12.5 distal lung bud expressing the Shh^GC/+; ROSA^mT/mG reporter (green epithelium, and red mesenchyme), imaged every 20 min (11 time steps). The white inset denotes the morphology of the lung bud, while the dotted area denotes the segmented cell patch. Corresponding visual provided in Figure 6—video 1. Cells on both the apical and basal domains were 2.5D segmented, and their morphology quantified. Corresponding visual provided in Figure 6—video 3. Scale bars 20 µm. (b) Cell neighbour frequencies for the apical and basal layers over time. (c) Observed average neighbour number and area coefficient of variation (CV) for the apical and basal layers over time. (d) Growing apical and basal layers follow the AW law (black line). Colour code applies to e-f. (e) The relative average apical and basal cell areas are linearly related to the number of neighbours (in all time points) and follow Lewis’ law (black line), or the quadratic relationship in the case of higher area variability (yellow line). (f) Observed fraction of hexagons versus area coefficient of variation (CV) on the apical and basal layers. The lines mark the theoretical prediction if polygonal cell layers follow either the linear Lewis’ law (black line) or the quadratic law (yellow line).

Figure 6—source data 1 2.5D neighbour number and area data over time.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig6-data1-v2.xls
Download elife-68135-fig6-data1-v2.xls
Figure 6—source data 2 2.5D area CV and avg. neighbour number data over time.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig6-data2-v2.xls
Download elife-68135-fig6-data2-v2.xls
Figure 6—source data 3 2.5D AW data over time.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig6-data3-v2.xls
Download elife-68135-fig6-data3-v2.xls

As the explant was growing, we readjusted the 2.5D segmented region such that the segmented surface area and cell numbers remained roughly constant over time (Figure 6—figure supplement 1a). Nonetheless, the segmented bud increased in volume as the thickness of the layer increased with time (Figure 6—figure supplement 1b). Much as in the static dataset, the neighbour number distributions (Figure 6b), and variability of cross-sectional areas (Figure 6c) differ between the apical and basal cell layers in all time points. However, for all time points, both the apical and basal layers conformed to Euler's polyhedron formula (Figure 6c), Aboav-Weaire's law (Figure 6d), and Lewis’ law (Figure 6e). Furthermore, the fraction of hexagons also followed from the variability of the cross-sectional areas, as predicted by the theory (Figure 6f).

We next sought to analyse the 3D dynamics of segmented epithelial cells. As the tracking of packed cells in growing pseudostratified epithelia is challenging, we focused on a small patch with 15 cells in total (Figure 7a). Sequential contour surfaces were drawn to follow cell membrane outlines on several planes along the apical-basal axis and interpolated to reconstruct 3D morphology for each time point (Figure 7—figure supplement 1, Figure 7—video 1). All planar segmentations were then pooled into five groups along the apical-basal axis to enable morphometric analysis in different tissue regions. Over the time course, the volume of individual cells varied between roughly 400 and 800 µm³ (Figure 7b), and the apical-basal length varied between roughly 20 and 30 µm (Figure 7c). We note that all layers conform to Aboav-Weaire's law (Figure 7d) and Lewis’ law (Figure 7e) in all time points. Moreover, consistent with our theory, the fraction of hexagons follows from the variability of the cross-sectional area, though more deviations are observed, given the small number of cells analysed (Figure 7f).

Figure 7 with 2 supplements see all

Download asset Open asset

3D cell organisation in growing epithelia.

(a) 3D segmentation of 15 epithelial cells from timelapse light-sheet microscopy imaging of a mouse E12.5 distal lung bud expressing the Shh^GC/+; ROSA^mT/mG reporter. The specimen was imaged every 20 min over 3 hr. Planar segmentations along the apical-basal axis were pooled into five groups to enable morphometric analysis in different tissue regions. A full timelapse panel is provided in Figure 7—figure supplement 1 and Figure 7—video 1. Scale bar 30 µm. (b) Epithelial cell volume, and (c) height over time (N=15). (d) Segmented cells in pooled layers along the apical-basal axis follow the AW law (black line) over time (left to right); see panel f for colour code. (e) The relative average cell area in each layer is linearly related to the number of neighbours for all time points (left to right) and follows Lewis’ law (black line), or the quadratic relationship in the case of higher area variability (yellow line); see panel f for colour code. (f) Temporal dynamics of observed fraction of hexagons versus area coefficient of variation (CV) along the apical-basal axis. Dotted lines denote variation per time point. Solid lines mark the theoretical prediction if polygonal cell layers follow either the linear Lewis’ law (black line) or the quadratic law (yellow line).

Figure 7—source data 1 3D neighbour number and area data over time.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig7-data1-v2.xls
Download elife-68135-fig7-data1-v2.xls
Figure 7—source data 2 3D area CV and hex percent data over time.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig7-data2-v2.xls
Download elife-68135-fig7-data2-v2.xls
Figure 7—source data 3 3D AW data over time.: https://cdn.elifesciences.org/articles/68135/elife-68135-fig7-data3-v2.xls
Download elife-68135-fig7-data3-v2.xls

Much as in the static dataset (Figures 3 and 4), we observe up to 14 neighbour number changes (T1L transitions) along the apical-basal axis (Figure 8a,b). The average number of T1L transitions is relatively constant over time (Figure 8b). The mean relative apical-basal position for T1L transitions is again roughly in the middle, but in this small dataset, we now observe more T1L transitions in the center of the cell than at the apical or basal boundaries (Figure 8c). By following a single cell over time, we can appreciate the dynamic cell shape changes, and how a change in the cross-sectional area correlates with a change in neighbour number (Figure 8d). The neighbour relationships are, of course, not determined by the local cell cross-section, but by the overall cross-sectional area distribution in that layer. Accordingly, the correlation between the cross-sectional area and the neighbour number is not perfect for a single cell. By considering a patch of cells, we can, however, see how those T1L transitions occur dynamically in developing tissues (Figure 8e).

Figure 8

Download asset Open asset

3D cell neighbour dynamics in growing epithelia.

(a) Apical, basal and lateral cross-sections from a light-sheet microscopy timelapse of a murine E12.5 distal lung bud expressing the Shh^GC/+; ROSA^mT/mG reporter. The specimen was imaged every 20 min over 3 hr. Cell membrane outlines illustrate fluid cell neighbour relationships along the apical-basal axis and over time. T1 transitions are marked with white stars; see panel e for cell colour code. Scale bar 14 µm. (b) Number of T1L transitions for all cells (N=15) over time. Diamonds represent the mean. Morphometric quantifications of planar segmentations along the apical-basal axis were used to examine T1L transition dynamics. (c) Spatial distribution along the apical-basal axis of T1L transitions for all cells and time points. (d) Temporal evolution of neighbour relationships along the apical-basal axis for a single cell. Schematic cell width corresponds to cross-sectional area. (e) (top row) 3D iso-surface segmentations of 15 epithelial cells. Scale bar 10 µm. (bottom row) Cross-sectional area and cell neighbour number along the apical-basal axis for a given cell (marked with a white star). Dotted lines indicate contact with a cell that was not segmented.

Discussion

Epithelial tissues remodel into complex geometries during morphogenesis. We used light-sheet microscopy and 3D cell segmentation to unravel the physical principles that define the 3D cell neighbour relationships in pseudostratified epithelial tissues. Our analysis reveals that pseudostratified epithelial layers adopt a far more complex packing solution than previously anticipated: the 3D epithelial cell shapes are highly irregular, and cell neighbour relationships change multiple times along the apical-basal axis, with some cells having up to 14 changes in their neighbour contacts along their apical-basal axis (Figure 3a). Curvature effects can result in neighbour changes, but the data does not show the dependency on cell neighbour numbers that would be expected if curvature effects played a dominating role (Figure 3g). There is also no apical-basal bias (Figure 3e), and the prevalence of contact remodelling is randomly distributed (Figure 4i).

Even though the neighbour relationships are uncorrelated between the apical and basal sides and appear random at first sight, they follow the same fundamental relationships that have previously been described for apical epithelial layers, that is, Euler's polyhedron formula, Lewis' law, and Aboav-Weaire's law across the entire tissue and at all times (Figures 2, 4, 6 and 7). This arrangement minimises the lateral cell-cell surface energy in each plane along the apical-basal axis, given the variability in the cell cross-sectional areas (Figure 1; Kokic et al., 2019; Vetter et al., 2019). Where present, the stiff nucleus determines the cell cross-sectional area, as is apparent from the strong correlation between the cell cross-sectional and the nuclear cross-sectional areas (Figure 5g). Accordingly, most changes in neighbour relationships occur at the apical and basal limits of the nucleus where cross-sectional areas change sharply (Figure 5i). As the nucleus moves along the apical-basal axis during the cell cycle, a process referred to as interkinetic nuclear migration (IKNM) (Meyer et al., 2011), cell neighbour relationships change continuously (Figure 8). We conclude that neighbour relationships in epithelia are fluidic, and the complex, dynamic 3D organisation of cells in growing epithelia follows simple physical principles.

Defining the physical principles behind cell neighbour relationships is only the first step in unravelling the determinants of epithelial 3D cell shapes. The second key aspect is the cell volume distribution along the apical-basal axis, which gives rise to the cell cross-sectional area distribution, which then determines the cell neighbour relationships (Figure 5e). The overall cell volume is determined by cell growth and division, but its distribution along the apical-basal axis depends on the nuclear dynamics (Figure 5g), and the epithelial cell heights. We find that the nucleus occupies, on average, 55% of the cell volume in the embryonic lung epithelia. As the cell nuclei move along the apical-basal axis during the cell cycle (Meyer et al., 2011), the cytoplasm fills the remaining space between the apical and basal surfaces, likely in a way that minimises the total surface area of all cells. The determinants of the epithelial thickness, that is, the distance between the apical and basal surfaces are still unknown, but signalling downstream of Fibroblastic Growth Factor (FGF), Sonic Hedgehog (SHH), Bone Morphogenetic Protein (BMP)/ transforming growth factor-beta (TGF-β), and WNT has been observed to affect cell height, presumably via an effect on cell tension and/or cell-cell adhesion (Gritli-Linde et al., 2002; Hirashima and Matsuda, 2021; Kadzik et al., 2014; Kondo and Hayashi, 2015; Widmann and Dahmann, 2009).

Cell-based modelling frameworks are heavily used to investigate epithelial processes and how they result in morphological changes such as tissue bending, folding, fusion, and anisotropic growth during morphogenesis (Fletcher et al., 2014; Tanaka, 2015). Our data confirms many underlying assumptions of cell-based modelling frameworks and provides quantitative data to calibrate parameters. Once calibrated to reproduce the here identified 3D cell shape distributions, such simulation frameworks will help to reveal the determinants of 3D cell shapes, and will be invaluable in providing insight into how local changes in cell growth, adhesion, tension, or in the basal lamina affect cell shapes locally and within the remaining epithelial layer.

In summary, this study offers a detailed view of 3D cell neighbour relationship dynamics and packing in growing epithelial tissues, and demonstrates that the 3D cell shapes are much more complex than previously anticipated, and that cell neighbour relationships are dynamic and change as result of cell growth and cell cycle-linked IKNM. The complex 3D cell neighbour relationships can nonetheless be understood based on simple physical principles. Although we recognize that tissue architecture is a multifactorial process, our work carries vast implications for the study of cell-cell signalling, epithelial cohesion, and energetic modelling of developing epithelial layers in both healthy and disease contexts.

Materials and methods

Ethical statement

Request a detailed protocol

Permission to use animals was obtained from the veterinary office of the Canton Basel-Stadt (license number 2777/26711). Experimental procedures were performed in accordance with the Guide for the Care and Use of Laboratory Animals and approved by the Ethics Committee for Animal Care of ETH Zurich. All animals were housed at the D-BSSE/UniBasel facility under standard water, chow, enrichment, and 12 hr light/dark cycles.

Animals

To investigate 3D cellular dynamics during mouse embryonic lung development, we used mouse lung rudiments from animals homozygous for the Rosa26^mTmG and heterozygous for the Shh^GFP-Cre allele (Shh^GFP-Cre/+; Rosa26^mTmG). The double-fluorescent Shh-controlled Cre reporter mouse expresses membrane-targeted tandem dimer Tomato (mT) before CRE-mediated excision and membrane-targeted green fluorescent protein (mG) after excision (Muzumdar et al., 2007). As a result, only epithelial cell membranes are labelled by GFP, while all adjacent mesenchymal tissue is labelled by tdTomato.

Immunofluorescence

Request a detailed protocol

E12.5 mouse lungs were fixated for 1 hr in 4% paraformaldehyde in PBS, and subsequently incubated with Lamin B1 (Thermo; Material No. 702972; 1:200) at 4°C for 3 days. As a structural component of the nuclear lamina, LaminB1 immunostaining makes crowded nuclei clearly distinguishable and easily segmentable. After washing in D-PBS, lungs were incubated with conjugated fluorescent secondary Alexa Fluor 555 donkey anti-mouse IgG (H+L) (Abcam; Material No. ab150106; 1:250) for 2 days at 4°C.

Optical clearing and lightsheet imaging

Request a detailed protocol

Optical clearing of embryonic lung rudiments enabled the 3D segmentation of numerous epithelial cells from single image stacks. To this extent, the whole-mount clearing of dissected E12.5 lung explants was performed with the Clear Unobstructed Brain/Body Imaging Cocktails and Computational Analysis (CUBIC) protocol (Susaki et al., 2015; Figure 2—figure supplement 3). Reagents for delipidation and refractive index (RI) matching were made as follows: CUBIC-1 [25% (w/w) urea, 25% ethylenediamine, 15% (w/w) Triton X-100 in distilled water], and CUBIC-2 [25% (w/w) urea, 50% (w/w) sucrose, 10% (w/w) nitrilotriethanol in distilled water], respectively. Following fixation and immunostaining, samples were incubated in 1/2 CUBIC-1 (CUBIC-1:H2O=1:1) for four days, and in 1X CUBIC-1 until they became transparent. All explants were subsequently washed several times in PBS and treated with 1/2 CUBIC-2 (CUBIC-2:PBS=1:1) for around four days. Lastly, incubation in 1X CUBIC-2 was done until the desired transparency was achieved. All solutions were changed daily, and CUBIC-1 steps were performed on a shaker at 37°C while CUBIC-2 steps at room temperature. Cleared samples were then embedded in 2% low melting point solid agarose cylinders and immersed in CUBIC-2 for two more days to increase the agarose refractive index. 3D image stacks were acquired on a Zeiss Lightsheet Z.1 microscope using a Zeiss 20x/1.0 clearing objective .

Timelapse light-sheet acquisitions

Request a detailed protocol

Light-sheet acquisitions of live epithelial cell morphology enabled the study of 3D organisation dynamics. Following dissection in DPBS at room temperature, E12.5 lung explants were cultured in sterile Dulbecco's modified Eagle's medium w/o phenol red (DMEM) (Life Technologies Europe BV; 11039021) containing 10% Fetal Bovine Serum (FBS) (Sigma-Aldrich Chemie GmbH; F9665-500ML), 1% Glutamax (Life Technologies Europe BV; A1286001), and 1% penicillin/streptomycin (Life Technologies Europe BV; 10378–016). All specimens were equilibrated at 37°C with 5% CO2 in a humidified incubator for 1 hr.

Following a 1 hr equilibration period, 1.5% LMP hollow agarose cylinders were prepared (Udan et al., 2014). Hollow cylinders, in contrast to solid ones, accommodate unencumbered 3D embryonic growth, provide boundaries to minimise tissue drift, enable imaging from multiple orientations, and allow for better perfusion of gasses and nutrients. All specimens were suspended within each hollow cylinder in undiluted Matrigel (VWR International GmbH; 734–1101), an ECM‐based optically clear hydrogel that provided a near-native 3D environment and supported cell growth and survival. All cylinders were kept at 37°C with 5% CO₂ in culture media for 1 hr before mounting.

For each overnight culture, the imaging chamber was prepared by sonication at 80°C, followed by ethanol and sterile PBS washes. After assembly, the chamber was filled with culture medium and allowed to equilibrate at 37°C with 5% CO₂ for at least 2 hr before a cylinder containing an explant was mounted for imaging. Furthermore, to compensate for evaporation and to maintain a fresh culture media environment, two peristaltic pumps were installed to supply 0.4 mL and extract 0.2 mL of culture medium per hour. Each lung explant was then aligned with the focal plane within the center of a thin light-sheet to enable fine optical sectioning with optimal lateral resolution. For this study, all live imaging was done with a 20x/1.0 Plan-APO water immersion objective.

Image processing

Request a detailed protocol

To efficiently process the resulting volumetric CZI datasets (10s-100s of GBs), all image stacks were transferred to a storage server and subsequently processed in remote workstations (Intel Xeon CPU E5-2650 with 512 GB memory). Deconvolution via Huygens Professional v19.04 (Scientific Volume Imaging, The Netherlands, http://svi.nl) improved overall contrast and resolution while Fiji (ImageJ v1.52t) (Schindelin et al., 2012) aided in accentuating cell membranes, enhancing local contrast, removing background fluorescence, and TIFF conversion.

Cell morphometric quantifications

Request a detailed protocol

Cell morphology on the apical and basal membranes of embryonic lung epithelia was investigated using the open-source software platform MorphoGraphX (MGX) (Barbier de Reuille et al., 2015). By meshing the curved boundaries of input 3D image stacks and projecting nearby signal onto it, MGX builds a curved 2.5D image projection that is distortion-free, unlike planar 2D projections that ignore curvature. We then proceeded to use a suitable implementation of the Watershed transform to extract individual cell geometries, with minimal manual curation, and quantify properties such as surface area and the number of cell neighbours. All border cells were excluded. Apical and basal cell meshes were exported as text files and traversed using the R Programming Environment to extract the neighbour relationships between cells as needed to generate Aboav-Weaire plots (Gómez et al., 2021).

To render time-lapse datasets and extract 3D volumetric surface reconstructions of entire epithelial cells, we employed Imaris v9.1.2 (Bitplane, South Windsor, CT, USA). By computationally interpolating between cell membrane contour surfaces from successive transverse frames into iso-surfaces, faithful cell and nuclear 3D volumes were obtained. Quantified volumetric features included cell and nuclear volume, total surface area, sphericity, and nuclear position along the apical-basal axis. Imaris was also used to generate high-resolution videos, which, despite being strongly downsampled to accommodate vast time-lapse datasets, presented little noticeable loss in image quality. Furthermore, to extract cell areas and the number of neighbours along the apical-basal axis, transverse image frames were imported into ImageJ and processed using the interactive plugin TissueAnalyzer (Aigouy et al., 2016). Like this, cell segmentation masks across layers could be generated, and cell geometry and neighbour topology quantified (Gómez et al., 2021).

Appendix 1

A simple theory for lateral T1 transitions in curved tissues

Drawing inspiration from the notion of ‘scutoids’ (Gómez-Gálvez et al., 2018), we develop here a simple geometrical theory to estimate the effect of tissue curvature, if any, on the occurrence and number of lateral T1 transitions of cells along the apical-basal axis in a pseudostratified epithelium. The purpose of this effort is to verify whether an idea brought forward in the aforementioned article, which is that tissue curvature could be responsible for the occurrence of up to one lateral T1 transition per cell between its apical and basal sides, is consistent with our observations in the developing, pseudostratified mouse lung epithelium.

We start by considering, like the authors of Gómez-Gálvez et al., 2018, the vicinity of an edge shared by two adjacent cells in a two-dimensional cross section of the tissue perpendicular to the apical-basal axis. This vicinity is defined as the quadrilateral (Appendix 1—figure 1A, orange) spanned by the four cell junctions (green) that are directly linked to the edge of interest (blue), covering a fraction of four cells in such a two-dimensional projection. As the projection plane is moved along the apical-basal axis, the edge length can shrink to zero and reappear with different orientation, leaving the cell neighbourhood of all four cells modified by one. It is this process which we refer to as lateral neighbourhood transition of type T1. Each such edge vicinity is characterised by an aspect ratio $ϵ = w / h$ , where $w = (w_{1} + w_{2}) / 2$ is the average width and $h = (h_{1} + h_{2}) / 2$ the average height. We now proceed to quantitatively estimate with a simple geometrical model how a change in tissue curvature translates into a change of aspect ratio of such motifs, and how that in turn would be expected to change the cell neighbourhood if all T1 transitions were effectively governed by the geometric effect of changing tissue curvature.

Appendix 1—figure 1

Download asset Open asset

Definition of the vicinity of a cell edge in a geometric projection on a tubular epithelium.

We consider two different tissue topologies: (i) a cylindrical tube and (ii) a spherical vesicle or hemispherical tube cap. The crucial difference between these two cases is that the former only has a mean curvature but no Gaussian curvature, whereas the latter has both. As the radius $r$ of these surfaces changes, the change of aspect ratio follows by the quotient rule of calculus as

\frac{\partial ϵ}{\partial r} = \frac{\partial (w / h)}{\partial r} = \frac{1}{h} (\frac{\partial w}{\partial r} - ϵ \frac{\partial h}{\partial r}) .

We now consider the general case where the edge motif is arbitrarily oriented on the surface, and will later average over all possible orientations. In general, the edge of interest spans an angle $θ$ with the first principal curvilinear axis of the surface, which we take as the main cylinder axis for the tubular epithelium (Appendix 1—figure 1), and as any great circle on the spherical epithelium. The change of width and height then read

\begin{array}{ll} \frac{\partial w}{\partial r} & = \cos θ \frac{\partial w_{x}}{\partial r} + \sin θ \frac{\partial w_{y}}{\partial r}, \\ \frac{\partial h}{\partial r} & = \sin θ \frac{\partial h_{x}}{\partial r} + \cos θ \frac{\partial h_{y}}{\partial r}, \end{array}

where w_x, w_y and h_x, h_y are the components of the width and height vectors, respectively (Appendix 1—figure 1B, gray):

\begin{array}{ll} w_{x} & = w \cos θ, w_{y} = w \sin θ, \\ h_{x} & = h \sin θ, h_{y} = h \cos θ . \end{array}

Their response to a change in tissue radius is given by

\begin{array}{ll} \frac{\partial w_{x}}{\partial r} & = κ_{1} w_{x}, \frac{\partial w_{y}}{\partial r} = κ_{2} w_{y}, \\ \frac{\partial h_{x}}{\partial r} & = κ_{1} h_{x}, \frac{\partial h_{y}}{\partial r} = κ_{2} h_{y}, \end{array}

where $κ_{1}$ and $κ_{2}$ are the principal curvatures of the surface in the respective directions. Substituting Equation 8 into Equations 7, 6 and 5 yields

\frac{\partial ϵ}{\partial r} = ϵ (κ_{1} - κ_{2}) \cos (2 θ) .

A cylinder with radius $r$ has the principal curvatures $κ_{1} = 0$ and $κ_{2} = 1 / r$ , so

\frac{\partial ϵ}{\partial r} = - \frac{ϵ}{r} \cos (2 θ) .

For a sphere with radius $r$ , on the other hand, the two principal curvatures are both equal, $κ_{1} = κ_{2} = 1 / r$ , and thus the aspect ratio does not change with varying radius,

\frac{\partial ϵ}{\partial r} = 0

independent of the orientation $θ$ .

From here onward, we consider the region along the apical-basal axis between two T1 transitions on a given cell, irrespective of whether these T1 transitions increase or decrease the neighbour number. In our experimental quantifications of lateral T1 transition in the mouse lung epithelium, they were counted in the same way. The sign of the T1 transition is encoded in the sign of $\partial ϵ / \partial r$ , hence we consider only the absolute value $| \partial ϵ / \partial r |$ hereafter, to quantify the change of edge aspect ratio between T1 transitions along the tissue radius.

Carrying on with the tubular geometry, we now average over all orientations $θ$ to calculate the mean change of aspect ratio for the entire cell ensemble in the tissue. The set of possible orientations is $θ \in [0, π / 2)$ , because all angles $θ \in [π / 2, π)$ just switch the roles of $w$ and $h$ . Therefore,

| \frac{\partial ϵ}{\partial r} | = \frac{2}{π} \int_{0}^{π / 2} | \frac{ϵ}{r} \cos (2 θ) | d θ = \frac{2 ϵ}{π r} .

This is an ordinary differential equation which can easily be solved by separation of variables. The result is

ϵ (r) = ϵ_{1} {(\frac{R_{1}}{r})}^{\pm 2 / π}

and therefore

| \frac{\partial ϵ}{\partial r} | = \frac{2 ϵ_{1}}{π r} {(\frac{R_{1}}{r})}^{\pm 2 / π}

where $ϵ_{1}$ is the aspect ratio at a given cylinder radius R₁. For the sphere, on the other hand, we trivially have a constant aspect ratio $ϵ (r) = ϵ_{1}$ .

To link this change of aspect ratio with lateral T1 transitions, a relationship between $ϵ$ and the number of neighbours a cell locally has, $n$ , needs to be established. A simple rough estimate can be found by recalling that cells typically tend toward a regular polygonal shape in a cross-sectional projection (Kokic et al., 2019; Vetter et al., 2019). In a honeycomb lattice, the aspect ratio is

{ϵ |}_{n = 6} = \frac{2}{\sqrt{3}} \approx 1.155 .

For general regular polygons with $n \geq 4$ , the relationship reads

ϵ (n) = \frac{1 + 2 \cos (2 π / n)}{2 \sin (2 π / n)}

as can be recognised from Appendix 1—figure 1B. Equation 10 holds locally in the region between two consecutive T1 transitions, where $n$ is constant. From this, we can estimate the change of aspect ratio between two lateral T1 transition as

\frac{\partial ϵ}{\partial n} = π \frac{2 + \cos (2 π / n)}{n^{2} \sin^{2} (2 π / n)} .

On average, for hexagonal cells ( $n = 6$ ), this evaluates to

{\frac{\partial ϵ}{\partial n} |}_{n = 6} = \frac{5 π}{54} \approx 0.291 .

So how many T1 transitions can be expected along the apical-basal axis? We can apply the chain rule of calculus to write

| \frac{d n}{d r} | = | \frac{\partial n}{\partial ϵ} \frac{\partial ϵ}{\partial r} | = | {(\frac{\partial ϵ}{\partial n})}^{- 1} \frac{\partial ϵ}{\partial r} | .

Substituting Equations 9 and 11 into Equation 12, we find

| \frac{d n}{d r} | = \frac{n^{2} \sin^{2} (2 π / n)}{2 + \cos (2 π / n)} \frac{2 ϵ_{1}}{π^{2} r} {(\frac{R_{1}}{r})}^{\pm 2 / π}

for a cylindrical epithelium, and $d n / d r = 0$ for a spherical one. A spherical topology is thus entirely compatible with cells being frusta that do not change their neighbour number from their apical to the basal side, whereas a cylindrical topology is not. For tubes, a change of radius $R_{1} \to R_{2}$ lets the neighbour number change by

| Δ n | = \int_{R_{1}}^{R_{2}} | \frac{d n}{d r} | d r = \frac{ϵ_{1}}{π} \frac{n^{2} \sin^{2} (2 π / n)}{2 + \cos (2 π / n)} [1 - {(\frac{R_{1}}{R_{2}})}^{\pm 2 / π}]

per cell edge. How far must R₁ and R₂ be apart to expect one T1 transition per cell? Over that distance, each of the cell’s $n$ edges will have undergone a change of $1 / n$ on average. The answer can therefore be found by setting $| Δ n | = 1 / n$ with $ϵ_{1} = ϵ (n)$ from Equation 10, and solving for the radius fold change:

\frac{R_{2}}{R_{1}} = {(1 - \frac{ϵ^{'} (n)}{n ϵ (n)})}^{\pm π / 2} = {(1 - \frac{2 π [2 + \cos (2 π / n)]}{n^{2} [\sin (2 π / n) + \sin (4 π / n)]})}^{\pm π / 2}

This is the expected fold change of mean curvature inside a tubular epithelium between two consecutive lateral T1 transitions per cell, if the T1 transitions were mainly driven by a change of curvature. Note that the two cases with different signs in the exponent reflect the $R_{1} \leftrightarrow R_{2}$ symmetry. As shown in Appendix 1—figure 2, Equation 13 predicts a strong dependency of the required change of curvature between adjacent lateral T1 transitions on the number of neighbours the cell has, $n$ . Cells with a large neighbour number are expected to transition at much smaller changes of curvature on average than cells with few neighbours.

Appendix 1—figure 2

Download asset Open asset

Predicted radius fold change in a tubular epithelium between a pair or consecutive T1 transitions per cell.

Equation 13 relates the fold change of tissue curvature on a tubular sheet between consecutive lateral T1 transitions along the apical-basal axis of a single cell, to the cell’s neighbour number $n$ in any cross section perpendicular to the radius between the two transitions. $n$ here is to be understood locally within the tissue, not as the cell connectivity on either the apical or basal surfaces of the cell. Equation 13 does not directly lend itself to estimating the number of lateral T1 transitions to be expected for cells in a curved tissue, although our theory might be adaptable to this case with some modifications. The sign of T1 transitions is ignored here, leaving the possibility of frequent back-and-forth transitions with zero net effect on $n$ —a phenomenon that we indeed do observe in the developing lung epithelium. Equation 13 quantifies only the net contribution of curvature on the number of T1 transitions a cell undergoes. This net effect strongly depends on the the local neighbour number $n$ , as shown in Appendix 1—figure 2.

Data availability

The source code and plotted data files are available as a git repository at https://git.bsse.ethz.ch/iber/Publications/2021_gomez_3d_cell_neighbour_dynamics.git. The raw data is publicly available as openBIS repository at https://openbis-data-repo.ethz.ch/openbis/webapp/eln-lims/?user=observer&pass=openbis under the name 3D Epithelium.

The following data sets were generated

1. Gómez HF
2. Dumond MS
3. Hodel L
4. Vetter R
5. Iber D
(2021) openBIS
ID 3D Epithelium. 3D Epithelium.

https://openbis-data-repo.ethz.ch/openbis/webapp/eln-lims/?user=observer&pass=openbis

References

1. Aboav DA
(1970) The arrangement of grains in a polycrystal
Metallography 3:383–390.

https://doi.org/10.1016/0026-0800(70)90038-8
- Google Scholar
(2016) Segmentation and quantitative analysis of epithelial tissues
Methods in Molecular Biology 1478:227–239.

https://doi.org/10.1007/978-1-4939-6371-3_13
- PubMed
- Google Scholar
1. Barbier de Reuille P
2. Routier-Kierzkowska AL
3. Kierzkowski D
4. Bassel GW
5. Schüpbach T
6. Tauriello G
7. Bajpai N
8. Strauss S
9. Weber A
10. Kiss A
11. Burian A
12. Hofhuis H
13. Sapala A
14. Lipowczan M
15. Heimlicher MB
16. Robinson S
17. Bayer EM
18. Basler K
19. Koumoutsakos P
20. Roeder AH
21. Aegerter-Wilmsen T
22. Nakayama N
23. Tsiantis M
24. Hay A
25. Kwiatkowska D
26. Xenarios I
27. Kuhlemeier C
28. Smith RS
(2015) MorphoGraphX: A platform for quantifying morphogenesis in 4D
eLife 4:05864.

https://doi.org/10.7554/eLife.05864
- PubMed
- Google Scholar
(2005) Hexagonal packing of Drosophila wing epithelial cells by the planar cell polarity pathway
Developmental Cell 9:805–817.

https://doi.org/10.1016/j.devcel.2005.10.016
- PubMed
- Google Scholar
1. Conrad L
2. Runser SVM
3. Fernando Gómez H
4. Lang CM
5. Dumond MS
6. Sapala A
7. Schaumann L
8. Michos O
9. Vetter R
10. Iber D
(2021) The biomechanical basis of biased epithelial tube elongation in lung and kidney development
Development 148:dev194209.

https://doi.org/10.1242/dev.194209
- PubMed
- Google Scholar
1. Drubin DG
2. Nelson WJ
(1996) Origins of cell polarity
Cell 84:335–344.

https://doi.org/10.1016/S0092-8674(00)81278-7
- PubMed
- Google Scholar
1. Escudero LM
2. Costa LF
3. Kicheva A
4. Briscoe J
5. Freeman M
6. Babu MM
(2011) Epithelial organisation revealed by a network of cellular contacts
Nature Communications 2:526.

https://doi.org/10.1038/ncomms1536
- PubMed
- Google Scholar
1. Etournay R
2. Popović M
3. Merkel M
4. Nandi A
5. Blasse C
6. Aigouy B
7. Brandl H
8. Myers G
9. Salbreux G
10. Jülicher F
11. Eaton S
(2015) Interplay of cell dynamics and epithelial tension during morphogenesis of the Drosophila pupal wing
eLife 4:e07090.

https://doi.org/10.7554/eLife.07090
- PubMed
- Google Scholar
(2007) The influence of cell mechanics, cell-cell interactions, and proliferation on epithelial packing
Current Biology 17:2095–2104.

https://doi.org/10.1016/j.cub.2007.11.049
- PubMed
- Google Scholar
(2014) Vertex models of epithelial morphogenesis
Biophysical Journal 106:2291–2304.

https://doi.org/10.1016/j.bpj.2013.11.4498
- PubMed
- Google Scholar
(2006) The emergence of geometric order in proliferating metazoan epithelia
Nature 442:1038–1041.

https://doi.org/10.1038/nature05014
- PubMed
- Google Scholar
Software
1. Gómez H
2. Dumond M
3. Hodel L
4. Vetter R
5. Iber D
(2021) 2021_Gomez_3D_Cell_Neighbour_Dynamics, version d3dcc601
GitLab.

https://git.bsse.ethz.ch/iber/Publications/2021_gomez_3d_cell_neighbour_dynamics
(2018) Scutoids are a geometrical solution to three-dimensional packing of epithelia
Nature Communications 9:2960.

https://doi.org/10.1038/s41467-018-05376-1
- PubMed
- Google Scholar
1. Gritli-Linde A
2. Bei M
3. Maas R
4. Zhang XM
5. Linde A
6. McMahon AP
(2002) Shh signaling within the dental epithelium is necessary for cell proliferation, growth and polarization
Development 129:5323–5337.

https://doi.org/10.1242/dev.00100
- PubMed
- Google Scholar
1. Gundersen GG
2. Worman HJ
(2013) Nuclear positioning
Cell 152:1376–1389.

https://doi.org/10.1016/j.cell.2013.02.031
- PubMed
- Google Scholar
(2016) Extracellular Forces Cause the Nucleus to Deform in a Highly Controlled Anisotropic Manner
Scientific Reports 6:21300.

https://doi.org/10.1038/srep21300
- PubMed
- Google Scholar
1. Heller D
2. Hoppe A
3. Restrepo S
4. Gatti L
5. Tournier AL
6. Tapon N
7. Basler K
8. Mao Y
(2016) EpiTools: An Open-Source Image Analysis Toolkit for Quantifying Epithelial Growth Dynamics
Developmental Cell 36:103–116.

https://doi.org/10.1016/j.devcel.2015.12.012
- PubMed
- Google Scholar
Preprint
1. Hirashima T
2. Matsuda M
(2021) ERK-mediated curvature feedback regulates branching morphogenesis in lung epithelial tissue
bioRxiv.

https://doi.org/10.1101/2021.07.11.451982
- Google Scholar
1. Kadzik RS
2. Cohen ED
3. Morley MP
4. Stewart KM
5. Lu MM
6. Morrisey EE
(2014) Wnt ligand/Frizzled 2 receptor signaling regulates tube shape and branch-point formation in the lung through control of epithelial cell shape
PNAS 111:12444–12449.

https://doi.org/10.1073/pnas.1406639111
- PubMed
- Google Scholar
Preprint
(2019) Minimisation of surface energy drives apical epithelial organisation and gives rise to Lewis' law
bioRxiv.

https://doi.org/10.1101/590729
- Google Scholar
1. Kondo T
2. Hayashi S
(2015) Mechanisms of cell height changes that mediate epithelial invagination
Development, Growth & Differentiation 57:313–323.

https://doi.org/10.1111/dgd.12224
- PubMed
- Google Scholar
1. Lammerding J
(2011) Mechanics of the nucleus
Comprehensive Physiology 1:783–807.

https://doi.org/10.1002/cphy.c100038
- PubMed
- Google Scholar
1. Lewis FT
(1928) The correlation between cell division and the shapes and sizes of prismatic cells in the epidermis of cucumis
The Anatomical Record 38:341–376.

https://doi.org/10.1002/ar.1090380305
- Google Scholar
(2011) Interkinetic nuclear migration is a broadly conserved feature of cell division in pseudostratified epithelia
Current Biology 21:485–491.

https://doi.org/10.1016/j.cub.2011.02.002
- PubMed
- Google Scholar
1. Muzumdar MD
2. Tasic B
3. Miyamichi K
4. Li L
5. Luo L
(2007) A global double-fluorescent Cre reporter mouse
Genesis 45:593–605.

https://doi.org/10.1002/dvg.20335
- PubMed
- Google Scholar
1. Neelam S
2. Hayes PR
3. Zhang Q
4. Dickinson RB
5. Lele TP
(2016) Vertical uniformity of cells and nuclei in epithelial monolayers
Scientific Reports 6:19689.

https://doi.org/10.1038/srep19689
- PubMed
- Google Scholar
1. Norden C
(2017) Pseudostratified epithelia - cell biology, diversity and roles in organ formation at a glance
Journal of Cell Science 130:1859–1863.

https://doi.org/10.1242/jcs.192997
- PubMed
- Google Scholar
(2019) Cell-Size pleomorphism drives aberrant clone dispersal in proliferating epithelia
Developmental Cell 51:49–61.

https://doi.org/10.1016/j.devcel.2019.08.005
- PubMed
- Google Scholar
1. Rivier N
2. Lissowski A
(1982) On the correlation between sizes and shapes of cells in epithelial mosaics
Journal of Physics A: Mathematical and General 15:L143–L148.

https://doi.org/10.1088/0305-4470/15/3/012
- Google Scholar
1. Rodriguez-Boulan E
2. Macara IG
(2014) Organization and execution of the epithelial polarity programme
Nature Reviews Molecular Cell Biology 15:225–242.

https://doi.org/10.1038/nrm3775
- PubMed
- Google Scholar
1. Rupprecht JF
2. Ong KH
3. Yin J
4. Huang A
5. Dinh HH
6. Singh AP
7. Zhang S
8. Yu W
9. Saunders TE
(2017) Geometric constraints alter cell arrangements within curved epithelial tissues
Molecular Biology of the Cell 28:3582–3594.

https://doi.org/10.1091/mbc.e17-01-0060
- PubMed
- Google Scholar
(2016) Fundamental physical cellular constraints drive self-organization of tissues
The EMBO Journal 35:77–88.

https://doi.org/10.15252/embj.201592374
- PubMed
- Google Scholar
1. Schindelin J
2. Arganda-Carreras I
3. Frise E
4. Kaynig V
5. Longair M
6. Pietzsch T
7. Preibisch S
8. Rueden C
9. Saalfeld S
10. Schmid B
11. Tinevez JY
12. White DJ
13. Hartenstein V
14. Eliceiri K
15. Tomancak P
16. Cardona A
(2012) Fiji: an open-source platform for biological-image analysis
Nature Methods 9:676–682.

https://doi.org/10.1038/nmeth.2019
- PubMed
- Google Scholar
1. Shah P
2. Hobson CM
3. Cheng S
4. Colville MJ
5. Paszek MJ
6. Superfine R
7. Lammerding J
(2021) Nuclear deformation causes DNA damage by increasing replication stress
Current Biology 31:753–765.

https://doi.org/10.1016/j.cub.2020.11.037
- PubMed
- Google Scholar
1. Shin K
2. Margolis B
(2006) ZOning out tight junctions
Cell 126:647–649.

https://doi.org/10.1016/j.cell.2006.08.005
- PubMed
- Google Scholar
1. Susaki EA
2. Tainaka K
3. Perrin D
4. Yukinaga H
5. Kuno A
6. Ueda HR
(2015) Advanced CUBIC protocols for whole-brain and whole-body clearing and imaging
Nature Protocols 10:1709–1727.

https://doi.org/10.1038/nprot.2015.085
- PubMed
- Google Scholar
1. Tanaka S
(2015) Simulation frameworks for morphogenetic problems
Computation 3:197–221.

https://doi.org/10.3390/computation3020197
- Google Scholar
(2014) Quantitative imaging of cell dynamics in mouse embryos using light-sheet microscopy
Development 141:4406–4414.

https://doi.org/10.1242/dev.111021
- PubMed
- Google Scholar
Preprint
1. Vetter R
2. Kokic M
3. Gomez H
4. Hodel L
5. Gjeta B
6. Iannini A
7. Villa-Fombuena G
8. Casares F
9. Iber D
(2019) Aboave-Weaire's law in epithelia results from an angle constraint in contiguous polygonal lattices
bioRxiv.

https://doi.org/10.1101/591461
- Google Scholar
1. Widmann TJ
2. Dahmann C
(2009) Dpp signaling promotes the cuboidal-to-columnar shape transition of Drosophila wing disc epithelia by regulating Rho1
Journal of Cell Science 122:1362–1373.

https://doi.org/10.1242/jcs.044271
- PubMed
- Google Scholar

Article and author information

Author details

Harold Fernando Gómez
1. Department of Biosystems Science and Engineering (D-BSSE), ETH Zürich, Basel, Switzerland
2. Swiss Institute of Bioinformatics (SIB), Basel, Switzerland
Contribution
Data curation, Software, Formal analysis, Supervision, Validation, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-5200-4933
Mathilde Sabine Dumond
1. Department of Biosystems Science and Engineering (D-BSSE), ETH Zürich, Basel, Switzerland
2. Swiss Institute of Bioinformatics (SIB), Basel, Switzerland
Contribution
Data curation, Software, Formal analysis, Funding acquisition, Validation, Investigation, Visualization, Methodology

Competing interests
No competing interests declared
Leonie Hodel

Department of Biosystems Science and Engineering (D-BSSE), ETH Zürich, Basel, Switzerland

Contribution
Formal analysis, Investigation

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-2541-4133
Roman Vetter
1. Department of Biosystems Science and Engineering (D-BSSE), ETH Zürich, Basel, Switzerland
2. Swiss Institute of Bioinformatics (SIB), Basel, Switzerland
Contribution
Software, Formal analysis, Supervision, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-2901-7036
Dagmar Iber
1. Department of Biosystems Science and Engineering (D-BSSE), ETH Zürich, Basel, Switzerland
2. Swiss Institute of Bioinformatics (SIB), Basel, Switzerland
Contribution
Conceptualization, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Writing - original draft, Project administration, Writing - review and editing

For correspondence
dagmar.iber@bsse.ethz.ch

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-8051-1035

Funding

Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (CRSII5_170930)

Dagmar Iber

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work has been supported through an SNF Sinergia grant to DI. We acknowledge the many discussions and varied input of Marco Kokic, Anđela Markovic, Odyssé Michos, and Erand Smakaj, and we thank Richard Smith for help with MorphographX.

Ethics

Animal experimentation: Permission to use animals was obtained from the veterinary office of the Canton Basel-Stadt (license number 2777/26711). Experimental procedures were performed in accordance with the Guide for the Care and Use of Laboratory Animals and approved by the Ethics Committee for Animal Care of ETH Zurich.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.