Positional information and pattern formation

Mechanical positional information guides the self-organized development of a polygonal network of creases in the skin of mammalian noses

The glabrous skin of the rhinarium (naked nose) of many mammalian species exhibits a polygonal pattern of grooves that retain physiological fluid, thereby keeping their nose wet and, among other effects, facilitating the collection of chemosensory molecules. Here, we perform volumetric imaging of whole-mount rhinaria from sequences of embryonic and juvenile cows, dogs, and ferrets. We demonstrate that rhinarial polygonal domains are not placode-derived skin appendages but arise through a self-organized mechanical process consisting of the constrained growth and buckling of the epidermal basal layer, followed by the formation of sharp epidermal creases exactly facing an underlying network of stiff blood vessels. Our numerical simulations show that the mechanical stress generated by excessive epidermal growth concentrates at the positions of vessels that form rigid base points, causing the epidermal layers to move outward and shape domes-akin to arches rising against stiff pillars. Remarkably, this gives rise to a larger length scale (the distance between the vessels) in the surface folding pattern than would otherwise occur in the absence of vessels. These results hint at a concept of "mechanical positional information" by which material properties of anatomical elements can impose local constraints on an otherwise globally self-organized mechanical pattern. In addition, our analyses of the rhinarial patterns in cow clones highlight a substantial level of stochasticity in the pre-pattern of vessels, while our numerical simulations also recapitulate the disruption of the folding pattern in cows affected by a hereditary disorder that causes hyperextensibility of the skin.

Mechanics of reproductive differentiation in the land plants: a paradigm shift?

This article addresses the physical mechanics of gametogenesis in vascular plants. The earliest events that initiate reproductive differentiation in the land plants are not well understood. How are the few cells that initiate reproductive differentiation specified and how is that information translated into action at the cellular level? In this article I propose a physical mechanism that resolves the problem of spatial targeting without invoking dependence on diffusible morphogens or other external signals. I suggest that the initiation of archesporial differentiation can instead be attributed to the confluence of organ geometry, surficial topography, and the physical mechanics of sporangial growth, resulting in the spontaneous emergence of an isotropic singularity that locates and precipitates archesporial differentiation. In discussing the logic of single-cell target selection and the limits of stochastic molecular signaling I propose that the sporangium would be better understood as a pressurized stress-mechanical lens that focuses turgor-generated growth forces on a central location, generating a physical singularity that locates and specifies the cell or cells that become the archesporium and initiates their transition from somatic proliferation to reproductive differentiation.

Evolutionary Implications of Self-Assembling Cybernetic Materials with Collective Problem-Solving Intelligence at Multiple Scales

In recent years, the scientific community has increasingly recognized the complex multi-scale competency architecture (MCA) of biology, comprising nested layers of active homeostatic agents, each forming the self-orchestrated substrate for the layer above, and, in turn, relying on the structural and functional plasticity of the layer(s) below. The question of how natural selection could give rise to this MCA has been the focus of intense research. Here, we instead investigate the effects of such decision-making competencies of MCA agential components on the process of evolution itself, using in silico neuroevolution experiments of simulated, minimal developmental biology. We specifically model the process of morphogenesis with neural cellular automata (NCAs) and utilize an evolutionary algorithm to optimize the corresponding model parameters with the objective of collectively self-assembling a two-dimensional spatial target pattern (reliable morphogenesis). Furthermore, we systematically vary the accuracy with which the uni-cellular agents of an NCA can regulate their cell states (simulating stochastic processes and noise during development). This allows us to continuously scale the agents' competency levels from a direct encoding scheme (no competency) to an MCA (with perfect reliability in cell decision executions). We demonstrate that an evolutionary process proceeds much more rapidly when evolving the functional parameters of an MCA compared to evolving the target pattern directly. Moreover, the evolved MCAs generalize well toward system parameter changes and even modified objective functions of the evolutionary process. Thus, the adaptive problem-solving competencies of the agential parts in our NCA-based in silico morphogenesis model strongly affect the evolutionary process, suggesting significant functional implications of the near-ubiquitous competency seen in living matter.

Periodic pattern formation during embryonic development

During embryonic development many organs and structures require the formation of series of repeating elements known as periodic patterns. Ranging from the digits of the limb to the feathers of the avian skin, the correct formation of these embryonic patterns is essential for the future form and function of these tissues. However, the mechanisms that produce these patterns are not fully understood due to the existence of several modes of pattern generation which often differ between organs and species. Here, we review the current state of the field and provide a perspective on future approaches to studying this fundamental process of embryonic development.

Morphogens enable interacting supracellular phases that generate organ architecture

During vertebrate organogenesis, increases in morphological complexity are tightly coupled to morphogen expression. In this work, we studied how morphogens influence self-organizing processes at the collective or "supra"-cellular scale in avian skin. We made physical measurements across length scales, which revealed morphogen-enabled material property differences that were amplified at supracellular scales in comparison to cellular scales. At the supracellular scale, we found that fibroblast growth factor (FGF) promoted "solidification" of tissues, whereas bone morphogenetic protein (BMP) promoted fluidity and enhanced mechanical activity. Together, these effects created basement membrane-less compartments within mesenchymal tissue that were mechanically primed to drive avian skin tissue budding. Understanding this multiscale process requires the ability to distinguish between proximal effects of morphogens that occur at the cellular scale and their functional effects, which emerge at the supracellular scale.

The Unreasonable Effectiveness of Reaction Diffusion in Vertebrate Skin Color Patterning

In 1952, Alan Turing published the reaction-diffusion (RD) mathematical framework, laying the foundations of morphogenesis as a self-organized process emerging from physicochemical first principles. Regrettably, this approach has been widely doubted in the field of developmental biology. First, we summarize Turing's line of thoughts to alleviate the misconception that RD is an artificial mathematical construct. Second, we discuss why phenomenological RD models are particularly effective for understanding skin color patterning at the meso/macroscopic scales, without the need to parameterize the profusion of variables at lower scales. More specifically, we discuss how RD models (a) recapitulate the diversity of actual skin patterns, (b) capture the underlying dynamics of cellular interactions, (c) interact with tissue size and shape, (d) can lead to ordered sequential patterning, (e) generate cellular automaton dynamics in lizards and snakes, (f) predict actual patterns beyond their statistical features, and (g) are robust to model variations. Third, we discuss the utility of linear stability analysis and perform numerical simulations to demonstrate how deterministic RD emerges from the underlying chaotic microscopic agents.

