Visualisation and quantification of morphogen gradient formation in the zebrafish

Temperature effects on development and lifelong behavior in zebrafish

In recent decades, global warming has intensified temperature changes, placing substantial pressure on organism survival. Understanding how temperature variations impact development and behavior is crucial for conservation strategies. This study examined how temperature affects zebrafish embryo development and behavior, focusing on mRNA expression changes under thermal challenges. Zebrafish embryos were reared at 27 °C (control), 22 °C, and 30 °C, monitored from 24 to 120 hpf for structural development, and tested for optomotor responses at 7 dpf. Juvenile (30 dpf) and adult (90 dpf) fish reared at 27 °C were subjected to acute temperature shifts (22 °C and 30 °C for 2 h), followed by behavioral assessments and brain sampling for hsp90a and hspb1 mRNA expression analysis. Survival rates were significantly lower at 22 °C, with higher hatching rates at 30 °C but decreased at 22 °C. Developmental abnormalities varied: head malformations were more common at 30 °C, pericardial and yolk sac edema at 22 °C, and tail malformations at both extremes. Optomotor responses were impaired in fish from 22 °C. Social and aggressive behaviors were mostly unaffected, but fish from extreme temperatures showed increased risk-taking and reduced response to alarm substances. hsp90a mRNA expression was elevated in fish raised at 30 °C and those exposed to the 30 °C challenge, while hspb1 mRNA expression remained stable across temperatures. Cooling environments detrimentally affected embryo growth and survival, while warmer conditions induced pronounced growth defects. Elevated temperatures posed greater risks, triggering heightened hsp90a expression crucial for stress adaptation. Understanding thermal variation impacts on embryo development is crucial for mitigating climate change effects on species' viability and reproduction.

Phenotypic consequences of logarithmic signaling in MAPK stress response

How cells respond to dynamic environmental changes is crucial for understanding fundamental biological processes and cell physiology. In this study, we developed an experimental and quantitative analytical framework to explore how dynamic stress gradients that change over time regulate cellular volume, signaling activation, and growth phenotypes. Our findings reveal that gradual stress conditions substantially enhance cell growth compared to conventional acute stress. This growth advantage correlates with a minimal reduction in cell volume dependent on the dynamic of stress. We explain the growth phenotype with our finding of a logarithmic signal transduction mechanism in the yeast mitogen-activated protein kinase (MAPK) osmotic stress response pathway. These insights into the interplay between gradual environments, cell volume change, dynamic cell signaling, and growth, advance our understanding of fundamental cellular processes in gradual stress environments.

Tissue-specific and endogenous protein labeling with split fluorescent proteins

The ability to label proteins by fusion with genetically encoded fluorescent proteins is a powerful tool for understanding dynamic biological processes. However, current approaches for expressing fluorescent protein fusions possess drawbacks, especially at the whole organism level. Expression by transgenesis risks potential overexpression artifacts while fluorescent protein insertion at endogenous loci is technically difficult and, more importantly, does not allow for tissue-specific study of broadly expressed proteins. To overcome these limitations, we have adopted the split fluorescent protein system mNeonGreen21-10/11 (split-mNG2) to achieve tissue-specific and endogenous protein labeling in zebrafish. In our approach, mNG21-10 is expressed under a tissue-specific promoter using standard transgenesis while mNG211 is inserted into protein-coding genes of interest using CRISPR/Cas-directed gene editing. Each mNG2 fragment on its own is not fluorescent, but when co-expressed the fragments self-assemble into a fluorescent complex. Here, we report successful use of split-mNG2 to achieve differential labeling of the cytoskeleton genes tubb4b and krt8 in various tissues. We also demonstrate that by anchoring the mNG21-10 component to specific cellular compartments, the split-mNG2 system can be used to manipulate protein localization. Our approach should be broadly useful for a wide range of applications.

Tissue-specific and endogenous protein labeling with split fluorescent proteins

The ability to label proteins by fusion with genetically encoded fluorescent proteins is a powerful tool for understanding dynamic biological processes. However, current approaches for expressing fluorescent protein fusions possess drawbacks, especially at the whole organism level. Expression by transgenesis risks potential overexpression artifacts while fluorescent protein insertion at endogenous loci is technically difficult and, more importantly, does not allow for tissue-specific study of broadly expressed proteins. To overcome these limitations, we have adopted the split fluorescent protein system mNeonGreen21-10/11 (split-mNG2) to achieve tissue-specific and endogenous protein labeling in zebrafish. In our approach, mNG21-10 is expressed under a tissue-specific promoter using standard transgenesis while mNG211 is inserted into protein-coding genes of interest using CRISPR/Cas-directed gene editing. Each mNG2 fragment on its own is not fluorescent, but when co-expressed the fragments self-assemble into a fluorescent complex. Here, we report successful use of split-mNG2 to achieve differential labeling of the cytoskeleton genes tubb4b and krt8 in various tissues. We also demonstrate that by anchoring the mNG21-10 component to specific cellular compartments, the split-mNG2 system can be used to manipulate protein function. Our approach should be broadly useful for a wide range of applications.

Phenotypic consequences of logarithmic signaling in MAPK stress response

How cells respond to dynamic environmental changes is crucial for understanding fundamental biological processes and cell physiology. In this study, we developed an experimental and quantitative analytical framework to explore how dynamic stress gradients that change over time regulate cellular volume, signaling activation, and growth phenotypes. Our findings reveal that gradual stress conditions substantially enhance cell growth compared to conventional acute stress. This growth advantage correlates with a minimal reduction in cell volume dependent on the dynamic of stress. We explain the growth phenotype with our finding of a logarithmic signal transduction mechanism in the yeast Mitogen-Activated Protein Kinase (MAPK) osmotic stress response pathway. These insights into the interplay between gradual environments, cell volume change, dynamic cell signaling, and growth, advance our understanding of fundamental cellular processes in gradual stress environments.

