Quantifying fat zonation in liver lobules: an integrated multiscale in silico model combining disturbed microperfusion and fat metabolism via a continuum biomechanical bi-scale, tri-phasic approach

Metabolic zonation refers to the spatial separation of metabolic functions along the sinusoidal axes of the liver. This phenomenon forms the foundation for adjusting hepatic metabolism to physiological requirements in health and disease (e.g., metabolic dysfunction-associated steatotic liver disease/MASLD). Zonated metabolic functions are influenced by zonal morphological abnormalities in the liver, such as periportal fibrosis and pericentral steatosis. We aim to analyze the interplay between microperfusion, oxygen gradient, fat metabolism and resulting zonated fat accumulation in a liver lobule. Therefore we developed a continuum biomechanical, tri-phasic, bi-scale, and multicomponent in silico model, which allows to numerically simulate coupled perfusion-function-growth interactions two-dimensionally in liver lobules. The developed homogenized model has the following specifications: (i) thermodynamically consistent, (ii) tri-phase model (tissue, fat, blood), (iii) penta-substances (glycogen, glucose, lactate, FFA, and oxygen), and (iv) bi-scale approach (lobule, cell). Our presented in silico model accounts for the mutual coupling between spatial and time-dependent liver perfusion, metabolic pathways and fat accumulation. The model thus allows the prediction of fat development in the liver lobule, depending on perfusion, oxygen and plasma concentration of free fatty acids (FFA), oxidative processes, the synthesis and the secretion of triglycerides (TGs). The use of a bi-scale approach allows in addition to focus on scale bridging processes. Thus, we will investigate how changes at the cellular scale affect perfusion at the lobular scale and vice versa. This allows to predict the zonation of fat distribution (periportal or pericentral) depending on initial conditions, as well as external and internal boundary value conditions.

From Corrosion Casting to Virtual Dissection: Contrast-Enhanced Vascular Imaging using Hafnium Oxide Nanocrystals

Vascular corrosion casting is a method used to visualize the three dimensional (3D) anatomy and branching pattern of blood vessels. A polymer resin is injected in the vascular system and, after curing, the surrounding tissue is removed. The latter often deforms or even fractures the fragile cast. Here, a method is proposed that does not require corrosion, and is based on in situ micro computed tomography (micro-CT) scans. To overcome the lack of CT contrast between the polymer cast and the animals' surrounding soft tissue, hafnium oxide nanocrystals (HfO2 NCs) are introduced as CT contrast agents into the resin. The NCs dramatically improve the overall CT contrast of the cast and allow for straightforward segmentation in the CT scans. Careful design of the NC surface chemistry ensures the colloidal stability of the NCs in the casting resin. Using only 5 m% of HfO2 NCs, high-quality cardiovascular casts of both zebrafish and mice can be automatically segmented using CT imaging software. This allows to differentiate even μ $\umu$ m-scale details without having to alter the current resin injection methods. This new method of virtual dissection by visualizing casts in situ using contrast-enhanced CT imaging greatly expands the application potential of the technique.

On a deterministic mathematical model which efficiently predicts the protective effect of a plant extract mixture in cirrhotic rats

In this work, we propose a mathematical model that describes liver evolution and concentrations of alanine aminotransferase and aspartate aminotransferase in a group of rats damaged with carbon tetrachloride. Carbon tetrachloride was employed to induce cirrhosis. A second groups damaged with carbon tetrachloride was exposed simultaneously a plant extract as hepatoprotective agent. The model reproduces the data obtained in the experiment reported in [Rev. Cub. Plant. Med. 22(1), 2017], and predicts that using the plants extract helps to get a better natural recovery after the treatment. Computer simulations show that the extract reduces the damage velocity but does not avoid it entirely. The present paper is the first report in the literature in which a mathematical model reliably predicts the protective effect of a plant extract mixture in rats with cirrhosis disease. The results reported in this manuscript could be used in the future to help in fighting cirrhotic conditions in humans, though more experimental and mathematical work is required in that case.

