Theranostics 2020; 10(15):6715-6727. doi:10.7150/thno.44900 This issue
1. The University of Queensland Diamantina Institute, The University of Queensland, Brisbane, QLD, 4102, Australia
2. Gallipoli Medical Research Institute, Greenslopes Private Hospital, Brisbane, QLD, 4120, Australia
3. Faculty of Medicine, The University of Queensland, Brisbane, QLD, 4072, Australia
4. Department of Biliary-pancreatic Surgery, Ren Ji Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai, 200127, China.
5. Preclinical Imaging Facility, Translational Research Institute, Brisbane, QLD, 4102, Australia
6. The University of Sydney, Concord Hospital, Concord, NSW, 2139, Australia
7. Department of General Surgery, Changzheng Hospital, The Second Military Medical University, Shanghai, 200003, China
Acetaminophen (APAP) is the foremost cause of drug-induced liver injury in the Western world. Most studies of APAP hepatotoxicity have focused on the hepatocellular injury, but current hepatocyte-related biomarkers have delayed presentation time and a lack of sensitivity. APAP overdose can induce hepatic microvascular congestion, which importantly precedes the injury of hepatocytes. However, the underlying molecular mechanisms remain unclear. It is imperative to discover and validate sensitive and specific translational biomarkers of APAP-induced liver injury.
Methods: In this study, we assessed APAP toxicity in sinusoidal endothelial cells and hepatocytes in mice treated with overdose APAP at different time points. The underlying mechanisms of APAP overdose induced sinusoidal endothelial cell injury were investigated by RT2 Profiler PCR arrays. The impact of APAP overdose on endothelial cell function was assessed by pseudovessel formation of endothelial cells in 2D Matrigel and in vivo hepatic vascular integrity using multiphoton microscopy. Finally, the effects of APAP overdose on oxygen levels in the liver and hepatic microcirculation were evaluated by contrast enhanced ultrasonography. Potential imaging-based vascular-related markers for early detection of APAP induced liver injury were assessed.
Results: Our study confirmed that hepatic endothelial cells are an early and direct target for APAP hepatotoxicity. ICAM1-related cellular adhesion pathways played a prominent role in APAP-induced endothelial cell injury, which was further validated in primary human sinusoidal endothelial cells and human livers after APAP overdose. APAP overdose impacted pseudovessel formation of endothelial cells and in vivo hepatic vascular integrity. Use of ultrasound to detect APAP-induced liver injury demonstrated that mean transit time, an imaging-based vascular-related biomarker, was more sensitive and precise for early detection of APAP hepatotoxicity and monitoring the treatment response in comparison with a conventional blood-based biomarker.
Conclusion: Imaging-based vascular-related biomarkers can identify early and mild liver injury induced by APAP overdose. With further development, such biomarkers may improve the assessment of liver injury and the efficacy of clinical decision-making, which can be extended to other microvascular dysfunction of deep organs.
Keywords: Acetaminophen, Liver Injury, Endothelial cells, Vascular-related biomarkers, Ultrasonography
Acetaminophen (APAP) is one of the most commonly-used drugs for analgesia and fever reduction . However, APAP hepatotoxicity remains one of the leading causes of liver injury, with about 26,000 hospitalizations and 500 deaths each year in the US alone [2, 3]. Current APAP toxicity assessment and the decision to start treatment are predominantly based on the ingested dose of APAP and a timed blood APAP concentration, and measurement of liver enzyme levels such as alanine aminotransaminase (ALT) [2, 3]. Clinical uncertainty regarding the presence of hepatocyte injury prevents treatment being individualized, resulting in patients either being over treated with a time-consuming potentially harmful antidote or being undertreated with the attendant risk of serious liver injury . Hence, novel biomarkers or panels of biomarkers are required for evaluating the severity and response to treatment of APAP induced liver injury at the earliest possible time.
Most investigations on biomarker development for APAP-induced liver injury have focussed on biomarkers based on mechanistic processes involved in hepatocyte damage [4, 5]. Recent studies have examined the potential use of new blood-based biomarkers such as microRNA-122, mitochondrial DNA, and keratin-18 from hepatocytes in human and murine models [2, 6]. We have also reported the use of the hepatic optical oxidative stress index for monitoring early treatment response . However, APAP can induce hepatic microvascular injury and congestion, which importantly precedes hepatocyte injury. APAP-induced hepatic congestion in humans was first reported in 1969 . Microvascular congestion has been found to precede hepatocyte injury [9-11]. APAP is also toxic to sinusoidal endothelial cells , resulting in collapse of the sinusoidal wall and the infiltration of blood elements into the space of Disse [10, 11]. Despite the evidence that liver sinusoidal endothelial cells are an early target of APAP injury to the liver , the mechanisms of APAP-induced microvascular injury remain to be elucidated. Furthermore, there are no vascular-related biomarkers for liver injury, mainly due to the absence of suitable technologies for detecting subtle alterations in hepatic microcirculation.
Only recently, advances in ultrasonography have led to the provision of non-invasive and quantitative imaging data on the hepatic microcirculation for diagnosis and evaluation of response to treatment in patients or animals with liver cirrhosis , fatty liver disease , and liver ischemia reperfusion injury . Two ultrasound-based modalities that have potential for reporting disease progression and evaluating treatment response are photoacoustic (PA) imaging [16, 17] and contrast-enhanced ultrasonography (CEU) [18-20]. PA imaging is a hybrid modality whereby short near infrared NIR laser pulses are used to generate thermoacoustic waves in tissues containing light-absorbing endogenous or exogenous contrast agents [21, 22]. PA imaging allows parallel collection of anatomical and functional information of tissue oxygen saturation which often changes in diseased conditions [23, 24]. CEU is another imaging technique with gas-filled microbubbles as contrast agents that infuse into the intravascular space of the body, which enables assessment of tissue perfusion and is already clinically approved [25-27]. Two perfusion models are commonly used to represent either bolus kinetics or replenishment kinetics following microbubble destruction. Microbubble velocity corresponds to the blood flow velocity . mTT is the average amount of time that it takes blood to transit through a given volume of an organ. In a destruction replenishment perfusion model, the amplitude of the plateau plus offset amplitude represents the relative blood volume in a quantified region. The product of microbubbles velocity, mTT, and tissue blood volume reflects regional perfusion . All these imaging-based vascular-related markers have the potential for early detection of APAP-induced liver injury.
Therefore, in present study, we investigated the mechanisms of APAP-induced hepatic microvascular injury in mice. We demonstrate that imaging-based vascular-related markers detected by ultrasonography are more sensitive and precise than conventional blood-based biomarkers for the early detection of APAP hepatotoxicity and for monitoring treatment response.
Animals. Male BALB/c mice (8 weeks old) were purchased from the Animal Resource Centre (Perth, Western Australia). All animal procedures were approved by the Animal Ethics Committee of the University of Queensland (SOM/TRI/211/15) and were carried out according to Australian Code for the Care and Use of Animals for Scientific Purposes 8th edition. APAP (Sigma-Aldrich) was dissolved in warm saline at the concentration of 15 mg/ml. The animals received a single gavage of 500 mg/kg APAP for 2 h or 12 h. Control animals were gavaged with saline. For n-acetylcysteine (NAC) treatment group (labelled with APAP 12h + NAC in all figures), animals were treated with 300 mg/kg NAC (Sigma-Aldrich, St Louis, MO, USA) by intraperitoneal injection 1 h after APAP administration and animals were assessed 12 h after APAP administration (n = 5 for each group), according to the reported procedure .
Histology and TUNEL assay. Liver tissues were collected, fixed, and embedded in paraffin. Paraffin sections were stained with hematoxylin and eosin for histological examination. TUNEL in situ assay was used for visualization of DNA stand breaks and apoptosis.
TEM imaging of liver tissues. Liver tissues were cut into approximately 1 mm3 cubes and fixed with 2.5% glutaraldehyde. Samples were embedded with epoxy resin, sectioned and imaged using a Philips CM10 electron microscope.
Serum biochemical measurements. A blood sample was collected and the plasma concentration of ALT was measured using a Hitachi 747 analyser (Hitachi Ltd., Tokyo, Japan) at Pathology Queensland, Princess Alexandra Hospital, Brisbane, Australia.
PCR array. Human umbilical vein endothelial cells, HUVECs, and human hepatic endothelial cells SK-HEP-1, were cultured in Endothelial Basal Medium (EBM-2) supplemented with EGM-2 SingleQuot supplements (Lonza, Basel, Switzerland) or DMEM containing 10% Fetal Bovine serum, respectively at 37 °C in 5% CO2.
Cells (1-2 X105) were seeded in triplicate into 12 well dishes and when 80% confluent, were treated with APAP (Sigma Chemical Company, 20mM) for 6 h. Cells were harvested into RLT buffer (Qiagen, Hilden, Germany). RNA was extracted using a RNeasy Micro Plus Kit (Qiagen, Hilden, Germany). Total RNA was quantified using a NanoDrop Spectrophotometer (Thermo Scientific). Reverse transcription was performed with a RT2 First Strand Kit (Qiagen, Hilden, Germany) and 250 ng total RNA. qPCR was carried out using RT² SYBR Green ROX qPCR Mastermix (Qiagen, Hilden, Germany) and a RT2 Profiler™ Human Endothelial Cell Biology Array (PAHS-015Z: Qiagen, Hilden, Germany) that contains 84 genes related to endothelial cell biology using the ABI Viia7 Real-Time PCR system (Thermofisher, Waltham, MA, USA).
Data was analysed using the RT2 Profiler PCR Array Data Analysis Webportal at GeneGlobe (http://www.qiagen.com/geneglobe). CT values were normalized using the ΔCt method based on an automatic selection from the house keeping gene panel of reference genes.
Genes that exhibited more than 1.5 fold change in expression from the untreated cells, with a p-value of ≤ 0.05, were further analyzed using Ingenuity Pathway Analysis (IPA) software (Qiagen, Hildan, Germany) to determine pathway enrichment and cellular context of the differentially expressed genes.
Patients and database. Samples of publicly available human datasets of APAP overdose from the Gene Expression Omnibus, GSE74000 were analysed using GEO2R at https://www.ncbi.nlm.nih.gov/geo/ [28, 29]. GSE74000 contains gene expression microarray data of liver biopsies from healthy humans (GSM1907918 and GSM1907919) and patients APAP-induced acute liver failure (Samples GSM1907915, GSM1907916 and GSM1907917). These data address differential gene expression in severe APAP-induced liver injury and were normalised and deposited by the original authors. Fold change (2^log-FC) and adjusted p-values were calculated relative to control samples.
Validation of relative gene expression. Mouse liver tissues were disrupted and lysed in RLT buffer containing 1% β-mercaptoethanol using a TissueLyser LT (Qiagen, Hilden, Germany). Lysates were further homogenised using a QIAshredder spin column (Qiagen, Hilden, Germany) followed by centrifugation to remove insoluble material. Primary human hepatic sinusoidal endothelial cells (Zenbio) were cultured in endothelial cell growth medium-1 and treated by 20 mM APAP for 6 h. Total RNA was extracted using an RNeasy Plus Mini Kit (Qiagen, Hilden, Germany). cDNA was prepared from purified RNA using a SensiFAST cDNA Synthesis Kit (Bioline, Memphis, TN). qRT-PCR was performed using a ViiA 7 Real-Time PCR System. Gene expression relative to HPRT (mouse liver samples, SK-HEP-1 cells and HSE cells) or B2M (beta-2-microglobulin) (HUVEC cells) was normalised to untreated samples.
Flow cytometric analysis for the specific cell marker. Liver tissues were dissociated into single cells in phosphate-buffered saline/bovine serum albumin/ethylenediaminetetraacetic acid. Dissociated single cells were further incubated with antibody combinations for multi-parameter flow acquisition and analysis. A Gallios flow cytometer was used for sample acquisition, whereas unbiased data analyses were performed with Kaluza analysis software (Beckman-Coulter, Brea, CA, USA). For each sample, 300,000 single cells were analyzed. The following combinations of antibodies were used to assess the endothelial population: Rat anti-mouse VE-Cadherin FITC, VEGFR2 PE, CD31 PE-Cy7, CD34 Alexa647, and Lineage PacBlue (Becton Dickinson, Franklin Lakes, NJ, USA).
Pseudovessels were formed by SK-HEP-1 cells mixed with green fluorescent protein (GFP) human endothelial colony forming cells (ECFC) in 2D Matrigel or HUVEC mixed with ECFC. SK-HEP-1 and HUVEC cells were stained with CellTracker™ Red CMTPX Dye (Thermofisher, C34552) at 37 °C for 1 h. 2.5x104 SK-HEP-1 cells and 2.5x104 ECFC cells were then mixed and seeded into 8-well chamber slide (Cellvis, C8-1-N) coated with 100 μl of Matrigel (CORNING, 356237) diluted in 100 μl of EGM-2 after 1hr of polymerization at 37 °C. Tube formation was monitored by Olympus IX81 Live Imaging Microscope. Cell were maintained at 37 °C, 5% CO2 during the whole imaging period. Time-lapse acquisitions were performed every 30 min up to 12 h by FITC and Cy3 channels to capture the signals of ECFC cells and SK-HEP-1 and HUVEC cells, respectively.
Multiphoton microscopy imaging of liver tissues. The settings of the multiphoton microscope were as previously described and established . A DermaInspect system (JenLab GmbH, Jena, Germany) equipped with an ultrashort pulsed mode-locked tunable Titanium: Sapphire laser (Mai Tai, USA) was used for imaging of mouse livers. Rhodamine B isothiocyanate/dextran 70 kDa (DXR70) was intravenously injected to mice and DXR70 fluorescence was detected at the excitation wavelength of 900 nm, which can be differentiated from autofluorescence of the liver. A bandpass filter was employed to detect fluorescence emission in the spectral range of 350 to 650 nm. High (40×, Zeiss) magnification objectives were used with the laser power set at 15 mW. DXR70 (250 μg/ml, 100 μl) was injected intravenously for labelling hepatic sinusoids.
In vivo ultrasound-based imaging. All PA and CEU imaging were performed on a Vevo 2100 LAZR imaging station (Visualsonics, Toronto, Canada). Depilated mice were anaesthetized by 1.5% isoflurane inhalation (1 L oxygen / min) delivered via a nose cone and kept warm on a heated stage. B-mode images were used to locate the same position of the left liver lobe on every mouse for functional analysis.
PA imaging and analysis. Photoacoustic data were acquired using a LZ400 linear array transducer, which has a centre frequency of 30 MHz. The laser energy was calibrated and optimised using an in-built energy meter prior to measurements. All mice were scanned with their surface 10 mm below the transducer, and the PA signals were amplified 40 dB. Tissue oxygen saturation can be quantified based on differences in the absorbance spectrum between oxygenated haemoglobin (Hboxy) and deoxygenated haemoglobin (Hbdeoxy). For oxygen saturation measurements, dual-wavelength photoacoustic imaging at 750 nm and 850 nm was performed. The relative amounts of Hboxy and Hbdeoxy were calculated using the OxyZated Tool within Vevolab software (Visualsonics, Toronto, Canada). Oxygen saturation was defined as [100% × Hboxy] / [Hboxy + Hbdeoxy]. The threshold for Hb was set at 20% of maximal intensity, a reasonable value determined empirically .
CEU imaging and analysis. Non-linear CEU imaging was performed using a MS250S transducer at a frequency of 18 MHz, power of 4%, and contrast gain of 30 dB, following intravenous bolus infusion of 50 ul of non-targeted microbubbles (Visualsonics, Toronto, Canada). Hepatic microvasculature perfusion kinetics were assessed by two models conducted in sequence; (1) bolus perfusion model, then (2) destruction-replenishment perfusion model . Bolus perfusion imaging is based on the wash-in and wash-out kinetics of microbubbles after bolus injection . For this, cine loop recordings at 11 frames/ sec in contrast mode commenced immediately after tail vein injection of a single bolus of microbubbles at a consistent speed. Calculations of WiR were made using the VevoCQ Tool within the Vevolab software (Visualsonics, Toronto, Canada), which is defined as
Where with the error function erf defined as . Where f(t) represents best-fit function of echo-power; t represents time; A is amplitude of the plateau and O represents offset amplitude. Immediately after bolus perfusion modeling, the destruction-replenishment model was employed. This involves first infusion of microbubbles to achieve a constant flow rate within an ROI, followed by a sudden application of a high powered ultrasound pulse (100% power) to burst the microbubbles, enabling calculation of other perfusion parameters such as mean transit time (mTT) and relative blood volume (rBV). Cine loop recordings in contrast mode were at 11 frame/sec. To determine the perfusion parameters, audio video interleave files for the ROIs were exported then calculated by a nonlinear analysis in GraphPad Prism 7.
Statistical analysis. Mixed-model regression analysis was performed using the R statistics package to fit the perfusion changes. Statistical significance between two groups was determined using a two-tailed Student t test. One-way ANOVA with unequal variances followed by Turkey's multiple comparison test was used to assess significance between multiple groups. All statistical analyses were performed using GraphPad Prism v 7 (GraphPad Software Inc., La Jolla, California). Results were considered statistically significant with a p-value < 0.05.
Endothelial cells are an early target of APAP-induced liver injury. To assess APAP toxicity in liver endothelial cells and hepatocytes, male BALB/c mice were treated with an overdose of APAP at different time points or with NAC treatment after APAP overdose, which is the first choice for treatment of APAP poisoning in the clinic [2, 32]. As shown in Figure 1A-B, mild to moderate sinusoidal congestions was observed 2 h after APAP overdose while there was no obvious hepatocyte injury. Similarly, hepatic DNA fragmentation was not detected using TUNEL staining (Figure 1C) at this time point, indicating there was no obvious hepatocyte apoptosis. Furthermore, ALT levels were not significantly increased compared to the control group (Figure 1E). However, sinusoidal endothelial cell swelling was detected by TEM 2 h after APAP administration (Figure 1D), and the endothelial cell injury became more severe at 12 h. There was a significant increase in the number of TUNEL-positive hepatocytes and ALT levels 12 h after APAP overdose. These results suggest that APAP induced sinusoidal endothelial cell injury precedes hepatocyte injury and sinusoidal endothelial cells are an early and direct target for APAP hepatotoxicity.
Mechanisms of APAP-induced sinusoidal endothelial cell injury. A dose-dependent toxicity of APAP was shown in HUVEC and SK-HEP-1 cells (Figure S1) and cell viability decreased as early as 6 h after APAP treatment. Both cell lines treated with 20 mM APAP for 6 and 24 h, reduced the viability of SK-HEP-1 cells to 79% and 35%, and the viability of HUVECs to 57% and 35%, respectively. To define cell signalling pathways involved in the APAP-induced endothelial cell injury, a total of 84 key genes of endothelial cell biology pathways (including vasoconstriction and vasodilation, inflammatory response, cell adhesion, apoptosis and coagulation) were simultaneously assayed using RT2 Profiler PCR arrays. SK-HEP-1 and HUVECs were treated with 20 mM APAP for 6 h. Scatter plots of relative mRNA expression for both cell lines before and after APAP treatment are shown in Figure S2. Seven genes were commonly upregulated and 5 genes downregulated in both cells lines after APAP overdose (Figure 2A and Figure S3).
APAP induced damage in the mouse liver. (A). H&E staining of liver sections. Necrosis is observed as pale eosinophilic staining and loss of hepatocyte nuclei. (B). High magnification view of area seen in box of (A) highlighted the sinusoidal area in liver (arrows indicated sinusoidal spaces). Sinusoidal congestion can be seen in livers after 2 h and 12 h APAP overdose as well as the NAC treated group. (C). Representative images of TUNEL stained liver sections in each group. Example of TUNEL positive cells are indicated by black arrows. (D). Transmission electron microscopic images of APAP-related acute sinusoidal endothelial cell injury. Sinusoidal endothelial cells are indicated by white arrows.
The expression of key genes was validated in both cell lines and primary human hepatic sinusoidal endothelial cells (Figure 2B and Figures S4 and S5A) and mouse liver tissues by qRT-PCR (Figure 2D). ICAM1 and SELPLG within the cellular adhesion molecular pathway, were both significantly up-regulated in two cell lines (Figure 2A-B, p ˂ 0.05). These results indicate that endothelial cells become activated after APAP treatment. APAP overdose also led to increased expressions of CCL5, a chemokine central to leukocyte recruitment, VEGFA, a key regulator of vascular permeability and EDN1 which is involved in regulation of vascular tone [33, 34]. Ingenuity Pathway analysis (IPA) was performed to ascertain the cellular context of the differentially expressed genes (Figure 2C) and indicated APAP induced leukocyte accumulation and endothelial cell damage. These results suggest ICAM plays a prominent role in the endothelial cells response to APAP. Importantly, analysis of gene expression of human livers using a public microarray database showed that there was a 7 fold increase in EDN1 expression and an overall 1.6 fold increase in ICAM expression after APAP overdose (Figure S5B), which were consistent with our in vitro and in vivo results.
Molecular mechanisms of APAP-induced sinusoidal endothelial cell injury. (A). The commonly up-regulated and down-regulated genes in SK-HEP-1 and HUVEC cells after APAP overdose for 6 h (fold ˃ 1.5, P ˂ 0.05). Data are presented as the mean ± SD (n = 3). (B). qPCR of VEGFA, EDN1, SELPLG, ICAM1, CCL5, and PTGS2 mRNAs in SK-HEP-1 cells after APAP overdose for 6 h. (C). IPA of differentially regulated gene networks after APAP overdose. ICAM1 is involved in the regulation of several genes related to apoptosis, necrosis, and activation of endothelial cells. (D). qRT-PCR of Icam1, Vegfa, Ccl5, and Edn1 mRNA in mouse liver tissues after APAP overdose for 2 and 12 h. (E). The expression of ICAM1 was demonstrated in mouse livers by immunofluorescence staining of ICAM1 (green color) with DAPI (blue color) labelling of the nuclei (CV: central vein; scale bar: 50 μm). (F). Mean number of VECAD+CD34+CD31+Lin- cells isolated from liver tissues were analysed and quantified by flow cytometry. Data are presented as the mean ± SD from 5 independent experiments. **P ˂ 0.01.
To visualise overall ICAM1 expression in APAP overdose animals, liver sections were stained with ICAM1 antibody (Figure 2E). Enhanced ICAM1 expression was found in the centrilobular regions of mouse livers after APAP overdose, further demonstrating the induction of endothelial adhesion molecules. We next examined the impact of APAP overdose on the endothelial cell-specific adhesion protein VE-cadherin (VECAD), which plays an important role in endothelial cell biology through control of the cohesion and organisation of the intercellular junctions . The number of VECAD positive endothelial cells decreased significantly 2 h after APAP overdose, and was further depleted at 12 h (Figure 2F, **P ˂ 0.01). Taken together, the results indicate that APAP overdose resulted in the activation of endothelial adhesion pathways and vascular damage evidenced by opening of the endothelial junctions, inflammatory cell homing, increased microvasculature permeability and vascular leakage.
APAP overdose effects on vascular tube formation and vascular integrity. As the level of ICAM1 is related to vascular tube formation , we further evaluated the effects of APAP on the formation of pseudo-vessels by SK-HEP-1 and HUVECs on 2D Matrigel. Pseudo-vessel networks were sufficiently formed by mixed SK-HEP-1 cells and human endothelial colony forming cells (ECFC) as well as HUVEC and ECFC (Figure 3A and Figure S6A, Video S1&2). Cells treated with 40 mM APAP did not form pseudo vessels but cells treated with 20mM APAP commenced vessel formation after 12 h (Figure 3A and Figure S6A, Video S3-6). To gain further insight into the vessel morphometry after APAP treatment, we quantified the number of pseudovessels, nodes and complete circular capillaries as well as the mean length of pseudovessels. Large numbers of vessels were lost after 20mM APAP treatment and the node number was significantly less than control groups (Figure 3B and Figure S6B). All these findings indicated APAP overdose impaired vascular tube formation of endothelial cells.
VECAD positive endothelial cells decreased in mouse livers after APAP treatment for 2h (Figure 2F), indicating opening of the endothelial junctions and vascular leakage. Then we examined the sinusoidal integrity by intravenous injection of DXR70 and imaged the liver by multiphoton microscopy according to the reported method . DXR70 was evenly distributed in sinusoids but was not taken up by hepatocytes in the control liver, while tapered sinusoids were shown in the APAP treated liver and DXR70 fluorescence was found in a few hepatocytes (Figure 3C). These results not only reflected compromised plasma membranes and sinusoidal vascular leak at the early phase of APAP overdose, but also indicated the injury of hepatocytes which do not show pathological changes.
APAP overdose effects on oxygen levels in the liver. Based on our findings above, we further explored potential imaging-based vascular-related biomarkers for early detection of APAP-induced liver injury using ultrasound. First, PA imaging was performed on the left lobe of the mouse liver, and relative oxygen saturation levels were plotted over time. Heat-map images of oxygen saturation of the liver, co-registered with greyscale B mode ultrasound, of representative mice from each group are shown in Figure 4A. Levels of hepatic oxygen saturation were associated with APAP-induced liver injury (Figure 4A). Oxygen saturation of the liver parenchyma was significantly decreased to a mean of 64 % after APAP overdose for 2 h and 56 % after APAP overdose for 12 h compared to control liver of over 75 % (Figure 4B). NAC treatment improved oxygen saturation to approximately 63 %.
We then assessed how oxygen saturation may reflect degrees of liver injury and found that oxygen saturation levels of livers correlated well with the ALT level from early stage of APAP overdose in the blood (R = 0.67, p = 0.0002) (Figure 4C) but not the late stage (this correlation only include control, 2 h APAP overdose, and NAC treated groups). Significant changes in hepatic oxygen saturation were detected as early as 2 h after APAP overdose, which preceded ALT changes. Thus, we suggest that oxygen saturation can be used as a biomarker for early or mild APAP-induced liver injury especially when ALT levels are normal.
APAP overdose effects on hepatic microcirculation. Immediately after photoacoustic imaging, contrast-enhanced ultrasound was performed to assess hepatic perfusion. Obvious effects of APAP overdose on perfusion parameters were observed using both bolus and destruction-replenishment perfusion models (Figure 5 and 6). High resolution parametric perfusion maps were developed to visualise and monitor microvascular perfusion changes in space and time, providing insights into the distributional changes in the microcirculation of the injured liver induced by APAP overdose. Real-time perfusion imaging clearly showed the decreased wash-in rate (WiR) of the initial contrast microbubbles in the APAP overdose treated liver (Figure 5A and Video S7-10), as presented in pseudo colour (the colour changed to blue and dark blue). The parametric maps revealed the WiR in the liver decreased as early as 2 h post APAP overdose and was more obviously reduced after 12 h. In NAC treated mice, delayed perfusion recovery persisted in the liver. These parametric perfusion images reflected rapid and marked reduction in the microvascular perfusion of the liver after APAP overdose. Representative single liver perfusion tracings were obtained for the control, 2 h and 12 h APAP overdose, and NAC treated animals, by monitoring initial filled-induced nonlinear signals of intact microbubbles (Figure 5B). WiR was significantly longer in the liver 12 h after APAP overdose compared to the control group (Figure 5C), indicating that APAP overdose decreased the blood flow velocity in the liver. In contrast, NAC treated livers showed obvious improvement in WiR compared to 12 h APAP overdose.
To further evaluate the hepatic microcirculation, the destruction-replenishment perfusion model was performed (Figure 6 and Video S11-14). Representative nonlinear contrast images of the liver after acutely collapsing microbubbles with an ultrasound burst and the refilled microbubbles were monitored (Figure 6A). By 0.5 s, the liver started to refill with microbubbles. After 1-2 s, the plateau phase of perfusion was reached in the control liver, whereas the hepatic perfusion 2 h after APAP overdose and in NAC treatment groups returned to levels close to baseline by about 5 s. The refilled microbubbles of 12 h after APAP overdose reached the flat phase over 10 s. The representative single hepatic perfusion tracings were obtained by monitoring refill induced nonlinear signals of intact microbubbles after immediately collapsing by ultrasound burst (Figure 6B). The imaging signals from slight breathing motion artefacts were digitally subtracted from the initial frames. The perfusion parameter, mTT, represents the microbubbles transit time in the regional liver and was determined by the destruction-replenishment perfusion model based on a curve fitting of data in Figure 6B. Compared to control liver, APAP overdose significantly decreased hepatic perfusion at 2 h, where mTT was significantly extended at 2 h and further increased at 12 h after APAP overdose (Figure 6C). After NAC treatment, hepatic perfusion was improved and mTT was significantly shorter compared to untreated 12 h APAP group. On the other hand, relative blood volume (rBV) quantified using this model (Figure 6D) showed that although rBV of the liver after APAP overdose for 12 h is less than normal liver, it is not statistically significant. This parameter is therefore not as sensitive as mTT and WiR. Furthermore, WiR correlated slightly with ALT levels (Figure 5D, R = 0.53, P = 0.0003), while mTT showed a significant correlation with ALT levels (Figure 6E, R = 0.9527, P ˂ 0.0001), suggesting these imaging biomarkers accurately detect the liver injury induced by APAP overdose at the early phase and reflect the treatment response of injured livers. Importantly, mTT is the optimal imaging marker for early detection of APAP induced liver injury.
APAP overdose affected vascular tube formation and vascular integrity in the liver. (A). 2D vascular tube formation assay using ECFC cells (green) and SK-HEP-1 cells (Red) either untreated and treated with various concentrations of APAP. Images represent t = 12 h, scale bar: 200 μm. (B). Quantitative analysis of vessel morphometry, including number of nodes, number of vessels, mean length of vessels, and number of complete circular capillaries. *P ˂ 0.05, **P ˂ 0.01; one-way ANOVA. (C) Multiphoton microscopy imaging of mouse liver after intravenous injection of DXR70. Red color represents DXR fluorescence in the liver and circled cells (dashed lines) indicate DXR uptake by hepatocytes. Scale: 40 μm.
APAP overdose and NAC treatment effect on oxygen saturation of the liver. (A). Representative PA images (right panel) co-registered with greyscale B mode imaging (left panel) of livers in each group. The heat map represents oxygen saturation levels ranging from 100 % (red) to 0% (dark blue). sO2 Avr total represents the oxygen saturation calculated by percentage of oxygenated haemoglobin in relation to total haemoglobin. (B). The average oxygen saturation of each group detected by PA imaging (n = 5). (C). Relationship between the levels of oxygen saturation of livers and ALT levels in control, APAP 2 h, and APAP+NAC groups (R2 = 0.67, P = 0.0002, n = 15).
It is imperative to discover and validate more sensitive and specific translational biomarkers of APAP-induced liver injury. Most studies of APAP hepatotoxicity have focused on the hepatocellular injury , and current hepatocyte-related biomarkers have delayed presentation time and a lack of sensitivity, remaining a critical impediment to the treatment of APAP overdose . Although studies have shown APAP is also toxic to endothelial cells and precedes the injury of hepatocytes [9-11, 39], the mechanisms involved are still unclear. This study revealed the detailed mechanisms of APAP-induced endothelial cell injury. Our data also provide evidence that imaging-based vascular-related markers are able to detect early and mild APAP-induced liver injury.
By profiling gene expression in endothelial cells after APAP overdose, we have shown that ICAM1 is at the apex of signalling pathways activated in endothelial cells during APAP hepatotoxicity. ICAM1 is an intercellular adhesion molecule continuously present in membranes of endothelial cells . ICAM1 and SELPLG, which were shown to be upregulated in both cells response to APAP, are involved in the activation of endothelial cells and binding of leukocytes to endothelial cells during leukocyte docking and transmigration through the endothelial cell layer . A few studies have shown that ICAM-1 expression was increased in mouse livers after APAP administration for 10 h and the marked protection against APAP hepatotoxicity was found in ICAM-1-deficient mice. But these studies demonstrated APAP overdose could induce the increase of ICAM-1 in tissues at the late or severe phase of liver injury, which was associated with the accumulation of neutrophils [41, 42]. This is the first study to show the upregulated ICAM-1 expression of human endothelial cells in direct response to APAP overdose at the early phase. Subsequent experiments showed that higher expression of ICAM1 was found in endothelial cells around the hepatic central vein as early as 2 h after APAP overdose, which is consistent with APAP induced hepatotoxicity initially occurring in the areas adjacent to the central vein . Along with these changes there was up-regulation of chemokines involved in leukocyte recruitment (CCL5) and growth factors involved with vascular tone (EDN1) . Analysis of human genes of primary human hepatic sinusoidal endothelial cells and patient tissues after APAP overdose induced acute liver failure showed that ICAM1 and EDN1 obviously increased, further supporting our findings. This is the only available public dataset of liver gene expression after APAP overdose in comparison to healthy liver. All these results suggest induction of vascular damage and an inflammatory response, leading to opening of the endothelial junctions and inflammatory cell recruitment.
APAP overdose and NAC treatment effect on wash-in rate (WiR) of the liver. (A). Left panel of images (B mode) represent mouse liver obtained using normal ultrasound to indicate liver anatomy (L: liver). High resolution parametric perfusion maps of the regional liver (Right panel). Livers of control, APAP overdose for 2 h and 12 h, and NAC treatment are displayed in pseudocolor. Red and orange colors represent the higher WiR, while blue color represents the lower WiR. After APAP overdose for 2 h, the red color changed to yellow and green in the liver. A noticeable blue band is present after APAP overdose for 12 h while NAC treated liver was less blue and more green and yellow. (B). Representative single liver perfusion tracings were obtained for the control, APAP overdose for 2 h and 12 h, and NAC treated animals, by monitoring initial filled nonlinear signals of intact microbubbles. (C). The average WiR of each group based on nonlinear contrast imaging (n = 5). (D). Relationship between the WiR and ALT levels in each group (R2 = 0.53, P = 0.0003, n = 20).
VECAD is an endothelial-specific cadherin that plays an important role in cohesion and organisation of intercellular junctions [35, 44]. The depletion of VECAD positive endothelial cells further suggests APAP overdose could result in the opening of the endothelial junctions, increased vascular permeability and vascular leakage. We also confirmed that APAP overdose impaired the functions of endothelial cells, such as vascular tube formation and vascular integrity, at the early stage of liver injury. Therefore, endothelial cells are an early target of APAP hepatotoxicity and vascular related biomarkers are potentially useful for early diagnosis of APAP induced liver injury. Hepatic perfusion may mirror the main consequences of APAP overdose induced liver injury at the early stage, which can be characterised by impairment of hepatic microcirculation and liver tissue oxygenation. Although standard imaging techniques allow the structural investigation of the vascular consequences of liver injury, the hepatic microcirculation is not easily assessed and monitored. Noninvasive measurement of hepatic blood flow by Doppler is not sensitive or accurate  and evaluation relies on indirect methods such as indocyanine green  and galactose clearance test. Taken together, new techniques capable of giving reliable information on liver microcirculation are still needed.
APAP overdose and NAC treatment effect on mean transit time (mTT) of the liver perfusion. (A). Nonlinear contrast imaging sequence at various time points of refilled microbubbles after “burst” at the regional area of the liver in each group. Top panel of images (B mode) represent the mouse liver obtained using normal ultrasound to indicate liver anatomy (L: liver). By 0.5 s, livers of control and APAP overdose for 2 h started to refill with microbubbles. After 1 s, control liver was almost fully refilled by microbubbles. By 5 s, livers of APAP overdose for 2 h and NAC treatment were fully refilled by microbubbles, while the plateau of liver perfusion was not reached in the liver of APAP overdose for 12 h. (B). Representative liver perfusion curves based on the sequence of refilled microbubbles after a high mechanical index “burst”. Bursts are indicated by arrows. The dots are linearized signals recorded by contrast imaging and the curve was fitted by the Destruction-replenishment perfusion model. (C). The average mTT of each group based on nonlinear contrast imaging and calculated by the destruction-replenishment model. Data are presented as the mean ± SD (n = 5, **P ˂ 0.01). (D). The average rBV of each group. Data are presented as the mean ± SD (n = 5). (E). Relationship between the mTT and ALT levels in each group (R2 = 0.9527, P ˂ 0.0001, n = 20).
Therefore, using advanced ultrasound technologies including PA and CEU imaging, we have demonstrated that APAP overdose resulted in a reproducible pattern of injury to the hepatic microcirculation. By PA imaging, oxygenation saturation of the liver could be evaluated based on measurements of oxygenated and deoxygenated haemoglobin of the tissue in vivo . Oxygenation saturation of the liver was found to be significantly reduced from 2 h after APAP overdose. The changes in oxygen content of blood vessels indicate the changes in blood flow and circulation . In recent years, CEU has been proposed as a tool for measuring real-time functional aspects of microvasculature in biological tissues [46-49]. Microbubbles can be used as contrast agents to generate nonlinear signals during CEU imaging, allowing real-time recording and analysis in all cases. We applied microbubbles as contrast agents to systematically investigate the effect of APAP overdose on hepatic microcirculation. Our results show that impaired hepatic microcirculation can be detected as early as 2 h after APAP overdose using CEU. We showed that not only was there a delay in reaching peak blood flow in the liver after APAP overdose, but overall reduced blood volume and blood flow (reflected by WiR) were also found in the liver. NAC treatment, which is shown to reduce the hepatotoxicity of APAP, significantly improves the conditions of hepatic microcirculation. These results suggest alterations in liver microperfusion occur in APAP-induced liver injury and potentially that these microcirculatory changes are important in injury progression and recovery. Importantly, the damage to hepatic microcirculation has been shown to precede evidence of hepatocyte necrosis and apoptosis. A further benefit of CEU is that there are no adverse effects related to the procedure . Thus, we believe that CEU is a feasible novel, quantitative, and noninvasive technique for evaluating hepatic microcirculation for assessing APAP toxicity in clinic.
The vascular related imaging biomarkers correlated well with a standard plasma biomarker, ALT . A parameter of oxygenation saturation correlated well with early and mild liver injury but may not reflect late or severe liver injury. This may be due to the oxygenation measurements, where haemoglobin is not present at a significant level once the tissue is necrotic or poorly vascularised . Hence, this parameter is more suitable for detecting early or mild liver injuries, especially those with normal ALT levels. WiR, which is determined by a bolus infusion model, showed a moderate correlation with ALT levels. WiR was affected by the injection rate, which may also limit the use of this parameter as a precision imaging biomarker. Our study showed that mTT had the best correlation with ALT levels compared to other parameters. It is noted that mTT was significantly prolonged 2 h after APAP overdose, while ALT levels did not significantly increase, indicating mTT can accurately evaluate the liver injury both at early phase and late phase. Overall, mTT was more reliable and gave a more comprehensive depiction of APAP induced liver injury at both early and mild as well as late phases of liver injury. We suggest this vascular related imaging biomarker also can be applied to monitor and evaluate microvascular dysfunction such as ischemia reperfusion [51, 52].
In conclusion, this study investigated the mechanisms of APAP overdose induced microcirculatory injury and endothelial cell damage for the first time. Systematic evaluation of the hepatic microcirculation after APAP overdose by ultrasonography showed that hepatic microvascular disturbance correlates well with the severity of APAP hepatotoxicity and precedes the necrosis of hepatocytes. Of all evaluated parameters, mTT by ultrasonography has better reproducibility and reliability, and could serve as a noninvasive biomarker for the diagnosis and treatment evaluation of APAP overdose. With further development, this novel technique can be translated to clinical use to improve the assessment of liver injury and the speed of clinical decision-making.
Supplementary figures and tables.
The study was supported by grants from National Health and Medical Research Council of Australia (APP1125794 and APP1055176), Gallipoli Medical Research Foundation, and Three-year Talent Program of Changzheng Hospital (2018). This research was carried out at the Translational Research Institute, Woolloongabba, QLD 4102, Australia. The authors also acknowledge the TRI for providing the excellent research environment and core facilities that enabled this research, particularly Preclinical Imaging, Biological Resources, and Microscopy. TRI is supported by grants from the Australian and Queensland Governments. The authors also acknowledge the Australian Research Council for funding the purchase of the Vevo LAZR through a LIEF grant (LE150100067).
H.W. designed and performed the study. L.J.B performed the PCR array and in vitro experiments. J.P. performed the flow cytometric experiments and data analysis. B.W.C. performed ultrasound imaging and helped with data analysis. K.R.B., D.H.G.C. and M.S.R. reviewed and edited the manuscript. V.C.C. performed and analysed TEM imaging. X. L. provided human primary HSE cells and performed validation experiments of PCR array. X.L. performed animal experiments. H.Y. performed functional assessments of endothelial cells. W.G. and X.L. developed the concept, designed the study and wrote the manuscript.
The authors have declared that no competing interest exists.
1. Barbier-Torres L, Iruzubieta P, Fernandez-Ramos D, Delgado TC, Taibo D, Guitierrez-de-Juan V. et al. The mitochondrial negative regulator MCJ is a therapeutic target for acetaminophen-induced liver injury. Nat Commun. 2017 8
2. Antoine DJ, Dear JW, Lewis PS, Platt V, Coyle J, Masson M. et al. Mechanistic biomarkers provide early and sensitive detection of acetaminophen-induced acute liver injury at first presentation to hospital. Hepatology. 2013;58:777-787
3. Andersson DA, Gentry C, Alenmyr L, Killander D, Lewis SE, Andersson A. et al. TRPA1 mediates spinal antinociception induced by acetaminophen and the cannabinoid Delta(9)-tetrahydrocannabiorcol. Nat Commun. 2011 2
4. Ghosh A, Sil PC. Protection of acetaminophen induced mitochondrial dysfunctions and hepatic necrosis via Akt-NF-kappa B pathway: Role of a novel plant protein. Chem-Biol Interact. 2009;177:96-106
5. Kass GEN. Paracetamol (acetaminophen)-induced hepatic injury: Apoptosis, necrosis and mitochondria. Toxicology. 2005;213:217-217
6. Wang K, Zhang SL, Marzolf B, Troisch P, Brightman A, Hu ZY. et al. Circulating microRNAs, potential biomarkers for drug-induced liver injury. P Natl Acad Sci USA. 2009;106:4402-4407
7. Wang HL, Zhang R, Bridle KR, Jayachandran A, Thomas JA, Zhang WZ. et al. Two-photon dual imaging platform for in vivo monitoring cellular oxidative stress in liver injury. Sci Rep-Uk. 2017 7
8. Rose PG. Paracetamol Overdose and Liver Damage. Brit Med J. 1969;1:381 -&
9. Lim SP, Andrews FJ, OBrien PE. Acetaminophen-induced microvascular injury in the rat liver: Protection with misoprostol. Hepatology. 1995;22:1776-1781
10. Ito Y, Bethea NW, Abril ER, McCuskey RS. Early hepatic microvascular injury in response to acetaminophen toxicity. Microcirculation. 2003;10:391-400
11. Ito Y, Machen NW, Abril ER, McCuskey RS. Effects of acetaminophen on hepatic microcirculation in mice. Comp Hepatol. 2004;3(Suppl 1):S33
12. McCuskey RS, Bethea NW, Wong J, McCuskey MK, Abril ER, Wang X. et al. Ethanol binging exacerbates sinusoidal endothelial and parenchymal injury elicited by acetaminophen. J Hepatol. 2005;42:371-377
13. Berzigotti A, Nicolau C, Bellot P, Abraldes JG, Gilabert R, Garcia-Pagan JC. et al. Evaluation of regional hepatic perfusion (RHP) by contrast-enhanced ultrasound in patients with cirrhosis. Journal of Hepatology. 2011;55:307-314
14. Kolachala VL, Jiang R, Abramowsky CR, Gupta NA. Contrast-Based Real-time assessment of microcirculatory changes in a fatty liver after ischemia reperfusion injury. J Pediatr Gastr Nutr. 2016;62:429-436
15. Li HY, Lu JN, Zhou XF, Pan DH, Guo DQ, Ling HY. et al. Quantitative analysis of hepatic microcirculation in rabbits after liver ischemia-reperfusion injury using contrast-enhanced ultrasound. Ultrasound Med Biol. 2017;43:2469-2476
16. Knox HJ, Hedhli J, Kim TW, Khalili K, Dobrucki LW, Chan J. A bioreducible N-oxide-based probe for photoacoustic imaging of hypoxia. Nat Commun. 2017 8
17. Chen YS, Yoon SJ, Frey W, Dockery M, Emelianov S. Dynamic contrast-enhanced photoacoustic imaging using photothermal stimuli-responsive composite nanomodulators. Nat Commun. 2017 8
18. St Clair JR, Ramirez D, Passman S, Benninger RKP. Contrast-enhanced ultrasound measurement of pancreatic blood flow dynamics predicts type 1 diabetes progression in preclinical models. Nat Commun. 2018 9
19. Zachs DP, Offutt SJ, Graham RS, Kim Y, Mueller J, Auger JL. et al. Noninvasive ultrasound stimulation of the spleen to treat inflammatory arthritis. Nat Commun. 2019 10
20. Mallidi S, Watanabe K, Timerman D, Schoenfeld D, Hasan T. Prediction of tumor recurrence and therapy monitoring using ultrasound-guided photoacoustic imaging. Theranostics. 2015;5:289-301
21. Wang HL, Li XX, Tse BWC, Yang HT, Thorling CA, Liu YX. et al. Indocyanine green-incorporating nanoparticles for cancer theranostics. Theranostics. 2018;8:1227-1242
22. Chaudhary Z, Khan GM, Abeer MM, Pujara N, Tse BW, McGuckin MA. et al. Efficient photoacoustic imaging using indocyanine green (ICG) loaded functionalized mesoporous silica nanoparticles. Biomater Sci. 2019;7:5002-5015
23. Gerling M, Zhao Y, Nania S, Norberg KJ, Verbeke CS, Englert B. et al. Real-time assessment of tissue hypoxia in vivo with combined photoacoustics and high-frequency ultrasound. Theranostics. 2014;4:604-613
24. Beard P. Biomedical photoacoustic imaging. Interface Focus. 2011;1:602-631
25. Albrecht T, Blomley MJK, Cosgrove DO, Taylor-Robinson SD, Jayaram V, Eckersley R. et al. Non-invasive diagnosis of hepatic cirrhosis by transit-time analysis of an ultrasound contrast agent. Lancet. 1999;353:1579-1583
26. Fischer K, Meral FC, Zhang YZ, Vangel MG, Jolesz FA, Ichimura T. et al. High-resolution renal perfusion mapping using contrast-enhanced ultrasonography in ischemia-reperfusion injury monitors changes in renal microperfusion. Kidney Int. 2016;89:1388-1398
27. Sirsi SR, Borden MA. Advances in ultrasound mediated gene therapy using microbubble contrast agents. Theranostics. 2012;2:1208-1222
28. James LP, McCullough SS, Lamps LW, Hinson JA. Effect of N-acetylcysteine on acetaminophen toxicity in mice: Relationship to reactive nitrogen and cytokine formation. Toxicological Sciences. 2003;75:458-467
29. Rodrigues RM, Govaere O, Roskams T, Vanhaecke T, Rogiers V, De Kock J. Gene expression data from acetaminophen-induced toxicity in human hepatic in vitro systems and clinical liver samples. Data Brief. 2016;7:1052-1057
30. Liang XW, Grice JE, Zhu Y, Liu D, Sanchez WY, Li Z. et al. Intravital multiphoton imaging of the intravital multiphoton imaging of the selective uptake of water-dispersible quantum dots into sinusoidal liver cells. Small. 2015;11:1711-1720
31. Wang HJ, Hristov D, Qin JL, Tian L, Willmann JK. Three-dimensional dynamic contrast-enhanced US imaging for early antiangiogenic treatment assessment in a mouse colon cancer model. radiology. 2015;277:424-434
32. Uzi D, Barda L, Scaiewicz V, Mills M, Mueller T, Gonzalez-Rodriguez A. et al. CHOP is a critical regulator of acetaminophen-induced hepatotoxicity. Journal of Hepatology. 2013;59:495-503
33. Sutton TA, Fisher CJ, Molitoris BA. Microvascular endothelial injury and dysfunction during ischemic acute renal failure. Kidney Int. 2002;62:1539-1549
34. Schon MP, Zollner TM, Boehncke WH. The molecular basis of lymphocyte recruitment to the skin: Clues for pathogenesis and selective therapies of inflammatory disorders. J Invest Dermatol. 2003;121:951-962
35. Sauteur L, Affolter M, Belting HG. Distinct and redundant functions of Esama and VE-cadherin during vascular morphogenesis. Development. 2017;144:1554-1565
36. Chabut D, Fischer AM, Hettey D, Colliec S. Low molecular weight fucoidan promotes FGF-2-induced vascular tube formation by human endothelial cells, with decreased PAI-1 release and ICAM-1 downregulation. Thromb Res. 2004;113:93-95
37. Li FC, Huang GT, Lin CJ, Wang SS, Sun TL, Lo SY. et al. Apical membrane rupture and backward bile flooding in acetaminophen-induced hepatocyte necrosis. Cell Death Dis. 2011;2:e183
38. Xie Y, McGill MR, Dorko K, Kumer SC, Schmitt TM, Forster J. et al. Mechanisms of acetaminophen-induced cell death in primary human hepatocytes. Int J Toxicol. 2015;34:80-81
39. Walker RM, Racz WJ, Mcelligott TF. Acetaminophen-induced hepatotoxic congestion in mice. Hepatology. 1985;5:233-240
40. Lawson C, Wolf S. ICAM-1 signaling in endothelial cells. Pharmacol Rep. 2009;61:22-32
41. Ishida Y, Kondo T, Ohshima T, Fujiwara H, Iwakura Y, Mukaida N. A pivotal involvement of IFN-gamma in the pathogenesis of acetaminophen-induced acute liver injury. FASEB J. 2002;16:1227-1236
42. Liu ZX, Han D, Gunawan B, Kaplowitz N. Neutrophil depletion protects against murine acetaminophen hepatotoxicity. Hepatology. 2006;43:1220-1230
43. Kennedy RC, Smith AK, Ropella GEP, McGill MR, Jaeschke H, Hunt CA. Propagation of pericentral necrosis during acetaminophen-induced liver injury: evidence for early interhepatocyte communication and information exchange. Toxicol Sci. 2019;169:151-166
44. Xiao KY, Allison DF, Kottke MD, Summers S, Sorescu GP, Faundez V. et al. Mechanisms of VE-cadherin processing and degradation in microvascular endothelial cells. J Biol Chem. 2003;278:19199-19208
45. Bolognesi M, Sacerdoti D, Merkel C, Gatta A. Relationship between portal blood-flow measured by image-Directed doppler ultrasonography and hepatic blood-flow measured by indocyanine green constant infusion in patients with cirrhosis. J Clin Ultrasound. 1995;23:297-303
46. Pessayre D, Lebrec D, Descatoire V, Peignoux M, Benhamou JP. Mechanism for Reduced Drug Clearance in Patients with Cirrhosis. Gastroenterology. 1978;74:566-571
47. Roddie IC, Shepherd JT, Whelan RF. Evidence from venous oxygen saturation measurements that the increase in forearm blood flow during body heating is confined to the Skin. J Physiol-London. 1956;134:444-450
48. Lim A, Cosgrove D. Functional studies. Eur Radiol. 2004;14:P110-P115
49. Lin F, Shelton SE, Espindola D, Rojas JD, Pinton G, Dayton PA. 3-D ultrasound localization microscopy for identifying microvascular morphology features of tumor angiogenesis at a resolution beyond the diffraction limit of conventional ultrasound. Theranostics. 2017;7:196-204
50. Jin G, Qiu G, Wu D, Hu Y, Qiao P, Fan C. et al. Allogeneic bone marrow-derived mesenchymal stem cells attenuate hepatic ischemia-reperfusion injury by suppressing oxidative stress and inhibiting apoptosis in rats. Int J Mol Med. 2013;31:1395-1401
51. Tsurui Y, Sho M, Kuzumoto Y, Hamada K, Akashi S, Kashizuka H. et al. Dual role of vascular endothelial growth factor in hepatic ischemia-reperfusion injury. Transplantation. 2005;79:1110-1115
52. Xu YM, Wang Y, Yan SY, Yang QH, Zhou YQ, Zeng XQ. et al. Regulation of endothelial intracellular adenosine via adenosine kinase epigenetically modulates vascular inflammation. Nat Commun. 2017 8
Corresponding authors: Wenchao Gao, MD, gaowenchao2007com; Xiaowen Liang, B.Sc., M.Pharm., Ph.D., Tel: 61-07-3443-7488; Email: x.liangedu.au