Dynamic scaling of morphogen gradients on growing domains

Brownian non-Gaussian polymer diffusion in non-static media

In nature, essentially, almost all the particles move irregularly in non-static media. With the advance of observation techniques, various kinds of new dynamical phenomena are detected, e.g., Brownian non-Gaussian diffusion. This paper focuses on the dynamical behavior of the center of mass (CM) of a polymer in non-static media and investigates the effect of polymer size fluctuations on the diffusion behavior. First, we establish a diffusing diffusivity model for polymer size fluctuations, linking the polymer size variation to the birth and death process, and introduce co-moving and physical coordinate systems to characterize the position of the CM for a polymer in non-static media. Next, the important statistical quantities for the CM diffusing diffusivity model in non-static media, such as mean square displacement (MSD) and kurtosis, are obtained by adopting the subordinate process approach, and the long-time asymptotic behavior of the MSD in the media of different types is specifically analyzed. Finally, the bivariate Fokker-Planck equation and the Feynman-Kac equation corresponding to the diffusing diffusivity model are detailedly derived and solved through the deep backward stochastic differential equation (BSDE) method to confirm the correctness of the derived equations.

Morphogens in the evolution of size, shape and patterning

Much of the striking diversity of life on Earth has arisen from variations in the way that the same molecules and networks operate during development to shape and pattern tissues and organs into different morphologies. However, we still understand very little about the potential for diversification exhibited by different, highly conserved mechanisms during evolution, or, conversely, the constraints that they place on evolution. With the aim of steering the field in new directions, we focus on morphogen-mediated patterning and growth as a case study to demonstrate how conserved developmental mechanisms can adapt during evolution to drive morphological diversification and optimise functionality, and to illustrate how evolution algorithms and computational tools can be used alongside experiments to provide insights into how these conserved mechanisms can evolve. We first introduce key conserved properties of morphogen-driven patterning mechanisms, before summarising comparative studies that exemplify how changes in the spatiotemporal expression and signalling levels of morphogens impact the diversification of organ size, shape and patterning in nature. Finally, we detail how theoretical frameworks can be used in conjunction with experiments to probe the role of morphogen-driven patterning mechanisms in evolution. We conclude that morphogen-mediated patterning is an excellent model system and offers a generally applicable framework to investigate the evolution of developmental mechanisms.

'Three signals - three body axes' as patterning principle in bilaterians

In vertebrates, the three orthogonal body axes, anteroposterior (AP), dorsoventral (DV) and left-right (LR) are determined at gastrula and neurula stages by the Spemann-Mangold organizer and its equivalents. A common feature of AP and DV axis formation is that an evolutionary conserved interplay between growth factors (Wnt, BMP) and their extracellular antagonists (e.g. Dkk1, Chordin) creates signaling gradients for axial patterning. Recent work showed that LR patterning in Xenopus follows the same principle, with R-spondin 2 (Rspo2) as an extracellular FGF antagonist, which creates a signaling gradient that determines the LR vector. That a triad of anti-FGF, anti-BMP, and anti-Wnt governs LR, DV, and AP axis formation reveals a unifying principle in animal development. We discuss how cross-talk between these three signals confers integrated AP-DV-LR body axis patterning underlying developmental robustness, size scaling, and harmonious regulation. We propose that Urbilateria featured three orthogonal body axes that were governed by a Cartesian coordinate system of orthogonal Wnt/AP, BMP/DV, and FGF/LR signaling gradients.

2D effects enhance precision of gradient-based tissue patterning

Robust embryonic development requires pattern formation with high spatial accuracy. In epithelial tissues that are patterned by morphogen gradients, the emerging patterns achieve levels of precision that have recently been explained by a simple one-dimensional reaction-diffusion model with kinetic noise. Here, we show that patterning precision is even greater if transverse diffusion effects are at play in such tissues. The positional error, a measure for spatial patterning accuracy, decreases in wider tissues but then saturates beyond a width of about ten cells. This demonstrates that the precision of gradient-based patterning in two- or higher-dimensional systems can be even greater than predicted by 1D models, and further attests to the potential of noisy morphogen gradients for high-precision tissue patterning.

A multiscale chemical-mechanical model predicts impact of morphogen spreading on tissue growth

