Switch and template pattern formation in a discrete reaction-diffusion system inspired by the Drosophila eye

Patterning of the Drosophila retina by the morphogenetic furrow

Pattern formation is the process by which cells within a homogeneous epithelial sheet acquire distinctive fates depending upon their relative spatial position to each other. Several proposals, starting with Alan Turing's diffusion-reaction model, have been put forth over the last 70 years to describe how periodic patterns like those of vertebrate somites and skin hairs, mammalian molars, fish scales, and avian feather buds emerge during development. One of the best experimental systems for testing said models and identifying the gene regulatory networks that control pattern formation is the compound eye of the fruit fly, Drosophila melanogaster. Its cellular morphogenesis has been extensively studied for more than a century and hundreds of mutants that affect its development have been isolated. In this review we will focus on the morphogenetic furrow, a wave of differentiation that takes an initially homogeneous sheet of cells and converts it into an ordered array of unit eyes or ommatidia. Since the discovery of the furrow in 1976, positive and negative acting morphogens have been thought to be solely responsible for propagating the movement of the furrow across a motionless field of cells. However, a recent study has challenged this model and instead proposed that mechanical driven cell flow also contributes to retinal pattern formation. We will discuss both models and their impact on patterning.

Ordered hexagonal patterns via notch-delta signaling

Many developmental processes in biology utilize notch-delta signaling to construct an ordered pattern of cellular differentiation. This signaling modality is based on nearest-neighbor contact, as opposed to the more familiar mechanism driven by the release of diffusible ligands. Here, exploiting this 'juxtacrine' property, we present an exact treatment of the pattern formation problem via a system of nine coupled ordinary differential equations. The possible patterns that are realized for realistic parameters can be analyzed by considering a co-dimension 2 pitchfork bifurcation of this system. This analysis explains the observed prevalence of hexagonal patterns with high delta at their center, as opposed to those with central high notch levels (referred to as anti-hexagons). We show that outside this range of parameters, in particular for lowcis-coupling, a novel kind of pattern is produced, where high delta cells have high notch as well. It also suggests that the biological system is only weakly first order, so that an additional mechanism is required to generate the observed defect-free patterns. We construct a simple strategy for producing such defect-free patterns.

Mathematical modeling of Notch dynamics in Drosophila neural development

Notch signalling is a well-conserved signalling pathway that regulates cell fate through cell-cell communication. A typical feature of Notch signalling is ‘lateral inhibition’, whereby two neighbouring cells of equivalent state of differentiation acquire different cell fates. Recently, mathematical and computational approaches have addressed the Notch dynamics in Drosophila neural development. Typical examples of lateral inhibition are observed in the specification of neural stem cells in the embryo and sensory organ precursors in the thorax. In eye disc development, Notch signalling cooperates with other signalling pathways to define the evenly spaced positioning of the photoreceptor cells. The interplay between Notch and epidermal growth factor receptor signalling regulates the timing of neural stem cell differentiation in the optic lobe. In this review, we summarize the theoretical studies that have been conducted to elucidate the Notch dynamics in these systems and discuss the advantages of combining mathematical models with biological experiments.

Hexagonal patterning of the Drosophila eye

A complex network of transcription factor interactions propagates across the larval eye disc to establish columns of evenly-spaced R8 precursor cells, the founding cells of Drosophila ommatidia. After the recruitment of additional photoreceptors to each ommatidium, the surrounding cells are organized into their stereotypical pattern during pupal development. These support cells - comprised of pigment and cone cells - are patterned to encapsulate the photoreceptors and separate ommatidia with an hexagonal honeycomb lattice. Since the proteins and processes essential for correct eye patterning are conserved, elucidating how these function and change during Drosophila eye patterning can substantially advance our understanding of transcription factor and signaling networks, cytoskeletal structures, adhesion complexes, and the biophysical properties of complex tissues during their morphogenesis. Our understanding of many of these aspects of Drosophila eye patterning is largely descriptive. Many important questions, especially relating to the regulation and integration of cellular events, remain.

Defect patterns on the curved surface of fish retinae suggest a mechanism of cone mosaic formation

