Search

International audience; Morphogenesis, wound healing and some cancer metastases depend upon migration of cell collectives that need to be guided to their destination as well as coordinated with other cell movements. During zebrafish gastrulation, extension of the embryonic axis is led by the mesendodermal polster that migrates towards the animal pole, followed by axial mesoderm that is undergoing convergence and extension. We here investigate how polster cells are guided towards the animal pole. Using a combination of precise laser ablations, advanced transplantations and functional as well as silico approaches, we establish that the directional information guiding polster cells is mechanical, and is provided by the anteriorward migration of the following cells. This information is detected by cell-cell contact through E-Cadherin/α-Catenin mechanotransduction and propagates from cell to cell over the whole tissue. Such guidance of migrating cells by followers ensures long-range coordination of movements and developmental robustness.

The mechanisms of organ size control remain poorly understood. A key question is how cells collectively sense the overall status of a tissue. We addressed this problem focusing on mouse liver regeneration. Using digital tissue reconstruction and quantitative image analysis, we found that the apical surface of hepatocytes forming the bile canalicular network expands concomitant with an increase in F‐actin and phospho‐myosin, to compensate an overload of bile acids. These changes are sensed by the Hippo transcriptional co‐activator YAP, which localizes to apical F‐actin‐rich regions and translocates to the nucleus in dependence of the integrity of the actin cytoskeleton. This mechanism tolerates moderate bile acid fluctuations under tissue homeostasis, but activates YAP in response to sustained bile acid overload. Using an integrated biophysical–biochemical model of bile pressure and Hippo signaling, we explained this behavior by the existence of a mechano‐sensory mechanism that activates YAP in a switch‐like manner. We propose that the apical surface of hepatocytes acts as a self‐regulatory mechano‐sensory system that responds to critical levels of bile acids as readout of tissue status., This study shows that the apical surface of hepatocytes acts as a self‐regulatory mechano‐sensory system that responds to critical levels of bile acids as readout of tissue status and activates YAP, by a mechanism dependent on the actin cytoskeleton.

Physiological liver cell replacement is central to maintaining the organ's high metabolic activity, although its characteristics are difficult to study in humans. Using retrospective 14 C birth dating of cells, we report that human hepatocytes show continuous and lifelong turnover, maintaining the liver a young organ (average age < 3 years ). Hepatocyte renewal is highly dependent on the ploidy level. Diploid hepatocytes show an seven-fold higher annual exchange rate than polyploid hepatocytes. These observations support the view that physiological liver cell renewal in humans is mainly dependent on diploid hepatocytes, whereas polyploid cells are compromised in their ability to divide. Moreover, cellular transitions between these two subpopulations are limited, with minimal contribution to the respective other ploidy class under homeostatic conditions. With these findings, we present a new integrated model of liver cell generation in humans that provides fundamental insights into liver cell turnover dynamics.

Simulations of tissue-specific effects of primary acute viral infections like COVID-19 are essential for understanding disease outcomes and optimizing therapies. Such simulations need to support continuous updating in response to rapid advances in understanding of infection mechanisms, and parallel development of components by multiple groups. We present an open-source platform for multiscale spatiotemporal simulation of an epithelial tissue, viral infection, cellular immune response and tissue damage, specifically designed to be modular and extensible to support continuous updating and parallel development. The base simulation of a simplified patch of epithelial tissue and immune response exhibits distinct patterns of infection dynamics from widespread infection, to recurrence, to clearance. Slower viral internalization and faster immune-cell recruitment slow infection and promote containment. Because antiviral drugs can have side effects and show reduced clinical effectiveness when given later during infection, we studied the effects on progression of treatment potency and time-of-first treatment after infection. In simulations, even a low potency therapy with a drug which reduces the replication rate of viral RNA greatly decreases the total tissue damage and virus burden when given near the beginning of infection. Many combinations of dosage and treatment time lead to stochastic outcomes, with some simulation replicas showing clearance or control (treatment success), while others show rapid infection of all epithelial cells (treatment failure). Thus, while a high potency therapy usually is less effective when given later, treatments at late times are occasionally effective. We illustrate how to extend the platform to model specific virus types (e.g., hepatitis C) and add additional cellular mechanisms (tissue recovery and variable cell susceptibility to infection), using our software modules and publicly-available software repository.

