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

LiverMap pipeline for 3D imaging of human liver reveals volumetric spatial dysregulation of cirrhotic vasculobiliary architecture

The liver contains an intricate microstructure that is critical for liver function. Architectural disruption of this spatial structure is pathologic. Unfortunately, 2D histopathology - the gold standard for pathological understanding of many liver diseases - can misrepresent or leave gaps in our understanding of complex 3D structural features. Here, we utilized immunostaining, tissue clearing, microscopy, and computational software to create 3D multilobular reconstructions of both non-fibrotic and cirrhotic human liver tissue. We found that spatial architecture in human cirrhotic liver samples with varying etiologies had sinusoid zonation dysregulation, reduction in glutamine synthetase-expressing pericentral hepatocytes, regression of central vein networks, disruption of hepatic arterial networks, and fragmentation of biliary networks, which together suggest a pro-portalization/decentralization phenotype in cirrhotic tissue. Further implementation of 3D pathological analyses may provide a deeper understanding of cirrhotic pathobiology and inspire novel treatments for liver disease.

In Vivo Tracking and 3D Mapping of Cell Death in Regeneration and Cancer Using Trypan Blue

Tracking cell death in vivo can enable a better understanding of the biological mechanisms underlying tissue homeostasis and disease. Unfortunately, existing cell death labeling methods lack compatibility with in vivo applications or suffer from low sensitivity, poor tissue penetration, and limited temporal resolution. Here, we fluorescently labeled dead cells in vivo with Trypan Blue (TBlue) to detect single scattered dead cells or to generate whole-mount three-dimensional maps of large areas of necrotic tissue during organ regeneration. TBlue effectively marked different types of cell death, including necrosis induced by CCl4 intoxication in the liver, necrosis caused by ischemia-reperfusion in the skin, and apoptosis triggered by BAX overexpression in hepatocytes. Moreover, due to its short circulating lifespan in blood, TBlue labeling allowed in vivo "pulse and chase" tracking of two temporally spaced populations of dying hepatocytes in regenerating mouse livers. Additionally, upon treatment with cisplatin, TBlue labeled dead cancer cells in livers with cholangiocarcinoma and dead thymocytes due to chemotherapy-induced toxicity, showcasing its utility in assessing anticancer therapies in preclinical models. Thus, TBlue is a sensitive and selective cell death marker for in vivo applications, facilitating the understanding of the fundamental role of cell death in normal biological processes and its implications in disease.

Multiscale transport and 4D time-lapse imaging in precision-cut liver slices (PCLS)

Background

Monitoring cellular processes across different levels of complexity, from the cellular to the tissue scale, is important for understanding tissue structure and function. However, it is challenging to monitor and estimate these structural and dynamic interactions within three-dimensional (3D) tissue models.

Objective

The aim of this study was to design a method for imaging, tracking, and quantifying 3D changes in cell morphology (shape and size) within liver tissue, specifically a precision-cut liver slice (PCLS). A PCLS is a 3D model of the liver that allows the study of the structure and function of liver cells in their native microenvironment.

Methods

Here, we present a method for imaging liver tissue during anisosmotic exposure in a multispectral four-dimensional manner. Three metrics of tissue morphology were measured to quantify the effects of osmotic stress on liver tissue. We estimated the changes in the volume of whole precision cut liver slices, quantified the changes in nuclei position, and calculated the changes in volumetric responses of tissue-embedded cells.

Results

During equilibration with cell-membrane-permeating and non-permeating solutes, the whole tissue experiences shrinkage and expansion. As nuclei showed a change in position and directional displacement under osmotic stress, we demonstrate that nuclei could be used as a probe to measure local osmotic and mechanical stress. Moreover, we demonstrate that cells change their volume within tissue slices as a result of osmotic perturbation and that this change in volume is dependent on the position of the cell within the tissue and the duration of the exposure.

Conclusion

The results of this study have implications for a better understanding of multiscale transport, mechanobiology, and triggered biological responses within complex biological structures.

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.

A framework for ultra-low-input spatial tissue proteomics

Spatial proteomics combining microscopy-based cell phenotyping with ultrasensitive mass-spectrometry-based proteomics is an emerging and powerful concept to study cell function and heterogeneity in (patho)physiology. However, optimized workflows that preserve morphological information for phenotype discovery and maximize proteome coverage of few or even single cells from laser microdissected tissue are currently lacking. Here, we report a robust and scalable workflow for the proteomic analysis of ultra-low-input archival material. Benchmarking in murine liver resulted in up to 2,000 quantified proteins from single hepatocyte contours and nearly 5,000 proteins from 50-cell regions. Applied to human tonsil, we profiled 146 microregions including T and B lymphocyte niches and quantified cell-type-specific markers, cytokines, and transcription factors. These data also highlighted proteome dynamics within activated germinal centers, illuminating sites undergoing B cell proliferation and somatic hypermutation. This approach has broad implications in biomedicine, including early disease profiling and drug target and biomarker discovery. A record of this paper's transparent peer review process is included in the supplemental information.

Periphery of porcine hepatic lobes has the smallest length density of hepatic sinusoids and bile canaliculi: A stereological histological study with implications for liver biopsies

BACKGROUND

Porcine liver is widely used in hepatologic research as a large animal model with many anatomical and physiological similarities with humans. However, only limited information on porcine liver spatial microstructure has been published, especially regarding the hepatic sinusoids and bile canaliculi. The aim of our study was to quantify the sinusoidal and bile canalicular network in healthy male and female porcine livers and to map the variability of these structures with heterogenous distribution to improve the evaluability of liver biopsy samples.

METHODS

