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

Improving Parameter Inference from FRAP Data: an Analysis Motivated by Pattern Formation in the Drosophila Wing Disc

Fluorescence recovery after photobleaching (FRAP) is used to obtain quantitative information about molecular diffusion and binding kinetics at both cell and tissue levels of organization. FRAP models have been proposed to estimate the diffusion coefficients and binding kinetic parameters of species for a variety of biological systems and experimental settings. However, it is not clear what the connection among the diverse parameter estimates from different models of the same system is, whether the assumptions made in the model are appropriate, and what the qualities of the estimates are. Here we propose a new approach to investigate the discrepancies between parameters estimated from different models. We use a theoretical model to simulate the dynamics of a FRAP experiment and generate the data that are used in various recovery models to estimate the corresponding parameters. By postulating a recovery model identical to the theoretical model, we first establish that the appropriate choice of observation time can significantly improve the quality of estimates, especially when the diffusion and binding kinetics are not well balanced, in a sense made precise later. Secondly, we find that changing the balance between diffusion and binding kinetics by changing the size of the bleaching region, which gives rise to different FRAP curves, provides a priori knowledge of diffusion and binding kinetics, which is important for model formulation. We also show that the use of the spatial information in FRAP provides better parameter estimation. By varying the recovery model from a fixed theoretical model, we show that a simplified recovery model can adequately describe the FRAP process in some circumstances and establish the relationship between parameters in the theoretical model and those in the recovery model. We then analyze an example in which the data are generated with a model of intermediate complexity and the parameters are estimated using models of greater or less complexity, and show how sensitivity analysis can be used to improve FRAP model formulation. Lastly, we show how sophisticated global sensitivity analysis can be used to detect over-fitting when using a model that is too complex.

Morphogene adsorption as a Turing instability regulator: Theoretical analysis and possible applications in multicellular embryonic systems

The Turing instability in the reaction-diffusion system is a widely recognized mechanism of the morphogen gradient self-organization during the embryonic development. One of the essential conditions for such self-organization is sharp difference in the diffusion rates of the reacting substances (morphogens). In classical models this condition is satisfied only for significantly different values of diffusion coefficients which cannot hold for morphogens of similar molecular size. One of the most realistic explanations of the difference in diffusion rate is the difference between adsorption of morphogens to the extracellular matrix (ECM). Basing on this assumption we develop a novel mathematical model and demonstrate its effectiveness in describing several well-known examples of biological patterning. Our model consisting of three reaction-diffusion equations has the Turing-type instability and includes two elements with equal diffusivity and immobile binding sites as the third reaction substance. The model is an extension of the classical Gierer-Meinhardt two-components model and can be reduced to it under certain conditions. Incorporation of ECM in the model system allows us to validate the model for available experimental parameters. According to our model introduction of binding sites gradient, which is frequently observed in embryonic tissues, allows one to generate more types of different spatial patterns than can be obtained with two-components models. Thus, besides providing an essential condition for the Turing instability for the system of morphogen with close values of the diffusion coefficients, the morphogen adsorption on ECM may be important as a factor that increases the variability of self-organizing structures.

Weak synchronization and large-scale collective oscillation in dense bacterial suspensions

Collective oscillatory behaviour is ubiquitous in nature, having a vital role in many biological processes from embryogenesis and organ development to pace-making in neuron networks. Elucidating the mechanisms that give rise to synchronization is essential to the understanding of biological self-organization. Collective oscillations in biological multicellular systems often arise from long-range coupling mediated by diffusive chemicals, by electrochemical mechanisms, or by biomechanical interaction between cells and their physical environment. In these examples, the phase of some oscillatory intracellular degree of freedom is synchronized. Here, in contrast, we report the discovery of a weak synchronization mechanism that does not require long-range coupling or inherent oscillation of individual cells. We find that millions of motile cells in dense bacterial suspensions can self-organize into highly robust collective oscillatory motion, while individual cells move in an erratic manner, without obvious periodic motion but with frequent, abrupt and random directional changes. So erratic are individual trajectories that uncovering the collective oscillations of our micrometre-sized cells requires individual velocities to be averaged over tens or hundreds of micrometres. On such large scales, the oscillations appear to be in phase and the mean position of cells typically describes a regular elliptic trajectory. We found that the phase of the oscillations is organized into a centimetre-scale travelling wave. We present a model of noisy self-propelled particles with strictly local interactions that accounts faithfully for our observations, suggesting that self-organized collective oscillatory motion results from spontaneous chiral and rotational symmetry breaking. These findings reveal a previously unseen type of long-range order in active matter systems (those in which energy is spent locally to produce non-random motion). This mechanism of collective oscillation may inspire new strategies to control the self-organization of active matter and swarming robots.

Establishment of the Vertebrate Germ Layers

The process of germ layer formation is a universal feature of animal development. The germ layers separate the cells that produce the internal organs and tissues from those that produce the nervous system and outer tissues. Their discovery in the early nineteenth century transformed embryology from a purely descriptive field into a rigorous scientific discipline, in which hypotheses could be tested by observation and experimentation. By systematically addressing the questions of how the germ layers are formed and how they generate overall body plan, scientists have made fundamental contributions to the fields of evolution, cell signaling, morphogenesis, and stem cell biology. At each step, this work was advanced by the development of innovative methods of observing cell behavior in vivo and in culture. Here, we take an historical approach to describe our current understanding of vertebrate germ layer formation as it relates to the long-standing questions of developmental biology. By comparing how germ layers form in distantly related vertebrate species, we find that highly conserved molecular pathways can be adapted to perform the same function in dramatically different embryonic environments.

Controlling contractile instabilities in the actomyosin cortex

