## Abstract

The purpose of this study is to predict preoperatively microvascular invasion (MVI) of solitary small hepatocellular cancer (sHCC) by using the kinetic parameters analysis on dynamic enhancement magnetic resonance imaging (MRI). Patients (${\displaystyle n=61}$) with known solitary sHCC(${\displaystyle \leq 3}$cm) were preoperatively examined with Gd-EOB-DTPA-enhanced MRI first before hepatic resection. The arterial peritumoral enhancement measured from the dynamic enhancement-MRI was analyzed by using quantitative kinetic parameters, including initial enhancement (${\displaystyle E_{1}}$), peak enhancement (${\displaystyle E_{\rm {peak}}}$), and enhancement ratio (${\displaystyle E_{R}}$) calculated. Correlations between quantitative kinetic parameters and MVI were evaluated and differences between MVI positive and negative groups were assessed. Histopathological analysis of liver resection confirmed that 19 patients had sHCC with MVI and that 42 patients had sHCC without MVI. Average (${\displaystyle \pm }$ standard deviation) ${\displaystyle E_{1}}$ is ${\displaystyle 0.36\pm 0.12}$ and ${\displaystyle 0.46\pm 0.09}$, ${\displaystyle E_{\rm {peak}}}$ is ${\displaystyle 0.78\pm 0.24}$ and ${\displaystyle 0.74\pm 0.18}$, and ${\displaystyle E_{R}}$ is ${\displaystyle 0.42\pm 0.20}$ and ${\displaystyle 0.56\pm 0.17}$ for negative and positive group, respectively. Statistical analysis showed that average ${\displaystyle E_{1}}$ and ER for the positive group were significantly higher (${\displaystyle p<0.05}$) than the negative group. The receiver operating characteristics (ROC) analysis between the two groups had area under the curve of 0.74 and 0.71 for ${\displaystyle E_{1}}$ and ${\displaystyle E_{R}}$, respectively. Quantitative kinetic parameters analysis for the arterial peritumoral enhancement is feasibility to the prediction and assist diagnosis of MVI in clinical practice.

Keywords: Microvascular invasion, solitary small hepatocellular cancer, quantitative analysis, kinetic parameters, dynamic enhancement magnetic resonance imaging

## 1. Introduction

Microvessel invasion (MVI) is a major prognostic factor in hepatocellular carcinoma (HCC) that influences the choice of treatment, but rarely can be evaluated preoperatively [1,2]. Several studies have reported that certain imaging findings on dynamic enhancement MRI are useful for predicting MVI of HCC, including tumor size [3], peritumoral enhancement [4], tumor margin [5], tumor hypointensity or peritumoral hypointensity on hepatobiliary phase (HBP) [6,7], radiological capsule on gadoxetic acid–enhanced MR imaging and so on [8].

Gd-EOB-DTPA enhanced MRI can provide more valuable information for the assessment of HCC and has been widely used in preoperative evaluation settings [9,10]. Some researches use the HBP images to describe the tumor margins in the axial and coronal hepatobiliary phase, including smooth [6], wedge-shaped peritumoral enhancement [11] and irregular circumferential peritumoral enhancement [12].

However, all of these studies are based on qualitative observations and lack the kinetic parameters analysis of arterial phase tumor peritumoral enhancement on MVI. Especially for the small hepatocellular carcinoma (sHCC) with the maximum tumor diameter ${\displaystyle \leq 3}$cm [13], there are limited studies that have investigated MR imaging finding for predicting MVI [11]. Therefore, techniques for quantitative analysis of dynamic enhancement MRI should be developed to directly predict MVI of sHCC, which results in early recurrence after hepatic resection and poor prognosis [13].

To our knowledge this is the first study to quantitatively preoperatively predict solitary sHCC with MVI based on kinetic parameters analysis of dynamic enhancement MRI. The kinetic analysis of pre-enhanced phase, arterial phase (20 seconds), portal venous phase (55 seconds), equilibrium phase (90 seconds) and delayed phase (180 seconds) enhancement curves of dynamic enhancement MRI were used in this study.

## 2. Materials and methods

This study was approved by the Institutional Review Board at our institution and the requirement for informed consent was waived. Patients were recruited from February 2018 to October 2019. The dynamic enhancement MRI was performed within one week before curative hepatectomy, and there is no macroscopic vascular invasion on MRI. A total of 61 patients (average age = 56 years old; 15 females and 46 males) with known single sHCC (${\displaystyle \pm 3}$cm) were enrolled in this study. Tumor size, number, and capsule condition were obtained at gross specimen. The entire tumor was examined for sHCC. MVI was defined as the presence of tumor in vascular space of surrounding hepatic tissue line by endothelial cells on microscopy [14].

