Abstract
Background:
The increased incidence of primary papillary thyroid microcarcinoma (PTMC) has led to increased research in the field of nonsurgical therapeutic options for those who refuse surgery or are at high risk for surgery. The study aimed at comprehensively evaluating the efficacy and safety of thermal ablation techniques for the treatment of PTMC via a systematic review and meta-analysis.
Methods:
The Pubmed MEDLINE and EMBASE databases were searched for studies reporting the efficacy and safety of thermal ablations (radiofrequency-, laser-, and microwave-ablations [RFA, LA, and MWA]) until August 10, 2019. A review of 105 potential papers identified 11 eligible papers, including 715 patients. The pooled proportions of complete disappearance and recurrence, and the pooled estimates of mean volume reduction and its rate of the treated PTMC were assessed by using random-effects modeling. The pooled proportions of overall and major complications were calculated. Subgroup analysis was performed according to the treatment modality. Between-study heterogeneity was explored by using χ 2 statistic for pooled estimates and inconsistency index I 2. Quality of the studies was evaluated by using the Risk of Bias Assessment Tool for Nonrandomized Studies.
Results:
The pooled proportions of complete disappearance and recurrence of PTMC were 57.6% [95% confidence interval (CI): 35.4–79.8] and 0.4% [95% CI: 0–1.1], respectively. The pooled estimates of mean volume reduction and its rate were 73.5 mm3 [52.4–94.6 mm3] and 98.1% [95% CI: 96.7–99.5], respectively. The pooled proportions of overall and major complications were 3.2% [95% CI: 1.1–5.2] and 0.7% [95% CI: 0–1.5], respectively. Significant between-study heterogeneity was observed for complete disappearance (p < 0.001, I 2: 99%), mean volume reduction (p < 0.001, I 2: 93%), and its rate (p < 0.001, I 2: 86%). Subgroup analysis revealed heterogeneity of the complete disappearance proportion among the treatment modality (I 2 range: 95–100%). RFA showed the highest mean volume reduction rate (99.3%), followed by MWA (95.3%) and LA (88.6%) (p < 0.001).
Conclusions:
All thermal ablation techniques were effective and safe for the treatment of PTMC. However, each treatment modality had significant heterogeneity with respect to complete disappearance of PTMC. Compared with RFA and MWA, LA was less effective in reducing the volume of PTMC.
Introduction
Among the subtypes of thyroid carcinomas, papillary thyroid carcinoma is the most common, accounting for more than 70% of the thyroid malignancies (1). Papillary thyroid microcarcinoma (PTMC) is defined as papillary thyroid carcinoma measuring less than 1.0 cm in size and has very favorable prognosis, with a disease-specific mortality rate of less than 1% (2,3). Among the treatment options for PTMC, surgical management is the first-line treatment (4). However, it is not suitable for patients for thyroidectomy may be high-risk surgery or for those refusing surgery (5). Active surveillance is also an alternative management approach recommended by the American Thyroid Association (4); however, there is a risk of locoregional metastasis (2–6%) (2,3).
Accumulating evidence suggests that thermal ablation techniques, including radiofrequency ablation (RFA), microwave ablation (MWA), and laser ablation (LA), are effective and safe treatment modalities for both benign thyroid nodules (6 –9) and recurrent thyroid cancers (10 –13). Compared with surgical treatment, similar efficacy and lower rates of complications have been reported for RFA (10,12) and MWA (14) for the treatment of recurrent papillary thyroid carcinoma and primary PTMC, respectively.
Recent studies have reported that thermal ablations are effective and safe for treating primary PTMC (14 –25). However, there are still concerns and comprehensive assessment of thermal ablation techniques in terms of the treatment efficacy and safety of PTMC is necessary to optimize personalized care in the face of uncertain evidence and to bridge the evidence gaps in the future. Therefore, this study aimed at comprehensively evaluating the efficacy and safety of thermal ablation techniques for the treatment of PTMC. To the best of our knowledge, this is the first systematic review and meta-analysis evaluating the efficacy and safety of different thermal ablations for the treatment of PTMC.
Materials and Methods
Literature search strategy
This meta-analysis was performed according to the Preferred Reporting Items for a Systematic Review and Meta-analysis (PRISMA) guidelines (26). The PubMed, MEDLINE, and EMBASE databases were searched to identify publications on thermal ablation in treating PTMC. No limitation was placed on the beginning date of literature, and the search was finalized on August 10, 2019. The following search terms were used: [(thyroid) AND (microcarcinoma OR micropapillary) AND (carcinoma OR cancer OR malignancy) AND (radiofrequency ablation OR RFA OR LA OR laser ablation OR microwave ablation OR thermal ablation)]. The bibliographies of the selected articles were further reviewed to search for other relevant articles. Two reviewers (blinded and blinded) independently performed the literature search; any disagreement in terms of inclusion was resolved by discussion until unanimity was attained by all authors. We were not blinded to authors, institutions, journals, or interventions while selecting studies or extracting data.
Inclusion and exclusion criteria
Studies investigating thermal ablation techniques (i.e., RFA, LA, and MWA) for the treatment of PTMC were eligible for inclusion. The included studies met the following criteria: (1) population and intervention: Patients were diagnosed with primary PTMC; the main treatment modality was either RFA, LA, or MWA; and sample size was ≥5 patients. (2) Study design: Both retrospective and prospective observational studies were included. (3) Outcomes: Studies reported details on treatment outcomes (i.e., volume reduction, complete disappearance, recurrence of nodules) and complications.
Any study meeting one of the following exclusion criteria was excluded: (1) case reports or series with <5 patients; (2) abstracts, reviews, comments, and editorials; (3) articles not written in English; (4) studies with partially overlapping cohorts. Two reviewers (blinded and blinded) independently reviewed the literature in consensus.
Data extraction
From the eligible studies, the following data were extracted into standardized formats: (a) study characteristics: authors, year of publication, affiliation, duration of study, study design, and sample size; (b) patients' demographic and clinical characteristics: age, sex, ablation technique (RFA, MWA, or LA), nodule characteristics (size and volume), operators' experience, and follow-up duration; (c) outcomes: mean volume reduction (mm3) of nodules between treatment sessions, volume reduction rate (%) calculated as [(initial volume − final volume) × 100/initial volume], complete disappearance (determined on the final follow-up ultrasound [US] scan), recurrence (determined pathologically), and volume reduction rate; and (d) major and minor complications per patient. One reviewer (blinded) extracted the data, and another reviewer (blinded) confirmed the data's validity.
Quality assessment
The quality of the included studies was then independently assessed in consensus by two authors (blinded and blinded) by using the structured criteria of the Risk of Bias Assessment Tool for Non-randomized Studies (RoBANS) (27).
Definitions of complications
Definitions of major and minor complications were based on definitions by the Society of Interventional Radiology (28,29). A major complication was defined as one that might be life-threatening, result in substantial morbidity or disability, or prolong hospitalization if left untreated. In this study, major complications included transient (lasting >1 month) or permanent voice changes and hypothyroidism. Minor complications included immediate hoarseness (lasting <1 week), hematoma, and moderate pain requiring medication to relieve it. Side effects—not requiring therapy or prescription medication—included mild, transient postprocedural pain, heat sensation, and neck discomfort.
Data synthesis for assessing the pooled proportions and estimates of treatment efficacy and complications
The primary outcomes of this meta-analysis were the pooled proportions of complete disappearance of nodule, recurrence of nodule, overall complications and major complications, and pooled estimates of mean volume reduction rates of nodules. Meta-analytic pooling was performed by using the inverse variance method to calculate weights. Pooled proportions with 95% confidence intervals (CIs) were obtained with the aid of Der Simonian-Laird random-effects modeling (30 –32). Subgroup meta-regression analysis was performed by the treatment modality (RFA, MWA, and LA).
Between-study heterogeneity was calculated by the χ 2 statistic for pooled estimates (p < 0.05, indicating significant heterogeneity) and the Higgins inconsistency index I 2 (0–40%, may have insignificant heterogeneity; 30–60%, 50–90%, and 75–100%, may have moderate, substantial, and considerable heterogeneity, respectively) (33). Publication bias was evaluated by using funnel plots and Egger's test, with a p-value <0.1 indicating significant bias (34). Publication bias-adjusted pooled proportions and estimates were also calculated via the trim-and-fill method (35), where agreement between the unadjusted and adjusted pooled proportions and estimates indicates minimal publication bias. All statistical meta-analyses were performed by one author (blinded) using the “meta” and “metafor” packages in R statistical and computing software (version 3.5.1, Vienna, Austria).
Results
Literature search
The selection process of studies from the literature search is depicted in Figure 1. The initial literature search identified 105 articles. After removal of duplicates, 90 were screened for eligibility. Of these, 77 were excluded after reviewing the abstracts (20 review articles, 1 case report, 4 letters/editorials/abstracts, 42 articles of irrelevant contents, and 10 non-English articles). The full texts of the remaining 13 articles were then thoroughly reviewed; two articles were further excluded due to a small sample size (<5 patients) (36) and having partially overlapping cohorts (15). Between the two studies with overlapping cohorts, the study with a larger sample size was chosen. Finally, 11 articles with a total of 715 patients with 748 primary PTMC were included (14 –20,22–25). No additional eligible study was identified after a search of bibliographies within these articles.

