Differential diffusivity of Nodal and Lefty underlies a reaction-diffusion patterning system

How a reaction-diffusion signal can control spinal cord regeneration in axolotls: A modeling study

Axolotls are uniquely able to completely regenerate the spinal cord after amputation. The underlying governing mechanisms of this regenerative response have not yet been fully elucidated. We previously found that spinal cord regeneration is mainly driven by cell-cycle acceleration of ependymal cells, recruited by a hypothetical signal propagating from the injury. However, the nature of the signal and its propagation remain unknown. In this theoretical study, we investigated whether the regeneration-inducing signal can follow a reaction-diffusion process. We developed a computational model, validated it with experimental data, and showed that the signal dynamics can be understood in terms of reaction-diffusion mechanism. By developing a theory of the regenerating outgrowth in the limit of fast reaction-diffusion, we demonstrate that control of regenerative response solely relies on cell-to-signal sensitivity and the signal reaction-diffusion characteristic length. This study lays foundations for further identification of the signal controlling regeneration of the spinal cord.

The evolution of developmental biology through conceptual and technological revolutions

Developmental biology-the study of the processes by which cells, tissues, and organisms develop and change over time-has entered a new golden age. After the molecular genetics revolution in the 80s and 90s and the diversification of the field in the early 21st century, we have entered a phase when powerful technologies provide new approaches and open unexplored avenues. Progress in the field has been accelerated by advances in genomics, imaging, engineering, and computational biology and by emerging model systems ranging from tardigrades to organoids. We summarize how revolutionary technologies have led to remarkable progress in understanding animal development. We describe how classic questions in gene regulation, pattern formation, morphogenesis, organogenesis, and stem cell biology are being revisited. We discuss the connections of development with evolution, self-organization, metabolism, time, and ecology. We speculate how developmental biology might evolve in an era of synthetic biology, artificial intelligence, and human engineering.

Maternal TGF-β ligand Panda breaks the radial symmetry of the sea urchin embryo by antagonizing the Nodal type II receptor ACVRII

In the highly regulative embryo of the sea urchin Paracentrotus lividus, establishment of the dorsal-ventral (D/V) axis critically depends on the zygotic expression of the TGF-β nodal in the ventral ectoderm. nodal expression is first induced ubiquitously in the 32-cell embryo and becomes progressively restricted to the presumptive ventral ectoderm by the early blastula stage. This early spatial restriction of nodal expression is independent of Lefty, and instead relies on the activity of Panda, a maternally expressed TGF-β ligand related to Lefty and Inhibins, which is required maternally for D/V axis specification. However, the mechanism by which Panda restricts the early nodal expression has remained enigmatic and it is not known if Panda works like a BMP ligand by opposing Nodal and antagonizing Smad2/3 signaling, or if it works like Lefty by sequestering an essential component of the Nodal signaling pathway. In this study, we report that Panda functions as an antagonist of the TGF-β type II receptor ACVRII (Activin receptor type II), which is the only type II receptor for Nodal signaling in the sea urchin and is also a type II receptor for BMP ligands. Inhibiting translation of acvrII mRNA disrupted D/V patterning across all 3 germ layers and caused acvrII morphants to develop with a typical Nodal loss-of-function phenotype. In contrast, embryos overexpressing acvrII displayed strong ectopic Smad1/5/8 signaling at blastula stages and developed as dorsalized larvae, a phenotype very similar to that caused by over activation of BMP signaling. Remarkably, embryos co-injected with acvrII mRNA and panda mRNA did not show ectopic Smad1/5/8 signaling and developed with a largely normal dorsal-ventral polarity. Furthermore, using an axis induction assay, we found that Panda blocks the ability of ACVRII to orient the D/V axis when overexpressed locally. Using co-immunoprecipitation, we showed that Panda physically interacts with ACVRII, as well as with the Nodal co-receptor Cripto, and with TBR3 (Betaglycan), which is a non-signaling receptor for Inhibins in mammals. At the molecular level, we have traced back the antagonistic activity of Panda to the presence of a single proline residue, conserved with all the Lefty factors, in the ACVRII binding motif of Panda, instead of a serine as in most of TGF-β ligands. Conversion of this proline to a serine converted Panda from an antagonist that opposed Nodal signaling and promoted dorsalization to an agonist that promoted Nodal signaling and triggered ventralization when overexpressed. Finally, using phylogenomics, we analyzed the emergence of the agonist and antagonist form of Panda in the course of evolution. Our data are consistent with the idea that the presence of a serine at that position, like in most TGF-β, was the ancestral condition and that the initial function of Panda was possibly in promoting and not in antagonizing Nodal signaling. These results highlight the existence of key functional and structural elements conserved between Panda and Lefty, allow to draw an intriguing parallel between sea urchin Panda and mammalian Inhibin α and raise the unexpected possibility that the original function of Panda may have been in activation of the Nodal pathway rather than in its inhibition.

Understanding and manipulating extracellular behaviors of Wnt ligands

Wnt, a family of secreted signaling proteins, serves diverse functions in embryogenesis, organogenesis, cancer, and stem cell functions. In the context of development, Wnt has been considered a representative morphogen, forming concentration gradients to give positional information to cells or tissues. However, although gradients are often illustrated in schemata, the reality of concentration gradients, or in other words, actual spatial distribution of Wnt ligands, and their behaviors in the extracellular space still remain poorly known. To understand extracellular behavior of Wnt ligands, quantitative analyses such as fluorescence correlation spectroscopy (FCS) and fluorescence recovery after photobleaching (FRAP) are highly informative because Wnt dispersal involves physical and biochemical processes, such as diffusion and binding to or dissociation from cell surface molecules, including heparan sulfate proteoglycans (HSPGs). Here, I briefly discuss representative methods to quantify morphogen dynamics. In addition, I discuss molecular manipulations of morphogens, mainly focusing on use of protein binders, and synthetic biology of morphogens as indicators of current and future directions in this field.

Optogenetic control of Nodal signaling patterns