Using FCS to accurately measure protein concentration in the presence of noise and photobleaching

Quantitative cell biology requires precise and accurate concentration measurements, resolved both in space and time. Fluorescence correlation spectroscopy (FCS) has been held as a promising technique to perform such measurements because the fluorescence fluctuations it relies on are directly dependent on the absolute number of fluorophores in the detection volume. However, the most interesting applications are in cells, where autofluorescence and confinement result in strong background noise and important levels of photobleaching. Both noise and photobleaching introduce systematic bias in FCS concentration measurements and need to be corrected for. Here, we propose to make use of the photobleaching inevitably occurring in confined environments to perform series of FCS measurements at different fluorophore concentration, which we show allows a precise in situ measurement of both background noise and molecular brightness. Such a measurement can then be used as a calibration to transform confocal intensity images into concentration maps. The power of this approach is first illustrated with in vitro measurements using different dye solutions, then its applicability for in vivo measurements is demonstrated in Drosophila embryos for a model nuclear protein and for two morphogens, Bicoid and Capicua.

The pattern of nodal morphogen signaling is shaped by co-receptor expression

Embryos must communicate instructions to their constituent cells over long distances. These instructions are often encoded in the concentration of signals called morphogens. In the textbook view, morphogen molecules diffuse from a localized source to form a concentration gradient, and target cells adopt fates by measuring the local morphogen concentration. However, natural patterning systems often incorporate numerous co-factors and extensive signaling feedback, suggesting that embryos require additional mechanisms to generate signaling patterns. Here, we examine the mechanisms of signaling pattern formation for the mesendoderm inducer Nodal during zebrafish embryogenesis. We find that Nodal signaling activity spans a normal range in the absence of signaling feedback and relay, suggesting that diffusion is sufficient for Nodal gradient formation. We further show that the range of endogenous Nodal ligands is set by the EGF-CFC co-receptor Oep: in the absence of Oep, Nodal activity spreads to form a nearly uniform distribution throughout the embryo. In turn, increasing Oep levels sensitizes cells to Nodal ligands. We recapitulate these experimental results with a computational model in which Oep regulates the diffusive spread of Nodal ligands by setting the rate of capture by target cells. This model predicts, and we confirm in vivo, the surprising observation that a failure to replenish Oep transforms the Nodal signaling gradient into a travelling wave. These results reveal that patterns of Nodal morphogen signaling are shaped by co-receptor-mediated restriction of ligand spread and sensitization of responding cells.

Generation of extracellular morphogen gradients: the case for diffusion

Cells within developing tissues rely on morphogens to assess positional information. Passive diffusion is the most parsimonious transport model for long-range morphogen gradient formation but does not, on its own, readily explain scaling, robustness and planar transport. Here, we argue that diffusion is sufficient to ensure robust morphogen gradient formation in a variety of tissues if the interactions between morphogens and their extracellular binders are considered. A current challenge is to assess how the affinity for extracellular binders, as well as other biophysical and cell biological parameters, determines gradient dynamics and shape in a diffusion-based transport system. Technological advances in genome editing, tissue engineering, live imaging and in vivo biophysics are now facilitating measurement of these parameters, paving the way for mathematical modelling and a quantitative understanding of morphogen gradient formation and modulation.

French flag gradients and Turing reaction-diffusion versus differentiation waves as models of morphogenesis

The Turing reaction-diffusion model and the French Flag Model are widely accepted in the field of development as the best models for explaining embryogenesis. Virtually all current attempts to understand cell differentiation in embryos begin and end with the assumption that some combination of these two models works. The result may become a bias in embryogenesis in assuming the problem has been solved by these two-chemical substance-based models. Neither model is applied consistently. We review the differences between the French Flag, Turing reaction-diffusion model, and a mechanochemical model called the differentiation wave/cell state splitter model. The cytoskeletal cell state splitter and the embryonic differentiation waves was first proposed in 1987 as a combined physics and chemistry model for cell differentiation in embryos, based on empirical observations on urodele amphibian embryos. We hope that the development of theory can be advanced and observations relevant to distinguishing the embryonic differentiation wave model from the French Flag model and reaction-diffusion equations will be taken up by experimentalists. Experimentalists rely on mathematical biologists for theory, and therefore depend on them for what parameters they choose to measure and ignore. Therefore, mathematical biologists need to fully understand the distinctions between these three models.

Zebrafish embryonic explants undergo genetically encoded self-assembly

Embryonic stem cell cultures are thought to self-organize into embryoid bodies, able to undergo symmetry-breaking, germ layer specification and even morphogenesis. Yet, it is unclear how to reconcile this remarkable self-organization capacity with classical experiments demonstrating key roles for extrinsic biases by maternal factors and/or extraembryonic tissues in embryogenesis. Here, we show that zebrafish embryonic tissue explants, prepared prior to germ layer induction and lacking extraembryonic tissues, can specify all germ layers and form a seemingly complete mesendoderm anlage. Importantly, explant organization requires polarized inheritance of maternal factors from dorsal-marginal regions of the blastoderm. Moreover, induction of endoderm and head-mesoderm, which require peak Nodal-signaling levels, is highly variable in explants, reminiscent of embryos with reduced Nodal signals from the extraembryonic tissues. Together, these data suggest that zebrafish explants do not undergo bona fide self-organization, but rather display features of genetically encoded self-assembly, where intrinsic genetic programs control the emergence of order.