The outer epithelial layer of zebrafish retinae contains a crystalline array of cone photoreceptors, called the cone mosaic. As this mosaic grows by mitotic addition of new photoreceptors at the rim of the hemispheric retina, topological defects, called "Y-Junctions", form to maintain approximately constant cell spacing. The generation of topological defects due to growth on a curved surface is a distinct feature of the cone mosaic not seen in other well-studied biological patterns like the R8 photoreceptor array in the Drosophila compound eye. Since defects can provide insight into cell-cell interactions responsible for pattern formation, here we characterize the arrangement of cones in individual Y-Junction cores as well as the spatial distribution of Y-junctions across entire retinae. We find that for individual Y-junctions, the distribution of cones near the core corresponds closely to structures observed in physical crystals. In addition, Y-Junctions are organized into lines, called grain boundaries, from the retinal center to the periphery. In physical crystals, regardless of the initial distribution of defects, defects can coalesce into grain boundaries via the mobility of individual particles. By imaging in live fish, we demonstrate that grain boundaries in the cone mosaic instead appear during initial mosaic formation, without requiring defect motion. Motivated by this observation, we show that a computational model of repulsive cell-cell interactions generates a mosaic with grain boundaries. In contrast to paradigmatic models of fate specification in mostly motionless cell packings, this finding emphasizes the role of cell motion, guided by cell-cell interactions during differentiation, in forming biological crystals. Such a route to the formation of regular patterns may be especially valuable in situations, like growth on a curved surface, where the resulting long-ranged, elastic, effective interactions between defects can help to group them into grain boundaries.

Regulation of Proneural Wave Propagation Through a Combination of Notch-Mediated Lateral Inhibition and EGF-Mediated Reaction Diffusion

Notch-mediated lateral inhibition regulates binary cell fate choice, resulting in salt-and-pepper pattern formation during various biological processes. In many cases, Notch signaling acts together with other signaling systems. However, it is not clear what happens when Notch signaling is combined with other signaling systems. Mathematical modeling and the use of a simple biological model system will be essential to address this uncertainty. A wave of differentiation in the Drosophila visual center, the "proneural wave," accompanies the activity of the Notch and EGF signaling pathways. Although all of the Notch signaling components required for lateral inhibition are involved in the proneural wave, no salt-and-pepper pattern is found during the progression of the proneural wave. Instead, Notch is activated along the wave front and regulates proneural wave progression. How does Notch signaling control wave propagation without forming a salt-and-pepper pattern? A mathematical model of the proneural wave, based on biological evidence, has demonstrated that Notch-mediated lateral inhibition is implemented within the proneural wave and that the diffusible action of EGF cancels salt-and-pepper pattern formation. The results from numerical simulation have been confirmed by genetic experiments in vivo and suggest that the combination of Notch-mediated lateral inhibition and EGF-mediated reaction diffusion enables a novel function of Notch signaling that regulates propagation of the proneural wave. Similar mechanisms may play important roles in diverse biological processes found in animal development and cancer pathogenesis.

The proneural wave in the Drosophila optic lobe is driven by an excitable reaction-diffusion mechanism

In living organisms, self-organised waves of signalling activity propagate spatiotemporal information within tissues. During the development of the largest component of the visual processing centre of the Drosophila brain, a travelling wave of proneural gene expression initiates neurogenesis in the larval optic lobe primordium and drives the sequential transition of neuroepithelial cells into neuroblasts. Here, we propose that this 'proneural wave' is driven by an excitable reaction-diffusion system involving epidermal growth factor receptor (EGFR) signalling interacting with the proneural gene l'sc. Within this framework, a propagating transition zone emerges from molecular feedback and diffusion. Ectopic activation of EGFR signalling in clones within the neuroepithelium demonstrates that a transition wave can be excited anywhere in the tissue by inducing signalling activity, consistent with a key prediction of the model. Our model illuminates the physical and molecular underpinnings of proneural wave progression and suggests a generic mechanism for regulating the sequential differentiation of tissues.

Notch-mediated lateral inhibition regulates proneural wave propagation when combined with EGF-mediated reaction diffusion