Volumetric Tissue Imaging of Surgical Tissue Specimens Using Micro-Computed Tomography: An Emerging Digital Pathology Modality for Nondestructive, Slide-Free Microscopy-Clinical Applications of Digital Pathology in 3 Dimensions

OBJECTIVES

Micro-computed tomography (micro-CT) is a novel, nondestructive, slide-free digital imaging modality that enables the acquisition of high-resolution, volumetric images of intact surgical tissue specimens. The aim of this systematic mapping review is to provide a comprehensive overview of the available literature on clinical applications of micro-CT tissue imaging and to assess its relevance and readiness for pathology practice.

METHODS

A computerized literature search was performed in the PubMed, Scopus, Web of Science, and CENTRAL databases. To gain insight into regulatory and financial considerations for performing and examining micro-CT imaging procedures in a clinical setting, additional searches were performed in medical device databases.

RESULTS

Our search identified 141 scientific articles published between 2000 and 2021 that described clinical applications of micro-CT tissue imaging. The number of relevant publications is progressively increasing, with the specialties of pulmonology, cardiology, otolaryngology, and oncology being most commonly concerned. The included studies were mostly performed in pathology departments. Current micro-CT devices have already been cleared for clinical use, and a Current Procedural Terminology (CPT) code exists for reimbursement of micro-CT imaging procedures.

CONCLUSIONS

Micro-CT tissue imaging enables accurate volumetric measurements and evaluations of entire surgical specimens at microscopic resolution across a wide range of clinical applications.

Predicting the risk of postsplenectomy thrombosis in patients with portal hypertension using computational hemodynamics models: A proof-of-concept study

BACKGROUND

The high incidence of thrombosis in the portal venous system following splenectomy (a frequently adopted surgery for treating portal hypertension in patients with splenomegaly and hypersplenism) is a critical clinical issue. The aim of this study was to address whether quantification of postsplenectomy hemodynamics has potential value for assessing the risk of postsplenectomy thrombosis.

METHODS

Computational models were constructed for three portal hypertensive patients treated with splenectomy based on their preoperative clinical data to quantify hemodynamics in the portal venous system before and after splenectomy, respectively. Each patient was followed up for three or five months after surgery and examined with CT to screen potential thrombosis.

FINDINGS

The area ratio of wall regions exposed to low wall shear stress was small before splenectomy in all patients, which increased markedly after splenectomy and exhibited enlarged inter-patient differences. The largest area ratio of low wall shear stress and most severe flow stagnation after splenectomy were predicted for the patient suffering from postsplenectomy thrombosis, with the wall regions exposed to low wall shear stress corresponding well with the CT-detected distribution of thrombus. Further analyses revealed that postoperative hemodynamic characteristics were considerably influenced by the anatomorphological features of the portal venous system.

INTERPRETATION

Postoperative hemodynamic conditions in the portal venous system are highly patient-specific and have a potential link to postsplenectomy thrombosis, which indicates that patient-specific hemodynamic studies may serve as a complement to routine clinical assessments for refining risk stratification and postoperative patient management.

Multiscale biomechanics and mechanotransduction from liver fibrosis to cancer

A growing body of multiscale biomechanical studies has been proposed to highlight the mechanical cues in the development of hepatic fibrosis and cancer. At the cellular level, changes in mechanical microenvironment induce phenotypic and functional alterations of hepatic cells, initiating a positive feedback loop that promotes liver fibrogenesis and hepatocarcinogenesis. Tumor mechanical microenvironment of hepatocellular carcinoma facilitates tumor cell growth and metastasis, and hinders the drug delivery and immunotherapy. At the molecular level, mechanical forces are sensed and transmitted into hepatic cells via allosteric activation of mechanoreceptors on the cell membrane, leading to the activation of various mechanotransduction pathways including integrin and YAP signaling and then regulating cell function. Thus, the application of mechanomedicine concept in the treatment of liver diseases is promising for rational design and cell-specific delivery of therapeutic drugs. This review mainly discusses the correlation between biomechanical cues and liver diseases from the viewpoint of mechanobiology.

