Search

Evolutionarily conserved signaling pathways are essential for early embryogenesis, and reducing or abolishing their activity leads to characteristic developmental defects. Classification of phenotypic defects can identify the underlying signaling mechanisms, but this requires expert knowledge and the classification schemes have not been standardized. Here we use a machine learning approach for automated phenotyping to train a deep convolutional neural network, EmbryoNet, to accurately identify zebrafish signaling mutants in an unbiased manner. Combined with a model of time-dependent developmental trajectories, this approach identifies and classifies with high precision phenotypic defects caused by loss of function of the seven major signaling pathways relevant for vertebrate development. Our classification algorithms have wide applications in developmental biology and robustly identify signaling defects in evolutionarily distant species. Furthermore, using automated phenotyping in high-throughput drug screens, we show that EmbryoNet can resolve the mechanism of action of pharmaceutical substances. As part of this work, we freely provide more than 2 million images that were used to train and test EmbryoNet. more...

Three dimensional (3D) geometrical models are not only a powerful tool for quantitatively characterizing complex tissues but also useful for probing structure-function relationships in a tissue. However, these models are generally incomplete due to experimental limitations in acquiring multiple (>4) fluorescent channels simultaneously. Indeed, predictive geometrical and functional models of the liver have been restricted to few tissue and cellular components, excluding important cellular populations such as hepatic stellate cells (HSCs) and Kupffer cells (KCs). Here, we performed deep-tissue immunostaining, multiphoton microscopy, deeplearning techniques, and 3D image processing to computationally expand the number of simultaneously reconstructed tissue structures. We then generated a spatio-temporal singlecell atlas of hepatic architecture (Hep3D), including all main tissue and cellular components at different stages of post-natal development in mice. We used Hep3D to quantitatively study 1) hepatic morphodynamics from early post-natal development to adulthood, and 2) the structural role of KCs in the murine liver homeostasis. In addition to a complete description of bile canaliculi and sinusoidal network remodeling, our analysis uncovered unexpected spatiotemporal patterns of non-parenchymal cells and hepatocytes differing in size, number of nuclei, and DNA content. Surprisingly, we found that the specific depletion of KCs alters the number and morphology of the HSCs. These findings reveal novel characteristics of liver heterogeneity and have important implications for both the structural organization of liver tissue and its function. Our next-gen 3D single-cell atlas is a powerful tool to understand liver tissue architecture, under both physiological and pathological conditions. more...

During vertebrate embryogenesis, the germ layers are patterned by secreted Nodal signals. In the classical model, Nodals elicit signaling by binding to a complex comprising Type I/II Activin receptors (Acvr) and the co-receptor Tdgf1. However, it is currently unclear whether receptor binding can also affect the distribution of Nodals themselves through the embryo, and it is unknown which of the putative Acvr paralogs mediate Nodal signaling in zebrafish. Here, we characterize three Type I (Acvr1) and four Type II (Acvr2) homologs and show that – except for Acvr1c – all receptor-encoding transcripts are maternally deposited and present during zebrafish embryogenesis. We generated mutants and used them together with combinatorial morpholino knockdown and CRISPR F0 knockout (KO) approaches to assess compound loss-of-function phenotypes. We discovered that the Acvr2 homologs function partly redundantly and partially independently of Nodal to pattern the early zebrafish embryo, whereas the Type I receptors Acvr1b-a and Acvr1b-b redundantly act as major mediators of Nodal signaling. By combining quantitative analyses with expression manipulations, we found that feedback-regulated Type I receptors and co-receptors can directly influence the diffusion and distribution of Nodals, providing a mechanism for the spatial restriction of Nodal signaling during germ layer patterning., This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement No 637840 (QUANTPATTERN) and grant agreement No 863952 (ACE-OF-SPACE)). This work was also funded by the Max Planck Society and the International Max Planck Research School “From Molecules to Organisms”. more...

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. more...

[This corrects the article DOI: 10.1371/journal.pcbi.1008412.].