The exact mechanism controlling cell growth remains a grand challenge in developmental biology and regenerative medicine. The Drosophila wing disc tissue serves as an ideal biological model to study mechanisms involved in growth regulation. Most existing computational models for studying tissue growth focus specifically on either chemical signals or mechanical forces. Here we developed a multiscale chemical-mechanical model to investigate the growth regulation mechanism based on the dynamics of a morphogen gradient. By comparing the spatial distribution of dividing cells and the overall tissue shape obtained in model simulations with experimental data of the wing disc, it is shown that the size of the domain of the Dpp morphogen is critical in determining tissue size and shape. A larger tissue size with a faster growth rate and more symmetric shape can be achieved if the Dpp gradient spreads in a larger domain. Together with Dpp absorbance at the peripheral zone, the feedback regulation that downregulates Dpp receptors on the cell membrane allows for further spreading of the morphogen away from its source region, resulting in prolonged tissue growth at a more spatially homogeneous growth rate.

The impact of cell size on morphogen gradient precision

Tissue patterning during embryonic development is remarkably precise. Here, we numerically determine the impact of the cell diameter, gradient length and the morphogen source on the variability of morphogen gradients. We show that the positional error increases with the gradient length relative to the size of the morphogen source, and with the square root of the cell diameter and the readout position. We provide theoretical explanations for these relationships, and show that they enable high patterning precision over developmental time for readouts that scale with expanding tissue domains, as observed in the Drosophila wing disc. Our analysis suggests that epithelial tissues generally achieve higher patterning precision with small cross-sectional cell areas. An extensive survey of measured apical cell areas shows that they are indeed small in developing tissues that are patterned by morphogen gradients. Enhanced precision may thus have led to the emergence of pseudostratification in epithelia, a phenomenon for which the evolutionary benefit had so far remained elusive.

Coupling cell proliferation rates to the duration of recruitment controls final size of the Drosophila wing

Organ growth driven by cell proliferation is an exponential process. As a result, even small variations in proliferation rates, when integrated over a relatively long developmental time, will lead to large differences in size. How organs robustly control their final size despite perturbations in cell proliferation rates throughout development is a long-standing question in biology. Using a mathematical model, we show that in the developing wing of the fruit fly, Drosophila melanogaster, variations in proliferation rates of wing-committed cells are inversely proportional to the duration of cell recruitment, a differentiation process in which a population of undifferentiated cells adopt the wing fate by expressing the selector gene, vestigial. A time-course experiment shows that vestigial-expressing cells increase exponentially while recruitment takes place, but slows down when recruitable cells start to vanish, suggesting that undifferentiated cells may be driving proliferation of wing-committed cells. When this observation is incorporated in our model, we show that the duration of cell recruitment robustly determines a final wing size even when cell proliferation rates of wing-committed cells are perturbed. Finally, we show that this control mechanism fails when perturbations in proliferation rates affect both wing-committed and recruitable cells, providing an experimentally testable hypothesis of our model.

Morphogen gradient scaling by recycling of intracellular Dpp

Morphogen gradients are fundamental to establish morphological patterns in developing tissues¹. During development, gradients scale to remain proportional to the size of growing organs^2,3. Scaling is a universal gear adjusting patterns to size in living organisms^3-8, yet its mechanisms remain unclear. Here, focusing on the Dpp gradient in the Drosophila wing disc, we unravel a cell biological basis behind scaling. From small to large discs, scaling of the Dpp gradient is achieved by increasing the contribution of the internalized Dpp molecules to Dpp transport: to expand the gradient, endocytosed molecules are re-exocytosed to spread extracellularly. To regulate the contribution of endocytosed Dpp to the spreading extracellular pool during tissue growth, it is the Dpp binding rates that are progressively modulated by the extracellular factor Pentagone, driving scaling. Thus, for some morphogens, evolution may act on endocytic trafficking to regulate the range of the gradient and its scaling, which could allow adaptation of shape and pattern to different sizes of organs in different species.

Scalable control of developmental timetables by epigenetic switching networks

