Patterning principles of morphogen gradients

Optogenetic control of pheromone gradients and mating behavior in budding yeast

During mating in budding yeast, cells use pheromones to locate each other and fuse. This model system has shaped our current understanding of signal transduction and cell polarization in response to extracellular signals. The cell populations producing extracellular signal landscapes themselves are, however, less well understood, yet crucial for functionally testing quantitative models of cell polarization and for controlling cell behavior through bioengineering approaches. Here we engineered optogenetic control of pheromone landscapes in mating populations of budding yeast, hijacking the mating-pheromone pathway to achieve spatial control of growth, cell morphology, cell-cell fusion, and distance-dependent gene expression in response to light. Using our tool, we were able to spatially control and shape pheromone gradients, allowing the use of a biophysical model to infer the properties of large-scale gradients generated by mating populations in a single, quantitative experimental setup, predicting that the shape of such gradients depends quantitatively on population parameters. Spatial optogenetic control of diffusible signals and their degradation provides a controllable signaling environment for engineering artificial communication and cell-fate systems in gel-embedded cell populations without the need for physical manipulation.

Origin, fate and function of extraembryonic tissues during mammalian development

Extraembryonic tissues have pivotal roles in morphogenesis and patterning of the early mammalian embryo. Developmental programmes mediated through signalling pathways and gene regulatory networks determine the sequence in which fate determination and lineage commitment of extraembryonic tissues take place, and epigenetic processes allow the memory of cell identity and state to be sustained throughout and beyond embryo development, even extending across generations. In this Review, we discuss the molecular and cellular mechanisms necessary for the different extraembryonic tissues to develop and function, from their initial specification up until the end of gastrulation, when the body plan of the embryo and the anatomical organization of its supporting extraembryonic structures are established. We examine the interaction between extraembryonic and embryonic tissues during early patterning and morphogenesis, and outline how epigenetic memory supports extraembryonic tissue development.

From genes to patterns: a framework for modeling the emergence of embryonic development from transcriptional regulation

Understanding embryonic patterning, the process by which groups of cells are partitioned into distinct identities defined by gene expression, is a central challenge in developmental biology. This complex phenomenon is driven by precise spatial and temporal regulation of gene expression across many cells, resulting in the emergence of highly organized tissue structures. While similar emergent behavior is well understood in other fields, such as statistical mechanics, the regulation of gene expression in development remains less clear, particularly regarding how molecular-level gene interactions lead to the large-scale patterns observed in embryos. In this study, we present a modeling framework that bridges the gap between molecular gene regulation and tissue-level embryonic patterning. Beginning with basic chemical reaction models of transcription at the single-gene level, we progress to model gene regulatory networks (GRNs) that mediate specific cellular functions. We then introduce phenomenological models of pattern formation, including the French Flag and Temporal Patterning/Speed Regulation models, and integrate them with molecular/GRN realizations. To facilitate understanding and application of our models, we accompany our mathematical framework with computer simulations, providing intuitive and simple code for each model. A key feature of our framework is the explicit articulation of underlying assumptions at each level of the model, from transcriptional regulation to tissue patterning. By making these assumptions clear, we provide a foundation for future experimental and theoretical work to critically examine and challenge them, thereby improving the accuracy and relevance of gene regulatory models in developmental biology. As a case study, we explore how different strategies for integrating enhancer activity affect the robustness and evolvability of GRNs that govern embryonic pattern formation. Our simulations suggest that a two-step regulation strategy, enhancer activation followed by competitive integration at the promoter, ensures more standardized integration of new enhancers into developmental GRNs, highlighting the adaptability of eukaryotic transcription. These findings shed new light on the transcriptional mechanisms underlying embryonic patterning, while the overall modeling framework serves as a foundation for future experimental and theoretical investigations.

Integrated multi-omics profiling characterizes the crucial role of human dental epithelium during tooth development

The development of early human tooth primordia is not well understood. Here, we linked single-cell RNA sequencing, spatial transcriptomics, and secretome analysis to characterize human fetal tooth development over time. A spatiotemporal atlas of human tooth development at multiple levels was mapped, identifying previously uncharacterized epithelial subpopulations with distinct gene expression profiles and spatial localization. Dynamic changes in epithelial-mesenchymal interactions across developmental stages were characterized. Secretome analysis confirmed the extensive paracrine signaling from the epithelial to mesenchymal compartments and uncovered signaling factors produced by dental epithelium (DE) that regulate mesenchymal cell fate and differentiation. Integration of these datasets highlighted the crucial role of the DE in orchestrating tooth morphogenesis. Our multi-omics approach not only provides unprecedented insights into the cellular and molecular mechanisms of ectoderm-derived tissue development but also serves as a valuable resource, which is publicly available online, for future studies on human tooth regeneration and related diseases.

