Abstract
PURPOSE:
Contrast enhanced magnetic resonance imaging (MRI) is able to assess liver function by characteristic changes of signal intensity (SI).
The aim was to evaluate dynamic contrast-enhanced SI-indices of the abdominal aorta, portal vein and liver.
METHODS:
72 patients underwent Gd-EOB-DTPA-enhanced MRI and a 13C-methacetin-based liver breath test (13C-MBT) for evaluation of liver function. Region-of-interest measurements in the liver, abdominal aorta and portal vein during native, arterial (AP), late arterial (LAP), portal venous (PVP) and hepatobiliary phase (HBP) were applied to analyze SI-indices in T1-weighted volume-interpolated breath-hold examination (VIBE) sequences with fat-suppression and relative enhancement (RE) analysis was performed.
RESULTS:
The liver (p < 0.001), the portal vein (p < 0.001) and abdominal aorta (p = 0.002) showed significant decrease of REs with decreasing liver function. An increasing trend between logarithmic values of 13C-MBT and REs of hepatic parenchyma (HBP; r = 0.662, p < 0.001), portal vein (PVP; r = 0.532, p < 0.001) and abdominal aorta (PVP; r = 0.421, p < 0.001) was observed.
CONCLUSIONS:
RE measurements of the hepatic parenchyma proofed to be a trustable evaluation method for liver function evaluation. In accordance with liver function, changes of REs in the portal vein and abdominal aorta occur.
Keywords
Introduction
Portal hypertension is one of the main implications of liver cirrhosis. It plays a key role in severe complications, such as rupturing of gastroesophageal varices, ascites, the occurrence of hepatorenal syndrome and hepatic encephalopathy. According to Ohm’s law, the blood pressure is dependent on flow and resistance, thus the portal pressure in portal hypertension is increased due to an increase in portal blood flow, an increase in vascular resistance or a combination of both. Hemodynamic changes are critical accessory phenomena of cirrhosis. The architectural distortion of the liver microcirculation caused by fibrosis, scarring and nodule formation as well as active contraction of the vascular smooth muscle cells, myofibroblasts and other contractile elements within or around the hepatic microcirculation have an impact on the intrahepatic vascular resistance. In cirrhotic patients, the development of portal hypertension will increase the risk of death and need for liver transplantation [3].
The liver has a unique dual blood supply which is based on venous blood from the portal vein and arterial blood from the hepatic artery. As there is communication between vessels, including transsinusoidal, transvasal and transplexal routes, vascular compromise will cause perfusion disorders. Established clinical tests to identify portal hypertension are based on different methods as the invasive balloon catheter method [2], analysis of xenobiotic-clearance [8, 15] or ultrasound measurements [4].
In addition to contrast enhanced ultrasound (CEUS)-based evaluation of parenchymal microcirculation [7], contrast enhanced MRI-based methods are able to assess hepatic flow parameters. Here, gadolinium-based contrast agent will be injected and changes in MR imaging will be monitored. Based on mathematical and physical hepatic perfusion models the data acquired during imaging will be fitted by algorithms and representative hepatic flow parameters, e.g. hepatic perfusion index or mean transit time will be assessed [10].
Gadoxetic acid (Gd-EOB-DTPA) is a hepatocyte-specific contrast agent that is capable to display liver function in an MRI-based test by analyzing the signal intensity caused by the interaction of gadolinium chelates with protons [13, 19]. As liver cirrhosis and fibrosis has an impact on the portal pressure, it may be possible to monitor changes using SI-based MR dynamic imaging of the portal vein and abdominal aorta.
The purpose of this study is to investigate the relationship between the liver function estimated by 13C-MBT and dynamic MRI-based SI values of the abdominal aorta, the portal vein and the liver parenchyma, respectively.
Material and methods
Patients
In this retrospective study, patients that received Gd-EOB-DTPA-enhanced T1-weighted volume-interpolated breath-hold examination (VIBE) MRI sequences with fat suppression at 3 T during native, arterial (AP), late arterial (LAP), portal venous (PVP) and hepatobiliary phase (HBP), were included. Further additional criteria were contemporary liver function evaluation by a 13C-MBT and absence of allergic reactions to liver-specific MRI contrast media, 13C-methacetin or specific contraindications to either MRI or Gd-EOB-DTPA administration (Primovist, Bayer Healthcare, Berlin, Germany). Informed consent was obtained from all participating patients and the study was conducted in full accordance with the ethical guidelines of the journal [1].
Magnetic resonance imaging
For T1-weighted VIBE sequences with fat suppression following conditions were applied: 3.09 ms repetition time, 1.17 and 2.49 ms echo time, 10° flip angle, parallel imaging factor 2 and 64 slices. The reconstructed voxel size was 1.25×1.25×3.0 m3 and the measured voxel size was 1.71×1.25×4.5 mm3. All images had an acquisition time of 14 sec and were acquired during breath-hold before Gd-EOB-DTPA administration (native phase) and during dynamic phases.
The patients received an i.v. bolus injected dose of Gd-EOB-DTPA (0.025 mmol/kg body weight; flow rate, 1 mL/sec) followed by 20 mL 0.9% sodium chloride for contrast-enhanced MRI.
Image analysis
The sequences were evaluated by region-of-interest (ROI) measurements. Circle shaped ROIs were manually placed at identical locations in every sequence, to obtain representative SI values while avoiding lesions and imaging artifacts. In total 8 ROIs (3 left liver lobe, 3 right liver lobe, 1 portal vein, 1 abdominal aorta) were defined per sequence.
The relative enhancement (RE) of the liver parenchyma (RE_LL), the portal vein (RE_PV) and the abdominal aorta (RE_AA) relative to the native phase was calculated for AP, LAP, PVP and HBP by following formulas:
a = AP, LAP, PVP, HBP
x = arterial, late arterial, portal venous, hepatobiliary
13C-methacetin breath test
The 13C-MBT was contemporary performed to the MRI scan (24 h) [16, 17]. Before the i.v. injection of 13C-methacetin, the 13CO2:12CO2 ratio control was recorded. After injection of 2 mg/kg body weight adapted 13C-methacetin and 20 mL 0.9% sodium chloride, a volatile analysis was performed by nondispersive isotope-selective infrared spectroscopy (FANci2-db16, Fischer Analysen Instrumente, Leipzig, Germany). According to the 13C-MBT readout values, the patients were classified: normal liver function (13C-MBT >315.0 [μg/kg/h]), intermediate liver function (13C-MBT: 140.0–315.0 [μg/kg/h]) and severely impaired liver function (13C-MBT <140.0 [μg/kg/h]) [17, 18].
Statistical analysis
The three 13C-MBT readout classes were compared as non-parametric independent samples by the Kruskal-Wallis-test. The predictive power of RE values were determined by Spearman-based bivariate correlation analysis, while in all tests the statistical significance level was set to 0.05 (two-sided). All analysis were performed using SPSS software (version 24; IBM, Chicago, IL, USA).
Results
Patient characteristics
In this retrospective study, patients received both a 13C-MBT and contrast enhanced dynamic T1-weighted MR imaging of the liver from native to hepatobiliary phase. 72 patients (51 men; 21 women) with a median age of 61.5 years were included (Table 1).
Patient characteristics of included cases
Patient characteristics of included cases
Arterial phase (AP), late arterial phase (LAP), portal venous phase (PVP) and hepatobiliary phase (HBP). Values indicate the mean±standard deviation. RE_LL: relative enhancement, a function of SI-based indices measured in the liver (LL), RE_PV: relative enhancement, a function of SI-based indices measured in the portal vein (PV), RE_AA: relative enhancement, a function of SI-based indices measured in the abdominal aorta, 13C-MBT: 13C-methacetin breath test.
The reasons for MR imaging were active monitoring of hepatocellular carcinoma (n = 17), follow-up reasons in the case of secondary malignancies and benign hepatic lesions (n = 13), preinterventional assessment in the case of HCC, known secondary liver malignancies or focal hepatic lesions (n = 11), postinterventional MR imaging (n = 17) and check-up for suspected liver disease or focal hepatic lesions (n = 14). The patients conditions were HCC with cirrhosis (n = 33), HCC without cirrhosis (n = 1), liver cirrhosis without HCC (n = 3), focal nodular hyperplasia (n = 1), carcinoma of the ileum (n = 1), carcinoma of the duodenum (n = 1), cholangiocarcinoma (n = 8), colon cancer (n = 1), hemangioma (n = 4), mamma carcinoma (n = 2), rectal cancer (n = 11), sigma carcinoma (n = 4), thymoma (n = 1) and uveal melanoma (n = 1).
A comparison of REs revealed that patients with normal liver function (13C-MBT >315.0 [μg/kg/h]) obtained the highest values in hepatic tissue, portal vein and abdominal aorta (Fig. 1). With dynamic MRI sequences, similar increase and decrease patterns, respecticvely, were observed in patients. While the RE of the hepatic tissue steadily increased with dynamic MRI sequences (Fig. 1A), the RE of the portal vein increased until late arterial phase before it started to decrease (Fig. 1B). In the abdominal aorta the highest values of the RE were obtained in the arterial phase (Fig. 1C).