Livers from 12 healthy piglets (6 females and 6 neutered males) were sampled into 36 tissue samples per organ, representing six hepatic lobes and three different regions related to the hepatic vasculature (peripheral, paracaval and paraportal region). Histological sections were processed with a random orientation of the cutting plane. The endothelium and the bile canaliculi were stained using Ricinus communis agglutinin I lectin histochemistry. The length densities of hepatic sinusoids LV(sinusoids,liver), of bile canaliculi LV(bile canaliculi,liver) and volume fraction VV(sinusoids,liver) and surface density SV(sinusoids,liver) of sinusoids were estimated using stereological methods. The newly acquired morphometric data were compared with previously published data on density of porcine hepatocytes and fractions of connective tissue.

RESULTS

The peripheral region had smallest LV(sinusoids,liver), smallest LV(bile canaliculi,liver) and greatest VV(sinusoids,liver). The six hepatic lobes had statistically comparable length densities of both sinusoids and bile canaliculi, but the left lateral lobe had smallest VV(sinusoids,liver). Regions with greater LV(sinusoids,liver) had also greater LV(bile canaliculi,liver) and SV(sinusoids,liver) and were accompanied by greater density of smaller hepatocytes. Regions with smaller LV(sinusoids,liver) and LV(bile canaliculi,liver) contained a greater fraction of interlobular connective tissue.

CONCLUSIONS

The length density of hepatic sinusoids is smaller in the peripheral regions of the porcine liver than in other regions related to the hepatic vasculature - paracaval and paraportal regions, and smaller in castrated males than in females. Greater length density of liver sinusoids was linked with greater local density of bile canaliculi, with local increase in the density of smaller hepatocytes and, simultaneously, with smaller fractions of hepatic connective tissue. The intrahepatic and inter-sexual variability of the porcine liver morphology needs to be taken into account when designing and interpreting experiments involving the histological quantification of the microvascular network. The complete primary morphometric data describing the distribution of morphometric parameters within porcine liver were made available in a form facilitating the power analysis to justify the minimal number of tissue samples or animals required when designing further histological evaluation studies. The macroscopic map of microvessels and bile canaliculi variability facilitates their assessment in liver biopsies in the pig.

Microneedle array facilitates hepatic sinusoid construction in a large-scale liver-acinus-chip microsystem

Hepatic sinusoids play a key role in maintaining high activities of liver cells in the hepatic acinus. However, the construction of hepatic sinusoids has always been a challenge for liver chips, especially for large-scale liver microsystems. Herein, we report an approach for the construction of hepatic sinusoids. In this approach, hepatic sinusoids are formed by demolding a self-developed microneedle array from a photocurable cell-loaded matrix in a large-scale liver-acinus-chip microsystem with a designed dual blood supply. Primary sinusoids formed by demolded microneedles and spontaneously self-organized secondary sinusoids can be clearly observed. Benefiting from significantly enhanced interstitial flows by formed hepatic sinusoids, cell viability is witnessed to be considerably high, liver microstructure formation occurs, and hepatocyte metabolism is enhanced. In addition, this study preliminarily demonstrates the effects of the resulting oxygen and glucose gradients on hepatocyte functions and the application of the chip in drug testing. This work paves the way for the biofabrication of fully functionalized large-scale liver bioreactors.

Bioengineering liver microtissues for modeling non-alcoholic fatty liver disease

Non-alcoholic fatty liver disease (NAFLD) has become the world's most common chronic liver disease. However, due to the lack of reliable in vitro NAFLD models, drug development studies have faced many limitations, and there is no food and drug administration-approved medicine for NAFLD treatment. A functional biomimetic in vitro human liver model requires an optimized natural microenvironment using appropriate cellular composition, to provide constructive cell-cell interactions, and niche-specific bio-molecules to supply crucial cues as cell-matrix interplay. Such a suitable liver model could employ appropriate and desired biochemical, mechanical, and physical properties similar to native tissue. Moreover, bioengineered three-dimensional tissues, specially microtissues and organoids, and more recently using infusion-based cultivation systems such as microfluidics can mimic natural tissue conditions and facilitate the exchange of nutrients and soluble factors to improve physiological function in the in vitro generated constructs. This review highlights the key players involved in NAFLD initiation and progression and discussed the available cells and matrices for in vitro NAFLD modeling. The strategies for optimizing the liver microenvironment to generate a powerful and biomimetic in vitro NAFLD model were described as well. Finally, the current challenges and future perospective for promotion in this subject were discussed.

In silico modeling for the hepatic circulation and transport: From the liver organ to lobules

The function of the liver depends critically on its blood supply. Numerous in silico models have been developed to study various aspects of the hepatic circulation, including not only the macro-hemodynamics at the organ level, but also the microcirculation at the lobular level. In addition, computational models of blood flow and bile flow have been used to study the transport, metabolism, and clearance of drugs in pharmacokinetic studies. These in silico models aim to provide insights into the liver organ function under both healthy and diseased states, and to assist quantitative analysis for surgical planning and postsurgery treatment. The purpose of this review is to provide an update on state-of-the-art in silico models of the hepatic circulation and transport processes. We introduce the numerical methods and the physiological background of these models. We also discuss multiscale frameworks that have been proposed for the liver, and their linkage with the large context of systems biology, systems pharmacology, and the Physiome project. This article is categorized under: Metabolic Diseases > Computational Models Metabolic Diseases > Biomedical Engineering Cardiovascular Diseases > Computational Models.

In vitro to in vivo acetaminophen hepatotoxicity extrapolation using classical schemes, pharmacodynamic models and a multiscale spatial-temporal liver twin

