Search

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. [ABSTRACT FROM AUTHOR]

Animal color patterns are widely studied in ecology, evolution, and through mathematical modeling. Patterns may vary among distinct body parts such as the head, trunk or tail. As large amounts of photographic data is becoming more easily available, there is a growing need for general quantitative methods for capturing and analyzing the full complexity and details of pattern variation. Detailed information on variation in color pattern elements is necessary to understand how patterns are produced and established during development, and which evolutionary forces may constrain such a variation. Here, we develop an approach to capture and analyze variation in melanistic color pattern elements in leopard geckos. We use this data to study the variation among different body parts of leopard geckos and to draw inferences about their development. We compare patterns using 14 different indices such as the ratio of melanistic versus total area, the ellipticity of spots, and the size of spots and use these to define a composite distance between two patterns. Pattern presence/absence among the different body parts indicates a clear pathway of pattern establishment from the head to the back legs. Together with weak within-individual correlation between leg patterns and main body patterns, this suggests that pattern establishment in the head and tail may be independent from the rest of the body. We found that patterns vary greatest in size and density of the spots among body parts and individuals, but little in their average shapes. We also found a correlation between the melanistic patterns of the two front legs, as well as the two back legs, and also between the head, tail and trunk, especially for the density and size of the spots, but not their shape or inter-spot distance. Our data collection and analysis approach can be applied to other organisms to study variation in color patterns between body parts and to address questions on pattern formation and establishment in animals. [ABSTRACT FROM AUTHOR]

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 MethodsModels of Systems Properties and Processes > Mechanistic ModelsDevelopmental Biology > Developmental Processes in Health and DiseaseAnalytical and Computational Methods > Analytical Methods [ABSTRACT FROM AUTHOR]

Genotypic variation, environmental variation, and their interaction may produce variation in the developmental process and cause phenotypic differences among individuals. Developmental noise, which arises during development from stochasticity in cellular and molecular processes when genotype and environment are fixed, also contributes to phenotypic variation. While evolutionary biology has long focused on teasing apart the relative contribution of genes and environment to phenotypic variation, our understanding of the role of developmental noise has lagged due to technical difficulties in directly measuring the contribution of developmental noise. The influence of developmental noise is likely underestimated in studies of phenotypic variation due to intrinsic mechanisms within organisms that stabilize phenotypes and decrease variation. Since we are just beginning to appreciate the extent to which phenotypic variation due to stochasticity is potentially adaptive, the contribution of developmental noise to phenotypic variation must be separated and measured to fully understand its role in evolution. Here, we show that variation in the component of the developmental process corresponding to environmental and genetic factors (here treated together as a unit called the LALI-type) versus the contribution of developmental noise, can be distinguished for leopard gecko (Eublepharis macularius) head color patterns using mathematical simulations that model the role of random variation (corresponding to developmental noise) in patterning. Specifically, we modified the parameters of simulations corresponding to variation in the LALI-type to generate the full range of phenotypic variation in color pattern seen on the heads of eight leopard geckos. We observed that over the range of these parameters, variation in color pattern due to LALI-type variation exceeds that due to developmental noise in the studied gecko cohort. However, the effect of developmental noise on patterning is also substantial. Our approach addresses one of the major goals of evolutionary biology: to quantify the role of stochasticity in shaping phenotypic variation. [ABSTRACT FROM AUTHOR]

