Abstract
Background
Repeated transarterial chemoembolization (TACE) can be associated with loss of its efficacy and subsequent tumor progression.
Purpose
To identify features of gadoxetic acid-enhanced magnetic resonance imaging (MRI) associated with TACE refractoriness and to develop a prediction model for estimating the risk of TACE refractoriness.
Material and Methods
Among 1025 patients with intermediate-stage hepatocellular carcinoma (HCC) who underwent TACE as a first-line treatment during 2010–2017, 427 patients who underwent preoperative gadoxetic acid-enhanced MRI were analyzed. According to the date of initial TACE, patients were divided into the development cohort (n = 211) and the test cohort (n = 216). TACE refractoriness was determined according to the Japan Society of Hepatology guidelines. Univariable and multivariable analyses were performed to investigate the association between clinical/MRI features and TACE refractoriness. The performance of the prediction model was internally and externally assessed using the C-index of discrimination and a Hosmer-Lemeshow goodness-of-fit test for calibration.
Results
By analyzing 427 patients, we constructed a prediction model with the following independent features associated with TACE refractoriness: maximum tumor size; tumor number; peritumoral hypointensity on hepatobiliary phase (HBP); and the presence of non-hypervascular hypointense nodule on HBP. This system enabled the prediction of TACE refractoriness in the development cohort (C-index, 0.796) and the test cohort (C-index, 0.738) with good discrimination and calibration abilities.
Conclusion
The prediction model based on gadoxetic acid-enhanced MRI features in addition to the known predictors including tumor size and number can be used to estimate the risk of TACE refractoriness in patients with intermediate-stage HCC.
Keywords
Introduction
Transarterial chemoembolization (TACE) is the standard treatment for patients with intermediate-stage hepatocellular carcinoma (HCC) according to the Barcelona Clinic Liver Cancer (BCLC) staging system (1,2). However, the efficacy of TACE is reduced upon repeated use, and the rate of local recurrence after initial TACE may be as high as 80% (3). Repeated TACE is associated with deterioration of liver function (4,5). Due to hypoxia caused by TACE (6,7), recurrent tumors after TACE have a significantly shorter doubling time compared with primary HCCs (3). In clinical practice, TACE is applied not only in patients with intermediate-stage HCCs but also in those with advanced-stage HCCs as well as those who already show TACE refractoriness, due to unsatisfactory results from systemic therapy (8,9).
The recent introduction of multikinase inhibitors such as lenvatinib, regorafenib, cabozantinib, and immune checkpoint inhibitors has changed the treatment landscapes of HCC (9–11). A recent ongoing prospective trial showed a possibility for a combination of TACE and systemic therapy in improving the outcomes of intermediate-stage HCC (12). Thus, it can be beneficial for patients with high risk of TACE refractoriness, if they are quickly identified so as to allow timely switching to combination therapy or systemic therapy. In addition, patients with intermediate-stage HCC (BCLC B) are highly heterogeneous in terms of tumor stages and prognosis (13,14). Previous attempts for developing a tailored subgroup stratification for the BCLC B stage include the Bolondi subclassification (13) and the Kinki criteria (14), which are mainly based on tumor factors such as tumor size and number, and liver function similar to the BCLC staging system.
Due to its excellent diagnostic performance, gadoxetic acid-enhanced liver magnetic resonance imaging (MRI) has been increasingly used to evaluate HCC (15,16). Gadoxetic acid provides information on not only hemodynamic profiles but also the biological behaviors of HCC by being taken up into hepatocytes via the organic-anion-transporting polypeptide B1 receptor during the hepatobiliary phase (HBP) (17–21). Imaging features on gadoxetic acid-enhanced MRI provide additional useful information beyond conventional staging systems for predicting post-treatment outcomes mainly after curative treatments, such as radiofrequency ablation (RFA), surgery, or liver transplantation (22–25). Thus, imaging features on gadoxetic acid-enhanced MRI has the potential to predict responses after TACE. To our knowledge, however, little has been known for its applicability to predict TACE refractoriness.
The aim of the present study was to identify the imaging features on gadoxetic acid-enhanced MRI associated with TACE refractoriness in patients with HCC and to develop an easily applicable prediction model for estimating the risk of TACE refractoriness.
Material and Methods
The present study was approved by the institutional review board of our center; the need for informed consent was waived due to the retrospective nature of the study. Analysis and reporting were carried out by adhering to the TRIPOD (Transparent Reporting of a multivariable prediction model for Individual Prognosis Or Diagnosis) guideline (26).
Study population
We reviewed the records of 4703 consecutive patients newly diagnosed with HCC and initially treated with TACE at our institution between January 2010 and December 2017. HCCs were diagnosed according to the guidelines of the American Association for the Study of Liver Diseases (AASLD) and the European Association for the Study of the Liver (EASL), or biopsy in cases with atypical imaging features (1,2). After identifying 1025 patients with BCLC B stage, we excluded the following: (i) those who did not undergo gadoxetic acid-enhanced MRI within one month before TACE (n = 535); (ii) those whose tumor response to TACE could not be assessed due to the lack of follow-up images or clinical information (n = 23); and (iii) those who were switched to or additionally underwent other treatment options for HCC (i.e. RFA, radiation therapy, surgery, or liver transplantation) immediately after TACE (n = 40). According to the date of the initial TACE, the remaining 427 patients were divided into the development cohort and the test cohort; the 211 patients who underwent TACE during 2010–2014 were allocated to the development cohort to derive the prediction model, and the 216 patients who underwent TACE during 2015–2017 were allocated to the test cohort (Fig. 1).