Hmgcr promotes a long-range signal to attract Drosophila germ cells independently of Hedgehog

During development, many cell types migrate along stereotyped routes determined through deployment of cell surface or secreted guidance molecules. Although we know the identity of many of these molecules, the distances over which they natively operate can be difficult to determine. Here, we have quantified the range of an attractive signal for the migration of Drosophila germ cells. Their migration is guided by an attractive signal generated by the expression of genes in the 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (Hmgcr) pathway, and by a repulsive signal generated by the expression of Wunens. We demonstrate that the attractive signal downstream of Hmgcr is cell-contact independent and acts at long range, the extent of which depends on Hmgcr levels. This range would be sufficient to reach all of the germ cells for their entire migration. Furthermore, Hmgcr-mediated attraction does not require Wunens but can operate simultaneously with Wunen-mediated repulsion. Finally, several papers posit Hedgehog (Hh) as being the germ cell attractant downstream of H mgcr Here, we provide evidence that this is not the case.

Setting up for gastrulation in zebrafish

Soon after fertilization the zebrafish embryo generates the pool of cells that will give rise to the germline and the three somatic germ layers of the embryo (ectoderm, mesoderm and endoderm). As the basic body plan of the vertebrate embryo emerges, evolutionarily conserved developmental signaling pathways, including Bmp, Nodal, Wnt, and Fgf, direct the nearly totipotent cells of the early embryo to adopt gene expression profiles and patterns of cell behavior specific to their eventual fates. Several decades of molecular genetics research in zebrafish has yielded significant insight into the maternal and zygotic contributions and mechanisms that pattern this vertebrate embryo. This new understanding is the product of advances in genetic manipulations and imaging technologies that have allowed the field to probe the cellular, molecular and biophysical aspects underlying early patterning. The current state of the field indicates that patterning is governed by the integration of key signaling pathways and physical interactions between cells, rather than a patterning system in which distinct pathways are deployed to specify a particular cell fate. This chapter focuses on recent advances in our understanding of the genetic and molecular control of the events that impart cell identity and initiate the patterning of tissues that are prerequisites for or concurrent with movements of gastrulation.

Homodimerization regulates an endothelial specific signature of the SOX18 transcription factor

During embryogenesis, vascular development relies on a handful of transcription factors that instruct cell fate in a distinct sub-population of the endothelium (1). The SOXF proteins that comprise SOX7, 17 and 18, are molecular switches modulating arterio-venous and lymphatic endothelial differentiation (2,3). Here, we show that, in the SOX-F family, only SOX18 has the ability to switch between a monomeric and a dimeric form. We characterized the SOX18 dimer in binding assays in vitro, and using a split-GFP reporter assay in a zebrafish model system in vivo. We show that SOX18 dimerization is driven by a novel motif located in the vicinity of the C-terminus of the DNA binding region. Insertion of this motif in a SOX7 monomer forced its assembly into a dimer. Genome-wide analysis of SOX18 binding locations on the chromatin revealed enrichment for a SOX dimer binding motif, correlating with genes with a strong endothelial signature. Using a SOX18 small molecule inhibitor that disrupts dimerization, we revealed that dimerization is important for transcription. Overall, we show that dimerization is a specific feature of SOX18 that enables the recruitment of key endothelial transcription factors, and refines the selectivity of the binding to discrete genomic locations assigned to endothelial specific genes.

Imaging and Quantification of P-Smad1/5 in Zebrafish Blastula and Gastrula Embryos

Spatiotemporal patterns of morphogen activity drive differential gene expression with a high degree of precision within a developing embryo and reproducibly between embryos. Understanding the formation and function of a morphogen gradient during development requires quantitative measurement of morphogen activity throughout an individual embryo and also between embryos within a population. Quantification of morphogen gradients in to presents unique challenges in imaging and image processing to minimize error and maximize the quality of the data so it may be used in computational models of development and in statistically testing hypotheses. Here we present methods for the preparation, immunostaining, imaging, and quantification of a bone morphogenetic protein (BMP) activity gradient in individual zebrafish embryos as well as methods for acquiring population-level statistics after embryo grouping and alignment. This quantitative approach can be extended to other morphogen systems, and the computational codes can be adapted to other imaging contexts in zebrafish and other organisms.

Functional equivalence of germ plasm organizers

The proteins Oskar (Osk) in Drosophila and Bucky ball (Buc) in zebrafish act as germ plasm organizers. Both proteins recapitulate germ plasm activities but seem to be unique to their animal groups. Here, we discover that Osk and Buc show similar activities during germ cell specification. Drosophila Osk induces additional PGCs in zebrafish. Surprisingly, Osk and Buc do not show homologous protein motifs that would explain their related function. Nonetheless, we detect that both proteins contain stretches of intrinsically disordered regions (IDRs), which seem to be involved in protein aggregation. IDRs are known to rapidly change their sequence during evolution, which might obscure biochemical interaction motifs. Indeed, we show that Buc binds to the known Oskar interactors Vasa protein and nanos mRNA indicating conserved biochemical activities. These data provide a molecular framework for two proteins with unrelated sequence but with equivalent function to assemble a conserved core-complex nucleating germ plasm.