A crucial step in early embryogenesis is the establishment of spatial patterns of signaling activity. Tools to perturb morphogen signals with high resolution in space and time can help reveal how embryonic cells decode these signals to make appropriate fate decisions. Here, we present new optogenetic reagents and an experimental pipeline for creaHng designer Nodal signaling patterns in live zebrafish embryos. Nodal receptors were fused to the light-sensitive heterodimerizing pair Cry2/CIB1N, and the Type II receptor was sequestered to the cytosol. The improved optoNodal2 reagents eliminate dark activity and improve response kinetics, without sacrificing dynamic range. We adapted an ultra-widefield microscopy platform for parallel light patterning in up to 36 embryos and demonstrated precise spatial control over Nodal signaling activity and downstream gene expression. Patterned Nodal activation drove precisely controlled internalization of endodermal precursors. Further, we used patterned illumination to generate synthetic signaling patterns in Nodal signaling mutants, rescuing several characteristic developmental defects. This study establishes an experimental toolkit for systematic exploration of Nodal signaling patterns in live embryos.

A diffusible small-RNA-based Turing system dynamically coordinates organ polarity

The formation of a flat and thin leaf presents a developmentally challenging problem, requiring intricate regulation of adaxial-abaxial (top-bottom) polarity. The patterning principles controlling the spatial arrangement of these domains during organ growth have remained unclear. Here we show that this regulation in Arabidopsis thaliana is achieved by an organ-autonomous Turing reaction-diffusion system centred on mobile small RNAs. The data illustrate how Turing dynamics transiently instructed by prepatterned information is sufficient to self-sustain properly oriented polarity in a dynamic, growing organ, presenting intriguing parallels to left-right patterning in the vertebrate embryo. Computational modelling demonstrates that this self-organizing system continuously adapts to coordinate the robust planar polarity of a flat leaf while affording flexibility to generate the tissue patterns of evolutionarily diverse organ shapes. Our findings identify a small-RNA-based Turing network as a dynamic regulator of organ polarity that accounts for leaf shape diversity at the level of the individual organ, plant or species.

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.

Long-range formation of the Bicoid gradient requires multiple dynamic modes that spatially vary across the embryo

Morphogen gradients provide essential positional information to gene networks through their spatially heterogeneous distribution, yet how they form is still hotly contested, with multiple models proposed for different systems. Here, we focus on the transcription factor Bicoid (Bcd), a morphogen that forms an exponential gradient across the anterior-posterior (AP) axis of the early Drosophila embryo. Using fluorescence correlation spectroscopy we find there are spatial differences in Bcd diffusivity along the AP axis, with Bcd diffusing more rapidly in the posterior. We establish that such spatially varying differences in Bcd dynamics are sufficient to explain how Bcd can have a steep exponential gradient in the anterior half of the embryo and yet still have an observable fraction of Bcd near the posterior pole. In the nucleus, we demonstrate that Bcd dynamics are impacted by binding to DNA. Addition of the Bcd homeodomain to eGFP::NLS qualitatively replicates the Bcd concentration profile, suggesting this domain regulates Bcd dynamics. Our results reveal how a long-range gradient can form while retaining a steep profile through much of its range.

Evolution of the gene regulatory network of body axis by enhancer hijacking in amphioxus

A central goal of evolutionary developmental biology is to decipher the evolutionary pattern of gene regulatory networks (GRNs) that control embryonic development, and the mechanism underlying GRNs evolution. The Nodal signaling that governs the body axes of deuterostomes exhibits a conserved GRN orchestrated principally by Nodal, Gdf1/3, and Lefty. Here we show that this GRN has been rewired in cephalochordate amphioxus. We found that while the amphioxus Gdf1/3 ortholog exhibited nearly no embryonic expression, its duplicate Gdf1/3-like, linked to Lefty, was zygotically expressed in a similar pattern as Lefty. Consistent with this, while Gdf1/3-like mutants showed defects in axial development, Gdf1/3 mutants did not. Further transgenic analyses showed that the intergenic region between Gdf1/3-like and Lefty could drive reporter gene expression as that of the two genes. These results indicated that Gdf1/3-like has taken over the axial development role of Gdf1/3 in amphioxus, possibly through hijacking Lefty enhancers. We finally demonstrated that, to compensate for the loss of maternal Gdf1/3 expression, Nodal has become an indispensable maternal factor in amphioxus and its maternal mutants caused axial defects as Gdf1/3-like mutants. We therefore demonstrated a case that the evolution of GRNs could be triggered by enhancer hijacking events. This pivotal event has allowed the emergence of a new GRN in extant amphioxus, presumably through a stepwise process. In addition, the co-expression of Gdf1/3-like and Lefty achieved by a shared regulatory region may have provided robustness during body axis formation, which provides a selection-based hypothesis for the phenomena called developmental system drift.

Quantitative Live Imaging of Zebrafish Scale Regeneration: From Adult Fish to Signaling Patterns and Tissue Flows

In regeneration, a damaged body part grows back to its original form. Understanding the mechanisms and physical principles underlying this process has been limited by the difficulties of visualizing cell signals and behaviors in regeneration. Zebrafish scales are emerging as a model system to investigate morphogenesis during vertebrate regeneration using quantitative live imaging. Scales are millimeter-sized dermal bone disks forming a skeletal armor on the body of the fish. The scale bone is deposited by an adjacent monolayer of osteoblasts that, after scale loss, regenerates in about 2 weeks. This intriguing regenerative process is accessible to live confocal microscopy, quantifications, and mathematical modeling. Here, I describe methods to image scale regeneration live, tissue-wide and at sub-cellular resolution. Furthermore, I describe methods to process the resulting images and quantify cell, tissue, and signal dynamics.

Uncovering developmental time and tempo using deep learning

During animal development, embryos undergo complex morphological changes over time. Differences in developmental tempo between species are emerging as principal drivers of evolutionary novelty, but accurate description of these processes is very challenging. To address this challenge, we present here an automated and unbiased deep learning approach to analyze the similarity between embryos of different timepoints. Calculation of similarities across stages resulted in complex phenotypic fingerprints, which carry characteristic information about developmental time and tempo. Using this approach, we were able to accurately stage embryos, quantitatively determine temperature-dependent developmental tempo, detect naturally occurring and induced changes in the developmental progression of individual embryos, and derive staging atlases for several species de novo in an unsupervised manner. Our approach allows us to quantify developmental time and tempo objectively and provides a standardized way to analyze early embryogenesis.

A chemo-mechanical model of endoderm movements driving elongation of the amniote hindgut