Belicova et al. report previously unrecognized apical membrane extensions, forming a pattern reminiscent of the bulkheads of boats, that determine the anisotropic expansion of hepatocyte lumina. Loss of the bulkheads caused by Rab35 silencing leads to re-engineering of epithelial polarity and liver tissue architecture., Lumen morphogenesis results from the interplay between molecular pathways and mechanical forces. In several organs, epithelial cells share their apical surfaces to form a tubular lumen. In the liver, however, hepatocytes share the apical surface only between adjacent cells and form narrow lumina that grow anisotropically, generating a 3D network of bile canaliculi (BC). Here, by studying lumenogenesis in differentiating mouse hepatoblasts in vitro, we discovered that adjacent hepatocytes assemble a pattern of specific extensions of the apical membrane traversing the lumen and ensuring its anisotropic expansion. These previously unrecognized structures form a pattern, reminiscent of the bulkheads of boats, also present in the developing and adult liver. Silencing of Rab35 resulted in loss of apical bulkheads and lumen anisotropy, leading to cyst formation. Strikingly, we could reengineer hepatocyte polarity in embryonic liver tissue, converting BC into epithelial tubes. Our results suggest that apical bulkheads are cell-intrinsic anisotropic mechanical elements that determine the elongation of BC during liver tissue morphogenesis., Graphical Abstract more...

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. more...

Can three-dimensional, microvasculature networks still ensure blood supply if individual links fail? We address this question in the sinusoidal network, a plexus-like microvasculature network, which transports nutrient-rich blood to every hepatocyte in liver tissue, by building on recent advances in high-resolution imaging and digital reconstruction of adult mice liver tissue. We find that the topology of the three-dimensional sinusoidal network reflects its two design requirements of a space-filling network that connects all hepatocytes, while using shortest transport routes: sinusoidal networks are sub-graphs of the Delaunay graph of their set of branching points, and also contain the corresponding minimum spanning tree, both to good approximation. To overcome the spatial limitations of experimental samples and generate arbitrarily-sized networks, we developed a network generation algorithm that reproduces the statistical features of 0.3-mm-sized samples of sinusoidal networks, using multi-objective optimization for node degree and edge length distribution. Nematic order in these simulated networks implies anisotropic transport properties, characterized by an empirical linear relation between a nematic order parameter and the anisotropy of the permeability tensor. Under the assumption that all sinusoid tubes have a constant and equal flow resistance, we predict that the distribution of currents in the network is very inhomogeneous, with a small number of edges carrying a substantial part of the flow—a feature known for hierarchical networks, but unexpected for plexus-like networks. We quantify network resilience in terms of a permeability-at-risk, i.e., permeability as function of the fraction of removed edges. We find that sinusoidal networks are resilient to random removal of edges, but vulnerable to the removal of high-current edges. Our findings suggest the existence of a mechanism counteracting flow inhomogeneity to balance metabolic load on the liver., Author summary The liver is the largest metabolic organ of the human body and pivotal for blood detoxification and the uptake of many medically relevant drugs. Inside liver tissue, a dense meshwork of blood micro-vessels, the sinusoidal network, contacts every hepatocyte cell. The architecture of these networks, and especially how their function responds to local damage is not well understood. Previous theoretical work addressed network resilience in simple two-dimensional networks such as vein networks of leaves but not in three-dimensional networks, partly due to the considerable difficulties of imaging three-dimensional tissues. Here, we build on unprecedented advances in imaging mouse liver tissue. By simulating arbitrarily-sized networks that faithfully reproduce the statistical features of spatially restricted experimental samples, we computationally characterize the relation between a weak alignment of the network along the flow direction and enhanced transport along this direction. A simple transport model predicts an inhomogeneous distribution of flow in the sinusoidal network. Concomitantly, these networks are resilient against the random removal of edges, but vulnerable against removal of ‘highway’ edges that carried a high flow in the unperturbed state. We speculate that yet unknown adaptive mechanisms balance the distribution of flows and thereby increase the resilience of this physiologically important network. more...

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. more...