In vitro to in vivo extrapolation represents a critical challenge in toxicology. In this paper we explore extrapolation strategies for acetaminophen (APAP) based on mechanistic models, comparing classical (CL) homogeneous compartment pharmacodynamic (PD) models and a spatial-temporal (ST), multiscale digital twin model resolving liver microarchitecture at cellular resolution. The models integrate consensus detoxification reactions in each individual hepatocyte. We study the consequences of the two model types on the extrapolation and show in which cases these models perform better than the classical extrapolation strategy that is based either on the maximal drug concentration (Cmax) or the area under the pharmacokinetic curve (AUC) of the drug blood concentration. We find that an CL-model based on a well-mixed blood compartment is sufficient to correctly predict the in vivo toxicity from in vitro data. However, the ST-model that integrates more experimental information requires a change of at least one parameter to obtain the same prediction, indicating that spatial compartmentalization may indeed be an important factor.

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

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.

Role of the Hippo pathway in liver regeneration and repair: recent advances

Although the signaling pathways involved in normal liver regeneration have been well characterized, less has been done for livers affected by chronic tissue damage. These "abnormal livers" have an impaired regenerative response that leads to liver repair and fibrosis. The tumor suppressor Hippo pathway plays a key role in liver regeneration and repair. On this basis, this review discusses recent studies focusing on the involvement of the Hippo signaling pathway during "normal healthy liver regeneration" (i.e., in a normal liver after 2/3 partial hepatectomy) and "abnormal liver regeneration" (i.e., in a liver damaged by chronic disease). This could be an important question to address with respect to new therapies aimed at improving impaired liver regenerative responses. The studies reported here have shown that activation of the Hippo coactivators YAP/TAZ during normal liver regeneration promotes the formation of a new bile duct network through direct BEC proliferation or/and hepatocyte dedifferentiation to HPCs which can trans-differentiate to BECs. Moreover, YAP/TAZ signaling interaction with other signaling pathways mediates the recruitment and activation of Kupffer cells, which release mitogenic cytokines for parenchymal and/or non-parenchymal cells and engage in phagocytosis of cellular debris. In addition, YAP-mediated activation of stellate cells (HSCs) promotes liver regeneration through the synthesis of extracellular matrix. However, in chronically diseased livers, where the predetermined threshold for proper liver regeneration is exceeded, YAP/TAZ activation results in a reparative process characterized by liver fibrosis. In this condition, YAP/TAZ activation in parenchymal and non-parenchymal cells results in (i) differentiation of quiescent HSCs into myofibroblastic HSCs; (ii) recruitment of macrophages releasing inflammatory cytokines; (iii) polarization of macrophages toward the M2 phenotype. Since accumulation of damaged hepatocytes in chronic liver injury represent a significant risk factor for the development of hepatocarcinoma, this review also discussed the involvement of the Hippo pathway in the clearance of damaged cells.

From Seeing to Simulating: A Survey of Imaging Techniques and Spatially-Resolved Data for Developing Multiscale Computational Models of Liver Regeneration

Liver regeneration, which leads to the re-establishment of organ mass, follows a specifically organized set of biological processes acting on various time and length scales. Computational models of liver regeneration largely focused on incorporating molecular and signaling detail have been developed by multiple research groups in the recent years. These modeling efforts have supported a synthesis of disparate experimental results at the molecular scale. Incorporation of tissue and organ scale data using noninvasive imaging methods can extend these computational models towards a comprehensive accounting of multiscale dynamics of liver regeneration. For instance, microscopy-based imaging methods provide detailed histological information at the tissue and cellular scales. Noninvasive imaging methods such as ultrasound, computed tomography and magnetic resonance imaging provide morphological and physiological features including volumetric measures over time. In this review, we discuss multiple imaging modalities capable of informing computational models of liver regeneration at the organ-, tissue- and cellular level. Additionally, we discuss available software and algorithms, which aid in the analysis and integration of imaging data into computational models. Such models can be generated or tuned for an individual patient with liver disease. Progress towards integrated multiscale models of liver regeneration can aid in prognostic tool development for treating liver disease.

Insight into Bile Duct Reaction to Obstruction from a Three-dimensional Perspective Using ex Vivo Phase-Contrast CT

Background Biliary obstruction leads to an increase in biliary pressure within the biliary system, which induces the morphologic adaptation of the biliary tree. Purpose To observe and to quantify the morphologic characteristics of the adaptation in a bile duct ligation rat model and verify it in patients with biliary atresia in a three-dimensional (3D) manner using x-ray phase-contrast CT. Materials and Methods A bile duct ligation model was induced in 40 male Sprague-Dawley rats, which were divided into five groups: the control group (no ligation) and groups 2, 4, 6, and 8 weeks after bile duct ligation (eight animals in each group). Liver tissue samples (approximately 1.8 cm in length and 1.3 cm in height) were imaged by using phase-contrast CT and compared with histologic analysis. With a combination of phase-contrast CT and 3D visualization technology, the entire biliary system and the intrahepatic vascular system were quantitatively analyzed according to downstream, midstream, and upstream domains based on bile duct volume, surface area, and other parameters. Additionally, liver explant tissues from 28 patients with biliary atresia were studied to determine the impact of biliary tract reconstruction. Results To offset the increased biliary pressure within the biliary system, the ductular reaction in the downstream, midstream, and upstream domains manifested as dilatation, spiderweb-like looping, and interconnected honeycomb-like patterns, respectively. The most severe ductular reaction occurred in the upstream domain, and the relative surface area (mean, 0.02 μm^-1 ± 0.01, 0.04 μm^-1 ± 0.01, 0.07 μm^-1 ± 0.02, and 0.10 μm^-1 ± 0.02 for the 2-8-week groups, respectively; P < .01 among the groups) and volume fraction of ductules (mean, 16.54% ± 4.62, 19.69% ± 6.41, 26.92% ± 5.82, and 38.34% ± 10.36 for the 2-8-week groups, respectively; P < .01 among the groups except between the 2- and 4-week groups [P = .062]) significantly increased over time. In patients with biliary atresia, it was observed that both fibrosis and proliferative ductules regressed after successful biliary tract reconstruction following Kasai portoenterostomy. Furthermore, ductular reaction was accompanied by a progressive increase in the arterial supply but a loss of portal blood supply. Conclusion X-ray phase-contrast CT with three-dimensional rendering of the biliary system in a bile duct ligation rat model provides key insights into ductular reaction or biliary self-adaptation triggered by increased biliary pressure. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Vannier and Wang in this issue.

