Intravital Dynamic and Correlative Imaging of Mouse Livers Reveals Diffusion-Dominated Canalicular and Flow-Augmented Ductular Bile Flux

Background and Aims Small‐molecule flux in tissue microdomains is essential for organ function, but knowledge of this process is scant due to the lack of suitable methods. We developed two independent techniques that allow the quantification of advection (flow) and diffusion in individual bile canaliculi and in interlobular bile ducts of intact livers in living mice, namely fluorescence loss after photoactivation and intravital arbitrary region image correlation spectroscopy. Approach and Results The results challenge the prevailing “mechano‐osmotic” theory of canalicular bile flow. After active transport across hepatocyte membranes, bile acids are transported in the canaliculi primarily by diffusion. Only in the interlobular ducts is diffusion augmented by regulatable advection. Photoactivation of fluorescein bis‐(5‐carboxymethoxy‐2‐nitrobenzyl)‐ether in entire lobules demonstrated the establishment of diffusive gradients in the bile canalicular network and the sink function of interlobular ducts. In contrast to the bile canalicular network, vectorial transport was detected and quantified in the mesh of interlobular bile ducts. Conclusions The liver consists of a diffusion‐dominated canalicular domain, where hepatocytes secrete small molecules and generate a concentration gradient and a flow‐augmented ductular domain, where regulated water influx creates unidirectional advection that augments the diffusive flux

Interruption of bile acid uptake by hepatocytes after acetaminophen overdose ameliorates hepatotoxicity.

Background & aimsAcetaminophen (APAP) overdose remains a frequent cause of acute liver failure, which is generally accompanied by increased levels of serum bile acids (BAs). However, the pathophysiological role of BAs remains elusive. Herein, we investigated the role of BAs in APAP-induced hepatotoxicity.MethodsWe performed intravital imaging to investigate BA transport in mice, quantified endogenous BA concentrations in the serum of mice and patients with APAP overdose, analyzed liver tissue and bile by mass spectrometry and MALDI-mass spectrometry imaging, assessed the integrity of the blood-bile barrier and the role of oxidative stress by immunostaining of tight junction proteins and intravital imaging of fluorescent markers, identified the intracellular cytotoxic concentrations of BAs, and performed interventions to block BA uptake from blood into hepatocytes.ResultsPrior to the onset of cell death, APAP overdose causes massive oxidative stress in the pericentral lobular zone, which coincided with a breach of the blood-bile barrier. Consequently, BAs leak from the bile canaliculi into the sinusoidal blood, which is then followed by their uptake into hepatocytes via the basolateral membrane, their secretion into canaliculi and repeated cycling. This, what we termed 'futile cycling' of BAs, led to increased intracellular BA concentrations that were high enough to cause hepatocyte death. Importantly, however, the interruption of BA re-uptake by pharmacological NTCP blockage using Myrcludex B and Oatp knockout strongly reduced APAP-induced hepatotoxicity.ConclusionsAPAP overdose induces a breach of the blood-bile barrier which leads to futile BA cycling that causes hepatocyte death. Prevention of BA cycling may represent a therapeutic option after APAP intoxication.Lay summaryOnly one drug, N-acetylcysteine, is approved for the treatment of acetaminophen overdose and it is only effective when given within ∼8 hours after ingestion. We identified a mechanism by which acetaminophen overdose causes an increase in bile acid concentrations (to above toxic thresholds) in hepatocytes. Blocking this mechanism prevented acetaminophen-induced hepatotoxicity in mice and evidence from patients suggests that this therapy may be effective for longer periods after ingestion compared to N-acetylcysteine

A predictive computational model shows that biomechanical cell cycle progression control can explain liver regeneration after partial hepatectomy

Partial hepatectomy (PHx) is a surgical intervention where a part of liver is removed. Due to its extraordinary capacity to regenerate, the organ is able to regenerate about two-thirds of its mass within a few weeks. Nevertheless, in some patients, regeneration fails. Understanding the principles and limitations underlying regeneration may permit to control this process. We established a computational model to mimic and prospectively improve regeneration of the liver lobe of mice. This model represents each hepatocyte individually and builds upon a previous computational model of regeneration of drug (CCl4)-induced damage in a single liver lobule, a small subunit of a lobe. The present study simulates entire liver lobes that can consist of hundreds to thousands of lobules. Different from that previous model it accounts for biomechanical control of cell cycle progression (Biomechanical Growth Control). The model reproduced the available experimental observations only if BGC was taken into account. Interestingly, the model predicted that BGC minimizes the number of proliferating neighbor cells of a proliferating cell resulting in a checkerboard-like proliferation pattern. Moreover, the model predicted different cell proliferation patterns in pigs and mice that corresponded to data obtained from regenerating tissue of the two species. In conclusion, the here established model suggests that biomechanical control mechanisms may play a significant role in liver regeneration after PHx