Notch-mediated lateral inhibition regulates binary cell fate choice, resulting in salt and pepper patterns during various developmental processes. However, how Notch signaling behaves in combination with other signaling systems remains elusive. The wave of differentiation in the Drosophila visual center or "proneural wave" accompanies Notch activity that is propagated without the formation of a salt and pepper pattern, implying that Notch does not form a feedback loop of lateral inhibition during this process. However, mathematical modeling and genetic analysis clearly showed that Notch-mediated lateral inhibition is implemented within the proneural wave. Because partial reduction in EGF signaling causes the formation of the salt and pepper pattern, it is most likely that EGF diffusion cancels salt and pepper pattern formation in silico and in vivo. Moreover, the combination of Notch-mediated lateral inhibition and EGF-mediated reaction diffusion enables a function of Notch signaling that regulates propagation of the wave of differentiation.

A two-step patterning process increases the robustness of periodic patterning in the fly eye

Complex periodic patterns can self-organize through dynamic interactions between diffusible activators and inhibitors. In the biological context, self-organized patterning is challenged by spatial heterogeneities ('noise') inherent to biological systems. How spatial variability impacts the periodic patterning mechanism and how it can be buffered to ensure precise patterning is not well understood. We examine the effect of spatial heterogeneity on the periodic patterning of the fruit fly eye, an organ composed of ∼800 miniature eye units (ommatidia) whose periodic arrangement along a hexagonal lattice self-organizes during early stages of fly development. The patterning follows a two-step process, with an initial formation of evenly spaced clusters of ∼10 cells followed by a subsequent refinement of each cluster into a single selected cell. Using a probabilistic approach, we calculate the rate of patterning errors resulting from spatial heterogeneities in cell size, position and biosynthetic capacity. Notably, error rates were largely independent of the desired cluster size but followed the distributions of signaling speeds. Pre-formation of large clusters therefore greatly increases the reproducibility of the overall periodic arrangement, suggesting that the two-stage patterning process functions to guard the pattern against errors caused by spatial heterogeneities. Our results emphasize the constraints imposed on self-organized patterning mechanisms by the need to buffer stochastic effects. Author summary Complex periodic patterns are common in nature and are observed in physical, chemical and biological systems. Understanding how these patterns are generated in a precise manner is a key challenge. Biological patterns are especially intriguing, as they are generated in a noisy environment; cell position and cell size, for example, are subject to stochastic variations, as are the strengths of the chemical signals mediating cell-to-cell communication. The need to generate a precise and robust pattern in this 'noisy' environment restricts the space of patterning mechanisms that can function in the biological setting. Mathematical modeling is useful in comparing the sensitivity of different mechanisms to such variations, thereby highlighting key aspects of their design.We use mathematical modeling to study the periodic patterning of the fruit fly eye. In this system, a highly ordered lattice of differentiated cells is generated in a two-dimensional cell epithelium. The pattern is first observed by the appearance of evenly spaced clusters of ∼10 cells that express specific genes. Each cluster is subsequently refined into a single cell, which initiates the formation and differentiation of a miniature eye unit, the ommatidium. We formulate a mathematical model based on the known molecular properties of the patterning mechanism, and use a probabilistic approach to calculate the errors in cluster formation and refinement resulting from stochastic cell-to-cell variations ('noise') in different quantitative parameters. This enables us to define the parameters most influencing noise sensitivity. Notably, we find that this error is roughly independent of the desired cluster size, suggesting that large clusters are beneficial for ensuring the overall reproducibility of the periodic cluster arrangement. For the stage of cluster refinement, we find that rapid communication between cells is critical for reducing error. Our work provides new insights into the constraints imposed on mechanisms generating periodic patterning in a realistic, noisy environment, and in particular, discusses the different considerations in achieving optimal design of the patterning network.

Periodic patterning of the Drosophila eye is stabilized by the diffusible activator Scabrous

Generation of periodic patterns is fundamental to the differentiation of multiple tissues during development. How such patterns form robustly is still unclear. The Drosophila eye comprises ∼750 units, whose crystalline order is set during differentiation of the eye imaginal disc: an activation wave sweeping across the disc is coupled to lateral inhibition, sequentially selecting pro-neural cells. Using mathematical modelling, here we show that this template-based lateral inhibition is highly sensitive to spatial variations in biochemical parameters and cell sizes. We reveal the basis of this sensitivity, and suggest that it can be overcome by assuming a short-range diffusible activator. Clonal experiments identify Scabrous, a previously implicated inhibitor, as the predicted activator. Our results reveal the mechanism by which periodic patterning in the fly eye is stabilized against spatial variations, highlighting how the need to maintain robustness shapes the design of patterning circuits.