Contact area and tissue growth dynamics shape synthetic juxtacrine signaling patterns

Cell-cell communication through direct contact, or juxtacrine signaling, is important in development, disease, and many areas of physiology. Synthetic forms of juxtacrine signaling can be precisely controlled and operate orthogonally to native processes, making them a powerful reductionist tool with which to address fundamental questions in cell-cell communication in vivo. Here, we investigate how cell-cell contact length and tissue growth dynamics affect juxtacrine signal responses through implementing a custom synthetic gene circuit in Drosophila wing imaginal discs alongside mathematical modeling to determine synthetic Notch (synNotch) activation patterns. We find that the area of contact between cells largely determines the extent of synNotch activation, leading to the prediction that the shape of the interface between signal-sending and signal-receiving cells will impact the magnitude of the synNotch response. Notably, synNotch outputs form a graded spatial profile that extends several cell diameters from the signal source, providing evidence that the response to juxtacrine signals can persist in cells as they proliferate away from source cells, or that cells remain able to communicate directly over several cell diameters. Our model suggests that the former mechanism may be sufficient, since it predicts graded outputs without diffusion or long-range cell-cell communication. Overall, we identify that cell-cell contact area together with output synthesis and decay rates likely govern the pattern of synNotch outputs in both space and time during tissue growth, insights that may have broader implications for juxtacrine signaling in general.

Intrinsically robust and scalable biofilm segmentation under diverse physical growth conditions

Developmental patterning is a shared feature across biological systems ranging from vertebrates to bacterial biofilms. While vertebrate patterning benefits from well-controlled homeostatic environments, bacterial biofilms can grow in diverse physical contexts. What mechanisms provide developmental robustness under diverse environments remains an open question. We show that a native clock-and-wavefront mechanism robustly segments biofilms in both solid-air and solid-liquid interfaces. Biofilms grown under these distinct physical conditions differ 4-fold in size yet exhibit robust segmentation. The segmentation pattern scaled with biofilm growth rate in a mathematically predictable manner independent of habitat conditions. We show that scaling arises from the coupling between wavefront speed and biofilm growth rate. In contrast to the complexity of scaling mechanisms in vertebrates, our data suggests that the minimal bacterial clock-and-wavefront mechanism is intrinsically robust and scales in real time. Consequently, bacterial biofilms robustly segment under diverse conditions without requiring cell-to-cell signaling to track system size.

Mobile signals, patterning, and positional information in root development

Multicellular organisms use mobile intercellular signals to generate spatiotemporal patterns of growth and differentiation. These signals, termed morphogens, arise from localized sources and move by diffusion or directional transport to be interpreted at target cells. The classical model for a morphogen is where a substance diffuses from a source to generate a concentration gradient that provides positional information across a field. This concept, presented by Wolpert and popularized as the "French Flag Model," remains highly influential, but other patterning models, which do not rely on morphogen gradients, also exist. Here, we review current evidence for mobile morphogenetic signals in plant root development and how they fit within existing conceptual frameworks for pattern formation. We discuss how the signals are formed, distributed, and interpreted in space and time, emphasizing the regulation of movement on the ability of morphogens to specify patterns. While significant advances have been made in the field since the first identification of mobile morphogenetic factors in plants, key questions remain to be answered, such as how morphogen movement is regulated, how these mechanisms allow scaling in different species, and how morphogens act to enable plant regeneration in response to damage.

Stable developmental patterns of gene expression without morphogen gradients

Gene expression patterns in developing organisms are established by groups of cross-regulating target genes that are driven by morphogen gradients. As development progresses, morphogen activity is reduced, leaving the emergent pattern without stabilizing positional cues and at risk of rapid deterioration due to the inherently noisy biochemical processes at the cellular level. But remarkably, gene expression patterns remain spatially stable and reproducible over long developmental time spans in many biological systems. Here we combine spatial-stochastic simulations with an enhanced sampling method (Non-Stationary Forward Flux Sampling) and a recently developed stability theory to address how spatiotemporal integrity of a gene expression pattern is maintained in developing tissue lacking morphogen gradients. Using a minimal embryo model consisting of spatially coupled biochemical reactor volumes, we study a prototypical stripe pattern in which weak cross-repression between nearest neighbor expression domains alternates with strong repression between next-nearest neighbor domains, inspired by the gap gene system in the Drosophila embryo. We find that tuning of the weak repressive interactions to an optimal level can prolong stability of the expression patterns by orders of magnitude, enabling stable patterns over developmentally relevant times in the absence of morphogen gradients. The optimal parameter regime found in simulations of the embryo model closely agrees with the predictions of our coarse-grained stability theory. To elucidate the origin of stability, we analyze a reduced phase space defined by two measures of pattern asymmetry. We find that in the optimal regime, intact patterns are protected via restoring forces that counteract random perturbations and give rise to a metastable basin. Together, our results demonstrate that metastable attractors can emerge as a property of stochastic gene expression patterns even without system-wide positional cues, provided that the gene regulatory interactions shaping the pattern are optimally tuned.