Digital twin demonstrates significance of biomechanical growth control in liver regeneration after partial hepatectomy

Summary: Partial liver removal is an important therapy option for liver cancer. In most patients within a few weeks, the liver is able to fully regenerate. In some patients, however, regeneration fails with often severe consequences. To better understand the control mechanisms of liver regeneration, experiments in mice were performed, guiding the creation of a spatiotemporal 3D model of the regenerating liver. The model represents cells and blood vessels within an entire liver lobe, a macroscopic liver subunit. The model could reproduce the experimental data only if a biomechanical growth control (BGC)-mechanism, inhibiting cell cycle entrance at high compression, was taken into account and predicted that BGC may act as a short-range growth inhibitor minimizing the number of proliferating neighbor cells of a proliferating cell, generating a checkerboard-like proliferation pattern. Model-predicted cell proliferation patterns in pigs and mice were found experimentally. The results underpin the importance of biomechanical aspects in liver growth control

Automated Detection of Portal Fields and Central Veins in Whole-Slide Images of Liver Tissue

Many physiological processes and pathological phenomena in the liver tissue are spatially heterogeneous. At a local scale, biomarkers can be quantified along the axis of the blood flow, from portal fields (PFs) to central veins (CVs), i.e., in zonated form. This requires detecting PFs and CVs. However, manually annotating these structures in multiple whole-slide images is a tedious task. We describe and evaluate a fully automated method, based on a convolutional neural network, for simultaneously detecting PFs and CVs in a single stained section. Trained on scans of hematoxylin and eosin-stained liver tissue, the detector performed well with an F1 score of 0.81 compared to annotation by a human expert. It does, however, not generalize well to previously unseen scans of steatotic liver tissue with an F1 score of 0.59. Automated PF and CV detection eliminates the bottleneck of manual annotation for subsequent automated analyses, as illustrated by two proof-of-concept applications: We computed lobulus sizes based on the detected PF and CV positions, where results agreed with published lobulus sizes. Moreover, we demonstrate the feasibility of zonated quantification of biomarkers detected in different stainings based on lobuli and zones obtained from the detected PF and CV positions. A negative control (hematoxylin and eosin) showed the expected homogeneity, a positive control (glutamine synthetase) was quantified to be strictly pericentral, and a plausible zonation for a heterogeneous F4/80 staining was obtained. Automated detection of PFs and CVs is one building block for automatically quantifying physiologically relevant heterogeneity of liver tissue biomarkers. Perspectively, a more robust and automated assessment of zonation from whole-slide images will be valuable for parameterizing spatially resolved models of liver metabolism and to provide diagnostic information

The transcription factor CHOP, a central component of the transcriptional regulatory network induced upon CCl4 intoxication in mouse liver, is not a critical mediator of hepatotoxicity

