Intravoxel incoherent motion model-based analysis of diffusion-weighted magnetic resonance imaging with 3 b-values for response assessment in locoregional therapy of hepatocellular carcinoma.

PURPOSE: The aim of this study was to evaluate an intravoxel incoherent motion (IVIM) model-based analysis of diffusion-weighted imaging (DWI) for assessing the response of hepatocellular carcinoma (HCC) to locoregional therapy. PATIENTS AND METHODS: Respiratory-gated DWI (b=0, 50, and 800 s/mm2) was retrospectively analyzed in 25 patients who underwent magnetic resonance imaging at 1.5 T before and 6 weeks following the first cycle of transarterial chemoembolization therapy, transarterial ethanol-lipiodol embolization therapy, and transarterial radioembolization therapy. In addition to the determination of apparent diffusion coefficient, ADC(0,800), an estimation of the diffusion coefficient, D', and the perfusion fraction, f', was performed by using a simplified IVIM approach. Parameters were analyzed voxel-wise. Tumor response was assessed in a central slice by using a region of interest (ROI) covering the whole tumor. HCCs were categorized into two groups, responders and nonresponders, according to tumor size changes on first and second follow ups (if available) and changes of contrast-enhanced region on the first follow up.
RESULTS: In total, 31 HCCs were analyzed: 17 lesions were assigned to responders and 14 were to nonresponders. In responders, ADC(0,800) and D' were increased after therapy by ~30% (P=0.00004) and ~42% (P=0.00001), respectively, whereas f' was decreased by ~37% (P=0.00094). No significant changes were found in nonresponders. Responders and nonresponders were better differentiated by changes in D' than by changes in ADC(0,800) (area under the curve =0.878 vs 0.819 or 0.714, respectively).
CONCLUSION: In patients with HCCs undergoing embolization therapy, diffusion changes were better reflected by D' than by conventional ADC(0,800), which is influenced by counteracting perfusion changes as assessed by f'.

Introduction

Despite routine clinical use of morphology-based tumor response evaluation, there is a growing consensus that size evaluation alone has certain limitations because therapy-induced changes in tumor volume often occur relatively late after treatment.1,2 Functional imaging parameters have been proposed to improve the accuracy of early response assessment.1 Diffusion-weighted imaging (DWI) is a functional imaging technique providing information on tissue cellularity, extracellular space tortuosity, and integrity of cell membranes by measuring random motion of water molecules without injection of exogenous contrast medium or radiation exposure.1,2 DWI is suited for early response evaluation by detecting loss of cellularity due to necrosis.1–3 Therapy-induced changes detected by using DWI in hepatocellular carcinomas (HCCs) after transcatheter arterial chemoembolization (TACE) showed stronger correlation with histopathological degree of necrosis than those from native or contrast-enhanced imaging according to response evaluation criteria in solid tumors (RECIST), European Association for the Study of the Liver (EASL), or modified RECIST (mRECIST).4 Furthermore, DWI improved the detection of viable tumor tissue in HCCs following radiotherapy, when added to conventional magnetic resonance imaging (MRI).5