During development, progenitor cells follow timetables for differentiation that span many cell generations. These developmental timetables are robustly encoded by the embryo, yet scalably adjustable by evolution, facilitating variation in organism size and form. Epigenetic switches, involving rate-limiting activation steps at regulatory gene loci, control gene activation timing in diverse contexts, and could profoundly impact the dynamics of gene regulatory networks controlling developmental lineage specification. Here, we develop a mathematical framework to model regulatory networks with genes controlled by epigenetic switches. Using this framework, we show that such epigenetic switching networks uphold developmental timetables that robustly span many cell generations, and enable the generation of differentiated cells in precisely defined numbers and fractions. Changes to epigenetic switching networks can readily alter the timing of developmental events within a timetable, or alter the overall speed at which timetables unfold, enabling scalable control over differentiated population sizes. With their robust, yet flexibly adjustable nature, epigenetic switching networks could represent central targets on which evolution acts to manufacture diversity in organism size and form.

Multiple morphogens and rapid elongation promote segmental patterning during development

The vertebrate hindbrain is segmented into rhombomeres (r) initially defined by distinct domains of gene expression. Previous studies have shown that noise-induced gene regulation and cell sorting are critical for the sharpening of rhombomere boundaries, which start out rough in the forming neural plate (NP) and sharpen over time. However, the mechanisms controlling simultaneous formation of multiple rhombomeres and accuracy in their sizes are unclear. We have developed a stochastic multiscale cell-based model that explicitly incorporates dynamic morphogenetic changes (i.e. convergent-extension of the NP), multiple morphogens, and gene regulatory networks to investigate the formation of rhombomeres and their corresponding boundaries in the zebrafish hindbrain. During pattern initiation, the short-range signal, fibroblast growth factor (FGF), works together with the longer-range morphogen, retinoic acid (RA), to specify all of these boundaries and maintain accurately sized segments with sharp boundaries. At later stages of patterning, we show a nonlinear change in the shape of rhombomeres with rapid left-right narrowing of the NP followed by slower dynamics. Rapid initial convergence improves boundary sharpness and segment size by regulating cell sorting and cell fate both independently and coordinately. Overall, multiple morphogens and tissue dynamics synergize to regulate the sizes and boundaries of multiple segments during development.

In segmental pattern formation, chemical gradients control gene expression in a concentration-dependent manner to specify distinct gene expression domains. Despite the stochasticity inherent to such biological processes, precise and accurate borders form between segmental gene expression domains. Previous work has revealed synergy between gene regulation and cell sorting in sharpening borders that are initially rough. However, it is still poorly understood how size and boundary sharpness of multiple segments are regulated in a tissue that changes dramatically in its morphology as the embryo develops. Here we develop a stochastic multiscale cell-base model to investigate these questions. Two novel strategies synergize to promote accurate segment formation, a combination of long- and short-range morphogens plus rapid tissue convergence, with one responsible for pattern initiation and the other enabling pattern refinement.

Diversity and robustness of bone morphogenetic protein pattern formation

Pattern formation by bone morphogenetic proteins (BMPs) demonstrates remarkable plasticity and utility in several contexts, such as early embryonic development, tissue patterning and the maintenance of stem cell niches. BMPs pattern tissues over many temporal and spatial scales: BMP gradients as short as 1-2 cell diameters maintain the stem cell niche of the Drosophila germarium over a 24-h cycle, and BMP gradients of several hundred microns establish dorsal-ventral tissue specification in Drosophila, zebrafish and Xenopus embryos in timescales between 30 min and several hours. The mechanisms that shape BMP signaling gradients are also incredibly diverse. Although ligand diffusion plays a dominant role in forming the gradient, a cast of diffusible and non-diffusible regulators modulate gradient formation and confer robustness, including scale invariance and adaptability to perturbations in gene expression and growth. In this Review, we document the diverse ways that BMP gradients are formed and refined, and we identify the core principles that they share to achieve reliable performance.

Scaling a Dpp Morphogen Gradient through Feedback Control of Receptors and Co-receptors