The actomyosin cell cortex is an active contractile material for driving cell- and tissue morphogenesis. The cortex has a tendency to form a pattern of myosin foci, which is a signature of potentially unstable behavior. How a system that is prone to such instabilities can rveliably drive morphogenesis remains an outstanding question. Here, we report that in the Caenorhabditis elegans zygote, feedback between active RhoA and myosin induces a contractile instability in the cortex. We discover that an independent RhoA pacemaking oscillator controls this instability, generating a pulsatory pattern of myosin foci and preventing the collapse of cortical material into a few dynamic contracting regions. Our work reveals how contractile instabilities that are natural to occur in mechanically active media can be biochemically controlled to robustly drive morphogenetic events.

Plasticity underlies tumor progression: role of Nodal signaling

The transforming growth factor beta (TGFβ) superfamily member Nodal is an established regulator of early embryonic development, with primary roles in endoderm induction, left-right asymmetry, and primitive streak formation. Nodal signals through TGFβ family receptors at the plasma membrane and induces signaling cascades leading to diverse transcriptional regulation. While conceptually simple, the regulation of Nodal and its molecular effects are profoundly complex and context dependent. Pioneering work by developmental biologists has characterized the signaling pathways, regulatory components, and provided detailed insight into the mechanisms by which Nodal mediates changes at the cellular and organismal levels. Nodal is also an important factor in maintaining pluripotency of embryonic stem cells through regulation of core transcriptional programs. Collectively, this work has led to an appreciation for Nodal as a powerful morphogen capable of orchestrating multiple cellular phenotypes. Although Nodal is not active in most adult tissues, its reexpression and signaling have been linked to multiple types of human cancer, and Nodal has emerged as a driver of tumor growth and cellular plasticity. In vitro and in vivo experimental evidence has demonstrated that inhibition of Nodal signaling reduces cancer cell aggressive characteristics, while clinical data have established associations with Nodal expression and patient outcomes. As a result, there is great interest in the potential targeting of Nodal activity in a therapeutic setting for cancer patients that may provide new avenues for suppressing tumor growth and metastasis. In this review, we evaluate our current understanding of the complexities of Nodal function in cancer and highlight recent experimental evidence that sheds light on the therapeutic potential of its inhibition.

Perspectives on Intra- and Intercellular Trafficking of Hedgehog for Tissue Patterning

Intercellular communication is a fundamental process for correct tissue development. The mechanism of this process involves, among other things, the production and secretion of signaling molecules by specialized cell types and the capability of these signals to reach the target cells in order to trigger specific responses. Hedgehog (Hh) is one of the best-studied signaling pathways because of its importance during morphogenesis in many organisms. The Hh protein acts as a morphogen, activating its targets at a distance in a concentration-dependent manner. Post-translational modifications of Hh lead to a molecule covalently bond to two lipid moieties. These lipid modifications confer Hh high affinity to lipidic membranes, and intense studies have been carried out to explain its release into the extracellular matrix. This work reviews Hh molecule maturation, the intracellular recycling needed for its secretion and the proposed carriers to explain Hh transportation to the receiving cells. Special focus is placed on the role of specialized filopodia, also named cytonemes, in morphogen transport and gradient formation.

Modelling Chemotactic Motion of Cells in Biological Tissues

Developmental processes in biology are underlined by proliferation, differentiation and migration of cells. The latter two are interlinked since cellular differentiation is governed by the dynamics of morphogens which, in turn, is affected by the movement of cells. Mutual effects of morphogenetic and cell movement patterns are enhanced when the movement is due to chemotactic response of cells to the morphogens. In this study we introduce a mathematical model to analyse how this interplay can result in a steady movement of cells in a tissue and associated formation of travelling waves in a concentration field of morphogen. Using the model we have identified four chemotactic scenarios for migration of single cell or homogeneous group of cells in a tissue. Such a migration can take place if moving cells are (1) repelled by a chemical produced by themselves or (2) attracted by a chemical produced by the surrounding cells in a tissue. Furthermore, the group of cells can also move if cells in surrounding tissue are (3) repelled by a chemical produced by moving cells or (4) attracted by a chemical produced by surrounding cells themselves. The proposed mechanisms can underlie migration of cells during embryonic development as well as spread of metastatic cells.

Self-organization in development, regeneration and organoids

Self-organization of cells is a fundamental design principle in biology, yet the inherent non-linearity of self-organizing systems often poses significant challenges in deciphering the underlying mechanisms. Here, we discuss recent progress in this respect, focusing on examples from development, regeneration and organoid differentiation. Together, these three paradigms emphasize the active material properties of tissues that result from the functional coupling between individual cells as active units. Further, we discuss the challenge of obtaining reproducible outcomes on the basis of self-organizing systems, which development and regeneration, but not the current organoid culture protocols, achieve.

Quantitative imaging reveals real-time Pou5f3-Nanog complexes driving dorsoventral mesendoderm patterning in zebrafish

Formation of the three embryonic germ layers is a fundamental developmental process that initiates differentiation. How the zebrafish pluripotency factor Pou5f3 (homologous to mammalian Oct4) drives lineage commitment is unclear. Here, we introduce fluorescence lifetime imaging microscopy and fluorescence correlation spectroscopy to assess the formation of Pou5f3 complexes with other transcription factors in real-time in gastrulating zebrafish embryos. We show, at single-cell resolution in vivo, that Pou5f3 complexes with Nanog to pattern mesendoderm differentiation at the blastula stage. Later, during gastrulation, Sox32 restricts Pou5f3–Nanog complexes to the ventrolateral mesendoderm by binding Pou5f3 or Nanog in prospective dorsal endoderm. In the ventrolateral endoderm, the Elabela / Aplnr pathway limits Sox32 levels, allowing the formation of Pou5f3–Nanog complexes and the activation of downstream BMP signaling. This quantitative model shows that a balance in the spatiotemporal distribution of Pou5f3–Nanog complexes, modulated by Sox32, regulates mesendoderm specification along the dorsoventral axis.

DOI:http://dx.doi.org/10.7554/eLife.11475.001

