Study design
This single-institutional study was conducted retrospectively and received institutional review board (IRB) approval from the Yale School of Medicine. Data collection and analysis were conducted in compliance with the Health Insurance Portability and Accountability Act (HIPPA). Study design was in agreement with the Standards for Reporting of Diagnostic Accuracy guidelines 31. Written informed consent from all subjects was waived by the IRB committee of Yale School of Medicine due to the retrospective nature of this study.
Study cohort
Initially, a total of 122 consecutive patients with biopsy-proven ICC lesions and history of systemic chemotherapy who underwent cTACE between October 2001 and February 2015 were evaluated. The following inclusion criteria were then applied: unresectable ICC (i.e. patients with mixed ICC and hepatocellular carcinoma were not included), naïve to locoregional treatment including percutaneous ablations or other intra-arterial therapies than cTACE as well as stereotactic body radiation therapy, and preprocedural multi-phase contrast-enhanced MRI with adequate/artifact-free image quality adequate for 3D quantitative tumor analysis.
Clinical and laboratory evaluation and staging
All patients had a complete clinical examination and baseline laboratory workup including bilirubin, albumin, and International Normalized Ratio (INR). Initial diagnosis of ICC was made on imaging and confirmed on pathology. Furthermore, the Child-Pugh classification for liver function and The Eastern Cooperative Oncology Group (ECOG) performance status were documented for each patient. The stage of disease was assessed using the Union for International Cancer Control (UICC) staging system.
cTACE protocol
All patients were discussed in the multidisciplinary tumor-board, and enrolled to undergo cTACE based on final consensus. The interventions were performed by the same interventional radiologist (with more than 20 years of experience) in a dedicated interventional radiology suite (Philips IR suites). After local anesthesia with lidocaine 1%, access was obtained through the common femoral artery via a 5 Fr vascular sheath, followed by a 0.035-in. guide wire in Seldinger technique. For orientation purposes, diagnostic angiography of the superior mesenteric artery and celiac trunk was obtained using a 5 Fr catheter to selectively advance into the tumor-supplying hepatic artery. Selective catheterization was achieved by placing a microcatheter and obtaining further imaging in order to more precisely target and spare healthy liver parenchyma. All patients were treated with a cTACE-protocol comprising of a 1:1 mixture of 50 mg doxorubicin (Adriamycin®, Pharmacia & Upjohn) and 10 mg of mitomycin-C with 10mL of iodized oil (Lipiodol®, Guerbet). This was followed by administration of 100–300 μm diameter microspheres (Embosphere®, Merit Medical); until complete stasis was reached, as the technical end point.
Imaging protocol
A standardized liver MRI protocol was performed in all patients enrolled using a 1.5T scanner (Magnetom Avanto, Siemens) with a phased array torso coil (repetition time ms/echo time ms, 5.77/2.77; field of view 320-400 mm; matrix 192 x 160; slice thickness 2.5 mm; receiver bandwidth 64 kHz; flip angle 10o). The protocol included single-shot breath-held gradient-echo diffusion weighted echo-planar images, axial T2-weighted fast spin-echo images, and unenhanced and contrast-enhanced (0.1 mmol/kg intravenous gadopentetate; Magnevist; Bayer) breath-held axial T1-weighted 3D fat-suppressed spoiled gradient-echo images in the hepatic arterial, portal venous and delayed phases (20s, 70s, and 180s, respectively)24,32,33.
Image analysis
All measurements were conducted using standard electronic calipers on Digital Imaging in Communications and Medicine (DICOM) files. Different sequences were assessed to distinguish between tumor enhancement and other hyperintense T1 signal abnormality (such as hemorrhage) in order to evaluate the true extent of tumor burden. If multiple lesions were present, the three largest lesions were assessed and the sum of total lesions in each patient was then processed in the analysis. For tumor burden analysis using each technique (1D, 2D and 3D), the most dominantly enhancing axial MRI sequence was used in each patient, since the enhancement pattern of ICC depends on tumor size and can vary accordingly 34. The largest overall tumor diameter (1D) and maximum cross-sectional area (2D), as well as the largest enhancing tumor diameter (1D) and maximum enhancing tumor area (2D) were measured by two radiological readers (XX with 5 years of experience and YY with 4 years of experience in abdominal MRI, respectively), using the RadiAnt™ DICOM Viewer (Medixant). The numbers used for the final analysis were the consensus of the simultaneous measurement and discussion of both readers; both readers were blinded to all clinical data.
For the 3D tumor assessment, a 3D quantitative semiautomatic tumor analysis software was used (IntelliSpacePortal V8, Philips ICAP) 19. The total tumor volume (TTV) and enhancing tumor volume (ETV) were assessed by an independent radiological reader with 1 year of experience with the software (ZZ). Three-dimensional segmentation-masks of the tumors were created to determine TTV and ETV using qEASL (Figure 2). Generally, the area within the segmentation-mask is considered the TTV and expressed in cubic centimeters (cm3). ETV was assessed using qEASL calculation in cubic centimeters (cm3)11. Initially, axial native MRI T1-weighted images were subtracted from axial enhancing phase images to remove false-positive background enhancement. After subtraction, 3D tumor segmentation-masks were used to select a region of interest (ROI) consisting of a 1 cm3 sized cube placed manually within non-tumor liver parenchyma as described previously in the literature 35. The ROI within the background liver parenchyma sets a cut-off value based on intensity that is used as a reference to calculate the ETV within the segmentation-masks of the tumors (Figure 2). After setting the ROI, the software automatically generates a color map of enhancing regions within the segmented 3D tumor mask. The non-enhancing and necrotic areas of the tumor were represented in blue, whereas enhancing and thus viable parts of the tumor were represented in red. The quantitative output resulted in volumetric values indicative of tumor enhancement.
To evaluate the enhancing tumor burden (ETB) within the liver, the total liver volume (TLV) was calculated using a software prototype (Medisys, Philips Research), that automatically generates a segmentation-mask of the entire liver. The software allows contraction and expansion of the segmentation-mask around control-points within the liver or its contour. Thereby, the mask can be manually adjusted by the reader to fully include all anatomical parts of the organ. The true volume of the segmented liver will be calculated and enunciated in cubic centimeters (cm3). The ETB (%) was calculated using the following formula:
This formula takes into account the ETV in relation to the TLV by calculating their ratio. For comparative purposes, patients were divided into low tumor burden (LTB) and high tumor burden (HTB) groups based on cut-off points defined for each 1D, 2D, ad 3D method.
Statistical analysis
Statistical analysis of the data was performed using SAS (SAS Institute, Version 9.4.3) and IBM SPSS Statistics (IBM, Version 23.0). Qualitative variables were presented as absolute numbers and percentages. Continuous variables were described by using mean ± standard deviation or median (range). Additionally, the Cox proportional hazard model was used to determine the predictive value of TTV, ETV, and ETB. Survival was calculated based on the interval between the date of embolization and death or last known alive date. The OS and cumulative survival analysis were calculated and represented using Kaplan-Meier curves, and the log-rank test was utilized to further contrast these survival curves. Q statistics was used to estimate the most significant cutoff values for each tumor assessment method, then the best area under the curve calculated by receiver operating characteristic (ROC) curve analysis was confirmed and utilized for every tumor assessment method to determine the optimal cutoff point for patient categorization into LTB and HTB groups, based on improved survival after cTACE. The demographic characteristics, child-Pugh classification, and tumor stage were considered for multivariate analysis. A p-value of less than 0.05 was considered statistically significant.