MRI examination was performed by using 3.0-T system (Ingenia, Philips Medical Systems, Eindhoven, and The Netherlands) in all patients. 32-channel phased-array coil was used. The scanning scale covered from the top to the lower edge of live. The protocol used consisted of the following sequences which were shown in Table 1.

Table 1. The protocol used consisted of the following sequences
 Sequences Index Value Axial fat-suppressed RT T2W single shot turbo spin echo TR/ TE ${\displaystyle 535/75}$ slice thickness /gap ${\displaystyle 7/1}$ mm FOV ${\displaystyle 350\times 392}$ cm matrix size ${\displaystyle 232\times 199}$ Coronal breath-hold T2W single shot turbo spin echo TR/TE ${\displaystyle 1100/80}$ slice thickness /gap ${\displaystyle 6/1}$ mm FOV ${\displaystyle 350\times 346}$ cm matrix size ${\displaystyle 292\times 253}$ Axial breath-hold dual-echo(in-phase and opposed-phase) T1W fast field-echo TR/TE1/TE2 ${\displaystyle 106/1.15/2.3}$ slice thickness /gap ${\displaystyle 7/1}$ mm FOV ${\displaystyle 400\times 322}$ cm matrix size ${\displaystyle 244\times 181}$ Fat-suppressed dynamic 3D volumetric interpolated breath-hold T1W sequence (before /after injection) TR/TE1/TE2 ${\displaystyle 3.6/1.32/2.3}$ slice thickness /gap ${\displaystyle 5/-2.5}$ mm FOV ${\displaystyle 320\times 427}$ cm matrix size ${\displaystyle 200\times 250}$ Injection of Gd-EOB-DTPA Primovist/ flow rate ${\displaystyle 0.1}$ mL/kg ${\displaystyle 1}$ mL/s saline/ flow rate ${\displaystyle 20}$ml ${\displaystyle 2}$ml/s The dynamic enhanced sequence pre-enhanced phase - arterial phase ${\displaystyle 20}$s after injection portal venous phase ${\displaystyle 55}$ s after injection equilibrium phase ${\displaystyle 90}$s after injection delayed phase ${\displaystyle 180}$s after injection Hepatobiliary phase images - ${\displaystyle 20}$min after injection TR, repetition time; TE, echo time; FOV, field of view; RT, respiratory-triggered

For each patient, the tumor ROI was evaluated by the radiologist (15-year abdominal diagnosis experience). Dynamic enhancement MRI data was analyzed using MATLAB (MathWorks, Natick, MA) with an in-house software package. The peritumoural region (dilated distance was analysed as dilated pixel counts multiplied by pixel size) of the labeled ROI was dilated at a radius of 1 cm automatically by topology algorithm. The peritumoural region in this paper is 8 pixels thick. Then, all the pixels on the peritumoural region were used for kinetic analysis.

ROI Topology Algorithm has follow steps:

Step 1: Search the ROI pixels (value=1);

Step 2: if the product of adjacent pixels value equal 0, that means the pixel is an edge;

Step 3: let the adjacent pixel equal 1 to enlarge the ROI with 1 pixel thick;

So, one loop can make the ROI enlarge 1 pixel thick, and in our research the algorithm need loop 8 times.