Summary: The Hedgehog (Hh) and Wnt/β-Catenin (Wnt) cascades are morphogen pathways whose pronounced influence on adult liver metabolism has been identified in recent years. How both pathways communicate and control liver metabolic functions are largely unknown. Detecting core components of Wnt and Hh signaling and mathematical modeling showed that both pathways in healthy liver act largely complementary to each other in the pericentral (Wnt) and the periportal zone (Hh) and communicate mainly by mutual repression. The Wnt/Hh module inversely controls the spatiotemporal operation of various liver metabolic pathways, as revealed by transcriptome, proteome, and metabolome analyses. Shifting the balance to Wnt (activation) or Hh (inhibition) causes pericentralization and periportalization of liver functions, respectively. Thus, homeostasis of the Wnt/Hh module is essential for maintaining proper liver metabolism and to avoid the development of certain metabolic diseases. With caution due to minor species-specific differences, these conclusions may hold for human liver as well. : Wnt/β-catenin and Hh signaling contribute to embryogenesis as well as to the maintenance of organ homeostasis through intensive crosstalk. Here, Kolbe et al. describe that both pathways act largely complementary to each other in the healthy liver and that this crosstalk is responsible for the maintenance of metabolic zonation.

SummaryBackgroundEndocytosis allows the import and distribution of cargo into a series of endosomes with distinct morphological and biochemical characteristics. Our current understanding of endocytic cargo trafficking is based on the kinetics of net cargo transport between endosomal compartments without considering individual endosomes. However, endosomes form a dynamic network of membranes undergoing fusion and fission, thereby continuously exchanging and redistributing cargo. The macroscopic kinetic properties, i.e., the properties of the endosomal network as a whole, result from the collective behaviors of many individual endosomes, a problem so far largely unaddressed.ResultsHere, we developed a general theoretical framework to describe the dynamics of cargo distributions in the endosomal network. We combined the theory with quantitative experiments to study how the macroscopic kinetic properties of the endosomal network emerge from microscopic processes at the level of individual endosomes. We compared our theory predictions to experimental data in which dynamic distributions of endocytosed low-density lipoprotein (LDL) were quantified.ConclusionsOur theory can quantitatively describe the observed cargo distributions as a function of time. Remarkably, the theory allows determining microscopic kinetic parameters such as the fusion rate between endosomes from still images of cargo distributions at different times of internalization. We show that this method is robust and sensitive because cargo distributions result from an average over many stochastic events in many cells. Our results provide theoretical and experimental support to the “funnel model” of endosome progression and suggest that the conversion of early to late endosomes is the major mode of LDL trafficking.

SummaryActivation of ErbB receptors by epidermal growth factor (EGF) or heregulin (HRG) determines distinct cell-fate decisions, although signals propagate through shared pathways. Using mathematical modeling and experimental approaches, we unravel how HRG and EGF generate distinct, all-or-none responses of the phosphorylated transcription factor c-Fos. In the cytosol, EGF induces transient and HRG induces sustained ERK activation. In the nucleus, however, ERK activity and c-fos mRNA expression are transient for both ligands. Knockdown of dual-specificity phosphatases extends HRG-stimulated nuclear ERK activation, but not c-fos mRNA expression, implying the existence of a HRG-induced repressor of c-fos transcription. Further experiments confirmed that this repressor is mainly induced by HRG, but not EGF, and requires new protein synthesis. We show how a spatially distributed, signaling-transcription cascade robustly discriminates between transient and sustained ERK activities at the c-Fos system level. The proposed control mechanisms are general and operate in different cell types, stimulated by various ligands.

