Search

Abstract Three-dimensional (3D) geometrical models are potent tools for quantifying complex tissue features and exploring structure–function relationships. However, these models are generally incomplete due to experimental limitations in acquiring multiple (> 4) fluorescent channels in thick tissue sections 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 combined deep-tissue immunostaining, multiphoton microscopy, deep-learning techniques, and 3D image processing to computationally expand the number of simultaneously reconstructed tissue structures. We then generated a spatial single-cell 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 effect on the liver's overall structure when changing the hepatic environment after removing KCs. 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 results in morphological changes in hepatocytes and 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, opening up avenues for in-depth investigations into tissue structure across both normal and pathological conditions.

Abstract 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 corrects the article DOI: 10.1371/journal.pcbi.1008412.].

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.

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.

Proteins can self-organize into spatial patterns via non-linear dynamic interactions on cellular membranes. Modelling and simulations have shown that small GTPases can generate patterns by coupling guanine nucleotide exchange factors (GEF) to effectors, generating a positive feedback of GTPase activation and membrane recruitment. Here, we reconstituted the patterning of the small GTPase Rab5 and its GEF/effector complex Rabex5/Rabaptin5 on supported lipid bilayers. We demonstrate a ‘handover’ of Rab5 from Rabex5 to Rabaptin5 upon nucleotide exchange. A minimal system consisting of Rab5, RabGDI and a complex of full length Rabex5/Rabaptin5 was necessary to pattern Rab5 into membrane domains. Rab5 patterning required a lipid membrane composition mimicking that of early endosomes, with PI(3)P enhancing membrane recruitment of Rab5 and acyl chain packing being critical for domain formation. The prevalence of GEF/effector coupling in nature suggests a possible universal system for small GTPase patterning involving both protein and lipid interactions.

Abstract Liver cells communicate with the extracellular environment to take up nutrients via endocytosis. Iron uptake is essential for metabolic activities and cell homeostasis. Here, we investigated the role of the endocytic system for maintaining iron homeostasis. We specifically depleted the small GTPase Rab5 in the mouse liver, causing a transient loss of the entire endo-lysosomal system. Strikingly, endosome depletion led to a fast reduction of hepatic iron levels, which was preceded by an increased abundance of the iron exporter ferroportin. Compensatory changes in livers of Rab5-depleted mice include increased expression of transferrin receptor 1 as well as reduced expression of the iron-regulatory hormone hepcidin. Serum iron indices (serum iron, free iron binding capacity and total iron binding capacity) in Rab5-KD mice were increased, consistent with an elevated splenic and hepatic iron export. Our data emphasize the critical importance of the endosomal compartments in hepatocytes to maintain hepatic and systemic iron homeostasis in vivo. The short time period (between day four and five) upon which these changes occur underscore the fast dynamics of the liver iron pool.

Intracellular trafficking depends on the function of Rab GTPases, whose activation is regulated by guanine exchange factors (GEFs). The Rab5 GEF, Rabex5, was previously proposed to be auto-inhibited by its C-terminus. Here, we studied full-length Rabex5 and Rabaptin5 proteins as well as domain deletion Rabex5 mutants using hydrogen deuterium exchange mass spectrometry. We generated a structural model of Rabex5, using chemical cross-linking mass spectrometry and integrative modeling techniques. By correlating structural changes with nucleotide exchange activity for each construct, we uncovered new auto-regulatory roles for the ubiquitin binding domains and the Linker connecting those domains to the catalytic core of Rabex5. We further provide evidence that enhanced dynamics in the catalytic core are linked to catalysis. Our results suggest a more complex auto-regulation mechanism than previously thought and imply that ubiquitin binding serves not only to position Rabex5 but to also control its Rab5 GEF activity through allosteric structural alterations.

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 of mouse liver tissue 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, 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.

Mitochondrial stress response is essential for cell survival, and damaged mitochondria are a hallmark of neurodegenerative diseases. Thus, it is fundamental to understand how mitochondria relay information within the cell. Here, by investigating mitochondrial-endosomal contact sites we made the surprising observation that the small GTPase Rab5 translocates from early endosomes to mitochondria upon oxidative stress. This process is reversible and accompanied by an increase in Rab5-positive endosomes in contact with mitochondria. Interestingly, activation of Rab5 on mitochondria depends on the Rab5-GEF ALS2/Alsin, encoded by a gene mutated in amyotrophic lateral sclerosis (ALS). Alsin-deficient human-induced pluripotent stem cell-derived spinal motor neurons are defective in relocating Rab5 to mitochondria and display increased susceptibility to oxidative stress. These findings define a novel pathway whereby Alsin catalyzes the assembly of the Rab5 endocytic machinery on mitochondria. Defects in stress-sensing by endosomes could be crucial for mitochondrial quality control during the onset of ALS.

The liver maintains glucose and lipid homeostasis by adapting its metabolic activity to the energy needs of the organism. Communication between hepatocytes and extracellular environment via endocytosis is key to such homeostasis. Here, we addressed the question of whether endosomes are required for gluconeogenic gene expression. We took advantage of the loss of endosomes in the mouse liver upon Rab5 silencing. Strikingly, we found hepatomegaly and severe metabolic defects such as hypoglycemia, hypercholesterolemia, hyperlipidemia, and glycogen accumulation that phenocopied those found in von Gierke’s disease, a glucose-6-phosphatase (G6Pase) deficiency. G6Pase deficiency alone can account for the reduction in hepatic glucose output and glycogen accumulation as determined by mathematical modeling. Interestingly, we uncovered functional alterations in the transcription factors, which regulate G6Pase expression. Our data highlight a requirement of Rab5 and the endosomal system for the regulation of gluconeogenic gene expression that has important implications for metabolic diseases.

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

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.