As an animal embryo develops, cells divide and establish three distinct layers called the ectoderm, mesoderm and endoderm. Proteins called transcription factors control this process by regulating the activity of particular genes. Two or more transcription factors may interact to modulate each other’s activity. Zebrafish embryos provide an ideal model system for monitoring how these embryonic layers form and the interactions between transcription factors in real-time because they are transparent and develop outside their parents. Pou5f3 and Nanog are two key transcription factors involved in this process in zebrafish. However, it is not clear how Pou5f3 and Nanog instruct cells to become ectoderm, mesoderm or endoderm.

Perez Camps et al. used imaging techniques to study Pou5f3 and Nanog. The experiments show that Pou5f3 and Nanog bind together to form complexes that instruct cells to form the temporary layer that later gives rise to both the mesoderm and endoderm. The cells in which there are less Pou5f3 and Nanog complexes form the ectoderm layer.

To develop the body shape of adult zebrafish, the embryos need to give individual cells information about their location in the body. For example, a signal protein called bone morphogenetic protein (BMP) accumulates on the side of the embryo that will become the underside of the fish. Perez Camps et al. show that once the endoderm, mesoderm and ectoderm have formed, Pou5f3–Nanog complexes regulate BMP signalling to specify the underside of the fish. Meanwhile, in the endoderm on the opposite side, another transcription factor called Sox32 binds to individual Pou5f3 and Nanog proteins. This prevents Pou5f3 and Nanog from forming complexes and determines which side of the embryo will make the topside of the fish. A future challenge is to explore other transcription factors that may prevent Pou5f1 and Nanog from binding in the mesoderm and ectoderm of the topside of the fish.

DOI:http://dx.doi.org/10.7554/eLife.11475.002

Lefty Glycoproteins in Human Embryonic Stem Cells: Extracellular Delivery Route and Posttranslational Modification in Differentiation

Lefty is a member of transforming growth factor-beta (TGF-β) superfamily and a potent antagonist of the TGF-β/Nodal/Activin signaling pathway. Lefty is critical in sustaining self-renewal/pluripotency status, and implicated in the differentiation of embryonic stem cells (ESCs). However, emerging studies depict Lefty as a multifaceted protein involved in myriad cellular events. Lefty proteins (human Lefty A and B) are secreted glycoproteins, but their mode of secretion and the significance of their "glycan" moiety remain mostly unexplored. By employing an in vitro system of human ESCs (hESCs), we observed that Lefty protein(s) are encased in exosomes for extracellular release. The exosomal- and cell-associated Lefty diverge in their proteolytic processing, and possess N-glycan structures of high mannose and complex nature. Differentiation of hESCs to mesenchymal cells (MSCs) or neuronal progenitor cells (NPCs) entails distinct changes in the Lefty A/Lefty B gene(s), and protein expression. Specifically, the proteolytic cleavage and N-glycan composition of the cell-associated and exosomal Lefty differ in the differentiated progenies. These modifications affected Lefty's inhibitory effect on Nodal signaling in aggressive melanoma cells. The microheterogeneity in the processing and glycosylation of Lefty protein(s) between hESCs, MSCs, and NPCs could present efficient means of diversifying the endogenous functions of Lefty. Whether Lefty's diverse functions in embryonic patterning, as well as its diffusion range in the extracellular environment, are similarly affected remains to be determined. Our studies underscore the potential relevance of Lefty-packaged exosomes for combating debilitating diseases such as cancer.

Pluripotent stem cells as a model for embryonic patterning: From signaling dynamics to spatial organization in a dish

In vivo studies have identified the signaling pathways and transcription factors involved in patterning the vertebrate embryo, but much remains unknown about how these are organized in space and time to orchestrate embryogenesis. Recently, embryonic stem cells have been established as a platform for studying spatial pattern formation and differentiation dynamics in the early mammalian embryo. The ease of observing and manipulating stem cell systems promises to fill gaps in our understanding of developmental dynamics and identify aspects that are uniquely human. Developmental Dynamics 245:976-990, 2016. © 2016 Wiley Periodicals, Inc.

Agonists and Antagonists of TGF-β Family Ligands

The discovery of the transforming growth factor β (TGF-β) family ligands and the realization that their bioactivities need to be tightly controlled temporally and spatially led to intensive research that has identified a multitude of extracellular modulators of TGF-β family ligands, uncovered their functions in developmental and pathophysiological processes, defined the mechanisms of their activities, and explored potential modulator-based therapeutic applications in treating human diseases. These studies revealed a diverse repertoire of extracellular and membrane-associated molecules that are capable of modulating TGF-β family signals via control of ligand availability, processing, ligand-receptor interaction, and receptor activation. These molecules include not only soluble ligand-binding proteins that were conventionally considered as agonists and antagonists of TGF-β family of growth factors, but also extracellular matrix (ECM) proteins and proteoglycans that can serve as "sink" and control storage and release of both the TGF-β family ligands and their regulators. This extensive network of soluble and ECM modulators helps to ensure dynamic and cell-specific control of TGF-β family signals. This article reviews our knowledge of extracellular modulation of TGF-β growth factors by diverse proteins and their molecular mechanisms to regulate TGF-β family signaling.

Precise Editing of the Zebrafish Genome Made Simple and Efficient

We present simple and efficient methods for creating heritable modifications of the zebrafish genome. Precisely modified alleles are generated by homologous recombination between the host genome and dsDNA donor molecules, stimulated by the induction of chromosomally targeted double-strand breaks. Several kilobase-long tracts of genome sequence can be replaced. Tagging donor sequences with reporter genes that can be subsequently excised improves recovery of edited alleles by an order of magnitude and facilitates recovery of recessive and phenotypically silent conditional mutations. We generate and demonstrate functionality of (1) alleles with a single codon change, (2) an allele encoding an epitope-tagged version of an endogenous protein, (3) alleles expressing reporter proteins, and (4) a conditional allele in which an exon is flanked by recombinogenic loxP sites. Our methods make recovery of a broad range of genome editing events very practicable, significantly advancing applicability of the zebrafish for studying normal biological processes and modeling diseases.