Axolotls are unique in their ability to regenerate the spinal cord. However, the mechanisms that underlie this phenomenon remain poorly understood. Previously, we showed that regenerating stem cells in the axolotl spinal cord revert to a molecular state resembling embryonic neuroepithelial cells and functionally acquire rapid proliferative divisions (Rodrigo Albors et al., 2015). Here, we refine the analysis of cell proliferation in space and time and identify a high- proliferation zone in the regenerating spinal cord that shifts posteriorly over time. By tracking sparsely-labeled cells, we also quantify cell influx into the regenerate. Taking a mathematical modeling approach, we integrate these quantitative datasets of cell proliferation, neural stem cell activation and cell influx, to predict regenerative tissue outgrowth. Our model shows that while cell influx and neural stem cell activation play a minor role, the acceleration of the cell cycle is the major driver of regenerative spinal cord outgrowth in axolotls., Facultad de Ciencias Exactas, Instituto de Física de Líquidos y Sistemas Biológicos

The polarisation of cells and tissues is fundamental for tissue morphogenesis during biological development and regeneration. A deeper understanding of biological polarity pattern formation can be gained from the consideration of pattern reorganisation in response to an opposing instructive cue, which we here consider by example of experimentally inducible body axis inversions in planarian flatworms. Our dynamically diluted alignment model represents three processes: entrainment of cell polarity by a global signal, local cell-cell coupling aligning polarity among neighbours and cell turnover inserting initially unpolarised cells. We show that a persistent global orienting signal determines the final mean polarity orientation in this stochastic model. Combining numerical and analytical approaches, we find that neighbour coupling retards polarity pattern reorganisation, whereas cell turnover accelerates it. We derive a formula for an effective neighbour coupling strength integrating both effects and find that the time of polarity reorganisation depends linearly on this effective parameter and no abrupt transitions are observed. This allows to determine neighbour coupling strengths from experimental observations. Our model is related to a dynamic $8$-Potts model with annealed site-dilution and makes testable predictions regarding the polarisation of dynamic systems, such as the planarian epithelium., Preprint as prior to first submission to Journal of the Royal Society Interface. 25 pages, 6 figures, plus supplement (18 pages, contains 1 table and 7 figures). A supplementary movie is available from https://dx.doi.org/10.6084/m9.figshare.c3887818

Bile, the central metabolic product of the liver, is transported by the bile canaliculi network. The impairment of bile flow in cholestatic liver diseases has urged a demand for insights into its regulation. Here, we developed a predictive 3D multi-scale model that simulates fluid dynamic properties successively from the subcellular to the tissue level. The model integrates the structure of the bile canalicular network in the mouse liver lobule, as determined by high-resolution confocal and serial block-face scanning electron microscopy, with measurements of bile transport by intravital microscopy. The combined experiment-theory approach revealed spatial heterogeneities of biliary geometry and hepatocyte transport activity. Based on this, our model predicts gradients of bile velocity and pressure in the liver lobule. Validation of the model predictions by pharmacological inhibition of Rho kinase demonstrated a requirement of canaliculi contractility for bile flow in vivo. Our model can be applied to functionally characterize liver diseases and quantitatively estimate biliary transport upon drug-induced liver injury. Graphical Abstract