Regeneration Defects in Yap and Taz Mutant Mouse Livers Are Caused by Bile Duct Disruption and Cholestasis

BACKGROUND AND AIMS

The Hippo pathway and its downstream effectors YAP and TAZ (YAP/TAZ) are heralded as important regulators of organ growth and regeneration. However, different studies provided contradictory conclusions about their role during regeneration of different organs, ranging from promoting proliferation to inhibiting it. Here we resolve the function of YAP/TAZ during regeneration of the liver, where Hippo's role in growth control has been studied most intensely.

METHODS

We evaluated liver regeneration after carbon tetrachloride toxic liver injury in mice with conditional deletion of Yap/Taz in hepatocytes and/or biliary epithelial cells, and measured the behavior of different cell types during regeneration by histology, RNA sequencing, and flow cytometry.

RESULTS

We found that YAP/TAZ were activated in hepatocytes in response to carbon tetrachloride toxic injury. However, their targeted deletion in adult hepatocytes did not noticeably impair liver regeneration. In contrast, Yap/Taz deletion in adult bile ducts caused severe defects and delay in liver regeneration. Mechanistically, we showed that Yap/Taz mutant bile ducts degenerated, causing cholestasis, which stalled the recruitment of phagocytic macrophages and the removal of cellular corpses from injury sites. Elevated bile acids activated pregnane X receptor, which was sufficient to recapitulate the phenotype observed in mutant mice.

CONCLUSIONS

Our data show that YAP/TAZ are practically dispensable in hepatocytes for liver development and regeneration. Rather, YAP/TAZ play an indirect role in liver regeneration by preserving bile duct integrity and securing immune cell recruitment and function.

Simulation of a detoxifying organ function: Focus on hemodynamics modeling and convection-reaction numerical simulation in microcirculatory networks

When modeling a detoxifying organ function, an important component is the impact of flow on the metabolism of a compound of interest carried by the blood. We here study the effects of red blood cells (such as the Fahraeus-Lindqvist effect and plasma skimming) on blood flow in typical microcirculatory components such as tubes, bifurcations and entire networks, with particular emphasis on the liver as important representative of detoxifying organs. In one of the plasma skimming models, under certain conditions, oscillations between states are found and analyzed in a methodical study to identify their causes and influencing parameters. The flow solution obtained is then used to define the velocity at which a compound would be transported. A convection-reaction equation is studied to simulate the transport of a compound in blood and its uptake by the surrounding cells. Different types of signal sharpness have to be handled depending on the application to address different temporal compound concentration profiles. To permit executing the studied models numerically stable and accurate, we here extend existing transport schemes to handle converging bifurcations, and more generally multi-furcations. We study the accuracy of different numerical schemes as well as the effect of reactions and of the network itself on the bolus shape. Even though this study is guided by applications in liver micro-architecture, the proposed methodology is general and can readily be applied to other capillary network geometries, hence to other organs or to bioengineered network designs.

Mutual Zonated Interactions of Wnt and Hh Signaling Are Orchestrating the Metabolism of the Adult Liver in Mice and Human

The Hedgehog (Hh) and Wnt/β-Catenin (Wnt) cascades are morphogen pathways whose pronounced influence on adult liver metabolism has been identified in recent years. How both pathways communicate and control liver metabolic functions are largely unknown. Detecting core components of Wnt and Hh signaling and mathematical modeling showed that both pathways in healthy liver act largely complementary to each other in the pericentral (Wnt) and the periportal zone (Hh) and communicate mainly by mutual repression. The Wnt/Hh module inversely controls the spatiotemporal operation of various liver metabolic pathways, as revealed by transcriptome, proteome, and metabolome analyses. Shifting the balance to Wnt (activation) or Hh (inhibition) causes pericentralization and periportalization of liver functions, respectively. Thus, homeostasis of the Wnt/Hh module is essential for maintaining proper liver metabolism and to avoid the development of certain metabolic diseases. With caution due to minor species-specific differences, these conclusions may hold for human liver as well.

Highlight report: New insights in liver physiology: Canalicular bile flux is diffusion dominated

One of the central functions of the liver is excretion of bile into the intestine. Currently, bile excretion is explained by the osmotic model, according to which bile acids are excreted by hepatocytes into the bile canaliculi and since bile acids are osmotically active they draw water into the canalicular lumen. Bile canaliculi are closed at the central side. Therefore, bile was postulated to flow to the open side into the ducts. However, bile flow in canaliculi has never been measured because of the small canalicular diameter which does not allow analysis of flux by conventional methods. Recently, methods have been developed that allow flow analysis in bile canaliculi and ducts. Interestingly, no measurable directed flow was observed in the canaliculi. Instead, small molecules in bile canaliculi reached the larger bile ducts by diffusion. Only there measurable flow sets in. The pathophysiological implications of this novel observation are discussed.

Tubular bile duct structure mimicking bile duct morphogenesis for prospective in vitro liver metabolite recovery