Although mechanical and biochemical descriptions of development are each essential, integration of upstream morphogenic cues with downstream tissue mechanics remains understudied during vertebrate morphogenesis. Here, we developed a two-dimensional chemo-mechanical model to investigate how mechanical properties of the endoderm and transport properties of fibroblast growth factor (FGF) regulate avian hindgut morphogenesis in a coordinated manner. Posterior endoderm cells convert a gradient of FGF ligands into a contractile force gradient, leading to a force imbalance that drives collective cell movements that elongate the forming hindgut tube. We formulated a 2D reaction-diffusion-advection model describing the formation of an FGF protein gradient as a result of posterior displacement of cells transcribing unstable Fgf8 mRNA during axis elongation, coupled with translation, diffusion and degradation of FGF protein. The endoderm was modeled as an active viscous fluid that generates contractile stresses in proportion to FGF concentration. With parameter values constrained by experimental data, the model replicates key aspects of hindgut morphogenesis, suggests that graded isotropic contraction is sufficient to generate large anisotropic cell movements, and provides new insight into how chemo-mechanical coupling across the mesoderm and endoderm coordinates hindgut elongation with axis elongation.

Diffusiophoresis-enhanced Turing patterns

Turing patterns are fundamental in biophysics, emerging from short-range activation and long-range inhibition processes. However, their paradigm is based on diffusive transport processes that yield patterns with shallower gradients than those observed in nature. A complete physical description of this discrepancy remains unknown. We propose a solution to this phenomenon by investigating the role of diffusiophoresis, which is the propulsion of colloids by a chemical gradient, in Turing patterns. Diffusiophoresis enables robust patterning of colloidal particles with substantially finer length scales than the accompanying chemical Turing patterns. A scaling analysis and a comparison to recent experiments indicate that chromatophores, ubiquitous in biological pattern formation, are likely diffusiophoretic and the colloidal Péclet number controls the pattern enhancement. This discovery suggests that important features of biological pattern formation can be explained with a universal mechanism that is quantified straightforwardly from the fundamental physics of colloids.

Molecular mechanisms controlling the biogenesis of the TGF-β signal Vg1

The TGF-beta signals Vg1 (Dvr1/Gdf3) and Nodal form heterodimers to induce vertebrate mesendoderm. The Vg1 proprotein is a monomer retained in the endoplasmic reticulum (ER) and is processed and secreted upon heterodimerization with Nodal, but the mechanisms underlying Vg1 biogenesis are largely elusive. Here, we clarify the mechanisms underlying Vg1 retention, processing, secretion, and signaling and introduce a Synthetic Processing (SynPro) system that enables the programmed cleavage of ER-resident and extracellular proteins. First, we find that Vg1 can be processed by intra- or extracellular proteases. Second, Vg1 can be processed without Nodal but requires Nodal for secretion and signaling. Third, Vg1-Nodal signaling activity requires Vg1 processing, whereas Nodal can remain unprocessed. Fourth, Vg1 employs exposed cysteines, glycosylated asparagines, and BiP chaperone-binding motifs for monomer retention in the ER. These observations suggest two mechanisms for rapid mesendoderm induction: Chaperone-binding motifs help store Vg1 as an inactive but ready-to-heterodimerize monomer in the ER, and the flexibility of Vg1 processing location allows efficient generation of active heterodimers both intra- and extracellularly. These results establish SynPro as an in vivo processing system and define molecular mechanisms and motifs that facilitate the generation of active TGF-beta heterodimers.

Loss of TJP1 disrupts gastrulation patterning and increases differentiation toward the germ cell lineage in human pluripotent stem cells

Biological patterning events that occur early in development establish proper tissue morphogenesis. Identifying the mechanisms that guide these patterning events is necessary in order to understand the molecular drivers of development and disease and to build tissues in vitro. In this study, we use an in vitro model of gastrulation to study the role of tight junctions and apical/basolateral polarity in modulating bone morphogenic protein-4 (BMP4) signaling and gastrulation-associated patterning in colonies of human pluripotent stem cells (hPSCs). Disrupting tight junctions via knockdown (KD) of the scaffolding tight junction protein-1 (TJP1, also known as ZO1) allows BMP4 to robustly and ubiquitously activate pSMAD1/5 signaling over time, resulting in loss of the patterning phenotype and marked differentiation bias of pluripotent stem cells to primordial germ cell-like cells (PGCLCs). These findings give important insights into how signaling events are regulated and lead to spatial emergence of diverse cell types in vitro.

The Activation of the LIMK/Cofilin Signaling Pathway via Extracellular Matrix-Integrin Interactions Is Critical for the Generation of Mature and Vascularized Cardiac Organoids

The generation of mature and vascularized human pluripotent stem cell-derived cardiac organoids (hPSC-COs) is necessary to ensure the validity of drug screening and disease modeling. This study investigates the effects of cellular aggregate (CA) stemness and self-organization on the generation of mature and vascularized hPSC-COs and elucidates the mechanisms underlying cardiac organoid (CO) maturation and vascularization. COs derived from 2-day-old CAs with high stemness (H-COs) and COs derived from 5-day-old CAs with low stemness (L-COs) were generated in a self-organized microenvironment via Wnt signaling induction. This study finds that H-COs exhibit ventricular, structural, metabolic, and functional cardiomyocyte maturation and vessel networks consisting of endothelial cells, smooth muscle cells, pericytes, and basement membranes compared to L-COs. Transcriptional profiling shows the upregulation of genes associated with cardiac maturation and vessel formation in H-COs compared with the genes in L-COs. Through experiments with LIMK inhibitors, the activation of ROCK-LIMK-pCofilin via ECM-integrin interactions leads to cardiomyocyte maturation and vessel formation in H-COs. Furthermore, the LIMK/Cofilin signaling pathway induces TGFβ/NODAL and PDGF pathway activation for the maturation and vascularization of H-COs. The study demonstrates for the first time that LIMK/Cofilin axis activation plays an important role in the generation of mature and vascularized COs.

A spatiotemporal barrier formed by Follistatin is required for left-right patterning