The muscle segments of fish have a folded shape, termed a chevron, which is thought to be optimal for the undulating body movements of swimming. However, the mechanism shaping the chevron during embryogenesis is not understood. Here, we use time-lapse microscopy of developing zebrafish embryos spanning the entire somitogenesis period to quantitate the dynamics of chevron shape development. Comparing such time courses with the start of movements in wildtype zebrafish and analyzing immobile mutants, we show that the previously implicated body movements do not play a role in chevron formation. Further, the monotonic increase of chevron angle along the anteroposterior axis revealed by our data constrains or rules out possible contributions by previously proposed mechanisms. In particular, we find that muscle pioneers are not required for chevron formation. We put forward a tension-and-resistance mechanism involving interactions between intra-segmental tension and segment boundaries. To evaluate this mechanism, we derive and analyze a mechanical model of a chain of contractile and resisting elements. The predictions of this model are verified by comparison to experimental data. Altogether, our results support the notion that a simple physical mechanism suffices to self-organize the observed spatiotemporal pattern in chevron formation.

In many animals, regenerative processes can replace lost body parts. Organ and tissue regeneration consequently also hold great medical promise. The regulation of regenerative processes is achieved through concerted actions of multiple organizational levels of the organism, from diffusing molecules and cellular gene expression patterns up to tissue mechanics. Our intuition is usually not adapted well to this degree of complexity and the quantitative aspects of the regulation of regenerative processes remain poorly understood. One way out of this dilemma lies in the combination of experimentation and mathematical modeling within an iterative process of model development/refinement, model predictions for novel experimental conditions, quantitative experiments testing these predictions, and subsequent model refinement. This interdisciplinary approach has already provided key insights into smaller scale processes during embryonic development and a so-far limited number of more complex regeneration processes. This review discusses selected theoretical and interdisciplinary studies and is structured along the three phases of regeneration: (1) initiation of a regeneration response, (2) tissue patterning during regenerate growth, (3) arresting the regeneration response. Moreover, we highlight the opportunities provided by extensions of mathematical models from developmental processes toward the study of related regenerative processes.

International audience; The Hedgehog (Hh) signalling cascade is highly conserved and involved in development and disease throughout evolution. Nevertheless, in comparison with other pathways, our mechanistic understanding of Hh signal transduction is remarkably incomplete. In the absence of ligand, the Hh receptor Patched (Ptc) represses the key signal transducer Smoothened (Smo) through an unknown mechanism. Hh binding to Ptc alleviates this repression, causing Smo redistribution to the plasma membrane, phosphorylation and opening of the Smo cytoplasmic tail, and Smo oligomerisation. However, the order and interdependence of these events is as yet poorly understood. We have mathematically modelled and simulated Smo activation for two alternative modes of pathway activation, with Ptc primarily affecting either Smo localisation or phosphorylation. Visualising Smo activation through a novel, fluorescence-based reporter allowed us to test these competing models. Here, we show that Smo localisation to the plasma membrane is sufficient for phosphorylation of the cytoplasmic tail in the presence of Ptc. Using fluorescence cross-correlation spectroscopy (FCCS), we also demonstrate that inactivation of Ptc by Hh induces Smo clustering irrespective of Smo phosphorylation. Our observations therefore support a model of Hh signal transduction whereby Smo subcellular localisation and not phosphorylation is the primary target of Ptc function.

The cell fate decision of multi-potent pancreatic progenitor cells between the exocrine and endocrine lineages is regulated by Notch signalling, mediated by cell–cell interactions. However, canonical models of Notch-mediated lateral inhibition cannot explain the scattered spatial distribution of endocrine cells and the cell-type ratio in the developing pancreas. Based on evidence from acinar-to-islet cell transdifferentiationin vitro, we propose that lateral stabilization, i.e. positive feedback between adjacent progenitor cells, acts in parallel with lateral inhibition to regulate pattern formation in the pancreas. A simple mathematical model of transcriptional regulation and cell–cell interaction reveals the existence of multi-stability of spatial patterns whose simultaneous occurrence causes scattering of endocrine cells in the presence of noise. The scattering pattern allows for control of the endocrine-to-exocrine cell-type ratio by modulation of lateral stabilization strength. These theoretical results suggest a previously unrecognized role for lateral stabilization in lineage specification, spatial patterning and cell-type ratio control in organ development.

