A computational framework for the morpho-elastic development of molluskan shells by surface and volume growth

A multiscale computational framework for the development of spines in molluscan shells

From mathematical models of growth to computer simulations of pigmentation, the study of shell formation has given rise to an abundant number of models, working at various scales. Yet, attempts to combine those models have remained sparse, due to the challenge of combining categorically different approaches. In this paper, we propose a framework to streamline the process of combining the molecular and tissue scales of shell formation. We choose these levels as a proxy to link the genotype level, which is better described by molecular models, and the phenotype level, which is better described by tissue-level mechanics. We also show how to connect observations on shell populations to the approach, resulting in collections of molecular parameters that may be associated with different populations of real shell specimens. The approach is as follows: we use a Quality-Diversity algorithm, a type of black-box optimization algorithm, to explore the range of concentration profiles emerging as solutions of a molecular model, and that define growth patterns for the mechanical model. At the same time, the mechanical model is simulated over a wide range of growth patterns, resulting in a variety of spine shapes. While time-consuming, these steps only need to be performed once and then function as look-up tables. Actual pictures of shell spines can then be matched against the list of existing spine shapes, yielding a potential growth pattern which, in turn, gives us matching molecular parameters. The framework is modular, such that models can be easily swapped without changing the overall working of the method. As a demonstration of the approach, we solve specific molecular and mechanical models, adapted from available theoretical studies on molluscan shells, and apply the multiscale framework to evaluate the characteristics of spines from three distinct populations of Turbo sazae.

Physical confinement and cell proximity increase cell migration rates and invasiveness: A mathematical model of cancer cell invasion through flexible channels

Cancer cell migration between different body parts is the driving force behind cancer metastasis, which is the main cause of mortality of patients. Migration of cancer cells often proceeds by penetration through narrow cavities in locally stiff, yet flexible tissues. In our previous work, we developed a model for cell geometry evolution during invasion, which we extend here to investigate whether leader and follower (cancer) cells that only interact mechanically can benefit from sequential transmigration through narrow micro-channels and cavities. We consider two cases of cells sequentially migrating through a flexible channel: leader and follower cells being closely adjacent or distant. Using Wilcoxon's signed-rank test on the data collected from Monte Carlo simulations, we conclude that the modelled transmigration speed for the follower cell is significantly larger than for the leader cell when cells are distant, i.e. follower cells transmigrate after the leader has completed the crossing. Furthermore, it appears that there exists an optimum with respect to the width of the channel such that cell moves fastest. On the other hand, in the case of closely adjacent cells, effectively performing collective migration, the leader cell moves 12% faster since the follower cell pushes it. This work shows that mechanical interactions between cells can increase the net transmigration speed of cancer cells, resulting in increased invasiveness. In other words, interaction between cancer cells can accelerate metastatic invasion.

A Numerical Investigation of Dimensionless Numbers Characterizing Meltpool Morphology of the Laser Powder Bed Fusion Process

Microstructure evolution in metal additive manufacturing (AM) is a complex multi-physics and multi-scale problem. Understanding the impact of AM process conditions on the microstructure evolution and the resulting mechanical properties of the printed component remains an active area of research. At the meltpool scale, the thermo-fluidic governing equations have been extensively modeled in the literature to understand the meltpool conditions and the thermal gradients in its vicinity. In many phenomena governed by partial differential equations, dimensional analysis and identification of important dimensionless numbers can provide significant insights into the process dynamics. In this context, we present a novel strategy using dimensional analysis and the linear least-squares regression method to numerically investigate the thermo-fluidic governing equations of the Laser Powder Bed Fusion AM process. First, the governing equations are solved using the Finite Element Method, and the model predictions are validated by comparing with experimentally estimated cooling rates, and with numerical results from the literature. Then, through dimensional analysis, an important dimensionless quantity interpreted as a measure of heat absorbed by the powdered material and the meltpool, is identified. This dimensionless measure of absorbed heat, along with classical dimensionless quantities such as Péclet, Marangoni, and Stefan numbers, are employed to investigate advective transport in the meltpool for different alloys. Further, the framework is used to study variations in the thermal gradients and the solidification cooling rate. Important correlations linking meltpool morphology and microstructure-evolution-related variables with classical dimensionless numbers are the key contribution of this work.

Comparison between a phenomenological approach and a morphoelasticity approach regarding the displacement of extracellular matrix

