1.Study characteristics
A total of 1430 relevant papers were retrieved from PubMed, SCI and Cochrane databases after removing duplicates, and another 21 papers were obtained from other relevant references. After reading the titles and abstracts, 59 papers were initially selected to meet the criteria. After carefully reading the literature, 25 papers were selected and the basic information and underlying data were extracted and presented (Table 1), and the specific screening process has been shown (Fig. 1), which included 5 RCT studies and 20 retrospective cohort studies, all were two-arm studies. These studies were published from 2005 to 2020, including 4025 patients.
Table 1
Author, year
|
Country
Region
|
Participants, n
|
Mean age
(years)
|
OS rates (%)
|
|
DFS rates (%)
|
|
|
|
|
|
1 year
|
3 years
|
1 year
|
3 years
|
Ng,2017[4]
|
Hong Kong
|
RFA:109
LR:109
|
57.0 (23.0–78.0)
55.0 (31.0–82.0)
|
95.4
94.5
|
82.3
80.6
|
70.6
74.1
|
46.6
50.9
|
Lee,2018[5]
|
Korea
|
RFA:34
LR:29
|
56.1 ± 7.4
55.6 ± 7.9
|
100
96.6
|
97.1
96.6
|
-
-
|
44.1
66.7
|
Chong,2020[12]
|
Hong Kong
|
RFA:46
MWA:47
|
64.5 (42.0–85.0)
63.0 (50.0–80.0)
|
93.5
97.9
|
72.7
67.1
|
58.7
51.1
|
22.7
24.1
|
Feng, 2012[20]
|
Chinese
mainland
|
RFA:84
LR:84
|
51.0 (24.0–83.0)
47.0 (18.0–76.0)
|
96
93.1
|
74.8
67.2
|
90.6
86.2
|
61.1
49.6
|
Wang,2015[21]
|
Chinese
mainland
|
RFA:180
CA:180
|
53.34 ± 8.91
53.87 ± 9.59
|
97
97
|
66
67
|
84
89
|
50
54
|
Ohmoto,2008[22]
|
Japan
|
RFA:34
MWA:49
|
67.0 (44.0–78.0)
64.0 (38.0–75.0)
|
100
89
|
70
49
|
63
52
|
26
15
|
Theodora,2016[23]
|
US
|
RFA:55
MWA:99
|
62.0 (23.0–88.0)
61.0 (44.0–82.0)
|
83.6
78.8
|
49.1
24
|
-
-
|
-
-
|
Lu, 2005[24]
|
Chinese
mainland
|
RFA:53
MWA:49
|
54.5 ± 11.7
50.1 ± 13.7
|
71.7
81.6
|
37.6
50.5
|
37.2
45.9
|
15.5
26.9
|
Liu,2018[25]
|
Chinese
mainland
|
RFA:436
MWA:126
|
56.0 (46.0–65.0)
54.0 (45.3–60.0)
|
96.4
99
|
80.7
94.6
|
60.4
78.9
|
22.8
46.9
|
Zhang, 2013[26]
|
Chinese
mainland
|
RFA:47/31
MWA:36/41
|
54.0 ± 10.5
54.0 ± 9.5
|
93.6/87.1
100/85.4
|
66/61.3
80.6/36.6
|
68.1/74.2
72.2/53.3
|
34/54.8
41.7/26.8
|
Thomas,2015[27]
|
Germany
|
RFA:25
MWA:28
|
57.0 ± 4.2
60.0 ± 4.2
|
100
100
|
72
79
|
96.9
97.2
|
90.6
91.7
|
Ding, 2013[28]
|
Chinese
mainland
|
RFA:85
MWA:113
|
58.64 ± 8.52
59.06 ± 11.68
|
98.7
98
|
82.7
77.6
|
80.3
75
|
39.5
32.1
|
Xu,2017[29]
|
Chinese
mainland
|
RFA:159
MWA:301
|
54.0 ± 11.0
54.2 ± 11.0
|
98.7
99.3
|
86.8
90.4
|
89.9
94.4
|
67.3
71.8
|
Casaccia,2017[30]
|
Italy
|
RFA:22
LR:24
|
60.82 ± 7.25
63.58 ± 9.55
|
81.8
95.8
|
50
58.3
|
40.9
70.8
|
18.2
25
|
Chong,2018[31]
|
Hong Kong
|
RFA:59
LR:59
|
59.3 ± 11.0
57.7 ± 10.5
|
96.6
94.9
|
78.7
88.2
|
59.3
86.3
|
25.3
68
|
Lai,2013[32]
|
Hong Kong
|
RFA:31
LR:80
|
63.1 ± 12.8
60.8 ± 9.9
|
100
92
|
92
75
|
76
76
|
40
60
|
Roberto, 2009[33]
|
Italy
|
RFA:74
LR:78
|
68.0 ± 7.0
68.0 ± 8.0
|
88
93
|
66
85
|
-
-
|
-
-
|
Shi,2014[34]
|
Chinese
mainland
|
MWA: 40/56
LR: 37/54
|
56.6 ± 9.2
54.5 ± 9.9
|
98/95
97/94
|
78/72
82/71
|
85/73
89/83
|
54/34
66/57
|
Sun,2020[35]
|
Chinese
mainland
|
MWA: 51/64
LR: 41/73
|
57.5 ± 9.6
54.5 ± 10.4
|
97.4/90.9
100/93.5
|
76.9/65.5
83.3/74.2
|
82.1/68.4
88.6/83.6
|
51.1/41.1
72.7/58.8
|
Zheng,2020[36]
|
Chinese
mainland
|
MWA:93
LR:300
|
58.0 (49.0–66.0)
57.0 (50.0–62.0)
|
90.1
95.5
|
82.9
86.1
|
66.3
71.8
|
46.1
44.3
|
Hu,2019[37]
|
Chinese
mainland
|
MWA:64
CA:56
|
55.2 ± 7.2
54.9 ± 11.3
|
85.8
92
|
63.5
87.4
|
77.8
81.4
|
49
58.5
|
Zhang, et al[38]
|
Chinese
mainland
|
RFA:67
LR:67
|
57.78 ± 10.97
57.51 ± 8.37
|
98.5
98.5
|
85
90
|
79.1
88.1
|
47
65.6
|
Cheng, et al[39]
|
Hong Kong
|
RFA:31
LR:99
|
65.48 ± 11.73
63.60 ± 9.86
|
96
97.9
|
79.2
96.2
|
60.6
91.2
|
37.9
66.7
|
Li, et al[40]
|
Taiwan
|
RFA:58
LR:58
|
61 (34–80)
61 (35–82)
|
91.4
98.2
|
77.2
88.8
|
91.2
93
|
62
71.7
|
Feng, et al[41]
|
Chinese
mainland
|
MWA: 121
LR: 121
|
57.1 ± 9.7
57.0 ± 8.4
|
98.3
96.5
|
84.7
81.8
|
81.8
85.4
|
54.4
67.8
|
2. Assessment of network
A complete triangular relationship existed between RFA, MWA and CA, the same as RFA, MWA and LR (Fig. 2), allowing for computational consistency and inconsistency analysis. Survival and disease-free survival outcomes for LR, RFA, MWA and CA combined are represented in statistical tables and forest plots (Table 2, Fig. 3). SUCRA curve plots reflected the different treatment modalities with different outcome rankings (Fig. 4). The largest number of studies was between RFA and MWA. Only one study addressed all outcomes between MWA and CA, the same situation arose between RFA and CA, and there was no direct study evidence between LR and CA.
Table 2
League table demonstrating the hazard ratio along with 95% CrI for OS rates at 1 year and 3 years and hazard ratio along with 95% CrI for DFS rates at 1 year and 3 years for different treatment comparison
|
Treatment
|
LR
|
RFA
|
MWA
|
CA
|
OS rates at 1 year
|
LR vs
|
NA
|
0.71(0.45,1.13)
|
0.81(0.49,1.35)
|
1.22(0.41,3.64)
|
|
RFA vs
|
1.40(0.88,2.22)
|
NA
|
1.14(0.74,1.76)
|
1.71(0.60,4.87)
|
|
MWA vs
|
1.23(0.74,2.04)
|
0.88(0.57,1.36)
|
NA
|
1.50(0.55,4.13)
|
|
CA vs
|
0.82(0.27,2.44)
|
0.58(0.21,1.66)
|
0.67(0.24,1.83)
|
NA
|
DFS rates at 1 year
|
LR vs
|
NA
|
0.54(0.38,0.76)
|
0.57(0.39,0.82)
|
0.76(0.31,1.87)
|
|
RFA vs
|
1.85(1.31,2.61)
|
NA
|
1.05(0.76,1.46)
|
1.41(0.60,3.33)
|
|
MWA vs
|
1.76(1.22,2.54)
|
0.95(0.69,1.32)
|
NA
|
1.35(0.57,3.16)
|
|
CA vs
|
1.31(0.53,3.21)
|
0.71(0.30,1.67)
|
0.74(0.32,1.74)
|
NA
|
OS rates at 3 year
|
LR vs
|
NA
|
0.79(0.53,1.17)
|
0.74(0.48,1.14)
|
1.52(0.53,4.34)
|
|
RFA vs
|
1.27(0.86,1.89)
|
NA
|
0.95(0.65,1.38)
|
1.52(0.53,4.34)
|
|
MWA vs
|
1.35(0.88,2.07)
|
1.06(0.73,1.54)
|
NA
|
2.04(0.75,5.58)
|
|
CA vs
|
0.66(0.23,1.89)
|
0.52(0.19,1.41)
|
0.49(0.18,1.34)
|
NA
|
DFS rates at 3 year
|
LR vs
|
NA
|
0.52(0.38,0.72)
|
0.53(0.37,0.76)
|
0.59(0.25,1.40)
|
|
RFA vs
|
1.92(1.38,2.66)
|
NA
|
1.02(0.74,1.42)
|
1.13(0.49,2.58)
|
|
MWA vs
|
1.88(1.31,2.68)
|
0.98(0.70,1.36)
|
NA
|
1.11(0.48,2.53)
|
|
CA vs
|
1.70(0.71,4.03)
|
0.89(0.39,2.02)
|
0.90(0.40,2.07)
|
NA
|
For one-year and three-year survival rates, HR values are all within the credible interval. In a conclusion, the one-year and three-year survival rates of the four treatment modalities were not statistically different between the groups.
For one-year disease-free survival, the HR was 1.85 (95% CrI: 1.31–2.61) for LR compared with RFA; 1.76 (95% CrI: 1.22–2.54) for LR compared with MWA. In a word, the one-year disease-free survival rate of LR was statistically higher than that of RFA and MWA, but there was no difference between the remaining groups compared.
For three-year disease-free survival, the HR was 1.92 (95% CrI: 1.38–2.66) for LR compared with RFA; 1.88 (95% CrI: 1.31–2.68) for LR compared with MW. As with one-year disease-free survival, LR had a higher disease-free survival rate than RFA and MWA at 3 years after the procedure, and there was no statistical difference in the three-year disease-free survival rates between LR and CA and the remaining three ablation modalities.
3. Rank probability of treatment regimen
The ranked probability ranking of the effect of all four treatment modalities has been reflected by the SUCRA graph (Fig. 4), for one-year survival, CA ranked first with the highest probability (0.621), followed by LR (0.316), then MWA and RFA; for one-year disease-free survival, LR ranked first (0.72), followed by CA (0.28), with both RFA and MWA ranking last; for three-year survival, again CA ranked first with the highest probability (0.773), followed by LR (0.206), then RFA (0.04) and MWA (0.01) with little difference; and for three-year disease-free survival, HR ranked first (0.882), followed by CA ( 0.118), with RFA and MWA ranking last both.
4. Consistency and inconsistency analysis
In order to develop the final model with all resulting parameters, the adaptation iterations were kept at 20,000, thus removing the initial simulation results of the MCMC analysis from the model. The number of simulation iterations was kept at 100,000. The sparsity factor was kept at 10 and the number of chains was kept at 4. For all four outcome parameters, the assessment of the adequacy of convergence of the Gelman Rubin diagnosis was below 1.05 (Fig. 5). The p value of inconsistency in Node-Splitting was not significant for all the outcomes in all the comparisons (Table 5).
5. Subgroup analysis
The subgroups of the study were divided into Chinese mainland and other regions according to the region where the study was conducted. Among the included literature, only Chinese mainland was analyzed for cryoablation efficacy, so the Chinese mainland subgroup included four different treatment modalities, while the other regions subgroup included only three treatment modalities such as partial hepatectomy, thermal radiofrequency ablation and microwave ablation.
In the subgroup from Chinese mainland, LR had an advantage in disease-free survival rates over MWA, but there was no statistical difference between other treatment modalities for both overall survival and disease-free survival. For LR compared with MWA in the one-year disease-free survival, the HR was 1.49 (95% CrI: 1.02–2.17); in the three-year disease-free survival, the HR was 1.63 (95% CrI: 1.04–2.56). In contrast, in the subgroups from regions other than Chinese mainland, there were no statistically significant differences between the three treatment modalities for one- and three-year overall survival rates. For the one-year disease-free survival, the HR was 2.33 (95% CrI: 1.34–4.06) for LR compared with RFA; 3.14 (95% CrI: 1.16–8.49) for LR compared with MWA (Table 3). For the three-year disease-free survival, the HR was 2.28 (95% CrI: 1.54–3.39); and for LR versus MWA, the HR was 2.74 (95% CrI: 1.26–5.59) (Table 4).
Table 3
League table of China mainland group
|
Treatment
|
LR
|
RFA
|
MWA
|
CA
|
OS rates at 1 year
|
LR vs
|
NA
|
0.80(0.42,1.55)
|
0.95(0.56,1.64)
|
1.47(0.54,4.01)
|
|
RFA vs
|
1.24(0.65,2.40)
|
NA
|
1.19(0.74,1.91)
|
1.83(0.72,4.65)
|
|
MWA vs
|
1.05(0.61,1.80)
|
0.84(0.52,1.36)
|
NA
|
1.54(0.65,3.65)
|
|
CA vs
|
0.68(0.25,1.86)
|
0.55(0.22,1.39)
|
0.65(0.27,1.54)
|
NA
|
DFS rates at 1 year
|
LR vs
|
NA
|
0.65(0.41,1.03)
|
0.67(0.46,0.98)
|
0.91(0.39,2.09)
|
|
RFA vs
|
1.55(0.97,2.46)
|
NA
|
1.04(0.73,1.49)
|
1.41(0.65,3.04)
|
|
MWA vs
|
1.49(1.02,2.17)
|
0.96(0.67,1.37)
|
NA
|
1.35(0.63,2.89)
|
|
CA vs
|
1.10(0.48,2.54)
|
0.71(0.33,1.54)
|
0.74(0.35,1.59)
|
NA
|
OS rates at 3 year
|
LR vs
|
NA
|
0.84(0.48,1.47)
|
0.88(0.56,1.40)
|
1.69(0.63,4.55)
|
|
RFA vs
|
1.19(0.68,2.08)
|
NA
|
1.05(0.67,1.65)
|
2.01(0.82,4.95)
|
|
MWA vs
|
1.13(0.72,1.79)
|
0.95(0.61,1.49)
|
NA
|
1.92(0.78,4.73)
|
|
CA vs
|
0.59(0.22,1.59)
|
0.50(0.20,1.22)
|
0.52(0.21,1.29)
|
NA
|
DFS rates at 3 year
|
LR vs
|
NA
|
0.63(0.36,1.09)
|
0.61(0.39,0.96)
|
0.69(0.26,1.84)
|
|
RFA vs
|
1.59(0.92,2.75)
|
NA
|
0.97(0.63,1.51)
|
1.10(0.45,2.70)
|
|
MWA vs
|
1.63(1.04,2.56)
|
1.03(0.66,1.60)
|
NA
|
1.13(0.46,2.77)
|
|
CA vs
|
1.44(0.54,3.83)
|
0.91(0.37,2.22)
|
0.88(0.36,2.16)
|
NA
|
Table 4
League table of other regions group
|
Treatment
|
LR
|
RFA
|
MWA
|
OS rates at 1 year
|
LR vs
|
NA
|
0.67(0.27,1.63)
|
0.50(0.10,2.40)
|
|
RFA vs
|
1.50(0.61,3.67)
|
NA
|
0.75(0.20,2.82)
|
|
MWA vs
|
1.99(0.42,9.52)
|
1.33(0.35,4.95)
|
NA
|
DFS rates at 1 year
|
LR vs
|
NA
|
0.43(0.25,0.75)
|
0.32(0.12,0.86)
|
|
RFA vs
|
2.33(1.34,4.06)
|
NA
|
0.74(0.33,1.70)
|
|
MWA vs
|
3.14(1.16,8.49)
|
1.35(0.59,3.08)
|
NA
|
OS rates at 3 year
|
LR vs
|
NA
|
0.75(0.41,1.36)
|
0.47(0.18,1.19)
|
|
RFA vs
|
1.33(0.74,2.41)
|
NA
|
0.62(0.30,1.29)
|
|
MWA vs
|
2.15(0.84,5.51)
|
1.61(0.77,3.35)
|
NA
|
DFS rates at 3 year
|
LR vs
|
NA
|
0.44(0.29,0.65)
|
0.37(0.17,0.80)
|
|
RFA vs
|
2.28(1.54,3.39)
|
NA
|
0.83(0.43,1.63)
|
|
MWA vs
|
2.74(1.26,5.95)
|
1.20(0.61,2.34)
|
NA
|
Table 5
Consistency inspection and Nodal analysis
Statistical approach
|
Inconsistent model
|
Nodel difference method
|
|
p-value
|
Minimum p-value
|
one-year OS rates
|
0. 48
|
0.480
|
three-year DFS rates
|
0.36
|
0.537
|
three-year OS rates
|
0.39
|
0.188
|
three-year DFS rates
|
0.79
|
0.491
|
6. Risk assessment
For the included RCT studies, we used the ROB 2 tool to assess the risk of bias, and the results are shown below (Table 6); while for the non-RCT category, we used the NOS scale for scoring, and the final assessment scores were above 6 (Table 7), suggesting less bias. In addition, analysis of the data using Stata revealed a symmetrical funnel plot distribution (Fig. 6), suggesting a small and acceptable publication bias.
Table 6
Assessment of the RCTs included in meta-analysis
SI. no.
|
RoB 2 parameters
|
Ng, et al[4]
|
Lee, et al[5]
|
Chong, et al[12]
|
Feng, et al[20]
|
Wang, et al[21]
|
1
|
Randomization process
|
LR
|
LR
|
LR
|
LR
|
LR
|
2
|
Concealment of allocation
|
LR
|
LR
|
LR
|
SC
|
LR
|
3
|
Assignment to intervention
|
LR
|
LR
|
LR
|
LR
|
LR
|
4
|
Adhering to intervention
|
LR
|
LR
|
LR
|
LR
|
LR
|
5
|
Missing of outcome data
|
LR
|
LR
|
LR
|
LR
|
LR
|
6
|
Measurements of outcome
|
LR
|
LR
|
LR
|
LR
|
LR
|
7
|
Selection of reported results
|
SC
|
SC
|
LR
|
LR
|
LR
|
8
|
Overall bias
|
SC
|
SC
|
LR
|
SC
|
LR
|
LR low risk, RR high risk, SC some concern |
Table 7
Assessment of the nRCTs included in meta-analysis
Study
|
Selection
|
|
|
|
Comparability
|
Exposure
|
|
|
Scores
|
|
Representativeness (expressed cohort)
|
Selection (non-exposed cohort)
|
Ascertainment (exposure)
|
no outcome events before study
|
Control for important
|
Assessment of outcome events
|
Adequacy of follow-up
|
Completeness of follow-up
|
Ohmoto, et al[22]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Theodora,et al[23]
|
1
|
1
|
1
|
0
|
1
|
1
|
1
|
1
|
7
|
Lu, et al[24]
|
1
|
1
|
1
|
0
|
1
|
1
|
1
|
1
|
7
|
Liu, et al[25]
|
1
|
1
|
1
|
0
|
1
|
1
|
1
|
1
|
7
|
Zhang, et al[26]
|
1
|
1
|
1
|
0
|
1
|
1
|
1
|
1
|
7
|
Thomas, et al[27]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Ding, et al[28]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Xu, et al[29]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Casaccia,et al[30]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Chong, et al[31]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Lai, et al[32]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Roberto, et al[33]
|
1
|
1
|
1
|
0
|
1
|
1
|
1
|
1
|
7
|
Shi, et al[34]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Sun, et al[35]
|
1
|
1
|
1
|
0
|
1
|
1
|
1
|
1
|
7
|
Zheng, et al[36]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Hu, et al[37]
|
1
|
1
|
1
|
0
|
1
|
1
|
1
|
1
|
7
|
Zhang, et al[38]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Cheng, et al[39]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Li, et al[40]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|
Feng, et al[41]
|
1
|
1
|
1
|
0
|
2
|
1
|
1
|
1
|
8
|