Gradients of decapentaplegic (Dpp) pattern Drosophila wing imaginal discs, establishing gene expression boundaries at specific locations. As discs grow, Dpp gradients expand, keeping relative boundary positions approximately stationary. Such scaling fails in mutants for Pentagone (pent), a gene repressed by Dpp that encodes a diffusible protein that expands Dpp gradients. Although these properties fit a recent mathematical model of automatic gradient scaling, that model requires an expander that spreads with minimal loss throughout a morphogen field. Here, we show that Pent's actions are confined to within just a few cell diameters of its site of synthesis and can be phenocopied by manipulating non-diffusible Pent targets strictly within the Pent expression domain. Using genetics and mathematical modeling, we develop an alternative model of scaling driven by feedback downregulation of Dpp receptors and co-receptors. Among the model's predictions is a size beyond which scaling fails-something we observe directly in wing discs.

Reaction-diffusion and reaction-subdiffusion equations on arbitrarily evolving domains

Reaction-diffusion equations are widely used as the governing evolution equations for modeling many physical, chemical, and biological processes. Here we derive reaction-diffusion equations to model transport with reactions on a one-dimensional domain that is evolving. The model equations, which have been derived from generalized continuous time random walks, can incorporate complexities such as subdiffusive transport and inhomogeneous domain stretching and shrinking. Inhomogeneously growing domains are frequently encountered in biological phenomena involving stochastic transport, such as tumor growth and morphogen gradient formation. A method for constructing analytic expressions for short-time moments of the position of the particles is developed and moments calculated from this approach are shown to compare favorably with results from random walk simulations and numerical integration of the reaction transport equation. The results show the important role played by the initial condition. In particular, it strongly affects the time dependence of the moments in the short-time regime by introducing additional drift and diffusion terms. We also discuss how our reaction transport equation could be applied to study the spreading of a population on an evolving interface. From a more general perspective, our findings help to mitigate the scarcity of analytic results for reaction-diffusion problems in geometries displaying nonuniform growth. They are also expected to pave the way for further results, including the treatment of first-passage problems associated with encounter-controlled reactions in such domains.

Modeling of Wnt-mediated tissue patterning in vertebrate embryogenesis

During embryogenesis, morphogens form a concentration gradient in responsive tissue, which is then translated into a spatial cellular pattern. The mechanisms by which morphogens spread through a tissue to establish such a morphogenetic field remain elusive. Here, we investigate by mutually complementary simulations and in vivo experiments how Wnt morphogen transport by cytonemes differs from typically assumed diffusion-based transport for patterning of highly dynamic tissue such as the neural plate in zebrafish. Stochasticity strongly influences fate acquisition at the single cell level and results in fluctuating boundaries between pattern regions. Stable patterning can be achieved by sorting through concentration dependent cell migration and apoptosis, independent of the morphogen transport mechanism. We show that Wnt transport by cytonemes achieves distinct Wnt thresholds for the brain primordia earlier compared with diffusion-based transport. We conclude that a cytoneme-mediated morphogen transport together with directed cell sorting is a potentially favored mechanism to establish morphogen gradients in rapidly expanding developmental systems.

Regulation of size and scale in vertebrate spinal cord development

All vertebrates have a spinal cord with dimensions and shape specific to their species. Yet how species‐specific organ size and shape are achieved is a fundamental unresolved question in biology. The formation and sculpting of organs begins during embryonic development. As it develops, the spinal cord extends in anterior–posterior direction in synchrony with the overall growth of the body. The dorsoventral (DV) and apicobasal lengths of the spinal cord neuroepithelium also change, while at the same time a characteristic pattern of neural progenitor subtypes along the DV axis is established and elaborated. At the basis of these changes in tissue size and shape are biophysical determinants, such as the change in cell number, cell size and shape, and anisotropic tissue growth. These processes are controlled by global tissue‐scale regulators, such as morphogen signaling gradients as well as mechanical forces. Current challenges in the field are to uncover how these tissue‐scale regulatory mechanisms are translated to the cellular and molecular level, and how regulation of distinct cellular processes gives rise to an overall defined size. Addressing these questions will help not only to achieve a better understanding of how size is controlled, but also of how tissue size is coordinated with the specification of pattern.

This article is categorized under:

Establishment of Spatial and Temporal Patterns > Regulation of Size, Proportion, and Timing

Signaling Pathways > Global Signaling Mechanisms

Nervous System Development > Vertebrates: General Principles

Evaluation of BMP-mediated patterning in a 3D mathematical model of the zebrafish blastula embryo

