Phenotypic characterization of liver tissue heterogeneity through a next-generation 3D single-cell atlas.

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.

Apical bulkheads accumulate as adaptive response to impaired bile flow in liver disease.

Hepatocytes form bile canaliculi that dynamically respond to the signalling activity of bile acids and bile flow. Little is known about their responses to intraluminal pressure. During embryonic development, hepatocytes assemble apical bulkheads that increase the canalicular resistance to intraluminal pressure. Here, we investigate whether they also protect bile canaliculi against elevated pressure upon impaired bile flow in adult liver. Apical bulkheads accumulate upon bile flow obstruction in mouse models and patients with primary sclerosing cholangitis (PSC). Their loss under these conditions leads to abnormally dilated canaliculi, resembling liver cell rosettes described in other hepatic diseases. 3D reconstruction reveals that these structures are sections of cysts and tubes formed by hepatocytes. Mathematical modelling establishes that they positively correlate with canalicular pressure and occur in early PSC stages. Using primary hepatocytes and 3D organoids, we demonstrate that excessive canalicular pressure causes the loss of apical bulkheads and formation of rosettes. Our results suggest that apical bulkheads are a protective mechanism of hepatocytes against impaired bile flow, highlighting the role of canalicular pressure in liver diseases.

The Rab5 effector FERRY links early endosomes with mRNA localization.

Localized translation is vital to polarized cells and requires precise and robust distribution of different mRNAs and ribosomes across the cell. However, the underlying molecular mechanisms are poorly understood and important players are lacking. Here, we discovered a Rab5 effector, the five-subunit endosomal Rab5 and RNA/ribosome intermediary (FERRY) complex, that recruits mRNAs and ribosomes to early endosomes through direct mRNA-interaction. FERRY displays preferential binding to certain groups of transcripts, including mRNAs encoding mitochondrial proteins. Deletion of FERRY subunits reduces the endosomal localization of transcripts in cells and has a significant impact on mRNA levels. Clinical studies show that genetic disruption of FERRY causes severe brain damage. We found that, in neurons, FERRY co-localizes with mRNA on early endosomes, and mRNA loaded FERRY-positive endosomes are in close proximity of mitochondria. FERRY thus transforms endosomes into mRNA carriers and plays a key role in regulating mRNA distribution and transport.

Structural basis of mRNA binding by the human FERRY Rab5 effector complex.

The pentameric FERRY Rab5 effector complex is a molecular link between mRNA and early endosomes in mRNA intracellular distribution. Here, we determine the cryo-EM structure of human FERRY. It reveals a unique clamp-like architecture that bears no resemblance to any known structure of Rab effectors. A combination of functional and mutational studies reveals that while the Fy-2 C-terminal coiled-coil acts as binding region for Fy-1/3 and Rab5, both coiled-coils and Fy-5 concur to bind mRNA. Mutations causing truncations of Fy-2 in patients with neurological disorders impair Rab5 binding or FERRY complex assembly. Thus, Fy-2 serves as a binding hub connecting all five complex subunits and mediating the binding to mRNA and early endosomes via Rab5. Our study provides mechanistic insights into long-distance mRNA transport and demonstrates that the particular architecture of FERRY is closely linked to a previously undescribed mode of RNA binding, involving coiled-coil domains.

Hepatocyte apical bulkheads provide a mechanical means to oppose bile pressure.

Hepatocytes grow their apical surfaces anisotropically to generate a 3D network of bile canaliculi (BC). BC elongation is ensured by apical bulkheads, membrane extensions that traverse the lumen and connect juxtaposed hepatocytes. We hypothesize that apical bulkheads are mechanical elements that shape the BC lumen in liver development but also counteract elevated biliary pressure. Here, by resolving their structure using STED microscopy, we found that they are sealed by tight junction loops, connected by adherens junctions, and contain contractile actomyosin, characteristics of mechanical function. Apical bulkheads persist at high pressure upon microinjection of fluid into the BC lumen, and laser ablation demonstrated that they are under tension. A mechanical model based on ablation results revealed that apical bulkheads double the pressure BC can hold. Apical bulkhead frequency anticorrelates with BC connectivity during mouse liver development, consistent with predicted changes in biliary pressure. Our findings demonstrate that apical bulkheads are load-bearing mechanical elements that could protect the BC network against elevated pressure.