Interplay between Notch signaling and mechanical forces during developmental patterning processes

The coordination between biochemical signals and cell mechanics has emerged in recent years as a crucial mechanism driving developmental patterning processes across a variety of developing and homeostatic systems. An important class of such developmental processes relies on local communication between neighboring cells through Notch signaling. Here, we review how the coordination between Notch-mediated differentiation and cell mechanics can give rise to unique cellular patterns. We discuss how global and local mechanical cues can affect, and be affected by, cellular differentiation and reorganization controlled by Notch signaling. We compare recent studies of such developmental processes, including the mammalian inner ear, Drosophila ommatidia, intestinal organoids, and zebrafish myocardium, to draw shared general concepts and their broader implications in biology.

MUM, a maternal unknown message, inhibits early establishment of the medio-lateral axis in the embryo of the kelp Saccharina latissima

Brown algae are multicellular photosynthetic organisms that have evolved independently of plants and other algae. Here, we have studied the determinism of body axis formation in the kelp Saccharina latissima. After microdissection of the embryo, we show that the stalk, an empty cell that retains the embryo on the maternal tissue, represses longitudinal cell divisions in the early embryo, thereby reinforcing the establishment of the initial apico-basal axis. In addition, it promotes cell growth and controls cell shape and arrangement in the flat oblong embryo composed of cells aligned in rows and columns. Although the stalk persists for several weeks until the embryo reaches at least 500 cells, proper embryogenesis requires connection to maternal tissue only during the first 4 days after fertilisation, i.e. before the embryo reaches the 8-cell stage. Transplantation experiments indicate that the maternal signal is not diffused in seawater, but requires contact between the embryo and the maternal tissue. This first global quantitative study of brown algal embryogenesis highlights the role of MUM, an unknown maternal message, in the control of growth axes and tissue patterning in kelp embryos.

'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.

The diversification of methods for studying cell-cell interactions and communication

No cell lives in a vacuum, and the molecular interactions between cells define most phenotypes. Transcriptomics provides rich information to infer cell-cell interactions and communication, thus accelerating the discovery of the roles of cells within their communities. Such research relies heavily on algorithms that infer which cells are interacting and the ligands and receptors involved. Specific pressures on different research niches are driving the evolution of next-generation computational tools, enabling new conceptual opportunities and technological advances. More sophisticated algorithms now account for the heterogeneity and spatial organization of cells, multiple ligand types and intracellular signalling events, and enable the use of larger and more complex datasets, including single-cell and spatial transcriptomics. Similarly, new high-throughput experimental methods are increasing the number and resolution of interactions that can be analysed simultaneously. Here, we explore recent progress in cell-cell interaction research and highlight the diversification of the next generation of tools, which have yielded a rich ecosystem of tools for different applications and are enabling invaluable discoveries.

Responses of organ precursors to correct and incorrect inductive signals

During embryonic development, the inductive molecules produced by local origins normally arrive at their target tissues in a nondirectional, diffusion manner. The target organ precursor cells must correctly interpret these inductive signals to ensure proper specification/differentiation, which is dependent on two prerequisites: (i) obtaining cell-intrinsic competence; and (ii) receiving correct inductive signals while resisting incorrect ones. Gain of intrinsic competence could avoid a large number of misinductions because the incompetent cells are nonresponsive to inductive signals. However, in cases of different precursor cells with similar competence and located in close proximity, resistance to incorrect inductive signals is essential for accurate determination of cell fate. Here we outline the mechanisms of how organ precursors respond to correct and incorrect inductive signals.

Contact area and tissue growth dynamics shape synthetic juxtacrine signaling patterns