Background

Liver metabolites are used to diagnose disease and examine drugs in clinical pharmacokinetics. Therefore, development of an in vitro assay system that reproduces liver metabolite recovery would provide important benefits to pharmaceutical research. However, liver models have proven challenging to develop because of the lack of an appropriate bile duct structure for the accumulation and transport of metabolites from the liver parenchyma. Currently available bile duct models, such as the bile duct cyst-embedded extracellular matrix (ECM), lack any morphological resemblance to the tubular morphology of the living bile duct. Moreover, these systems cannot overcome metabolite recovery issues because they are established in isolated culture systems. Here, we successfully established a non-continuous tubular bile duct structure model in an open-culture system, which closely resembled an in vivo structure. This system was utilized to effectively collect liver metabolites separately from liver parenchymal cells.

Results

Triple-cell co-culture of primary rat hepatoblasts, rat biliary epithelial cells, and mouse embryonic fibroblasts was grown to mimic the morphogenesis of the bile duct during liver development. Overlaying the cells with ECM containing a Matrigel and collagen type I gel mixture promoted the development of a tubular bile duct structure. In this culture system, the expression of specific markers and signaling molecules related to biliary epithelial cell differentiation was highly upregulated during the ductal formation process. This bile duct structure also enabled the separate accumulation of metabolite analogs from liver parenchymal cells.

Conclusions

A morphogenesis-based culture system effectively establishes an advanced bile duct structure and improves the plasticity of liver models feasible for autologous in vitro metabolite-bile collection, which may enhance the performance of high-throughput liver models in cell-based assays.

A quantitative high-resolution computational mechanics cell model for growing and regenerating tissues

Mathematical models are increasingly designed to guide experiments in biology, biotechnology, as well as to assist in medical decision making. They are in particular important to understand emergent collective cell behavior. For this purpose, the models, despite still abstractions of reality, need to be quantitative in all aspects relevant for the question of interest. This paper considers as showcase example the regeneration of liver after drug-induced depletion of hepatocytes, in which the surviving and dividing hepatocytes must squeeze in between the blood vessels of a network to refill the emerged lesions. Here, the cells' response to mechanical stress might significantly impact the regeneration process. We present a 3D high-resolution cell-based model integrating information from measurements in order to obtain a refined and quantitative understanding of the impact of cell-biomechanical effects on the closure of drug-induced lesions in liver. Our model represents each cell individually and is constructed by a discrete, physically scalable network of viscoelastic elements, capable of mimicking realistic cell deformation and supplying information at subcellular scales. The cells have the capability to migrate, grow, and divide, and the nature and parameters of their mechanical elements can be inferred from comparisons with optical stretcher experiments. Due to triangulation of the cell surface, interactions of cells with arbitrarily shaped (triangulated) structures such as blood vessels can be captured naturally. Comparing our simulations with those of so-called center-based models, in which cells have a largely rigid shape and forces are exerted between cell centers, we find that the migration forces a cell needs to exert on its environment to close a tissue lesion, is much smaller than predicted by center-based models. To stress generality of the approach, the liver simulations were complemented by monolayer and multicellular spheroid growth simulations. In summary, our model can give quantitative insight in many tissue organization processes, permits hypothesis testing in silico, and guide experiments in situations in which cell mechanics is considered important.

ECM1 Prevents Activation of Transforming Growth Factor β, Hepatic Stellate Cells, and Fibrogenesis in Mice

BACKGROUND & AIMS

Activation of TGFB (transforming growth factor β) promotes liver fibrosis by activating hepatic stellate cells (HSCs), but the mechanisms of TGFB activation are not clear. We investigated the role of ECM1 (extracellular matrix protein 1), which interacts with extracellular and structural proteins, in TGFB activation in mouse livers.

METHODS

We performed studies with C57BL/6J mice (controls), ECM1-knockout (ECM1-KO) mice, and mice with hepatocyte-specific knockout of EMC1 (ECM1Δhep). ECM1 or soluble TGFBR2 (TGFB receptor 2) were expressed in livers of mice after injection of an adeno-associated virus vector. Liver fibrosis was induced by carbon tetrachloride (CCl4) administration. Livers were collected from mice and analyzed by histology, immunohistochemistry, in situ hybridization, and immunofluorescence analyses. Hepatocytes and HSCs were isolated from livers of mice and incubated with ECM1; production of cytokines and activation of reporter genes were quantified. Liver tissues from patients with viral or alcohol-induced hepatitis (with different stages of fibrosis) and individuals with healthy livers were analyzed by immunohistochemistry and in situ hybridization.

RESULTS

ECM1-KO mice spontaneously developed liver fibrosis and died by 2 months of age without significant hepatocyte damage or inflammation. In liver tissues of mice, we found that ECM1 stabilized extracellular matrix-deposited TGFB in its inactive form by interacting with αv integrins to prevent activation of HSCs. In liver tissues from patients and in mice with CCl4-induced liver fibrosis, we found an inverse correlation between level of ECM1 and severity of fibrosis. CCl4-induced liver fibrosis was accelerated in ECM1Δhep mice compared with control mice. Hepatocytes produced the highest levels of ECM1 in livers of mice. Ectopic expression of ECM1 or soluble TGFBR2 in liver prevented fibrogenesis in ECM1-KO mice and prolonged their survival. Ectopic expression of ECM1 in liver also reduced the severity of CCl4-induced fibrosis in mice.

CONCLUSIONS

ECM1, produced by hepatocytes, inhibits activation of TGFB and its activation of HSCs to prevent fibrogenesis in mouse liver. Strategies to increase levels of ECM1 in liver might be developed for treatment of fibrosis.