Patterning in the Drosophila eye is achieved by a series of signalling cascades over several cell distances. Here Gavish et al. model lateral inhibition in the developing eye to understand how developmental noise refines such patterning, identifying a novel activator required for buffering spatial variability.

Modeling Notch signaling: a practical tutorial

Theoretical and computational approaches for understanding different aspects of Notch signaling and Notch dependent patterning are gaining popularity in recent years. These in silico methodologies can provide dynamic insights that are often not intuitive and may help guide experiments aimed at elucidating these processes. This chapter is an introductory tutorial intended to allow someone with basic mathematical and computational knowledge to explore new mathematical models of Notch-mediated processes and perform numerical simulations of these models. In particular, we explain how to define and simulate models of lateral inhibition patterning processes. We provide a Matlab code for simulating various lateral inhibition models in a simple and intuitive manner, and show how to present the results from the computational models. This code can be used as a starting point for exploring more specific models that include additional aspects of the Notch pathway and its regulation.

Regulation of neuronal differentiation at the neurogenic wavefront

Signaling mediated by the Delta/Notch system controls the process of lateral inhibition, known to regulate neurogenesis in metazoans. Lateral inhibition takes place in equivalence groups formed by cells having equal capacity to differentiate, and it results in the singling out of precursors, which subsequently become neurons. During normal development, areas of active neurogenesis spread through non-neurogenic regions in response to specific morphogens, giving rise to neurogenic wavefronts. Close contact of these wavefronts with non-neurogenic cells is expected to affect lateral inhibition. Therefore, a mechanism should exist in these regions to prevent disturbances of the lateral inhibitory process. Focusing on the developing chick retina, we show that Dll1 is widely expressed by non-neurogenic precursors located at the periphery of this tissue, a region lacking Notch1, lFng, and differentiation-related gene expression. We investigated the role of this Dll1 expression through mathematical modeling. Our analysis predicts that the absence of Dll1 ahead of the neurogenic wavefront results in reduced robustness of the lateral inhibition process, often linked to enhanced neurogenesis and the presence of morphological alterations of the wavefront itself. These predictions are consistent with previous observations in the retina of mice in which Dll1 is conditionally mutated. The predictive capacity of our mathematical model was confirmed further by mimicking published results on the perturbation of morphogenetic furrow progression in the eye imaginal disc of Drosophila. Altogether, we propose that Notch-independent Delta expression ahead of the neurogenic wavefront is required to avoid perturbations in lateral inhibition and wavefront progression, thus optimizing the neurogenic process.

A dynamical model of ommatidial crystal formation

The crystalline photoreceptor lattice in the Drosophila eye is a paradigm for pattern formation during development. During eye development, activation of proneural genes at a moving front adds new columns to a regular lattice of R8 photoreceptors. We present a mathematical model of the governing activator-inhibitor system, which indicates that the dynamics of positive induction play a central role in the selection of certain cells as R8s. The "switch and template" patterning mechanism we observe is mathematically very different from the well-known Turing instability. Unlike a standard lateral inhibition model, our picture implies that R8s are defined before the appearance of the complete group of proneural cells. The model reproduces the full time course of proneural gene expression and accounts for specific features of the refinement of proneural groups that had resisted explanation. It moreover predicts that perturbing the normal template can lead to eyes containing stripes of R8 cells. We observed these stripes experimentally after manipulation of the Notch and scabrous genes. Our results suggest an alternative to the generally assumed mode of operation for lateral inhibition during development; more generally, they hint at a broader role for bistable switches in the initial establishment of patterns as well as in their maintenance.

Modeling bistable cell-fate choices in the Drosophila eye: qualitative and quantitative perspectives

A major goal of developmental biology is to understand the molecular mechanisms whereby genetic signaling networks establish and maintain distinct cell types within multicellular organisms. Here, we review cell-fate decisions in the developing eye of Drosophila melanogaster and the experimental results that have revealed the topology of the underlying signaling circuitries. We then propose that switch-like network motifs based on positive feedback play a central role in cell-fate choice, and discuss how mathematical modeling can be used to understand and predict the bistable or multistable behavior of such networks.