For the kinetic analysis, the signal intensity curve (${\displaystyle S(t)}$) were utilized from the dynamic enhancement MRI data [15]. For each patient, average peritumoural region enhancement curve was used to calculate the kinetic parameters. The baseline signal intensity value (${\displaystyle S_{0}}$) of each pixel was the signal intensity (${\displaystyle S(t)}$) of the precontrast time point. Three quantitative kinetic parameters, namely initial enhancement (${\displaystyle E_{1}}$), peak enhancement (${\displaystyle E_{\rm {peak}}}$), and enhancement ratio (${\displaystyle E_{R}}$) were calculated using ${\displaystyle S(t)}$ as follows [16]:

 ${\displaystyle E_{1}=(S_{1}-S_{0}/S_{0}}$,
(1)
 ${\displaystyle E_{\rm {peak}}=(S_{\rm {peak}}-S_{0})/S_{0}}$,
(2)
 ${\displaystyle E_{R}=(S_{1}-S_{0})/(S_{\rm {last}}-S_{0})}$,
(3)

where ${\displaystyle S_{0}}$, ${\displaystyle S_{1}}$, ${\displaystyle S_{\rm {peak}}}$, and ${\displaystyle S_{\rm {last}}}$ are the signal intensity at the baseline, the 1st enhancement time point, the peak point and the last time point, respectively. The detailed data analysis steps are shown in Figure 1 flowchart.

 Figure 1. The flowchart of detailed steps used in the kinetic analysis

Student t-tests were performed to examine whether there were significant differences between MVI positive and negative groups for calculated quantitative parameters. Receiver operating characteristics (ROC) analysis was performed to evaluate whether calculated parameters could be used to distinguish between MVI positive and negative groups. A p-value less than 0.05 was considered as statistically significant.

## 3. Results

### 3.1. Patients

For all 61 patients, they were diagnosed as a MVI positive group (${\displaystyle n=19}$) if the invasion of tumor cells within a vascular space lined by endothelium that is visible only on microscopy, and otherwise as an MVI negative group (${\displaystyle n=42}$).

Figure 2 shows the peritumoral region and the corresponding signal intensity curve ${\displaystyle S(t)}$ of a MVI positive patient and a negative one for example. The curve demonstrated that there was clearly higher enhancement for the MVI positive group compared to the MVI negative group.

 Figure 2. The peritumoural region and the corresponding signal intensity curve ${\displaystyle S(t)}$ of MVI positive patient, male, 47 years old (a and b) and negative patient, male, 44 years old (c and d). The curve demonstrated that there was clearly faster enhancement (E1) and washout (ER) for the MVI positive group compared to the MVI negative group

### 3.2. Kinetic analysis

Figure 3 shows boxplots of kinetic parameters for both MVI positive and negative groups: (a) initial enhancement ${\displaystyle E_{1}}$, (b) peak enhancement ${\displaystyle E_{\rm {peak}}}$, and (c) enhancement ratio ${\displaystyle E_{R}}$. Average (${\displaystyle \pm SD}$) ${\displaystyle E_{1}}$ is ${\displaystyle 0.36\pm 0.12}$ and ${\displaystyle 0.46\pm 0.09}$, ${\displaystyle E_{\rm {peak}}}$ is ${\displaystyle 0.78\pm 0.24}$ and ${\displaystyle 0.74\pm 0.18}$, and ${\displaystyle E_{R}}$ is ${\displaystyle 0.42\pm 0.20}$ and ${\displaystyle 0.56\pm 0.17}$ for negative and positive group, respectively. Statistical analysis showed that average ${\displaystyle E_{1}}$ and ${\displaystyle E_{R}}$ for the positive group were significantly higher (${\displaystyle p<0.05}$) than the negative group.The ${\displaystyle E_{\rm {peak}}}$ for positive group was slightly lower (${\displaystyle p=0.49}$) than the negative group.

 Figure 3. Boxplots of parameters for both negative and positive groups: (a) ${\displaystyle E_{1}}$ and (b) ${\displaystyle E_{R}}$. The square (□) indicates mean of the data. The ${\displaystyle E_{1}}$ value for negative group and positive group were ${\displaystyle 0.36\pm 0.12}$ and ${\displaystyle 0.46\pm 0.09}$, the ${\displaystyle E_{R}}$ value ${\displaystyle 0.42\pm 0.20}$ and ${\displaystyle 0.56\pm 0.17}$ for negative and positive group, respectively. Statistical analysis showed that average ${\displaystyle E_{1}}$ and ${\displaystyle E_{R}}$ for the positive group were significantly higher (${\displaystyle p<0.05}$) than the negative group

Finally, Figure 4 shows the ROC analysis results for the parameters ${\displaystyle E_{1}}$ and ${\displaystyle E_{R}}$ obtained from kinetic analysis parameters. By selecting the Youden index cut-off point, the calculated sensitivity is 0.57 and 0.52; and specificity is 0.75 and 0.75 for the ${\displaystyle E_{1}}$ and ${\displaystyle E_{R}}$, respectively. The corresponding cut-off value is 0.44 and 0.58 for the ${\displaystyle E_{1}}$ and ${\displaystyle E_{R}}$, respectively.

 Figure 4. The receiver operating characteristics (ROC) analysis results between MVI negative and MVI positive group for the parameters ${\displaystyle E_{1}}$ and ${\displaystyle E_{R}}$. The same color dot on the same color line is indicating the Youden index cut-off point

## 4. Discussion

We investigated the feasibility of using dynamic enhancement MRI to directly assess the peritumoral enhancement for prediction the MVI of sHCC. The dynamic enhancement MRI included pre-enhanced phase, arterial phase, portal venous phase, equilibrium phase and delayed phase, which can help us track the blood supply. The kinetic analysis for dynamic enhancement MRI could be more sensitive to imperceptible changes. Furthermore, the enhancement could be analyzed with the physiological compartment model or mathematical model to obtain the qualitative results of slight changes in the lesion hemodynamics, so as to help doctors make accurate diagnosis. The downside of using enhancement curve is that the dynamic enhancement MRI needs to be acquired with more time points, which may add to the burden of clinical work. However, dynamic enhanced liver examination has become a routine method for clinical diagnosis. We believe that kinetic analysis can bring more reliable information for clinical diagnosis.

Our results demonstrated that the kinetic analysis of peritumoral enhancement measured from the dynamic enhancement MRI could be useful for quantitatively evaluating the MVI of sHCC. The parameters ${\displaystyle E_{1}}$ and ${\displaystyle R}$ are higher in the positive group than negative group, which means there are faster blood supply and faster washout in the liver parenchyma around the tumor. The feature indicates there is the same kinetic characteristic in the peritumoral region as intratumoral of HCC which increases the probability of MVI. These agreements were very encouraging to inspire more studies in the future. We believe that this is the first study to quantitatively evaluate peritumoral enhancement measured from the dynamic enhancement MRI for prediction the MVI of sHCC.

Previous studies of predicting the MVI of sHCC were mainly on tumor morphology by evaluating the smooth of tumor margin [6,17], the regularity of tumor shape [18], the completeness of tumor capsule [19]. There is a lack of relevant quantitative research.

The tumor morphology may be not equal to the MVI of sHCC for several reasons. First, the judgment of sHCC from morphology is limited [20,21]. Second, some image features due to individual blood supply differences may not always appear in the corresponding acquisition phase. Third, these features are mainly evaluated by the naked eye of experienced radiologist. Our quantitative analysis of dynamic enhancement MRI could be used to overcome above problems to some extent by quantifying the prediction the MVI of sHCC.

A potential advantage of dynamic enhancement MRI would be the possibility to perform one-stop examination just before HBP image. As a routine clinical examination, dynamic enhancement data are relatively easy to obtain.

There were several limitations to this study. First, the scale of our research was relatively small, especially the positive group. Second, the study did not include comparisons of blood supply changes between different contrast agents. Third, the sHCC in the study is restricted to diameter not bigger than 3cm. The HCC with diameter between 3 and 5cm should also be predicted in the future. Finally, our developed techniques were not tested on different datasets to evaluate the sensitivity and specificity obtained in this study.

There is potential clinical application using kinetic analysis for dynamic enhancement MRI in assisting the MVI of sHCC. Quantitative analysis parameters could be used as objective measures in diagnosis. Our study would inspire more research on the prediction the MVI for the solitary sHCC patients in the future.

## Funding

This research was funded by Harbin youth reserve talent project (2017RAQXJ102) and Harbin Medical University Cancer Hospitals Haiyan Funds (JJZD2020-17).

## Conflicts of Interest

The authors declare no conflict of interest.

## References

[1] Imai K., Yamashita Y.I., Yusa T., et al. Microvascular invasion in small-sized hepatocellular carcinoma: significance for outcomes following hepatectomy and radiofrequency ablation. Anticancer Res., 38(2):1053-1060, 2018.

[2] Ahn S.J., Kim J.H., Park S.J., Kim S.T., Han J.K. Hepatocellular carcinoma: preoperative gadoxetic acid-enhanced MR imaging can predict early recurrence after curative resection using image features and texture analysis. Abdom. Radiol., 44(2):539-548, 2019.

[3] Yamashita Y.I., Imai K., Yusa T., et al. Microvascular invasion of single small hepatocellular carcinoma ≤3 cm: Predictors and optimal treatments. Ann. Gastroenterol. Surg., 2(3):197-203, 2018.

[4] Kim K.A., Kim M.J., Jeon H.M., et al. Prediction of microvascular invasion of hepatocellular carcinoma: usefulness of peritumoral hypointensity seen on gadoxetate disodium-enhanced hepatobiliary phase images. J Magn. Reson. Imaging, 35(3):629-34, 2012.

[5] Ariizumi S., Kitagawa K., Kotera Y., et al. A non-smooth tumor margin in the hepatobiliary phase of gadoxetic acid disodium (Gd-EOB-DTPA)-enhanced magnetic resonance imaging predicts microscopic portal vein invasion, intrahepatic metastasis, and early recurrence after hepatectomy in patients with hepatocellular carcinoma. J. Hepatobiliary Pancreat. Sci., 18(4):575-85, 2011.

[6] An C., Rhee H., Han K., et al. Added value of smooth hypointense rim in the hepatobiliary phase of gadoxetic acid-enhanced MRI in identifying tumour capsule and diagnosing hepatocellular carcinoma. Eur. Radiol., 27(6):2610-2618, 2017.

[7] Choi Y.S., Rhee H., Choi J.Y., et al. Histological characteristics of small hepatocellular carcinomas showing atypical enhancement patterns on gadoxetic acid-enhanced MR imaging. J. Magn. Reson. Imaging, 37(6):1384-1391, 2013.

[8] Grazioli L., Olivetti L., Fugazzola C., et al. The pseudocapsule in hepatocellular carcinoma: correlation between dynamic MR imaging and pathology. Eur. Radiol., 9(1):62-7, 1999.

[9] Choi J.W., Lee J.M., Kim S.J., et al. Hepatocellular carcinoma: imaging patterns on gadoxetic acid-enhanced MR Images and their value as an imaging biomarker. Radiology, 267(3):776-86, 2013.

[10] Joo I., Lee J.M. Recent Advances in the imaging diagnosis of hepatocellular carcinoma: value of gadoxetic acid-enhanced MRI. Liver Cancer, 5(1):67-87, 2016.

[11] Kim H., Park M.S., Choi J.Y., et al. Can microvessel invasion of hepatocellular carcinoma be predicted by pre-operative MRI. Eur. Radiol., 19(7):1744-1751, 2009.

[12] Huang M., Liao B., Xu P., et al. Prediction of microvascular invasion in hepatocellular carcinoma: preoperative Gd-EOB-DTPA-dynamic enhanced MRI and histopathological correlation. Contrast Media Mol. Imaging, 9674565, pp. 9, 2018.

[13] Du M., Chen L., Zhao J., et al. Microvascular invasion (MVI) is a poorer prognostic predictor for small hepatocellular carcinoma. BMC Cancer, 14:38, 2014.

[14] Feng L.H., Dong H., Lau W.Y., et al. Novel microvascular invasion-based prognostic nomograms to predict survival outcomes in patients after R0 resection for hepatocellular carcinoma. J. Cancer Res. Clin. Oncol., 143(2):293-303, 2017.

[15] Speidel Ma, Bateman C.l., Tao Y., et al. Reduction of image noise in low tube current dynamic CT myocardial perfusion imaging using HYPR processing: a time-attenuation curve analysis. Med. Phys., 40:011904, 2013.

[16] Zhou X.Y., Zhang D.M. et al. Quantitative analysis of lower leg muscle enhancement measured from dynamic computed tomographic angiography for diagnosis of peripheral arterial occlusive disease. Journal of Computer Assisted Tomography, 44(1):20-25, 2020.

[17] Chen Z.H., Zhang X.P., Wang H., et al. Effect of microvascular invasion on the postoperative long-term prognosis of solitary small HCC: a systematic review and meta-analysis. HPB (Oxford), 21(8):935-944, 2019.

[18] Lee S., Kim S.H., Lee J.E., Sinn D.H., Park C.K. Preoperative gadoxetic acid-enhanced MRI for predicting microvascular invasion in patients with single hepatocellular carcinoma. J. Hepatol., 67(3):526-534, 2017.

[19] Wu D., Tan M., Zhou M., et al. Liver computed tomographic perfusion in the assessment of microvascular invasion in patients with small hepatocellular carcinoma. Invest. Radiol., 50(4):188-94, 2015.

[20] Yang C.B., Zhang S., Jia Y.J., et al. Dual energy spectral CT imaging for the evaluation of small hepatocellular carcinoma microvascular invasion. Eur. J. Radiol., 95:222-227, 2017.

[21] El-Gendi A., El-Shafei M., Abdel-Aziz F., Bedewy E. Intraoperative ablation for small HCC not amenable for percutaneous radiofrequency ablation in Child A cirrhotic patients. J. Gastrointest. Surg., 17(4):712-718, 2013.

### Document information

Published on 24/03/20
Accepted on 23/03/20
Submitted on 05/03/20

Volume 36, Issue 1, 2020
DOI: 10.23967/j.rimni.2020.03.003

### Document Score

0

Views 17
Recommendations 0