The Hydra stem cell system - Revisited

Cnidarians are >600 million years old and are considered the sister group of Bilateria based on numerous molecular phylogenetic studies. Apart from Hydra, the genomes of all major clades of Cnidaria have been uncovered (e.g. Aurelia, Clytia, Nematostella and Acropora) and they reveal a remarkable completeness of the metazoan genomic toolbox. Of particular interest is Hydra, a model system of aging research, regenerative biology, and stem cell biology. With the knowledge gained from scRNA research, it is now possible to characterize the expression profiles of all cell types with great precision. In functional studies, our picture of the Hydra stem cell biology has changed, and we are in the process of obtaining a clear picture of the homeostasis and properties of the different stem cell populations. Even though Hydra is often compared to plant systems, the new data on germline and regeneration, but also on the dynamics and plasticity of the nervous system, show that Hydra with its simple body plan represents in a nutshell the prototype of an animal with stem cell lineages, whose properties correspond in many ways to Bilateria. This review provides an overview of the four stem cell lineages, the two epithelial lineages that constitute the ectoderm and the endoderm, as well as the multipotent somatic interstitial lineage (MPSC) and the germline stem cell lineage (GSC), also known as the interstitial cells of Hydra.

Embryonic Tissue and Blastema Transplantations

Embryo grafts have been an experimental pillar in developmental biology, and particularly, in amphibian biology. Grafts have been essential in constructing fate maps of different cell populations and migratory patterns. Likewise, autografts and allografts in older larvae or adult salamanders have been widely used to disentangle mechanisms of regeneration. The combination of transgenesis and grafting has widened even more the application of this technique.In this chapter, we provide a detailed protocol for embryo transplants in the axolotl (Ambystoma mexicanum ). The location and stages to label connective tissue, muscle, or blood vessels in the limb and blood cells in the whole animal. However, the potential of embryo transplants is enormous and impossible to cover in one chapter. Furthermore, we provide a protocol for blastema transplantation as an example of allograft in older larvae.

Reaction-diffusion phenomena in antagonistic bipolar diffusion fields

Operating natural or artificial chemical systems requires nonequilibrium conditions under which temporal and spatial control of the process is realizable. Open reaction-diffusion systems provide a general way to create such conditions. A key issue is the proper design of reactors in which the nonequilibrium conditions can be maintained. A hydrogel with flow-through channels is a simple, flexible, and easy-to-make device in which chemical reactions are performed in the diffusion field of localized separated sources of reactants. Two reactants separated in two channels create a bipolar antagonistic diffusion field, where the reaction intermediates firmly separate in space. Numerical simulations and corresponding experiments are performed to present this inhomogeneous diffusion field-induced chemical separation in sequential reactions. A remarkable result of this bipolar spatial control is localized wave phenomena in a nonlinear activatory-inhibitory reaction. These findings may help design functioning artificial nonequilibrium systems with the desired spatial separation of chemicals.

Modern perspectives on near-equilibrium analysis of Turing systems

In the nearly seven decades since the publication of Alan Turing's work on morphogenesis, enormous progress has been made in understanding both the mathematical and biological aspects of his proposed reaction-diffusion theory. Some of these developments were nascent in Turing's paper, and others have been due to new insights from modern mathematical techniques, advances in numerical simulations and extensive biological experiments. Despite such progress, there are still important gaps between theory and experiment, with many examples of biological patterning where the underlying mechanisms are still unclear. Here, we review modern developments in the mathematical theory pioneered by Turing, showing how his approach has been generalized to a range of settings beyond the classical two-species reaction-diffusion framework, including evolving and complex manifolds, systems heterogeneous in space and time, and more general reaction-transport equations. While substantial progress has been made in understanding these more complicated models, there are many remaining challenges that we highlight throughout. We focus on the mathematical theory, and in particular linear stability analysis of 'trivial' base states. We emphasize important open questions in developing this theory further, and discuss obstacles in using these techniques to understand biological reality. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Mathematical modeling of the eyespots in butterfly wings

Butterfly wing color patterns are a representative model system for studying biological pattern formation, due to their two-dimensional simple structural and high inter- and intra-specific variabilities. Moreover, butterfly color patterns have demonstrated roles in mate choice, thermoregulation, and predator avoidance via disruptive coloration, attack deflection, aposematism, mimicry, and masquerade. Because of the importance of color patterns to many aspects of butterfly biology and their apparent tractability for study, color patterns have been the subjects of many attempts to model their development. Early attempts focused on generalized mechanisms of pattern formation such as reaction-diffusion, diffusion gradient, lateral inhibition, and threshold responses, without reference to any specific gene products. As candidate genes with expression patterns that resembled incipient color patterns were identified, genetic regulatory networks were proposed for color pattern formation based on gene functions inferred from other insects with wings, such as Drosophila. Particularly detailed networks incorporating the gene products, Distal-less (Dll), Engrailed (En), Hedgehog (Hh), Cubitus interruptus (Ci), Transforming growth factor-β (TGF-β), and Wingless (Wg), have been proposed for butterfly border ocelli (eyespots) which helps the investigation of the formation of these patterns. Thus, in this work, we develop a mathematical model including the gene products En, Hh, Ci, TGF-β, and Wg to mimic and investigate the eyespot formation in butterflies. Our simulations show that the level of En has peaks in the inner and outer rings and the level of Ci has peaks in the inner and middle rings. The interactions among these peaks activate cells to produce white, black, and yellow pigments in the inner, middle, and outer rings, respectively, which captures the eyespot pattern of wild type Bicyclus anynana butterflies. Additionally, our simulations suggest that lack of En generates a single black spot and lack of Hh or Ci generates a single white spot, and a deficiency of TGF-β or Wg will cause the loss of the outer yellow ring. These deficient patterns are similar to those observed in the eyespots of Vanessa atalanta, Vanessa altissima, and Chlosyne nycteis. Thus, our model also provides a hypothesis to explain the mechanism of generating the deficient patterns in these species.