Multicellular organisms use gametes for their propagation. Gametes are formed from germ cells, which are specified during embryogenesis in some animals by the inheritance of RNP granules known as germ plasm. Transplantation of germ plasm induces extra germ cells, whereas germ plasm ablation leads to the loss of gametes and sterility. Therefore, germ plasm is key for germ cell formation and reproduction. However, the molecular mechanisms of germ cell specification by germ plasm in the vertebrate embryo remain an unsolved question. Proteins, which assemble the germ plasm, are known as germ plasm organizers. Here, we show that the two germ plasm organizers Oskar from the fly and Bucky ball from the fish show similar functions by using a cross species approach. Both are intrinsically disordered proteins, which rapidly changed their sequence during evolution. Moreover, both proteins still interact with conserved components of the germ cell specification pathway. These data might provide a first example of two proteins with the same biological role, but distinct sequence.

TGF-β Family Signaling in Early Vertebrate Development

TGF-β family ligands function in inducing and patterning many tissues of the early vertebrate embryonic body plan. Nodal signaling is essential for the specification of mesendodermal tissues and the concurrent cellular movements of gastrulation. Bone morphogenetic protein (BMP) signaling patterns tissues along the dorsal-ventral axis and simultaneously directs the cell movements of convergence and extension. After gastrulation, a second wave of Nodal signaling breaks the symmetry between the left and right sides of the embryo. During these processes, elaborate regulatory feedback between TGF-β ligands and their antagonists direct the proper specification and patterning of embryonic tissues. In this review, we summarize the current knowledge of the function and regulation of TGF-β family signaling in these processes. Although we cover principles that are involved in the development of all vertebrate embryos, we focus specifically on three popular model organisms: the mouse Mus musculus, the African clawed frog of the genus Xenopus, and the zebrafish Danio rerio, highlighting the similarities and differences between these species.

Timing of gene expression in a cell-fate decision system

During development, morphogens provide extracellular cues allowing cells to select a specific fate by inducing complex transcriptional programs. The mating pathway in budding yeast offers simplified settings to understand this process. Pheromone secreted by the mating partner triggers the activity of a MAPK pathway, which results in the expression of hundreds of genes. Using a dynamic expression reporter, we quantified the kinetics of gene expression in single cells upon exogenous pheromone stimulation and in the physiological context of mating. In both conditions, we observed striking differences in the timing of induction of mating‐responsive promoters. Biochemical analyses and generation of synthetic promoter variants demonstrated how the interplay between transcription factor binding and nucleosomes contributes to determine the kinetics of transcription in a simplified cell‐fate decision system.

Yeast-based assays for detecting protein-protein/drug interactions and their inhibitors

Understanding cellular processes at molecular levels in health and disease requires the knowledge of protein-protein interactions (PPIs). In line with this, identification of PPIs at genome-wide scale is highly valuable to understand how different cellular pathways are interconnected, and it eventually facilitates designing effective drugs against certain PPIs. Furthermore, investigating PPIs at a small laboratory scale for deciphering certain biochemical pathways has been demanded for years. In this regard, yeast two hybrid system (Y2HS) has proven an extremely useful tool to discover novel PPIs, while Y2HS derivatives and novel yeast-based assays are contributing significantly to identification of protein-drug/inhibitor interaction at both large- and small-scale set-ups. These methods have been evolving over time to provide more accurate, reproducible and quantitative results. Here we briefly describe different yeast-based assays for identification of various protein-protein/drug/inhibitor interactions and their specific applications, advantages, shortcomings, and improvements. The broad range of yeast-based assays facilitates application of the most suitable method(s) for each specific need.

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.

Nuclear movement regulated by non-Smad Nodal signaling via JNK is associated with Smad signaling during zebrafish endoderm specification

Although asymmetric nuclear positioning is observed during animal development, the regulation and significance of this nuclear positioning in cell differentiation remains poorly understood. Using zebrafish blastulae, we provide evidence that nuclear movement toward the yolk syncytial layer, which comprises extraembryonic tissue, occurs in the first endoderm specified cells during endoderm specification. Nodal signaling is essential for nuclear movement, whereas nuclear envelope proteins are involved in the movement through the microtubule formation. The positioning of the microtubule organizing center, which is proposed to be critical for nuclear movement, is regulated by Nodal signaling and nuclear envelope proteins. The non-Smad JNK signaling pathway, which is downstream of Nodal signaling, regulates nuclear movement independent of the Smad pathway, and this nuclear movement is associated with Smad signal transduction toward the nucleus. Our study provides insights into the function of nuclear movement in Smad signaling toward the nucleus, and could be applied to the control of Transforming Growth Factor-β signaling.

Zebrafish Pronephros Development

The pronephros is the first kidney type to form in vertebrate embryos. The first step of pronephrogenesis in the zebrafish is the formation of the intermediate mesoderm during gastrulation, which occurs in response to secreted morphogens such as BMPs and Nodals. Patterning of the intermediate mesoderm into proximal and distal cell fates is induced by retinoic acid signaling with downstream transcription factors including wt1a, pax2a, pax8, hnf1b, sim1a, mecom, and irx3b. In the anterior intermediate mesoderm, progenitors of the glomerular blood filter migrate and fuse at the midline and recruit a blood supply. More posteriorly localized tubule progenitors undergo epithelialization and fuse with the cloaca. The Notch signaling pathway regulates the formation of multi-ciliated cells in the tubules and these cells help propel the filtrate to the cloaca. The lumenal sheer stress caused by flow down the tubule activates anterior collective migration of the proximal tubules and induces stretching and proliferation of the more distal segments. Ultimately these processes create a simple two-nephron kidney that is capable of reabsorbing and secreting solutes and expelling excess water-processes that are critical to the homeostasis of the body fluids. The zebrafish pronephric kidney provides a simple, yet powerful, model system to better understand the conserved molecular and cellular progresses that drive nephron formation, structure, and function.

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.

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.