Quantitative analysis of DWI data was conducted from a set of images with varying degrees of diffusion weighting as selected by the diffusion weighting factors (b-values) of the sequence.1,2 For oncologic imaging, the apparent diffusion coefficient (ADC) is commonly analyzed from ≥2 b-values between 0 and 500–1,000 s/mm2, assuming a monoexponential dependence of signal intensity on the b-value.1,2 An increase in the ADC is related to therapy-induced necrosis.1–3 On the other hand, it is well known that DWI is sensitive not only to molecular diffusion but also to pseudorandom movements such as blood flow in the randomly oriented capillary network.2 Perfusion leads to additional signal intensity attenuation at low b-values, whereas the signal intensity attenuation at higher b-values is mainly caused by true diffusion.2 Thus, the ADC is a composite coefficient affected by both perfusion and diffusion. For DWI acquired with >2 b-values, a refined DWI analysis based on intravoxel incoherent motion (IVIM) theory was proposed by Le Bihan.6 In this analysis, a biexponential behavior of the signal intensity was assumed, which enables the separation of the effects of diffusion and perfusion by determining true diffusion coefficient, D, pseudodiffusion coefficient, D*, and perfusion fraction, f.7 The true diffusion coefficient D represents the mobility of water in tissues, the perfusion fraction f reflects the relative contribution of microvascular perfusion to the DWI signal, and the pseudodiffusion coefficient D* is related to the blood flow velocity and vessel architecture.7,8 In a normal liver tissue, an accurate determination of D* by nonlinear least squares fitting procedures with numerous b-values is required, if changes in blood flow velocity should be detected as in diffuse liver disease.9–13 However, in the liver, a simultaneous determination of D, f, and D* by using nonlinear biexponential fitting is challenging due to motion-induced signal variations and low signal-to-noise ratio (SNR).7,14 In malignant liver lesions, the application of such fitting procedures is even more problematic, as D* and f values in malignant liver lesions are lower than in normal liver tissue, leading to numerical instabilities and unreliable results.15–17 Thus, for oncological applications, dedicated IVIM approaches should be developed, which allow a numerically stable estimation of IVIM parameters, at least D and f, in settings with low D* and/or low f values. The estimation of these parameters should be performed voxel-wise for spatially resolved assessment and within clinically acceptable acquisition times in the interest of patient. Even if such estimates may deviate from exact D, D*, and f values, they serve as useful empirical biomarkers for diffusion and perfusion changes.

Therefore, the aim of the present study was to evaluate a simplified IVIM model–based DWI analysis based on 3 b-values for response assessment in patients with HCCs undergoing locoregional therapy.

Patients and methods

Study population

This retrospective study includes data from clinical routine examinations. The DWI sequence was part of the standardized abdominal imaging protocol since May 2008. This retrospective study was approved by the local institutional review board of the University Hospital Bonn, which also waived the need for informed patient consent because it was a retrospective re-analysis study of clinically routine images. Patients with HCC (confirmed histologically or by using typical MRI findings according to EASL criteria), who underwent TACE, transarterial ethanol–lipiodol embolization (TAELE), or radioembolization (RE) and who underwent MRI before and 6 weeks following the first cycle of therapy (the first follow up), were included in the study. Patients were excluded if they were treated with any other local or systemic chemotherapy or radiotherapy before baseline or the first follow up MRI or if the quality of diffusion-weighted image was not sufficient due to motion artifacts.

The final study population is 25 patients (5 female and 20 male mean age =66±11 years; range =40–83 years), treated with TACE (n=12), TAELE (n=2), and RE (n=11). For TACE, Doxorubicin (AdriblastinHL™; Pharmacia, Erlangen, Germany) in combination with ethiodized oil (Lipiodol Ultra-Fluid, Guerbet, France), degradable starch microspheres (EmboCept® S; PharmaCept, Berlin, Germany), or tandem beads (Embozene TANDEM™ Microspheres; CeloNova BioSciences, San Antonio, TX, USA) were used. For TAELE, ethanol 95% and ethiodized oil (Lipiodol Ultra Fluid) were injected. For RE, SIR-Spheres (SIRTEX Medical Limited, North Sydney, Australia) were used in 7 patients and TheraSphere (BTG International, Canada) in 4 patients. In 18 patients the right liver lobe was treated, in 4 patients the left liver lobe, and in 3 patients both the liver lobes. Mean time between pretherapeutic MRI and therapy was 13±15 days (range =0–47 days) and that between therapy and first follow up MRI was 44±21 days (range =21–113 days). In 22 of 25 patients, a second follow up MRI was available (mean time between therapy and second follow up MRI was 104±28 days; range =51–177 days), 3 patients had survived <3 months or were too ill to perform a second follow up. Posttherapeutic survival was monitored in 20 patients; 4 patients were alive at the date of analysis (>1,000 days since intervention).

Magnetic resonance imaging