Placozoan fiber cells: mediators of innate immunity and participants in wound healing

Placozoa is a phylum of non-bilaterian marine animals. These small, flat organisms adhere to the substrate via their densely ciliated ventral epithelium, which mediates mucociliary locomotion and nutrient uptake. They have only six morphological cell types, including one, fiber cells, for which functional data is lacking. Fiber cells are non-epithelial cells with multiple processes. We used electron and light microscopic approaches to unravel the roles of fiber cells in Trichoplax adhaerens, a representative member of the phylum. Three-dimensional reconstructions of serial sections of Trichoplax showed that each fiber cell is in contact with several other cells. Examination of fiber cells in thin sections and observations of live dissociated fiber cells demonstrated that they phagocytose cell debris and bacteria. In situ hybridization confirmed that fiber cells express genes involved in phagocytic activity. Fiber cells also are involved in wound healing as evidenced from microsurgery experiments. Based on these observations we conclude that fiber cells are multi-purpose macrophage-like cells. Macrophage-like cells have been described in Porifera, Ctenophora, and Cnidaria and are widespread among Bilateria, but our study is the first to show that Placozoa possesses this cell type. The phylogenetic distribution of macrophage-like cells suggests that they appeared early in metazoan evolution.

The Emerging Mechanisms of Wnt Secretion and Signaling in Development

Wnts are highly-conserved lipid-modified secreted proteins that activate multiple signaling pathways. These pathways regulate crucial processes during various stages of development and maintain tissue homeostasis in adults. One of the most fascinating aspects of Wnt protein is that despite being hydrophobic, they are known to travel several cell distances in the extracellular space. Research on Wnts in the past four decades has identified several factors and uncovered mechanisms regulating their expression, secretion, and mode of extracellular travel. More recently, analyses on the importance of Wnt protein gradients in the growth and patterning of developing tissues have recognized the complex interplay of signaling mechanisms that help in maintaining tissue homeostasis. This review aims to present an overview of the evidence for the various modes of Wnt protein secretion and signaling and discuss mechanisms providing precision and robustness to the developing tissues.

COUP-TFI specifies the medial entorhinal cortex identity and induces differential cell adhesion to determine the integrity of its boundary with neocortex

Development of cortical regions with precise, sharp, and regular boundaries is essential for physiological function. However, little is known of the mechanisms ensuring these features. Here, we show that determination of the boundary between neocortex and medial entorhinal cortex (MEC), two abutting cortical regions generated from the same progenitor lineage, relies on COUP-TFI (chicken ovalbumin upstream promoter-transcription factor I), a patterning transcription factor with graded expression in cortical progenitors. In contrast with the classical paradigm, we found that increased COUP-TFI expression expands MEC, creating protrusions and disconnected ectopic tissue. We further developed a mathematical model that predicts that neuronal specification and differential cell affinity contribute to the emergence of an instability region and boundary sharpness. Correspondingly, we demonstrated that high expression of COUP-TFI induces MEC cell fate and protocadherin 19 expression. Thus, we conclude that a sharp boundary requires a subtle interplay between patterning transcription factors and differential cell affinity.

Lewis Wolpert (1929-2021)

Lewis Wolpert, who died on 28 January 2021, was an inspirational figure to generations of developmental biologists, and a man whose influence extended far beyond his subject. By seeing significance in the commonplace, and by identifying and re-framing research questions, he inspired new ways of thinking about embryonic development. As a teacher and as a lecturer he was nonpareil, and his ideas were brilliantly distilled in his textbook Principles of Development, which will continue to imbue students with his clarity of thought. In writing about his depression, Wolpert reduced the stigma associated with the disease and provided scientific understanding, solace and hope. An open, kind and generous man, he went out of his way to support young scientists.

Emergence of self-organizational patterning at the mesoscopic scale

In this issue of Developmental Cell, Chuyen et al. (2021) suggest that the Scf/Kit pathway controls mutual repulsion of multiciliated cells and their affinity for epidermal cell junctions through soluble and membrane-associated Scf ligands, respectively. Effective self-organizational patterning emerges at the mesoscopic scale as a small set of effective behaviors.

On the evolution and development of morphological complexity: A view from gene regulatory networks

How does morphological complexity evolve? This study suggests that the likelihood of mutations increasing phenotypic complexity becomes smaller when the phenotype itself is complex. In addition, the complexity of the genotype-phenotype map (GPM) also increases with the phenotypic complexity. We show that complex GPMs and the above mutational asymmetry are inevitable consequences of how genes need to be wired in order to build complex and robust phenotypes during development. We randomly wired genes and cell behaviors into networks in EmbryoMaker. EmbryoMaker is a mathematical model of development that can simulate any gene network, all animal cell behaviors (division, adhesion, apoptosis, etc.), cell signaling, cell and tissues biophysics, and the regulation of those behaviors by gene products. Through EmbryoMaker we simulated how each random network regulates development and the resulting morphology (i.e. a specific distribution of cells and gene expression in 3D). This way we obtained a zoo of possible 3D morphologies. Real gene networks are not random, but a random search allows a relatively unbiased exploration of what is needed to develop complex robust morphologies. Compared to the networks leading to simple morphologies, the networks leading to complex morphologies have the following in common: 1) They are rarer; 2) They need to be finely tuned; 3) Mutations in them tend to decrease morphological complexity; 4) They are less robust to noise; and 5) They have more complex GPMs. These results imply that, when complexity evolves, it does so at a progressively decreasing rate over generations. This is because as morphological complexity increases, the likelihood of mutations increasing complexity decreases, morphologies become less robust to noise, and the GPM becomes more complex. We find some properties in common, but also some important differences, with non-developmental GPM models (e.g. RNA, protein and gene networks in single cells).

Appendage Regeneration in Vertebrates: What Makes This Possible?