How left-right (LR) asymmetry emerges in a patterning field along the anterior-posterior axis remains an unresolved problem in developmental biology. Left-biased Nodal emanating from the LR organizer propagates from posterior to anterior (PA) and establishes the LR pattern of the whole embryo. However, little is known about the regulatory mechanism of the PA spread of Nodal and its asymmetric activation in the forebrain. Here, we identify bilaterally expressed Follistatin (Fst) as a regulator blocking the propagation of the zebrafish Nodal ortholog Southpaw (Spaw) in the right lateral plate mesoderm (LPM), and restricting Spaw transmission in the left LPM to facilitate the establishment of a robust LR asymmetric Nodal patterning. In addition, Fst inhibits the Activin-Nodal signaling pathway in the forebrain thus preventing Nodal activation prior to the arrival, at a later time, of Spaw emanating from the left LPM. This contributes to the orderly propagation of asymmetric Nodal activation along the PA axis. The LR regulation function of Fst is further confirmed in chick and frog embryos. Overall, our results suggest that a robust LR patterning emerges by counteracting a Fst barrier formed along the PA axis.

Recording morphogen signals reveals origins of gastruloid symmetry breaking

When cultured in three dimensional spheroids, mammalian stem cells can reproducibly self-organize a single anterior-posterior axis and sequentially differentiate into structures resembling the primitive streak and tailbud. Whereas the embryo's body axes are instructed by spatially patterned extra-embryonic cues, it is unknown how these stem cell gastruloids break symmetry to reproducibly define a single anterior-posterior (A-P) axis. Here, we use synthetic gene circuits to trace how early intracellular signals predict cells' future anterior-posterior position in the gastruloid. We show that Wnt signaling evolves from a homogeneous state to a polarized state, and identify a critical 6-hour time period when single-cell Wnt activity predicts future cellular position, prior to the appearance of polarized signaling patterns or morphology. Single-cell RNA sequencing and live-imaging reveal that early Wnt-high and Wnt-low cells contribute to distinct cell types and suggest that axial symmetry breaking is driven by sorting rearrangements involving differential cell adhesion. We further extend our approach to other canonical embryonic signaling pathways, revealing that even earlier heterogeneity in TGFβ signaling predicts A-P position and modulates Wnt signaling during the critical time period. Our study reveals a sequence of dynamic cellular processes that transform a uniform cell aggregate into a polarized structure and demonstrates that a morphological axis can emerge out of signaling heterogeneity and cell movements even in the absence of exogenous patterning cues.

A chemo-mechanical model of endoderm movements driving elongation of the amniote hindgut

While mechanical and biochemical descriptions of development are each essential, integration of upstream morphogenic cues with downstream tissue mechanics remains understudied in many contexts during vertebrate morphogenesis. A posterior gradient of Fibroblast Growth Factor (FGF) ligands generates a contractile force gradient in the definitive endoderm, driving collective cell movements to form the hindgut. Here, we developed a two-dimensional chemo-mechanical model to investigate how mechanical properties of the endoderm and transport properties of FGF coordinately regulate this process. We began by formulating a 2-D reaction-diffusion-advection model that describes the formation of an FGF protein gradient due to posterior displacement of cells transcribing unstable Fgf8 mRNA during axis elongation, coupled with translation, diffusion, and degradation of FGF protein. This was used together with experimental measurements of FGF activity in the chick endoderm to inform a continuum model of definitive endoderm as an active viscous fluid that generates contractile stresses in proportion to FGF concentration. The model replicated key aspects of hindgut morphogenesis, confirms that heterogeneous - but isotropic - contraction is sufficient to generate large anisotropic cell movements, and provides new insight into how chemo-mechanical coupling across the mesoderm and endoderm coordinates hindgut elongation with outgrowth of the tailbud.

Direct BMP signaling to chordoblasts is required for the initiation of segmented notochord sheath mineralization in zebrafish vertebral column development

The vertebral column, with the centra as its iteratively arranged building blocks, represents the anatomical key feature of the vertebrate phylum. In contrast to amniotes, where vertebrae are formed from chondrocytes and osteoblasts deriving from the segmentally organized neural crest or paraxial sclerotome, teleost vertebral column development is initiated by chordoblasts of the primarily unsegmented axial notochord, while sclerotomal cells only contribute to later steps of vertebrae formation. Yet, for both mammalian and teleostean model systems, unrestricted signaling by Bone Morphogenetic Proteins (BMPs) or retinoic acid (RA) has been reported to cause fusions of vertebral elements, while the interplay of the two signaling processes and their exact cellular targets remain largely unknown. Here, we address this interplay in zebrafish, identifying BMPs as potent and indispensable factors that, as formerly shown for RA, directly signal to notochord epithelial cells/chordoblasts to promote entpd5a expression and thereby metameric notochord sheath mineralization. In contrast to RA, however, which promotes sheath mineralization at the expense of further collagen secretion and sheath formation, BMP defines an earlier transitory stage of chordoblasts, characterized by sustained matrix production/col2a1 expression and concomitant matrix mineralization/entpd5a expression. BMP-RA epistasis analyses further indicate that RA can only affect chordoblasts and their further progression to merely mineralizing cells after they have received BMP signals to enter the transitory col2a1/entpd5a double-positive stage. This way, both signals ensure consecutively for proper mineralization of the notochord sheath within segmented sections along its anteroposterior axis. Our work sheds further light onto the molecular mechanisms that orchestrate early steps of vertebral column segmentation in teleosts. Similarities and differences to BMP's working mechanisms during mammalian vertebral column formation and the pathomechanisms underlying human bone diseases such as Fibrodysplasia Ossificans Progressiva (FOP) caused by constitutively active BMP signaling are discussed.

EmbryoNet: using deep learning to link embryonic phenotypes to signaling pathways

Evolutionarily conserved signaling pathways are essential for early embryogenesis, and reducing or abolishing their activity leads to characteristic developmental defects. Classification of phenotypic defects can identify the underlying signaling mechanisms, but this requires expert knowledge and the classification schemes have not been standardized. Here we use a machine learning approach for automated phenotyping to train a deep convolutional neural network, EmbryoNet, to accurately identify zebrafish signaling mutants in an unbiased manner. Combined with a model of time-dependent developmental trajectories, this approach identifies and classifies with high precision phenotypic defects caused by loss of function of the seven major signaling pathways relevant for vertebrate development. Our classification algorithms have wide applications in developmental biology and robustly identify signaling defects in evolutionarily distant species. Furthermore, using automated phenotyping in high-throughput drug screens, we show that EmbryoNet can resolve the mechanism of action of pharmaceutical substances. As part of this work, we freely provide more than 2 million images that were used to train and test EmbryoNet.

Gradient to sectioning CUBE workflow for the generation and imaging of organoids with localized differentiation