Molecular specification of germ layers in vertebrate embryos

In order to generate the tissues and organs of a multicellular organism, different cell types have to be generated during embryonic development. The first step in this process of cellular diversification is the formation of the three germ layers: ectoderm, endoderm and mesoderm. The ectoderm gives rise to the nervous system, epidermis and various neural crest-derived tissues, the endoderm goes on to form the gastrointestinal, respiratory and urinary systems as well as many endocrine glands, and the mesoderm will form the notochord, axial skeleton, cartilage, connective tissue, trunk muscles, kidneys and blood. Classic experiments in amphibian embryos revealed the tissue interactions involved in germ layer formation and provided the groundwork for the identification of secreted and intracellular factors involved in this process. We will begin this review by summarising the key findings of those studies. We will then evaluate them in the light of more recent genetic studies that helped clarify which of the previously identified factors are required for germ layer formation in vivo, and to what extent the mechanisms identified in amphibians are conserved across other vertebrate species. Collectively, these studies have started to reveal the gene regulatory network (GRN) underlying vertebrate germ layer specification and we will conclude our review by providing examples how our understanding of this GRN can be employed to differentiate stem cells in a targeted fashion for therapeutic purposes.

Spatiotemporal analysis of different mechanisms for interpreting morphogen gradients

During development, multicellular organisms must accurately control both temporal and spatial aspects of tissue patterning. This is often achieved using morphogens, signaling molecules that form spatially varying concentrations and so encode positional information. Typical analysis of morphogens assumes that spatial information is decoded in steady state by measuring the value of the morphogen concentration. However, recent experimental work suggests that both pre-steady-state readout and measurement of spatial and temporal derivatives of the morphogen concentration can play important roles in defining boundaries. Here, we undertake a detailed theoretical and numerical study of the accuracy of patterning-both in space and time-in models where readout is provided not by the morphogen concentration but by its spatial and temporal derivatives. In both cases we find that accurate patterning can be achieved, with sometimes even smaller errors than directly reading the morphogen concentration. We further demonstrate that such models provide other potential benefits to the system, such as the ability to switch on and off gene response with a high degree of spatiotemporal accuracy. Finally, we discuss how such derivatives might be calculated biologically and examine these models in relation to Sonic Hedgehog signaling in the vertebrate central nervous system. We show that, when coupled to a downstream transcriptional network, pre-steady-state measurement of the temporal change in the Shh morphogen is a plausible mechanism for determining precise gene boundaries in both space and time.

Automated detection and quantification of single RNAs at cellular resolution in zebrafish embryos

Analysis of differential gene expression is crucial for the study of cell fate and behavior during embryonic development. However, automated methods for the sensitive detection and quantification of RNAs at cellular resolution in embryos are lacking. With the advent of single-molecule fluorescence in situ hybridization (smFISH), gene expression can be analyzed at single-molecule resolution. However, the limited availability of protocols for smFISH in embryos and the lack of efficient image analysis pipelines have hampered quantification at the (sub)cellular level in complex samples such as tissues and embryos. Here, we present a protocol for smFISH on zebrafish embryo sections in combination with an image analysis pipeline for automated transcript detection and cell segmentation. We use this strategy to quantify gene expression differences between different cell types and identify differences in subcellular transcript localization between genes. The combination of our smFISH protocol and custom-made, freely available, analysis pipeline will enable researchers to fully exploit the benefits of quantitative transcript analysis at cellular and subcellular resolution in tissues and embryos.

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.

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.

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

Activin/Nodal signalling in stem cells

Activin/Nodal growth factors control a broad range of biological processes, including early cell fate decisions, organogenesis and adult tissue homeostasis. Here, we provide an overview of the mechanisms by which the Activin/Nodal signalling pathway governs stem cell function in these different stages of development. We describe recent findings that associate Activin/Nodal signalling to pathological conditions, focusing on cancer stem cells in tumorigenesis and its potential as a target for therapies. Moreover, we will discuss future directions and questions that currently remain unanswered on the role of Activin/Nodal signalling in stem cell self-renewal, differentiation and proliferation.

Late Alk4/5/7 signaling is required for anterior skeletal patterning in sea urchin embryos

Skeletal patterning in the sea urchin embryo requires a conversation between the skeletogenic primary mesenchyme cells (PMCs) and the overlying pattern-dictating ectoderm; however, our understanding of the molecular basis for this process remains incomplete. Here, we show that TGF-β-receptor signaling is required during gastrulation to pattern the anterior skeleton. To block TGF-β signaling, we used SB431542 (SB43), a specific inhibitor of the TGF-β type I receptor Alk4/5/7. Treatment with SB43 during gastrulation blocks anterior PMC positioning and the formation of the anterior skeleton, but does not perturb general ectoderm specification or development. This is the first example of a signaling event required for patterning of a specific part of the skeleton. Alk4/5/7 inhibition does not prevent the formation of a mouth, although SB43-treated plutei display reduced feeding ability, presumably due to the loss of the structural support for the mouth conferred by the anterior skeleton. Both Univin and Nodal are potential ligands for Alk4/5/7; however, Nodal is unilaterally expressed on only the right side, whereas Univin is bilaterally expressed in the ectoderm adjacent to the anterior skeleton during the relevant time period. Our results demonstrate that Univin is both necessary and sufficient for secondary skeletal development in a control background, consistent with the hypothesis that Univin is a relevant Alk4/5/7 ligand for anterior skeletal patterning. Taken together, our data demonstrate that Alk4/5/7 signaling during gastrulation is required to direct PMCs to the oral hood, and suggest that Univin is a relevant ligand for this signaling event.