The ability to regenerate amputated or injured tissues and organs is a fascinating property shared by several invertebrates and, interestingly, some vertebrates. The mechanism of evolutionary loss of regeneration in mammals is not understood, yet from the biomedical and clinical point of view, it would be very beneficial to be able, at least partially, to restore that capability. The current availability of new experimental tools, facilitating the comparative study of models with high regenerative ability, provides a powerful instrument to unveil what is needed for a successful regeneration. The present review provides an updated overview of multiple aspects of appendage regeneration in three vertebrates: lizard, salamander, and zebrafish. The deep investigation of this process points to common mechanisms, including the relevance of Wnt/β-catenin and FGF signaling for the restoration of a functional appendage. We discuss the formation and cellular origin of the blastema and the identification of epigenetic and cellular changes and molecular pathways shared by vertebrates capable of regeneration. Understanding the similarities, being aware of the differences of the processes, during lizard, salamander, and zebrafish regeneration can provide a useful guide for supporting effective regenerative strategies in mammals.

Salient brain entities labelled in P2rx7-EGFP reporter mouse embryos include the septum, roof plate glial specializations and circumventricular ependymal organs

The purinergic system is one of the oldest cell-to-cell communication mechanisms and exhibits relevant functions in the regulation of the central nervous system (CNS) development. Amongst the components of the purinergic system, the ionotropic P2X7 receptor (P2X7R) stands out as a potential regulator of brain pathology and physiology. Thus, P2X7R is known to regulate crucial aspects of neuronal cell biology, including axonal elongation, path-finding, synapse formation and neuroprotection. Moreover, P2X7R modulates neuroinflammation and is posed as a therapeutic target in inflammatory, oncogenic and degenerative disorders. However, the lack of reliable technical and pharmacological approaches to detect this receptor represents a major hurdle in its study. Here, we took advantage of the P2rx7-EGFP reporter mouse, which expresses enhanced green fluorescence protein (EGFP) immediately downstream of the P2rx7 proximal promoter, to conduct a detailed study of its distribution. We performed a comprehensive analysis of the pattern of P2X7R expression in the brain of E18.5 mouse embryos revealing interesting areas within the CNS. Particularly, strong labelling was found in the septum, as well as along the entire neural roof plate zone of the brain, except chorioidal roof areas, but including specialized circumventricular roof formations, such as the subfornical and subcommissural organs (SFO; SCO). Moreover, our results reveal what seems a novel circumventricular organ, named by us postarcuate organ (PArcO). Furthermore, this study sheds light on the ongoing debate regarding the specific presence of P2X7R in neurons and may be of interest for the elucidation of additional roles of P2X7R in the idiosyncratic histologic development of the CNS and related systemic functions.

Mechanics of Development

Mechanical forces are integral to development-from the earliest stages of embryogenesis to the construction and differentiation of complex organs. Advances in imaging and biophysical tools have allowed us to delve into the developmental mechanobiology of increasingly complex organs and organisms. Here, we focus on recent work that highlights the diversity and importance of mechanical influences during morphogenesis. Developing tissues experience intrinsic mechanical signals from active forces and changes to tissue mechanical properties as well as extrinsic mechanical signals, including constraint and compression, pressure, and shear forces. Finally, we suggest promising avenues for future work in this rapidly expanding field.

Transcriptional Enhancers in Drosophila

Key discoveries in Drosophila have shaped our understanding of cellular "enhancers." With a special focus on the fly, this chapter surveys properties of these adaptable cis-regulatory elements, whose actions are critical for the complex spatial/temporal transcriptional regulation of gene expression in metazoa. The powerful combination of genetics, molecular biology, and genomics available in Drosophila has provided an arena in which the developmental role of enhancers can be explored. Enhancers are characterized by diverse low- or high-throughput assays, which are challenging to interpret, as not all of these methods of identifying enhancers produce concordant results. As a model metazoan, the fly offers important advantages to comprehensive analysis of the central functions that enhancers play in gene expression, and their critical role in mediating the production of phenotypes from genotype and environmental inputs. A major challenge moving forward will be obtaining a quantitative understanding of how these cis-regulatory elements operate in development and disease.

Multistable and dynamic CRISPRi-based synthetic circuits

Gene expression control based on CRISPRi (clustered regularly interspaced short palindromic repeats interference) has emerged as a powerful tool for creating synthetic gene circuits, both in prokaryotes and in eukaryotes; yet, its lack of cooperativity has been pointed out as a potential obstacle for dynamic or multistable synthetic circuit construction. Here we use CRISPRi to build a synthetic oscillator (“CRISPRlator”), bistable network (toggle switch) and stripe pattern-forming incoherent feed-forward loop (IFFL). Our circuit designs, conceived to feature high predictability and orthogonality, as well as low metabolic burden and context-dependency, allow us to achieve robust circuit behaviors in Escherichia coli populations. Mathematical modeling suggests that unspecific binding in CRISPRi is essential to establish multistability. Our work demonstrates the wide applicability of CRISPRi in synthetic circuits and paves the way for future efforts towards engineering more complex synthetic networks, boosted by the advantages of CRISPR technology.

Synthetic circuits based on CRISPRi have not achieved multistable and dynamic behaviors. Here the authors build an oscillator, a toggle switch and an incoherent feed-forward loop using CRISPRi, and provide a mathematical model suggesting that unspecific binding in CRISPRi enables multistability.

Can the French flag and reaction-diffusion models explain flower patterning? Celebrating the 50th anniversary of the French flag model

It has been 50 years since Lewis Wolpert introduced the French flag model proposing the patterning of different cell types based on threshold concentrations of a morphogen diffusing in the tissue. Sixty-seven years ago, Alan Turing introduced the idea of patterns initiating de novo from a reaction-diffusion network. Together these models have been used to explain many patterning events in animal development, so here we take a look at their applicability to flower development. First, although many plant transcription factors move through plasmodesmata from cell to cell, in the flower there is little evidence that they specify fate in a concentration-dependent manner, so they cannot yet be described as morphogens. Secondly, the reaction-diffusion model appears to be a reasonably good description of the formation of spots of pigment on petals, although additional nuances are present. Thirdly, aspects of both of these combine in a new fluctuation-based patterning system creating the scattered pattern of giant cells in Arabidopsis sepals. In the future, more precise imaging and manipulations of the dynamics of patterning networks combined with mathematical modeling will allow us to better understand how the multilayered complex and beautiful patterns of flowers emerge de novo.