Modulating Portal Hemodynamics With Vascular Ring Allows Efficient Regeneration After Partial Hepatectomy in a Porcine Model

OBJECTIVE

To investigate safety and efficacy of temporary portal hemodynamics modulation with a novel percutaneously adjustable vascular ring (MID-AVR) onto a porcine model of 75% hepatectomy.

BACKGROUND

Postoperative liver failure is a leading cause of mortality after major hepatectomy. Portal flow modulation is an increasingly accepted concept to prevent postoperative liver failure. Nonetheless, the current strategies have shortcomings.

METHODS

Resection was performed under hemodynamic monitoring in 17 large, white pigs allocated into 2 groups. Eight pigs had ring around the portal vein for 3 days with the aim of reducing changes in hemodynamics due to hepatectomy. Analysis of hemodynamics, laboratory, and histopathological parameters was performed.

RESULTS

Percutaneous inflation, deflation, and removal of the MID-AVR were safe. Two (25%) pigs in the MID-AVR group and 4 (45%) controls died before day 3 (P = NS). A moderate increase of portal flow rate per liver mass after resection was associated with better survival (P = 0.017). The portocaval pressure gradient was lower after hepatectomy in the MID-AVR group (P = 0.001). Postoperative serum bilirubin levels were lower in the MID-AVR group (P = 0.007 at day 5). In the MID-AVR group, the Ki67 index was significantly higher on day 3 (P = 0.043) and the architectural derangement was lower (P < 0.05). Morphometric quantification of the bile canaliculi revealed a significantly lower number of intersection branches (P < 0.05) and intersection nodes (P < 0.001) on day 7 compared with the preoperative specimen, in the control group. These differences were not found in the ring group.

CONCLUSIONS

MID-AVR is safe for portal hemodynamics modulation. It might improve liver regeneration by protecting liver microarchitecture.

Prediction of human drug-induced liver injury (DILI) in relation to oral doses and blood concentrations

Drug-induced liver injury (DILI) cannot be accurately predicted by animal models. In addition, currently available in vitro methods do not allow for the estimation of hepatotoxic doses or the determination of an acceptable daily intake (ADI). To overcome this limitation, an in vitro/in silico method was established that predicts the risk of human DILI in relation to oral doses and blood concentrations. This method can be used to estimate DILI risk if the maximal blood concentration (C_max) of the test compound is known. Moreover, an ADI can be estimated even for compounds without information on blood concentrations. To systematically optimize the in vitro system, two novel test performance metrics were introduced, the toxicity separation index (TSI) which quantifies how well a test differentiates between hepatotoxic and non-hepatotoxic compounds, and the toxicity estimation index (TEI) which measures how well hepatotoxic blood concentrations in vivo can be estimated. In vitro test performance was optimized for a training set of 28 compounds, based on TSI and TEI, demonstrating that (1) concentrations where cytotoxicity first becomes evident in vitro (EC₁₀) yielded better metrics than higher toxicity thresholds (EC₅₀); (2) compound incubation for 48 h was better than 24 h, with no further improvement of TSI after 7 days incubation; (3) metrics were moderately improved by adding gene expression to the test battery; (4) evaluation of pharmacokinetic parameters demonstrated that total blood compound concentrations and the 95%-population-based percentile of C_max were best suited to estimate human toxicity. With a support vector machine-based classifier, using EC₁₀ and C_max as variables, the cross-validated sensitivity, specificity and accuracy for hepatotoxicity prediction were 100, 88 and 93%, respectively. Concentrations in the culture medium allowed extrapolation to blood concentrations in vivo that are associated with a specific probability of hepatotoxicity and the corresponding oral doses were obtained by reverse modeling. Application of this in vitro/in silico method to the rat hepatotoxicant pulegone resulted in an ADI that was similar to values previously established based on animal experiments. In conclusion, the proposed method links oral doses and blood concentrations of test compounds to the probability of hepatotoxicity.

Second exposure to acetaminophen overdose is associated with liver fibrosis in mice

Repeated administration of hepatotoxicants is usually accompanied by liver fibrosis. However, the difference in response as a result of repeated exposures of acetaminophen (APAP) compared to a single dose is not well-studied. Therefore, in the current study, the liver response after a second dose of APAP was investigated. Adult fasted Balb/C mice were exposed to two toxic doses of 300 mg/kg APAP, which were administered 72 h apart from each other. Subsequently, blood and liver from the treated mice were collected 24 h and 72 h after both APAP administrations. Liver transaminase, i.e. alanine amino transferase (ALT) and aspartate amino transferase (AST) levels revealed that the fulminant liver damage was reduced after the second APAP administration compared to that observed at the same time point after the first treatment. These results correlated with the necrotic areas as indicated by histological analyses. Surprisingly, Picro Sirius Red (PSR) staining showed that the accumulation of extracellular matrix after the second dose coincides with the upregulation of some fibrogenic signatures, e.g., alpha smooth muscle actin. Non-targeted liver tissue metabolic profiling indicates that most alterations occur 24 h after the first dose of APAP. However, the levels of most metabolites recover to basal values over time. This organ adaptation process is also confirmed by the upregulation of antioxidative systems like e.g. superoxide dismutase and catalase. From the results, it can be concluded that there is a different response of the liver to APAP toxic doses, if the liver has already been exposed to APAP. A necroinflammatory process followed by a liver regeneration was observed after the first APAP exposure. However, fibrogenesis through the accumulation of extracellular matrix is observed after a second challenge. Therefore, further studies are required to mechanistically understand the so called "liver memory".

Bile Microinfarcts in Cholestasis Are Initiated by Rupture of the Apical Hepatocyte Membrane and Cause Shunting of Bile to Sinusoidal Blood