Application of Magnetic Resonance Imaging in Liver Biomechanics: A Systematic Review

MRI-based biomechanical studies can provide a deep understanding of the mechanisms governing liver function, its mechanical performance but also liver diseases. In addition, comprehensive modeling of the liver can help improve liver disease treatment. Furthermore, such studies demonstrate the beginning of an engineering-level approach to how the liver disease affects material properties and liver function. Aimed at researchers in the field of MRI-based liver simulation, research articles pertinent to MRI-based liver modeling were identified, reviewed, and summarized systematically. Various MRI applications for liver biomechanics are highlighted, and the limitations of different viscoelastic models used in magnetic resonance elastography are addressed. The clinical application of the simulations and the diseases studied are also discussed. Based on the developed questionnaire, the papers' quality was assessed, and of the 46 reviewed papers, 32 papers were determined to be of high-quality. Due to the lack of the suitable material models for different liver diseases studied by magnetic resonance elastography, researchers may consider the effect of liver diseases on constitutive models. In the future, research groups may incorporate various aspects of machine learning (ML) into constitutive models and MRI data extraction to further refine the study methodology. Moreover, researchers should strive for further reproducibility and rigorous model validation and verification.

Assessment of Vascular Network Connectivity of Hepatocellular Carcinoma Using Graph-Based Approach

BACKGROUND

The angiogenesis of liver cancer is a key condition for its growth, invasion, and metastasis. This study aims to investigate vascular network connectivity of hepatocellular carcinoma (HCC) using graph-based approach.

METHODS

Orthotopic HCC xenograft models (n=10) and the healthy controls (n=10) were established. After 21 days of modeling, hepatic vascular casting and Micro-CT scanning were performed for angiography, followed by blood vessels automatic segmentation and vascular network modeling. The topologic parameters of vascular network, including clustering coefficient (CC), network structure entropy (NSE), and average path length (APL) were quantified. Topologic parameters of the tumor region, as well as the background liver were compared between HCC group and normal control group.

RESULTS

Compared with normal control group, the tumor region of HCC group showed significantly decreased CC [(0.046 ± 0.005) vs. (0.052 ± 0.006), P=0.026], and NSE [(0.9894 ± 0.0015) vs. (0.9927 ± 0.0010), P<0.001], and increased APL [(0.433 ± 0.138) vs. (0.188 ± 0.049), P<0.001]. Compared with normal control group, the background liver of HCC group showed significantly decreased CC [(0.047 ± 0.004) vs. (0.052 ± 0.006), P=0.041] and increased NSE [0.9938 (0.9936~0.9940) vs. (0.9927 ± 0.0010), P=0.035]. No significant difference was identified for APL between the two groups.

CONCLUSION

Graph-based approach allows quantification of vascular connectivity of HCC. Disrupted vascular topological connectivity exists in the tumor region, as well as the background liver of HCC.

Influences of Anatomorphological Features of the Portal Venous System on Postsplenectomy Hemodynamic Characteristics in Patients With Portal Hypertension: A Computational Model-Based Study