The tetrapod appendicular skeleton is initiated as spatially patterned mesenchymal condensations. The size and spacing of these condensations in avian limb buds are mediated by a reaction-diffusion-adhesion network consisting of galectins Gal-1A, Gal-8 and their cell surface receptors. In cell cultures, the appearance of condensations is synchronized across distances greater than the characteristic wavelength of their spatial pattern. We explored the possible role of observed oscillations of the transcriptional co-regulator Hes1 in this phenomenon. Treatment of micromass cultures with DAPT, a γ-secretase inhibitor, damped Hes1 oscillations, elevated Gal-1A and -8 mRNA levels, and led to irregularly-sized proto-condensations that subsequently fused. In developing limb buds, DAPT led to spatially non-uniform Hes1 expression and fused, truncated and misshapen digits. Periodicity in adhesive response to Gal-1A, a plausible Hes1-dependent function, was added to a previously tested mathematical model for condensation patterning by the two-galectin network. The enhanced model predicted regularization of patterning due to synchronization of Hes1 oscillations and resulting spatiotemporal coordination of its expression. The model also predicted changes in galectin expression and patterning in response to suppression of Hes1 expression, which were confirmed in in vitro experiments. Our results indicate that the two-galectin patterning network is regulated by Hes1 dynamics, the synchronization of which refines and regularizes limb skeletogenesis. • hes1 oscillates during limb precartilage condensation formation. • DAPT treatment damps hes1 expression and alters condensation and digital patterning. • Modeling suggests synchronized hes1 oscillations coordinate condensation formation. [ABSTRACT FROM AUTHOR]

The paired appendages (fins or limbs) of jawed vertebrates contain an endoskeleton consisting of nodules, bars and, in some groups, plates of cartilage, or bone arising from replacement of cartilaginous templates. The generation of the endoskeletal elements occurs by processes involving production and diffusion of morphogens, with, variously, positive and negative feedback circuits, adhesion, and receptor dynamics with similarities to the mechanism for chemical pattern formation proposed by Alan Turing. This review presents a unified interpretation of the evolution and functioning of these mechanisms. Studies are described indicating that protocondensations, compacted mesenchymal cell aggregates that prefigure the appendicular skeleton, arise through the adhesive activity of galectin-1, a matricellular protein with skeletogenic homologs in all jawed vertebrates. In the cartilaginous and lobe-finned fishes (and to a variable extent in ray-finned fishes) it additionally cooperates with an isoform of galectin-8 to constitute a self-organizing network capable of generating arrays of preskeletal nodules, bars and plates. Further, in the tetrapods, a putative galectin-8 control module was acquired that may have enabled proximodistal increase in the number of protocondensations. In parallel to this, other self-organizing networks emerged that acted, via Bmp, Wnt, Sox9 and Runx2, as well as transforming factor-β and fibronectin, to convert protocondensations into skeletal tissues. The progressive appearance and integration of these skeletogenic networks over evolution occurred in the context of an independently evolved system of Hox protein and Shh gradients that interfaced with them to tune the spatial wavelengths and refine the identities of the resulting arrays of elements. [ABSTRACT FROM AUTHOR]

We investigate the stability of Turing-type patterns in one spatial dimension in a system of reaction-diffusion equations with a term depending linearly on the spatial position. The system is a generic model of two interacting chemical species where production rates are dependent on a linear external gradient. This is motivated by mathematical models in developmental biology. In a previous paper, we found analytic approximations of Turing-like steady state patterns. In the present article, we derive conditions for the stability of these patterns and show bifurcation diagrams in two small parameters related to the slope of the external gradient and the ratio of the diffusion coefficients. [ABSTRACT FROM AUTHOR]

Background: A multiscale network of two galectins Galectin-1 (Gal-1) and Galectin-8 (Gal-8) patterns the avian limb skeleton. Among vertebrates with paired appendages, chondrichthyan fins typically have one or more cartilage plates and many repeating parallel endoskeletal elements, actinopterygian fins have more varied patterns of nodules, bars and plates, while tetrapod limbs exhibit tandem arrays of few, proximodistally increasing numbers of elements. We applied a comparative genomic and protein evolution approach to understand the origin of the galectin patterning network. Having previously observed a phylogenetic constraint on Gal-1 structure across vertebrates, we asked whether evolutionary changes of Gal-8 could have critically contributed to the origin of the tetrapod pattern. Results: Translocations, duplications, and losses of Gal-8 genes in Actinopterygii established them in different genomic locations from those that the Sarcopterygii (including the tetrapods) share with chondrichthyans. The sarcopterygian Gal-8 genes acquired a potentially regulatory non-coding motif and underwent purifying selection. The actinopterygian Gal-8 genes, in contrast, did not acquire the non-coding motif and underwent positive selection. Conclusion: These observations interpreted through the lens of a reaction-diffusion-adhesion model based on avian experimental findings can account for the distinct endoskeletal patterns of cartilaginous, ray-finned, and lobe-finned fishes, and the stereotypical limb skeletons of tetrapods. [ABSTRACT FROM AUTHOR]