Tissue culture on a chip: Developmental biology applications of self-organized capillary networks in microfluidic devices

Organ culture systems are used to elucidate the mechanisms of pattern formation in developmental biology. Various organ culture techniques have been used, but the lack of microcirculation in such cultures impedes the long-term maintenance of larger tissues. Recent advances in microfluidic devices now enable us to utilize self-organized perfusable capillary networks in organ cultures. In this review, we will overview past approaches to organ culture and current technical advances in microfluidic devices, and discuss possible applications of microfluidics towards the study of developmental biology.

Scaling of morphogenetic patterns in reaction-diffusion systems

Development of multicellular organisms is commonly associated with the response of individual cells to concentrations of chemical substances called morphogens. Concentration fields of morphogens form a basis for biological patterning and ensure its properties including ability to scale with the size of the organism. While mechanisms underlying the formation of morphogen gradients are reasonably well understood, little is known about processes responsible for their scaling. Here, we perform a formal analysis of scaling for chemical patterns forming in continuous systems. We introduce a quantity representing the sensitivity of systems to changes in their size and use it to analyse scaling properties of patterns forming in a few different systems. Particularly, we consider how scaling properties of morphogen gradients forming in diffusion-decay systems depend on boundary conditions and how the scaling can be improved by passive modulation of morphogens or active transport in the system. We also analyse scaling of morphogenetic signal caused by two opposing gradients and consider scaling properties of patterns forming in activator-inhibitor systems. We conclude with a few possible mechanisms which allow scaling of morphogenetic patterns.

Theoretical analysis of degradation mechanisms in the formation of morphogen gradients

Fundamental biological processes of development of tissues and organs in multicellular organisms are governed by various signaling molecules, which are called morphogens. It is known that spatial and temporal variations in the concentration profiles of signaling molecules, which are frequently referred as morphogen gradients, lead to a cell differentiation via activating specific genes in a concentration-dependent manner. It is widely accepted that the establishment of the morphogen gradients involves multiple biochemical reactions and diffusion processes. One of the critical elements in the formation of morphogen gradients is a degradation of signaling molecules. We develop a new theoretical approach that provides a comprehensive description of the degradation mechanisms. It is based on the idea that the degradation works as an effective potential that drives the signaling molecules away from the source region. Utilizing the method of first-passage processes, the dynamics of the formation of morphogen gradients for various degradation mechanisms is explicitly evaluated. It is found that linear degradation processes lead to a dynamic behavior specified by times to form the morphogen gradients that depend linearly on the distance from the source. This is because the effective potential due to the degradation is quite strong. At the same time, nonlinear degradation mechanisms yield a quadratic scaling in the morphogen gradients formation times since the effective potentials are much weaker. Physical-chemical explanations of these phenomena are presented.

Where to Go: Breaking the Symmetry in Cell Motility

Cell migration in the “correct” direction is pivotal for many biological processes. Although most work is devoted to its molecular mechanisms, the cell’s preference for one direction over others, thus overcoming intrinsic random motility, epitomizes a profound principle that underlies all complex systems: the choice of one axis, in structure or motion, from a uniform or symmetric set of options. Explaining directional motility by an external chemo-attractant gradient does not solve but only shifts the problem of causation: whence the gradient? A new study in PLOS Biology shows cell migration in a self-generated gradient, offering an opportunity to take a broader look at the old dualism of extrinsic instruction versus intrinsic symmetry-breaking in cell biology.

Directional cell motility is enabled by chemoattractant gradient and symmetry-breaking. This Primer argues that the recent observation of cells generating a gradient in a uniformly distributed nutrient reveals the multilayered nature of symmetry-breaking in cell locomotion.

Turing Patterning Using Gene Circuits with Gas-Induced Degradation of Quorum Sensing Molecules

The Turing instability was proposed more than six decades ago as a mechanism leading to spatial patterning, but it has yet to be exploited in a synthetic biology setting. Here we characterize the Turing instability in a specific gene circuit that can be implemented in vitro or in populations of clonal cells producing short-range activator N-Acyl homoserine lactone (AHL) and long-range inhibitor hydrogen peroxide (H2O2) gas. Slowing the production rate of the AHL-degrading enzyme, AiiA, generates stable fixed states, limit cycle oscillations and Turing patterns. Further tuning of signaling parameters determines local robustness and controls the range of unstable wavenumbers in the patterning regime. These findings provide a roadmap for optimizing spatial patterns of gene expression based on familiar quorum and gas sensitive E. coli promoters. The circuit design and predictions may be useful for (re)programming spatial dynamics in synthetic and natural gene expression systems.

Abstracting the principles of development using imaging and modeling

Here we look at modern developmental biology with a focus on the relationship between different approaches of investigation. We argue that direct imaging is a powerful approach not only for obtaining descriptive information but also for model generation and testing that lead to mechanistic insights. Modeling, on the other hand, conceptualizes imaging data and provides guidance to perturbations. The inquiry progresses most efficiently when a trinity of approaches—quantitative imaging (measurement), modeling (theory) and perturbation (test)—are pursued in concert, but not when one approach is dominant. Using recent studies of the zebrafish system, we show how this combination has effectively advanced classic topics in developmental biology compared to a perturbation-centric approach. Finally, we show that interdisciplinary expertise and perhaps specialization are necessary for carrying out a systematic approach, and discuss the technical hurdles.

Extracellular interactions and ligand degradation shape the nodal morphogen gradient