Splenectomy, as an effective surgery for relieving complications caused by portal hypertension, is often accompanied by a significantly increased incidence of postoperative thrombosis in the portal venous system (PVS). While the underlying mechanisms remain insufficiently understood, the marked changes in hemodynamic conditions in the PVS following splenectomy have been suggested to be a potential contributing factor. The aim of this study was to investigate the influences of the anatomorphological features of the PVS on hemodynamic characteristics before and after splenectomy, with emphasis on identifying the specific anatomorphological features that make postoperative hemodynamic conditions more clot-promoting. For this purpose, idealized computational hemodynamics models of the PVS were constructed based on general anatomical structures and population-averaged geometrical parameters of the PVS. In the models, we incorporated various anatomorphological variations to represent inter-patient variability. The analyses of hemodynamic data were focused on the spatial distribution of wall shear stress (WSS) and the area ratio of wall regions exposed to low WSS (ALS). Obtained results showed that preoperative hemodynamic conditions were comparable among different models in terms of space-averaged WSS and ALS (all were small) irrespective of the considerable differences in spatial distribution of WSS, whereas, the inter-model differences in ALS were significantly augmented after splenectomy, with the value of ALS reaching up to over 30% in some models, while being smaller than 15% in some other models. Postoperative ALS was mainly determined by the anatomical structure of the PVS, followed by some morphogeometrical parameters, such as the diameter and curvature of the splenic vein, and the distance between the inferior mesenteric vein and splenoportal junction. Relatively, the angles between tributary veins and trunk veins only had mild influences on ALS. In addition, a marked increase in blood viscosity was predicted after splenectomy, especially in regions with low WSS, which may play an additive role to low WSS in initiating thrombosis. These findings suggest that the anatomical structure and some morphogeometrical features of the PVS are important determinants of hemodynamic conditions following splenectomy, which may provide useful clues to assessing the risk of postsplenectomy thrombosis based on medical imaging data.

Doppler ultrasonography and contrast-enhanced ultrasonography to evaluate liver allograft discard: A pilot prospective study

BACKGROUND

Broad hemodynamic changes, is believed to have a profoundly damaging effect on donor livers after brain death (DBD) or cardiac death (DCD). It remains unclear whether Doppler ultrasonography (DUS) and contrast-enhanced ultrasonography (CEUS), the imaging modalities to evaluate perfusion, could provide more information of liver discarded.

OBJECTIVE

To evaluate the ability of DUS and CEUS to predict the risk of DBD or DCD liver discarded.

METHODS

The consecutive DBD or DCD donors with DUS/CEUS examinations before surgical procurement from February 2016 to June 2018 at our institution were included. The US and CEUS images of each donor liver were analyzed and the parameters were recorded.

RESULTS

Among the 67 eligible donor livers, 15 (22.4%) were discarded and 52 (77.6%) were used. The discarded livers showed prolonged SAT of hepatic artery (0.08s vs 0.06s, OR = 2.169, P = 0.008) on DUS, less cases with homogeneous enhancement (40.0% vs 73.1%, OR = 0.243, P = 0.028) on CEUS, more cases with decreased enhancement (53.3% vs 19.2%, OR = 4.800, P = 0.009), and less difference of the peak time between portal vein and liver parenchymal (0.5s vs 6.7s, OR = 0.917, P = 0.034). The multivariable analysis showed that donor liver with prolonged SAT of hepatic artery (OR = 7.304, 95% CI: 1.195-44.655, P = 0.031) and decreased enhancement (OR = 2.588, 95% CI: 1.234-5.426, P = 0.012) were independent factors of liver discarded.

CONCLUSIONS

DUS/CEUS could be applied as a promising predictive tool to screen high-risk liver donors. The prolonged SAT of hepatic artery on DUS and the decrease of liver donor in enhancement on CEUS, indicating hemodynamic changes in DBD and DCD donor livers, were risk factors of liver discarded.

Calibration of an Electrical Analog Model of Liver Hemodynamics in Fontan Patients

Fontan associated liver disease is a common complication in patients with Fontan circulation, who were born with a single functioning heart ventricle. The hepatic venous pressure gradient (HVPG) is used to assess liver health and is a surrogate measure of the pressure gradient across the entire liver (portal pressure gradient (PPG)). However, it is thought to be inaccurate in Fontan patients. The main objectives of this study were (1) to apply an existing detailed lumped parameter model (LPM) of the liver to Fontan patients using patient-specific clinical data and (2) to determine whether HVPG is a suitable measurement of PPGs in these patients. An existing LPM of the liver blood circulation was applied and tuned to simulate patient-specific liver hemodynamics. Geometries were collected from seven adult Fontan patients and used to evaluate model parameters. The model was solved and tuned using waveform measurements of flows, inlet and outlet pressures. The predicted ratio of portal to hepatic venous pressures is comparable to in vivo measurements. The results confirmed that HVPG is not suitable for Fontan patients, as it would underestimate the portal pressures gradient by a factor of 3 to 4. Our patient-specific liver model provides an estimate of the pressure drop across the liver, which differs from the clinically used metric HVPG. This work represents a first step toward models suitable to assess liver health in Fontan patients and improve its long-term management.