Phase-dependent enhancement of relative enhancement (RE): hepatic parenchyma (RE LL; A), portal vein (RE PV; B) and abdominal aorta (RE AA; C). Cases were stratified in the distinct liver function groups normal liver function (light grey), intermediate liver function (dark grey) and severely impaired liver function (black-white stripes). arterial phase (AP), late arterial phase (LAP), portal venous phase (PVP) and hepatobiliary phase (HBP). Values indicate the mean±standard deviation. RE_LL: relative enhancement, measured in the liver parenchyma, RE_PV: relative enhancement measured in the portal vein (PV), RE_AA: relative enhancement measured in the abdominal aorta (AA), 13C-MBT: 13C-methacetin breath test, Level of significance: *p < 0.05; **p < 0.01; ***p < 0.001.
REs of the specific MRI sequences were tested against the 13C-MBT readout classes (Table 2). For the hepatic tissue, severely impaired liver function and normal liver function were significantly different (p < 0.002) in the late arterial phase (RE_LL_LAP, p < 0.001) and portal venous phase (RE_LL_PVP; p = 0.003). Additionally, in the hepatobiliary phase (RE_LL_HBP, p < 0.001) severely impaired liver function and intermediate liver function were also significant different (p < 0.001; Fig. 2A). For the portal vein, REs obtained in the arterial phase showed a significant difference between severely impaired and intermediate liver function (p < 0.001) as well as normal liver function (p < 0.001; Fig. 2B). Further significant differences were observed in the late arterial phase (RE_PV_LAP), portal venous (RE_PV_PVP) and hepatobiliary phase (RE_PV_HBP). The abdominal aorta was significant different between severely impaired liver function and normal liver function in the arterial phase (RE_AA_AP; p < 0.05), the late arterial phase (RE_AA_LAP; p < 0.01) and portal venous phase (RE_AA_PVP; p < 0.001; Fig. 2C).
Non-parametric test of SI-based indices between distinct liver function groups.
Non-parametric test of SI-based indices between distinct liver function groups.
REs of hepatic parenchyma, portal vein and abdominal aorta are tested based on 13C-MBT readout classification in a non-parametric Kruskal-Wallis test. The post-hoc tests reveal significant differences between the classes. Arterial phase (AP), late arterial phase (LAP), portal venous phase (PVP) and hepatobiliary phase (HBP). Values indicate the mean±standard deviation. RE_LL: relative enhancement, measured in the liver (LL), RE_PV: relative enhancement measured in the portal vein (PV), RE_AA: relative enhancement measured in the abdominal aorta (AA), 13C-MBT: 13C-methacetin breath test.

