Synchronization of Hes1 oscillations coordinates and refines condensation formation and patterning of the avian limb skeleton

Form, function, mind: What doesn't compute (and what might)

The applicability of computational and dynamical systems models to organisms is scrutinized, using examples from developmental biology and cognition. Developmental morphogenesis is dependent on the inherent material properties of developing animal (metazoan) tissues, a non-computational modality, but cell differentiation, which utilizes chromatin-based revisable memory banks and program-like function-calling, via the developmental gene co-expression system unique to the metazoans, has a quasi-computational basis. Multi-attractor dynamical models are argued to be misapplied to global properties of development, and it is suggested that along with computationalism, classic forms of dynamicism are similarly unsuitable to accounting for cognitive phenomena. Proposals are made for treating brains and other nervous tissues as novel forms of excitable matter with inherent properties which enable the intensification of cell-based basal cognition capabilities present throughout the tree of life. Finally, some connections are drawn between the viewpoint described here and active inference models of cognition, such as the Free Energy Principle.

A Cellular Potts Model of the interplay of synchronization and aggregation

We investigate the behavior of systems of cells with intracellular molecular oscillators ("clocks") where cell-cell adhesion is mediated by differences in clock phase between neighbors. This is motivated by phenomena in developmental biology and in aggregative multicellularity of unicellular organisms. In such systems, aggregation co-occurs with clock synchronization. To account for the effects of spatially extended cells, we use the Cellular Potts Model (CPM), a lattice agent-based model. We find four distinct possible phases: global synchronization, local synchronization, incoherence, and anti-synchronization (checkerboard patterns). We characterize these phases via order parameters. In the case of global synchrony, the speed of synchronization depends on the adhesive effects of the clocks. Synchronization happens fastest when cells in opposite phases adhere the strongest ("opposites attract"). When cells of the same clock phase adhere the strongest ("like attracts like"), synchronization is slower. Surprisingly, the slowest synchronization happens in the diffusive mixing case, where cell-cell adhesion is independent of clock phase. We briefly discuss potential applications of the model, such as pattern formation in the auditory sensory epithelium.

A two-galectin network establishes mesenchymal condensation phenotype in limb development

Previous work showed that Gal-1A and Gal-8, two proteins belonging to the galactoside-binding galectin family, are the earliest determinants of the patterning of the skeletal elements of embryonic chicken limbs, and further, that their experimentally determined interactions in the embryonic limb bud can be interpreted via a reaction-diffusion-adhesion (2GL: two galectin plus ligands) model. Here, we use an ordinary differential equation-based approach to analyze the intrinsic switching modality of the 2GL network and characterize the network behavior independent of the diffusive and adhesive arms of the patterning mechanism. We identify two states: where the concentrations of both the galectins are respectively, negligible, and very high. This bistable switch-like system arises via a saddle-node bifurcation from a monostable state. For the case of mass-action production terms, we provide an explicit Lyapunov function for the system, which shows that it has no periodic solutions. Our model therefore predicts that the galectin network may exist in low expression and high expression states separated in space or time, without any intermediate states. We test these predictions in experiments performed with high density cultures of chick limb mesenchymal cells and observe that cells inside precartilage protocondensations express Gal-1A at a much higher rate than those outside, for which it was negligible. The Gal-1A and -8-based patterning network is therefore sufficient to partition the mesenchymal cell population into two discrete cell states with different developmental (chondrogenic vs. non-chondrogenic) fates. When incorporated into an adhesion and diffusion-enabled framework this system can generate a spatially patterned limb skeleton.

The vertebrate Embryo Clock: Common players dancing to a different beat

Vertebrate embryo somitogenesis is the earliest morphological manifestation of the characteristic patterned structure of the adult axial skeleton. Pairs of somites flanking the neural tube are formed periodically during early development, and the molecular mechanisms in temporal control of this early patterning event have been thoroughly studied. The discovery of a molecular Embryo Clock (EC) underlying the periodicity of somite formation shed light on the importance of gene expression dynamics for pattern formation. The EC is now known to be present in all vertebrate organisms studied and this mechanism was also described in limb development and stem cell differentiation. An outstanding question, however, remains unanswered: what sets the different EC paces observed in different organisms and tissues? This review aims to summarize the available knowledge regarding the pace of the EC, its regulation and experimental manipulation and to expose new questions that might help shed light on what is still to unveil.

Signaling oscillations in embryonic development

