Figures
Abstract
Objective
To investigate the intravoxel incoherent motion diffusion weighted imaging (IVIM-DWI) as a potential valuable marker to monitor the therapy responses of VX2 to radiofrequency ablation (RF Ablation).
Methods
The institutional animal care and use committee approved this study. In 10 VX2 tumor–bearing rabbits, IVIM-DWI examinations were performed with a 3.0T imaging unit by using 16 b values from 0 to 800 sec/mm2. The true diffusion coefficient (D), pseudodiffusion coefficient (D*) and perfusion fraction (f) of tumors were compared between before and instantly after RF Ablation treatment. The differences of D, D* and f and conventional perfusion parameters (from perfusion CT and dynamic enhanced magnetic resonance imaging, DCE-MRI) in the coagulation necrosis area, residual unablated area, untreated area, and normal control had been calculated by compared t- test. The correlation between f or D* with perfusion weighted CT including blood flow, BF (milliliter per 100 mL/min), blood volume, BV (milliliter per 100 mL/min), and capillary permeability–surface area, PMB (as a fraction) or from DCE-MRI: transfer constant (Ktrans), extra-vascular extra-cellular volume fraction (Ve) and reflux constant (Kep) values had been analyzed by region-of-interest (ROI) methods to calculate Pearson’s correlation coefficients.
Results
In the ablated necrosis areas, f and D* significantly decreased and D significantly increased, compared with residual unblazed areas or untreated control groups and normal control groups (P < 0.001). The IVIM-DWI derived f parameters showed significant increases in the residual unablated tumor area. There was no significant correlations between f or D* and conventional perfusion parameters.
Citation: Guo Z, Zhang Q, Li X, Jing Z (2015) Intravoxel Incoherent Motion Diffusion Weighted MR Imaging for Monitoring the Instantly Therapeutic Efficacy of Radiofrequency Ablation in Rabbit VX2 Tumors without Evident Links between Conventional Perfusion Weighted Images. PLoS ONE 10(5): e0127964. https://doi.org/10.1371/journal.pone.0127964
Academic Editor: Johan Pallud, Sainte-Anne Hospital Center, FRANCE
Received: November 21, 2014; Accepted: April 21, 2015; Published: May 28, 2015
Copyright: © 2015 Guo et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: All relevant data are within the paper.
Funding: This work was supported by the Haikou Municipal Bureau of Industry and Information Bureau of the key projects (contract/grant numbers 2012075) and Hainan Natural Science Foundation Project (contract/grant numbers 813261).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Image-guided percutaneous ablative therapies using thermal energy sources, such as radiofrequency (RF) [1], microwave [2], laser-induced thermal therapy [3], and high-intensity focused ultrasound [4], are rapidly evolving as nonsurgical minimally invasive techniques for the treatment of malignant neoplasms. Imaging of solid tumors after RF Ablation is crucial to judge the completeness of the ablation and later, to detect therapy response [5]. This assessment is usually performed with contrast-enhanced computed tomography [6] or magnetic resonance imaging [7]. When the use of contrast agents is contraindicated, assessment of the residual tumor can be a problem. In addition, assessment of the accuracy of ablation completeness by immediate CT or MR imaging has inevitable limitations, because postablation changes make an accurate diagnosis difficult [8][9]. The development of imaging tests without contrast agent is a challenge to immediately assess the treatment effects of RF Ablation.
Diffusion-weighted imaging (DWI) has been applied to measure tumor changes to assess the therapy response [10][11]. Normally, the apparent diffusion coefficient (ADC) of water molecules within living tissues is derived from diffusion images with an assumption that the water molecular diffusion is a random process [12]. This type of diffusion is also known as Gaussian diffusion and is characterized by a simple mono-exponential decay model. However, tumor tissue is a complex and restricted environment, which hinders the distribution of water molecules, resulting in distributions that are far from Gaussian diffusion [13]. DWI reflects the microscopic translational motions from the influence of both capillary perfusion and molecular diffusion that occurs within each image voxel. Increasing diffusion gradients lead to a multi-exponential signal decrease within each voxel, in agreement the intravoxel incoherent motion (IVIM) theory [14][15]. IVIM-DWI provides methods to measure the respective influence of true molecular diffusion and microcirculation within a voxel. The parameters derived from IVIM including: diffusivity (D), perfusivity (D*) and its volume fraction (f) values, which reflect the extent of diffusion and blood perfusion at the capillary level.
However, whether the perfusion parameters D* and f correlate with standard perfusion CT and MRI measurement has not been fully illuminated. The purpose of this study was to investigate the value of IVIM-DWI, DCE-MRI, and perfusion CT techniques in assessing the response of VX2 to RF Ablation therapy, with particular emphasis on the role of IVIM-DWI as a potential valuable imaging marker of therapy response.
Material and Methods
VX2 tumor model
Ten New Zealand white rabbits weighing from 2.5 to 3.0 kg were used. The protocol was approved by the Committee on the Ethics of Animal Experiments of Beijing Union Medical College Hospital (Permit Number: SYXK20100028), and the study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Research Council. All surgery was performed under anesthesia using sodium pentobarbital administered into the marginal ear vein. The VX2 tumor strain (National Infrastructure of Cell Line Resources, Beijing, China) was transplanted into the hind thigh muscles of tumor-carrier rabbits. The tumor-carrier rabbit was then sacrificed by an overdose of anesthesia, and the VX2 tumor was excised and placed in saline. Necrotic tissue was discarded. The tumor tissue was cut into 1-mm3 fragments. Two to three tumor fragments were symmetrically implanted into each side of the erector spinae of 10 rabbits under CT-guidance. After 10–14 d, both plain and contrast enhancement CT scans were performed on the 10 rabbits to evaluate the growth of VX2 tumor. The implanted tumor mass had grown to an average size of 2.0 cm.
CT guided RF Ablation for the treatment of VX2 tumors
To prepare for RF Ablation, the hair was shaved and the surgical site was disinfected using povidone-iodine. A scout scan through the tumor was obtained. Axial CT scans both plain and perfusion CT imaging had been used to select the target treatment level of the tumor. VX2 tumor ablations were achieved using a single 15-gauge internally cooled bipolar electrodes (Celon ProSurge, Berlin, Germany) with an active tip of 2.0 cm length. The RF Ablation needle (CelonProSurge, Berlin, Germany) was injected at a depth two-thirds of the target tumor’s diameter under the gauidence of CT. The applied parameters as follow: power 25W, about 5 min, total 5 KJ. The electrode was internally cooled with chilled saline circulated by a peristaltic pump. The RF system is constituted by one 200-watt power generator supplying pulsed energy automatically adjusted based upon tissue impedance during a full 5-minute ablation session. Each animal was treated once in a single setting. The diameter of RF Ablation zone was about 1cm. At least 60 min apart, allowing the contrast material from the preceding acquisition to clear, perfusion CT imaging and MRI examination was performed again to evaluate the curative effects of the RF Ablation as shown in Fig 1. The rabbits were fastened to a plastic board to ensure the same postion during the whole experiment.
Perfusion CT Examination and Image Analysis
For lesion detection, an unenhanced scan of the entire tumor was obtained before beginning cine CT scanning. The native scan parameters were as follows: 80 kVp, 120 mAs, slice thickness 5 mm, matrix = 512×512, FOV = 150 mm, 30 slices. This was followed by a cine scan (1 phase / second) during the whole proccess of intravenous administration of contrast medium, iopromide (350 mg of iodine per milliliter, 2 mL/kg) at a rate of 2mL/sec into an ear vein via an 18-gauge catheter. Contrast medium was injected into an ear vein by using a power injector (Stellant; MEDRAD, Pittsburgh, PA). All CT examinations were performed with a multisection spiral CT scanner (Somatom Perspective; Siemens Medical Solutions, Forchheim, Germany). Cine CT scanning parameters were as follows: 120 mAs/rotation, 80 kV, 512 × 512 matrix, 150 mm FOV, section thickness/spacing of 5/5 mm, and an acquisition series of 90 time frames at 1-second intervals with 10 phases of blank scans before the injection of contrast agent. Cine CT was performed again after the CT guided RF Ablation therapy.
All CT perfusion images were transferred to an image processing workstation (Syngo via; Siemens Healthcare, Erlangen, Germany). Tumor perfusion software (Siemens Healthcare Systems) was used for CT perfusion analysis. The software used curve-fitting by least mean squares to obtain mathematical descriptions of the time–density curves for each pixel [16]. Regions of interest were placed on tumors (RF Ablation treated necrosis areas, residual unablated areas and contralateral untreated, namely TRF Ablation, Tresidual, and Tcontrol), and normal tissue (Normal), as shown in Fig 2B to 2E. Subsequently, functional maps had been automatically generated, including blood flow BF (milliliter per 100 mL/min), blood volume BV (milliliter per 100 mL/min), and capillary permeability–surface area PMB (as a fraction). Each pixel value had been represented with a color scale. The BF, BV, and PMB were measured in the tumors and in an area with normal muscle tissue (no tumor area) in each case by 2 senior radiologists (G.Z.Y. and Z.Q.).
(A) Arterial input function for the simulation curve; (B) and (C). Raw images for CT perfusion mapping after and before RF Ablation treatment. (D to G) before RF Ablation treatment (H to K) after RF Ablation treatment: (D) fused and (E) pure BF maps (F) Blood volume (unit in mL/100 mg), (G) permeability of capillary vessel surface (scale in 0.5 mL/min/100 mg) show greater enhancement values in the tumor than that of the adjacent muscle. Higher blood flow, blood volume and permeability had been observed for the tumor (arrow) compared to the adjacent muscle tissues and almost equal to contralateral tumor. (I to K) Significantly reduced values of BF, BV, and PMB in the RF Ablation treated tumor (arrow) than that of the contralateral untreated tumors.
The absence of early arterial hyperperfusion at the margins of the tumor was interpreted as complete remission and grouped as TRF Ablation. Peripheral hyperperfusion of the lesion in the early arterial hyperperfusion phase was considered to indicate residual vital unablated tumor tissue, grouped as Tresidual. After RF Ablation has been completed, rabbits were moved back to the MR scanner for post-treatment MRI.
MRI Examination
All MR imaging examinations were performed at our institution with the same 3.0-T MR imaging system Skyra (Siemens Healthcare, Erlangen, Germany) and two 8-channel surface radiofrequency receiving coils. The MR imaging protocol included transverse T1-weighted turbo spin echo and transverse T2-weighted turbo spin echo sequences. DW images were acquired prior to gadolinium chelate injection. Finally, a time-resolved MR angiography with Stochastic Trajectories (TWIST) T1-weighted perfusion weighted sequence was repeated 75 times during the whole process of contrast injection with 5 phases of blank scans before the injection of contrast agent. The total of 0.2 mL per kilogram of body weight gadoteric acid (Omniscan TM; General Electric healthcare, USA) had been injected both before and after RF Ablation treatment. The time interval of two injections was over 1 hour to eliminate the residual effects of contrast agent. Protocol details are shown in Table 1.
Image Analysis of DCE-MRI
All MRI perfusion images were transferred to an image processing workstation, tissue 4D software (Siemens Healthcare Systems) was used for time resolved DCE-MRI analysis. A parameter related to vessel permeability and tissue blood flow transfer constant (Ktrans), extra-vascular extra-cellular volume fraction (Ve) and reflux constant (Kep) values had been calcualted. Calculation is based on the Tofts model and pharmacokinetic calculation using a 2-compartment model on a pixel-by-pixel basis [17]. Regions of interest (ROIs) were manually drawn as the same way as CT perfusion analysis. The Ktrans, Kep, Ve had been measured in the tumors (TRF Ablation, Tresidual, Tcontrol), and normal muscle tissue (Normal) in each case. Tumor lesions were identified by two senior radiologists (GZY and ZQ) for consensus.
Image Analysis of DWI
The bi-exponential model was mathematically expressed as follows to analysis IVIM-DWI:
1
where f is the unit of perfusion volume fraction. D* is the pseudodiffusion coefficient influenced by tissue microcapillary perfusion, and D is the true diffusion coefficient. The diffusion quantification and analysis of IVIM were conducted by the MITK-Diffusion software (http://www.mitk.org/) [18].
DWI datasets were analyzed in consensus by two senior radiology researchers. Each lesion observed in the DW images had also been defined on the T2W, T1W with post-contrast images, all serving as a roadmap for regions of interest (ROIs) demarcation. ROIs were drawn nearly the same methods as perfusion CT. D, D*, f, were calculated in each of the ROIs, as shown in Fig 3B to 3E). All the differences before and after ablation were compared to confirm any recorded changes. The typical image processing is demonstrated in Fig 4A to 4C).
A. Arterial input function for the simulation curve; and (B-E) ROIs selection for TRF Ablation, Tresidual, Tcontrol and normal regions for claculation the Ktrans, Ve, and Kep parameters before RF Ablation treatment; F-J after RF Ablation treatment corresponding to A to E. All the mapping from A to J show reduced perfusion values in the RF Ablation treated tumor than that of the contralateral untreated tumors and before treatment in color coded mappings. K to L show the T2, T1 post contrast enhancement MRI before RF Ablation treatment and after treatment (M to N). All the mapping from K to N show instantly morphological chages of routine MRI in the RF Ablation treated tumor. The T2WI and T1WI showed the necrosic focus in tumor after RF Ablation (M, N).
(A) Top row demonstrates the sequencial diffusion weighted images b value from 50 to 800, low row corresponding to the before treatment as control. (B to C) the IVIM-DWI parameters f, D and D* derived from above multiple diffusion weighted images (A). The IVIM analysis derived f, D, D* map were obtained with pixel-by-pixel basis from the each pixel's signal intensity in the IVIM-DWI data. It was shown that reduced f value and increased D value indicated the RF Ablation therapy necrosis responses.
Trypan blue injection
Trypan blue (TB; Sigma-Aldrich Co., LLC, St. Louis, MO, USA) was dissolved in 10 mM phosphate buffer (pH 7.4) and sterilized by passing through 0.2 μm pore membrane filters. The final dye was 2.0% TB. Animals were injected in the ear vein with ~1.0 mL/kg body wt and then sacrificed by anesthetic overdose. The animals were visually inspected for dye uptake into the tumor, indicated by blue.
Histopathologic analysis
The tumor and surrounding tissue were orderly harvested according to the distance from the central necrosis area. The necrosis area was confirmed by in vivo TB injection intravenously into the ear margin. All histopathologic analyses were used to gather evidence of the curative effects of thermal ablations. The specimen tissues were cut into 5.0-μm sections and prepared for histological assessment using hematoxylin and eosin staining. The histopathologic analysis was performed using microscopy by one senior pathologist.
Statistical analyses
Comparisons between before and after treatment were analyzed using the compared t-tests for parametric continuous data of BF, BV, PMB, Ktrans, Kep, Ve, D, D*, and f in the RF Ablation treated necrosis, residual unablated, untreatment control, and normal groups. Pearson's correlation analysis were used to evaluate the correlations between the f or D* values obtained with IVIM and BF, BV, PMB, Ktrans, Kep, Ve values obtained with perfusion CT and MRI. All these data had been analyzed by the ROI-based approach. Statistical analyses were performed using SPSS for Windows Version 16.0 (SPSS Inc., Chicago, IL, USA), and P < 0.05 was considered significant.
Results
Establishment of rabbit VX2 tumor model
Direct mass implantation tumor had 100% successful established tumors and produced two similar tumors at same level for the following tests. This direct mass implantation was suitable for evaluating RF Ablation efficacy, the implanted tumor mass grew to an average size of 2.0 cm (ranging from 1.2 to 3.0 cm) by 10 to 12 days after its introduction.
Therapeutic effect of RF Ablation revealed by perfusion CT and DCE-MRI
Using perfusion CT and MRI technique, we evaluated a total of 10 VX2 tumor models before and immediately after treatment with RF Ablation. These lesions were located in both sides of erector spinae. RF Ablation treatment of these VX2 tumors produced an average of a 1.0 cm coagulation necrosis as shown by contrast enhancement CT and MRI. The coagulation necrosis either remained unchanged or had no contrast agent perfusion on the perfusion images, as shown in Figs 2 and 3. If the perfusion image showed a weighted portion, it meant that the tumor tissue was not completely ablated and is an indication of residual tumors. The descriptive analysis of the perfusion parameters measured for all the VX2 tumor models showed the following results (Table 2 and 3). The conventional perfusion parameters including BF, BV, PMB, Ktrans, Ve, and Kep were all significantly reduced (P < 0.001) in the TRF Ablation group, which indicated the therapy effects of RF Ablation on tumor. And the Ve, in the T residual significantly increased (P < 0.001) after RF Ablation treatment indicated the alterations of vascular exchage fuction alterations instantly after RF Ablation therapy.
Therapeutic effect of RF Ablation revealed by IVIM-Diffusion weighted imaging
A summary of three IVIM parameters for each group there were measured before and immediately after treatment with RF Ablation is presented in Table 3. The D values had significantly increased, f and D* value measured in the tumor were significantly reduced in the TRF Ablation group induced by RF Ablation therapy (P<0.001). The f value in the Tresidual was significantly increased in the Tresidual group after RF Ablation therapy (P<0.001) (Table 3). There were not significant changes in the other groups including Tcontrol and Normal before and after RF Ablation treatment.
IVIM-Diffusion and perfusion CT and MRI Correlation
Correlation between perfusion parameters derived by perfusion CT and DCE-MR measurements and IVIM-based parameters are summarized in Table 4. No significantly correlations of f or D* to conventional perfusion parameters had been found.
Histopathologic results
The rabbits in the treatment group were sacrificed immediately after RF Ablation treatment for histologic studies. For hemotoxylin-eosin–stained tumor specimens, ablation zones were observed to encompass the transitions from necrosis—residual unablated- normal tumor within samples collected immediately after RF Ablation. The ablated zones within tumor were pale relative to adjacent residual unablated tissue regions, indicating necrosis, away from the necrosis zone (Fig 5A). Partial loss of tissue was found at the center of the electrode insertion for the rabbit showing complete tumor ablation, and this was surrounded by a pale coagulated necrotic area and a darker portion with venous TB stain. The VX2 cells appeared to have undergone coagulation necrosis in a large area, and nucleus condensation, karyorrhexis, and karyolysis could be seen in individual cells. Coagulation necrosis is characterized by preservation of the basic cell outlines of necrotic cells (Fig 5B to 5E).
(A) Viable cells exclude Trypan blue (TB), while dead cells stain blue from TB uptake. TB staining of VX2 tumor after in vivo treatment shows damage to the tissue structure within the area occupied by the RF Ablation. VX2 tumor treated with RF Ablation in vivo, its curative effects had been verified by pathology Hematoxylin/eosin staining. (B) Region near the right side demonstrates areas near the RF Ablation needle and shows coagulation necrosis areas, X 100. (C) Region demonstrates nucleus of tumor cells was partially dissolved and necrotic, X 200. (D to E) Regions demonstrate transition from total coagulation necrosis to residual unablated tumor (partial necrosis), X 200.
Discussion
The major goals of this work were to evaluate the usefulness of IVIM-diffusion weighted, imaging-based methods to evaluate the therapeutic response of RF Ablation on VX2 tumors model. In our study, direct tumor mass implantation forms a spherical tumor, which is similar to the naturally occurring tumor. For RF Ablation treatment in this experiment, we set the RF energy at total 5KJ (25 W, about 5 minutes) and to ensure high consistency and repeatability of the therapeutic range. An average 1.0 cm focus is formed using these parameters. CT guidance is needed during RF Ablation treatment to verify tumor location and avoid tissue damage. Because the electrode had been inserted from the dorsal lateral side to the target tumor, the postion remain stable and without motion artifacts make the study of perfusion CT or MRI, IVIM-DWI and their relationship simple and repeatable.
We demonstrated that IVIM-DWI is able to depict therapy responses in a rabbit VX2 tumor model after administration of RF Ablation treatments, as shown in Fig 4. Among the IVIM-DW imaging parameters, relative changes in D*, D, and f can be utilized as predictors of tumor response immediately after RF Ablation. In the RF Ablation treated group, perfusion-related IVIM parameters, including D* and f immediately decreased, whereas diffusion related parameter D increased at the same time. Our findings in this study are consistent with previous study about the early evaluation of the treatment response after RF Ablation of unresectable lung tumours by DWI [19]. These results suggest that pre- and post-RF Ablation IVIM-DW can be used to predict the therapeutic outcome before the change in tumor size becomes detectable on CT or MRI images. In addition, perfusion-related IVIM parameters immediately increased in the residual tumor tissue, which is consistent with the changes of Ve and BV from conventional perfusion-weighted imaging.
One retrospective cohort study of MRI data included both pre- and post- laser interstitial thermal therapy (LITT) T2WI and DWI acquisitions, which yielded an Apparent Diffusion Co-efficient (ADC) map. Their initial results indicate that the ADC maps demonstrated highly sensitive, albeit non-specific, response to LITT-related changes in vivo, due to the significant change in water permeability in tissue as a result of treatment [20]. One in vivo canine study indicated that ADC may be a better tool to predict microwave-induced tissue necrosis than contrast-enhanced imaging, including DCE-MRI, and contrast enhanced ultrasound (CEUS) [21]. There is a substantial change in perfusion after MR-guided high-intensity focused ultrasound ablation of uterine fibroids due to vessel destruction. Low b-values ADC might be the good choice to emphasis perfusion changes in ADC maps [22]. Most recently, intravoxel incoherent motion (IVIM) diffusion-weighted (DW), magnetic resonance (MR) imaging has been used to evaluate the diagnostic value of quantitative assessment of the therapeutic efficacy of a vascular disrupting agent (VDA) (CKD-516) in rabbit VX2 liver tumors [23]. In this study, the authors derived the f and D*, in addition to D, which reflects the perfusion-insensitive water diffusion in tissues by performing biexponential IVIM model fitting to the multiple–b value DW imaging data [24]. Our current results and the conclusions from previous studies indicated that IVIM-DW could be used to obtained tumor perfusion information on the capillary level relevant to RF Ablation treatments and may be useful to evalue the reaction of tumor after RF Ablation treatments.
Another purpose of this study was to explore whether the IVIM-based D* and f metrics measured tissue perfusion by perfusion CT or DCE-MRI, and thus may deliver diffusion and perfusion information simultaneously without the need for intravenous, contrast media based on this VX2 tumor model. Many authors have indicated that IVIM-DW may be an alternative (or even superior) to DCE-MRI and could either link or correlate D* and f values with relative cerebral blood volume in brain [25–27] and modeling of the tumor vasculature in prostate tracer [28]. For this reason, we had examined VX2 tumor models by means of IVIM, DCE-MRI, and perfusion CT imaging, both before and after RF Ablation treatment, in the same experiment. Our current results indicate that the correlation coefficient between perfusion parameters derived from perfusion CT or DCE-MRI and IVIM-based perfusion parameters are low (r never above 0.4), irrelevant (r less than 0.2), or negative in most instances, as shown in Table 4. Our findings are in agreement with other authors’ conclusion regarding no clearly evident link between D*, f, and DSC- and DCE-derived metrics [29][30].
A cause of these findings may have been the low signal-to-noise ratio of D* estimates in VX2 tumor, the extensively altered hemodynamics and challenge to model capillary properties in tumor tissue, and the more complex than assumed nature of origin for the biexponential signal decay in the IVIM-DW measurements; this does not share the same principles as the intravascular tracer behavior [31]. In addition, it is worth noting that the kinetic model used in DCE-MRI or perfusion CT is only valid for tissues that contain small blood volumes [32], so the parameters might not be accurate for VX2 tumors because VX2 tumors have a larger blood volume?. It was believed that the IVIM method to measure perfusion remains technically challenging; for example, the optimal fitting algorithm and the optimal distribution of b values, as well as the accuracy and robustness of the measurements [33]. Overall, the experimental verification of the relation between the IVIM perfusion parameters and the standard perfusion parameters are still lacking.
This study had several limitations. First, we only evaluated the immediately changes of RF Ablation treatment. Thus, we could not assess late histologic changes or tumor recurrence that might have explained the late changes in IVIM parameters and its relation to conventional perfusion parameters. Second, to measure IVIM parameters, we selected one ROI on one axial MR image of the tumor. However, considering the three-dimensional structure of the tumors, measurement on cross-sectional images might be less representative of the entire tumor than volumetric measurements. Third, due to the lower spatial resolution and lower signal noise ratio in IVIM-DWI and its parameter mappings, only ROIs based analysis had been performed. At last, only the intramuscular tumors were evaluated in this study. Further research is needed to determine the depth and location of tumors (subcutaneous, intramuscular or intraperitoneal), size of tumors, and size of the ablative lesions that may have any influence on IVIM-DWI monitoring.
Conclusion
The therapeutic effect induced by RF Ablation treatment can be effectively evaluated with IVIM-derived parameters. The IVIM-DW derived perfusion parameters have the potential to indicate therapy effects immediately after RF Ablation treatment without significant association with classical tumor perfusion metrics derived from DCE-MRI and perfusion CT measurements.
Acknowledgments
This work was supported by the Haikou Municipal Bureau of Industry and Information Bureau of the key projects (contract/grant numbers 2012075) and Hainan Natural Science Foundation Project (contract/grant numbers 813261).
Author Contributions
Conceived and designed the experiments: ZYG QZ XGL. Analyzed the data: ZYG XGL ZYJ QZ. Wrote the paper: ZYG XGL ZYJ QZ.
References
- 1. Wistuba II, Gelovani JG, Jacoby JJ, Davis SE, Herbst RS. Methodological and practical challenges for personalized cancer therapies. Nat Rev Clin Oncol. 2011; 8:135–41. pmid:21364686
- 2. Nishikawa H, Kimura T, Kita R, Osaki Y. Radiofrequency ablation for hepatocellular carcinoma. Int J Hyperthermia. 2013; 29:558–68. pmid:23937321
- 3. Kitchin D, Lubner M, Ziemlewicz T, Hinshaw JL, Alexander M, Brace CL, et al. Microwave ablation of malignant hepatic tumours: intraperitoneal fluid instillation prevents collateral damage and allows more aggressive case selection. Int J Hyperthermia. 2014; 30:299–305. pmid:25144819
- 4. Vogl TJ, Straub R, Zangos S, Mack MG, Eichler K. MR-guided laser-induced thermotherapy (LITT) of liver tumours: experimental and clinical data. Int J Hyperthermia. 2004; 20:713–24. pmid:15675667
- 5. Hennedige T, Venkatesh SK. Imaging of hepatocellular carcinoma: diagnosis, staging and treatment monitoring. Cancer Imaging. 2013; 12:530–47. pmid:23400006
- 6. Kim YS, Rhim H, Lim HK, Park CK, Lee WJ, Do YS, et al. Completeness of treatment in hepatocellular carcinomas treated with image-guided tumor therapies: Evaluation of positive predictive value of contrast-enhanced CT with histopathologic correlation in the explanted liver specimen. J Comput Assist Tomogr. 2006; 30:578–82. pmid:16845287
- 7. Park MS, Kim S, Patel J, Hajdu CH, Do RK, Mannelli L, et al. Hepatocellular carcinoma: detection with diffusion-weighted versus contrast-enhanced magnetic resonance imaging in pretransplant patients. Hepatology. 2012; 56:140–8. pmid:22370974
- 8. Miki K, Shimomura T, Yamada H, Kishimoto K, Ohishi Y, Harada J, et al. Percutaneous cryoablation of renal cell carcinoma guided by horizontal open magnetic resonance imaging. Int J Urol. 2006; 13:880–4. pmid:16882047
- 9. Schima W, Ba-Ssalamah A, Kurtaran A, Schindl M, Gruenberger T. Post-treatment imaging of liver tumours. Cancer Imaging. 2007; 7:S28–36. pmid:17921098
- 10. Bains LJ, Zweifel M, Thoeny HC. Therapy response with diffusion MRI: an update. Cancer Imaging. 2012; 12:395–402. pmid:23022595
- 11. Hoang JK, Choudhury KR, Chang J, Craciunescu OI, Yoo DS, Brizel DM. Diffusion-weighted imaging for head and neck squamous cell carcinoma: quantifying repeatability to understand early treatment-induced change. AJR Am J Roentgenol. 2014; 203:1104–8. pmid:25341151
- 12. Padhani AR, Liu G, Koh DM, Chenevert TL, Thoeny HC, Takahara T, et al. Diffusion-weighted magnetic resonance imaging as a cancer biomarker: consensus and recommendations. Neoplasia. 2009; 11:102–25. pmid:19186405
- 13. Koh DM, Scurr E, Collins D, Kanber B, Norman A, Leach MO, et al. Predicting response of colorectal hepatic metastasis: value of pretreatment apparent diffusion coefficients. AJR Am J Roentgenol. 2007; 188:1001–8. pmid:17377036
- 14. Le Bihan D, Breton E, Lallemand D, Grenier P, Cabanis E, Laval-Jeantet M. MR imaging of intravoxel incoherent motions: application to diffusion and perfusion in neurologic disorders. Radiology. 1986; 161:401–7. pmid:3763909
- 15. Le Bihan D, Breton E, Lallemand D, Aubin ML, Vignaud J, Laval-Jeantet M. Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. Radiology. 1988; 168:497–505. pmid:3393671
- 16. Sanelli PC, Lev MH, Eastwood JD, Gonzalez RG, Lee TY. The effect of varying user-selected input parameters on quantitative values in CT perfusion maps. Acad Radiol. 2004; 11:1085–92. pmid:15530801
- 17. Mazzetti S, Gliozzi AS, Bracco C, Russo F, Regge D, Stasi M. Comparison between PUN and Tofts models in the quantification of dynamic contrast-enhanced MR imaging. Phys Med Biol. 2012; 57:8443–53. pmid:23202297
- 18. Fritzsche KH, Neher PF, Reicht I, van Bruggen T, Goch C, Reisert M, et al. MITK diffusion imaging. Methods Inf Med. 2012;51:441–8. pmid:23038239
- 19. Okuma T, Matsuoka T, Yamamoto A, Hamamoto S, Nakamura K, Inoue Y. Assessment of early treatment response after CT-guided radiofrequency ablation of unresectable lung tumours by diffusion-weighted MRI: a pilot study. Br J Radiol. 2009; 82:989–94. pmid:19470575
- 20. Viswanath S, Toth R, Rusu M, Sperling D, Lepor H, Futterer J, et al. Quantitative Evaluation of Treatment Related Changes on Multi-Parametric MRI after Laser Interstitial Thermal Therapy of Prostate Cancer. Proc Soc Photo Opt Instrum Eng. 2013;8671:86711F. pmid:24817802
- 21. Cheng HL, Haider MA, Dill-Macky MJ, Sweet JM, Trachtenberg J, Gertner MR. MRI and contrast-enhanced ultrasound monitoring of prostate microwave focal thermal therapy: an in vivo canine study. J Magn Reson Imaging. 2008; 28:136–43. pmid:18581405
- 22. Ikink ME, Voogt MJ, van den Bosch MA, Nijenhuis RJ, Keserci B, Kim YS, et al. Diffusion-weighted magnetic resonance imaging using different b-value combinations for the evaluation of treatment results after volumetric MR-guided high-intensity focused ultrasound ablation of uterine fibroids. Eur Radiol. 2014; 24:2118–27. pmid:24962829
- 23. Joo I, Lee JM, Han JK, Choi BI. Intravoxel incoherent motion diffusion-weighted MR imaging for monitoring the therapeutic efficacy of the vascular disrupting agent CKD-516 in rabbit VX2 liver tumors. Radiology. 2014; 272:417–26. pmid:24697148
- 24. Koh DM. Science to practice: can intravoxel incoherent motion diffusion-weighted MR imaging be used to assess tumor response to antivascular drugs? Radiology. 2014; 272:307–8. pmid:25058129
- 25. Wirestam R, Borg M, Brockstedt S, Lindgren A, Holtås S, Ståhlberg F. Perfusion-related parameters in intravoxel incoherent motion MR imaging compared with CBV and CBF measured by dynamic susceptibility-contrast MR technique. Acta Radiol. 2001; 42:123–8. pmid:11281143
- 26. Federau C, Meuli R, O'Brien K, Maeder P, Hagmann P. Perfusion measurement in brain gliomas with intravoxel incoherent motion MRI. AJNR Am J Neuroradiol. 2014; 35:256–62. pmid:23928134
- 27. Kim HS, Suh CH, Kim N, Choi CG, Kim SJ. Histogram analysis of intravoxel incoherent motion for differentiating recurrent tumor from treatment effect in patients with glioblastoma: initial clinical experience. AJNR Am J Neuroradiol. 2014; 35:490–7. pmid:23969343
- 28. Pang Y, Turkbey B, Bernardo M, Kruecker J, Kadoury S, Merino MJ, et al. Intravoxel incoherent motion MR imaging for prostate cancer: an evaluation of perfusion fraction and diffusion coefficient derived from different b-value combinations. Magn Reson Med. 2013; 69:553–62. pmid:22488794
- 29. Patel J, Sigmund EE, Rusinek H, Oei M, Babb JS, Taouli B. Diagnosis of cirrhosis with intravoxel incoherent motion diffusion MRI and dynamic contrast-enhanced MRI alone and in combination: preliminary experience. J Magn Reson Imaging. 2010; 31:589–600. pmid:20187201
- 30. Bisdas S, Braun C, Skardelly M, Schittenhelm J, Teo TH, Thng CH, et al. Correlative assessment of tumor microcirculation using contrast-enhanced perfusion MRI and intravoxel incoherent motion diffusion-weighted MRI: is there a link between them? NMR Biomed. 2014; 27:1184–91. pmid:25088433
- 31. Sourbron SP, Buckley DL. On the scope and interpretation of the Tofts models for DCE-MRI. Magn Reson Med. 2011; 66:735–45. pmid:21384424
- 32. Lemke A, Stieltjes B, Schad LR, Laun FB. Toward an optimal distribution of b values for intravoxel incoherent motion imaging. Magn Reson Imaging. 2011; 29:766–76. pmid:21549538
- 33. King MD, van Bruggen N, Busza AL, Houseman J, Williams SR, Gadian DG. Perfusion and diffusion MR imaging. Magn Reson Med. 1992; 24:288–301. pmid:1569868