Design centering enables robustness screening of pattern formation models.

MOTIVATION: Access to unprecedented amounts of quantitative biological data allows us to build and test biochemically accurate reaction-diffusion models of intracellular processes. However, any increase in model complexity increases the number of unknown parameters and, thus, the computational cost of model analysis. To efficiently characterize the behavior and robustness of models with many unknown parameters remains, therefore, a key challenge in systems biology. RESULTS: We propose a novel computational framework for efficient high-dimensional parameter space characterization of reaction-diffusion models in systems biology. The method leverages the Lp-Adaptation algorithm, an adaptive-proposal statistical method for approximate design centering and robustness estimation. Our approach is based on an oracle function, which predicts for any given point in parameter space whether the model fulfills given specifications. We propose specific oracles to efficiently predict four characteristics of Turing-type reaction-diffusion models: bistability, instability, capability of spontaneous pattern formation and capability of pattern maintenance. We benchmark the method and demonstrate that it enables global exploration of a model's ability to undergo pattern-forming instabilities and to quantify robustness for model selection in polynomial time with dimensionality. We present an application of the framework to pattern formation on the endosomal membrane by the small GTPase Rab5 and its effectors, and we propose molecular mechanisms underlying this system. AVAILABILITY AND IMPLEMENTATION: Our code is implemented in MATLAB and is available as open source under https://git.mpi-cbg.de/mosaic/software/black-box-optimization/rd-parameter-space-screening. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Active APPL1 sequestration by Plasmodium favors liver-stage development.

Intracellular pathogens manipulate host cells to survive and thrive. Cellular sensing and signaling pathways are among the key host machineries deregulated to favor infection. In this study, we show that liver-stage Plasmodium parasites compete with the host to sequester a host endosomal-adaptor protein (APPL1) known to regulate signaling in response to endocytosis. The enrichment of APPL1 at the parasitophorous vacuole membrane (PVM) involves an atypical Plasmodium Rab5 isoform (Rab5b). Depletion of host APPL1 alters neither the infection nor parasite development; however, upon overexpression of a GTPase-deficient host Rab5 mutant (hRab5_Q79L), the parasites are smaller and their PVM is stripped of APPL1. Infection with the GTPase-deficient Plasmodium berghei Rab5b mutant (PbRab5b_Q91L) in this case rescues the PVM APPL1 signal and parasite size. In summary, we observe a robust correlation between the level of APPL1 retention at the PVM and parasite size during exoerythrocytic development.

Endosomal escape of delivered mRNA from endosomal recycling tubules visualized at the nanoscale.

Delivery of exogenous mRNA using lipid nanoparticles (LNPs) is a promising strategy for therapeutics. However, a bottleneck remains in the poor understanding of the parameters that correlate with endosomal escape versus cytotoxicity. To address this problem, we compared the endosomal distribution of six LNP-mRNA formulations of diverse chemical composition and efficacy, similar to those used in mRNA-based vaccines, in primary human adipocytes, fibroblasts, and HeLa cells. Surprisingly, we found that total uptake is not a sufficient predictor of delivery, and different LNPs vary considerably in endosomal distributions. Prolonged uptake impaired endosomal acidification, a sign of cytotoxicity, and caused mRNA to accumulate in compartments defective in cargo transport and unproductive for delivery. In contrast, early endocytic/recycling compartments have the highest probability for mRNA escape. By using super-resolution microscopy, we could resolve a single LNP-mRNA within subendosomal compartments and capture events of mRNA escape from endosomal recycling tubules. Our results change the view of the mechanisms of endosomal escape and define quantitative parameters to guide the development of mRNA formulations toward higher efficacy and lower cytotoxicity.