Bone Morphogenetic Proteins (BMPs) play an important role in dorsal-ventral (DV) patterning of the early zebrafish embryo. BMP signaling is regulated by a network of extracellular and intracellular factors that impact the range and signaling of BMP ligands. Recent advances in understanding the mechanism of pattern formation support a source-sink mechanism, however it is not clear how the source-sink mechanism shapes patterns in 3D, nor how sensitive the pattern is to biophysical rates and boundary conditions along both the anteroposterior (AP) and DV axes of the embryo. We propose a new three-dimensional growing Partial Differential Equation (PDE)-based model to simulate the BMP patterning process during the blastula stage. This model provides a starting point to elucidate how different mechanisms and components work together in 3D to create and maintain the BMP gradient in the embryo. We also show how the 3D model fits the BMP signaling gradient data at multiple time points along both axes. Furthermore, sensitivity analysis of the model suggests that the spatiotemporal patterns of Chordin and BMP ligand gene expression are dominant drivers of shape in 3D and more work is needed to quantify the spatiotemporal profiles of gene and protein expression to further refine the models.

Positional information and tissue scaling during development and regeneration

In order to contribute to the appropriate tissues during development, cells need to know their position within the embryo. This positional information is conveyed by gradients of signaling molecules, termed morphogens, that are produced in specific regions of the embryo and induce concentration-dependent responses in target tissues. Positional information is remarkably robust, and embryos often develop with the correct proportions even if large parts of the embryo are removed. In this Review, we discuss classical embryological experiments and modern quantitative analyses that have led to mechanistic insights into how morphogen gradients adapt, scale and properly pattern differently sized domains. We analyze these experimental findings in the context of mathematical models and synthesize general principles that apply to multiple systems across species and developmental stages.

Standard and fractional Ornstein-Uhlenbeck process on a growing domain

We study normal diffusive and subdiffusive processes in a harmonic potential (Ornstein-Uhlenbeck process) on a uniformly growing or contracting domain. Our starting point is a recently derived fractional Fokker-Planck equation, which covers both the case of Brownian diffusion and the case of a subdiffusive continuous-time random walk (CTRW). We find a high sensitivity of the random walk properties to the details of the domain growth rate, which gives rise to a variety of regimes with extremely different behaviors. At the origin of this rich phenomenology is the fact that the walkers still move while they wait to jump, since they are dragged by the deterministic drift arising from the domain growth. Thus, the increasingly long waiting times associated with the aging of the subdiffusive CTRW imply that, in the time interval between two consecutive jumps, the walkers might travel over much longer distances than in the normal diffusive case. This gives rise to seemingly counterintuitive effects. For example, on a static domain, both Brownian diffusion and subdiffusive CTRWs yield a stationary particle distribution with finite width when a harmonic potential is at play, thus indicating a confinement of the diffusing particle. However, for a sufficiently fast growing or contracting domain, this qualitative behavior breaks down, and differences between the Brownian case and the subdiffusive case are found. In the case of Brownian particles, a sufficiently fast exponential domain growth is needed to break the confinement induced by the harmonic force; in contrast, for subdiffusive particles such a breakdown may already take place for a sufficiently fast power-law domain growth. Our analytic and numerical results for both types of diffusion are fully confirmed by random walk simulations.

Scale invariance of BMP signaling gradients in zebrafish

In both vertebrates and invertebrates, spatial patterning along the Dorsal-ventral (DV) embryonic axis depends on a morphogen gradient of Bone Morphogenetic Protein signaling. Scale invariance of DV patterning by BMPs has been found in both vertebrates and invertebrates, however the mechanisms that regulate gradient scaling remain controversial. To obtain quantitative data that can be used to address core questions of scaling, we introduce a method to tune the size of zebrafish embryos by reducing varying amounts of vegetal yolk. We quantified the BMP signaling gradient in wild-type and perturbed embryos and found that the system scales for reductions in cross-sectional perimeter of up to 30%. Furthermore, we found that the degree of scaling for intraspecies scaling within zebrafish is greater than that between Danioninae species.

Growth and size control during development