Fluorescently targeted proteins are widely used for studies of intracellular organelles dynamic. Peripheral proteins are transiently associated with organelles and a significant fraction of them are located at the cytosol. Image analysis of peripheral proteins poses a problem on properly discriminating membrane-associated signal from the cytosolic one. In most cases, signals from organelles are compact in comparison with diffuse signal from cytosol. Commonly used methods for background estimation depend on the assumption that background and foreground signals are separable by spatial frequency filters. However, large non-stained organelles (e.g., nuclei) result in abrupt changes in the cytosol intensity and lead to errors in the background estimation. Such mistakes result in artifacts in the reconstructed foreground signal. We developed a new algorithm that estimates background intensity in fluorescence microscopy images and does not produce artifacts on the borders of nuclei. more...

Plasticity is an essential condition for cancer cells to invade surrounding tissues. The nucleus is the most rigid cellular organelle and it undergoes substantial deformations to get through environmental constrictions. Nuclear stiffness mostly depends on the nuclear lamina and chromatin, which in turn might be affected by nuclear architectural proteins. Among these is the HMGA1 (High Mobility Group A1) protein, a factor that plays a causal role in neoplastic transformation and that is able to disentangle heterochromatic domains by H1 displacement. Here we made use of atomic force microscopy to analyze the stiffness of breast cancer cellular models in which we modulated HMGA1 expression to investigate its role in regulating nuclear plasticity. Since histone H1 is the main modulator of chromatin structure and HMGA1 is a well-established histone H1 competitor, we correlated HMGA1 expression and cellular stiffness with histone H1 expression level, post-translational modifications, and nuclear distribution. Our results showed that HMGA1 expression level correlates with nuclear stiffness, is associated to histone H1 phosphorylation status, and alters both histone H1 chromatin distribution and expression. These data suggest that HMGA1 might promote chromatin relaxation through a histone H1-mediated mechanism strongly impacting on the invasiveness of cancer cells. more...

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. more...

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. more...

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. more...

A prerequisite for the systems biology analysis of tissues is an accurate digital three-dimensional reconstruction of tissue structure based on images of markers covering multiple scales. Here, we designed a flexible pipeline for the multi-scale reconstruction and quantitative morphological analysis of tissue architecture from microscopy images. Our pipeline includes newly developed algorithms that address specific challenges of thick dense tissue reconstruction. Our implementation allows for a flexible workflow, scalable to high-throughput analysis and applicable to various mammalian tissues. We applied it to the analysis of liver tissue and extracted quantitative parameters of sinusoids, bile canaliculi and cell shapes, recognizing different liver cell types with high accuracy. Using our platform, we uncovered an unexpected zonation pattern of hepatocytes with different size, nuclei and DNA content, thus revealing new features of liver tissue organization. The pipeline also proved effective to analyse lung and kidney tissue, demonstrating its generality and robustness. DOI: http://dx.doi.org/10.7554/eLife.11214.001, eLife digest Understanding how individual cells interact to form tissues in animals and plants is a key problem in cell and developmental biology. To be able to answer this question researchers need to use microscopy to observe the cells in a tissue, extract structural information from the images, and then generate three-dimensional digital models of the tissue. However, the software solutions that are currently available are limited, and reconstructing three-dimensional tissue from microscopy images remains problematic. To meet this challenge, Morales-Navarrete et al. extended the free software platform called MotionTracking, which had been used previously for two-dimensional work. The software now combines a series of new and established algorithms for analysing fluorescence microscopy images that make it possible to identify the different structures that make up a tissue and then create and analyse a three-dimensional model. Morales-Navarrete et al. used the software to analyse liver tissue from mice. The resulting model revealed that liver cells called hepatocytes are arranged in particular zones within the tissue according to their size and DNA content. The software was also applied successfully to analyse lung and kidney tissue, which demonstrates that the approach can be used to create three-dimensional models of a variety of tissues. Morales-Navarrete et al.’s approach can rapidly generate accurate models of larger tissues than were previously possible. Therefore, it provides researchers with a powerful tool to analyse the different features of tissues. This tool will be useful for many areas of research: from understanding of how cells form tissues, to diagnosing diseases based on the changes to features in particular tissues. DOI: http://dx.doi.org/10.7554/eLife.11214.002 more...

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 more...

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. more...