Transarterial drug delivery for liver cancer: numerical simulations and experimental validation of particle distribution in patient-specific livers

Background: Transarterial therapies are routinely used for the locoregional treatment of unresectable hepatocellular carcinoma (HCC). However, the impact of clinical parameters (i.e. injection location, particle size, particle density etc.) and patient-specific conditions (i.e. hepatic geometry, cancer burden) on the intrahepatic particle distribution (PD) after transarterial injection of embolizing microparticles is still unclear. Computational fluid dynamics (CFD) may help to better understand this impact.Methods: Using CFD, both the blood flow and microparticle mass transport were modeled throughout the 3D-reconstructed arterial vasculature of a patient-specific healthy and cirrhotic liver. An experimental feasibility study was performed to simulate the PD in a 3D-printed phantom of the cirrhotic arterial network.Results: Axial and in-plane injection locations were shown to be effective parameters to steer particles toward tumor tissue in both geometries. Increasing particle size or density made it more difficult for particles to exit the domain. As cancer burden increased, the catheter tip location mattered less. The in vitro study and numerical results confirmed that PD largely mimics flow distribution, but that significant differences are still possible.Conclusions: Our findings highlight that optimal parameter choice can lead to selective targeting of tumor tissue, but that targeting potential highly depends on patient-specific conditions.

Quantification of Hemodynamic Changes in Chronic Liver Disease: Correlation of Perfusion-CT Data with Histopathologic Staging of Fibrosis

RATIONALE AND OBJECTIVES

To noninvasively estimate the severity of liver fibrosis using perfusion-CT (PCT)-based quantification of dual liver blood supply prior to liver transplantation or liver resections and to correlate results with histological grading of fibrosis stages and AST-platelet ratio index.

MATERIALS AND METHODS

Institutional review board approved this retrospective study. We analysed 41 consecutive patients (19 classified as Child-Pugh A, 17 as Child-Pugh B, and 5 as Child-Pugh C; MELD score ranged from 7 to 28) who underwent PCT prior to liver transplantation/liver resections between 2013 and 2016. The examination protocol included a scan time of 40 s, 80 kV, 100/120 mAs. Arterial liver perfusion, portal-venous perfusion and hepatic perfusion index (HPI) were registered in liver parenchyma by three readers. Fibrosis was histological graded according to Ishak scoring system as liver fibrosis (F3, n = 10), incomplete liver cirrhosis (F5, n = 5), and complete liver cirrhosis (F6, n = 26).

RESULTS