The correct distribution and activity of secreted signaling proteins called morphogens is required for many developmental processes. Nodal morphogens play critical roles in embryonic axis formation in many organisms. Models proposed to generate the Nodal gradient include diffusivity, ligand processing, and a temporal activation window. But how the Nodal morphogen gradient forms in vivo remains unclear. Here, we have measured in vivo for the first time, the binding affinity of Nodal ligands to their major cell surface receptor, Acvr2b, and to the Nodal inhibitor, Lefty, by fluorescence cross-correlation spectroscopy. We examined the diffusion coefficient of Nodal ligands and Lefty inhibitors in live zebrafish embryos by fluorescence correlation spectroscopy. We also investigated the contribution of ligand degradation to the Nodal gradient. We show that ligand clearance via degradation shapes the Nodal gradient and correlates with its signaling range. By computational simulations of gradient formation, we demonstrate that diffusivity, extra-cellular interactions, and selective ligand destruction collectively shape the Nodal morphogen gradient.

DOI:http://dx.doi.org/10.7554/eLife.13879.001

Animals develop from a single fertilized egg cell into multicellular organisms. This process requires chemical signals called “morphogens” that instruct the cells how to behave during development. The morphogens move across cells and tissues to form gradients of the signal. Cells then respond in different ways depending on how much of the signal they receive. This, in turn, depends on several factors: first, how quickly or slowly the signal moves; second, how well the morphogen binds to responding cells and other molecules in its path; and third, how much signal is lost or destroyed during the movement.

Many researchers study morphogen gradients in the transparent zebrafish, since it grows quickly and it is easy to see developmental changes. However, until now it was not fully clear how the well-known morphogen called Nodal moves in live zebrafish as they develop.

Wang, Wang et al. have now investigated how well Nodal signals bind to the surface of cells that receive the signal and to a molecule called “Lefty”, which is present in the same path and interferes with Nodal signals. Advanced techniques called fluorescence correlation and cross-correlation spectroscopy were used to measure Nodal signals at the level of single molecules in growing zebrafish. The experiments gave insights into how far Nodal signals move and remain active. The results showed that, in addition to Nodal diffusing and binding to receiving cells, one of the most important factors determining how far and quickly Nodal moves is its inactivation and destruction. Lastly, Wang, Wang et al. built computational models to test their observations from live zebrafish.

The current work was based on forcing zebrafish to produce molecules including Nodal at locations within the fish that normally do not make them. Therefore future experiments will aim to examine these molecules and their interactions when they are produced at their normal locations in the animal over time.

DOI:http://dx.doi.org/10.7554/eLife.13879.002

Visualizing retinoic acid morphogen gradients

Morphogens were originally defined as secreted signaling molecules that diffuse from local sources to form concentration gradients, which specify multiple cell fates. More recently morphogen gradients have been shown to incorporate a range of mechanisms including short-range signal activation, transcriptional/translational feedback, and temporal windows of target gene induction. Many critical cell-cell signals implicated in both embryonic development and disease, such as Wnt, fibroblast growth factor (Fgf), hedgehog (Hh), transforming growth factor beta (TGFb), and retinoic acid (RA), are thought to act as morphogens, but key information on signal propagation and ligand distribution has been lacking for most. The zebrafish provides unique advantages for genetics and imaging to address gradients during early embryonic stages when morphogens help establish major body axes. This has been particularly informative for RA, where RA response elements (RAREs) driving fluorescent reporters as well as Fluorescence Resonance Energy Transfer (FRET) reporters of receptor binding have provided evidence for gradients, as well as regulatory mechanisms that attenuate noise and enhance gradient robustness in vivo. Here we summarize available tools in zebrafish and discuss their utility for studying dynamic regulation of RA morphogen gradients, through combined experimental and computational approaches.

Reaction-diffusion processes at the nano- and microscales

The bottom-up fabrication of nano- and microscale structures from primary building blocks (molecules, colloidal particles) has made remarkable progress over the past two decades, but most research has focused on structural aspects, leaving our understanding of the dynamic and spatiotemporal aspects at a relatively primitive stage. In this Review, we draw inspiration from living cells to argue that it is now time to move beyond the generation of structures and explore dynamic processes at the nanoscale. We first introduce nanoscale self-assembly, self-organization and reaction-diffusion processes as essential features of cells. Then, we highlight recent progress towards designing and controlling these fundamental features of life in abiological systems. Specifically, we discuss examples of reaction-diffusion processes that lead to such outcomes as self-assembly, self-organization, unique nanostructures, chemical waves and dynamic order to illustrate their ubiquity within a unifying context of dynamic oscillations and energy dissipation. Finally, we suggest future directions for research on reaction-diffusion processes at the nano- and microscales that we find hold particular promise for a new understanding of science at the nanoscale and the development of new kinds of nanotechnologies for chemical transport, chemical communication and integration with living systems.

A Theoretical Model of Jigsaw-Puzzle Pattern Formation by Plant Leaf Epidermal Cells

Plant leaf epidermal cells exhibit a jigsaw puzzle-like pattern that is generated by interdigitation of the cell wall during leaf development. The contribution of two ROP GTPases, ROP2 and ROP6, to the cytoskeletal dynamics that regulate epidermal cell wall interdigitation has already been examined; however, how interactions between these molecules result in pattern formation remains to be elucidated. Here, we propose a simple interface equation model that incorporates both the cell wall remodeling activity of ROP GTPases and the diffusible signaling molecules by which they are regulated. This model successfully reproduces pattern formation observed in vivo, and explains the counterintuitive experimental results of decreased cellulose production and increased thickness. Our model also reproduces the dynamics of three-way cell wall junctions. Therefore, this model provides a possible mechanism for cell wall interdigitation formation in vivo.

Noggin4 is a long-range inhibitor of Wnt8 signalling that regulates head development in Xenopus laevis