All the examinations were performed on clinical 1.5-T MRI scanners (Gyroscan Intera and ACS-NT, 1.5 T; Philips Healthcare, Best, the Netherlands; gradient system: maximum amplitude =30 mT/m, maximum slew rate =150 T/m/s) by using a commercially available 4-element phased array surface receiver coil with the patient in supine position. In each patient, baseline and follow up MRI examinations were performed by using the same MRI scanner. The DWI sequence (Table 1), a respiratory-triggered single-shot spin-echo echo-planar imaging variant with 3 b-values (b0=0, b1=50, and b2=800 s/mm2) and motion-probing gradients in 3 orthogonal directions, was always acquired before contrast agent injection. Isotropic diffusion-weighted images were reconstructed online on the MRI system. The protocol also included a T2-weighted sequence with and without fat suppression, a T1-weighted sequence, and a T1-weighted dynamic contrast-enhanced sequence.

Table 1

Diffusion-weighted imaging sequence parameters

Name	Value
FOV (RL × AP)/orientation	380×326 mm/transversal
Slice number/thickness/gap	28/7.0 mm/0.7 mm
Matrix/resolution	112×93/3.4×3.5 mm
Echo time	63 ms
Repetition time	1 respiratory cycle
Imaging time per respiration	1,648 ms
EPI/half-Fourier/SENSE factor	51/0.6/2
Diffusion gradients	3 orthogonal directions
Duration/distance	22.3 ms/32.0 ms
b-values	0 (2 NSA), 50 s/mm2 (2 NSA ×3 directions), and 800 s/mm2 (6 NSA ×3 directions)
Fat-suppression method	SPIR
Water–fat shift/BW	9.2 pixel/23.6 Hz
BW in EPI frequency direction	1,680.3 Hz
Acquisition time	~3 min (1:12 min without gating)

Abbreviations: AP, anterior–posterior; BW, bandwidth; EPI, echo-planar imaging; FOV, field of view; NSA, number of averages; RL, right–left; SENSE, parallel imaging with sensitivity encoding; SPIR, spectral presaturation by inversion recovery.

Image analysis

The IVIM concept assumes a 2-compartment model of extravascular and intravascular space leading to a biexponential approach of the signal attenuation. where D is the true diffusion coefficient, D* the pseudodiffusion coefficient, f the perfusion fraction, and S(b) and S(0) the signal intensities with and without motion-probing gradients, respectively. D and f were approximated as introduced by Le Bihan.6 Recently, it was applied to abdominal imaging with the 3 b-values: b0=0, b1=50, and b2=800 s/mm2.18

For comparison, the ADC that was conventionally obtained from the 2 b-values was also calculated:

Parameter maps were calculated offline by using MATLAB (MathWorks, Natick, MA, USA).

All HCCs were analyzed in consensus by a board-certified radiologist with >11 years of clinical experience in MRI and a physicist with >16 years of experience in DWI under consideration of all the available anatomical images. Tumor response was assessed only for HCCs >2 cm in diameter in order to prevent partial volume effects. An ROI covering the whole tumor was placed in a central slice without visible motion artifacts, pixel misalignments, or susceptibility artifacts. For each HCC, one ROI was placed on the diffusion-weighted image with b=800 s/mm2 as large as possible, excluding the areas close to the rim, with initially very low SNR or containing intralesional hemorrhage. The anatomical position of each ROI was visually cross-checked between the images with b=0, 50, and 800 s/mm2; then, the ROI was copied into the parameter maps, and mean values and standard deviations were determined.

For the clinical response assessment, changes in tumor size were determined by using the longest diameter of the whole tumor (LDT) on the first and the second (if available) follow up MRI measured on 2D slices. Moreover, the change in the longest diameter of regions with contrast enhancement (LDCE) at the time point of the first follow up was measured in the arterial phase images. In a case in which the HCC appeared hyperintense on the native images due to lipiodol deposition following TACE, subtraction images (contrast-enhanced minus native images) were used for measuring LDCE.

