Search

Early disease diagnosis is key to the effective treatment of diseases. Histopathological analysis of human biopsies is the gold standard to diagnose tissue alterations. However, this approach has low resolution and overlooks 3D (three-dimensional) structural changes resulting from functional alterations. Here, we applied multiphoton imaging, 3D digital reconstructions and computational simulations to generate spatially resolved geometrical and functional models of human liver tissue at different stages of non-alcoholic fatty liver disease (NAFLD). We identified a set of morphometric cellular and tissue parameters correlated with disease progression, and discover profound topological defects in the 3D bile canalicular (BC) network. Personalized biliary fluid dynamic simulations predicted an increased pericentral biliary pressure and micro-cholestasis, consistent with elevated cholestatic biomarkers in patients’ sera. Our spatially resolved models of human liver tissue can contribute to high-definition medicine by identifying quantitative multiparametric cellular and tissue signatures to define disease progression and provide new insights into NAFLD pathophysiology. A combination of high-resolution imaging and modeling approaches facilitates the study of the mechanisms and clinical progression of non-alcoholic fatty liver disease in humans.

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.

Research software has become a central asset in academic research. It optimizes existing and enables new research methods, implements and embeds research knowledge, and constitutes an essential research product in itself. Research software must be sustainable in order to understand, replicate, reproduce, and build upon existing research or conduct new research effectively. In other words, software must be available, discoverable, usable, and adaptable to new needs, both now and in the future. Research software therefore requires an environment that supports sustainability. Hence, a change is needed in the way research software development and maintenance are currently motivated, incentivized, funded, structurally and infrastructurally supported, and legally treated. Failing to do so will threaten the quality and validity of research. In this paper, we identify challenges for research software sustainability in Germany and beyond, in terms of motivation, selection, research software engineering personnel, funding, infrastructure, and legal aspects. Besides researchers, we specifically address political and academic decision-makers to increase awareness of the importance and needs of sustainable research software practices. In particular, we recommend strategies and measures to create an environment for sustainable research software, with the ultimate goal to ensure that software-driven research is valid, reproducible and sustainable, and that software is recognized as a first class citizen in research. This paper is the outcome of two workshops run in Germany in 2019, at deRSE19 - the first International Conference of Research Software Engineers in Germany - and a dedicated DFG-supported follow-up workshop in Berlin.

How epithelial cells coordinate their polarity to form functional tissues is an open question in cell biology. Here, we characterize a unique type of polarity found in liver tissue, nematic cell polarity, which is different from vectorial cell polarity in simple, sheet-like epithelia. We propose a conceptual and algorithmic framework to characterize complex patterns of polarity proteins on the surface of a cell in terms of a multipole expansion. To rigorously quantify previously observed tissue-level patterns of nematic cell polarity (Morales-Navarrete et al., eLife 2019), we introduce the concept of co-orientational order parameters, which generalize the known biaxial order parameters of the theory of liquid crystals. Applying these concepts to three-dimensional reconstructions of single cells from high-resolution imaging data of mouse liver tissue, we show that the axes of nematic cell polarity of hepatocytes exhibit local coordination and are aligned with the biaxially anisotropic sinusoidal network for blood transport. Our study characterizes liver tissue as a biological example of a biaxial liquid crystal. The general methodology developed here could be applied to other tissues and in-vitro organoids., Author summary Cell polarity enables cells to carry out specific functions. Cell polarity is characterized by the formation of different plasma membrane domains, each with specific composition of proteins, phospholipids and cytoskeletal components. In simple epithelial sheets, or tube-like tissues such as kidney, epithelial cells are known to display a single apical domain, facing a lumenal cavity, and a single basal domain on the opposite side of the cell, facing a basal layer of extracellular matrix. This apico-basal polarity defines a vector of cell polarity, which provides a direction of fluid transport, e.g., from the basal side of the sheet to the lumen-facing side. In more complex, three-dimensional epithelial tissues, such as liver tissue with its complex network of blood-transporting sinusoids, the membrane domains of hepatocyte cells display more intricate patterns, including rings and antipodal pairs of apical membrane. Here, we develop a mathematical framework to precisely characterize and quantify complex polarity patterns. Thereby, we reveal ordered patterns of cell polarity that span across a liver lobule. Our new method builds on physical concepts originally developed for ordered phases of liquid crystals. It provides a versatile tool to characterize the spatial organization of a complex three-dimensional tissue.

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.Author summaryThis study presents an open-source, extensible, multiscale platform for simulating viral immune interactions in epithelial tissues, which enables the rapid development and deployment of sophisticated models of viruses, infection mechanisms and tissue types. The model is used to investigate how potential treatments influence disease progression. Simulation results suggest that drugs which interfere with virus replication (e.g., remdesivir) yield substantially better infection outcomes when administered prophylactically even at very low doses than when used at high doses as treatment for an infection that has already begun.

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.

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, which is controlled by Hippo signalling. 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 of F-actin and phospho-Myosin, to compensate an overload of bile acids. Interestingly, these changes are sensed by the Hippo transcriptional co-activator YAP, which localizes to the apical F-actin-rich region and translocates to the nucleus in dependence of the acto-myosin system. 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 signalling, we explained this behaviour 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.

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.

Functional tissue architecture originates by self-assembly of distinct cell types, following tissue-specific rules of cell-cell interactions. In the liver, a structural model of the lobule was pioneered by Elias in 1949. This model, however, is in contrast with the apparent random 3D arrangement of hepatocytes. Since then, no significant progress has been made to derive the organizing principles of liver tissue. To solve this outstanding problem, we computationally reconstructed 3D tissue geometry from microscopy images and analyzed it applying soft-condensed-matter-physics concepts. Surprisingly, analysis of the spatial organization of cell polarity revealed that hepatocytes are not randomly oriented but follow a long-range liquid-crystal order. This does not depend exclusively on hepatocytes receiving instructive signals by endothelial cells as generally assumed, since silencing Integrin-ß1 disrupted both liquid-crystal order and organization of the sinusoidal network. Our results suggest that bi-directional communication between hepatocytes and sinusoids underlies the self-organization of liver tissue.

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.

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

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

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.

Early disease diagnosis is key for the effective treatment of diseases. It relies on the identification of biomarkers and morphological inspection of organs and tissues. Histopathological analysis of human biopsies is the gold standard to diagnose tissue alterations. However, this approach has low resolution and overlooks 3D structural changes that are consequence of functional alterations. Here, we applied multiphoton imaging, 3D digital reconstructions and computational simulations to generate spatially-resolved geometrical and functional models of human liver tissue at different stages of non-alcoholic fatty liver disease (NAFLD). We identified a set of new morphometric cellular parameters correlated with disease progression. Moreover, we found profound topological defects in the 3D bile canaliculi (BC) network. Personalized biliary fluid dynamic simulations predicted an increased pericentral biliary pressure and zonated cholestasis, consistent with elevated cholestatic biomarkers in patients’ sera. Our spatially-resolved models of human liver tissue can contribute to high-definition medicine by identifying quantitative multi-parametric cellular and tissue signatures to define disease progression and provide new insights into NAFLD pathophysiology.