Since xenobiotics enter the organism via the liver, hepatocytes must cope with numerous perturbations, including modifications of proteins leading to endoplasmic reticulum stress (ER-stress). This triggers a signaling pathway termed unfolded protein response (UPR) that aims to restore homeostasis or to eliminate disturbed hepatocytes by apoptosis. In the present study, we used the well-established CCl4 hepatotoxicity model in mice to address the questions whether CCl4 induces ER-stress and, if so, whether the well-known ER-stress effector CHOP is responsible for CCl4-induced apoptosis. For this purpose, we treated mice with a high dose of CCl4 injected i.p. and followed gene expression profile over time using Affymetrix gene array analysis. This time resolved gene expression analysis allowed the identification of gene clusters with overrepresented binding sites for the three most important ER-stress induced transcription factors, CHOP, XBP1 and ATF4. Such result was confirmed by the demonstration of CCl4-induced XBP1 splicing, upregulation of CHOP at mRNA and protein levels, and translocation of CHOP to the nucleus. Two observations indicated that CHOP may be responsible for CCl4-induced cell death: (1) Nuclear translocation of CHOP was exclusively observed in the pericentral fraction of hepatocytes that deteriorate in response to CCl4 and (2) CHOP-regulated genes with previously reported pro-apoptotic function such as GADD34, TRB3 and ERO1L were induced in the pericentral zone as well. Therefore, we compared CCl4 induced hepatotoxicity in CHOP knockout versus wild-type mice. Surprisingly, genetic depletion of CHOP did not afford protection against CCl4-induced damage as evidenced by serum GOT and GPT as well as quantification of dead tissue areas. The negative result was obtained at several time points (8, 24 and 72 h) and different CCl4 doses (1.6 and 0.132 g/kg). Overall, our results demonstrate that all branches of the UPR are activated in mouse liver upon CCl4 treatment. However, CHOP does not play a critical role in CCl4-induced cell death and cannot be considered as a biomarker strictly linked to hepatotoxicity. The role of alternative UPR effectors such as XBP1 remains to be investigated

Intravital Dynamic and Correlative Imaging of Mouse Livers Reveals Diffusion-Dominated Canalicular and Flow-Augmented Ductular Bile Flux

Background and Aims: Small-molecule flux in tissue microdomains is essential for organ function, but knowledge of this process is scant due to the lack of suitable methods. We developed two independent techniques that allow the quantification of advection (flow) and diffusion in individual bile canaliculi and in interlobular bile ducts of intact livers in living mice, namely fluorescence loss after photoactivation and intravital arbitrary region image correlation spectroscopy. Approach and Results: The results challenge the prevailing "mechano-osmotic" theory of canalicular bile flow. After active transport across hepatocyte membranes, bile acids are transported in the canaliculi primarily by diffusion. Only in the interlobular ducts is diffusion augmented by regulatable advection. Photoactivation of fluorescein bis-(5-carboxymethoxy-2-nitrobenzyl)-ether in entire lobules demonstrated the establishment of diffusive gradients in the bile canalicular network and the sink function of interlobular ducts. In contrast to the bile canalicular network, vectorial transport was detected and quantified in the mesh of interlobular bile ducts. Conclusions: The liver consists of a diffusion-dominated canalicular domain, where hepatocytes secrete small molecules and generate a concentration gradient and a flow-augmented ductular domain, where regulated water influx creates unidirectional advection that augments the diffusive flux

Transcriptomic Cross-Species Analysis of Chronic Liver Disease Reveals Consistent Regulation Between Humans and Mice

Mouse models are frequently used to study chronic liver diseases (CLDs). To assess their translational relevance, we quantified the similarity of commonly used mouse models to human CLDs based on transcriptome data. Gene-expression data from 372 patients were compared with data from acute and chronic mouse models consisting of 227 mice, and additionally to nine published gene sets of chronic mouse models. Genes consistently altered in humans and mice were mapped to liver cell types based on single-cell RNA-sequencing data and validated by immunostaining. Considering the top differentially expressed genes, the similarity between humans and mice varied among the mouse models and depended on the period of damage induction. The highest recall (0.4) and precision (0.33) were observed for the model with 12-months damage induction by CCl4 and by a Western diet, respectively. Genes consistently up-regulated between the chronic CCl4 model and human CLDs were enriched in inflammatory and developmental processes, and mostly mapped to cholangiocytes, macrophages, and endothelial and mesenchymal cells. Down-regulated genes were enriched in metabolic processes and mapped to hepatocytes. Immunostaining confirmed the regulation of selected genes and their cell type specificity. Genes that were up-regulated in both acute and chronic models showed higher recall and precision with respect to human CLDs than exclusively acute or chronic genes. CONCLUSION: Similarly regulated genes in human and mouse CLDs were identified. Despite major interspecies differences, mouse models detected 40% of the genes significantly altered in human CLD. The translational relevance of individual genes can be assessed at https://saezlab.shinyapps.io/liverdiseaseatlas/

Colitis ameliorates cholestatic liver disease via suppression of bile acid synthesis