Bile duct ligation (BDL) is an experimental procedure that mimics obstructive cholestatic disease. One of the early consequences of BDL in rodents is the appearance of so-called bile infarcts that correspond to Charcot-Gombault necrosis in human cholestasis. The mechanisms causing bile infarcts and their pathophysiological relevance are unclear. Therefore, intravital two photon-based imaging of BDL mice was performed with fluorescent bile salts (BS) and non-BS organic anion analogues. Key findings were followed up by matrix-assisted laser desorption ionization imaging, clinical chemistry, immunostaining, and gene expression analyses. In the acute phase, 1-3 days after BDL, BS concentrations in bile increased and single-cell bile microinfarcts occurred in dispersed hepatocytes throughout the liver caused by the rupture of the apical hepatocyte membrane. This rupture occurred after loss of mitochondrial membrane potential, followed by entry of bile, cell death, and a "domino effect" of further death events of neighboring hepatocytes. Bile infarcts provided a trans-epithelial shunt between bile canaliculi and sinusoids by which bile constituents leaked into blood. In the chronic phase, ≥21 days after BDL, uptake of BS tracers at the sinusoidal hepatocyte membrane was reduced. This contributes to elevated concentrations of BS in blood and decreased concentrations in the biliary tract. Conclusion: Bile microinfarcts occur in the acute phase after BDL in a limited number of dispersed hepatocytes followed by larger infarcts involving neighboring hepatocytes, and they allow leakage of bile from the BS-overloaded biliary tract into blood, thereby protecting the liver from BS toxicity; in the chronic phase after BDL, reduced sinusoidal BS uptake is a dominant protective factor, and the kidney contributes to the elimination of BS until cholemic nephropathy sets in.

Pipe-3D: A Pipeline Based on Immunofluorescence, 3D Confocal Imaging, Reconstructions, and Morphometry for Biliary Network Analysis in Cholestasis

Cholestasis, the impairment of bile flux out of the liver, is a common complication of many pathological liver disorders, such as cholangiopathies, primary biliary sclerosis, and primary biliary cirrhosis. Besides accumulation of bile acids in the liver and blood, it leads to a proliferative response of the biliary tree termed as a ductular reaction. The ductular reaction is characterized by enhanced proliferation of cholangiocytes, which form the epithelial lining of bile ducts. This strong reaction of the biliary tree has been reported to generate a source of progenitor cells that can differentiate to hepatocytes or cholangiocytes during regeneration. On the other hand, it can cause periportal fibrosis eventually progressing to cirrhosis and death. In 2D histology, this leads to the appearance of an increased number of duct lumina per area of tissue. Yet, the biliary tree is a 3D vstructure and the appearance of lumina in thin slices may be explained by the appearance of novel ducts or by ramification or convolution of existing ducts in 3D. In many such aspects, traditional 2D histology on thin slices limits our understanding of the response of the biliary tree. A comprehensive understanding of architecture remodeling of the biliary network in cholestasis depends on robust 3D sample preparation and analysis methods. To that end, we describe pipe-3D, a processing and analysis pipeline visualization based on immunofluorescence, confocal imaging, surface reconstructions, and automated morphometry of the biliary network in 3D at subcellular resolution. This pipeline has been used to discover extensive remodeling of interlobular bile ducts in cholestasis, wherein elongation, branching, and looping create a dense ductular mesh around the portal vein branch. Surface reconstructions generated by Pipe-3D from confocal data also show an approximately fivefold enhancement of the luminal duct surface through corrugation of the epithelial lamina, which may increase bile reabsorption and alleviate cholestasis. The response of interlobular ducts in cholestasis was shown to be in sharp contrast to that of large bile ducts, de novo duct formation during embryogenesis. It is also distinct from ductular response in other models of hepatic injury such as choline-deficient, ethionine-supplemented diet, where parenchymal tissue invasion by ducts and their branches is observed. Pipe-3D is applicable to any model of liver injury, and optionally integrates tissue clearing techniques for 3D analysis of thick (>500 μm) tissue sections.

Modeling of xenobiotic transport and metabolism in virtual hepatic lobule models

Computational models of normal liver function and xenobiotic induced liver damage are increasingly being used to interpret in vitro and in vivo data and as an approach to the de novo prediction of the liver’s response to xenobiotics. The microdosimetry (dose at the level of individual cells) of xenobiotics vary spatially within the liver because of both compound-independent and compound-dependent factors. In this paper, we build model liver lobules to investigate the interplay between vascular structure, blood flow and cellular transport that lead to regional variations in microdosimetry. We then compared simulation results obtained using this complex spatial model with a simpler linear pipe model of a sinusoid and a very simple single box model. We found that variations in diffusive transport, transporter-mediated transport and metabolism, coupled with complex liver sinusoid architecture and blood flow distribution, led to three essential patterns of xenobiotic exposure within the virtual liver lobule: (1) lobular-wise uniform, (2) radially varying and (3) both radially and azimuthally varying. We propose to use these essential patterns of exposure as a reference for selection of model representations when a computational study involves modeling detailed hepatic responses to xenobiotics.

Shotgun Lipidomics Combined with Laser Capture Microdissection: A Tool To Analyze Histological Zones in Cryosections of Tissues