Two MYB Proteins in a Self-Organizing Activator-Inhibitor System Produce Spotted Pigmentation Patterns

Many organisms exhibit visually striking spotted or striped pigmentation patterns. Developmental models predict that such spatial patterns can form when a local autocatalytic feedback loop and a long-range inhibitory feedback loop interact. At its simplest, this self-organizing network only requires one self-activating activator that also activates a repressor, which inhibits the activator and diffuses to neighboring cells. However, the molecular activators and inhibitors fully fitting this versatile model remain elusive in pigmentation systems. Here, we characterize an R2R3-MYB activator and an R3-MYB repressor in monkeyflowers (Mimulus). Through experimental perturbation and mathematical modeling, we demonstrate that the properties of these two proteins correspond to an activator-inhibitor pair in a two-component, reaction-diffusion system, explaining the formation of dispersed anthocyanin spots in monkeyflower petals. Notably, disrupting this pattern impacts pollinator visitation. Thus, subtle changes in simple activator-inhibitor systems are likely essential contributors to the evolution of the remarkable diversity of pigmentation patterns in flowers.

MicroRNAs organize intrinsic variation into stem cell states

Understanding how mammalian organisms achieve the full diversity of cell types in the adult organism is a central goal of developmental cell biology. Recent work has shown that some embryonic precursor cells can self-organize into developmental structures but the mechanisms of gene regulation that contribute to this process remain unknown. Here we show embryonic stem cells self-organize into distinct gene expression states that resemble developmental gene programs. We find that microRNAs, small noncoding regulators of gene expression, play a critical role in organizing fluctuations across gene networks to help achieve this organization into distinct expression states.

Pluripotent embryonic stem cells (ESCs) contain the potential to form a diverse array of cells with distinct gene expression states, namely the cells of the adult vertebrate. Classically, diversity has been attributed to cells sensing their position with respect to external morphogen gradients. However, an alternative is that diversity arises in part from cooption of fluctuations in the gene regulatory network. Here we find ESCs exhibit intrinsic heterogeneity in the absence of external gradients by forming interconverting cell states. States vary in developmental gene expression programs and display distinct activity of microRNAs (miRNAs). Notably, miRNAs act on neighborhoods of pluripotency genes to increase variation of target genes and cell states. Loss of miRNAs that vary across states reduces target variation and delays state transitions, suggesting variable miRNAs organize and propagate variation to promote state transitions. Together these findings provide insight into how a gene regulatory network can coopt variation intrinsic to cell systems to form robust gene expression states. Interactions between intrinsic heterogeneity and environmental signals may help achieve developmental outcomes.

Remembering where we are: Positional information in salamander limb regeneration

Fifty years ago, Lewis Wolpert defined an important question in developmental biology: how are cell fates determined by the positions of cells within a system? He proposed that cells retain positional values as if they lie within a coordinate system and that the interpretation of these values produces patterns in development. He referred to this concept as positional information. Though initially controversial, this concept of positional information has proven to be profoundly influential in developmental biology. One area in which the influence of Wolpert's theoretical work can be clearly demonstrated is the study of limb regeneration in salamanders. Here, we review the work in limb regeneration leading up to Wolpert defining the concept of positional information and how his theory has guided regeneration research over the subsequent 50 years.

Distinct Mesenchymal Cell Populations Generate the Essential Intestinal BMP Signaling Gradient

Intestinal stem cells (ISCs) are confined to crypt bottoms and their progeny differentiate near crypt-villus junctions. Wnt and bone morphogenic protein (BMP) gradients drive this polarity, and colorectal cancer fundamentally reflects disruption of this homeostatic signaling. However, sub-epithelial sources of crucial agonists and antagonists that organize this BMP gradient remain obscure. Here, we couple whole-mount high-resolution microscopy with ensemble and single-cell RNA sequencing (RNA-seq) to identify three distinct PDGFRA⁺ mesenchymal cell types. PDGFRA(hi) telocytes are especially abundant at the villus base and provide a BMP reservoir, and we identified a CD81⁺ PDGFRA(lo) population present just below crypts that secretes the BMP antagonist Gremlin1. These cells, referred to as trophocytes, are sufficient to expand ISCs in vitro without additional trophic support and contribute to ISC maintenance in vivo. This study reveals intestinal mesenchymal structure at fine anatomic, molecular, and functional detail and the cellular basis for a signaling gradient necessary for tissue self-renewal.

Mathematical modeling of chondrogenic pattern formation during limb development: Recent advances in continuous models

The phenomenon of chondrogenic pattern formation in the vertebrate limb is one of the best studied examples of organogenesis. Many different models, mathematical as well as conceptual, have been proposed for it in the last fifty years or so. In this review, we give a brief overview of the fundamental biological background, then describe in detail several models which aim to describe qualitatively and quantitatively the corresponding biological phenomena. We concentrate on several new models that have been proposed in recent years, taking into account recent experimental progress. The major mathematical tools in these approaches are ordinary and partial differential equations. Moreover, we discuss models with non-local flux terms used to account for cell-cell adhesion forces and a structured population model with diffusion. We also include a detailed list of gene products and potential morphogens which have been identified to play a role in the process of limb formation and its growth.

Dynamics of Spaetzle morphogen shuttling in the Drosophila embryo shapes gastrulation patterning

Establishment of morphogen gradients in the early Drosophila embryo is challenged by a diffusible sextracellular milieu, and by rapid nuclear divisions that occur at the same time. To understand how a sharp gradient is formed within this dynamic environment, we followed the generation of graded nuclear Dorsal protein, the hallmark of pattern formation along the dorso-ventral axis, in live embryos. The dynamics indicate that a sharp extracellular gradient is formed through diffusion-based shuttling of the Spaetzle (Spz) morphogen that progresses through several nuclear divisions. Perturbed shuttling in wntD mutant embryos results in a flat activation peak and aberrant gastrulation. Re-entry of Dorsal into the nuclei at the final division cycle plays an instructive role, as the residence time of Dorsal in each nucleus is translated to the amount of zygotic transcript that will be produced, thereby guiding graded accumulation of specific zygotic transcripts that drive patterned gastrulation. We conclude that diffusion-based ligand shuttling, coupled with dynamic readout, establishes a refined pattern within the diffusible environment of early embryos.