Noggin4 is a Noggin family secreted protein whose molecular and physiological functions remain unknown. In this study, we demonstrate that in contrast to other Noggins, Xenopus laevis Noggin4 cannot antagonise BMP signalling; instead, it specifically binds to Wnt8 and inhibits the Wnt/β -catenin pathway. Live imaging demonstrated that Noggin4 diffusivity in embryonic tissues significantly exceeded that of other Noggins. Using the Fluorescence Recovery After Photobleaching (FRAP) assay and mathematical modelling, we directly estimated the affinity of Noggin4 for Wnt8 in living embryos and determined that Noggin4 fine-tune the Wnt8 posterior-to-anterior gradient. Our results suggest a role for Noggin4 as a unique, freely diffusing, long-range inhibitor of canonical Wnt signalling, thus explaining its ability to promote head development.

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

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

Self-Organization of Spatial Patterning in Human Embryonic Stem Cells

The developing embryo is a remarkable example of self-organization, where functional units are created in a complex spatiotemporal choreography. Recently, human embryonic stem cells (ESCs) have been used to recapitulate in vitro the self-organization programs that are executed in the embryo in vivo. This represents an unique opportunity to address self-organization in humans that is otherwise not addressable with current technologies. In this chapter, we review the recent literature on self-organization of human ESCs, with a particular focus on two examples: formation of embryonic germ layers and neural rosettes. Intriguingly, both activation and elimination of TGFβ signaling can initiate self-organization, albeit with different molecular underpinnings. We discuss the mechanisms underlying the formation of these structures in vitro and explore future challenges in the field.

High-throughput mathematical analysis identifies Turing networks for patterning with equally diffusing signals