Tight spatiotemporal control of cellular behavior and cell fate decisions is paramount to the formation of multicellular organisms during embryonic development. Intercellular communication via signaling pathways mediates this control. Interestingly, these signaling pathways are not static, but dynamic and change in activity over time. Signaling oscillations as a specific type of dynamics are found in various signaling pathways and model systems. Functions of oscillations include the regulation of periodic events or the transmission of information by encoding signals in the dynamic properties of a signaling pathway. For instance, signaling oscillations in neural or pancreatic progenitor cells modulate their proliferation and differentiation. Oscillations between neighboring cells can also be synchronized, leading to the emergence of waves traveling through the tissue. Such population-wide signaling oscillations regulate for example the consecutive segmentation of vertebrate embryos, a process called somitogenesis. Here, we outline our current understanding of signaling oscillations in embryonic development, how signaling oscillations are generated, how they are studied and how they contribute to the regulation of embryonic development.

Geometric analysis of chondrogenic self-organisation of embryonic limb bud cells in micromass culture

Spatial and temporal control of chondrogenesis generates precise, species-specific patterns of skeletal structures in the developing vertebrate limb. The pattern-template is laid down when mesenchymal cells at the core of the early limb bud condense and undergo chondrogenic differentiation. Although the mechanisms involved in organising such complex patterns are not fully understood, the interplay between BMP and Wnt signalling pathways is fundamental. Primary embryonic limb bud cells grown under high-density micromass culture conditions spontaneously create a simple cartilage nodule pattern, presenting a model to investigate pattern generation. We describe a novel analytical approach to quantify geometric properties and spatial relationships between chondrogenic condensations, utilizing the micromass model. We follow the emergence of pattern in live cultures with nodules forming at regular distances, growing and changing shape over time. Gene expression profiling supports rapid chondrogenesis and transition to hypertrophy, mimicking the process of endochondral ossification within the limb bud. Manipulating the signalling environment through addition of BMP or Wnt ligands, as well as the BMP pathway antagonist Noggin, altered the differentiation profile and nodule pattern. BMP2 addition increased chondrogenesis while WNT3A or Noggin had the opposite effect, but with distinct pattern outcomes. Titrating these pro- and anti-chondrogenic factors and examining the resulting patterns support the hypothesis that regularly spaced cartilage nodules formed by primary limb bud cells in micromass culture are influenced by the balance of Wnt and BMP signalling under a Turing-like mechanism. This study demonstrates an approach for investigating the mechanisms governing chondrogenic spatial organization using simple micromass culture.

Spatial waves and temporal oscillations in vertebrate limb development

The mesenchymal tissue of the developing vertebrate limb bud is an excitable medium that sustains both spatial and temporal periodic phenomena. The first of these is the outcome of general Turing-type reaction-diffusion dynamics that generate spatial standing waves of cell condensations. These condensations are transformed into the nodules and rods of the cartilaginous, and eventually (in most species) the bony, endoskeleton. In the second, temporal periodicity results from intracellular regulatory dynamics that generate oscillations in the expression of one or more gene whose products modulate the spatial patterning system. Here we review experimental evidence from the chicken embryo, interpreted by a set of mathematical and computational models, that the spatial wave-forming system is based on two glycan-binding proteins, galectin-1A and galectin-8 in interaction with each other and the cells that produce them, and that the temporal oscillation occurs in the expression of the transcriptional coregulator Hes1. The multicellular synchronization of the Hes1 oscillation across the limb bud serves to coordinate the biochemical states of the mesenchymal cells globally, thereby refining and sharpening the spatial pattern. Significantly, the wave-forming reaction-diffusion-based mechanism itself, unlike most Turing-type systems, does not contain an oscillatory core, and may have evolved to this condition as it came to incorporate the cell-matrix adhesion module that enabled its pattern-forming capability.

Interplay of mesoscale physics and agent-like behaviors in the parallel evolution of aggregative multicellularity

Myxobacteria and dictyostelids are prokaryotic and eukaryotic multicellular lineages, respectively, that after nutrient depletion aggregate and develop into structures called fruiting bodies. The developmental processes and resulting morphological outcomes resemble one another to a remarkable extent despite their independent origins, the evolutionary distance between them and the lack of traceable homology in molecular mechanisms. We hypothesize that the morphological parallelism between the two lineages arises as the consequence of the interplay within multicellular aggregates between generic processes, physical and physicochemical processes operating similarly in living and non-living matter at the mesoscale (~10^-3-10^-1 m) and agent-like behaviors, unique to living systems and characteristic of the constituent cells, considered as autonomous entities acting according to internal rules in a shared environment. Here, we analyze the contributions of generic and agent-like determinants in myxobacteria and dictyostelid development and their roles in the generation of their common traits. Consequent to aggregation, collective cell-cell contacts mediate the emergence of liquid-like properties, making nascent multicellular masses subject to novel patterning and morphogenetic processes. In both lineages, this leads to behaviors such as streaming, rippling, and rounding-up, as seen in non-living fluids. Later the aggregates solidify, leading them to exhibit additional generic properties and motifs. Computational models suggest that the morphological phenotypes of the multicellular masses deviate from the predictions of generic physics due to the contribution of agent-like behaviors of cells such as directed migration, quiescence, and oscillatory signal transduction mediated by responses to external cues. These employ signaling mechanisms that reflect the evolutionary histories of the respective organisms. We propose that the similar developmental trajectories of myxobacteria and dictyostelids are more due to shared generic physical processes in coordination with analogous agent-type behaviors than to convergent evolution under parallel selection regimes. Insights from the biology of these aggregative forms may enable a unified understanding of developmental evolution, including that of animals and plants.