The size and shape of organs are characteristic for each species. Even when organisms develop to different sizes due to varying environmental conditions, such as nutrition, organ size follows species-specific rules of proportionality to the rest of the body, a phenomenon referred to as allometry. Therefore, for a given environment, organs stop growth at a predictable size set by the species's genotype. How do organs stop growth? How can related species give rise to organs of strikingly different size? No definitive answer has been given to date. One of the major models for the studies of growth termination is the vinegar fly Drosophila melanogaster. Therefore, this review will focus mostly on work carried out in Drosophila to try to tease apart potential mechanisms and identify routes for further investigation. One general rule, found across the animal kingdom, is that the rate of growth declines with developmental time. Therefore, answers to the problem of growth termination should explain this seemingly universal fact. In addition, growth termination is intimately related to the problems of robustness (i.e. precision) and plasticity in organ size, symmetric and asymmetric organ development, and of how the ‘target’ size depends on extrinsic, environmental factors.

Critical Point in Self-Organized Tissue Growth

We present a theory of pattern formation in growing domains inspired by biological examples of tissue development. Gradients of signaling molecules regulate growth, while growth changes these graded chemical patterns by dilution and advection. We identify a critical point of this feedback dynamics, which is characterized by spatially homogeneous growth and proportional scaling of patterns with tissue length. We apply this theory to the biological model system of the developing wing of the fruit fly Drosophila melanogaster and quantitatively identify signatures of the critical point.

Simulation of Morphogen and Tissue Dynamics

Morphogenesis, the process by which an adult organism emerges from a single cell, has fascinated humans for a long time. Modeling this process can provide novel insights into development and the principles that orchestrate the developmental processes. This chapter focuses on the mathematical description and numerical simulation of developmental processes. In particular, we discuss the mathematical representation of morphogen and tissue dynamics on static and growing domains, as well as the corresponding tissue mechanics. In addition, we give an overview of numerical methods that are routinely used to solve the resulting systems of partial differential equations. These include the finite element method and the Lattice Boltzmann method for the discretization as well as the arbitrary Lagrangian-Eulerian method and the Diffuse-Domain method to numerically treat deforming domains.

Mechanisms and Measurements of Scale Invariance of Morphogen Gradients

Morphogen gradients provide positional information to underlying cells that translate the information into differential gene expression and eventually different cell fates. Scale invariance is the property where the gradients of the morphogen adjust proportionately to the size of the domain. Scale invariance of morphogen gradients or patterns of differentiation is a common phenomenon observed between individuals within the same species and between homologous tissues or structures in different species. To determine whether or not a pattern is scale invariant, others and we have developed definitions and measurements of gradient scaling. These include point-wise and global scaling errors as well as global scaling power. Furthermore, there are a number of mathematical conditions for scale invariance of advection-diffusion-reaction models that inform mechanisms of scaling. Herein we provide a deeper perspective on modeling and measurement of scale invariance of morphogen gradients.

On the importance of protein diffusion in biological systems: The example of the Bicoid morphogen gradient

Morphogens are proteins that form concentration gradients in embryos and developing tissues, where they act as postal codes, providing cells with positional information and allowing them to behave accordingly. Bicoid was the first discovered morphogen, and remains one of the most studied. It regulates segmentation in flies, forming a striking exponential gradient along the anterior-posterior axis of early Drosophila embryos, and activating the transcription of multiple target genes in a concentration-dependent manner. In this review, the work done by us and by others to characterize the mobility of Bicoid in D. melanogaster embryos is presented. The central role played by the diffusion of Bicoid in both the establishment of the gradient and the activation of target genes is discussed, and placed in the context of the need for these processes to be all at once rapid, precise and robust. The Bicoid system, and morphogen gradients in general, remain amongst the most amazing examples of the coexistence, often observed in living systems, of small-scale disorder and large-scale spatial order. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Continuous-time random-walk model for anomalous diffusion in expanding media