The Turing reaction-diffusion model explains how identical cells can self-organize to form spatial patterns. It has been suggested that extracellular signaling molecules with different diffusion coefficients underlie this model, but the contribution of cell-autonomous signaling components is largely unknown. We developed an automated mathematical analysis to derive a catalog of realistic Turing networks. This analysis reveals that in the presence of cell-autonomous factors, networks can form a pattern with equally diffusing signals and even for any combination of diffusion coefficients. We provide a software (available at http://www.RDNets.com) to explore these networks and to constrain topologies with qualitative and quantitative experimental data. We use the software to examine the self-organizing networks that control embryonic axis specification and digit patterning. Finally, we demonstrate how existing synthetic circuits can be extended with additional feedbacks to form Turing reaction-diffusion systems. Our study offers a new theoretical framework to understand multicellular pattern formation and enables the wide-spread use of mathematical biology to engineer synthetic patterning systems.

Green-to-red primed conversion of Dendra2 using blue and red lasers

Primed photoconversion of Dendra2 in commonly available confocal and super-resolution microscopy setups equipped with blue and red lasers is described.

Using Confocal Analysis of Xenopus laevis to Investigate Modulators of Wnt and Shh Morphogen Gradients

This protocol describes a method to visualise ligands distributed across a field of cells. The ease of expressing exogenous proteins, together with the large size of their cells in early embryos, make Xenopus laevis a useful model for visualising GFP-tagged ligands. Synthetic mRNAs are efficiently translated after injection into early stage Xenopus embryos, and injections can be targeted to a single cell. When combined with a lineage tracer such as membrane tethered RFP, the injected cell (and its descendants) that are producing the overexpressed protein can easily be followed. This protocol describes a method for the production of fluorescently tagged Wnt and Shh ligands from injected mRNA. The methods involve the micro dissection of ectodermal explants (animal caps) and the analysis of ligand diffusion in multiple samples. By using confocal imaging, information about ligand secretion and diffusion over a field of cells can be obtained. Statistical analyses of confocal images provide quantitative data on the shape of ligand gradients. These methods may be useful to researchers who want to test the effects of factors that may regulate the shape of morphogen gradients.

A Temporal Window for Signal Activation Dictates the Dimensions of a Nodal Signaling Domain

Morphogen signaling is critical for the growth and patterning of tissues in embryos and adults, but how morphogen signaling gradients are generated in tissues remains controversial. The morphogen Nodal was proposed to form a long-range signaling gradient via a reaction-diffusion system, on the basis of differential diffusion rates of Nodal and its antagonist Lefty. Here we use a specific zebrafish Nodal biosensor combined with immunofluorescence for phosphorylated Smad2 to demonstrate that endogenous Nodal is unlikely to diffuse over a long range. Instead, short-range Nodal signaling activation in a temporal window is sufficient to determine the dimensions of the Nodal signaling domain. The size of this temporal window is set by the differentially timed production of Nodal and Lefty, which arises mainly from repression of Lefty translation by the microRNA miR-430. Thus, temporal information is transformed into spatial information to define the dimensions of the Nodal signaling domain and, consequently, to specify mesendoderm.

Tissue patterning and cellular mechanics

In development, cells organize into biological tissues through cell growth, migration, and differentiation. Globally, this process is dictated by a genetically encoded program in which secreted morphogens and cell-cell interactions prompt the adoption of unique cell fates. Yet, at its lowest level, development is achieved through the modification of cell-cell adhesion and actomyosin-based contractility, which set the level of tension within cells and dictate how they pack together into tissues. The regulation of tension within individual cells and across large groups of cells is a major driving force of tissue organization and the basis of all cell shape change and cell movement in development.

In vivo dynamics of skeletal muscle Dystrophin in zebrafish embryos revealed by improved FRAP analysis

Dystrophin forms an essential link between sarcolemma and cytoskeleton, perturbation of which causes muscular dystrophy. We analysed Dystrophin binding dynamics in vivo for the first time. Within maturing fibres of host zebrafish embryos, our analysis reveals a pool of diffusible Dystrophin and complexes bound at the fibre membrane. Combining modelling, an improved FRAP methodology and direct semi-quantitative analysis of bleaching suggests the existence of two membrane-bound Dystrophin populations with widely differing bound lifetimes: a stable, tightly bound pool, and a dynamic bound pool with high turnover rate that exchanges with the cytoplasmic pool. The three populations were found consistently in human and zebrafish Dystrophins overexpressed in wild-type or dmd(ta222a/ta222a) zebrafish embryos, which lack Dystrophin, and in Gt(dmd-Citrine)(ct90a) that express endogenously-driven tagged zebrafish Dystrophin. These results lead to a new model for Dystrophin membrane association in developing muscle, and highlight our methodology as a valuable strategy for in vivo analysis of complex protein dynamics.

Critical waves and the length problem of biology

It is pointed out that the mystery of how biological systems measure their lengths vanishes away if one premises that they have discovered a way to generate linear waves analogous to compressional sound. These can be used to detect length at either large or small scales using echo timing and fringe counting. It is shown that suitable linear chemical potential waves can, in fact, be manufactured by tuning to criticality conventional reaction-diffusion with a small number substance. Min oscillations in Escherichia coli are cited as precedent resonant length measurement using chemical potential waves analogous to laser detection. Mitotic structures in eukaryotes are identified as candidates for such an effect at higher frequency. The engineering principle is shown to be very general and functionally the same as that used by hearing organs.

Molecular bioelectricity: how endogenous voltage potentials control cell behavior and instruct pattern regulation in vivo

In addition to biochemical gradients and transcriptional networks, cell behavior is regulated by endogenous bioelectrical cues originating in the activity of ion channels and pumps, operating in a wide variety of cell types. Instructive signals mediated by changes in resting potential control proliferation, differentiation, cell shape, and apoptosis of stem, progenitor, and somatic cells. Of importance, however, cells are regulated not only by their own V_mem but also by the V_mem of their neighbors, forming networks via electrical synapses known as gap junctions. Spatiotemporal changes in V_mem distribution among nonneural somatic tissues regulate pattern formation and serve as signals that trigger limb regeneration, induce eye formation, set polarity of whole-body anatomical axes, and orchestrate craniofacial patterning. New tools for tracking and functionally altering V_mem gradients in vivo have identified novel roles for bioelectrical signaling and revealed the molecular pathways by which V_mem changes are transduced into cascades of downstream gene expression. Because channels and gap junctions are gated posttranslationally, bioelectrical networks have their own characteristic dynamics that do not reduce to molecular profiling of channel expression (although they couple functionally to transcriptional networks). The recent data provide an exciting opportunity to crack the bioelectric code, and learn to program cellular activity at the level of organs, not only cell types. The understanding of how patterning information is encoded in bioelectrical networks, which may require concepts from computational neuroscience, will have transformative implications for embryogenesis, regeneration, cancer, and synthetic bioengineering.

Molecular Framework of a Regulatory Circuit Initiating Two-Dimensional Spatial Patterning of Stomatal Lineage

Stomata, valves on the plant epidermis, are critical for plant growth and survival, and the presence of stomata impacts the global water and carbon cycle. Although transcription factors and cell-cell signaling components regulating stomatal development have been identified, it remains unclear as to how their regulatory interactions are translated into two-dimensional patterns of stomatal initial cells. Using molecular genetics, imaging, and mathematical simulation, we report a regulatory circuit that initiates the stomatal cell-lineage. The circuit includes a positive feedback loop constituting self-activation of SCREAMs that requires SPEECHLESS. This transcription factor module directly binds to the promoters and activates a secreted signal, EPIDERMAL PATTERNING FACTOR2, and the receptor modifier TOO MANY MOUTHS, while the receptor ERECTA lies outside of this module. This in turn inhibits SPCH, and hence SCRMs, thus constituting a negative feedback loop. Our mathematical model accurately predicts all known stomatal phenotypes with the inclusion of two additional components to the circuit: an EPF2-independent negative-feedback loop and a signal that lies outside of the SPCH•SCRM module. Our work reveals the intricate molecular framework governing self-organizing two-dimensional patterning in the plant epidermis.

Models for patterning primary embryonic body axes: The role of space and time

Models for the generation and interpretation of spatial patterns are discussed. Crucial for these processes is an intimate link between self-enhancing and antagonistic reactions. For spatial patterning, long-ranging antagonistic reactions are required that restrict the self-enhancing reactions to generate organizing regions. Self-enhancement is also required for a permanent switch-like activation of genes. This self-enhancement is antagonized by the mutual repression of genes, making sure that in a particular cell only one gene of a set of possible genes become activated - a long range inhibition in the 'gene space'. The understanding how the main body axes are initiated becomes more straightforward if the evolutionary ancestral head/brain pattern and the trunk pattern is considered separately. To activate a specific gene at particular concentration of morphogenetic gradient, observations are compatible with a systematic and time-requiring 'promotion' from one gene to the next until the local concentration is insufficient to accomplish a further promotion. The achieved determination is stable against a fading of the morphogen, as required to allow substantial growth. Minor modifications lead to a purely time-dependent activation of genes; both mechanisms are involved to pattern the anteroposterior axis. A mutual activation of cell states that locally exclude each other accounts for many features of the segmental patterning of the trunk. A possible scenario for the evolutionary invention of segmentation is discussed that is based on a reemployment of interactions involved in asexual reproduction.

Temporally coordinated signals progressively pattern the anteroposterior and dorsoventral body axes

The vertebrate body plan is established through the precise spatiotemporal coordination of morphogen signaling pathways that pattern the anteroposterior (AP) and dorsoventral (DV) axes. Patterning along the AP axis is directed by posteriorizing signals Wnt, fibroblast growth factor (FGF), Nodal, and retinoic acid (RA), while patterning along the DV axis is directed by bone morphogenetic proteins (BMP) ventralizing signals. This review addresses the current understanding of how Wnt, FGF, RA and BMP pattern distinct AP and DV cell fates during early development and how their signaling mechanisms are coordinated to concomitantly pattern AP and DV tissues.

In vivo single-cell labeling by confined primed conversion

Spatially confined green-to-red photoconversion of fluorescent proteins with high-power, pulsed laser illumination is negligible, thus precluding optical selection of single cells in vivo. We report primed conversion, in which low-power, dual-wavelength, continuous-wave illumination results in pronounced photoconversion. With a straightforward addition to a conventional confocal microscope, we show confined primed conversion in living zebrafish and reveal the complex anatomy of individual neurons packed between neighboring cells.

The role of source delocalization in the development of morphogen gradients

Successful biological development via spatial regulation of cell differentiation relies on the action of multiple signaling molecules that are known as morphogens. It is now well-established that signaling molecules create non-uniform concentration profiles, morphogen gradients, that activate different genes, leading to patterning in the developing embryos. The current view of the formation of morphogen gradients is that it is a result of complex reaction-diffusion processes that include the strongly localized production, diffusion and uniform degradation of signaling molecules. However, multiple experimental studies also suggest that the production of morphogen in many cases is delocalized. We develop a theoretical method that allows us to investigate the role of the delocalization in the formation of morphogen gradients. The approach is based on discrete-state stochastic models that can be solved exactly for arbitrary production lengths and production rates of morphogen molecules. Our analysis shows that the delocalization might have a strong effect on mechanisms of the morphogen gradient formation. The physical origin of this effect is discussed.

Response to Nodal morphogen gradient is determined by the kinetics of target gene induction

Morphogen gradients expose cells to different signal concentrations and induce target genes with different ranges of expression. To determine how the Nodal morphogen gradient induces distinct gene expression patterns during zebrafish embryogenesis, we measured the activation dynamics of the signal transducer Smad2 and the expression kinetics of long- and short-range target genes. We found that threshold models based on ligand concentration are insufficient to predict the response of target genes. Instead, morphogen interpretation is shaped by the kinetics of target gene induction: the higher the rate of transcription and the earlier the onset of induction, the greater the spatial range of expression. Thus, the timing and magnitude of target gene expression can be used to modulate the range of expression and diversify the response to morphogen gradients.

DOI:http://dx.doi.org/10.7554/eLife.05042.001

How a cell can tell where it is in a developing embryo has fascinated scientists for decades. The pioneering computer scientist and mathematical biologist Alan Turing was the first person to coin the term ‘morphogen’ to describe a protein that provides information about locations in the body. A morphogen is released from a group of cells (called the ‘source’) and as it moves away its activity (called the ‘signal’) declines gradually. Cells sense this signal gradient and use it to detect their position with respect to the source. Nodal is an important morphogen and is required to establish the correct identity of cells in the embryo; for example, it helps determine which cells should become a brain or heart or gut cell and so on.

The zebrafish is a widely used model to study animal development, in part because its embryos are transparent; this allows cells and proteins to be easily observed under a microscope. When Nodal acts on cells, another protein called Smad2 becomes activated, moves into the cell's nucleus, and then binds to specific genes. This triggers the expression of these genes, which are first copied into mRNA molecules via a process known as transcription and are then translated into proteins. The protein products of these targeted genes control cell identity and movement.

Several models have been proposed to explain how different concentrations of Nodal switch on the expression of different target genes; that is to say, to explain how a cell interprets the Nodal gradient. Dubrulle et al. have now measured factors that underlie how this gradient is interpreted. Individual cells in zebrafish embryos were tracked under a microscope, and Smad2 activation and gene expression were assessed. Dubrulle et al. found that, in contradiction to previous models, the amount of Nodal present on its own was insufficient to predict the target gene response. Instead, their analysis suggests that the size of each target gene's response depends on its rate of transcription and how quickly it is first expressed in response to Nodal.

These findings of Dubrulle et al. suggest that timing and transcription rate are important in determining the appropriate response to Nodal. Further work will be now needed to find out whether similar mechanisms regulate other processes that rely on the activity of morphogens.

DOI:http://dx.doi.org/10.7554/eLife.05042.002

Organ-level quorum sensing directs regeneration in hair stem cell populations

Coordinated organ behavior is crucial for an effective response to environmental stimuli. By studying regeneration of hair follicles in response to patterned hair plucking, we demonstrate that organ-level quorum sensing allows coordinated responses to skin injury. Plucking hair at different densities leads to a regeneration of up to five times more neighboring, unplucked resting hairs, indicating activation of a collective decision-making process. Through data modeling, the range of the quorum signal was estimated to be on the order of 1 mm, greater than expected for a diffusible molecular cue. Molecular and genetic analysis uncovered a two-step mechanism, where release of CCL2 from injured hairs leads to recruitment of TNF-α-secreting macrophages, which accumulate and signal to both plucked and unplucked follicles. By coupling immune response with regeneration, this mechanism allows skin to respond predictively to distress, disregarding mild injury, while meeting stronger injury with full-scale cooperative activation of stem cells.

Cooperativity to increase Turing pattern space for synthetic biology

It is hard to bridge the gap between mathematical formulations and biological implementations of Turing patterns, yet this is necessary for both understanding and engineering these networks with synthetic biology approaches. Here, we model a reaction-diffusion system with two morphogens in a monostable regime, inspired by components that we recently described in a synthetic biology study in mammalian cells.1 The model employs a single promoter to express both the activator and inhibitor genes and produces Turing patterns over large regions of parameter space, using biologically interpretable Hill function reactions. We applied a stability analysis and identified rules for choosing biologically tunable parameter relationships to increase the likelihood of successful patterning. We show how to control Turing pattern sizes and time evolution by manipulating the values for production and degradation relationships. More importantly, our analysis predicts that steep dose-response functions arising from cooperativity are mandatory for Turing patterns. Greater steepness increases parameter space and even reduces the requirement for differential diffusion between activator and inhibitor. These results demonstrate some of the limitations of linear scenarios for reaction-diffusion systems and will help to guide projects to engineer synthetic Turing patterns.