Shotgun analysis provides a quantitative snapshot of the lipidome composition of cells, tissues, or model organisms; however, it does not elucidate the spatial distribution of lipids. Here we demonstrate that shotgun analysis could quantify low-picomole amounts of lipids isolated by laser capture microdissection (LCM) of hundred micrometer-sized histological zones visualized at the cryosections of tissues. We identified metabolically distinct periportal (pp) and pericentral (pc) zones by immunostaining of 20 μm thick cryosections of a healthy mouse liver. LCM was used to ablate, catapult, and collect the tissue material from 10 to 20 individual zones covering a total area of 0.3-0.5 mm² and containing ca. 500 cells. Top-down shotgun profiling relying upon computational stitching of 61 targeted selective ion monitoring ( t-SIM) spectra quantified more than 200 lipid species from 17 lipid classes including glycero- and glycerophospholipids, sphingolipids, cholesterol esters, and cholesterol. Shotgun LCM revealed the overall commonality of the full lipidome composition of pp and pc zones along with significant ( p < 0.001) difference in the relative abundance of 13 lipid species. Follow-up proteomics analyses of pellets recovered from an aqueous phase saved after the lipid extraction identified 13 known and 7 new protein markers exclusively present in pp or in pc zones and independently validated the specificity of their visualization, isolation, and histological assignment.

Confounding influence of tamoxifen in mouse models of Cre recombinase-induced gene activity or modulation

Tamoxifen (TAM) is commonly used for cell type specific Cre recombinase-induced gene inactivation and in cell fate tracing studies. Inducing a gene knockout by TAM and using non-TAM exposed mice as controls lead to a situation where differences are interpreted as consequences of the gene knockout but in reality result from TAM-induced changes in hepatic metabolism. The degree to which TAM may compromise the interpretation of animal experiments with inducible gene expression still has to be elucidated. Here, we report that TAM strongly attenuates CCl₄-induced hepatotoxicity in male C57Bl/6N mice, even after a 10 days TAM exposure-free period. TAM decreased (p < 0.0001) the necrosis index and the level of aspartate- and alanine transaminases in CCl₄-treated compared to vehicle-exposed mice. TAM pretreatment also led to the downregulation of CYP2E1 (p = 0.0045) in mouse liver tissue, and lowered its activity in CYP2E1 expressing HepG2 cell line. Furthermore, TAM increased the level of the antioxidant ascorbate, catalase, SOD2, and methionine, as well as phase II metabolizing enzymes GSTM1 and UGT1A1 in CCl₄-treated livers. Finally, we found that TAM increased the presence of resident macrophages and recruitment of immune cells in necrotic areas of the livers as indicated by F4/80 and CD45 staining. In conclusion, we reveal that TAM increases liver resistance to CCl₄-induced toxicity. This finding is of high relevance for studies using the tamoxifen-inducible expression system particularly if this system is used in combination with hepatotoxic compounds such as CCl₄.

Galunisertib modifies the liver fibrotic composition in the Abcb4Ko mouse model

Transforming growth factor (TGF)-β stimulates extracellular matrix (ECM) deposition during development of liver fibrosis and cirrhosis, the most important risk factor for the onset of hepatocellular carcinoma. In liver cancer, TGF-β is responsible for a more aggressive and invasive phenotype, orchestrating remodeling of the tumor microenvironment and triggering epithelial-mesenchymal transition of cancer cells. This is the scientific rationale for targeting the TGF-β pathway via a small molecule, galunisertib (intracellular inhibitor of ALK5) in clinical trials to treat liver cancer patients at an advanced disease stage. In this study, the hypothesis that galunisertib modifies the tissue microenvironment via inhibition of the TGF-β pathway is tested in an experimental preclinical model. At the age of 6 months, Abcb4ko mice-a well-established model for chronic liver disease development and progression-are treated twice daily with galunisertib (150 mg/kg) via oral gavage for 14 consecutive days. Two days after the last treatment, blood plasma and livers are harvested for further assessment, including fibrosis scoring and ECM components. The reduction of Smad2 phosphorylation in both parenchymal and non-parenchymal liver cells following galunisertib administration confirms the treatment effectiveness. Damage-related galunisertib does not change cell proliferation, macrophage numbers and leucocyte recruitment. Furthermore, no clear impact on the amount of fibrosis is evident, as documented by PicroSirius red and Gomori-trichome scoring. On the other hand, several fibrogenic genes, e.g., collagens (Col1α1 and Col1α2), Tgf-β1 and Timp1, mRNA levels are significantly downregulated by galunisertib administration when compared to controls. Most interestingly, ECM/stromal components, fibronectin and laminin-332, as well as the carcinogenic β-catenin pathway, are remarkably reduced by galunisertib-treated Abcb5ko mice. In conclusion, TGF-β inhibition by galunisertib interferes, to some extent, with chronic liver progression, not by reducing the stage of liver fibrosis as measured by different scoring systems, but rather by modulating the biochemical composition of the deposited ECM, likely affecting the fate of non-parenchymal cells.

A comprehensive model for assessment of liver stage therapies targeting Plasmodium vivax and Plasmodium falciparum

Malaria liver stages represent an ideal therapeutic target with a bottleneck in parasite load and reduced clinical symptoms; however, current in vitro pre-erythrocytic (PE) models for Plasmodium vivax and P. falciparum lack the efficiency necessary for rapid identification and effective evaluation of new vaccines and drugs, especially targeting late liver-stage development and hypnozoites. Herein we report the development of a 384-well plate culture system using commercially available materials, including cryopreserved primary human hepatocytes. Hepatocyte physiology is maintained for at least 30 days and supports development of P. vivax hypnozoites and complete maturation of P. vivax and P. falciparum schizonts. Our multimodal analysis in antimalarial therapeutic research identifies important PE inhibition mechanisms: immune antibodies against sporozoite surface proteins functionally inhibit liver stage development and ion homeostasis is essential for schizont and hypnozoite viability. This model can be implemented in laboratories in disease-endemic areas to accelerate vaccine and drug discovery research.