HCCs were categorized into two groups, Group A (responders) with definite LDT reduction (at least 5%) on the first follow up and stable or further reduced LDT on the second follow up (if available) and also reduced LDCE (at least 5%) and Group B (nonresponders) with stable (LDT reduction <5%) or increased LDT on the first follow up and stable or further increased LDT on the second follow up (if no other therapies were performed meanwhile) and also unchanged (LDCE reduction <5%) or increased LDCE.

Statistical analysis

Statistically significant differences (P<0.05) between groups were tested by using nonparametric Mann–Whitney U-test (independent samples) and with Student’s t-test (paired samples) after proving Gaussian distribution by using Kolmogorov–Smirnov test. Discrimination between groups was analyzed by using receiver operating characteristics. Statistics were performed by using Statistical Package for Social Sciences (Version 21.0; IBM, Armonk, NY, USA).

Results

A total of 31 HCCs were eligible for analysis, including 1 HCC per patient in 21 patients, 2 HCCs in 3 patients, and 4 HCCs in 1 patient. The pretherapeutically measured mean lesion size was 54±27 mm (range =19–143 mm). Seventeen HCCs were assigned to Group A with mean LDT change (postvalue minus prevalue) of −16% (range =−37% to −5%, n=17) at the first follow up and −29% (range =−60% to −10%, n=16) at the second follow up and a mean LDCE reduction of −58% (range =−100% to −14%, n=17) at the first follow up; 14 HCCs were assigned to Group B with mean LDT change of +5% (range =−4% to +26%, n=14) at the first follow up and +2% (range =−16% to +31%, n=11) at the second follow up (decreased LDT values in 4 cases can be explained by further treatments). The mean LDCE change in Group B was +2% (range =−2% to +19%, n=9); in 5 cases, the quality of the images with contrast enhancement was not sufficient for mRECIST analysis due to motion artifacts (n=3) or due to hypovascularity in all dynamic phases (n=2). The overall mean survival time was longer in Group A (mean =744±484 days, range =81–1,637 days) than in Group B (mean =654±643 days, range =129–1,861 days), but the difference was not statistically significant.

Figure 1 shows the examples of diffusion-weighted images and maps of ADC(0,800), D′, and f′. Several HCCs of Group A showed a homogeneous response: loss of hyperintensity on DWI b800 measured posttherapeutically and an increase in ADC(0,800) with an enhanced increase in D′ and a decrease in f′ in the entire tumor (Figure 1A). In some of these cases, an additional loss of hyperintensity on DWI b0 was found (Figure 1B). Other HCCs showed more heterogeneous response (Figure 1C). The values represent a net effect measured in the ROI. Some other HCCs showed homogeneous hyperintense appearance after therapy, but with a decrease in size or signal intensity. In one HCC, there was no increase in ADC(0,800) and D′, but f′ was markedly decreased. Table 2 shows the results of DWI analysis. For the pretherapeutically measured parameter values, no differences were found between Groups A and B. Posttherapeutically, for Group A, an increase in ADC(0,800) and D′ by ~30% (P=0.00004) and ~42% (P=0.00001), respectively, and a decrease in f′ by ~37% (P=0.00094) were measured, whereas for Group B no significant changes were detected. For Group A compared to Group B, posttherapeutic D′ values were significantly larger (P=0.029), f′ values were significantly lower (P=0.017), and the ADC(0,800) showed larger values only by tendency. The changes in the parameter values (differences between post- and pretreatment values) showed the following significant between-group differences: for Group A compared to Group B, a significant larger increase in ADC(0,800) and D′ (P=0.00190 and P=0.00016, respectively) and a significant larger decrease in f′ (P=0.04384) were found. Responders and nonresponders were better differentiated by changes in D′ than by changes in ADC(0,800) and f′ or by posttherapeutical D′ and f′ values as indicated by the results of receiver operating characteristic analysis (area under the curve =0.878 compared to 0.819, 0.714, 0.731, and 0.750, respectively).

Figure 1

Diffusion-weighted images with maps of ADC(0,800), D′, and f′ pre- and post-TACE (A–C, E) and RE (D).