Advancements in organoid culture have led to various in vitro mini-organs that mimic native tissues in many ways. Yet, the bottleneck remains to generate complex organoids with body axis patterning, as well as keeping the orientation of organoids during post-experiment analysis processes. Here, we present a workflow for culturing organoids with morphogen gradient using a CUBE culture device, followed by sectioning samples with the CUBE to retain information on gradient direction. We show that hiPSC spheroids cultured with two separated differentiation media on opposing ends of the CUBE resulted in localized expressions of the respective differentiation markers, in contrast to homogeneous distribution of markers in controls. We also describe the processes for cryo and paraffin sectioning of spheroids in CUBE to retain gradient orientation information. This workflow from gradient culture to sectioning with CUBE can provide researchers with a convenient tool to generate increasingly complex organoids and study their developmental processes in vitro.

YAP1 is essential for self-organized differentiation of pluripotent stem cells

Induced pluripotent stem cells (iPSCs) form aggregates that recapitulate aspects of the self-organization in early embryogenesis. Within few days, cells undergo a transition from epithelial-like structures to organized three-dimensional embryoid bodies (EBs) with upregulation of germ layer-specific genes. However, it is largely unclear, which signaling cascades regulate self-organized differentiation. The Yes-associated protein 1 (YAP1) is a downstream effector of the Hippo pathway and essential mechanotransducer. YAP1 has been suggested to play a crucial role for early embryo development, but the relevance for early germ layer commitment of human iPSCs remains to be elucidated. To gain insights into the function of YAP1 in early cell-fate decisions, we generated YAP1 knockout (YAP^-/-) iPSC lines with CRISPR/Cas9 technology and analyzed transcriptomic and epigenetic modifications. YAP^-/- iPSCs showed increased expression of several YAP1 targets and of NODAL, an important regulator of cell differentiation. Furthermore, YAP1 deficiency evoked global DNA methylation changes. Directed differentiation of adherent iPSC colonies towards endoderm, mesoderm, and ectoderm could be induced, albeit endodermal and ectodermal differentiation showed transcriptomic and epigenetic changes in YAP^-/- lines. Notably, in undirected self-organized YAP^-/- EBs germ layer specification was clearly impaired. This phenotype was rescued via lentiviral overexpression of YAP1 and also by NODAL inhibitors. Our results demonstrate that YAP1 plays an important role during early germ layer specification of iPSCs, particularly for the undirected self-organization of EBs, and this is at least partly attributed to activation of the NODAL signaling.

Tissue interplay during morphogenesis

The process by which biological systems such as cells, tissues and organisms acquire shape has been named as morphogenesis and it is central to a plethora of biological contexts including embryo development, wound healing, or even cancer. Morphogenesis relies in both self-organising properties of the system and in environmental inputs (biochemical and biophysical). The classical view of morphogenesis is based on the study of external biochemical molecules, such as morphogens. However, recent studies are establishing that the mechanical environment is also used by cells to communicate within tissues, suggesting that this mechanical crosstalk is essential to synchronise morphogenetic transitions and self-organisation. In this article we discuss how tissue interaction drive robust morphogenesis, starting from a classical biochemical view, to finalise with more recent advances on how the biophysical properties of a tissue feedback with their surroundings to allow form acquisition. We also comment on how in silico models aid to integrate and predict changes in cell and tissue behaviour. Finally, considering recent advances from the developmental biomechanics field showing that mechanical inputs work as cues that promote morphogenesis, we invite to revisit the concept of morphogen.

The logistics of Wnt production and delivery

Wnts are secreted proteins that control stem cell maintenance, cell fate decisions, and growth during development and adult homeostasis. Wnts carry a post-translational modification not seen in any other secreted protein: during biosynthesis, they are appended with a palmitoleoyl moiety that is required for signaling but also impairs solubility and hence diffusion in the extracellular space. In some contexts, Wnts act only in a juxtacrine manner but there are also instances of long range action. Several proteins and processes ensure that active Wnts reach the appropriate target cells. Some, like Porcupine, Wntless, and Notum are dedicated to Wnt function; we describe their activities in molecular detail. We also outline how the cell infrastructure (secretory, endocytic, and retromer pathways) contribute to the progression of Wnts from production to delivery. We then address how Wnts spread in the extracellular space and form a signaling gradient despite carrying a hydrophobic moiety. We highlight particularly the role of lipid-binding Wnt interactors and heparan sulfate proteoglycans. Finally, we briefly discuss how evolution might have led to the emergence of this unusual signaling pathway.

The present and future of Turing models in developmental biology

The Turing model (or reaction-diffusion model), first published in 1952, is a mathematical model that can account for autonomy in the morphogenesis of organisms. Although initially controversial, the model has gradually gained wider acceptance among experimental embryologists due to the accumulation of experimental data to support it. More recently, this model and others based on it have been used not only to explain biological phenomena conceptually but also as working hypotheses for molecular-level experiments and as internal components of more-complex 3D models. In this Spotlight, I will provide a personal perspective from an experimental biologist on some of the recent developments of the Turing model.

Synthetic gene circuits for cell state detection and protein tuning in human pluripotent stem cells

During development, cell state transitions are coordinated through changes in the identity of molecular regulators in a cell type‐ and dose‐specific manner. The ability to rationally engineer such transitions in human pluripotent stem cells (hPSC) will enable numerous applications in regenerative medicine. Herein, we report the generation of synthetic gene circuits that can detect a desired cell state using AND‐like logic integration of endogenous miRNAs (classifiers) and, upon detection, produce fine‐tuned levels of output proteins using an miRNA‐mediated output fine‐tuning technology (miSFITs). Specifically, we created an “hPSC ON” circuit using a model‐guided miRNA selection and circuit optimization approach. The circuit demonstrates robust PSC‐specific detection and graded output protein production. Next, we used an empirical approach to create an “hPSC‐Off” circuit. This circuit was applied to regulate the secretion of endogenous BMP4 in a state‐specific and fine‐tuned manner to control the composition of differentiating hPSCs. Our work provides a platform for customized cell state‐specific control of desired physiological factors in hPSC, laying the foundation for programming cell compositions in hPSC‐derived tissues and beyond.

Regulatory factor identification for nodal genes in zebrafish by causal inference