Boxplot analysis of SI-based indices separated by 13C-MBT readout categories. RE values obtained during the hepatobiliary phase of the liver (RE_LL_HBP; A) show significant differences between severely impaired liver function (median: 0.437) and intermediate (median: 0.901) as well as normal liver function (median: 1.069). RE values obtained in the portal vein during the late arterial phase (RE_PV_AP; B) are also significant different between severely impaired liver function (median: 1.036) and intermediate (median: 2.633) as well as normal liver function (median: 2.482). In contrast, RE values obtained in the abdominal aorta during the portal venous phase (RE_AA_PVP; C) show only between severely impaired (median: 1.537) and normal liver function (median: 2.336) a significant differentiation. 13C-MBT >315.0 [μg/kg/h]: normal liver function, 13C-MBT: 140.0–315.0 [μg/kg/h]: intermediate liver function, 13C-MBT <140.0 [μg/kg/h]: severely impaired liver function, 13C-MBT: liver function-dependent 13C-methacetin metabolism measured in a breath test, Level of significance: *p < 0.05; **p < 0.01; ***p < 0.001.
In a bivariate Spearman correlation analysis the REs of the tissues at specific time points were tested against logarithmic values of 13C-MBT readouts (Table 3). According to the Spearman correlation analysis, REs of the hepatic tissue had the highest correlation with logarithmic 13C-MBT readout values in the hepatobiliary phase (r = 0.648; p < 0.001; Fig. 3A). For portal vein, best correlations were reveived the portal venous phase, (r = 0.528; p < 0.001; Fig. 3B) and for abdominal aorta in the arterial phase (r = 0.405; p < 0.001; Fig. 3C).
Linear spearman correlation analysis of RE-values to logarithmic values of 13C-MBT
Linear spearman correlation analysis of RE-values to logarithmic values of 13C-MBT
Arterial phase (AP), late arterial phase (LAP), portal venous phase (PVP) and hepatobiliary phase (HBP). Values indicate the mean±standard deviation. RE_LL: relative enhancement, measured in the liver (LL), RE_PV: relative enhancement measured in the portal vein (PV), RE_AA: relative enhancement measured in the abdominal aorta (AA), 13C-MBT: 13C-methacetin breath test, r: correlation coefficient, n.s.: not significant.