Quantitative intracellular retention of delivered RNAs through optimized cell fixation and immunostaining.

Detection of nucleic acids within subcellular compartments is key to understanding their function. Determining the intracellular distribution of nucleic acids requires quantitative retention and estimation of their association with different organelles by immunofluorescence microscopy. This is particularly important for the delivery of nucleic acid therapeutics, which depends on endocytic uptake and endosomal escape. However, the current protocols fail to preserve the majority of exogenously delivered nucleic acids in the cytoplasm. To solve this problem, by monitoring Cy5-labeled mRNA delivered to primary human adipocytes via lipid nanoparticles (LNP), we optimized cell fixation, permeabilization, and immunostaining of a number of organelle markers, achieving quantitative retention of mRNA and allowing visualization of levels that escape detection using conventional procedures. The optimized protocol proved effective on exogenously delivered siRNA, miRNA, as well as endogenous miRNA. Our protocol is compatible with RNA probes of single molecule fluorescence in situ hybridization (smFISH) and molecular beacon, thus demonstrating that it is broadly applicable to study a variety of nucleic acids in cultured cells.

Correction: Quantification of nematic cell polarity in three-dimensional tissues.

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

Profiling Structural Alterations During Rab5 Nucleotide Exchange by HDX-MS.

Hydrogen deuterium exchange mass spectrometry (HDX-MS) gives insight into the structure of proteins. By monitoring the rate of exchange of the amide hydrogens on the protein backbone with deuterium atoms in the solvent, one can determine if a given region is highly structured or dynamic, map binding sites of interacting molecules or determine if a binding event is associated with allosteric structural alterations in a protein. Herein, we illustrate the use of this technique to monitor the nucleotide exchange process in Rab5, using the guanine nucleotide exchange factor (GEF)-effector complex, Rabex5:Rabaptin5. By simultaneously monitoring the HDX in Rab5, Rabex5 and Rabaptin5, we can directly visualize nucleotide exchange in Rab5, gain mechanistic insights into the exchange reaction and, by witnessing the transfer of Rab5 from Rabex5 to Rabaptin5, provide direct evidence for the positive feedback loop generated by a GEF-effector complex. HDX-MS can be used to monitor a variety of Rab protein-effector and -regulator interactions and be widely applied to other enzymatic processes as well.

Anisotropic expansion of hepatocyte lumina enforced by apical bulkheads.

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.

Quantification of nematic cell polarity in three-dimensional tissues.

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.

Resilience of three-dimensional sinusoidal networks in liver tissue.

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.

A non-linear system patterns Rab5 GTPase on the membrane.

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.

Bile canaliculi remodeling activates YAP via the actin cytoskeleton during liver regeneration.

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.

A drug discovery platform to identify compounds that inhibit EGFR triple mutants.

Receptor tyrosine kinases (RTKs) are transmembrane receptors of great clinical interest due to their role in disease. Historically, therapeutics targeting RTKs have been identified using in vitro kinase assays. Due to frequent development of drug resistance, however, there is a need to identify more diverse compounds that inhibit mutated but not wild-type RTKs. Here, we describe MaMTH-DS (mammalian membrane two-hybrid drug screening), a live-cell platform for high-throughput identification of small molecules targeting functional protein-protein interactions of RTKs. We applied MaMTH-DS to an oncogenic epidermal growth factor receptor (EGFR) mutant resistant to the latest generation of clinically approved tyrosine kinase inhibitors (TKIs). We identified four mutant-specific compounds, including two that would not have been detected by conventional in vitro kinase assays. One of these targets mutant EGFR via a new mechanism of action, distinct from classical TKI inhibition. Our results demonstrate how MaMTH-DS is a powerful complement to traditional drug screening approaches.

Three-dimensional spatially resolved geometrical and functional models of human liver tissue reveal new aspects of NAFLD progression.

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.