Activation of nodal genes is critical for mesoderm and endoderm induction. Our previous study reported that zebrafish nodal genes ndr1/squint and ndr2/cyclops are coordinately regulated by maternal Eomesa, Hwa-activated β-catenin (Hwa/β-catenin) signaling, and Nodal autoregulation (Nodal/Smad2) signaling. However, the exact contribution and underlying mechanisms are still elusive. Here, we applied “causal inference” to evaluate the causal between the independent and dependent variables, and we found that Hwa/β-catenin and Smad2 are the cause of ndr1 activation, while Eomesa is the cause of ndr2 activation. Mechanistically, the different cis-regulatory regions of ndr1 and ndr2 bound by Eomesa, β-catenin, and Smad2 were screened out via ChIP-qPCR and verified by the transgene constructs. The marginal GFP expression driven by ndr1 transgenesis could be diminished without both maternal Eomesa and Hwa/β-catenin, while Eomesa, not β-catenin, could bind and activate ndr2 demonstrated by ndr2 transgenesis. Thus, the distinct regulation of ndr1/ndr2 relies on different cis-regulatory regions.

Single-molecule tracking of Nodal and Lefty in live zebrafish embryos supports hindered diffusion model

The hindered diffusion model postulates that the movement of a signaling molecule through an embryo is affected by tissue geometry and binding-mediated hindrance, but these effects have not been directly demonstrated in vivo. Here, we visualize extracellular movement and binding of individual molecules of the activator-inhibitor signaling pair Nodal and Lefty in live developing zebrafish embryos using reflected light-sheet microscopy. We observe that diffusion coefficients of molecules are high in extracellular cavities, whereas mobility is reduced and bound fractions are high within cell-cell interfaces. Counterintuitively, molecules nevertheless accumulate in cavities, which we attribute to the geometry of the extracellular space by agent-based simulations. We further find that Nodal has a larger bound fraction than Lefty and shows a binding time of tens of seconds. Together, our measurements and simulations provide direct support for the hindered diffusion model and yield insights into the nanometer-to-micrometer-scale mechanisms that lead to macroscopic signal dispersal.

Static and Dynamic Gradient Based Directional Transportation of Neutral Molecules in Swollen Polymer Films

Materials which selectively transport molecules offer powerful opportunities for concentrating and separating chemical agents. Here, utilizing static and dynamic chemical gradients, transport of molecules within swollen crosslinked polymers is demonstrated. Using an ≈200 μm static hydroxyl to hexyl gradient, the neutral ambipolar nerve agent surrogate diethyl (cyanomethyl)phosphonate (DECP) is directionally transported and concentrated 60-fold within 4 hours. To accelerate transport kinetics, a dynamic gradient (a "travelling wave") is utilized. Here, the non-polar dye pyrene was transported. The dynamic gradient is generated by an ion exchange process triggered by the localized introduction of an aqueous NaCl solution, which converts the gel from hydrophobic to hydrophilic. As the hydrophilic region expands, associated water enters the gel, and pyrene is pushed ahead of the expansion front. The dynamic gradient provides about 10-fold faster transport kinetics than the static gradient.

Patterning principles of morphogen gradients

Metazoan embryos develop from a single cell into three-dimensional structured organisms while groups of genetically identical cells attain specialized identities. Cells of the developing embryo both create and accurately interpret morphogen gradients to determine their positions and make specific decisions in response. Here, we first cover intellectual roots of morphogen and positional information concepts. Focusing on animal embryos, we then provide a review of current understanding on how morphogen gradients are established and how their spans are controlled. Lastly, we cover how gradients evolve in time and space during development, and how they encode information to control patterning. In sum, we provide a list of patterning principles for morphogen gradients and review recent advances in quantitative methodologies elucidating information provided by morphogens.

A self-generated Toddler gradient guides mesodermal cell migration

The sculpting of germ layers during gastrulation relies on the coordinated migration of progenitor cells, yet the cues controlling these long-range directed movements remain largely unknown. While directional migration often relies on a chemokine gradient generated from a localized source, we find that zebrafish ventrolateral mesoderm is guided by a self-generated gradient of the initially uniformly expressed and secreted protein Toddler/ELABELA/Apela. We show that the Apelin receptor, which is specifically expressed in mesodermal cells, has a dual role during gastrulation, acting as a scavenger receptor to generate a Toddler gradient, and as a chemokine receptor to sense this guidance cue. Thus, we uncover a single receptor-based self-generated gradient as the enigmatic guidance cue that can robustly steer the directional migration of mesoderm through the complex and continuously changing environment of the gastrulating embryo.

A conserved molecular template underlies color pattern diversity in estrildid finches

The color patterns that adorn animals' coats not only exhibit extensive diversity linked to various ecological functions but also display recurrences in geometry, orientation, or body location. How processes of pattern formation shape such phenotypic trends remains a mystery. Here, we surveyed plumage color patterns in passerine finches displaying extreme apparent variation and identified a conserved set of color domains. We linked these domains to putative embryonic skin regions instructed by early developmental tissues and outlined by the combinatory expression of few genetic markers. We found that this embryonic prepattern is largely conserved in birds displaying drastic color differences in the adult, interspecies variation resulting from the masking or display of each domain depending on their coloration. This work showed that a simple molecular landscape serves as common spatial template to extensive color pattern variation in finches, revealing that early conserved landmarks and molecular pathways are a major cause of phenotypic trends.

Long-range morphogen gradient formation by cell-to-cell signal propagation

Morphogen gradients are a central concept in developmental biology. Their formation often involves the secretion of morphogens from a local source, that spread by diffusion in the cell field, where molecules eventually get degraded. This implies limits to both the time and length scales over which morphogen gradients can form which are set by diffusion coefficients and degradation rates. Towards the goal of identifying plausible mechanisms capable of extending the gradient range, we here use theory to explore properties of a cell-to-cell signaling relay. Inspired by the millimeter-scale Wnt-expression and signaling gradients in flatworms, we consider morphogen-mediated morphogen production in the cell field. We show that such a relay can generate stable morphogen and signaling gradients that are oriented by a local, morphogen-independent source of morphogen at a boundary. This gradient formation can be related to an effective diffusion and an effective degradation that result from morphogen production due to signaling relay. If the secretion of morphogen produced in response to the relay is polarized, it further gives rise to an effective drift. We find that signaling relay can generate long-ranged gradients in relevant times without relying on extreme choices of diffusion coefficients or degradation rates, thus exceeding the limits set by physiological diffusion coefficients and degradation rates. A signaling relay is hence an attractive principle to conceptualize long-range gradient formation by slowly diffusing morphogens that are relevant for patterning in adult contexts such as regeneration and tissue turn-over.