Correlation analysis of SI-based indices to logarithmic values of 13C-MBT readout in scatterplots. The RE values of the hepatic tissue show the strongest correlation to logarithmic values of 13C-MBT readout (r = 0.648; A) compared to RE values of the portal vein (r = 0.528; B) and abdominal aorta (r = 0.405; C). RE_LL_HBP: RE values obtained during the hepatobiliary phase of the liver, RE_PV_PVP: RE values obtained during the portal venous phase of the portal vein, RE_AA_PVP: RE values obtained during the portal venous phase of the abdominal aorta, 13C-MBT: liver function-dependent 13C-methacetin metabolism measured in a breath test.
Increased portal hypertension is associated with severe clinical complications of liver cirrhosis like variceal bleeding, ascites or portal hypertensive gastropathy. According to anatomical location, the aetiology of increased portal resistance can be categorized into pre-hepatic, intra-hepatic and post-hepatic causes [2]. While cirrhotic portal hypertension is associated with an increased hepatic venous pressure gradient, in non-cirrhotic patients suffering from portal hypertension, the hepatic venous pressure gradient is normal or only mildly elevated but lower than the portal venous pressure [9]. Architectural changes influence the sinusoidal blood flow and cause an increase of hepatic resistance. This results in an increased splanchnic vasodilation and increased blood flow into the portal venous system and initial changes regarding hemodynamic and microcirculation in liver parenchyma might be detected [6]. Similar to other hepatic disorders, the initial presence of portal hypertension occurs almost unremarkable. Nevertheless, the assessment of the hepatic venous pressure gradient can be a strong indication for portal hypertension. Accurate measurements of the hepatic venous pressure gradient rely on invasive techniques. Therefore a balloon catheter needs to be introduced into the hepatic vein to allow a direct measurement [2]. Another option to access hepatic perfusion parameters is by analyzing the flow limited clearance of xenobiotics like indocyanine green or sorbitol [8, 15]. But, as in cases of advanced liver disease the perfusion as well as intracellular enzymatic and transport processes are altered, so the application of xenobiotics for hepatic perfusion assessment may not represent the perfusion status accurate enough. Another option is to assess hepatic perfusion parameters by a hepatocyte model for pharmacokinetic analysis which is based on dynamic liver MR imaging over time [10].
The radiological imaging of Gadoxetic-acid in MRI-based functional imaging is an established tool for assessing liver function and displaying hepatic lesions. One convincing benefit of the contrast agent is the resistance to enzymatic metabolism. In pioneer-like studies, the use of Gadolinum-based dynamic MRI for the noninvasive assessment of hepatic perfusion parameters was evaluated [11]. They showed that MRI-based compartmental analysis of signal intensity versus time curves can be applied as alternative method to measure hemodynamic parameters, as hepatic flow parameters correlated with the severity of cirrhosis and portal hypertension.
Liver cirrhosis is associated with a loss of functioning hepatocytes, bridging of portal spaces or nodular regeneration of the liver parenchyma, which hinders hepatocyte contrast agent uptake and causes decreased SI-based values [12, 14]. The estimation of SI-based indices provides insights into regional hepatocyte-specific function, efficiency, functionality, condition and perfusion.
In order to reduce the number of different examinations and to reduce ionizing radiation dose, researchers evaluate established tests in order to extract further information about patients’ health status. Different studies proofed that hepatic flow parameters achieved by contrast agent-based MR images are able to display liver function. We hypothesized that T1-weighted contrast-enhanced SI-based MR imaging of the portal vein and the abdominal aorta will be able to display liver function. Therefore we analyzed time-dependent RE values and compared them to liver function, represented by 13C-MBT readout values.
Our results demonstrate that SI-based MR-imaging of the portal vein and abdominal aorta show a characteristic pattern which is determined by liver function and time-dependent image acquisition. Similar to Nilsson et al. [13], we observed that patients with liver disease achieve maximum RE after a shorter time than patients with better liver function. Our data permit a crude assumption on hepatic function already in the portal venous phase, based on MR imaging of liver parenchyma, as patients with severely impaired liver function show now further enhancement. In accordance to the finding that hepatic perfusion index is dependent on fibrosis severity made by Leporq et al. [10], we observed that dynamic SI-indices of the liver, portal vein and abdominal aorta are decreased with increasing progression of liver damage. In contrast to our measurements, Leporq et al. used a dual-input three-compartment hepatocyte model for pharmacokinetic analysis. We in contrast expected to see an effect on SI-values without applying a perfusion model, as SI-indices depend indirect on contrast agent concentration, which is influenced by aortic flow and pressure. Our results revealed a significant correlation between SI-based indices of hepatic tissue, portal vein (r = 0.528; p < 0.001) and abdominal aorta (r = 0.405; p < 0.001) and 13C-MBT readout values, which correlates to Child-Pugh classes (r = –0.707; p < 0.001) [5].
A recent published study by Wildner et al. [20] demonstrated a significant dependence of angiogenesis parameters assessed by dynamic contrast-enhanced ultrasound (DCEUS) to specific malignant liver lesions. Their data showed that benign lesions (focal-nodular-hyperplasia) and malignant focal liver lesions (hepatocellular carcinoma, cholangiocellular carcinoma, pancreatic adenocarcinoma, breast cancer, colorectal cancer and melanoma) experience specific absolute signal intensity values and characteristic washout. While the focus of our study was set on MRI-based signal intensity values of portal vein and abdominal aorta, an additional analysis based on specific hepatic lesion and DCEUS assessed angiogenesis parameters would improve the informative value.
Our study has several limitations. In order to decrease the influence of nonhomogeneous distributed parenchymal changes, analysis were done with the mean of six distinct ROI measurements in the case of hepatic parenchyma. Another limitation is the absence of healthy patients without hepatic diseases. Also, an analysis based on specific hepatic diseases and hepatic venous pressure gradient would be beneficial, as they might be crucial parameters for portal hypertension which now are not taken into account. An analysis of relaxometry values might add further information, which may improve the evaluation of the abdominal aorta and portal vein.
In conclusion, evaluation of dynamic SI-based indices of liver parenchyma, portal vein and abdominal aorta can give a strong hint on liver function with best results for liver parenchyma evaluation.
