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.

Validation of AI-based software for objectification of conjunctival provocation test.

Background: Provocation tests are widely used in allergology to objectively reveal patients' sensitivity to specific allergens. The objective quantification of an allergic reaction is a crucial characteristic of these tests. Because of the absence of objective quantitative measurements, the conjunctival provocation test (CPT) is a less frequently used method despite its sensitivity and simplicity. Objective: We developed a new artificial intelligence (AI)-based method, called AllergoEye, for quantitative evaluation of conjunctival allergic reactions and validated it in a clinical study. Methods: AllergoEye was implemented as a 2-component system. The first component is based on an Android smartphone camera for screening and imaging the patient's eye, and the second is personal computer-based for image analysis and quantification. For the validation of AllergoEye, an open-label, prospective, monocentric study was carried out on 41 patients. Standardized CPT was performed with sequential titration of grass allergens in 4 dilutions, with the reaction evaluated by subjective/qualitative symptom scores and by quantitative AllergoEye scores. Results: AllergoEye demonstrated high sensitivity (98%) and specificity (90%) as compared with human estimation of allergic reaction. Tuning cutoff thresholds allowed us to increase the specificity of AllergoEye to 97%, at which point the correlation between detected sensitivity to allergen and specific IgE carrier-polymer system class becomes obvious. Strikingly, such correlation was not found with sensitivity to allergen detected on the basis of subjective and qualitative symptom scores. Conclusion: The clinical validation demonstrated that AllergoEye is a sensitive and efficient instrument for objective measurement of allergic reactions in CPT for clinical studies as well as for routine therapy control.

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.

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.

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.

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

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

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.

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.

SNX27-retromer assembly recycles MT1-MMP to invadopodia and promotes breast cancer metastasis.

A variety of metastatic cancer cells use actin-rich membrane protrusions, known as invadopodia, for efficient ECM degradation, which involves trafficking of proteases from intracellular compartments to these structures. Here, we demonstrate that in the metastatic breast cancer cell line MDA-MB-231, retromer regulates the matrix invasion activity by recycling matrix metalloprotease, MT1-MMP. We further found that MT2-MMP, another abundantly expressed metalloprotease, is also invadopodia associated. MT1- and MT2-MMP showed a high degree of colocalization but were located on the distinct endosomal domains. Retromer and its associated sorting nexin, SNX27, phenocopied each other in matrix degradation via selectively recycling MT1-MMP but not MT2-MMP. ITC-based studies revealed that both SNX27 and retromer could directly interact with MT1-MMP. Analysis from a publicly available database showed SNX27 to be overexpressed or frequently altered in the patients having invasive breast cancer. In xenograft-based studies, SNX27-depleted cell lines showed prolonged survival of SCID mice, suggesting a possible implication for overexpression of the sorting nexin in tumor samples.

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.

Auto-regulation of Rab5 GEF activity in Rabex5 by allosteric structural changes, catalytic core dynamics and ubiquitin binding.

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.

Liquid-crystal organization of liver tissue.

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.

Correlative single-molecule localization microscopy and electron tomography reveals endosome nanoscale domains.

Many cellular organelles, including endosomes, show compartmentalization into distinct functional domains, which, however, cannot be resolved by diffraction-limited light microscopy. Single molecule localization microscopy (SMLM) offers nanoscale resolution but data interpretation is often inconclusive when the ultrastructural context is missing. Correlative light electron microscopy (CLEM) combining SMLM with electron microscopy (EM) enables correlation of functional subdomains of organelles in relation to their underlying ultrastructure at nanometer resolution. However, the specific demands for EM sample preparation and the requirements for fluorescent single-molecule photo-switching are opposed. Here, we developed a novel superCLEM workflow that combines triple-color SMLM (dSTORM & PALM) and electron tomography using semi-thin Tokuyasu thawed cryosections. We applied the superCLEM approach to directly visualize nanoscale compartmentalization of endosomes in HeLa cells. Internalized, fluorescently labeled Transferrin and EGF were resolved into morphologically distinct domains within the same endosome. We found that the small GTPase Rab5 is organized in nanodomains on the globular part of early endosomes. The simultaneous visualization of several proteins in functionally distinct endosomal sub-compartments demonstrates the potential of superCLEM to link the ultrastructure of organelles with their molecular organization at nanoscale resolution.

Retrograde transport of Akt by a neuronal Rab5-APPL1 endosome.

Long-distance axonal trafficking plays a critical role in neuronal function and transport defects have been linked to neurodegenerative disorders. Various lines of evidence suggest that the small GTPase Rab5 plays a role in neuronal signaling via early endosomal transport. Here, we characterized the motility of Rab5 endosomes in primary cultures of mouse hippocampal pyramidal cells by live-cell imaging and showed that they exhibit bi-directional long-range motility in axons, with a strong bias toward retrograde transport. Characterization of key Rab5 effectors revealed that endogenous Rabankyrin-5, Rabenosyn-5 and APPL1 are all present in axons. Further analysis of APPL1-positive endosomes showed that, similar to Rab5-endosomes, they display more frequent long-range retrograde than anterograde movement, with the endosomal levels of APPL1 correlated with faster retrograde movement. Interestingly, APPL1-endosomes transport the neurotrophin receptor TrkB and mediate retrograde axonal transport of the kinase Akt1. FRET analysis revealed that APPL1 and Akt1 interact in an endocytosis-dependent manner. We conclude that Rab5-APPL1 endosomes exhibit the hallmarks of axonal signaling endosomes to transport Akt1 in hippocampal pyramidal cells.