Abstract Primary sclerosing cholangitis (PSC) is a chronic cholestatic liver disease characterized by chronic inflammation and progressive fibrosis of the biliary tree. The majority of PSC patients suffer from concomitant inflammatory bowel disease (IBD), which has been suggested to promote disease development and progression. However, the molecular mechanisms by which intestinal inflammation may aggravate cholestatic liver disease remain incompletely understood. Here, we employ an IBD-PSC mouse model to investigate the impact of colitis on bile acid metabolism and cholestatic liver injury. Unexpectedly, intestinal inflammation and barrier impairment improve acute cholestatic liver injury and result in reduced liver fibrosis in a chronic colitis model. This phenotype is independent of colitis-induced alterations of microbial bile acid metabolism but mediated via hepatocellular NF-κB activation by lipopolysaccharide (LPS), which suppresses bile acid metabolism in-vitro and in-vivo. This study identifies a colitis-triggered protective circuit suppressing cholestatic liver disease and encourages multi-organ treatment strategies for PSC

Protocols for staining of bile canalicular and sinusoidal networks of human, mouse and pig livers, three-dimensional reconstruction and quantification of tissue microarchitecture by image processing and analysis

Histological alterations often constitute a fingerprint of toxicity and diseases. The extent to which these alterations are cause or consequence of compromised organ function, and the underlying mechanisms involved is a matter of intensive research. In particular, liver disease is often associated with altered tissue microarchitecture, which in turn may compromise perfusion and functionality. Research in this field requires the development and orchestration of new techniques into standardized processing pipelines that can be used to reproducibly quantify tissue architecture. Major bottlenecks include the lack of robust staining, and adequate reconstruction and quantification techniques. To bridge this gap, we established protocols employing specific antibody combinations for immunostaining, confocal imaging, three-dimensional reconstruction of approximately 100-μm-thick tissue blocks and quantification of key architectural features. We describe a standard procedure termed ‘liver architectural staining’ for the simultaneous visualization of bile canaliculi, sinusoidal endothelial cells, glutamine synthetase (GS) for the identification of central veins, and DAPI as a nuclear marker. Additionally, we present a second standard procedure entitled ‘S-phase staining’, where S-phase-positive and S-phase-negative nuclei (stained with BrdU and DAPI, respectively), sinusoidal endothelial cells and GS are stained. The techniques include three-dimensional reconstruction of the sinusoidal and bile canalicular networks from the same tissue block, and robust capture of position, size and shape of individual hepatocytes, as well as entire lobules from the same tissue specimen. In addition to the protocols, we have also established image analysis software that allows relational and hierarchical quantifications of different liver substructures (e.g. cells and vascular branches) and events (e.g. cell proliferation and death). Typical results acquired for routinely quantified parameters in adult mice (C57Bl6/N) include the hepatocyte volume (5,128.3 ± 837.8 μm3) and the fraction of the hepatocyte surface in contact with the neighbouring hepatocytes (67.4 ± 6.7 %), sinusoids (22.1 ± 4.8 %) and bile canaliculi (9.9 ± 3.8 %). Parameters of the sinusoidal network that we also routinely quantify include the radius of the sinusoids (4.8 ± 2.25 μm), the branching angle (32.5 ± 11.2°), the length of intersection branches (23.93 ± 5.9 μm), the number of intersection nodes per mm3 (120.3 × 103 ± 42.1 × 103), the average length of sinusoidal vessel per mm3 (5.4 × 103 ± 1.4 × 103mm) and the percentage of vessel volume in relation to the whole liver volume (15.3 ± 3.9) (mean ± standard deviation). Moreover, the provided parameters of the bile canalicular network are: length of the first-order branches (7.5 ± 0.6 μm), length of the second-order branches (10.9 ± 1.8 μm), length of the dead-end branches (5.9 ± 0.7 μm), the number of intersection nodes per mm3 (819.1 × 103 ± 180.7 × 103), the number of dead-end branches per mm3 (409.9 × 103 ± 95.6 × 103), the length of the bile canalicular network per mm3 (9.4 × 103 ± 0.7 × 103 mm) and the percentage of the bile canalicular volume with respect to the total liver volume (3.4 ± 0.005). A particular strength of our technique is that quantitative parameters of hepatocytes and bile canalicular as well as sinusoidal networks can be extracted from the same tissue block. Reconstructions and quantifications performed as described in the current protocols can be used for quantitative mathematical modelling of the underlying mechanisms. Furthermore, protocols are presented for both human and pig livers. The technique is also applicable for both vibratome blocks and conventional paraffin slices