Flow chart of the study population. BCLC, the Barcelona Clinic Liver Cancer staging system; HCC, hepatocellular carcinoma; MRI, magnetic resonance imaging; TACE, tranarterial chemoembolization.
The clinical and demographic characteristics of the study population, such as age, sex, etiology of chronic liver disease, and laboratory data including total bilirubin, albumin, Child-Pugh class, and alpha-fetoprotein (AFP), were collected from electronic medical records.
TACE procedures
All TACE procedures were performed by highly experienced interventional radiologists. Before chemoembolization, superior mesenteric and common hepatic angiographies were performed to assess the vascular anatomy, tumor vascularity, and portal vein patency. In all patients, cisplatin-based TACE was performed. Detailed description of procedure is presented in the Supplementary Material and Methods.
Acquisition of gadoxetic acid-enhanced MRI
MRI examinations were performed using 1.5-T (Magnetom Avanto, Siemens Medical Solutions, Erlangen, Germany) or 3-T (Magnetom Skyra, Siemens Medical Solutions, Erlangen, Germany) systems. Detailed MRI techniques and parameters for MRI sequences are summarized in the Supplementary Materials and Methods and Supplementary Table 1. The average time interval between the acquisition of pre-treatment MRI and TACE for each patient was eight days (range = 1–30 days).
Pre-procedural MRI data analysis
All pre-procedural MRIs were retrospectively reviewed by two abdominal radiologists (with >6 years of experience in abdominal MRI interpretation) in consensus, who were blinded to the clinical characteristics and follow-up imaging results after TACE. The reviewers evaluated the following MRI features that potentially reflected pathologic characteristics of HCC (17,20,27–29): (i) maximum tumor size; (ii) tumor number; (iii) morphology (nodular or multinodular vs. infiltrative); (iv) intratumoral fat (tumor area with decreased signal intensity on opposed-phase images compared with in-phase images (17)); (v) arterial peritumoral enhancement (detectable portion of crescent- or polygonal-shaped enhancement outside the tumor margin with broad contact with the tumor border on arterial phase, becoming isointense with background liver parenchyma on delayed phase (27)); (vi) enhancing capsule(peripheral rim of smooth hyperenhancement in the portal venous or transitional phase images (17,28)); (vii) peritumoral hypointensity on HBP (wedge-shaped or flame-like hypointense area of hepatic parenchyma located outside of the tumor margin on HBP images (20)); and (viii) presence of non-hypervascular hypointense nodules (NHHN; nodules showing unequivocal hypointensity on HBP images compared with the surrounding hepatic parenchyma without arterial phase hyperenhancement (29)).
Assessment of TACE refractoriness
Two abdominal radiologists that were not involved in pre-procedural MRI analyses evaluated all available follow-up CT or MRI scans performed 1–3 months after TACE to determine the radiological treatment responses of the HCCs. TACE refractoriness was defined according to the Japan Society of Hepatology (JSH) Consensus Guidelines (5): (i) two or more consecutive insufficient responses of the treated tumor (viable lesion >50%); (ii) more than two consecutive new intrahepatic lesions; (iii) continuous elevation of tumor markers immediately after TACE even with slight transient decrease; (iv) appearance of vascular invasion; and (v) appearance extrahepatic spread. For determining the elevation of tumor markers, we evaluated the level of AFP at baseline and 1–3 months after TACE and defined a continuous elevation of the tumor marker as an increase of >20% from baseline.
Statistical analysis
Continuous variables were expressed as mean ± SD or median (interquartile range); categorical variables were expressed as number or frequency. Data were analyzed with the Student’s t-test or the Mann–Whitney U test for continuous variables, and the chi-square test or Fisher’s exact test for categorical variables.
A logistic regression analysis was performed to assess the powers of clinical characteristics and MRI features for the prediction of TACE refractoriness. In the development cohort, univariable and multivariable analyses were performed to investigate the association between TACE refractoriness and various clinical and MRI factors. To develop a statistically robust scoring model predicting the refractoriness, we used a bootstrap resampling method. The variable selection process in the system was conducted as follows. First, we generated 1000 bootstrap resamples and a backward step-down variable selection procedure, based on a Wald chi-square statistic, was performed in the logistic regression model for each bootstrap sample. We then counted how many of each candidate variables remained in the 1000 bootstrap samples: if a variable appeared ≥500 times in the final model (i.e. bootstrap reliability criterion ≥50%), the variable was included in the scoring system. To allocate risk points in the scoring system, the bootstrapping method was used again to obtain bias-corrected regression coefficients. Risk points were obtained by bias-corrected regression coefficients multiplied by a reference value in the corresponding category. A reference risk factor profile was chosen by selecting a base category for each risk factor, which was assigned with 0 points in the scoring system. The total point scores based on patient profiles were converted into a risk estimate by using logit-transformation of the total risk score: 1/[1+exp(-
P values < 0.05 were considered statistically significant. Statistical analyses were performed using SAS software version 9.1.3 (SAS Inc., Cary, NC, USA) and R 2.8.1 (The R Foundation for Statistical Computing, Vienna, Austria).
Results
Baseline characteristics of the development and test cohorts
During the study period, 4703 patients with HCC were treated with TACE at our institution; after applying the exclusion criteria, 427 patients (mean age = 59.3 ± 9.7 years) were finally selected for analysis. The baseline clinical and MRI features of the development (n = 211) and the test cohort (n = 216) are summarized in Table 1. The development cohort consisted of 133 (63.0%) TACE-responsive patients and 78 (37.0%) TACE-refractory patients, and the test cohort consisted of 136 (63.0%) TACE-responsive patients and 80 (37.0%) TACE-refractory patients. The prevalence of TACE refractoriness was not significantly different between the two cohorts (P = 0.988). The median levels of bilirubin and albumin were significantly higher in the development cohort (P < 0.001 for both). The prevalence of enhancing capsule on MRI was higher in the test cohort (P = 0.001). Other characteristics showed similar results between the two cohorts.
Baseline characteristics of the study population.
Values are given as n (%), mean ± SD, or median (IQR).
AFP, alpha-fetoprotein; HBP, hepatobiliary phase; HBV, hepatitis B virus; HCV, hepatitis C virus; MRI, magnetic resonance imaging; NHHN, non-hypervascular hypointense nodule; TACE, transarterial chemoembolization.
Development of the predictive scoring system
In the development cohort, univariable analysis showed that AFP ≥ 200 (odds ratio [OR] = 2.375; P = 0.004), maximum tumor size ≥7 cm (OR = 5.099; P < 0.001), tumor number ≥1 (OR = 3.961; P < 0.001), infiltrative morphology (OR = 3.125; P = 0.016), arterial peritumoral enhancement (OR = 3.216; P < 0.001), peritumoral hypointensity on HBP (OR = 5.3; P < 0.001), and the presence of NHHN (OR = 2.527; P = 0.002) on HBP were significantly associated with TACE refractoriness (Table 2).
Clinical features and pre-procedural gadoxetic acid-enhanced MRI findings in the development cohort and the results of univariate analysis for factors associated with TACE refractoriness.
Values are given as n (%) or mean ± SD unless otherwise specified.
AFP, alpha-fetoprotein; CI, confidence interval; HBP, hepatobiliary phase; MRI, magnetic resonance imaging; NHHN, non-hypervascular hypointense nodule; OR, odds ratio; TACE, transarterial chemoembolization.
The following four variables met the bootstrap reliability criterion of ≥50%: maximum tumor size ≥7 cm; tumor number ≥11; peritumoral hypointensity on HBP; and the presence of NHHN on HBP (Supplementary Table 2). A multivariable model for predicting TACE refractoriness was developed with these four variables, and a scoring point was derived on the basis of the adjusted coefficients of the predictors in the final model. Each risk point was rescaled to designate the presence of NHHN as 1 (Table 3). Fig. 2a shows the nomogram for prediction of TACE refractoriness; the total score was in the range of 0–5 points, with a higher score reflecting a higher probability of TACE refractoriness. The highest predicted risk of being TACE refractoriness was 90.2% with a score of 5 if the following predictors were present: maximum tumor size ≥7 cm; tumor number ≥11; peritumoral hypointensity on HBP; and the presence of the NHHN on HBP (Fig. 3).
Predictors selected by multivariate logistic regression model for TACE refractoriness.
*These are shrunk coefficients (by factor 0.935) to correct for over-fitting.
†The point for each predictor was determined on the basis of the β regression coefficient of each predictor.
CI, confidence interval; HBP, hepatobiliary phase; NHHN, non-hypervascular hypointense nodule; OR, odds ratio; TACE, transarterial chemoembolization.

Prediction model for TACE refractoriness. (a) Nomogram for predicting the probability of TACE. The normogram was developed in the development cohort refractoriness based on clinical and pre-treatment MRI findings. Top, predictor points are found on the uppermost point scale that corresponds to each variable. Bottom, points of all variables are added up, and the total point projected on the bottom scale indicates the probability of TACE-refractory HCC. (b–e) ROC curve (b) and calibration curve (c) for TACE refractoriness prediction model in the development cohort show good discrimination with a C-index of 0.796 and adequate fitting and reliability. In the validation cohort, a C-index of ROC analysis (d) was 0.738. Calibration curve (e) yields good agreement between predicted and actual probability of TACE refractoriness. HCC, hepatocellular carcinoma; MRI, magnetic resonance imaging; ROC, receiver operating characteristic; TACE, tranarterial chemoembolization.

A 67-year-old man with three HCCs. (a) Arterial phase images of gadoxetic acid-enhanced MRI show a large HCC (arrows) with a maximal diameter of 8.5 cm (1.5 points) in the right lobe of the liver. (b) In the right hepatic dome, a 1-cm-sized non-hypervascular nodule (arrowhead) is seen adjacent to another 1.2-cm-sized HCC (arrow). (c) Hepatobiliary phase image reveals peritumoral hypointensity (arrowheads, 1.5 points) around the large HCC. (d) On hepatobiliary phase imaging, two non-hypervascular hypointense nodules (arrowheads, 1 point) are present adjacent to the HCC (arrow) in the right hepatic dome. According to the system developed in the present study, the total score was 4 points with an estimated risk for TACE refractoriness of 79.7%. Follow-up CT (e–g) after two sessions of TACE showing residual viable tumor around the lipodol-laden HCC (arrows, e) with rapidly progressing intrahepatic metastases (arrowheads, e), mediastinal metastatic lymph nodes (arrows, f), bone metastasis with extraosseous soft-tissue mass (arrowheads, f) and lung metastasis (arrow, g). CT, computed tomography; HCC, hepatocellular carcinoma; MRI, magnetic resonance imaging; TACE, tranarterial chemoembolization.
Validation of the predictive scoring system
In the internal validation analysis, the final prediction model had a bias-corrected C-index of 0.796 (95% confidence interval [CI] = 0.733–0.860; Fig 2b). The calibration curve for the predicted and actual probabilities of TACE refractoriness yielded good calibration abilities based on the results of the Homer-Lemeshow test (χ2 = 1.031; df = 2; P = 0.597; Fig 2c). The external validation analysis confirmed that the prediction model well-performed in terms of discrimination (C-index = 0.738; 95% CI = 0.670–0.807; Fig. 2d) and calibration (Hosmer-Lemeshow goodness-of-fit test, χ2 = 2.817; df = 2; P = 0.245; Fig 2e).
Discussion
Our results show that gadoxetic acid-enhanced MRI features provide additional information to the conventional staging system in predicting TACE refractoriness in patients with HCC. In addition to the known predictors, including maximum tumor size and tumor number, multivariable analysis showed that peritumoral hypointensity and the presence of NHHN on HBP on gadoxetic acid-enhanced MRI were associated with TACE refractoriness. We developed a prediction model derived from these independent predictors, which showed good performances in the development cohort (C-index = 0.796) and was well-validated in the temporally distinct test cohort (C-index = 0.738).
The maximum tumor size and tumor number are well-validated predictors of early HCC recurrence after TACE (31–33). In line with these studies, we also found that increased tumor burden was significantly related to a higher risk of TACE refractoriness. Due to the heterogeneity of intermediate-stage HCC with respect to tumor burden, previous attempts for further stratification of such patients including the Bolondi subclassification (13) and the Kinki criteria (14) also included tumor size and tumor number as important factors. Our results support these modified systems for selecting a subset of patients with higher risk of TACE refractoriness in whom alternative remedies may be needed.
The present study unveiled the potential of tailored staging using imaging features on gadoxetic acid-enhanced MRI including peritumoral hypointensity and the presence of NHHN on HBP. Among the clinical and imaging features investigated in the study, peritumoral hypointensity on HBP showed the highest OR of 4.187 and was allocated 1.5 points in the scoring system for predicting TACE refractoriness. The clinical implications of peritumoral hypointensity on HBP has been actively studied: recent studies showed that peritumoral hypointensity on HBP is strongly associated with the presence of microvascular invasion in pathologic analysis (20,24), which is a strong predictor of early recurrence after curative hepatectomy or RFA (34,35) as well as TACE-refractory HCC (36). Although several studies showed that the presence of peritumoral hypointensity on HBP can predict early recurrence after curative treatment (23,24,37), only few studies focused on its clinical implications following TACE. Our results suggest that the peritumoral hypointensity on HBP can also be considered as a predictive feature of poor treatment outcomes after TACE. As the use of gadoxetic acid-enhanced MRI in evaluating HCC expands, NHHN on HBP continues to gain attention. NHHN may be considered as a borderline lesion that may transition into HCC (17) and ∼30% of NHHNs eventually transform into hypervascular HCC (38). Thus, the presence of NHHN is considered to indicate enhanced hepatocarcinogenic potential in background liver (29). The concurrence of NHHN on HBP is a significant risk factor for the recurrence of HCC after hepatectomy and RFA, especially for multicentric recurrences (23,29,39). Our results show that NHHN is also a prognostic factor for poor outcomes after TACE. We assume that NHHN has low responsiveness to TACE due to its non-hypervascular nature and the high risk of transformation into overt HCC after TACE.
By using the aforementioned tumor-related factors and imaging features from gadoxetic acid-enhanced MRI, we developed a scoring system that estimates the risk for TACE refractoriness with good discrimination and calibration, and validated its performance on a separate test cohort. Each total score presents the probability of TACE refractoriness and supposedly provides reliable guidance for predicting TACE refractoriness. Gadoxetic acid-enhanced MRI is expected to improve the accuracy of conventional staging systems based on tumor size and number by adding new information on tumor biology. In line with our results, a recent study demonstrated that additional gadoxetic acid-enhanced MRI before HCC treatment was associated with better survival by assisting patient selection (40). Our new staging system for intermediate-stage HCC may be useful for improving patient selection as well by identifying those at risk for TACE refractoriness.
The present study has some limitations. First, there may have been a potential selection bias due to the retrospective nature of this study. Approximately half of the initial eligible patients were excluded from the study as they did not undergo pre-procedural gadoxetic acid-enhanced MRI. This reflects the current practice pattern of pre-procedural imaging in which gadoxetic acid-enhanced MRI is more frequently performed before curative treatments than before TACE. Our results suggest that gadoxetic acid-enhanced MRI may be encouraged before determining treatment options, not only due to its excellent diagnostic performance but also its ability for better patient selection and prognostic implications. Second, although we used robust statistical methods to enhance the reliability of our results, the present study was performed using data from a single center. Therefore, our results should be further validated through prospective studies in multiple institutions.
In conclusion, we developed a simple and easily applicable scoring system based on gadoxetic acid-enhanced pre-procedural MRI for estimating the risk of TACE refractoriness in patients with HCC. By combining conventional staging information with detailed information from gadoxetic acid-enhanced MRI, HCC patients at risk for TACE refractoriness may be identified before TACE and allocated to optimal treatment options.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