Notes: Examples with different changes in LDT and LDCE can be seen: (A) HCC with a decrease of 24% in LDT and 100% in LDCE (Group A, responders). The homogeneous increase in ADC(0,800) might indicate necrosis over the entire tumor area. This is supported by the stronger increase in D′ accompanied by the decrease in f′. The degree of necrosis might be underestimated by ADC(0,800). (B) The decrease in LDT and LDCE was 10% and 100%, respectively, for both lesions (Group A). After therapy, both HCCs show an increase in ACD(0,800), and they appear isointense in the postimage but were slightly hypointense in the preimage. In the DWI images, both HCCs appeared hypointense after therapy, but also in the b0 image. This might be due to coagulative necrosis and is supported by the even stronger increase in D′ and the decrease in f′. (C) This HCC with a decrease of 5% in LDT and 67% in LDCE (Group A) shows heterogeneous response. The main part remains hyperintense in the b800 postimage with no changes in ADC(0,800), D′ and f′ indicating vital tumor. A small area (star) shows clear increase in ADC(0,800) and D′ and decrease in f′, which are in accordance with necrosis. Another small area (round circle), which appears hypointense in all the DWI images, shows no increase in ADC(0,800) and D′, but a decrease in f′. This might indicate coagulative necrosis or hemorrhage or embolization. (D) HCC with a decrease of 10% in LDT and 27% in LDCE (Group A). The intensity on contrast-enhanced images (not shown) is also decreased. On DWI, the area of hyperintensity is decreased, but the ADC(0,800) is slightly decreased. However, D′ and f′ show that the decrease in ADC(0,800) is caused by a decrease in f′, whereas D′ is rather unchanged. This might be associated with embolization or necrosis of low degree. (E) HCC with an increase of 11% in LDT and 19% in LDCE (Group B, nonresponders). On DWI, the area of hyperintensity is also slightly increased. ADC(0,800) and D′ values are rather unchanged, indicating vital tumor. The slight decrease of f′ may be caused by embolization.

Abbreviations: ADC, apparent diffusion coefficient; DWI, diffusion-weighted imaging; HCC, hepatocellular carcinoma; LDCE, longest diameter of the region with contrast enhancement; LDT, longest diameter of the whole tumor on morphological images; RE, radioembolization; TACE, transcatheter arterial chemoembolization.

Table 2

Apparent diffusion coefficient ADC(0,800) and parameters D′ and f′ estimated by using the IVIM model before and after therapy with changes

	LDT	LDCE	LDT	ADC(0,800)			D′			f′
	(%)			(10−6 mm2/s)			(10−6 mm2/s)
	rdiff	rdiff	rdiff2	pre	post	diff	pre	post	diff	pre	post	diff
Group A
Mean	−16	−58	−29	1,067	1,358	291	953	1,323	370	0.087	0.052	−0.035
SD	9	32	15	262	258	213	218	263	229	0.042	0.037	0.036
N	17	17	16	17	17	17	17	17	17	17	17	17
Min	−37	−100	−60	691	928	−61	592	864	1	0.032	0.002	−0.136
Max	−5	−14	−10	1,591	1,822	641	1,317	1,823	693	0.192	0.152	0.026
Group B
Mean	5	2	2	1,190	1,207	17	1,067	1,103	37	0.094	0.088	−0.006
SD	7	6	12	260	285	205	233	259	167	0.036	0.044	0.034
N	14	9	11	14	14	14	14	14	14	14	14	14
Min	−4	−2	−16	775	748	−249	731	685	−205	0.045	0.020	−0.033
Max	26	19	31	1,774	1,721	310	1,628	1,633	280	0.170	0.170	0.072
P-value				ns	ns	0.00190	ns	0.02902	0.00016	ns	0.01706	0.04384

Notes: P-values are related to differences between Groups A (responders) and B (nonresponders). Group A includes all HCCs with definite reduction of LDT and Group B all other cases. Pre- and postvalues of ADC(0,800), D′, and f′ differed only within Group A (P=0.00004, P=0.00001, and P=0.00094, respectively). Moreover, the LDT changes at first follow up differed from those at second follow up and the ADC(0,800) changes differed from the D′ changes only within Group A (P=0.00008 and P=0.00018, respectively).

