Cholecystectomy is associated with higher risk of recurrence after microwave ablation of hepatocellular carcinoma: a propensity score matching analysis ======================================================================================================================================================= * Hongcai Yang * Yi Yang * Jianping Dou * Rui Cui * Zhigang Cheng * Zhiyu Han * Fangyi Liu * Xiaoling Yu * Xiang Zhou * Jie Yu * Ping Liang ## Abstract **Objective:** To explore the association between cholecystectomy and the prognostic outcomes of patients with hepatocellular carcinoma (HCC) who underwent microwave ablation (MWA). **Methods:** Patients with HCC (*n* = 921) who underwent MWA were included and divided into cholecystectomy (*n* = 114) and non-cholecystectomy groups (*n* = 807). After propensity score matching (PSM) at a 1:2 ratio, overall survival (OS) and disease-free survival (DFS) rates were analyzed to compare prognostic outcomes between the cholecystectomy (*n* = 114) and non-cholecystectomy groups (*n* = 228). Univariate and multivariate Cox analyses were performed to assess potential risk factors for OS and DFS. Major complications were also compared between the groups. **Results:** After matching, no significant differences between groups were observed in baseline characteristics. The 1-, 3-, and 5-year OS rates were 96.5%, 82.1%, and 67.1% in the cholecystectomy group, and 97.4%, 85.2%, and 74.4% in the non-cholecystectomy group (*P* = 0.396); the 1-, 3-, and 5-year DFS rates were 58.4%, 34.5%, and 26.6% in the cholecystectomy group, and 73.6%, 44.7%, and 32.2% in the non-cholecystectomy group (*P* = 0.026), respectively. The intrahepatic distant recurrence rate in the cholecystectomy group was significantly higher than that in the non-cholecystectomy group (*P* = 0.026), and the local tumor recurrence and extrahepatic recurrence rates did not significantly differ between the groups (*P* = 0.609 and *P* = 0.879). Multivariate analysis revealed that cholecystectomy (HR = 1.364, 95% CI 1.023–1.819, *P* = 0.035), number of tumors (2 *vs.* 1: HR = 2.744, 95% CI 1.925–3.912, *P* < 0.001; 3 *vs.* 1: HR = 3.411, 95% CI 2.021–5.759, *P* < 0.001), and γ-GT levels (HR = 1.003, 95% CI 1.000–1.006, *P* < 0.024) were independent risk factors for DFS. The best γ-GT level cut-off value for predicting median DFS was 39.6 U/L (area under the curve = 0.600, *P* < 0.05). A positive correlation was observed between cholecystectomy and γ-GT level (*r* = 0.108, 95% CI −0.001–0.214, *P* = 0.047). Subgroup analysis showed that the DFS rates were significantly higher in the non-cholecystectomy group than the cholecystectomy group when γ-GT ≥39.6 U/L (*P* = 0.044). The 5-, 10-, 15-, 20-, and 25-year recurrence rates from the time of cholecystectomy were 2.63%, 21.93%, 42.11%, 58.77%, and 65.79%, respectively. A significant positive correlation was observed between cholecystectomy and the time from cholecystectomy to recurrence (*r* = 0.205, 95% CI 0.016–0.379, *P* = 0.029). There were no significant differences in complications between groups (*P* = 0.685). **Conclusions:** Patients with HCC who underwent cholecystectomy were more likely to develop intrahepatic distant recurrence after MWA, an outcome probably associated with increased γ-GT levels. Moreover, the recurrence rates increased with time. * Cholecystectomy * microwave ablation * hepatocellular carcinoma * propensity score matching ## Introduction Hepatocellular carcinoma (HCC) is one of the most common malignancies in China, with 364,800 new cases and 318,800 death cases in 2015, which accounted for more than 50% of the newly diagnosed cases and deaths worldwide1,2. Surgical resection and orthotopic liver transplantation, which are accepted as radical therapies for HCC, are limited by cirrhosis, poor hepatic reserves, advanced tumor progression, and the scarcity of donated organs3. Minimally invasive thermal ablation of HCC has become common since the advent of modern imaging and is recommended by the American Association for the Study of Liver Diseases (AASLD), European Association for the Study of the Liver (EASL), Asian Pacific Association for the Study of the Liver (APASL), and National Comprehensive Cancer Network (NCCN) guidelines as an effective therapy for HCC4–7. Microwave ablation (MWA) has been established as a frequently used and acceptable therapeutic method for HCC treatment, according to the Milan criteria, and has been associated with many theoretical advantages, such as higher intra-tumoral temperature, greater ablation volume, and shorter operation time than those of radiofrequency ablation (RFA). Cholecystectomy has been used as a therapeutic surgery to treat gallbladder-related diseases for more than a century, and it is thought to be associated with the declining morbidity rates for gallbladder cancer8. Cholecystectomy can also be conducted simultaneously with HCC curative resection, because the tumor is adjacent to the gallbladder. Typical clinical manifestations followed by cholecystectomy include dilation of the common bile duct, and increased bile duct pressure and cholestasis, which may cause chronic infections in peripheral liver tissue and even systemic inflammation9–12. Chronic inflammation frequently induces necrosis and hepatocyte regeneration, thus resulting in DNA mutations and hepatocyte cancerization13,14. Recently, many epidemiological investigations and systematic reviews have addressed the relationship between cholecystectomy and HCC, and most of these studies have indicated that cholecystectomy is associated with a high risk of HCC15–17. Moreover, one study by Li et al.18 has found that cholecystectomy is associated with a higher risk of early recurrence after curative resection for early-stage HCC. The relationship between cholecystectomy and the prognostic outcomes of patients with HCC after thermal ablation has not been clarified. Therefore, we conducted a retrospective study with long-term follow-up of consecutive patients with HCC who underwent MWA, by using propensity score matching (PSM) to strengthen causal arguments by reducing potential bias. To our knowledge, this is the first study to consider cholecystectomy as a potential risk factor in analyzing the prognostic outcomes of patients with HCC who underwent thermal ablation. ## Materials and methods ### Patients This retrospective study was approved by the Institutional Review Board of the Chinese PLA General Hospital and was exempt from requirements to obtain informed consent. Between January 2005 and December 2017, 1,023 consecutive patients with HCC underwent MWA. Of these, 921 patients (754 men and 167 women) were included in this study after meeting the following inclusion criteria: (a) refusal of, or intolerance to, surgical therapy (e.g., cirrhosis or insufficient liver function); (b) a solitary HCC nodule of 5 cm or less, or as many as 3 nodules 3 cm or less, no extrahepatic metastasis or major vessel invasion (Milan criteria); (c) Child-Pugh class of 0 or A; (d) Eastern Cooperative Oncology Group (ECOG) performance status of 0 or 1; and (e) more than 12 months of follow-up time. The baseline variables were selected according to clinical relevance and the results of previous studies19,20, and included age, sex, maximum tumor diameter, number of tumors, cirrhosis type, Child-Pugh class, a-fetoprotein (AFP) concentration, aspartic transaminase (AST) level, alanine aminotransferase (ALT) level, γ-glutamyl transpeptidase (γ-GT) level, albumin (ALB) level, total bilirubin (TBil) level, creatinine (Cre) level, blood platelet (PLT) count, international normalized ratio (INR), neutrophil count, and lymphocyte count. The upper limits of the normal values in our hospital were used as the cut-off values. Preoperative serum tests were performed 1 week before MWA. ### MWA procedure All MWA procedures were performed by experienced doctors (P.L. and J.Y.) in our department. Color Doppler and grayscale ultrasound were performed to choose the safest access for the antenna. The antenna was connected to the MW generator and inserted into the target area of the tumor. When there were 2 antennas, they were usually deployed in parallel 1–2.5 cm apart. In addition, the antennas were usually deployed parallel to vessels to decrease the damage to large vessels. In general, the microwave energy application was set to 50–60 W for 5–10 min per session. The region of ablation was monitored with real-time ultrasound. MWA emission was stopped when the hyperechoic zone covered the entire tumor, including a safety margin. If necessary, owing to tumor size, multiple overlapping ablations were usually needed to envelope the entire tumor with at least a 5 mm safety margin. ### Follow-up Contrast-enhanced (CE) imaging (CE-magnetic resonance imaging, CE-computed tomography, or CE ultrasound) was performed within 3 days after the ablation to assess whether complete ablation was achieved. After a residual tumor was recognized, additional ablation was performed to achieve complete coagulation. Patients were followed up at 1 month, every 3 months in the first year after ablation, and then every 6 months in the following years. Patients were also monitored prospectively for recurrence with a standard protocol that included serum AFP and CE imaging. An MR scan of the abdomen was performed every 6 months. Recurrent tumors identified during follow-up were treated with optimal treatment, such as MWA, transcatheter arterial chemoembolization, radiation therapy, targeted therapy, or liver transplantation, according to the recommendations of a multidisciplinary tumor board, on the basis of liver function, patient performance, and recurrent tumor characteristics. ### Definition HCC diagnosis was confirmed according to the pathologic examination findings or the current practice guidelines of the AASLD21,22. Standardized terminology and reporting criteria were used for imaging-guided tumor ablation23. Overall survival (OS) was recorded from the local hepatic therapy to the date of death or the date of the last follow-up. Disease-free survival (DFS) was calculated from the beginning of therapy to the time of tumor recurrence or to the time of death or last patient visit. The recurrence of tumors was divided into 3 types: local tumor recurrence, intrahepatic distant recurrence, and extrahepatic recurrence. Local tumor recurrence was defined as the appearance of new tumor foci at the ablative margin after the local eradication of all tumor cells with ablation. Intrahepatic distant recurrence was defined as the detection of the tumor distally from the site of ablation within the liver. Extrahepatic recurrence was defined as the appearance of tumors in extrahepatic organs. Ultrasound-guided fine-needle biopsy was required to confirm the diagnosis when imaging findings were atypical. ### Statistical analysis All statistical analyses were performed in SPSS Version 22 (College Station, TX, USA). The *t*-test or Mann-Whitney U test was used to assess continuous variables including the maximum tumor diameter, AST level, ALT level, γ-GT level, ALB level, TBil level, Cre level, PLT count, INR, neutrophil count, and lymphocyte count. Fisher’s exact test and χ2-test were used to assess the categorical variables including age, sex, number of tumors, AFP level, type of cirrhosis, and Child-Pugh class. Binary logistic regression was performed to achieve 1:2 matching between the cholecystectomy group and the non-cholecystectomy group. A 1:2 nearest-neighbor algorithm was used with a maximum tolerated difference between propensity scores in the range of 0.05. The Kaplan-Meier method was used to evaluate cumulative incidence rates of OS and DFS, and the log-rank test was used to evaluate differences between the cholecystectomy group and the non-cholecystectomy group. Univariable and multivariable Cox proportional hazard models were used to assess the association between possible risk factors for OS and DFS between the non-cholecystectomy group and the cholecystectomy group. The best cut-off value of the γ-GT level was calculated with time-dependent receiver operating characteristic (ROC) curve analysis. Spearman’s correlation analysis was used to assess potential correlations between risk factors. A two-tailed 95% confidence interval (CI) was applied to reveal the accuracy of the hazard ratio (HR), and a *P* value less than 0.05 was considered to indicate a statistically significant difference. ## Results ### Baseline characteristics A total of 1,023 patients received MWA for HCC between January 2005 and December 2017 in our department; among them, 102 patients were excluded. A total of 921 patients were included in this study, of whom 114 (12.4%) were in the cholecystectomy group, and 807 (87.2%) were in the non-cholecystectomy group. The reasons for cholecystectomy were as follows: cholelithiasis in 72 patients, gallbladder polyps in 26 patients, acute cholecystitis in 9 patients, and gallbladder adenomyoma in 7 patients. Before PSM, the maximum tumor diameter, number of tumors, AST level, ALT level, and lymphocyte count were significantly different between the 2 groups (**Table 1**). In view of nonrandomized study design and confounding variables, we used PSM for the analysis of observational data to reduce potential bias and performed PSM at a 1:2 ratio. All 114 patients in the cholecystectomy group were successfully matched with 228 patients in the non-cholecystectomy group. No statistically significant differences were observed in the baseline characteristics between the groups after PSM. A flow diagram of the patient selection process is shown in **Figure 1**. The baseline characteristics of patients before and after matching are listed in **Table 1**. ![Figure 1](http://www.cancerbiomed.org/https://www.cancerbiomed.org/content/cbm/17/2/478/F1.medium.gif) [Figure 1](http://www.cancerbiomed.org/content/17/2/478/F1) Figure 1 Flow diagram of the patient selection process. View this table: [Table 1](http://www.cancerbiomed.org/content/17/2/478/T1) Table 1 Baseline characteristics of patients before and after matching ### Comparison of survival The median follow-up times was 45.8 months (range, 0.6–96.0 months) in the cholecystectomy group and 41.3 months (range, 1.3–96.0 months) in the non-cholecystectomy group (*P* = 0.697). Until the follow-up ended, the death incidence rates were 28.95% and 23.25% for the cholecystectomy group and non-cholecystectomy group, respectively. The cumulative 1-, 3-, and 5-year OS rates were 96.5%, 82.1%, and 67.1% in the cholecystectomy group, and 97.4%, 85.2%, and 74.4% in the non-cholecystectomy group, respectively (*P* = 0.396, **Figure 2A**). ![Figure 2](http://www.cancerbiomed.org/https://www.cancerbiomed.org/content/cbm/17/2/478/F2.medium.gif) [Figure 2](http://www.cancerbiomed.org/content/17/2/478/F2) Figure 2 Cumulative OS rates (A) and DFS rates (B) between the non-cholecystectomy group and the cholecystectomy group. The median DFS time was 16.7 months (range, 0.6–96.0 months) in the cholecystectomy group and 21.4 months (range, 1.3–94.6 months) in the non-cholecystectomy group (*P* = 0.358). Until the follow-up ended, the recurrence rates were 70.2% and 58.8% for the cholecystectomy group and non-cholecystectomy group, respectively. The cumulative 1-, 3-, and 5-year DFS rates were 58.4%, 34.5%, and 26.6% in the cholecystectomy group, and 73.6%, 44.7%, and 32.2% in the non-cholecystectomy group, respectively. The cumulative DFS rates of the non-cholecystectomy group were significantly higher than those of the cholecystectomy group (*P* = 0.026, **Figure 2B**). ### Cox analysis of risk factors associated with OS and DFS All baseline characteristics of patients with HCC and cholecystectomy were included in univariate and multivariate Cox analyses of potential risk factors associated with survival. For OS, univariate analysis revealed that age (HR = 1.678, 95% CI 1.089–2.585, *P* = 0.010), number of tumors (2 *vs.* 1: HR = 2.261, 95% CI 1.426–3.587, *P* < 0.001; 3 *vs.* 1: HR = 1.678, 95% CI 1.089–2.585, *P* = 0.003), Child-Pugh class (HR = 2.318, 95% CI 1.158–4.637,* P* = 0.018), and lymphocyte count (HR = 0.074, 95% CI 0.009–0.625,* P* = 0.017) were associated risk factors (**Table 2**). Multivariate analysis revealed that age (HR = 1.778, 95% CI 1.130–2.829, *P* = 0.013) and number of tumors (2 *vs.* 1: HR = 3.010, 95% CI 1.658–5.463, *P* < 0.001; 3 *vs.* 1: HR = 5.981, 95% CI 2.556–13.997, *P* < 0.001) were independent risk factors associated with OS (**Table 2**). View this table: [Table 2](http://www.cancerbiomed.org/content/17/2/478/T2) Table 2 Cox analyses of risk factors associated with OS For DFS, univariate analysis revealed that cholecystectomy (HR = 1.367, 95% CI 1.036–1.804, *P* = 0.027), sex (HR = 1.486, 95% CI 1.015–2.176, *P* = 0.041), number of tumors (2 *vs.* 1: HR = 2.332, 95% CI 1.745–3.118, *P* < 0.001; 3 *vs.* 1: HR = 2.522, 95% CI 1.647–3.862, *P* < 0.001), γ-GT level (HR = 1.004, 95% CI 1.002–1.006, *P* < 0.001), and Cre level (HR = 1.002, 95% CI 1.000–1.004, *P* = 0.012) were associated risk factors (**Table 3**). Multivariate analysis revealed that cholecystectomy (HR = 1.364, 95% CI 1.023–1.819, *P* = 0.035), number of tumors (2 *vs.* 1: HR = 2.744, 95% CI 1.925–3.912, *P* < 0.001; 3 *vs.* 1: HR = 3.411, 95% CI 2.021–5.759, *P* < 0.001), and γ-GT level (HR = 1.003, 95% CI 1.000–1.006, *P* = 0.024) were independent risk factors associated with DFS (**Table 3**). View this table: [Table 3](http://www.cancerbiomed.org/content/17/2/478/T3) Table 3 Cox analyses of risk factors associated with DFS ### Subgroup analysis according to the risk factors The cumulative OS rates of the group < 60 years of age were significantly higher than those of the group ≥ 60 years of age (*P* = 0.018, **Figure 3A**). The cumulative OS rates of the 1-tumor group were significantly higher than those of the 2-tumor group (*P* < 0.001, **Figure 3B**) and the 3-tumor group (*P* = 0.003, **Figure 3B**), but there was no statistically significant difference between the 2-tumor group and the 3-tumor group (*P* = 0.612, **Figure 3B**). ![Figure 3](http://www.cancerbiomed.org/https://www.cancerbiomed.org/content/cbm/17/2/478/F3.medium.gif) [Figure 3](http://www.cancerbiomed.org/content/17/2/478/F3) Figure 3 Cumulative OS rates in groups < 60 years and ≥ 60 years of age (A), and in the 1-tumor, 2-tumor, and 3-tumor groups (B). Furthermore, the cumulative DFS rates of the 1-tumor group were significantly higher than those of the 2-tumor group (*P* < 0.001, **Figure 4A**) and 3-tumor group (*P* < 0.001, **Figure 4A**), but there was no statistically significant difference between the 2-tumor group and the 3-tumor group (*P* = 0.716, **Figure 4A**). The cumulative DFS rates of the γ-GT < 50 U/L group were significantly higher than those of the γ-GT ≥ 50 U/L group (*P* = 0.006, **Figure 4B**). ![Figure 4](http://www.cancerbiomed.org/https://www.cancerbiomed.org/content/cbm/17/2/478/F4.medium.gif) [Figure 4](http://www.cancerbiomed.org/content/17/2/478/F4) Figure 4 Cumulative DFS rates in the γ-GT < 50 U/L and γ-GT ≥ 50 U/L groups (A) and in the 1-tumor, 2-tumor, and 3-tumor groups (B). In the subgroup of patients with 1 tumor, the 1-, 3-, and 5-year DFS rates were 63.8%, 46.9%, and 34.8% in the cholecystectomy group, and 84.1%, 61.9%, and 45.0% in the non-cholecystectomy group, respectively (*P* = 0.012, **Figure 5A**). In the subgroup of patients with 2 or 3 tumors, no statistically significant differences in DFS rates were observed between the groups (*P* = 0.723 and *P* = 0.783, respectively, **Figure 5A**). ![Figure 5](http://www.cancerbiomed.org/https://www.cancerbiomed.org/content/cbm/17/2/478/F5.medium.gif) [Figure 5](http://www.cancerbiomed.org/content/17/2/478/F5) Figure 5 Subgroup analysis according to the number of tumors (A) and γ-GT levels (B) in the non-cholecystectomy group and the cholecystectomy group. When the cut-off value of γ-GT was 50 U/L (the upper limit of normal values in our hospital), subgroup analysis revealed no statistically significant differences in DFS rates between the cholecystectomy group and the non-cholecystectomy group, not only for patients with γ-GT < 50 U/L (*P* = 0.140, **Figure 5B**) but also for patients with γ-GT ≥ 50 U/L (*P* = 0.084, **Figure 5B**). Given that the γ-GT level was an independent risk factor associated with DFS, time-dependent ROC curve analysis indicated that 39.6 U/L was the best cut-off value of the γ-GT level for predicting median DFS at 20.2 months (area under the curve = 0.600, *P* < 0.05) (**Figure 6A**). Then 342 patients were divided into 2 groups: a low γ-GT group (< 39.6 U/L, *n* = 161, 47.1%) and a high γ-GT group (≥ 39.6 U/L, *n* = 181, 52.9%). The DFS rates of the low γ-GT group were significantly higher than those of the high γ-GT group (*P* < 0.001, **Figure 6B**). In the subgroup of patients with γ-GT ≥ 39.6 U/L, the DFS rates of the non-cholecystectomy group were significantly higher than those of the cholecystectomy group (*P* = 0.044, **Figure 6B**); however, in the subgroup of patients with γ-GT < 39.6 U/L, there were no statistically significant differences in the DFS rates between the groups (*P* = 0.465, **Figure 6B**). ![Figure 6](http://www.cancerbiomed.org/https://www.cancerbiomed.org/content/cbm/17/2/478/F6.medium.gif) [Figure 6](http://www.cancerbiomed.org/content/17/2/478/F6) Figure 6 Time-dependent ROC curve for the cut-off value of the γ-GT level (A); subgroup analysis according to the γ-GT level (B) of the non-cholecystectomy and cholecystectomy groups. ### Comparison of recurrence As mentioned above, the total recurrence rates were 70.2% (80/114) in the cholecystectomy group and 58.8% (134/228) in the non-cholecystectomy group, and a higher recurrence rate was observed in the cholecystectomy group (*P* = 0.040); the local tumor recurrence rates were 8.8% (10/114) and 10.5% (24/228), respectively, and there was no statistically significant difference between the groups (*P* = 0.609). The intrahepatic distant recurrence rates were 54.4% (62/114) and 41.7% (95/228) in the cholecystectomy group and non-cholecystectomy group, respectively, and higher recurrence rates were observed in the cholecystectomy group (*P* = 0.026). The extrahepatic recurrence rates were 7.0% (8/80) and 4.4% (15/134) in the cholecystectomy group and non-cholecystectomy group, respectively, and no statistically significant difference was observed between the groups (*P* = 0.879) (**Figure 7**). ![Figure 7](http://www.cancerbiomed.org/https://www.cancerbiomed.org/content/cbm/17/2/478/F7.medium.gif) [Figure 7](http://www.cancerbiomed.org/content/17/2/478/F7) Figure 7 Forest plot for the comparison of recurrence. ### Analysis of correlation between cholecystectomy and ***γ***-GT level As mentioned above, both cholecystectomy and γ-GT levels were independent risk factors for DFS and were associated with bile metabolism and liver inflammation. Therefore, we performed Spearman’s correlation analysis and found a significantly positive correlation between cholecystectomy and γ-GT levels (*r* = 0.108, 95% CI −0.001–0.214, *P* = 0.047, **Figure 8**). ![Figure 8](http://www.cancerbiomed.org/https://www.cancerbiomed.org/content/cbm/17/2/478/F8.medium.gif) [Figure 8](http://www.cancerbiomed.org/content/17/2/478/F8) Figure 8 Correlation between cholecystectomy and γ-GT level. ### Analysis of time from cholecystectomy to recurrence In the cholecystectomy group (*n* = 114), patients were divided into a recurrence subgroup (*n* = 80) and non-recurrence subgroup (*n* = 34). The time from cholecystectomy to recurrence was 172.8 ± 84.5 months (median 144.6 months; range 31.5–307.0 months) in the recurrence subgroup, and the mean time from cholecystectomy to death or the ending of follow-up was 132.5 ± 72.6 months (median 138.6 months; range 29.2–288.6 months) in the non-recurrence subgroup; a statistically significant difference was observed between the groups (*P* = 0.017, **Figure 9A**). The 5-, 10-, 15-, 20-, and 25-year recurrence rates from the time of cholecystectomy were 2.6%, 21.9%, 42.1%, 58.8%, and 65.8%, respectively (**Figure 9B**). By Spearman’s correlation analysis, we found a significantly positive correlation between recurrence and the time from cholecystectomy to recurrence (*r* = 0.205, 95% CI 0.016–0.379, *P* = 0.029, **Figure 9C**). ![Figure 9](http://www.cancerbiomed.org/https://www.cancerbiomed.org/content/cbm/17/2/478/F9.medium.gif) [Figure 9](http://www.cancerbiomed.org/content/17/2/478/F9) Figure 9 T test for time from cholecystectomy to recurrence in recurrence group and time from cholecystectomy to death or the ending of follow-up in non-recurrence group (A); proportion of patient recurrence from the time of cholecystectomy (B); correlation between recurrence and the time from cholecystectomy to recurrence (C). ### Complications The complication rates associated with MWA were 3.3% (30/921) in the overall data (11 in the cholecystectomy group; 19 in the non-cholecystectomy group, *P* = 0.685). Complications that may be associated with cholecystectomy (such as liver abscess and bile leakage) were not significantly more frequent in the cholecystectomy group than in the non-cholecystectomy group (1.8% in the cholecystectomy group *vs.* 1.3% in the non-cholecystectomy group, *P* = 0.750). ## Discussion HCC is an inflammation-related cancer in which more than 90% of cases have a background of hepatic injury or chronic inflammation, such as hepatitis B or C virus infection, cirrhosis, cholestasis, and cholecystitis24–26. Cholecystectomy-induced cholestasis is the major cause of the intrahepatic chronic inflammatory environment, which is considered to increase the risk and recurrence of HCC24. However, whether cholecystectomy influences the prognostic outcomes of patients with HCC undergoing thermal ablation has not yet been clarified. Cholecystectomy was found to be an independent risk factor associated with DFS in this study (*P* = 0.035), and the recurrence rates in the cholecystectomy group were significantly higher than those in the non-cholecystectomy group (*P* = 0.034). High recurrence rate is considered the primary reason for the dismal oncological outcomes in HCC27. Some studies28–30 have demonstrated that removal of the gallbladder is always accompanied by dysfunction in the sphincter of Oddi, thus potentially leading to dilation of the common bile duct, elevated bile duct pressure, subsequent cholestasis, and chronic inflammation. The intrahepatic chronic inflammatory microenvironment is the most critical factor associated with HCC, and inflammation-induced neoangiogenesis, inflammation-related cytokines, chemokines, and growth factors all play critical roles in cancer recurrence and metastasis. Moreover, the mean time from cholecystectomy to recurrence was 172.8 ± 84.5 months, and a statistically significant difference was observed in this group compared with time from cholecystectomy to death or the ending of follow-up in non-recurrence group (*P* = 0.017). The 5-, 10-, 15-, 20-, and 25-year recurrence rates from the time of cholecystectomy were 2.6%, 21.9%, 42.1%, 58.8%, and 65.8%, respectively, and significantly positive correlations were observed between recurrence and the time from cholecystectomy to recurrence (*r* = 0.205, 95% CI 0.016–0.379, *P* = 0.029), thus indicating an increasing trend of recurrence over time. Lagergren et al.16 have reported a significantly greater overall risk of HCC in patients who have undergone cholecystectomy than in the corresponding background population (standardized incidence ratio = 1.24, 95% CI 1.11–1.38), and the risk gradually increases with longer follow-up time (*P* for trend = 0.003). Our findings revealed that the chronic inflammation introduced by cholecystectomy plays an important role in the carcinogenic process. In addition to cholecystectomy, the γ-GT level was found to be an independent risk factor associated with DFS (*P* = 0.024). Elevated γ-GT levels result from both alcoholic and non-alcoholic fatty liver disease, cholestatic liver disease, and induction by drugs such as phenytoin. γ-GT has been widely applied for detection and diagnosis of liver disease, which is a nearly ubiquitous epithelial enzyme that is responsible for the catabolism of extracellular glutathione. In addition to its association with liver disease, γ-GT has been associated with high all-cause mortality, cardiovascular disease incidence and death, diabetes, and cancer incidence and death31–33. Some studies34–36 have reported that γ-GT levels have significant prognostic value as a risk factor for HCC, in agreement with our findings. The synthesis of γ-GT is excessive in patients with HCC; moreover, cholecystectomy-related cholestasis may block the excretion of γ-GT, and the blocked γ-GT then returns to the blood and is accompanied by bile, thus increasing serum γ-GT levels. Our findings revealed a positive correlation between cholecystectomy and γ-GT levels (*r* = 0.108, 95% CI −0.001–0.214, *P* = 0.047). The above results suggested that patients with HCC who underwent cholecystectomy and had high γ-GT levels were prone to recurrence after thermal ablation. In this study, the median follow-up time in the cholecystectomy group and non-cholecystectomy group was 45.8 and 41.3 months, respectively (*P* = 0.697). However, no statistically significant difference in OS rates was observed between the cholecystectomy group and the non-cholecystectomy group (*P* = 0.396). The reason for this finding may be that most patients with recurrence probably benefitted from our multimodal therapy and timely comprehensive management. A phase III randomized controlled trial examining MWA for patients with HCC37 has reported 1-, 3-, and 5-year OS rates of 96.4%, 81.9%, and 67.3%, respectively, findings similar to our oncological outcomes. Several limitations in our study must be noted. First, our retrospective study design may potentially have introduced selection bias. Although we attempted to simulate randomization by using PSM analysis, there remains a possibility of uncontrolled confounding factors. Second, because this was a single-center study, and there was no universal consensus of MWA, our results might not be consistently reproducible in other settings. Despite these limitations, to our knowledge, no previous study has addressed a possible relationship between cholecystectomy and increased risk of recurrence in patients with HCC after thermal ablation. Our findings may help validate the association between cholecystectomy and the prognostic outcomes of patients with HCC who underwent thermal ablation. In summary, cholecystectomy is an independent prognostic factor associated with the DFS of patients with HCC who underwent thermal ablation. Though it is too early to consider any potential clinical recommendations, we will continue to focus on patients with HCC who underwent cholecystectomy and perform further medical research to explore the relevant mechanisms. ## Acknowledgements This work was supported by the National Key R&D Program of China (Grant No. 2017YFC0112000), and the National Natural Science Foundation of China (Grant Nos. 81627803 and 91859201). ## Footnotes * *These authors contributed equally to this work. * **Conflict of interest statement** No potential conflicts of interest are disclosed. * Received July 28, 2019. * Accepted October 7, 2019. * Copyright: © 2020, Cancer Biology & Medicine [https://creativecommons.org/licenses/by/4.0/](https://creativecommons.org/licenses/by/4.0/) This is an open access article distributed under the terms of the [Creative Commons Attribution License (CC BY) 4.0](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium, provided the original author and source are credited. ## References 1. 1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019; 69: 7–34. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.3322/caac.21551&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=30620402&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) 2. 2. Zheng R, Qu C, Zhang S, Zeng H, Sun K, Gu X, et al. Liver cancer incidence and mortality in China: temporal trends and projections to 2030. Chin J Cancer Res. 2018; 30(6):571–9. 3. 3. Moris D, Vernadakis S, Papalampros A, Petrou A, Dimitroulis D, Spartalis E, et al. The effect of Guidelines in surgical decision making: the paradigm of hepatocellular carcinoma. J BUON. 2016; 21: 1332–6. 4. 4. Omata M, Lesmana LA, Tateishi R, Chen PJ, Lin SM, Yoshida H, et al. Asian Pacific Association for the Study of the Liver consensus recommendations on hepatocellular carcinoma. Hepatol Int. 2010; 4: 439–74. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1007/s12072-010-9165-7&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=20827404&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) 5. 5. Marrero JA, Kulik LM, Sirlin CB, Zhu AX, Finn RS, Abecassis MM, et al. Diagnosis, staging, and management of hepatocellular carcinoma: 2018 practice guidance by the American Association for the Study of Liver Diseases. Hepatology. 2018; 68: 723–50. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1002/hep.29913&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) 6. 6.European Association for the Study of the Liver, European Organisation for Research and Treatment of Cancer. EASLEORTC clinical practice guidelines: management of hepatocellular carcinoma. J Hepatol. 2012; 56: 908–43. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1016/j.jhep.2011.12.001&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=22424438&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) [Web of Science](http://www.cancerbiomed.org/lookup/external-ref?access_num=000302272600024&link_type=ISI) 7. 7. Benson AB 3rd., D’Angelica MI, Abbott DE, Abrams TA, Alberts SR, Saenz DA, et al. NCCN guidelines insights: hepatobiliary cancers, Version 1.2017. J Natl Compr Canc Netw. 2017; 15: 563–73. [Abstract/FREE Full Text](http://www.cancerbiomed.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiam5jY24iO3M6NToicmVzaWQiO3M6ODoiMTUvNS81NjMiO3M6NDoiYXRvbSI7czoxODoiL2NibS8xNy8yLzQ3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 8. 8. Wood R, Fraser LA, Brewster DH, Garden OJ. Epidemiology of gallbladder cancer and trends in cholecystectomy rates in Scotland, 1968–1998. Eur J Cancer. 2003; 39: 2080–6. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1016/S0959-8049(03)00370-8&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=12957463&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) [Web of Science](http://www.cancerbiomed.org/lookup/external-ref?access_num=000185510400023&link_type=ISI) 9. 9. Li Y, He R, Chen S, Qu Y. Effect of dexmedetomidine on early postoperative cognitive dysfunction and peri-operative inflammation in elderly patients undergoing laparoscopic cholecystectomy. Exp Ther Med. 2015; 10: 1635–42. 10. 10. Csendes A, Kruse A, Funch-Jensen P, Oster MJ, Ornsholt J, Amdrup E. Pressure measurements in the biliary and pancreatic duct systems in controls and in patients with gallstones, previous cholecystectomy, or common bile duct stones. Gastroenterology. 1979; 77: 1203–10. [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=499707&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) [Web of Science](http://www.cancerbiomed.org/lookup/external-ref?access_num=A1979HU67000006&link_type=ISI) 11. 11. Tanaka M, Ikeda S, Nakayama F. Change in bile duct pressure responses after cholecystectomy: loss of gallbladder as a pressure reservoir. Gastroenterology. 1984; 87: 1154–9. [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=6479537&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) 12. 12. Madureira FA, Manso JE, Madureira Filho D, Iglesias AC. Inflammation in laparoendoscopic single-site surgery versus laparoscopic cholecystectomy. Surg Innov. 2014; 21: 263–8. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1177/1553350613499454&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=23945842&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) 13. 13. El-Serag HB, Rudolph KL. Hepatocellular carcinoma: epidemiology and molecular carcinogenesis. Gastroenterology. 2007; 132: 2557–76. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1053/j.gastro.2007.04.061&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=17570226&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) [Web of Science](http://www.cancerbiomed.org/lookup/external-ref?access_num=000247146900032&link_type=ISI) 14. 14. Tarao K, Hoshino H, Shimizu A, Ohkawa S, Nakamura Y, Harada M, et al. Role of increased DNA synthesis activity of hepatocytes in multicentric hepatocarcinogenesis in residual liver of hepatectomized cirrhotic patients with hepatocellular carcinoma. Jpn J Cancer Res. 1994; 85: 1040–4. 15. 15. Nogueira L, Freedman ND, Engels EA, Warren JL, Castro F, Koshiol J. Gallstones, cholecystectomy, and risk of digestive system cancers. Am J Epidemiol. 2014; 179: 731–9. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1093/aje/kwt322&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=24470530&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) 16. 16. Lagergren J, Mattsson F, El-Serag H, Nordenstedt H. Increased risk of hepatocellular carcinoma after cholecystectomy. Br J Cancer. 2011; 105: 154–6. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1038/bjc.2011.181&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=21610710&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) [Web of Science](http://www.cancerbiomed.org/lookup/external-ref?access_num=000292182400023&link_type=ISI) 17. 17. Guo L, Mao J, Li Y, Jiao Z, Guo J, Zhang J, et al. Cholelithiasis, cholecystectomy and risk of hepatocellular carcinoma: a metaanalysis. J Cancer Res Ther. 2014; 10: 834–8. 18. 18. Li T, Wang SK, Zhi XT, Zhou J, Dong ZR, Zhang ZL, et al. Cholecystectomy is associated with higher risk of early recurrence and poorer survival after curative resection for early stage hepatocellular carcinoma. Sci Rep. 2016; 6: 28229. 19. 19. Berber E, Siperstein A. Local recurrence after laparoscopic radiofrequency ablation of liver tumors: an analysis of 1032 tumors. Ann Surg Oncol. 2008; 15: 2757–64. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1245/s10434-008-0043-7&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=18618182&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) [Web of Science](http://www.cancerbiomed.org/lookup/external-ref?access_num=000259186300020&link_type=ISI) 20. 20. Li T, Fan J, Qin LX, Zhou J, Sun HC, Qiu SJ, et al. Risk factors, prognosis, and management of early and late intrahepatic recurrence after resection of primary clear cell carcinoma of the liver. Ann Surg Oncol. 2011; 18: 1955–63. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1245/s10434-010-1540-z&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=21240562&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) 21. 21. Bruix J, Sherman M, American Association for the Study of Liver Disease. Management of hepatocellular carcinoma: an update. Hepatology. 2011; 53: 1020–2. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1002/hep.24199&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=21374666&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) [Web of Science](http://www.cancerbiomed.org/lookup/external-ref?access_num=000288211200032&link_type=ISI) 22. 22. Bruix J, Sherman M, Llovet JM, Beaugrand M, Lencioni R, Burroughs AK, et al. Clinical management of hepatocellular carcinoma. Conclusions of the Barcelona – 2000 EASL conference. European Association for the Study of the Liver. J Hepatol. 2001; 35: 421–30. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1016/S0168-8278(01)00130-1&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=11592607&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) [Web of Science](http://www.cancerbiomed.org/lookup/external-ref?access_num=000171210600017&link_type=ISI) 23. 23. Ahmed M, Solbiati L, Brace CL, Breen DJ, Callstrom MR, Charboneau JW, et al. Image-guided tumor ablation: standardization of terminology and reporting criteria – a 10-year update. Radiology. 2014; 273: 241–60. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1148/radiol.14132958&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=24927329&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) 24. 24. Kao WY, Hwang CY, Su CW, Chang YT, Luo JC, Hou MC, et al. Risk of hepato-biliary cancer after cholecystectomy: a nationwide cohort study. J Gastrointest Surg. 2013; 17: 345–51. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1007/s11605-012-2090-4&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=23188223&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) 25. 25. Triolo M, Della Corte C, Colombo M. Impact of HBV therapy on the incidence of hepatocellular carcinoma. Liver Int. 2014; 34: 139–45. 26. 26. Yu L, Cheng YJ, Cheng ML, Yao YM, Zhang Q, Zhao XK, et al. Quantitative assessment of common genetic variations in HLA-DP with hepatitis B virus infection, clearance and hepatocellular carcinoma development. Sci Rep. 2015; 5: 14933. 27. 27. Yang Y, Si T. Yttrium-90 transarterial radioembolization versus conventional transarterial chemoembolization for patients with hepatocellular carcinoma: a systematic review and meta-analysis. Cancer Biol Med. 2018; 15: 299–310. [Abstract/FREE Full Text](http://www.cancerbiomed.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiY2JtIjtzOjU6InJlc2lkIjtzOjg6IjE1LzMvMjk5IjtzOjQ6ImF0b20iO3M6MTg6Ii9jYm0vMTcvMi80NzguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 28. 28. Wang Q, Blank S, Fiel MI, Kadri H, Luan W, Warren L, et al. The severity of liver fibrosis influences the prognostic value of inflammation-based scores in hepatitis B-associated hepatocellular carcinoma. Ann Surg Oncol. 2015; 22: S1125–32. 29. 29. Aravalli RN, Greten TF. FoxC1: novel regulator of inflammationinduced metastasis in hepatocellular carcinoma. Gastroenterology. 2015; 149: 861–3. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1053/j.gastro.2015.08.032&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=26302490&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) 30. 30. Poon RT, Ng IO, Lau C, Yu WC, Yang ZF, Fan ST, et al. Tumor microvessel density as a predictor of recurrence after resection of hepatocellular carcinoma: a prospective study. J Clin Oncol. 2002; 20: 1775–85. [Abstract/FREE Full Text](http://www.cancerbiomed.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjk6IjIwLzcvMTc3NSI7czo0OiJhdG9tIjtzOjE4OiIvY2JtLzE3LzIvNDc4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 31. 31. Kazemi-Shirazi L, Endler G, Winkler S, Schickbauer T, Wagner O, Marsik C. Gamma glutamyltransferase and long-term survival: is it just the liver? Clin Chem. 2007; 53: 940–6. [Abstract/FREE Full Text](http://www.cancerbiomed.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiY2xpbmNoZW0iO3M6NToicmVzaWQiO3M6ODoiNTMvNS85NDAiO3M6NDoiYXRvbSI7czoxODoiL2NibS8xNy8yLzQ3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 32. 32. Strasak AM, Rapp K, Brant LJ, Hilbe W, Gregory M, Oberaigner W, et al. Association of gamma-glutamyltransferase and risk of cancer incidence in men: a prospective study. Cancer Res. 2008; 68: 3970–7. [Abstract/FREE Full Text](http://www.cancerbiomed.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjEwOiI2OC8xMC8zOTcwIjtzOjQ6ImF0b20iO3M6MTg6Ii9jYm0vMTcvMi80NzguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 33. 33. Sakuta H, Suzuki T, Yasuda H, Ito T. Gamma-glutamyl transferase and metabolic risk factors for cardiovascular disease. Intern Med. 2005; 44: 538–41. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.2169/internalmedicine.44.538&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=16020876&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) [Web of Science](http://www.cancerbiomed.org/lookup/external-ref?access_num=000230096900006&link_type=ISI) 34. 34. Ruhl CE, Everhart JE. Elevated serum alanine aminotransferase and gamma-glutamyltransferase and mortality in the United States population. Gastroenterology. 2009; 136: 477–85.e411. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.1053/j.gastro.2008.10.052&link_type=DOI) [PubMed](http://www.cancerbiomed.org/lookup/external-ref?access_num=19100265&link_type=MED&atom=%2Fcbm%2F17%2F2%2F478.atom) [Web of Science](http://www.cancerbiomed.org/lookup/external-ref?access_num=000262967700019&link_type=ISI) 35. 35. Le Y, Shen JX, Zhang YF, He MK, Kan A, Chen HL, et al. Transarterial chemoembolization related to good survival for selected patients with advanced hepatocellular carcinoma. J Cancer. 2019; 10: 665–71. 36. 36. Xu XS, Wan Y, Song SD, Chen W, Miao RC, Zhou YY, et al. Model based on gamma-glutamyltransferase and alkaline phosphatase for hepatocellular carcinoma prognosis. World J Gastroenterol. 2014; 20: 10944–52. [CrossRef](http://www.cancerbiomed.org/lookup/external-ref?access_num=10.3748/wjg.v20.i31.10944&link_type=DOI) 37. 37. Yu J, Yu XL, Han ZY, Cheng ZG, Liu FY, Zhai HY, et al. Percutaneous cooled-probe microwave versus radiofrequency ablation in early-stage hepatocellular carcinoma: a phase III randomised controlled trial. Gut. 2017; 66: 1172–3. [FREE Full Text](http://www.cancerbiomed.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjk6IjY2LzYvMTE3MiI7czo0OiJhdG9tIjtzOjE4OiIvY2JtLzE3LzIvNDc4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==)