A flow diagram of PRISMA.
Characteristics of included studies
The detailed characteristics of the 11 eligible studies are shown in Table 1. Among the 11 original articles included, there were 7 retrospective studies (14,15,17 –20,25), 2 prospective studies (22,24), and 2 studies with unclear study design (16,23). All included studies were carried out in single centers. The cohort sizes ranged from 6 to 185 patients, with numbers of PTMC ranging from 6 to 206 total nodules. The mean ages of included studies' patients ranged from 42.2 to 59.5 years. The mean largest diameters of PTMC ranged from 1.8 to 7.3 mm. Nine of 11 studies had a mean volume reduction rate of more than 90% (range, 56.3–100%). The mean follow-up periods of included studies ranged from 7.8 to 25.7 months, with six studies reporting more than 1 year of mean follow-up periods (14,15,18 –20,25). The type of thermal ablation techniques included four MWA (14,18,22,23), four RFA (16,17,20,24), and three LA (15,19,25). The experience of operators was reported in seven studies (range, 2–20 years of experience). All study settings were based in Asian countries—10 in China (14 –19,22–25) and 1 in Korea (20). A single thermal ablation session was performed in four studies (14,15,20,22), while one study indicated a maximum of two LA sessions (19). Six studies did not describe the number of sessions performed (16 –18,23–25). The quality assessment by RoBANS showed moderate overall quality of the included studies, with all studies satisfying at least five of the eight criteria (Fig. 2). All studies described the thermal ablation techniques and equipment clearly.