Mathematical models for cell migration: a non-local perspective

We provide a review of recent advancements in non-local continuous models for migration, mainly from the perspective of its involvement in embryonal development and cancer invasion. Particular emphasis is placed on spatial non-locality occurring in advection terms, used to characterize a cell's motility bias according to its interactions with other cellular and acellular components in its vicinity (e.g. cell-cell and cell-tissue adhesions, non-local chemotaxis), but we also briefly address spatially non-local source terms. Following a short introduction and description of applications, we give a systematic classification of available PDE models with respect to the type of featured non-localities and review some of the mathematical challenges arising from such models, with a focus on analytical aspects. This article is part of the theme issue 'Multi-scale analysis and modelling of collective migration in biological systems'.

Multiscale modeling of vertebrate limb development

We review the current state of mathematical modeling of cartilage pattern formation in vertebrate limbs. We place emphasis on several reaction-diffusion type models that have been proposed in the last few years. These models are grounded in more detailed knowledge of the relevant regulatory processes than previous ones but generally refer to different molecular aspects of these processes. Considering these models in light of comparative phylogenomics permits framing of hypotheses on the evolutionary order of appearance of the respective mechanisms and their roles in the fin-to-limb transition. This article is categorized under: Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Mechanistic Models Developmental Biology > Developmental Processes in Health and Disease Analytical and Computational Methods > Analytical Methods.

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

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

Why a Constant Number of Vertebrae? Digital Control of Segmental Identity during Vertebrate Development: The Somite Cycle Controls a Digital, Chromatin-Based Counter That Defines Segmental Identity and Body Plans in Vertebrate Animals

It is not understood how the numbers and identities of vertebrae are controlled during mammalian development. The remarkable robustness and conservation of segmental numbers may suggest the digital nature of the underlying process. The study proposes a mechanism that allows cells to obtain and store the segmental information in digital form, and to produce a pattern of chromatin accessibility that in turn regulates Hox gene expression specific to the metameric segment. The model requires that a regulatory element be present such that the number of occurrences of the motif between two consecutive Hox genes equals the number of segments under the control of the anterior gene. This is true for the recently discovered hydroxyl radical cleavage 3bp-periodic (HRC3) motif, associated with histone modifications and developmental genes. The finding not only allows the correct prediction of the numbers of segments using only sequence information, but also resolves the 40-year-old enigma of the function of temporal and spatial collinearity of Hox genes. The logic of the mechanism is illustrated in the attached animated video. How different aspects of the proposed mechanism can be tested experimentally is also discussed.

Inherency of Form and Function in Animal Development and Evolution

I discuss recent work on the origins of morphology and cell-type diversification in Metazoa – collectively the animals – and propose a scenario for how these two properties became integrated, with the help of a third set of processes, cellular pattern formation, into the developmental programs seen in present-day metazoans. Inherent propensities to generate familiar forms and cell types, in essence a parts kit for the animals, are exhibited by present-day organisms and were likely more prominent in primitive ones. The structural motifs of animal bodies and organs, e.g., multilayered, hollow, elongated and segmented tissues, internal and external appendages, branched tubes, and modular endoskeletons, can be accounted for by the properties of mesoscale masses of metazoan cells. These material properties, in turn, resulted from the recruitment of “generic” physical forces and mechanisms – adhesion, contraction, polarity, chemical oscillation, diffusion – by toolkit molecules that were partly conserved from unicellular holozoan antecedents and partly novel, distributed in the different metazoan phyla in a fashion correlated with morphological complexity. The specialized functions of the terminally differentiated cell types in animals, e.g., contraction, excitability, barrier function, detoxification, excretion, were already present in ancestral unicellular organisms. These functions were implemented in metazoan differentiation in some cases using the same transcription factors as in single-celled ancestors, although controlled by regulatory mechanisms that were hybrids between earlier-evolved processes and regulatory innovations, such as enhancers. Cellular pattern formation, mediated by released morphogens interacting with biochemically responsive and excitable tissues, drew on inherent self-organizing processes in proto-metazoans to transform clusters of holozoan cells into animal embryos and organs.