We consider the geometric optics problem of finding a system of two reflectors that transform a spherical wavefront into a beam of parallel rays with prescribed intensity distribution. Using techniques from optimal transportation theory, it has been shown previously that this problem is equivalent to an infinite-dimensional linear programming (LP) problem. Here we investigate techniques for constructing the two reflectors numerically by considering the finite-dimensional LP problems which arise as approximations to the infinite-dimensional problem. A straightforward discretization has the disadvantage that the number of constraints increases rapidly with the mesh size, so only very coarse meshes are practical. To address this well-known issue we propose an iterative solution scheme. In each step, an LP problem is solved, where information from the previous iteration step is used to reduce the number of necessary constraints. As an illustration, we apply our proposed scheme to solve a problem with synthetic data, demonstrating that the method allows for much finer meshes than a simple discretization. We also give evidence that the scheme converges. There exists a growing literature for the application of optimal transportation theory to other beam shaping problems, and our proposed scheme is easy to adapt for these problems as well. [ABSTRACT FROM AUTHOR]

We consider three-dimensional Turing patterns and their isoconcentration surfaces corresponding to the equilibrium concentration of the reaction kinetics. We call these surfaces equilibrium concentration surfaces (EC surfaces). They are the interfaces between the regions of "high" and "low" concentrations in Turing patterns. We give alternate characterizations of EC surfaces by means of two variational principles, one of them being that they are optimal for diffusive transport. Several examples of EC surfaces are considered. Remarkably, they are often very well approximated by certain minimal surfaces. We give a dynamical explanation for the emergence of Scherk's surface in certain cases, a structure that has been observed numerically previously in [De Wit et al., 1997]. [ABSTRACT FROM AUTHOR]

In this article the authors discuss their software product CompuCell3D. The software is designed to allow researchers in biology and physics to set up various environments where experiments can be run and the conditions can be manipulated. For convenience the software has been written in a format with greater flexibility than other problem solving environment software. The authors discuss why the product is an advancement and also what they hope to upgrade in future versions of the software.

An abstract of the article "Multiscale Models for Vertebrate Limb Development," by Stuart A. Newman and colleagues is presented.

Ovarian cancer is amongst the most morbid of gynecological malignancies due to its diagnosis at an advanced stage, a transcoelomic mode of metastasis, and rapid transition to chemotherapeutic resistance. Like all other malignancies, the progression of ovarian cancer may be interpreted as an emergent outcome of the conflict between metastasizing cancer cells and the natural defense mounted by microenvironmental barriers to such migration. Here, we asked whether senescence in coelom-lining mesothelia, brought about by drug exposure, affects their interaction with disseminated ovarian cancer cells. We observed that cancer cells adhered faster on senescent human and murine mesothelial monolayers than on non-senescent controls. Time-lapse epifluorescence microscopy showed that mesothelial cells were cleared by a host of cancer cells that surrounded the former, even under sub-confluent conditions. A multiscale computational model predicted that such colocalized mesothelial clearance under sub-confluence requires greater adhesion between cancer cells and senescent mesothelia. Consistent with the prediction, we observed that senescent mesothelia expressed an extracellular matrix with higher levels of fibronectin, laminins and hyaluronan than non-senescent controls. On senescent matrix, cancer cells adhered more efficiently, spread better, and moved faster and persistently, aiding the spread of cancer. Inhibition assays using RGD cyclopeptides suggested the adhesion was predominantly contributed by fibronectin and laminin. These findings led us to propose that the senescence-associated matrisomal phenotype of peritoneal barriers enhances the colonization of invading ovarian cancer cells contributing to the metastatic burden associated with the disease., (© 2023. The Author(s).)

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., (© 2020 Wiley Periodicals, Inc.)