Expanding media are typical in many different fields, e.g., in biology and cosmology. In general, a medium expansion (contraction) brings about dramatic changes in the behavior of diffusive transport properties such as the set of positional moments and the Green's function. Here, we focus on the characterization of such effects when the diffusion process is described by the continuous-time random-walk (CTRW) model. As is well known, when the medium is static this model yields anomalous diffusion for a proper choice of the probability density function (pdf) for the jump length and the waiting time, but the behavior may change drastically if a medium expansion is superimposed on the intrinsic random motion of the diffusing particle. For the case where the jump length and the waiting time pdfs are long-tailed, we derive a general bifractional diffusion equation which reduces to a normal diffusion equation in the appropriate limit. We then study some particular cases of interest, including Lévy flights and subdiffusive CTRWs. In the former case, we find an analytical exact solution for the Green's function (propagator). When the expansion is sufficiently fast, the contribution of the diffusive transport becomes irrelevant at long times and the propagator tends to a stationary profile in the comoving reference frame. In contrast, for a contracting medium a competition between the spreading effect of diffusion and the concentrating effect of contraction arises. In the specific case of a subdiffusive CTRW in an exponentially contracting medium, the latter effect prevails for sufficiently long times, and all the particles are eventually localized at a single point in physical space. This "big crunch" effect, totally absent in the case of normal diffusion, stems from inefficient particle spreading due to subdiffusion. We also derive a hierarchy of differential equations for the moments of the transport process described by the subdiffusive CTRW model in an expanding medium. From this hierarchy, the full time evolution of the second-order moment is obtained for some specific types of expansion. In the case of an exponential expansion, exact recurrence relations for the Laplace-transformed moments are obtained, whence the long-time behavior of moments of arbitrary order is subsequently inferred. Our analytical and numerical results for both Lévy flights and subdiffusive CTRWs confirm the intuitive expectation that the medium expansion hinders the mixing of diffusive particles occupying separate regions. In the case of Lévy flights, we quantify this effect by means of the so-called "Lévy horizon."

Collective Space-Sensing Coordinates Pattern Scaling in Engineered Bacteria

Scale invariance refers to the maintenance of a constant ratio of developing organ size to body size. Although common, its underlying mechanisms remain poorly understood. Here, we examined scaling in engineered Escherichia coli that can form self-organized core-ring patterns in colonies. We found that the ring width exhibits perfect scale invariance to the colony size. Our analysis revealed a collective space-sensing mechanism, which entails sequential actions of an integral feedback loop and an incoherent feedforward loop. The integral feedback is implemented by the accumulation of a diffusive chemical produced by a colony. This accumulation, combined with nutrient consumption, sets the timing for ring initiation. The incoherent feedforward is implemented by the opposing effects of the domain size on the rate and duration of ring maturation. This mechanism emphasizes a role of timing control in achieving robust pattern scaling and provides a new perspective in examining the phenomenon in natural systems.

A Model of the Spatio-temporal Dynamics of Drosophila Eye Disc Development

Patterning and growth are linked during early development and have to be tightly controlled to result in a functional tissue or organ. During the development of the Drosophila eye, this linkage is particularly clear: the growth of the eye primordium mainly results from proliferating cells ahead of the morphogenetic furrow (MF), a moving signaling wave that sweeps across the tissue from the posterior to the anterior side, that induces proliferating cells anterior to it to differentiate and become cell cycle quiescent in its wake. Therefore, final eye disc size depends on the proliferation rate of undifferentiated cells and on the speed with which the MF sweeps across the eye disc. We developed a spatio-temporal model of the growing eye disc based on the regulatory interactions controlled by the signals Decapentaplegic (Dpp), Hedgehog (Hh) and the transcription factor Homothorax (Hth) and explored how the signaling patterns affect the movement of the MF and impact on eye disc growth. We used published and new quantitative data to parameterize the model. In particular, two crucial parameter values, the degradation rate of Hth and the diffusion coefficient of Hh, were measured. The model is able to reproduce the linear movement of the MF and the termination of growth of the primordium. We further show that the model can explain several mutant phenotypes, but fails to reproduce the previously observed scaling of the Dpp gradient in the anterior compartment.

How domain growth is implemented determines the long-term behavior of a cell population through its effect on spatial correlations