Encoding of temporal signals by the TGF-β pathway and implications for embryonic patterning

Genetics and biochemistry have defined the components and wiring of the signaling pathways that pattern the embryo. Among them, the transforming growth factor β (TGF-β) pathway has the potential to behave as a morphogen: in vitro experiments established that it can dictate cell fate in a concentration-dependent manner. How morphogens convey positional information in a developing embryo, when signal levels change with time, is less understood. Using integrated microfluidic cell culture and time-lapse microscopy, we demonstrate here that the speed of ligand presentation has a key and previously unexpected influence on TGF-β signaling outcomes. The response to a TGF-β concentration step is transient and adaptive: slowly increasing the ligand concentration diminishes the response, and well-spaced pulses of ligand combine additively, resulting in greater pathway output than with constant stimulation. Our results suggest that in an embryonic context, the speed of change of ligand concentration is an instructive signal for patterning.

The integrator complex subunit 6 (Ints6) confines the dorsal organizer in vertebrate embryogenesis

Dorsoventral patterning of the embryonic axis relies upon the mutual antagonism of competing signaling pathways to establish a balance between ventralizing BMP signaling and dorsal cell fate specification mediated by the organizer. In zebrafish, the initial embryo-wide domain of BMP signaling is refined into a morphogenetic gradient following activation dorsally of a maternal Wnt pathway. The accumulation of β-catenin in nuclei on the dorsal side of the embryo then leads to repression of BMP signaling dorsally and the induction of dorsal cell fates mediated by Nodal and FGF signaling. A separate Wnt pathway operates zygotically via Wnt8a to limit dorsal cell fate specification and maintain the expression of ventralizing genes in ventrolateral domains. We have isolated a recessive dorsalizing maternal-effect mutation disrupting the gene encoding Integrator Complex Subunit 6 (Ints6). Due to widespread de-repression of dorsal organizer genes, embryos from mutant mothers fail to maintain expression of BMP ligands, fail to fully express vox and ved, two mediators of Wnt8a, display delayed cell movements during gastrulation, and severe dorsalization. Consistent with radial dorsalization, affected embryos display multiple independent axial domains along with ectopic dorsal forerunner cells. Limiting Nodal signaling or restoring BMP signaling restores wild-type patterning to affected embryos. Our results are consistent with a novel role for Ints6 in restricting the vertebrate organizer to a dorsal domain in embryonic patterning.

Morphogenetic action through flux-limited spreading

A central question in biology is how secreted morphogens act to induce different cellular responses within a group of cells in a concentration-dependent manner. Modeling morphogenetic output in multicellular systems has so far employed linear diffusion, which is the normal type of diffusion associated with Brownian processes. However, there is evidence that at least some morphogens, such as Hedgehog (Hh) molecules, may not freely diffuse. Moreover, the mathematical analysis of such models necessarily implies unrealistic instantaneous spreading of morphogen molecules, which are derived from the assumptions of Brownian motion in its continuous formulation. A strict mathematical model considering Fick's diffusion law predicts morphogen exposure of the whole tissue at the same time. Such a strict model thus does not describe true biological patterns, even if similar and attractive patterns appear as results of applying such simple model. To eliminate non-biological behaviors from diffusion models we introduce flux-limited spreading (FLS), which implies a restricted velocity for morphogen propagation and a nonlinear mechanism of transport. Using FLS and focusing on intercellular Hh-Gli signaling, we model a morphogen gradient and highlight the propagation velocity of morphogen particles as a new key biological parameter. This model is then applied to the formation and action of the Sonic Hh (Shh) gradient in the vertebrate embryonic neural tube using our experimental data on Hh spreading in heterologous systems together with published data. Unlike linear diffusion models, FLS modeling predicts concentration fronts and the evolution of gradient dynamics and responses over time. In addition to spreading restrictions by extracellular binding partners, we suggest that the constraints imposed by direct bridges of information transfer such as nanotubes or cytonemes underlie FLS. Indeed, we detect and measure morphogen particle velocity in such cell extensions in different systems.

Quantitative imaging of morphogen gradients in Drosophila imaginal discs

Cells at different positions in a developing tissue receive different concentrations of signaling molecules, called morphogens, and this influences their cell fate. Morphogen concentration gradients have been proposed to control patterning as well as growth in many developing tissues. Some outstanding questions about tissue patterning by morphogen gradients are the following: What are the mechanisms that regulate gradient formation and shape? Is the positional information encoded in the gradient sufficiently precise to determine the positions of target gene domain boundaries? What are the temporal dynamics of gradients and how do they relate to patterning and growth? These questions are inherently quantitative in nature and addressing them requires measuring morphogen concentrations in cells, levels of downstream signaling activity, and kinetics of morphogen transport. Here we first present methods for quantifying morphogen gradient shape in which the measurements can be calibrated to reflect actual morphogen concentrations. We then discuss using fluorescence recovery after photobleaching to study the kinetics of morphogen transport at the tissue level. Finally, we present particle tracking as a method to study morphogen intracellular trafficking.

A role for Vg1/Nodal signaling in specification of the intermediate mesoderm

