Home
JournalsCollections
For Authors For Reviewers For Editorial Board Members
Article Processing Charges Open Access
Ethics Advertising Policy
Editorial Policy Resource Center
Company Information Contact Us
OPEN ACCESS

Predicting Infiltrative Hepatocellular Carcinoma Patient Outcome Post-TACE: MR Bias Field Correction Effect on 3D-quantitative Image Analysis

  • Cuihong Liu1,2,
  • Susanne Smolka2,3,
  • Xenophon Papademetris4,
  • Duc Do Minh2,3,
  • Geliang Gan5,
  • Yanhong Deng5,
  • MingDe Lin2,
  • Julius Chapiro2,*,
  • Ximing Wang*,6,7,
  • Christos Georgiades8 and
  • Kelvin Hong8
Journal of Clinical and Translational Hepatology   2020;8(3):292-298

doi: 10.14218/JCTH.2020.00054

Received:

Revised:

Accepted:

Published online:

 Author information

Citation: Liu C, Smolka S, Papademetris X, Minh DD, Gan G, Deng Y, et al. Predicting Infiltrative Hepatocellular Carcinoma Patient Outcome Post-TACE: MR Bias Field Correction Effect on 3D-quantitative Image Analysis. J Clin Transl Hepatol. 2020;8(3):292-298. doi: 10.14218/JCTH.2020.00054.

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).

Image processing workflow.
Fig. 1.  Image processing workflow.

Whole liver 3D segmentation (red outline) on MR image (A) and its volumetric rendering (B). Before (C) and after (D) MR bias field correction. (E) Enhancing tumor volume colormap overlay. Note the heterogeneous tumor enhancement.

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

Variablesn (%)
Demographics
Total number of patients101
Age at baseline
 <65 years72 (71)
 ≥65 years29 (29)
Sex
 Male76 (75)
 Female25 (25)
Vital status at time of analysis
 Alive23 (22)
 Deceased78 (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 virus16 (16)
 Hepatitis C virus45 (45)
 Hepatitis B virus+hepatitis C virus2 (2)
 Human immunodeficiency virus5 (5)
 Alcohol35 (35)
 Non-alcoholic steatohepatitis7 (7)
 Cryptogenic7 (7)
Cirrhosis36 (36)
Ethnicity
 White65 (64)
 African American16 (16)
 Hispanic1 (1)
 Other8 (8)
 Unknown11 (11)
Barcelona Clinic Liver Cancer stage
 B29 (29)
 C69 (68)
 D3 (3)
Child-Pugh class
 A68 (67)
 B32 (32)
 C1 (1)
ECOG performance score
 056 (55)
 139 (39)
 24 (4)
 32 (2)
Tumor characteristics
Number of lesions
 118 (18)
 219 (19)
 311 (11)
 4 or more53 (52)
Tumor size in cm
 Mean (range)8.6 (1.5-19.5)
 Median (interquartile range)8.0 (5.0-11.5)
Portal invasion
 Yes42 (42)
 No59 (58)
Extrahepatic invasion
 Yes22 (22)
 No79 (78)
TACE type
 Lipiodol®67 (66)
 Drug-eluting beads34 (34)
Number of TACE sessions
 148 (48)
 219 (18)
 313 (13)
 4 or more21 (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 cm3p
Without MR bias field correctionWith MR bias field correction
Baseline ETV823.98 (62.52-3560.86)502.71 (141.41-1294.27)<0.001
ETV after TACE629.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

VariablesWithout MR bias field correctionWith MR bias field correction
Hazard ratio (95% CI)PHazard ratio (95% CI)P
Baseline ETV1.061 (1.029, 1.093)<0.0011.161 (1.066, 1.265)<0.001
Absolute ETV change1.010 (0.982, 1.039)0.491.038 (0.932, 1.157)0.50
ETV%1.003 (0.979, 1.028)0.790.663 (0.535, 0.822)<0.001
Cut-off points for MR bias field corrected baseline ETV and ETV%.
Fig. 2.  Cut-off points for MR bias field corrected baseline ETV and ETV%.

(A) ETV was 415 cm3 (p=0.003) and (B) ETV% was 41% (p=0.010), respectively.

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).

Kaplan-Meier survival curves for baseline ETV and ETV%.
Fig. 3.  Kaplan-Meier survival curves for baseline ETV and ETV%.

The cut-offs were 415 cm3 and 41% for choosing TACE candidates and predicting treatment response, respectively.

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).