Portal-venous perfusion was significantly higher in liver fibrosis (F3 69.5±23.7 ml/100 ml/min) compared to incomplete liver cirrhosis (F5, 52.9±25.7 ml/100 ml/min) and complete liver cirrhosis (F6, 46.4±24.8 ml/100 ml/min (range 6.3-112.0 ml/100 ml/min; F = 15, p < 0.0001). HPI showed the same group differences (F = 20, p < 0.0001; HPI F3: 19.1±10.7%, HPI F5: 38.5±24.3%, HPI F6: 43.4±25.8%). Group comparisons were not significant for arterial liver perfusion (F = 3, p = 0.15). PCT parameters as well as histological fibrosis grading did neither correlate with laboratory findings including AST-platelet ratio index and MELD-Score, nor with Child-Pugh-Score.

CONCLUSION

Quantitative data from perfusion-CT can be used to differentiate between liver fibrosis (F3) and liver cirrhosis (F5/F6).

Corrosion casting in anatomy: Visualizing the architecture of hollow structures and surface details

Corrosion casting is the technique by which a solid, negative replica is created from a hollow anatomical structure and liberated from its surrounding tissues. For centuries, different types of hardening substances have been developed to create such casts, but nowadays, thermosetting polymers are mostly used as casting medium. Although the principle and initial set-up are relatively easy, producing high-quality casts that serve their intended purpose can be quite challenging. This paper evaluates some of the more popular casting resins that are currently available and provides a step-by-step overview of the corrosion casting procedure, including surface casts of anatomical structures. Hurdles and pitfalls are discussed, along with possible solutions to circumvent them, based on personal experience by the authors.

Pathological Features of Mitochondrial Ultrastructure Predict Susceptibility to Post-TIPS Hepatic Encephalopathy

BACKGROUND

Post-TIPS hepatic encephalopathy (PSE) is a complex process involving numerous risk factors; the root cause is unclear, but an elevation of blood ammonia due to portosystemic shunt and metabolic disorders in hepatocytes has been proposed as an important risk factor.

AIMS

The aim of this study was to investigate the impact of pathological features of mitochondrial ultrastructure on PSE via transjugular liver biopsy at TIPS implantation.

METHODS

We evaluated the pathological damage of mitochondrial ultrastructure on recruited patients by the Flameng classification system. A score ≤2 (no or low damage) was defined as group A, and a score >2 (high damage level) was defined as group B; routine follow-up was required at 1 and 2 years; the incidence of PSE and multiple clinical data were recorded.

RESULTS

A total of 78 cases in group A and 42 in group B completed the study. The incidence of PSE after 1 and 2 years in group B (35.7% and 45.2%, respectively) was significantly higher than that in group A (16.7% and 24.4%, respectively); the 1- and 2-year OR (95% CI) were 2.778 (1.166-6.615) and 2.565 (1.155-5.696), respectively, for groups A and B. Importantly, group B had worse incidence of PSE than group A [P=0.014, hazard ratio (95%CI): 2.172 (1.190-4.678)].

CONCLUSION

Aggressive damage to mitochondrial ultrastructure in liver shunt predicts susceptibility to PSE. The registration number is NCT02540382.

Computational Modeling in Liver Surgery

The need for extended liver resection is increasing due to the growing incidence of liver tumors in aging societies. Individualized surgical planning is the key for identifying the optimal resection strategy and to minimize the risk of postoperative liver failure and tumor recurrence. Current computational tools provide virtual planning of liver resection by taking into account the spatial relationship between the tumor and the hepatic vascular trees, as well as the size of the future liver remnant. However, size and function of the liver are not necessarily equivalent. Hence, determining the future liver volume might misestimate the future liver function, especially in cases of hepatic comorbidities such as hepatic steatosis. A systems medicine approach could be applied, including biological, medical, and surgical aspects, by integrating all available anatomical and functional information of the individual patient. Such an approach holds promise for better prediction of postoperative liver function and hence improved risk assessment. This review provides an overview of mathematical models related to the liver and its function and explores their potential relevance for computational liver surgery. We first summarize key facts of hepatic anatomy, physiology, and pathology relevant for hepatic surgery, followed by a description of the computational tools currently used in liver surgical planning. Then we present selected state-of-the-art computational liver models potentially useful to support liver surgery. Finally, we discuss the main challenges that will need to be addressed when developing advanced computational planning tools in the context of liver surgery.

Partial hepatectomy hemodynamics changes: Experimental data explained by closed-loop lumped modeling

The liver function may be degraded after partial liver ablation surgery. Adverse liver hemodynamics have been shown to be associated to liver failure. The link between these hemodynamics changes and ablation size is however poorly understood. This article proposes to explain with a closed-loop lumped model the hemodynamics changes observed during twelve surgeries in pigs. The portal venous tree is modeled with a pressure-dependent variable resistor. The variables measured, before liver ablation, are used to tune the model parameters. Then, the liver partial ablation is simulated with the model and the simulated pressures and flows are compared with post-operative measurements. Fluid infusion and blood losses occur during the surgery. The closed-loop model presented accounts for these blood volume changes. Moreover, the impact of blood volume changes and the liver lobe mass estimations on the simulated variables is studied. The typical increase of portal pressure, increase of liver pressure loss, slight decrease of portal flow and major decrease in arterial flow are quantitatively captured by the model for a 75% hepatectomy. It appears that the 75% decrease in hepatic arterial flow can be explained by the resistance increase induced by the surgery, and that no hepatic arterial buffer response (HABR) mechanism is needed to account for this change. The different post-operative states, observed in experiments, are reproduced with the proposed model. Thus, an explanation for inter-subjects post-operative variability is proposed. The presented framework can easily be adapted to other species circulations and to different pathologies for clinical hepatic applications.

Stereological analysis of size and density of hepatocytes in the porcine liver

The porcine liver is frequently used as a large animal model for verification of surgical techniques, as well as experimental therapies. Often, a histological evaluation is required that include measurements of the size, nuclearity or density of hepatocytes. Our aims were to assess the mean number-weighted volume of hepatocytes, the numerical density of hepatocytes, and the fraction of binuclear hepatocytes (BnHEP) in the porcine liver, and compare the distribution of these parameters among hepatic lobes and macroscopic regions of interest (ROIs) with different positions related to the liver vasculature. Using disector and nucleator as design-based stereological methods, the morphometry of hepatocytes was quantified in seven healthy piglets. The samples were obtained from all six hepatic lobes and three ROIs (peripheral, paracaval and paraportal) within each lobe. Histological sections (thickness 16 μm) of formalin-fixed paraffin-embedded material were stained with the periodic acid-Schiff reaction to indicate the cell outlines and were assessed in a series of 3-μm-thick optical sections. The mean number-weighted volume of mononuclear hepatocytes (MnHEP) in all samples was 3670 ± 805 μm³ (mean ± SD). The mean number-weighted volume of BnHEP was 7050 ± 2550 μm³ . The fraction of BnHEP was 4 ± 2%. The numerical density of all hepatocytes was 146 997 ± 15 738 cells mm^-3 of liver parenchyma. The porcine hepatic lobes contained hepatocytes of a comparable size, nuclearity and density. No significant differences were identified between the lobes. The peripheral ROIs of the hepatic lobes contained the largest MnHEP with the smallest numerical density. The distribution of a larger MnHEP was correlated with a larger volume of BnHEP and a smaller numerical density of all hepatocytes. Practical recommendations for designing studies that involve stereological evaluations of the size, nuclearity and density of hepatocytes in porcine liver are provided.

A multilevel framework to reconstruct anatomical 3D models of the hepatic vasculature in rat livers

The intricate (micro)vascular architecture of the liver has not yet been fully unravelled. Although current models are often idealized simplifications of the complex anatomical reality, correct morphological information is instrumental for scientific and clinical purposes. Previously, both vascular corrosion casting (VCC) and immunohistochemistry (IHC) have been separately used to study the hepatic vasculature. Nevertheless, these techniques still face a number of challenges such as dual casting in VCC and limited imaging depths for IHC. We have optimized both techniques and combined their complementary strengths to develop a framework for multilevel reconstruction of the hepatic circulation in the rat. The VCC and micro-CT scanning protocol was improved by enabling dual casting, optimizing the contrast agent concentration, and adjusting the viscosity of the resin (PU4ii). IHC was improved with an optimized clearing technique (CUBIC) that extended the imaging depth for confocal microscopy more than five-fold. Using in-house developed software (DeLiver), the vascular network - in both VCC and IHC datasets - was automatically segmented and/or morphologically analysed. Our methodological framework allows 3D reconstruction and quantification of the hepatic circulation, ranging from the major blood vessels down to the intertwined and interconnected sinusoids. We believe that the presented framework will have value beyond studies of the liver, and will facilitate a better understanding of various parenchymal organs in general, in physiological and pathological circumstances.