The intermediate mesoderm (IM) is the embryonic source of all kidney tissue in vertebrates. The factors that regulate the formation of the IM are not yet well understood. Through investigations in the chick embryo, the current study identifies and characterizes Vg1/Nodal signaling (henceforth referred to as 'Nodal-like signaling') as a novel regulator of IM formation. Excess Nodal-like signaling at gastrulation stages resulted in expansion of the IM at the expense of the adjacent paraxial mesoderm, whereas inhibition of Nodal-like signaling caused repression of IM gene expression. IM formation was sensitive to levels of the Nodal-like pathway co-receptor Cripto and was inhibited by a truncated form of the secreted molecule cerberus, which specifically blocks Nodal, indicating that the observed effects are specific to the Nodal-like branch of the TGFβ signaling pathway. The IM-promoting effects of Nodal-like signaling were distinct from the known effects of this pathway on mesoderm formation and left-right patterning, a finding that can be attributed to specific time windows for the activities of these Nodal-like functions. Finally, a link was observed between Nodal-like and BMP signaling in the induction of IM. Activation of IM genes by Nodal-like signaling required an active BMP signaling pathway, and Nodal-like signals induced phosphorylation of Smad1/5/8, which is normally associated with activation of BMP signaling pathways. We postulate that Nodal-like signaling regulates IM formation by modulating the IM-inducing effects of BMP signaling.

Morphogen transport

The graded distribution of morphogens underlies many of the tissue patterns that form during development. How morphogens disperse from a localized source and how gradients in the target tissue form has been under debate for decades. Recent imaging studies and biophysical measurements have provided evidence for various morphogen transport models ranging from passive mechanisms, such as free or hindered extracellular diffusion, to cell-based dispersal by transcytosis or cytonemes. Here, we analyze these transport models using the morphogens Nodal, fibroblast growth factor and Decapentaplegic as case studies. We propose that most of the available data support the idea that morphogen gradients form by diffusion that is hindered by tortuosity and binding to extracellular molecules.

Diversity in genetic in vivo methods for protein-protein interaction studies: from the yeast two-hybrid system to the mammalian split-luciferase system

The yeast two-hybrid system pioneered the field of in vivo protein-protein interaction methods and undisputedly gave rise to a palette of ingenious techniques that are constantly pushing further the limits of the original method. Sensitivity and selectivity have improved because of various technical tricks and experimental designs. Here we present an exhaustive overview of the genetic approaches available to study in vivo binary protein interactions, based on two-hybrid and protein fragment complementation assays. These methods have been engineered and employed successfully in microorganisms such as Saccharomyces cerevisiae and Escherichia coli, but also in higher eukaryotes. From single binary pairwise interactions to whole-genome interactome mapping, the self-reassembly concept has been employed widely. Innovative studies report the use of proteins such as ubiquitin, dihydrofolate reductase, and adenylate cyclase as reconstituted reporters. Protein fragment complementation assays have extended the possibilities in protein-protein interaction studies, with technologies that enable spatial and temporal analyses of protein complexes. In addition, one-hybrid and three-hybrid systems have broadened the types of interactions that can be studied and the findings that can be obtained. Applications of these technologies are discussed, together with the advantages and limitations of the available assays.

Genome-wide mechanisms of Smad binding

A dual role of transforming growth factor β (TGF-β), to both suppress and promote tumor progression and metastasis, has been well established, but its molecular basis has remained elusive. In this review, we focus on Smad proteins, which are central mediators of the signal transduction of TGF-β family members. We describe current knowledge of cell-type-specific binding patterns of Smad proteins and mechanisms of transcriptional regulation, obtained from recent studies on genome-wide binding sites of Smad molecules. We also discuss potential application of the genome-wide analyses for cancer research, which will allow clarification of the complex mechanisms occurring during cancer progression, and the identification of potential biomarkers for future cancer diagnosis, prognosis and therapy.

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

Biological systems involving short-range activators and long-range inhibitors can generate complex patterns. Reaction-diffusion models postulate that differences in signaling range are caused by differential diffusivity of inhibitor and activator. Other models suggest that differential clearance underlies different signaling ranges. To test these models, we measured the biophysical properties of the Nodal/Lefty activator/inhibitor system during zebrafish embryogenesis. Analysis of Nodal and Lefty gradients revealed that Nodals have a shorter range than Lefty proteins. Pulse-labeling analysis indicated that Nodals and Leftys have similar clearance kinetics, whereas fluorescence recovery assays revealed that Leftys have a higher effective diffusion coefficient than Nodals. These results indicate that differential diffusivity is the major determinant of the differences in Nodal/Lefty range and provide biophysical support for reaction-diffusion models of activator/inhibitor-mediated patterning.

Gene regulatory logic for reading the Sonic Hedgehog signaling gradient in the vertebrate neural tube

Secreted signals, known as morphogens, provide the positional information that organizes gene expression and cellular differentiation in many developing tissues. In the vertebrate neural tube, Sonic Hedgehog (Shh) acts as a morphogen to control the pattern of neuronal subtype specification. Using an in vivo reporter of Shh signaling, mouse genetics, and systems modeling, we show that a spatially and temporally changing gradient of Shh signaling is interpreted by the regulatory logic of a downstream transcriptional network. The design of the network, which links three transcription factors to Shh signaling, is responsible for differential spatial and temporal gene expression. In addition, the network renders cells insensitive to fluctuations in signaling and confers hysteresis--memory of the signal. Our findings reveal that morphogen interpretation is an emergent property of the architecture of a transcriptional network that provides robustness and reliability to tissue patterning.

Brain regionalization: of signaling centers and boundaries