International audience; Many internal epithelial organs derive from cysts, which are tissues comprised of bent epithelial cell layers enclosing a lumen. Ion accumulation in the lumen drives water influx and consequently water accumulation and cyst expansion. Lumen-size recognition is important for the regulation of organ size. When lumen size and cyst size are not controlled, diseases can result; for instance, renal failure of the kidney. We develop a mechanistic mathematical model of lumen expansion in order to investigate the mechanisms for saturation of cyst growth. We include fluid accumulation in the lumen, osmotic and elastic pressure, ion transport and stretch-induced cell division. We find that the lumen volume increases in two phases: first, due to fluid accumulation stretching the cells, then in the second phase, the volume increase follows the increase in cell number until proliferation ceases as stretch forces relax. The model is quantitatively fitted to published data of in vitro cyst growth and predicts steady state lumen size as a function of the model parameters.

Key cellular functions and developmental processes rely on cascades of GTPases. GTPases of the Rab family provide a molecular ID code to the generation, maintenance and transport of intracellular compartments. Here, we addressed the molecular design principles of endocytosis by focusing on the conversion of early endosomes into late endosomes, which entails replacement of Rab5 by Rab7. We modelled this process as a cascade of functional modules of interacting Rab GTPases. We demonstrate that intermodule interactions share similarities with the toggle switch described for the cell cycle. However, Rab5-to-Rab7 conversion is rather based on a newly characterized ‘cut-out switch' analogous to an electrical safety-breaker. Both designs require cooperativity of auto-activation loops when coupled to a large pool of cytoplasmic proteins. Live cell imaging and endosome tracking provide experimental support to the cut-out switch in cargo progression and conversion of endosome identity along the degradative pathway. We propose that, by reconciling module performance with progression of activity, the cut-out switch design could underlie the integration of modules in regulatory cascades from a broad range of biological processes.

We evaluate different mechanisms for spatial domain formation of guanosine triphosphatases (GTPases) on cellular membranes. A kinetic model of the basic guanine-nucleotide cycle common to all GTPases is developed and coupled along a one-dimensional axis by diffusion of inactive and activated GTPases. We ask whether a parameter set exists such that domain formation is possible by Turing’s mechanism, i.e., purely by reactions and diffusion, and show that the Turing instability does not occur in this model for any parameter combination. But, as revealed by stability and bifurcation analysis, domain formation is reproduced after augmenting the model with combinations of two spatial interaction mechanisms: 1. attraction and 2. adhesion among active GTPases. These interactions can be mediated by effector proteins that bind active GTPases, and the model therefore predicts domains to disintegrate if effector binding is inhibited.

The modulated amplitude waves (MAWs) in the complex Ginzburg-Landau equation (CGLE) were systematically studied. These periodic coherent structures were found to originate in supercritical bifurcations due to the BF instability of the CGLE. Indications were given in favor of the existence of the phase turbulent regime even in the thermodynamic limit., EU under Contract No. ERBFMBICT972554.

Cell fate reprogramming, such as the generation of insulin-producing β cells from other pancreas cells, can be achieved by external modulation of key transcription factors. However, the known gene regulatory interactions that form a complex network with multiple feedback loops make it increasingly difficult to design the cell reprogramming scheme because the linear regulatory pathways as schemes of causal influences upon cell lineages are inadequate for predicting the effect of transcriptional perturbation. However, sufficient information on regulatory networks is usually not available for detailed formal models. Here we demonstrate that by using the qualitatively described regulatory interactions as the basis for a coarse-grained dynamical ODE (ordinary differential equation) based model, it is possible to recapitulate the observed attractors of the exocrine and β, δ, α endocrine cells and to predict which gene perturbation can result in desired lineage reprogramming. Our model indicates that the constraints imposed by the incompletely elucidated regulatory network architecture suffice to build a predictive model for making informed decisions in choosing the set of transcription factors that need to be modulated for fate reprogramming.