Domain growth plays an important role in many biological systems, and so the inclusion of domain growth in models of these biological systems is important to understanding how these systems function. In this work we present methods to include the effects of domain growth on the evolution of spatial correlations in a continuum approximation of a lattice-based model of cell motility and proliferation. We show that, depending on the way in which domain growth is implemented, different steady-state densities are predicted for an agent population. Furthermore, we demonstrate that the way in which domain growth is implemented can result in the evolution of the agent density depending on the size of the domain. Continuum approximations that ignore spatial correlations cannot capture these behaviors, while those that account for spatial correlations do. These results will be of interest to researchers in developmental biology, as they suggest that the nature of domain growth can determine the characteristics of cell populations.

Developmental Pattern Formation in Phases

Cells in developing organs undergo a series of changes in their transcriptional state until a complete repertoire of cell types is specified. These changes in cell identity, together with the control of tissue growth, determine the pattern of gene expression in the tissue. Recent studies explore the dynamics of pattern formation during development and provide new insights into the control mechanisms. Changes in morphogen signalling and transcriptional networks control the specification of cell types. This is often followed by a distinct second phase, where pattern is elaborated by tissue growth. Here, we discuss the transitions between distinct phases in pattern formation. We consider the implications of the underlying mechanisms for understanding how reproducible patterns form during development.

Temporal and spatial dynamics of scaling-specific features of a gene regulatory network in Drosophila

A widely appreciated aspect of developmental robustness is pattern formation in proportion to size. But how such scaling features emerge dynamically remains poorly understood. Here we generate a data set of the expression profiles of six gap genes in Drosophila melanogaster embryos that differ significantly in size. Expression patterns exhibit size-dependent dynamics both spatially and temporally. We uncover a dynamic emergence of under-scaling in the posterior, accompanied by reduced expression levels of gap genes near the middle of large embryos. Simulation results show that a size-dependent Bicoid gradient input can lead to reduced Krüppel expression that can have long-range and dynamic effects on gap gene expression in the posterior. Thus, for emergence of scaled patterns, the entire embryo may be viewed as a single unified dynamic system where maternally derived size-dependent information interpreted locally can be propagated in space and time as governed by the dynamics of a gene regulatory network.

How pattern formation is regulated relative to the size of an organism is unclear. Here, Wu et al. take data from gap gene expression in flies of different sizes together with simulations, identifying how scaling emerges dynamically and that local patterning influences global gene regulatory networks.

Read-Out of Dynamic Morphogen Gradients on Growing Domains

Quantitative data from the Drosophila wing imaginal disc reveals that the amplitude of the Decapentaplegic (Dpp) morphogen gradient increases continuously. It is an open question how cells can determine their relative position within a domain based on a continuously increasing gradient. Here we show that pre-steady state diffusion-based dispersal of morphogens results in a zone within the growing domain where the concentration remains constant over the patterning period. The position of the zone that is predicted based on quantitative data for the Dpp morphogen corresponds to where the Dpp-dependent gene expression boundaries of spalt (sal) and daughters against dpp (dad) emerge. The model also suggests that genes that are scaling and are expressed at lateral positions are either under the control of a different read-out mechanism or under the control of a different morphogen. The patterning mechanism explains the extraordinary robustness that is observed for variations in Dpp production, and offers an explanation for the dual role of Dpp in controlling patterning and growth. Pre-steady-state dynamics are pervasive in morphogen-controlled systems, thus making this a probable general mechanism for the scaled read-out of morphogen gradients in growing developmental systems.

An interplay of geometry and signaling enables robust lung branching morphogenesis

Early branching events during lung development are stereotyped. Although key regulatory components have been defined, the branching mechanism remains elusive. We have now used a developmental series of 3D geometric datasets of mouse embryonic lungs as well as time-lapse movies of cultured lungs to obtain physiological geometries and displacement fields. We find that only a ligand-receptor-based Turing model in combination with a particular geometry effect that arises from the distinct expression domains of ligands and receptors successfully predicts the embryonic areas of outgrowth and supports robust branch outgrowth. The geometry effect alone does not support bifurcating outgrowth, while the Turing mechanism alone is not robust to noisy initial conditions. The negative feedback between the individual Turing modules formed by fibroblast growth factor 10 (FGF10) and sonic hedgehog (SHH) enlarges the parameter space for which the embryonic growth field is reproduced. We therefore propose that a signaling mechanism based on FGF10 and SHH directs outgrowth of the lung bud via a ligand-receptor-based Turing mechanism and a geometry effect.