Dynamical systems in which geometrically extended model cells produce and interact with diffusible (morphogen) and nondiffusible (extracellular matrix) chemical fields have proved very useful as models for developmental processes. The embryonic vertebrate limb is an apt system for such mathematical and computational modeling since it has been the subject of hundreds of experimental studies, and its normal and variant morphologies and spatiotemporal organization of expressed genes are well known. Because of its stereotypical proximodistally generated increase in the number of parallel skeletal elements, the limb lends itself to being modeled by Turing-type systems which are capable of producing periodic, or quasiperiodic, arrangements of spot- and stripe-like elements. This chapter describes several such models, including, (i) a system of partial differential equations in which changing cell density enters into the dynamics explicitly, (ii) a model for morphogen dynamics alone, derived from the latter system in the "morphostatic limit" where cell movement relaxes on a much slower time-scale than cell differentiation, (iii) a discrete stochastic model for the simplified pattern formation that occurs when limb cells are placed in planar culture, and (iv) several hybrid models in which continuum morphogen systems interact with cells represented as energy-minimizing mesoscopic entities. Progress in devising computational methods for handling 3D, multiscale, multimodel simulations of organogenesis is discussed, as well as for simulating reaction-diffusion dynamics in domains of irregular shape.

The cellular Potts model (CPM) has been used for simulating various biological phenomena such as differential adhesion, fruiting body formation of the slime mold Dictyostelium discoideum, angiogenesis, cancer invasion, chondrogenesis in embryonic vertebrate limbs, and many others. We derive a continuous limit of a discrete one-dimensional CPM with the chemotactic interactions between cells in the form of a Fokker-Planck equation for the evolution of the cell probability density function. This equation is then reduced to the classical macroscopic Keller-Segel model. In particular, all coefficients of the Keller-Segel model are obtained from parameters of the CPM. Theoretical results are verified numerically by comparing Monte Carlo simulations for the CPM with numerics for the Keller-Segel model.

We present COMPUCELL3D, a software framework for three-dimensional simulation of morphogenesis in different organisms. COMPUCELL3D employs biologically relevant models for cell clustering, growth, and interaction with chemical fields. COMPUCELL3D uses design patterns for speed, efficient memory management, extensibility, and flexibility to allow an almost unlimited variety of simulations. We have verified COMPUCELL3D by building a model of growth and skeletal pattern formation in the avian (chicken) limb bud. Binaries and source code are available, along with documentation and input files for sample simulations, at http:// compucell.sourceforge.net.

We describe a 'reactor-diffusion' mechanism for precartilage condensation based on recent experiments on chondrogenesis in the early vertebrate limb and additional hypotheses. Cellular differentiation of mesenchymal cells into subtypes with different fibroblast growth factor (FGF) receptors occurs in the presence of spatio-temporal variations of FGFs and transforming growth factor-betas (TGF-betas). One class of differentiated cells produces elevated quantities of the extracellular matrix protein fibronectin, which initiates adhesion-mediated preskeletal mesenchymal condensation. The same class of cells also produces an FGF-dependent laterally acting inhibitor that keeps condensations from expanding beyond a critical size. We show that this 'reactor-diffusion' mechanism leads naturally to patterning consistent with skeletal form, and describe simulations of spatio-temporal distribution of these differentiated cell types and the TGF-beta and inhibitor concentrations in the developing limb bud.