Stripe and spot selection in cusp patterning of mammalian molar formation

Tooth development is governed largely by epithelial-mesenchymal interactions and is mediated by numerous signaling pathways. This type of morphogenetic processes has been explained by reaction-diffusion systems, especially in the framework of a Turing model. Here we focus on morphological and developmental differences between upper and lower molars in mice by modeling 2D pattern formation in a Turing system. Stripe vs. spot patterns are the primary types of variation in a Turing model. We show that the complexity of the cusp cross-sections can distinguish between stripe vs. spot patterns, and mice have stripe-like upper and spot-like lower molar morphologies. Additionally, our computational modeling that incorporates empirical data on tooth germ growth traces the order of cusp formation and relative position of the cusps in upper and lower molars in mice. We further propose a hypothetical framework of developmental mechanism that could help us understand the evolution of the highly variable nature of mammalian molars associated with the acquisition of the hypocone and the increase of lophedness.

Maternal Factors and Nodal Autoregulation Orchestrate Nodal Gene Expression for Embryonic Mesendoderm Induction in the Zebrafish

Nodal proteins provide crucial signals for mesoderm and endoderm induction. In zebrafish embryos, the nodal genes ndr1/squint and ndr2/cyclops are implicated in mesendoderm induction. It remains elusive how ndr1 and ndr2 expression is regulated spatiotemporally. Here we investigated regulation of ndr1 and ndr2 expression using Mhwa mutants that lack the maternal dorsal determinant Hwa with deficiency in β-catenin signaling, Meomesa mutants that lack maternal Eomesodermin A (Eomesa), Meomesa;Mhwa double mutants, and the Nodal signaling inhibitor SB431542. We show that ndr1 and ndr2 expression is completely abolished in Meomesa;Mhwa mutant embryos, indicating an essential role of maternal eomesa and hwa. Hwa-activated β-catenin signaling plays a major role in activation of ndr1 expression in the dorsal blastodermal margin, while eomesa is mostly responsible for ndr1 expression in the lateroventral margin and Nodal signaling contributes to ventral expansion of the ndr1 expression domain. However, ndr2 expression mainly depends on maternal eomesa with minor or negligible contribution of maternal hwa and Nodal autoregulation. These mechanisms may help understand regulation of Nodal expression in other species.

The role of cell geometry and cell-cell communication in gradient sensing

Cells can measure shallow gradients of external signals to initiate and accomplish a migration or a morphogenetic process. Recently, starting from mathematical models like the local-excitation global-inhibition (LEGI) model and with the support of empirical evidence, it has been proposed that cellular communication improves the measurement of an external gradient. However, the mathematical models that have been used have over-simplified geometries (e.g., they are uni-dimensional) or assumptions about cellular communication, which limit the possibility to analyze the gradient sensing ability of more complex cellular systems. Here, we generalize the existing models to study the effects on gradient sensing of cell number, geometry and of long- versus short-range cellular communication in 2D systems representing epithelial tissues. We find that increasing the cell number can be detrimental for gradient sensing when the communication is weak and limited to nearest neighbour cells, while it is beneficial when there is long-range communication. We also find that, with long-range communication, the gradient sensing ability improves for tissues with more disordered geometries; on the other hand, an ordered structure with mostly hexagonal cells is advantageous with nearest neighbour communication. Our results considerably extend the current models of gradient sensing by epithelial tissues, making a step further toward predicting the mechanism of communication and its putative mediator in many biological processes.

Precision of morphogen gradients in neural tube development

Morphogen gradients encode positional information during development. How high patterning precision is achieved despite natural variation in both the morphogen gradients and in the readout process, is still largely elusive. Here, we show that the positional error of gradients in the mouse neural tube has previously been overestimated, and that the reported accuracy of the central progenitor domain boundaries in the mouse neural tube can be achieved with a single gradient, rather than requiring the simultaneous readout of opposing gradients. Consistently and independently, numerical simulations based on measured molecular noise levels likewise result in lower gradient variabilities than reported. Finally, we show that the patterning mechanism yields progenitor cell numbers with even greater precision than boundary positions, as gradient amplitude changes do not affect interior progenitor domain sizes. We conclude that single gradients can yield the observed developmental precision, which provides prospects for tissue engineering.

Establishment and interpretation of NODAL and BMP signaling gradients in early vertebrate development

Transforming growth factor β (TGF-β) family ligands play crucial roles in orchestrating early embryonic development. Most significantly, two family members, NODAL and BMP form signaling gradients and indeed in fish, frogs and sea urchins these two opposing gradients are sufficient to organize a complete embryonic axis. This review focuses on how these gradients are established and interpreted during early vertebrate development. The review highlights key principles that are emerging, in particular the importance of signaling duration as well as ligand concentration in both gradient generation and their interpretation. Feedforward and feedback loops involving other signaling pathways are also essential for providing spatial and temporal information downstream of the NODAL and BMP signaling pathways. Finally, new data suggest the existence of buffering mechanisms, whereby early signaling defects can be readily corrected downstream later in development, suggesting that signaling gradients do not have to be as precise as previously thought.

Waves in Embryonic Development

Embryonic development hinges on effective coordination of molecular events across space and time. Waves have recently emerged as constituting an ubiquitous mechanism that ensures rapid spreading of regulatory signals across embryos, as well as reliable control of their patterning, namely, for the emergence of body plan structures. In this article, we review a selection of recent quantitative work on signaling waves and present an overview of the theory of waves. Our aim is to provide a succinct yet comprehensive guiding reference for the theoretical frameworks by which signaling waves can arise in embryos. We start, then, from reaction-diffusion systems, both static and time dependent; move to excitable dynamics; and conclude with systems of coupled oscillators. We link these theoretical models to molecular mechanisms recently elucidated for the control of mitotic waves in early embryos, patterning of the vertebrate body axis, micropattern cultures, and bone regeneration. Our goal is to inspire experimental work that will advance theory in development and connect its predictions to quantitative biological observations. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Nodal is a short-range morphogen with activity that spreads through a relay mechanism in human gastruloids