Cell-cell communication through direct contact, or juxtacrine signaling, is important in development, disease, and many areas of physiology. Synthetic forms of juxtacrine signaling can be precisely controlled and operate orthogonally to native processes, making them a powerful reductionist tool with which to address fundamental questions in cell-cell communication in vivo. Here we investigate how cell-cell contact length and tissue growth dynamics affect juxtacrine signal responses through implementing a custom synthetic gene circuit in Drosophila wing imaginal discs alongside mathematical modeling to determine synthetic Notch (synNotch) activation patterns. We find that the area of contact between cells largely determines the extent of synNotch activation, leading to the prediction that the shape of the interface between signal-sending and signal-receiving cells will impact the magnitude of the synNotch response. Notably, synNotch outputs form a graded spatial profile that extends several cell diameters from the signal source, providing evidence that the response to juxtacrine signals can persist in cells as they proliferate away from source cells, or that cells remain able to communicate directly over several cell diameters. Our model suggests the former mechanism may be sufficient, since it predicts graded outputs without diffusion or long-range cell-cell communication. Overall, we identify that cell-cell contact area together with output synthesis and decay rates likely govern the pattern of synNotch outputs in both space and time during tissue growth, insights that may have broader implications for juxtacrine signaling in general.

Oscillatory control of embryonic development

Proper embryonic development depends on the timely progression of a genetic program. One of the key mechanisms for achieving precise control of developmental timing is to use gene expression oscillations. In this Review, we examine how gene expression oscillations encode temporal information during vertebrate embryonic development by discussing the gene expression oscillations occurring during somitogenesis, neurogenesis, myogenesis and pancreas development. These oscillations play important but varied physiological functions in different contexts. Oscillations control the period of somite formation during somitogenesis, whereas they regulate the proliferation-to-differentiation switch of stem cells and progenitor cells during neurogenesis, myogenesis and pancreas development. We describe the similarities and differences of the expression pattern in space (i.e. whether oscillations are synchronous or asynchronous across neighboring cells) and in time (i.e. different time scales) of mammalian Hes/zebrafish Her genes and their targets in different tissues. We further summarize experimental evidence for the functional role of their oscillations. Finally, we discuss the outstanding questions for future research.

Validation of a mouse 3D gastruloid-based embryotoxicity assay in reference to the ICH S5(R3) guideline chemical exposure list

There is growing interest in establishing alternative methods in place of conventional animal tests to assess the developmental and reproductive toxicity (DART) of chemicals. Gastruloids are 3D aggregates of pluripotent stem cells that spontaneously exhibit axial elongation morphogenesis similar to gastrulation. They have been explored as in vitro embryogenesis models for developmental and toxicological studies. Here, a mouse gastruloid-based assay was validated for DART assessment in accordance with the ICH S5(R3) guideline, which provides the plasma concentration data of various reference drugs in rodents, specifically Cmax and AUC for NOAEL and LOAEL. First, adverse effect concentrations of the reference drugs and their known metabolites on gastruloid development were determined based on morphological impact, namely reduced growth or aberrant elongation. Then, the NOAEL to LOAEL concentration range obtained from the gastruloid assay was compared with that in rodents to examine similarities in sensitivity between the in vitro and in vivo assays for each chemical. For 18 out of the 24 reference drugs that have both NOAEL and LOAEL information in rodents, the sensitivity of the gastruloid assay was comparable to the in vivo assay within an 8-fold concentration margin. For 7 out of the 8 additional reference drugs that have only NOAEL or LOAEL information in rodents, the gastruloid assay was in line with the in vivo data. Altogether, these results support the effectiveness of the gastruloid assay, which may be exploited as a non-animal alternative method for DART assessment.

Bone morphogenetic protein signaling: the pathway and its regulation

In the mid-1960s, bone morphogenetic proteins (BMPs) were first identified in the extracts of bone to have the remarkable ability to induce heterotopic bone. When the Drosophila gene decapentaplegic (dpp) was first identified to share sequence similarity with mammalian BMP2/BMP4 in the late-1980s, it became clear that secreted BMP ligands can mediate processes other than bone formation. Following this discovery, collaborative efforts between Drosophila geneticists and mammalian biochemists made use of the strengths of their respective model systems to identify BMP signaling components and delineate the pathway. The ability to conduct genetic modifier screens in Drosophila with relative ease was critical in identifying the intracellular signal transducers for BMP signaling and the related transforming growth factor-beta/activin signaling pathway. Such screens also revealed a host of genes that encode other core signaling components and regulators of the pathway. In this review, we provide a historical account of this exciting time of gene discovery and discuss how the field has advanced over the past 30 years. We have learned that while the core BMP pathway is quite simple, composed of 3 components (ligand, receptor, and signal transducer), behind the versatility of this pathway lies multiple layers of regulation that ensures precise tissue-specific signaling output. We provide a sampling of these discoveries and highlight many questions that remain to be answered to fully understand the complexity of BMP signaling.

Spatiotemporal control of pattern formation during somitogenesis