Robust Wnt signaling is maintained by a Wg protein gradient and Fz2 receptor activity in the developing Drosophila wing

Wnts are secreted proteins that regulate cell fate during development of all metazoans. Wnt proteins were proposed to spread over several cells to activate signaling directly at a distance. In the Drosophila wing epithelium, an extracellular gradient of the Wnt1 homolog Wingless (Wg) was observed extending over several cells away from producing cells. Surprisingly, however, it was also shown that a membrane-tethered Neurotactin-Wg fusion protein (NRT-Wg) can largely replace endogenous Wg, leading to proper patterning of the wing. Therefore, the functional range of Wg and whether Wg spreading is required for correct tissue patterning remains controversial. Here, by capturing secreted Wg on cells away from the source, we show that Wg acts over a distance of up to 11 cell diameters to induce signaling. Furthermore, cells located outside the reach of extracellular Wg depend on the Frizzled2 receptor to maintain signaling. Frizzled2 expression is increased in the absence of Wg secretion and is required to maintain signaling and cell survival in NRT-wg wing discs. Together, these results provide insight into the mechanisms by which robust Wnt signaling is achieved in proliferating tissues.

Agent-based modeling of morphogenetic systems: Advantages and challenges

The complexity of morphogenesis poses a fundamental challenge to understanding the mechanisms governing the formation of biological patterns and structures. Over the past century, numerous processes have been identified as critically contributing to morphogenetic events, but the interplay between the various components and aspects of pattern formation have been much harder to grasp. The combination of traditional biology with mathematical and computational methods has had a profound effect on our current understanding of morphogenesis and led to significant insights and advancements in the field. In particular, the theoretical concepts of reaction–diffusion systems and positional information, proposed by Alan Turing and Lewis Wolpert, respectively, dramatically influenced our general view of morphogenesis, although typically in isolation from one another. In recent years, agent-based modeling has been emerging as a consolidation and implementation of the two theories within a single framework. Agent-based models (ABMs) are unique in their ability to integrate combinations of heterogeneous processes and investigate their respective dynamics, especially in the context of spatial phenomena. In this review, we highlight the benefits and technical challenges associated with ABMs as tools for examining morphogenetic events. These models display unparalleled flexibility for studying various morphogenetic phenomena at multiple levels and have the important advantage of informing future experimental work, including the targeted engineering of tissues and organs.

Modeling Edar expression reveals the hidden dynamics of tooth signaling center patterning

When patterns are set during embryogenesis, it is expected that they are straightly established rather than subsequently modified. The patterning of the three mouse molars is, however, far from straight, likely as a result of mouse evolutionary history. The first-formed tooth signaling centers, called MS and R2, disappear before driving tooth formation and are thought to be vestiges of the premolars found in mouse ancestors. Moreover, the mature signaling center of the first molar (M1) is formed from the fusion of two signaling centers (R2 and early M1). Here, we report that broad activation of Edar expression precedes its spatial restriction to tooth signaling centers. This reveals a hidden two-step patterning process for tooth signaling centers, which was modeled with a single activator-inhibitor pair subject to reaction-diffusion (RD). The study of Edar expression also unveiled successive phases of signaling center formation, erasing, recovering, and fusion. Our model, in which R2 signaling center is not intrinsically defective but erased by the broad activation preceding M1 signaling center formation, predicted the surprising rescue of R2 in Edar mutant mice, where activation is reduced. The importance of this R2-M1 interaction was confirmed by ex vivo cultures showing that R2 is capable of forming a tooth. Finally, by introducing chemotaxis as a secondary process to RD, we recapitulated in silico different conditions in which R2 and M1 centers fuse or not. In conclusion, pattern formation in the mouse molar field relies on basic mechanisms whose dynamics produce embryonic patterns that are plastic objects rather than fixed end points.

Engineered Living Materials: Prospects and Challenges for Using Biological Systems to Direct the Assembly of Smart Materials

Vast potential exists for the development of novel, engineered platforms that manipulate biology for the production of programmed advanced materials. Such systems would possess the autonomous, adaptive, and self-healing characteristics of living organisms, but would be engineered with the goal of assembling bulk materials with designer physicochemical or mechanical properties, across multiple length scales. Early efforts toward such engineered living materials (ELMs) are reviewed here, with an emphasis on engineered bacterial systems, living composite materials which integrate inorganic components, successful examples of large-scale implementation, and production methods. In addition, a conceptual exploration of the fundamental criteria of ELM technology and its future challenges is presented. Cradled within the rich intersection of synthetic biology and self-assembling materials, the development of ELM technologies allows the power of biology to be leveraged to grow complex structures and objects using a palette of bio-nanomaterials.

CONSERVATION OF MOLECULAR PREPATTERNS DURING THE EVOLUTION OF CUTICLE MORPHOLOGY IN DROSOPHILA LARVAE

We are using patterns of cuticle specialization in Drosophila larvae as models to investigate the molecular, genetic, and developmental bases of morphological evolution. Members of the virilis species group differ markedly from one another in the distribution of hairs on the dorsal surface of first instar larvae. In particular, characteristic bands of hairs cover about 20% of each trunk segment in some species but about 70% in others. These major types do not correlate with recently proposed phylogenetic relationships, suggesting that similar phenotypes have arisen independently in different lineages. The patterns of expression of several genes that control or reflect intrasegmental patterning are indistinguishable in species with very different cuticle morphologies. We conclude that, in this case, morphology probably has evolved via altered response to a conserved molecular prepattern.

Vertebrate Axial Patterning: From Egg to Asymmetry