An outstanding question is how receptor tyrosine kinases (RTKs) determine different cell-fate decisions despite sharing the same signalling cascades. Here, we uncovered an unexpected mechanism of RTK trafficking in this process. By quantitative high-resolution FRET microscopy, we found that phosphorylated epidermal growth factor receptor (p-EGFR) is not randomly distributed but packaged at constant mean amounts in endosomes. Cells respond to higher EGF concentrations by increasing the number of endosomes but keeping the mean p-EGFR content per endosome almost constant. By mathematical modelling, we found that this mechanism confers both robustness and regulation to signalling output. Different growth factors caused specific changes in endosome number and size in various cell systems and changing the distribution of p-EGFR between endosomes was sufficient to reprogram cell-fate decision upon EGF stimulation. We propose that the packaging of p-RTKs in endosomes is a general mechanism to ensure the fidelity and specificity of the signalling response.

Functional genomics screens using multi-parametric assays are powerful approaches for identifying genes involved in particular cellular processes. However, they suffer from problems like noise, and often provide little insight into molecular mechanisms. A bottleneck for addressing these issues is the lack of computational methods for the systematic integration of multi-parametric phenotypic datasets with molecular interactions. Here, we present Integrative Multi Profile Analysis of Cellular Traits (IMPACT). The main goal of IMPACT is to identify the most consistent phenotypic profile among interacting genes. This approach utilizes two types of external information: sets of related genes (IMPACT-sets) and network information (IMPACT-modules). Based on the notion that interacting genes are more likely to be involved in similar functions than non-interacting genes, this data is used as a prior to inform the filtering of phenotypic profiles that are similar among interacting genes. IMPACT-sets selects the most frequent profile among a set of related genes. IMPACT-modules identifies sub-networks containing genes with similar phenotype profiles. The statistical significance of these selections is subsequently quantified via permutations of the data. IMPACT (1) handles multiple profiles per gene, (2) rescues genes with weak phenotypes and (3) accounts for multiple biases e.g. caused by the network topology. Application to a genome-wide RNAi screen on endocytosis showed that IMPACT improved the recovery of known endocytosis-related genes, decreased off-target effects, and detected consistent phenotypes. Those findings were confirmed by rescreening 468 genes. Additionally we validated an unexpected influence of the IGF-receptor on EGF-endocytosis. IMPACT facilitates the selection of high-quality phenotypic profiles using different types of independent information, thereby supporting the molecular interpretation of functional screens.

Rab GTPases constitute the largest subfamily of the Ras protein superfamily. Rab proteins regulate organelle biogenesis and transport, and display distinct binding preferences for effector and activator proteins, many of which have not been elucidated yet. The underlying molecular recognition motifs, binding partner preferences and selectivities are not well understood.Comparative analysis of the amino acid sequences and the three-dimensional electrostatic and hydrophobic molecular interaction fields of 62 human Rab proteins revealed a wide range of binding properties with large differences between some Rab proteins. This analysis assists the functional annotation of Rab proteins 12, 14, 26, 37 and 41 and provided an explanation for the shared function of Rab3 and 27. Rab7a and 7b have very different electrostatic potentials, indicating that they may bind to different effector proteins and thus, exert different functions. The subfamily V Rab GTPases which are associated with endosome differ subtly in the interaction properties of their switch regions, and this may explain exchange factor specificity and exchange kinetics.We have analysed conservation of sequence and of molecular interaction fields to cluster and annotate the human Rab proteins. The analysis of three dimensional molecular interaction fields provides detailed insight that is not available from a sequence-based approach alone. Based on our results, we predict novel functions for some Rab proteins and provide insights into their divergent functions and the determinants of their binding partner selectivity.

The small GTPase Rab5 is a key regulator of clathrin-mediated endocytosis. On early endosomes, within a spatially restricted domain enriched in phosphatydilinositol-3-phosphate [PI(3)P], Rab5 coordinates a complex network of effectors that functionally cooperate in membrane tethering, fusion, and organelle motility. Here we discovered a novel PI(3)P-binding Rab5 effector, Rabankyrin-5, which localises to early endosomes and stimulates their fusion activity. In addition to early endosomes, however, Rabankyrin-5 localises to large vacuolar structures that correspond to macropinosomes in epithelial cells and fibroblasts. Overexpression of Rabankyrin-5 increases the number of macropinosomes and stimulates fluid-phase uptake, whereas its downregulation inhibits these processes. In polarised epithelial cells, this function is primarily restricted to the apical membrane. Rabankyrin-5 localises to large pinocytic structures underneath the apical surface of kidney proximal tubule cells, and its overexpression in polarised Madin-Darby canine kidney cells stimulates apical but not basolateral, non-clathrin-mediated pinocytosis. In demonstrating a regulatory role in endosome fusion and (macro)pinocytosis, our studies suggest that Rab5 regulates and coordinates different endocytic mechanisms through its effector Rabankyrin-5. Furthermore, its active role in apical pinocytosis in epithelial cells suggests an important function of Rabankyrin-5 in the physiology of polarised cells.