Abbreviations: ADC, apparent diffusion coefficient; diff, absolute difference of postvalues (at first follow up) minus prevalues; HCC, hepatocellular carcinoma; IVIM, intravoxel incoherent motion; LDCE, longest diameter of the region with contrast enhancement; LDT, longest diameter of the whole tumor; ns, not significant; post, at the first follow up; pre, before therapy; rdiff, relative difference at the first follow up; rdiff2, relative difference at the second follow up; SD, standard deviation.

Discussion

In the present study, a simplified IVIM approach was used to obtain parameter maps reflecting intralesional diffusion and perfusion changes in patients with HCC undergoing locoregional tumor therapy. As a main result, responders showed an increase in the estimated diffusion coefficient D′ accompanied by a decrease in the estimated perfusion fraction f′. In contrast, no significant changes were found in nonresponders. Responders could be better differentiated from nonresponders by changes in D′ than by changes in the conventionally determined ADC(0,800). Moreover, for differentiation, the changes in D′ were better suited than the changes in f′ and the posttherapeutical D′ and f′ values. The pretherapeutically measured parameter values showed no between-group differences, and thus, they were nonpredictive of response.

Previous studies on locoregional therapy response evaluation by DWI in HCCs focused on the assessment of the conventional ADC and typically implemented only a low maximum b-value of 500 s/mm2 for the assessment of RE19,20 and TACE.21–26 Because low maximum b-values have the disadvantage of higher perfusion sensitivity of the corresponding ADC values, recently published studies used a higher maximum b-value of 1,000 s/mm2, which reduces (but not eliminates) the effects of perfusion on ADC.4 One of these studies reported a correlation of the ADC change with the histopathological degree of necrosis.4 The results of the present work suggest that changes in D′ might better reflect necrosis than changes in ADC, as ADC is influenced by counteracting perfusion changes. Stable or decreased D′ found in nonresponders may indicate persisting viable tumor.2,27

In the IVIM model, information about perfusion is provided by the perfusion fraction f and the pseudodiffusion coefficient D*.7,8 In nude mice implanted with the HT29 colorectal cancer cells, f correlated with histological microvessel density.28 In the liver, a voxel-wise simultaneous determination of D, f, and D* by using nonlinear biexponential fitting is challenging due to motion-induced signal variations and low SNR.7,14 For data with low SNR, fitted f values were known to be systematically overestimated.29,30 On the other hand, the acquisition of numerous b-values with sufficient signal averages leads to long acquisition times, especially at 1.5 T.9,31 Thus, D* was kept constant,32 ROI-wise fitting instead of voxel-based analysis was performed,9–12 and/or advanced fitting procedures such as Bayesian fitting33 or segmented fitting3,9–13 were used. In patients with diffuse liver diseases such as isolated steatosis, steatohepatitis, fibrosis, and cirrhosis, D and especially D* were found to be pathologically altered, whereby for accurate D* determination the acquisition of b-values between 0 and 50 s/mm2 was essential.9–13 However, in regions with low f values (small IVIM effects), problems with unreliable fit results occurred.32 In malignant liver lesions, IVIM analysis was further complicated by the fact that D* and f values are lower than those in healthy or surrounding liver tissues,15–17 especially in hypovascular lesions34 and in areas with necrosis.35 Slow or stagnant blood flow through small damaged tumor vessels may lead to low D* values (and indirectly also to lower f values). This hypothesis is in line with a mouse mammary carcinoma model in which D* negatively correlated with the interstitial fluid pressure (IFP) that influences blood flow.36 In case of low D* and f values, the use of nonlinear biexponential fitting procedures is hampered by signal decay curves with only weak biexponential behavior and small IVIM effects, leading to fitting failure.37 For HCCs and liver metastases, unacceptably poor reproducibility especially of fitted D* was found.15,17 The determination of low D* values is challenging, and it cannot be improved by the acquisition of b-values between 0 and 50 s/mm2 due to the flatter curve at low b-values. Thus, dedicated analysis methods are required for the clinical investigation of malignant liver lesions. The determination of estimated diffusion coefficients D′ and perfusion fractions f′ from a low number of b-values with sufficient signal averages by using numerically stable voxel-wise analysis methods such as approximation formula6 or asymptotic fitting30 is promising.18,27 Although f′ (and also D′) values depend on the b-values used and are influenced not only by f but also by D* (not directly determined), f′ may serve as an additional empirical parameter that enables the assessment of perfusion effects. Its voxel-wise determination is an important advantage for treatment response evaluation, enabling a spatially resolved assessment.