The emergence of the bilateral embryonic body axis from a symmetrical egg has been a long-standing question in developmental biology. Historical and modern experiments point to an initial symmetry-breaking event leading to localized Wnt and Nodal growth factor signaling and subsequent induction and formation of a self-regulating dorsal "organizer." This organizer forms at the site of notochord cell internalization and expresses primarily Bone Morphogenetic Protein (BMP) growth factor antagonists that establish a spatiotemporal gradient of BMP signaling across the embryo, directing initial cell differentiation and morphogenesis. Although the basics of this model have been known for some time, many of the molecular and cellular details have only recently been elucidated and the extent that these events remain conserved throughout vertebrate evolution remains unclear. This chapter summarizes historical perspectives as well as recent molecular and genetic advances regarding: (1) the mechanisms that regulate symmetry-breaking in the vertebrate egg and early embryo, (2) the pathways that are activated by these events, in particular the Wnt pathway, and the role of these pathways in the formation and function of the organizer, and (3) how these pathways also mediate anteroposterior patterning and axial morphogenesis. Emphasis is placed on comparative aspects of the egg-to-embryo transition across vertebrates and their evolution. The future prospects for work regarding self-organization and gene regulatory networks in the context of early axis formation are also discussed.

On the developmental self-regulatory dynamics and evolution of individuated multicellular organisms

Changes in gene expression are thought to regulate the cell differentiation process intrinsically through complex epigenetic mechanisms. In fundamental terms, however, this assumed regulation refers only to the intricate propagation of changes in gene expression or else leads to non-explanatory regresses. The developmental self-regulatory dynamics and evolution of individuated multicellular organisms also lack a unified and falsifiable description. To fill this gap, I computationally analyzed publicly available high-throughput data of histone H3 post-translational modifications and mRNA abundance for different Homo sapiens, Mus musculus, and Drosophila melanogaster cell-type/developmental-period samples. My analysis of genomic regions adjacent to transcription start sites generated a profile from pairwise partial correlations between histone modifications controlling for the respective mRNA levels for each cell-type/developmental-period dataset. I found that these profiles, while explicitly uncorrelated with the respective transcriptional "identities" by construction, associate strongly with cell differentiation states. This association is not expected if cell differentiation is, in effect, regulated by epigenetic mechanisms. Based on these results, I propose a general, falsifiable theory of individuated multicellularity, which relies on the synergistic coupling across the extracellular space of two explicitly uncorrelated "self-organizing" systems constraining histone modification states at the same sites. This theory describes how the simplest multicellular individual-understood as an intrinsic, higher-order constraint-emerges from proliferating undifferentiated cells, and could explain the intrinsic regulation of gene transcriptional changes for cell differentiation and the evolution of individuated multicellular organisms.

Positional information and reaction-diffusion: two big ideas in developmental biology combine

One of the most fundamental questions in biology is that of biological pattern: how do the structures and shapes of organisms arise? Undoubtedly, the two most influential ideas in this area are those of Alan Turing's 'reaction-diffusion' and Lewis Wolpert's 'positional information'. Much has been written about these two concepts but some confusion still remains, in particular about the relationship between them. Here, we address this relationship and propose a scheme of three distinct ways in which these two ideas work together to shape biological form.

Sea Urchin Morphogenesis

In the sea urchin morphogenesis follows extensive molecular specification. The specification controls the many morphogenetic events and these, in turn, precede patterning steps that establish the larval body plan. To understand how the embryo is built it was necessary to understand those series of molecular steps. Here an example of the historical sequence of those discoveries is presented as it unfolded over the last 50 years, the years during which major progress in understanding development of many animals and plants was documented by CTDB. In sea urchin development a rich series of experimental studies first established many of the phenomenological components of skeletal morphogenesis and patterning without knowledge of the molecular components. The many discoveries of transcription factors, signals, and structural proteins that contribute to the shape of the endoskeleton of the sea urchin larva then followed as molecular tools became available. A number of transcription factors and signals were discovered that were necessary for specification, morphogenesis, and patterning. Perturbation of the transcription factors and signals provided the means for assembling models of the gene regulatory networks used for specification and controlled the subsequent morphogenetic events. The earlier experimental information informed perturbation experiments that asked how patterning worked. As a consequence it was learned that ectoderm provides a series of patterning signals to the skeletogenic cells and as a consequence the skeletogenic cells secrete a highly patterned skeleton based on their ability to genotypically decode the localized reception of several signals. We still do not understand the complexity of the signals received by the skeletogenic cells, nor do we understand in detail how the genotypic information shapes the secreted skeletal biomineral, but the current knowledge at least outlines the sequence of events and provides a useful template for future discoveries.

Dual Monitoring of Secretion and ATP Levels during Chondrogenesis Using Perfusion Culture-Combined Bioluminescence Monitoring System

Skeletal pattern formation in limb development depends on prechondrogenic condensation which prefigures the cartilage template. However, although morphogens such as TGF-βs and BMPs have been known to play essential roles in skeletal patterning, how the morphogens induce prechondrogenic cells to aggregate and determine patterns of cartilage elements has remained unclear. Our previous study reported that ATP oscillations are induced during chondrogenesis. This result suggests the possibility that ATP oscillations lead to the oscillatory secretion of morphogens, due to the fact that secretion process requires ATP. To examine the correlation between ATP oscillations and secretion levels of morphogens, we have developed perfusion culture-combined bioluminescence monitoring system to simultaneously monitor intracellular ATP levels and secretion levels. Using this system, we found that secretory activity oscillates in phase with ATP oscillations and that secretion levels of TGF-β1 and BMP2 oscillate during chondrogenesis. The oscillatory secretion of the morphogens would contribute to amplifying the fluctuation of the morphogens, underlie the spatial patterning of morphogens, and consequently lead to skeletal pattern formation.

A newly established culture method highlights regulatory roles of retinoic acid on morphogenesis and calcification of mammalian limb cartilage