Morphogens are signaling molecules that convey positional information and dictate cell fates during development. Although ectopic expression in model organisms suggests that morphogen gradients form through diffusion, little is known about how morphogen gradients are created and interpreted during mammalian embryogenesis due to the combined difficulties of measuring endogenous morphogen levels and observing development in utero. Here we take advantage of a human gastruloid model to visualize endogenous Nodal protein in living cells, during specification of germ layers. We show that Nodal is extremely short range so that Nodal protein is limited to the immediate neighborhood of source cells. Nodal activity spreads through a relay mechanism in which Nodal production induces neighboring cells to transcribe Nodal. We further show that the Nodal inhibitor Lefty, while biochemically capable of long-range diffusion, also acts locally to control the timing of Nodal spread and therefore of mesoderm differentiation during patterning. Our study establishes a paradigm for tissue patterning by an activator-inhibitor pair.

Size matters: tissue size as a marker for a transition between reaction-diffusion regimes in spatio-temporal distribution of morphogens

The reaction-diffusion model constitutes one of the most influential mathematical models to study distribution of morphogens in tissues. Despite its widespread use, the effect of finite tissue size on model-predicted spatio-temporal morphogen distributions has not been completely elucidated. In this study, we analytically investigated the spatio-temporal distributions of morphogens predicted by a reaction-diffusion model in a finite one-dimensional domain, as a proxy for a biological tissue, and compared it with the solution of the infinite-domain model. We explored the reduced parameter, the tissue length in units of a characteristic reaction-diffusion length, and identified two reaction-diffusion regimes separated by a crossover tissue size estimated in approximately three characteristic reaction-diffusion lengths. While above this crossover the infinite-domain model constitutes a good approximation, it breaks below this crossover, whereas the finite-domain model faithfully describes the entire parameter space. We evaluated whether the infinite-domain model renders accurate estimations of diffusion coefficients when fitted to finite spatial profiles, a procedure typically followed in fluorescence recovery after photobleaching (FRAP) experiments. We found that the infinite-domain model overestimates diffusion coefficients when the domain is smaller than the crossover tissue size. Thus, the crossover tissue size may be instrumental in selecting the suitable reaction-diffusion model to study tissue morphogenesis.

Front speed and pattern selection of a propagating chemical front in an active fluid

Spontaneous pattern formation in living systems is driven by reaction-diffusion chemistry and active mechanics. The feedback between chemical and mechanical forces is often essential to robust pattern formation, yet it remains poorly understood in general. In this analytical and numerical paper, we study an experimentally motivated minimal model of coupling between reaction-diffusion and active matter: a propagating front of an autocatalytic and stress-generating species. In the absence of activity, the front is described by the well-studied Kolmogorov, Petrovsky, and Piskunov equation. We find that front propagation is maintained even in active systems, with crucial differences: an extensile stress increases the front speed beyond a critical magnitude of the stress, while a contractile stress has no effect on the front speed but can generate a periodic instability in the high-concentration region behind the front. We expect our results to be useful in interpreting pattern formation in active systems with mechanochemical coupling in vivo and in vitro.

Turing pattern design principles and their robustness

Turing patterns have morphed from mathematical curiosities into highly desirable targets for synthetic biology. For a long time, their biological significance was sometimes disputed but there is now ample evidence for their involvement in processes ranging from skin pigmentation to digit and limb formation. While their role in developmental biology is now firmly established, their synthetic design has so far proved challenging. Here, we review recent large-scale mathematical analyses that have attempted to narrow down potential design principles. We consider different aspects of robustness of these models and outline why this perspective will be helpful in the search for synthetic Turing-patterning systems. We conclude by considering robustness in the context of developmental modelling more generally. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

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

Systems for intricate patterning of the vertebrate anatomy

Periodic patterns form intricate arrays in the vertebrate anatomy, notably the hair and feather follicles of the skin, but also internally the villi of the gut and the many branches of the lung, kidney, mammary and salivary glands. These tissues are composite structures, being composed of adjoined epithelium and mesenchyme, and the patterns that arise within them require interaction between these two tissue layers. In embryonic development, cells change both their distribution and state in a periodic manner, defining the size and relative positions of these specialized structures. Their placement is determined by simple spacing mechanisms, with substantial evidence pointing to a variety of local enhancement/lateral inhibition systems underlying the breaking of symmetry. The nature of the cellular processes involved, however, has been less clear. While much attention has focused on intercellular soluble signals, such as protein growth factors, experimental evidence has grown for contributions of cell movement or mechanical forces to symmetry breaking. In the mesenchyme, unlike the epithelium, cells may move freely and can self-organize into aggregates by chemotaxis, or through generation and response to mechanical strain on their surrounding matrix. Different modes of self-organization may coexist, either coordinated into a single system or with hierarchical relationships. Consideration of a broad range of distinct biological processes is required to advance understanding of biological pattern formation. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Computational modelling of nephron progenitor cell movement and aggregation during kidney organogenesis

During early kidney organogenesis, nephron progenitor (NP) cells move from the tip to the corner region of the ureteric bud (UB) branches in order to form the pretubular aggregate, the early structure giving rise to nephron formation. NP cells derive from metanephric mesenchymal cells and physically interact with them during the movement. Chemotaxis and cell-cell adhesion differences are believed to drive the cell patterning during this critical period of organogenesis. However, the effect of these forces to the cell patterns and their respective movements are known in limited details. We applied a Cellular Potts Model to explore how these forces and organizations contribute to directed cell movement and aggregation. Model parameters were estimated based on fitting to experimental data obtained in ex vivo kidney explant and dissociation-reaggregation organoid culture studies. Our simulations indicated that optimal enrichment and aggregation of NP cells in the UB corner niche requires chemoattractant secretion from both the UB epithelial cells and the NP cells themselves, as well as differences in cell-cell adhesion energies. Furthermore, NP cells were observed, both experimentally and by modelling, to move at higher speed in the UB corner as compared to the tip region where they originated. The existence of different cell speed domains along the UB was confirmed using self-organizing map analysis. In summary, we saw faster NP cell movements near aggregation. The applicability of Cellular Potts Model approach to simulate cell movement and patterning was found to be good during for this early nephrogenesis process. Further refinement of the model should allow us to recapitulate the effects of developmental changes of cell phenotypes and molecular crosstalk during further organ development.