Plastic (permanent) deformations were earlier, modeled by a phenomenological model in Peng and Vermolen (Biomech Model Mechanobiol 19(6):2525–2551, 2020). In this manusctipt, we consider a more physics-based formulation that is based on morphoelasticity. We firstly introduce the morphoelasticity approach and investigate the impact of various input variables on the output parameters by sensitivity analysis. A comparison of both model formulations shows that both models give similar computational results. Furthermore, we carry out Monte Carlo simulations of the skin contraction model containing the morphoelasticity approach. Most statistical correlations from the two models are similar, however, the impact of the collagen density on the severeness of contraction is larger for the morphoelasticity model than for the phenomenological model.

The physical basis of mollusk shell chiral coiling

Snails are model organisms for studying the genetic, molecular, and developmental bases of left-right asymmetry in Bilateria. However, the development of their typical helicospiral shell, present for the last 540 million years in environments as different as the abyss or our gardens, remains poorly understood. Conversely, ammonites typically have a bilaterally symmetric, planispiraly coiled shell, with only 1% of 3,000 genera displaying either a helicospiral or a meandering asymmetric shell. A comparative analysis suggests that the development of chiral shells in these mollusks is different and that, unlike snails, ammonites with asymmetric shells probably had a bilaterally symmetric body diagnostic of cephalopods. We propose a mathematical model for the growth of shells, taking into account the physical interaction during development between the soft mollusk body and its hard shell. Our model shows that a growth mismatch between the secreted shell tube and a bilaterally symmetric body in ammonites can generate mechanical forces that are balanced by a twist of the body, breaking shell symmetry. In gastropods, where a twist is intrinsic to the body, the same model predicts that helicospiral shells are the most likely shell forms. Our model explains a large diversity of forms and shows that, although molluscan shells are incrementally secreted at their opening, the path followed by the shell edge and the resulting form are partly governed by the mechanics of the body inside the shell, a perspective that explains many aspects of their development and evolution.

A formalism for modelling traction forces and cell shape evolution during cell migration in various biomedical processes

The phenomenological model for cell shape deformation and cell migration Chen (BMM 17:1429–1450, 2018), Vermolen and Gefen (BMM 12:301–323, 2012), is extended with the incorporation of cell traction forces and the evolution of cell equilibrium shapes as a result of cell differentiation. Plastic deformations of the extracellular matrix are modelled using morphoelasticity theory. The resulting partial differential differential equations are solved by the use of the finite element method. The paper treats various biological scenarios that entail cell migration and cell shape evolution. The experimental observations in Mak et al. (LC 13:340–348, 2013), where transmigration of cancer cells through narrow apertures is studied, are reproduced using a Monte Carlo framework.

Learning from beautiful monsters: phylogenetic and morphogenetic implications of left-right asymmetry in ammonoid shells

BACKGROUND

Many pathologies that modify the shell geometry and ornamentation of ammonoids are known from the fossil record. Since they may reflect the developmental response of the organism to a perturbation (usually a sublethal injury), their study is essential for exploring the developmental mechanisms of these extinct animals. Ammonoid pathologies are also useful to assess the value of some morphological characters used in taxonomy, as well as to improve phylogenetic reconstructions and evolutionary scenarios.

RESULTS

We report on the discovery of an enigmatic pathological middle Toarcian (Lower Jurassic) ammonoid specimen from southern France, characterized by a pronounced left-right asymmetry in both ornamentation and suture lines. For each side independently, the taxonomic interpretations of ornamentation and suture lines are congruent, suggesting a Hildoceras semipolitum species assignment for the left side and a Brodieia primaria species assignment for the right side. The former exhibits a lateral groove whereas the second displays sinuous ribs. This specimen, together with the few analogous cases reported in the literature, lead us to erect a new forma-type pathology herein called "forma janusa" for specimens displaying a left-right asymmetry in the absence of any clear evidence of injury or parasitism, whereby the two sides match with the regular morphology of two distinct, known species.

CONCLUSIONS

Since "forma janusa" specimens reflect the underlying developmental plasticity of the ammonoid taxa, we hypothesize that such specimens may also indicate unsuspected phylogenetic closeness between the two displayed taxa and may even reveal a direct ancestor-descendant relationship. This hypothesis is not, as yet, contradicted by the stratigraphical data at hand: in all studied cases the two distinct taxa correspond to contemporaneous or sub-contemporaneous taxa. More generally, the newly described specimen suggests that a hitherto unidentified developmental link may exist between sinuous ribs and lateral grooves. Overall, we recommend an integrative approach for revisiting aberrant individuals that illustrate the intricate links among shell morphogenesis, developmental plasticity and phylogeny.