During mammalian embryogenesis, sclerotome-derived chondrocytes in the limb bud are arranged into a complicated bone shape with specific areas undergoing hypertrophy and calcification, creating a region-specific mineralized pattern in the cartilage. To follow chondrogenesis progression in vitro, we isolated limb cartilage from mice on embryonic day 13 (E13) and cultured it at the air-liquid interface after microsurgical removal of the ectoderm/epidermis. Explants underwent proper morphogenesis, giving rise to complete templates for limb bones in vitro. We found that region-specific calcification patterns resembling limbs of prepartum mature embryos could be induced in explants using culture medium containing high concentrations of CaCl2 (Ca), ascorbic acid (AA), and β-glycerophosphoric acid (BGP). In this culture system, excess amounts of all-trans retinoic acid (RA) severely disrupted morphogenesis and calcification patterns in limb cartilage. These effects were more pronounced in forearms than in phalanges. Although dissociated, the nascent chondrocytes in culture did not give rise to cartilage units even though augmented calcification was induced in these cell aggregates in the presence of RA. Taken together, our newly established culture system revealed that RA independently regulates three-dimensional morphogenesis and calcification.

The chemical nature of the factor responsible for embryonic induction: an historical overview

The discovery by Hans Spemann of the "organizer" tissue and its ability to induce the formation of the amphibian embryo's neural tube inspired leading embryologists to attempt to elucidate embryonic induction's underlying mechanisms. Since then several studies have described several developmental model system to better understand the role of specific signaling molecules, the interplay of different signals and tissue interactions in regulating tissue induction and patterning events. Different groups of workers set out to subject embryonic amphibian tissues and inductive adult organs to various extraction methods in the hope that the active agents could be isolated and chemically identified. In addition, a large number of well characterized chemical compounds were tested.

Mathematical modeling of vertebrate limb development

In this paper, we review the major mathematical and computational models of vertebrate limb development and their roles in accounting for different aspects of this process. The main aspects of limb development that have been modeled include outgrowth and shaping of the limb bud, establishment of molecular gradients within the bud, and formation of the skeleton. These processes occur interdependently during development, although (as described in this review), there are various interpretations of the biological relationships among them. A wide range of mathematical and computational methods have been used to study these processes, including ordinary and partial differential equation systems, cellular automata and discrete, stochastic models, finite difference methods, finite element methods, the immersed boundary method, and various combinations of the above. Multiscale mathematical modeling and associated computational simulation have become integrated into the study of limb morphogenesis and pattern formation to an extent with few parallels in the field of developmental biology. These methods have contributed to the design and analysis of experiments employing microsurgical and genetic manipulations, evaluation of hypotheses for limb bud outgrowth, interpretation of the effects of natural mutations, and the formulation of scenarios for the origination and evolution of the limb skeleton.

Morphogenetic fields in embryogenesis, regeneration, and cancer: non-local control of complex patterning

Establishment of shape during embryonic development, and the maintenance of shape against injury or tumorigenesis, requires constant coordination of cell behaviors toward the patterning needs of the host organism. Molecular cell biology and genetics have made great strides in understanding the mechanisms that regulate cell function. However, generalized rational control of shape is still largely beyond our current capabilities. Significant instructive signals function at long range to provide positional information and other cues to regulate organism-wide systems properties like anatomical polarity and size control. Is complex morphogenesis best understood as the emergent property of local cell interactions, or as the outcome of a computational process that is guided by a physically encoded map or template of the final goal state? Here I review recent data and molecular mechanisms relevant to morphogenetic fields: large-scale systems of physical properties that have been proposed to store patterning information during embryogenesis, regenerative repair, and cancer suppression that ultimately controls anatomy. Placing special emphasis on the role of endogenous bioelectric signals as an important component of the morphogenetic field, I speculate on novel approaches for the computational modeling and control of these fields with applications to synthetic biology, regenerative medicine, and evolutionary developmental biology.

A new method for choosing the computational cell in stochastic reaction-diffusion systems

How to choose the computational compartment or cell size for the stochastic simulation of a reaction-diffusion system is still an open problem, and a number of criteria have been suggested. A generalized measure of the noise for finite-dimensional systems based on the largest eigenvalue of the covariance matrix of the number of molecules of all species has been suggested as a measure of the overall fluctuations in a multivariate system, and we apply it here to a discretized reaction-diffusion system. We show that for a broad class of first-order reaction networks this measure converges to the square root of the reciprocal of the smallest mean species number in a compartment at the steady state. We show that a suitably re-normalized measure stabilizes as the volume of a cell approaches zero, which leads to a criterion for the maximum volume of the compartments in a computational grid. We then derive a new criterion based on the sensitivity of the entire network, not just of the fastest step, that predicts a grid size that assures that the concentrations of all species converge to a spatially-uniform solution. This criterion applies for all orders of reactions and for reaction rate functions derived from singular perturbation or other reduction methods, and encompasses both diffusing and non-diffusing species. We show that this predicts the maximal allowable volume found in a linear problem, and we illustrate our results with an example motivated by anterior-posterior pattern formation in Drosophila, and with several other examples.

Genomic differences distinguish the myofibroblast phenotype of distal lung fibroblasts from airway fibroblasts

Primary human distal lung/parenchymal fibroblasts (DLFs) exhibit a different phenotype from airway fibroblasts (AFs), including the expression of high levels of α-smooth muscle actin (α-SMA). The scope of the differences between these anatomically differentiated fibroblasts, or the mechanisms driving them, has remained unknown. To determine whether the different characteristics of regional fibroblasts are predicted by distinct genomic differences in AFs versus DLFs, matched human fibroblast pairs were isolated from proximal and distal lung tissue and evaluated. Microarray analysis was performed on 12 matched fibroblast pairs (four normal and eight asthmatic samples) and validated by quantitative real-time PCR. The potential functional implications of these differences were analyzed using computational approaches. Four hundred seventy-four transcripts were up-regulated in AFs, and 611 were up-regulated in DLFs via microarray analysis. No differences in normal and asthmatic fibroblasts were evident, and the data were combined for subsequent analyses. Gene ontology and network analyses suggested distinct patterns of pathway activation between AFs and DLFs. The up-regulation of extracellular matrix-associated molecules in AFs was observed, whereas genes associated with actin binding and cytoskeletal organization were up-regulated in DLFs. The up-regulation of activated/total SMAD3 and c-Jun N-terminal kinase in DLFs may partly explain these myofibroblast-like characteristics in DLFs. Thus, marked genomic differences exist between these two populations of regional lung fibroblasts. These striking differences may help identify potential mechanisms by which AFs and DLFs differ in their responses to injury, regeneration, and remodeling in the lung.