Spatiotemporal patterns widely occur in biological, chemical, and physical systems. Particularly, embryonic development displays a diverse gamut of repetitive patterns established in many tissues and organs. Branching treelike structures in lungs, kidneys, livers, pancreases, and mammary glands as well as digits and bones in appendages, teeth, and palates are just a few examples. A fascinating instance of repetitive patterning is the sequential segmentation of the primary body axis, which is conserved in all vertebrates and many arthropods and annelids. In these species, the body axis elongates at the posterior end of the embryo containing an unsegmented tissue. Meanwhile, segments sequentially bud off from the anterior end of the unsegmented tissue, laying down an exquisite repetitive pattern and creating a segmented body plan. In vertebrates, the paraxial mesoderm is sequentially divided into somites. In this review, we will discuss the most prominent models, the most puzzling experimental data, and outstanding questions in vertebrate somite segmentation.

Information integration during bioelectric regulation of morphogenesis of the embryonic frog brain

Spatiotemporal patterns of cellular resting potential regulate several aspects of development. One key aspect of the bioelectric code is that transcriptional and morphogenetic states are determined not by local, single-cell, voltage levels but by specific distributions of voltage across cell sheets. We constructed and analyzed a minimal dynamical model of collective gene expression in cells based on inputs of multicellular voltage patterns. Causal integration analysis revealed a higher-order mechanism by which information about the voltage pattern was spatiotemporally integrated into gene activity, as well as a division of labor among and between the bioelectric and genetic components. We tested and confirmed predictions of this model in a system in which bioelectric control of morphogenesis regulates gene expression and organogenesis: the embryonic brain of the frog Xenopus laevis. This study demonstrates that machine learning and computational integration approaches can advance our understanding of the information-processing underlying morphogenetic decision-making, with a potential for other applications in developmental biology and regenerative medicine.

The salt-and-pepper pattern in mouse blastocysts is compatible with signaling beyond the nearest neighbors

Embryos develop in a concerted sequence of spatiotemporal arrangements of cells. In the preimplantation mouse embryo, the distribution of the cells in the inner cell mass evolves from a salt-and-pepper pattern to spatial segregation of two distinct cell types. The exact properties of the salt-and-pepper pattern have not been analyzed so far. We investigate the spatiotemporal distribution of NANOG- and GATA6-expressing cells in the ICM of the mouse blastocysts with quantitative three-dimensional single-cell-based neighborhood analyses. A combination of spatial statistics and agent-based modeling reveals that the cell fate distribution follows a local clustering pattern. Using ordinary differential equations modeling, we show that this pattern can be established by a distance-based signaling mechanism enabling cells to integrate information from the whole inner cell mass into their cell fate decision. Our work highlights the importance of longer-range signaling to ensure coordinated decisions in groups of cells to successfully build embryos.

Closing the loop on morphogenesis: a mathematical model of morphogenesis by closed-loop reaction-diffusion

Morphogenesis, the establishment and repair of emergent complex anatomy by groups of cells, is a fascinating and biomedically-relevant problem. One of its most fascinating aspects is that a developing embryo can reliably recover from disturbances, such as splitting into twins. While this reliability implies some type of goal-seeking error minimization over a morphogenic field, there are many gaps with respect to detailed, constructive models of such a process. A common way to achieve reliability is negative feedback, which requires characterizing the existing body shape to create an error signal-but measuring properties of a shape may not be simple. We show how cells communicating in a wave-like pattern could analyze properties of the current body shape. We then describe a closed-loop negative-feedback system for creating reaction-diffusion (RD) patterns with high reliability. Specifically, we use a wave to count the number of peaks in a RD pattern, letting us use a negative-feedback controller to create a pattern with N repetitions, where N can be altered over a wide range. Furthermore, the individual repetitions of the RD pattern can be easily stretched or shrunk under genetic control to create, e.g., some morphological features larger than others. This work contributes to the exciting effort of understanding design principles of morphological computation, which can be used to understand evolved developmental mechanisms, manipulate them in regenerative-medicine settings, or engineer novel synthetic morphology constructs with desired robust behavior.

Zebrafish as a Model to Study Retinoic Acid Signaling in Development and Disease

Retinoic acid (RA) is a metabolite of vitamin A (retinol) that plays various roles in development to influence differentiation, patterning, and organogenesis. RA also serves as a crucial homeostatic regulator in adult tissues. The role of RA and its associated pathways are well conserved from zebrafish to humans in both development and disease. This makes the zebrafish a natural model for further interrogation into the functions of RA and RA-associated maladies for the sake of basic research, as well as human health. In this review, we explore both foundational and recent studies using zebrafish as a translational model for investigating RA from the molecular to the organismal scale.