Abstract
Background and Aims: To investigate the impact of MR bias field correction on response determination and survival prediction using volumetric tumor enhancement analysis in patients with infiltrative hepatocellular carcinoma, after transcatheter arterial chemoembolization (TACE).
Methods: This study included 101 patients treated with conventional or drug-eluting beads TACE between the years of 2001 and 2013. Semi-automated 3D quantification software was used to segment and calculate the enhancing tumor volume (ETV) of the liver with and without bias-field correction on multi-phasic contrast-enhanced MRI before and 1-month after initial TACE. ETV (expressed as cm3) at baseline imaging and the relative change in ETV (as % change, ETV%) before and after TACE were used to predict response and survival, respectively. Statistical survival analyses included Kaplan-Meier curve generation and Cox proportional hazards modeling. Q statistics were calculated and used to identify the best cut-off value for ETV to separate responders and non-responders (ETV cm3). The difference in survival was evaluated between responders and non-responders using Kaplan-Meier and Cox models.
Results: MR bias field correction correlated with improved response calculation from baseline MR as well as survival after TACE; using a 415 cm3 cut-off for ETV at baseline (hazard ratio: 2.00, 95% confidence interval: 1.23-3.26, p=0.01) resulted in significantly improved response prediction (median survival in patients with baseline ETV <415 cm3: 19.66 months vs. ≥415 cm3: 9.21 months, p<0.001, log-rank test). A ≥41% relative decrease in ETV (hazard ratio: 0.58, 95%confidence interval: 0.37-0.93, p=0.02) was significant in predicting survival (ETV ≥41%: 19.20 months vs. ETV <41%: 8.71 months, p=0.008, log-rank test). Without MR bias field correction, response from baseline ETV could be predicted but survival after TACE could not.
Conclusions: MR bias field correction improves both response assessment and accuracy of survival prediction using whole liver tumor enhancement analysis from baseline MR after initial TACE in patients with infiltrative hepatocellular carcinoma.
Keywords
3D volume measurement,
Infiltrative HCC,
TACE,
Segmentation,
Bias field correction
Introduction
Hepatocellular carcinoma (HCC) is the fifth most common cancer and the second leading cause of cancer-related death worldwide.1 Patients in the United States and Europe with this type of cancer have the third lowest 5-year survival rate of all other cancer types.2 The morphology of HCC can be classified into three groups: multi-nodular, massive, and infiltrative.3 Infiltrative HCC represents 7-20% of all HCC cases4 and it is characterized by multifocal growth pattern, inhomogeneous expansion into liver parenchyma, as well as portal vein thrombosis in most cases, resulting in an advanced stage classification at initial diagnosis and a higher mortality.5,6 Thus, most patients with infiltrative HCC are ineligible for curative therapies, like ablation, resection, or transplantation.7 Catheter-based intra-arterial therapies, such as trans-arterial chemoembolization (TACE), are commonly used as palliative treatment in patients with infiltrative HCC.
In our previous work,8 we demonstrated that, for patients with unresectable HCC treated with TACE, 3D quantitative enhancement-based assessment (quantitative European Association for the Study of the Liver, qEASL) could define treatment response and non-response better than traditional uni- and bi-dimensional methods, such as RECIST, mRECIST, World Health Organization and EASL. Moreover, a significant reduction of enhancing tumor volume (ETV) was found to correspond with a significant improvement in overall survival (OS) according to qEASL.9 Hence, the early and reliable identification of responder (R) and non-responder (NR) by qEASL after TACE enables treatment decisions to be made sooner.10,11
MR imaging is the reference standard and routinely used in clinical practice for the diagnosis of HCC and the assessment of TACE therapy. However, due to imperfections in the radio frequency coil, eddy currents or patient anatomy, bias field effects, presenting as a slowly varying intensity variation across the image, are common during MR imaging. The appearance of bias field is often seen as varying bands of artificial image hyper- and hypo-enhancement across the MR image.12 Put another way, an image with uniform intensity would appear to have artificial regions of increased and/or decreased image enhancement. Such image artifact could limit the accuracy of image interpretation, segmentation, and quantification.13,14 A characteristic presentation of infiltrative HCC is its diffuse appearance/ill-defined morphology that typically encompasses a large portion of the liver with absence of discrete tumor-liver parenchyma boundaries. As such, the MR bias field artifact can interfere with and exacerbate the already challenging visualization, characterization, and quantification of infiltrative HCC.11,14 For example, regions of hyperenhancement could be due to the contrast agent uptake in the HCC, the bias field artifact, or a combination of the two.
The purpose of our study was to investigate the utility for MR bias field correction in improving response calculation among patients with infiltrative HCC and predict their survival after initial TACE using volumetric tumor enhancement analysis.
Methods
Study cohort
This retrospective, single-institution study (The Johns Hopkins Hospital, Baltimore, MD, USA) was compliant with the Health Insurance Portability and Accountability Act and approved by the institutional review board. Informed consent was waived. A total of 397 patients diagnosed with infiltrative HCC and who underwent conventional (Lipiodol®) TACE or drug-eluting beads TACE in our institution between August 2001 and November 2013 were identified. Patients were included in this study if they underwent contrast-enhanced multi-phasic MR imaging within 1 month before and 6 weeks after TACE. Eighty-six patients were excluded because of inadequate imaging, including missing phases or imaging artifact due to motion or metal. Furthermore, 210 patients were excluded because of missing contrast agent information or use of different contrast agents between the baseline and follow-up MR. At baseline, all of the included patients were therapy-naïve with respect to systemic therapy with tyrosine kinase inhibitors and none received systemic therapy during the initial course of TACE treatment until the time point of first imaging follow-up. Finally, a total of 101 patients were included into the final analysis.
TACE procedure
The decision to treat using TACE was determined by a multidisciplinary tumor board according to our standard institutional protocol. For conventional TACE, an emulsion containing 50 mg doxorubicin (Adriamycin; Pharmacia & Upjohn, Peapack, NJ, USA) and 10 mg mitomycin C in a 1:>1 mixture with iodized oil (Lipiodol®; Guerbet, Aulnay-sous-Bois, France) was infused and followed by 100-300-μm microspheres (Embosphere Microspheres; Biosphere Medical, Rockland, MA, USA) till arterial inflow was substantially reduced. For the drug-eluting beads TACE, 100–300-μm DC Beads (Biocompatible/BTG, Surrey, England) were loaded with up to 100 mg of doxorubicin hydrochloride (25 mg/mL) and mixed with 4 mL of Oxilan (300 mg iodine per mL; Guerbet, Bloomington, IN, USA) to form a 4 mL solution. The technical end-point was reduction of tumor arterial inflow to between 2-5 heart beats. Complete occlusion was avoided in order to maintain arterial patency for repeat treatment.
MR imaging protocol
All patients underwent MR imaging within 4 weeks before and 6 weeks after the first TACE, on a 1.5-Tesla whole-body MR scanner (Magnetom Avanto; Siemens Medical Solutions, Erlangen, Germany). A standardized MR liver protocol was used with the same imaging parameters to ensure consistency in image acquisition and timing. The imaging protocol included: (1) axial T2-weighted fast spin-echo images (TR/TE, 5000/100 ms; matrix size, 256×256; slice thickness, 8 mm; interslice gap, 2 mm; receiver bandwidth, 32 kHz); (2) axial single-shot breath-hold gradient echo diffusion-weighted echo-planar images (TR/TE, 5000-6500/110 ms; matrix size, 128×128; slice thickness, 8 mm; interslice gap, 2 mm; b-value, 500 s/mm2; receiver bandwidth, 64 kHz); (3) axial breath-hold unenhanced and contrast-enhanced (0.1 mmol/kg of body weight, intravenous injection of gadodiamide [Omniscan, General Electric Healthcare, Princeton, NJ, USA]) T1-weighted three-dimensional fat suppressed spoiled gradient-echo images (TR/TE, 5.1/1.2 ms; field of view, 320-400 mm2; matrix size, 192×160; slice thickness, 4-6 mm; receiver bandwidth, 64kHz; flip angle, 15°) in the pre-contrast injection, post-contrast injection arterial, portal venous and delayed phases (20, 70 and 180 s after intravenous contrast administration, respectively; a fixed time-delay technique was used). Finally, the pre-contrast and 20 s post-contrast scans were used for MR image processing and 3D enhancement quantification.
MR image processing: Registration, segmentation and bias field correction
All images were presented as Digital Imaging in Communications and Medicine (“DICOM”) files. Image assessment was conducted independently by two radiologists (CL and SS, with 8 and 2 years of experience in abdominal MR) who did not participate in the TACE procedures and were blinded to survival outcome.
An automated 3D deformable registration (BioImage Suite, Yale University, New Haven, CT, USA) was used to register the pre- and 20 second post-contrast scans as a preprocessing step,15,16 in order to reduce inter-phase hepatic displacement and improve image quality. The whole liver from the arterial phase MR image was segmented in 3D semi-automatically using a prototype software which generated a 3D segmentation mask.8,17 Care was taken to exclude the major hepatic vessels, such as the hepatic artery and portal vein (Fig. 1A). Finally, the MR bias field18,19 correction using the software (BioImage Suite) to remove the intensity variation was applied to the 3D segmented liver in order to mitigate the effects of shading artifacts (Fig. 1C). This is a two-step procedure. First, the slice-by-slice inhomogeneity (shading artifacts inside and outside of the liver) in all three (x, y, z) directions is normalized for the entire MR image by estimating a linear regression between any adjacent slices and removing any changes.19 Assuming that at macro-resolution (∼20 mm) the liver intensity is approximately homogeneous, the liver (as defined by the liver segmentation above) in the resultant image above is divided by a heavily blurred (macro-resolution) version of itself in order to eliminate additional intensity inhomogeneity (Fig. 1D).
3D enhancement quantification
From the whole liver 3D segmentation, the ETV of the entire liver was calculated (qEASL) three times, and the results averaged for before and after bias field correction. This validated method has high inter-reader reproducibility and radiology-pathology validation.8,17,20 Briefly, the procedure was as follows:
The 3D registered pre-contrast T1-weighted MR image was subtracted from the contrast-enhanced arterial phase image, in order to remove contributions of background enhancement.21,22 Second, the 3D segmentation mask of the entire liver was transposed to the subtracted images. Third, a 1-cm3 3D region of interest23 was manually placed in healthy liver parenchyma that would have a coefficient of variation less than 25%. Care was taken to place the region of interest away from tumors, blood vessels, gallbladder, and liver capsule during region of interest placement. Viable tumor tissue was defined as enhancement that was more than twice the standard deviation of the region of interest average signal intensity.23 A color map visualization in which the brightness and color indicating the strength and location of the enhanced signal (Fig. 1E) was also made. To improve results, we repeated the ETV calculation three times based on three different region of interest placements for each patient and the three ETV measurements were averaged. Based on the ETV results, the absolute and relative change in ETV was calculated as:
Absolute ETV change = Baseline ETV - ETV after TACE
Relative ETV change (ETV%) = Absolute ETV change / Baseline ETV
The ETV measurements/calculations were carried out with and without MR bias field correction.
Statistical analysis
The baseline characteristics were summarized using descriptive statistics. Paired t-test was used to compare differences with and without MR bias field correction within the same patient. OS, measured from the time of TACE to death, was estimated using the Kaplan-Meier method. Cut-off points were calculated to identify Rs and NRs.24 Differences between response groups were assessed using log-rank test. Univariate and multivariate Cox modelings were performed to identify prognostic factors. A two-sided p-value of less than 0.05 was considered statistically significant. All statistical analyses were carried out using SAS 9.4 (SAS Institute, Cary, NC, USA).
Results
Patient characteristics
The mean age at TACE treatment was 60.3 years for the 101 patients included in our study and 76 of them (75%) were male. Among them, 42 patients (42%) had portal venous invasion and 22 patients (22%) had extra-hepatic metastasis. Also, 67 patients were treated using Lipiodol® and 34 patients using drug-eluting beads. The mean time (standard deviation) between baseline imaging and the first TACE procedure was 13.36 (13.95) days, and the mean time between TACE and follow-up imaging was 26.88 (9.23) days. None of the patients died before the first follow-up. Patients were followed up for a median (interquartile range) of 10.8 (4.0-20.7) months. The median OS of the study cohort was 12.8 months (95% confidence interval: 10.29-16.08) and by the end of observation date (February 2016), 78 patients (77%) were deceased. Additional baseline characteristics are summarized in Table 1.
Table 1.Baseline patients and tumor characteristics
| Variables | n (%) |
| Demographics | |
| Total number of patients | 101 |
| Age at baseline | |
| <65 years | 72 (71) |
| ≥65 years | 29 (29) |
| Sex | |
| Male | 76 (75) |
| Female | 25 (25) |
| Vital status at time of analysis | |
| Alive | 23 (22) |
| Deceased | 78 (78) |
| Follow-up time in months | |
| Mean (range) | 17.1 (0.1-140.4) |
| Median (interquartile range) | 10.8 (4.0-20.7) |
| Etiology | |
| Hepatitis B virus | 16 (16) |
| Hepatitis C virus | 45 (45) |
| Hepatitis B virus+hepatitis C virus | 2 (2) |
| Human immunodeficiency virus | 5 (5) |
| Alcohol | 35 (35) |
| Non-alcoholic steatohepatitis | 7 (7) |
| Cryptogenic | 7 (7) |
| Cirrhosis | 36 (36) |
| Ethnicity | |
| White | 65 (64) |
| African American | 16 (16) |
| Hispanic | 1 (1) |
| Other | 8 (8) |
| Unknown | 11 (11) |
| Barcelona Clinic Liver Cancer stage | |
| B | 29 (29) |
| C | 69 (68) |
| D | 3 (3) |
| Child-Pugh class | |
| A | 68 (67) |
| B | 32 (32) |
| C | 1 (1) |
| ECOG performance score | |
| 0 | 56 (55) |
| 1 | 39 (39) |
| 2 | 4 (4) |
| 3 | 2 (2) |
| Tumor characteristics | |
| Number of lesions | |
| 1 | 18 (18) |
| 2 | 19 (19) |
| 3 | 11 (11) |
| 4 or more | 53 (52) |
| Tumor size in cm | |
| Mean (range) | 8.6 (1.5-19.5) |
| Median (interquartile range) | 8.0 (5.0-11.5) |
| Portal invasion | |
| Yes | 42 (42) |
| No | 59 (58) |
| Extrahepatic invasion | |
| Yes | 22 (22) |
| No | 79 (78) |
| TACE type | |
| Lipiodol® | 67 (66) |
| Drug-eluting beads | 34 (34) |
| Number of TACE sessions | |
| 1 | 48 (48) |
| 2 | 19 (18) |
| 3 | 13 (13) |
| 4 or more | 21 (21) |
Imaging analysis and volumetric calculation
Paired t-test was used to compare the ETV with and without MR bias field correction. As shown in Table 2, the mean baseline ETV was 823.98 cm3 without bias field correction and 502.71 cm3 with bias field correction (p<0.001); similarly, after TACE treatment, mean ETV was 629.91 cm3 and 322.05 cm3 (p<0.001) with and without MR bias field correction, respectively, denoting that the application of MR bias field could impact the evaluation and the measurement of the tumor. Furthermore, we also calculated absolute ETV change and relative ETV change (ETV%). Without MR bias field correction, mean absolute ETV decrease was 194.06 (-1012.29-1671.84) cm3 and mean ETV% was 0.10 (-5.35-0.83), whereas with MR bias field correction, mean absolute ETV decrease was 180.66 (27.08-627.70) cm3 and mean ETV% was 0.37 (0.06-0.62), respectively.
Table 2.Comparison of ETV with and without MR bias field correction
| Mean (range) in cm3 | p |
| Without MR bias field correction | With MR bias field correction |
| Baseline ETV | 823.98 (62.52-3560.86) | 502.71 (141.41-1294.27) | <0.001 |
| ETV after TACE | 629.91 (46.70-3380.53) | 322.05 (109.50-1016.39) | <0.001 |
With the results above, the Cox proportional hazards model was used to determine whether ETV or volume change would influence prognosis of patients with infiltrative HCC. As shown in Table 3, both baseline ETV and ETV% with MR bias field correction correlated with prognosis. Without MR bias field correction, only baseline ETV was correlated with prognosis. Therefore, subsequent analysis was done with MR bias field corrected images only. First, the cut-off points for baseline ETV and ETV% for prognosis prediction were identified. To determine such cut-off points for baseline ETV and ETV% with MR bias field correction, the Q statistics test was used. The LOESS smoothed line of Martingale residuals was not straight nor close to zero, which implied that it would be better to covert ETV or ETV% into an indicator with appropriate cut-off points instead of keeping them linear in the Cox model. Therefore, cut-off points were calculated with the method from Mandrekar’s paper21 to identify Rs and NRs. For baseline ETV, the cut-off point was 415 cm3 (p=0.003) (Fig. 2A), and for ETV% the cut-off point was 41% (p=0.01) (Fig. 2B).
Table 3.Hazard ratio for baseline ETV, absolute ETV change and ETV% with and without MR bias field correction
| Variables | Without MR bias field correction | | With MR bias field correction |
| Hazard ratio (95% CI) | P | Hazard ratio (95% CI) | P |
| Baseline ETV | 1.061 (1.029, 1.093) | <0.001 | | 1.161 (1.066, 1.265) | <0.001 |
| Absolute ETV change | 1.010 (0.982, 1.039) | 0.49 | | 1.038 (0.932, 1.157) | 0.50 |
| ETV% | 1.003 (0.979, 1.028) | 0.79 | | 0.663 (0.535, 0.822) | <0.001 |
Survival analysis
With the cut-off points determined, the patients were separated into different groups to check whether the cut-off points could define Rs and NRs to the treatment. As illustrated in Fig. 3A, 43 patients (42.6%) were in the ‘Baseline ETV <415 cm3’ group and 58 patients (57.4%) were in the ‘Baseline ETV ≥415 cm3’ group; median survival time for them was 19.66 (15.58-31.79) months and 9.21 (6.21-11.64) months, respectively (p<0.001). Fifty-six patients (55.4%) were in the ‘ETV% <41%’ group and forty-five patients (44.6%) were in the ‘ETV%≥41%’ group; median survival time for them was 8.71 (6.12-11.21) months and 19.20 (14.53-21.17) months, respectively (p=0.008) (Fig. 3B).
Univariate and multivariate analyses
Univariate and multivariate analyses were carried out following adjustment for factors such as age, gender, ethnicity, TACE type (Lipiodol® or drug-eluting beads), Barcelona Clinic Liver Cancer stage and Child–Pugh score. Baseline ETV and ETV% were identified to be prognostic factors for patients with infiltrative HCC. According to the univariate analysis, the mortality rate for patients in the ‘Baseline ETV ≥415 cm3’ group was twice as high as those in the ‘Baseline ETV ≥415 cm3’ group (p=0.01) (Supplementary Table 1). The mortality rate for patients in the ‘ETV% <41%’ group was 42% lower than patients in group ‘ETV% <41%’ (p=0.02) (Supplementary Table 2). Of note, univariate analysis did not demonstrate the type of TACE (Lipiodol® or drug-eluting beads) to be a statistically significant prognostic factor.
Discussion
This study concludes that with MR bias field correction, it is possible to use whole liver tumor enhancement analysis to predict treatment response from baseline MR and survival after initial TACE in patients with infiltrative HCC. Specifically, baseline ETV could be used to identify candidates who may benefit from TACE before treatment, and ETV% change after first TACE could be used to predict survival.
Characterizing infiltrative HCC is challenging because it frequently presents as a large, diffuse process.3–5 As such, image assessment of infiltrative HCC is typically done through visual observation, without objective quantitative method and standard. In a previous study,25 the image intensity of the lesion was measured in 2D and they were unable to predict survival and treatment response. In our work, the measurement is based on the whole liver in 3D and we were able to find statistical significance in predicting survival and therapy response even though we used an imaging method (i.e. MR) that has non-calibrated image intensities (as opposed to calibrated intensities like Hounsfield Units in computed tomography). The main feature is that our assessment was done on a whole liver basis using MR bias field correction. The intensity inhomogeneity observed in MR images due to the bias field produced by imperfections in the radio frequency coil, eddy currents, as well as object susceptibility make quantitative image analysis challenging.12–14 Overcoming this limitation was essential to mitigating the influence of the bias field and in doing so enabled statistically significant response and survival predictions. Our hypothesis was tested by comparing the survival and response predictive ability of volumetric enhancement with and without MR bias field correction. It was found that with bias correction, both ETV at baseline and ETV% change after TACE could be used as predictive factors for prognosis; whereas, without MR bias field correction, only baseline ETV could predict prognosis.
In our previous work,8,9 qEASL was used for volumetric measurement of unifocal and multifocal HCC and it was determined that a volumetric ETV threshold of 65 cm3 for unifocal HCC and 45 cm3 for multifocal HCC could predict response. However, in the current study, the previous cut-offs could not be applied due to the characteristics of infiltrative HCC. Thus, a unique cut-point determination method using Q-statistics24 was employed to determine the cut-off for baseline ETV and ETV% that could predict treatment response in patients with infiltrative HCC. It was found that patients with a baseline ETV less than 415 cm3 of infiltrative HCC tumor burden could be better candidates for TACE, and for after TACE treatment, that patients with ETV% over 41% were found to have a longer survival time. This information can help to identify patients with infiltrative HCC who could benefit from treatments other than TACE before initiating therapy. It could further help identify NRs to TACE so that changes in therapy management could happen sooner, especially for those patients who failed to respond to the initial and second TACE.26 Additionally, these findings could help identify patients that would potentially benefit from earlier onset of systemic therapy (e.g., with sorafenib or checkpoint inhibitors) despite being formally classified as Barcelona Clinic Liver Cancer stage B.
There are several limitations to the current study. First, it was a single-institutional retrospective study with a limited cohort size. Second, the whole tumor signal is summarized into one result. Although the main branches of the portal vein and hepatic artery were excluded, smaller vessels with intensity similar to the tumor were inevitably included, possibly leading to unspecific background signal.
In conclusion, MR bias field correction has the potential to improve 3D quantitative tumor assessment techniques and may add to the utility of these techniques in predicting both response and survival outcomes after initial TACE in patients with infiltrative HCC using whole liver tumor enhancement analysis. Having this predictive information, especially at time of staging/baseline imaging, could potentially improve allocation of patients to a more appropriate TACE schedule or therapy selection early during the clinical decision making process.
Supplementary information
Supplementary Table 1
Univariate and multivariate analyses results for baseline ETV with MR bias field correction
(DOCX)
Supplementary Table 2
Univariate and multivariate analyses results for ETV% with MR bias field correction
(DOCX)
Abbreviations
- ETV:
enhancing tumor volume
- HCC:
hepatocellular carcinoma
- NR:
non-responder
- OS:
overall survival
- R:
responder
- TACE:
transcatheter arterial chemoembolization
Declarations
Funding
This research was funded by grant from NIH/NCI (R01 CA160771).
Conflict of interest
MingDe Lin is a former employee of Philips Research North America and is now an employee of Visage Imaging, Inc. Julius Chapiro declares relation to Koninklijke Philips, Guerbet SA, Eisai Co. Xenophon Papademetris declares relation to Brain Electrophysiology Laboratory Company (no relation to this work). The other authors have no conflict of interests related to this publication.
Authors’ contributions
Study design (CL, SS, ML, JC), performance of experiments (CL, SS, XP, DDM), analysis and interpretation of data (CL, SS, GG, YD, ML, JC), manuscript writing (CL, SS, XP, GG, YD, ML, JC, CG, KH), critical revision of the manuscript (CL, ML, JC), statistical analysis (GG, YD), critical funding (ML, JC), administration (XW, CG, KH), technical or material support (ML, XP).