The RoBANS criteria of the included studies. RoBANS, Risk of Bias Assessment Tool for Nonrandomized Studies.
Characteristics of Selected Studies
LA, laser ablation; MWA, microwave ablation; N/A, not available; RFA, radiofrequency ablation; SD, standard deviation.
Pooled estimates of mean volume reduction rates
The pooled estimates of mean volume reduction and its rate were 73.5 mm3 (95% CI: 52.4–94.6 mm3) and 98.1% (95% CI: 96.7–99.5), respectively (Fig. 3A, B). There were significant heterogeneities for both (mean volume reduction, I 2 = 93%; mean volume reduction rate, I 2 = 86%; Cochran's Q method, both, p < 0.001). The funnel plot and Egger's test showed a significant publication bias (Fig. 4A, B, p < 0.001). After the trim-and-fill method adjusting publication bias, the adjusted pooled estimates for mean volume reduction and its rate were 61.2 mm3 (95% CI: 48.8–73.7 mm3) and 99.2% (95% CI: 97.6–100), respectively.

Forest plots of pooled estimates of (

Funnel plots of (
Pooled proportions of complete disappearance and recurrence of nodules
Among 728 nodules, 415 nodules showed complete disappearance after treatment. Among eight studies reporting mean follow-up periods, the mean follow-up of completely disappeared PTCM (n = 298) was 21.1 months (range, 11–25.7 months). The pooled proportion of complete disappearance was 57.6% (95% CI: 35.4–79.8; I 2 = 99%) (Table 2 and Fig. 3C). The funnel plot showed symmetric distribution of complete disappearances of nodules, and publication bias was not significant in Egger's test (p = 0.273) (Fig. 4C). Significant heterogeneity was observed in the proportion of complete disappearance of nodules (p < 0.001). Nine recurrences among 728 nodules were reported in three studies (14,15,18). The pooled proportion of recurrence was 0.4% (95% CI: 0–1.1; I 2 = 0%) (Table 2 and Fig. 3D). In both the funnel plot and Egger's test, significant publication bias was observed for recurrences (p = 0.055) (Fig. 4D). However, the trim-and-fill estimate of recurrence (0.6%) was inconsistent with the pooled proportion. No heterogeneity was observed in the reported recurrence rate (p = 0.704).
Summary Estimates of the Meta-Analytic Pooled Proportions and Estimates for Efficacy and Safety of Thermal Ablations
p-Value by the Cochran's Q method to test the heterogeneity of the pooled data (p < 0.05 indicating significant heterogeneity).
Higgins' index for heterogeneity (0–40% may have insignificant heterogeneity; 30–60%, 50–90%, and 75–100% may indicate moderate, substantial, and considerable heterogeneity, respectively).
Egger's test (p < 0.10 indicating significant publication bias).
Pooled estimates for mean volume reduction with 95% CI.
Pooled estimates for mean volume reduction rate calculated as [(initial volume − final volume) × 100/initial volume].
CI, confidence interval.
Pooled proportions of complications
Thirty-four complications among 728 nodules in 715 patients were reported. The pooled proportion of major complications was 0.7% (95% CI: 0–1.5; I 2 = 0%) (Table 2 and Fig. 5A). The funnel plot showed asymmetric distribution of major complications, but publication bias was not significant in Egger's test (p = 0.199) (Fig. 6A). The pooled proportion of overall complications was 3.2% (95% CI: 1.1–5.2; I 2 = 34%) (Table 2 and Fig. 5B). In both the funnel plot and Egger's test, no significant publication bias was seen for overall complications (p = 0.247) (Fig. 6B). No heterogeneity was noted in both overall (p = 0.123) and major (p = 0.8) complications. The trim-and-fill estimates of overall and major complications were consistent with pooled proportions, indicating little to no publication bias (Table 2).

Forest plots of pooled proportions of (

Funnel plots for assessing publication bias of (
Details on major and minor complications
All 11 major complications were voice changes—10 transient voice changes lasting for at least 1 month (14,18) and 1 permanent voice change (14). Twenty-three minor complications were reported, including 3 instances of moderate pain, 2 abnormal thyroid hormone levels, 11 hematomas, and 7 instances of immediate hoarseness. All immediate hoarseness occurred soon after the completion of the procedure and resolved within 10 minutes to 1 week (16,18,23,24). Other than the moderate pain relieved by medication, all minor complications resolved spontaneously within 3 months. The frequency of major and minor complications was 1.5% (11/715) and 3.2% (23/715), respectively.
Subgroup analysis according to the treatment modality
Subgroup analysis according to the treatment modality was performed for the treatment efficacy outcomes showing significant heterogeneity (i.e., complete disappearance and mean volume reduction rate) (Table 3). There were no significant differences among the treatment modalities for complete disappearance of nodules (p = 0.351). However, all modalities showed large heterogeneity (I 2: 95–100%). For the mean volume reduction rate, RFA showed the highest volume reduction rate (99.3%) followed by MWA (95.3%) and LA (88.6%), and the differences were statistically significant (p < 0.001). However, no significant differences were found among the treatment modalities in the subgroup analysis of overall (p = 0.286) and major (p = 0.406) complications (Table 4).
Subgroup Analysis of Complete Disappearance, Mean Volume Reduction, and Mean Volume Reduction Rate of Nodules According to the Treatment Modality
Higgins' index for heterogeneity (0–40%, may have insignificant heterogeneity; 30–60%, 50–90%, and 75–100% may indicate moderate, substantial, and considerable heterogeneity, respectively).
Egger's test for publication and reporting bias (p < 0.10 indicating significant publication bias).
Calculated as [(initial volume − final volume) × 100/initial volume].
LA, laser ablation; MWA, microwave ablation; RFA, radiofrequency ablation.
Subgroup Analysis of Overall and Major Complications According to the Treatment Modality
Higgins' index for heterogeneity (0–40%, may have insignificant heterogeneity; 30–60%, 50–90%, and 75–100% may indicate moderate, substantial, and considerable heterogeneity, respectively).
Egger's test (p < 0.10 indicating significant publication bias).
Discussion
In the current systemic review and meta-analysis evaluating the efficacy and safety of thermal ablation techniques (i.e., RFA, MWA, LA), treatment efficacy was favorable, as noted by the pooled proportion of complete disappearance of PTMC (57.6% [95% CI: 35.4–79.8]) and that of recurrence (0.4% [95% CI: 0–1.1]). The pooled estimate for mean volume reduction rate of PTMC after thermal ablation was 98.1% [95% CI: 96.7–99.5]. In addition, the pooled proportion of major complications was very low (0.7% [95% CI: 0–1.5]), and all of them were non-life-threatening voice changes. Based on these results, all three thermal ablation techniques demonstrated efficacy and safety for the treatment of PTMC.
In subgroup analysis by treatment modality, the mean volume reduction rate of PTMC was comparable between RFA (99.3%) and MWA (95.3%), whereas that of LA was lower (88.6%). This finding is consistent with a previous study by Korkusuz et al., who found comparable volume reductions between RFA and MWA in the treatment of benign thyroid nodules (37). Further, the relatively larger volume reduction of RFA compared with LA is consistent with a previous meta-analysis comparing the efficacy of RFA and LA, which reported significantly higher volume reduction after RFA (77.8% [95% CI: 67.7–88.0]) than after LA (49.5% [95% CI: 26.7–72.4]) (38). We speculate the superior efficacy of RFA to the use of the moving-shot technique, which allows multi-dimensional ablation of the entire tumor safely and effectively (38,39). Interestingly, six of eight included studies of either RFA or MWA applied the moving-shot technique (14,16,18,20,23,24). However, none of the LA studies used it due to the procedure-specific steps in LA, which simultaneously deploys multiple fibers—each with an individual energy emission—in an array around the tumor (40). Our findings were consistent with a previous study by Baek et al., who reported the efficacy of RFA to be slightly superior to that of LA, and that attributed this to the use of the moving-shot technique and an internal cooling system (40). However, our findings were in PTMC and not benign thyroid nodules reported in these studies and that were greater in size (37,38,40).
Overall, thermal ablation techniques were found to be safe in the treatment of PTMC. The pooled proportions of both overall (3.2% [95% CI: 1.1–5.2]) and major complications (0.7% [95% CI: 0–1.5]) were similar to those previously reported in a meta-analysis evaluating the safety of RFA for the treatment of benign thyroid nodules and recurrent thyroid cancers (pooled proportions of overall and major complications: 2.38% and 1.35%, respectively) (41). However, in a meta-analysis evaluating the safety and efficacy of MWA for benign thyroid nodules and PTMC, the pooled proportions of major and minor complications were 11.5% and 5.1%, respectively (42), higher than those reported in this study. Indeed, in the subgroup analysis of safety by treatment modality, we also found MWA to have higher pooled proportions of overall and major complications than the other two modalities, although the difference was not statistically significant. This slight difference may be due to variation in operators' experience or procedural techniques, and thus may not necessarily imply that MWA is safer than the other two modalities. Aside from subtle differences in the frequency of complications by treatment modality, no major complications were reported to be life-threatening, confirming that minimal risks are associated with all thermal ablations.
With respect to study heterogeneity, very little heterogeneity was present for the recurrence of PTMC (p = 0.702; I 2 = 0%). However, significant between-study heterogeneity was noted for the complete disappearance of PTMC (p < 0.001; I 2 = 99%). After subgroup analysis by treatment modality, large heterogeneity was noted for each modality with respect to complete disappearance of PTMC, with I 2 of more than 90% for all modalities. We assume that a small number of included studies within each subgroup (<4 studies) might be attributable to this finding. Further, significant between-study heterogeneity was noted in the mean volume reduction rate of PTMC (p < 0.001; I 2 = 86%). The MWA accounted for most of the heterogeneity (I 2 = 92%), whereas LA showed no heterogeneity (I 2 = 0%), possibly suggesting that LA may have more consistent treatment effects than MWA.
All three modalities are based on thermal ablation as the fundamental mechanism of delivering therapeutic effect. The application of localized heating on thyroid tissue creates irreversible cellular damage in LA (36,43), RFA (40,44), and MWA (45). In addition to the widely known benefits of thermal ablations over surgery, such as avoiding surgical scarring and postoperative risks, US-guided delivery of heat specifically to the lesion allows preservation of surrounding normal thyroid tissues and, consequently, of thyroid function. One issue to consider is that lesion-specific targeting may have a risk of overlooking other very tiny microcarcinomas not visible on US; however, the same risk exists for lobectomy when tiny microcarcinomas are in the contralateral lobe (36). Among the included studies, the longest mean follow-up period was two years, which may not be long enough to confidently determine the efficacy of thermal ablations. Although not included in our data synthesis due to a small sample size (consisting of small papillary thyroid carcinomas not strictly defined as PTMC), a study by Kim et al. included four patients with PTMC who were treated with RFA (46); over a long follow-up period ranging from 36 to 65 months, all four PTMC showed complete disappearance without recurrence. Future studies are needed to further validate the efficacy and safety of thermal ablations with a longer follow-up period.
There are a few limitations to our study. First, the included studies were regionally concentrated in Asia, and our meta-analytic results may not be similar in other regions where clinical practice settings may differ. Second, the number of included studies was relatively small because treatment of PTMC by thermal ablations was introduced recently, as evidenced from the years of publication of included studies (8 of 11 studies published between 2018 and 2019). In addition, the operators' experiences were not included in the subgroup analysis due to insufficient information and variability of reported data.
In conclusion, the current systematic review and meta-analysis provides a comprehensive review of results on the efficacy and safety of thermal ablation techniques—RFA, MWA, and LA—for the treatment of PTMC. Our results suggest that thermal ablation techniques have reasonable safety and efficacy profiles. However, future research needs to investigate studies with longer follow-up and more heterogeneous populations.
Footnotes
Acknowledgment
The authors appreciate Dr. Eun Ju Ha for her time and effort in reviewing this article.
Author Disclosure Statement
No competing financial interests exist.
Funding Information
No funding was received for this article.