In the present study, the same simplified IVIM approach and b-values were used as by Penner et al18. The pretherapeutically measured f′ values of Group A (0.087±0.042) and Group B (0.094±0.036) in the present study are in agreement with the previous work (0.089±0.045); f values found in studies using other IVIM approaches were higher.15,16,27 However, it was expected that D′ values overestimate D and f′ values underestimate f, especially for tissues with low D*.18 In general, results from different studies may vary according to not only the b-values used but also the localization of the HCC within the liver,13,18 gating techniques used,14 contrast agent injection,16 and echo and repetition times.38 For therapy monitoring, it is important to use the same MRI scanner hardware and DWI protocol in the follow up MRI as performed in the present study.

Current knowledge about IVIM analysis for response assessment in the liver is limited. The results of increased D′ and decreased f′ in responders are in agreement with increased D and decreased f and D* in histologically confirmed coagulation necrosis area of rabbit VX2 tumors after radiofrequency ablation.39 In a recent study, the early response of liver metastases from colorectal cancer to a combined therapy of antiangiogenic plus cytotoxic agents (FOLFIRI plus bevacizumab) was assessed.40 Responders showed unchanged D and D* values, but decreased f values as a sign of reduced microvessel density. This was explained by the short time interval from drug administration (2 weeks), where only perfusion but no apoptosis took place. Tumor vasculature normalized before its destruction. In a rabbit VX2 liver tumor model, perfusion changes caused by a vascular disruption agent were observed by D* and f.41 On the other hand, in a response analysis of vital tumor regions of neuroendocrine liver metastases to RE by fractioned ADC analysis, responders showed lower pretherapeutical ADC(0, 50) values and a larger therapy-induced increase in this perfusion-sensitive parameter.42 Lower pretherapeutical values might be related to higher pretherapeutical IFPs, the increase after therapy may indicate a decrease in IFP due to embolization, which may improve drug penetration.43 In a study addressing the response of advanced HCCs following antiangiogenic therapy with sorafenib, D remained stable, indicating absence of significant necrotic effects, whereas f differentiated responders from nonresponders by lower pretherapeutic values and increased values after therapy.27 This finding may be explained by higher pretherapeutical IFP in responders that decreased after therapy due to less microvessel leakage. It has to be mentioned that observed changes in f may be partly caused by changes in relaxation times.38 This means that the measured reduction in f′ in this work, which is probably related to vessel destruction due to necrosis or embolization, might also be in part an effect of an increase in tissue relaxation time T2 due to necrosis.

In this basic retrospective study, changes in functional IVIM parameters were investigated 6 weeks postinterventionally on the first follow up MRI. However, DWI has the potential to detect changes earlier after therapy when changes in tumor size are still absent, as demonstrated for HCCs after TACE.44 Thus, in future studies, the optimal time point for early response assessment by D′ and f′ measurements should be determined. Moreover, the reproducibility of the parameters determined with the simplified IVIM approach should be analyzed. In future studies, therapy response analyses should also be investigated separately for different types of interventions in order to provide therapy-specific response criteria.

Conclusion

The present work provided spatially resolved information on diffusion and perfusion changes of HCCs after locoregional therapy from 3 b-values DWI by using a simplified IVIM approach. D′ was better suited than conventional ADC(0,800) to assess therapeutic diffusion changes due to counteracting perfusion influences, which were separately addressed by f′. The results indicate that an IVIM-based simplified analysis has the potential to improve the accuracy of HCC response assessment compared to conventional ADC analysis.