Our knowledge of the general mechanisms controlling the formation of the vertebrate central nervous system has advanced tremendously in the last decade. Here, we discuss the impact of the combined use of cell manipulation, in vivo imaging and genetics in the zebrafish on recent progress in understanding how signaling processes progressively control regionalization of the central nervous system. We highlight the unresolved issues and speculate upon the fundamental role the zebrafish will continue having in answering them.

Proteins on the move: insights gained from fluorescent protein technologies

Proteins are always on the move, and this may occur through diffusion or active transport. The realization that the regulation of signal transduction is highly dynamic in space and time has stimulated intense interest in the movement of proteins. Over the past decade, numerous new technologies using fluorescent proteins have been developed, allowing us to observe the spatiotemporal dynamics of proteins in living cells. These technologies have greatly advanced our understanding of protein dynamics, including protein movement and protein interactions.

Morphogen gradients: from generation to interpretation

Morphogens are long-range signaling molecules that pattern developing tissues in a concentration-dependent manner. The graded activity of morphogens within tissues exposes cells to different signal levels and leads to region-specific transcriptional responses and cell fates. In its simplest incarnation, a morphogen signal forms a gradient by diffusion from a local source and clearance in surrounding tissues. Responding cells often transduce morphogen levels in a linear fashion, which results in the graded activation of transcriptional effectors. The concentration-dependent expression of morphogen target genes is achieved by their different binding affinities for transcriptional effectors as well as inputs from other transcriptional regulators. Morphogen distribution and interpretation are the result of complex interactions between the morphogen and responding tissues. The response to a morphogen is dependent not simply on morphogen concentration but also on the duration of morphogen exposure and the state of the target cells. In this review, we describe the morphogen concept and discuss the mechanisms that underlie the generation, modulation, and interpretation of morphogen gradients.

Extracellular movement of signaling molecules

Extracellular signaling molecules have crucial roles in development and homeostasis, and their incorrect deployment can lead to developmental defects and disease states. Signaling molecules are released from sending cells, travel to target cells, and act over length scales of several orders of magnitude, from morphogen-mediated patterning of small developmental fields to hormonal signaling throughout the organism. We discuss how signals are modified and assembled for transport, which routes they take to reach their targets, and how their range is affected by mobility and stability.

Steady-state invariant genetics: probing the role of morphogen gradient dynamics in developmental patterning

Morphogen-mediated patterning is the predominant mechanism by which positional information is established during animal development. In the classical view, the interpretation of positional signals depends on the equilibrium distribution of a morphogen, regardless of the dynamics of gradient formation. The problem of whether or not morphogen dynamics contribute to developmental patterning has not been explored in detail, partly because genetic experiments, which selectively affect signalling dynamics while maintaining unchanged the steady-state morphogen profile, are difficult to design and interpret. Here, I present a modelling-based approach to identify genetic mutations in developmental patterning that may affect the transient, but leave invariant the steady-state signalling gradient. As a case study, this approach is used to explore the dynamic properties of Hedgehog (Hh) signalling in the developing wing of the fruitfly, Drosophila melanogaster. This analysis provides insights into how different properties of the Hh gradient dynamics, such as the duration of exposure to the signal or the maximum width of the transient gradient, can be genetically perturbed without affecting the steady-state distribution of the Hh concentration profile. I propose that this method can be used as an experimental design tool to investigate the role of transient morphogen gradients in developmental patterning and discuss the generality of these ideas in other problems.

Integration of Hedgehog and BMP signalling by the engrailed2a gene in the zebrafish myotome

Different levels and timing of Hedgehog (Hh) signalling activity have been proposed to specify three distinct cell types in the zebrafish myotome. Two of these, the medial fast-twitch fibres (MFFs) and the slow-twitch muscle pioneers (MPs) are characterised by expression of eng1a, -1b and -2a and require the highest levels of Hh for their specification. We have defined a minimal eng2a element sufficient to drive reporter expression specifically in MPs and MFFs. This element binds both Gli2a, a mediator of Hh signalling, and activated Smads (pSmads), mediators of bone morphogenic protein (BMP) signalling, in vivo. We found a strict negative correlation between nuclear accumulation of pSmad, and eng2a expression in myotomal cells and show that abrogation of pSmad accumulation results in activation of eng2a, even when Hh signalling is attenuated. Conversely, driving nuclear accumulation of pSmad suppresses the induction of eng expression even when Hh pathway activity is maximal. Nuclear accumulation of pSmads is depleted by maximal Hh pathway activation. We show that a synthetic form of the Gli2 repressor interacts with Smad1 specifically in the nuclei of myotomal cells in the developing embryo and that this interaction depends upon BMP signalling activity. Our results demonstrate that the eng2a promoter integrates repressive and activating signals from the BMP and Hh pathways, respectively, to limit its expression to MPs and MFFs. We suggest a novel basis for crosstalk between the Hh and BMP pathways, whereby BMP-mediated repression of Hh target genes is promoted by a direct interaction between Smads and truncated Glis, an interaction that is abrogated by Hh induced depletion of the latter.

What's all the noise about developmental stochasticity?

In October 2010, researchers from diverse backgrounds collided at the historic Cumberland Lodge (Windsor, UK) to discuss the role of randomness in cell and developmental biology. Organized by James Briscoe and Alfonso Marinez-Arias, The Company of Biologists' workshop was the latest in a series of meetings aimed at encouraging interdisciplinary interactions between biologists. This aim was reflected in talks at this workshop that ranged from the tissue to the cellular scale, and that integrated experimental and theoretical approaches to examining stochastic behavior in diverse systems.