References

  1. EASL Clinical Practice Guidelines: Management of hepatocellular carcinoma. J Hepatol 2018;69:182-236 View Article PubMed/NCBI
  2. Forner A, Llovet JM, Bruix J. Hepatocellular carcinoma. Lancet 2012;379:1-7 View Article PubMed/NCBI
  3. Sereni CP, Rodgers SK, Kirby CL, Goykhman I. Portal vein thrombus and infiltrative HCC: a pictoral review. Abdom Radiol (NY) 2017;42:159-170 View Article PubMed/NCBI
  4. Reynolds AR, Furlan A, Fetzer DT, Sasatomi E, Borhani AA, Heller MT. Infiltrative hepatocellular carcinoma: what radiologists need to know. Radiographics 2015;35:371-386 View Article PubMed/NCBI
  5. Kneuertz PJ, Demirjian A, Firoozmand A, Corona-Villalobos C, Bhagat N, Herman J. Diffuse infiltrative hepatocellular carcinoma: assessment of presentation, treatment, and outcomes. Ann Surg Oncol 2012;19:2897-2907 View Article PubMed/NCBI
  6. Demirjian A, Peng P, Geschwind JF, Cosgrove D, Schutz J, Kamel IR. Infiltrating hepatocellular carcinoma: seeing the tree through the forest. J Gastrointest Surg 2011;15:2089-2097 View Article PubMed/NCBI
  7. Dai QS, Gu HL, Ye S, Zhang YJ, Lin XJ, Lau WY. Transarterial chemoembolization vs. conservative treatment for unresectable infiltrating hepatocellular carcinoma: A retrospective comparative study. Mol Clin Oncol 2014;2:1047-1054 View Article PubMed/NCBI
  8. Fleckenstein FN, Schernthaner RE, Duran R, Sohn JH, Sahu S, Zhao Y. 3D Quantitative tumour burden analysis in patients with hepatocellular carcinoma before TACE: comparing single-lesion vs. multi-lesion imaging biomarkers as predictors of patient survival. Eur Radiol 2016;26:3243-3252 View Article PubMed/NCBI
  9. Chapiro J, Duran R, Lin M, Schernthaner RE, Wang Z, Gorodetski B. Identifying staging markers for hepatocellular carcinoma before transarterial chemoembolization: Comparison of three-dimensional quantitative versus non-three-dimensional imaging markers. Version 2. Radiology 2015;275:438-447 View Article PubMed/NCBI
  10. Tacher V, Lin M, Duran R, Yarmohammadi H, Lee H, Chapiro J. Comparison of existing response criteria in patients with hepatocellular carcinoma treated with transarterial chemoembolization using a 3D quantitative approach. Version 2. Radiology 2016;278:275-284 View Article PubMed/NCBI
  11. Zhao Y, Duran R, Bai W, Sahu S, Wang W, Kabus S. Which criteria applied in multi-phasic CT can predict early tumor response in patients with hepatocellular carcinoma treated using conventional TACE: RECIST, mRECIST, EASL or qEASL?. Cardiovasc Intervent Radiol 2018;41:433-442 View Article PubMed/NCBI
  12. Ding H, Johnson T, Lin M, Le HQ, Ducote JL, Su MY. Breast density quantification using magnetic resonance imaging (MRI) with bias field correction: a postmortem study. Med Phys 2013;40:122305 View Article PubMed/NCBI
  13. Zhang K, Liu Q, Song H, Li X. A variational approach to simultaneous image segmentation and bias correction. IEEE Trans Cybern 2015;45:1426-1437 View Article PubMed/NCBI
  14. Despotović I, Goossens B, Philips W. MRI segmentation of the human brain: challenges, methods, and applications. Comput Math Methods Med 2015;2015:450341 View Article PubMed/NCBI
  15. Dura E, Domingo J, Ayala G, Martí-Bonmatí L. Evaluation of the registration of temporal series of contrast-enhanced perfusion magnetic resonance 3D images of the liver. Comput Methods Programs Biomed 2012;108:932-945 View Article PubMed/NCBI
  16. Chitphakdithai N, Chiang VL, Duncan JS. Non-rigid registration of longitudinal brain tumor treatment MRI. Conf Proc IEEE Eng Med Biol Soc 2011;2011:4893-4896 View Article PubMed/NCBI
  17. Sahu S, Schernthaner R, Ardon R, Chapiro J, Zhao Y, Sohn JH. Imaging biomarkers of tumor response in neuroendocrine liver metastases treated with transarterial chemoembolization: Can enhancing tumor burden of the whole liver help predict patient survival? Version 2. Radiology 2017;283:883-894 View Article PubMed/NCBI
  18. Zhang Y, Brady M, Smith S. Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans Med Imaging 2001;20:45-57 View Article PubMed/NCBI
  19. Petersen KF, Dufour S, Savage DB, Bilz S, Solomon G, Yonemitsu S. The role of skeletal muscle insulin resistance in the pathogenesis of the metabolic syndrome. Version 2. Proc Natl Acad Sci U S A 2007;104:12587-12594 View Article PubMed/NCBI
  20. Chapiro J, Wood LD, Lin M, Duran R, Cornish T, Lesage D. Radiologic-pathologic analysis of contrast-enhanced and diffusion-weighted MR imaging in patients with HCC after TACE: diagnostic accuracy of 3D quantitative image analysis. Version 2. Radiology 2014;273:746-758 View Article PubMed/NCBI
  21. Duran R, Chapiro J, Frangakis C, Lin M, Schlachter TR, Schernthaner RE. Uveal melanoma metastatic to the liver: The role of quantitative volumetric contrast-enhanced MR imaging in the assessment of early tumor response after transarterial chemoembolization. Transl Oncol 2014;7:447-455 View Article PubMed/NCBI
  22. Lin M, Pellerin O, Bhagat N, Rao PP, Loffroy R, Ardon R. Quantitative and volumetric European Association for the Study of the Liver and Response Evaluation Criteria in Solid Tumors measurements: feasibility of a semiautomated software method to assess tumor response after transcatheter arterial chemoembolization. J Vasc Interv Radiol 2012;23:1629-1637 View Article PubMed/NCBI
  23. Chockalingam A, Duran R, Sohn JH, Schernthaner R, Chapiro J, Lee H. Radiologic-pathologic analysis of quantitative 3D tumour enhancement on contrast-enhanced MR imaging: a study of ROI placement. Eur Radiol 2016;26:103-113 View Article PubMed/NCBI
  24. Mandrekar JN, Mandrekar SJ, Cha SS. Cutpoint determination methods in survival analysis using SAS®. Available from: http://support.sas.com/resources/papers/proceedings/proceedings/sugi28/261-28.pdf View Article PubMed/NCBI
  25. Koulakian H, Allaham W, Vilgrain V, Ronot M. Non-measurable infiltrative HCC: is post-contrast attenuation on CT a sign of tumor response?. Eur Radiol 2019;29:4389-4399 View Article PubMed/NCBI
  26. Georgiades C, Geschwind JF, Harrison N, Hines-Peralta A, Liapi E, Hong K. Lack of response after initial chemoembolization for hepatocellular carcinoma: does it predict failure of subsequent treatment? Version 2. Radiology 2012;265:115-123 View Article PubMed/NCBI
  27. Treilhard J, Smolka S, Staib L, Chapiro J, Lin M, Shakirin G. Liver tissue classification in patients with hepatocellular carcinoma by fusing structured and rotationally invariant context representation. Med Image Comput Comput Assist Interv 2017;10435:81-88 View Article PubMed/NCBI
  28. Zhang F, Yang J, Nezami N, Laage-Gaupp F, Chapiro J, De Lin M. Liver tissue classification using an auto-context-based deep neural network with a multi-phase training framework. Patch Based Tech Med Imaging (2018) 2018;11075:59-66 View Article PubMed/NCBI