Abstract
The presence of distant metastasis is associated with an adverse outcome in papillary thyroid cancer. We performed a meta-analysis to investigate the role of molecular markers as predictors for distant metastasis in papillary thyroid cancer. Four electronic databases including PubMed, Web of Science, Scopus, and Virtual Health Library were searched, and odds ratio and its 95% confidence interval concerning the association of
Introduction
Thyroid cancer is the most common among endocrine malignancies, and its incidence has been increasing rapidly over the years. 1 Papillary thyroid cancer (PTC) is the most common histological type, accounting for about 90% of thyroid cancer.2,3 Generally, the prognosis of PTC patients is excellent with a survival rate of more than 95% and usually a good response to initial treatment. 3 Nonetheless, there is a small group of PTCs that behave aggressively at presentation, develop distant metastasis (DM), and are associated with an increased rate of mortality. 4 DM is not a frequent event in PTCs, accounting about 2%–5% of cases,5–8 but the presence of DM has been well known as an adverse prognostic factor in PTCs.5,9–11 In the latest guidelines of the American Thyroid Association (ATA), differentiated thyroid cancers (DTCs) with DM are classified as ATA high risk. 4
Several clinicopathological factors such as old age, large size, vascular invasion, and extrathyroidal extension have been demonstrated to be risk factors for DM in DTC.5,12,13 However, there are ongoing debates in other studies on this topic.14–17 There has been rapid development in understanding the pathogenesis and genetic profiles of thyroid cancer in recent years. A few genetic events have been shown to be associated with aggressive behaviors and poor outcomes in PTC.6,18,19 In addition to having worse impact on patients’ survival, the survival rate of thyroid cancer patients with DM has not improved over the last two decades despite the fact that there have been huge improvements in patient management and treatment strategies. 20 Therefore, it is really needful to identify reliable molecular markers that are predictive of DM in PTCs.
In this study, we aimed to perform a systematic review and meta-analysis of observational studies to investigate the role of molecular markers in predicting DM in PTCs.
Methods
Literature search
We searched for relevant studies in four electronic databases including PubMed, Scopus, Web of Science, and Virtual Health Library from inception to November 2016. We used the following search terms: “
Selection criteria and abstract screening
We imported all search results from four libraries into Endnote (Thomson Reuters) and deleted duplicates. Two reviewers independently screened the titles and abstracts of all included articles using the following inclusion criteria: (1) studies on PTC population and (2) studies in which the authors reported the association of DM with mutational status of any of the following genetic alterations (
Full-text screening and data extraction
We downloaded corresponding full-text of all included articles following abstract screening steps, and two reviewers independently screened these full-texts. Data were extracted into a predefined extraction form. Any discrepancies were resolved again by discussion and consensus. The following details were extracted: institution, city, country, publication year, surgical period of included patients, types of genetic events, detection method, gender, age, follow-up duration, number of overall DM, DM at the time of diagnosis, and DM during follow-up. Those studies with a number of cases with DM less than five were further excluded in this step.
Risk of bias assessment
The methodologic quality of included studies was evaluated using the Newcastle–Ottawa Scale (NOS) for quality of cohort and case-control studies in our meta-analysis. 22 Stars were given for each study (maximum nine stars) based on a developed checklist. 22 Studies with at least six stars were considered as moderate- to high-quality studies and those with less than six stars were regarded as low-quality studies.
Data analysis
We used Review Manager version 5.3 (Cochrane Collaboration, Oxford, UK) for data analysis. Random-effects model was used to calculate pooled estimates of odds ratio (OR) and 95% confidence interval (CI). We assessed the heterogeneity across the included studies by the I2 statistic. 23 An I2 statistic <25% indicates low amount of heterogeneity and >50% indicates high amount of heterogeneity. 24 We used subgroup analyses and meta-regression to determine whether study-level covariates accounted for the heterogeneity among studies. Meta-regression analysis was calculated by Meta-Essentials: Workbook for meta-analysis. 25 Egger’s regression test and funnel plot observation were used to evaluate the presence of publication bias, investigated by using MAVIS version 1.1.2—a R package. 26 A p value less than 0.05 was considered to be statistically significant publication bias.
Results
A total of 7505 and 2647 articles were found after searching four electronic databases and deleting duplicates, respectively. We identified 136 potential studies following title and abstract screening. In total, 94 articles were excluded after reading full-text and leaving 42 studies with 11,253 patients for final analyses (Figure 1). Characteristics of all included studies are described in Table 1.

Study flowchart.
Characteristics of included studies.
NOS: Newcastle–Ottawa Scale; DM: distant metastasis; ND: no description; DS: direct sequencing; RT-PCR: reverse transcription polymerase chain reaction; PCR-RFLP: PCR–restriction fragment length polymorphism; PCR-SSCP: PCR–single-strand conformation polymorphism; IHC: immunohistochemistry; MASA-PCR: mutant allele–specific amplification–PCR; NGS: next generation sequencing; Pyroseq: pyrosequencing; FMCA: fluorescent melting curve analysis; OH: oligonucleotide hybridization; qPCR: quantitative real-time PCR; LH-MSA: loop-hybrid mobility shift assay.
Median value of age.
TERT promoter mutations
Overall, nine studies with 2434 PTC patients had sufficient data to be included for meta-analysis.6,12,33,34,42,44,49,52,59 Our analysis showed that 38.2% and 8.8% of PTCs with and without DM, respectively, had

Forest plots concerning the association of distant metastasis with (a)
BRAF mutations
Relevant data were found in 36 studies,6–8,12,19,27–36,38–41,43,45–51,53–57,59–62 PTC population from the following studies7,8,30–32,50,57,59 could overlap with PTCs from the multicenter study by Xing et al. 18 In addition, patients from the studies by Gandolfi et al., 33 Zheng et al., 61 and Alzahrani et al. 62 possibly overlapped with patients from the studies by Sancisi et al., 54 Zheng et al., 60 and Abubaker et al., 27 respectively. We selected studies with higher number of cases for meta-analysis.
Finally, 25 studies comprising 6884 PTC patients were included for this meta-analysis.
RAS mutations
Overall, five studies including 629 patients were included for analysis28,37,45,55,56
RET/PTC rearrangements
There were nine studies comprising 1200 PTCs included for meta-analysis.28,41,45,47,53,55,56,58,63 The population from the two studies47,63 could duplicate with each other, and the study with higher number of cases was selected. The prevalence of
Coexisting mutations
We performed an additional meta-analysis to examine the impact of PTCs with coexisting mutations on DM as compared with PTCs with each mutation alone. There were very limited data on this topic, and we could only find data regarding the association of coexisting mutations with DM in only four studies which provided data on interaction between
Quality of studies
The NOS tool was used to assess the quality of included studies. The number of stars awarded to each study ranged from five to seven stars. Details of given stars within each domain of NOS are described in Table 1.
Subgroup analyses and meta-regression
We performed subgroup analyses according to the region of origin (Caucasian and Asian), time of DM (at presentation and during follow-up), and mutation detection methods.
Subgroup analyses.
OR: odds ratio; CI: confidence interval; NA: not available; DM: distant metastasis; PCR: polymerase chain reaction.
In the meta-regression of the 24 studies on the association between
Publication bias
Funnel plot observation and Egger’s regression test did not show strong evidence of publication bias among the set of studies except for the meta-analysis of
Discussion
DM is a well-known prognostic factor in thyroid cancer, and its presence is associated with an increased risk of mortality. 4 Several studies have investigated the clinicopathological risk factors for DM in DTCs but there are inconsistencies among various studies.14,15,20 Recently, in three decades, novel genetic alterations in thyroid cancer have been discovered in thyroid cancer with the development of translational medicine. It is very essential to investigate the usefulness of these molecular markers to discriminate aggressive tumors from those with an indolent course to further tailor the initial therapeutic decisions and the long-term surveillance as well as to avoid overtreatment. There have been a number of studies reporting the association of molecular biomarkers and DM with PTCs but there are still conflicting results in the literature. As a result, there is a definite need to perform a meta-analysis to clarify the results.
The
Over the last two decades, the association of
The recently discovered genetic event in thyroid cancer,
This study is the first meta-analysis to investigate the use of various genetic markers to predict the most unfavorable prognostic factor in PTCs, DM. Our study can help clinicians and thyroid oncologists to choose reliable molecular markers to assess patient prognosis and select appropriate management. However, there are a few limitations in our study that need to be addressed. The first limitation is that the majority of included studies are retrospective studies so selection bias is inevitable. Another concern is that considerable heterogeneity exists among the studies assessing the association of
In conclusion,
Supplemental Material
Supplemtary_Figures – Supplemental material for Role of molecular markers to predict distant metastasis in papillary thyroid carcinoma: Promising value of TERT promoter mutations and insignificant role of BRAF mutations—a meta-analysis
Supplemental material, Supplemtary_Figures for Role of molecular markers to predict distant metastasis in papillary thyroid carcinoma: Promising value of TERT promoter mutations and insignificant role of BRAF mutations—a meta-analysis by Huy Gia Vuong, Ahmed MA Altibi, Uyen NP Duong, Hanh TT Ngo, Thong Quang Pham, Hung Minh Tran, Naoki Oishi, Kunio Mochizuki, Tadao Nakazawa, Lewis Hassell, Ryohei Katoh and Tetsuo Kondo in Tumor Biology
Footnotes
Acknowledgements
H.G.V., A.M.A.A., and L.H. were responsible for conceptualization; methodology; software; validation; formal analysis; investigation; resources; data curation; original draft preparation, reviewing, and editing; visualization, supervision, and project administration. U.N.P.D. was responsible for methodology, validation, formal analysis, investigation, resources, data curation, and original draft preparation, reviewing, and editing. H.T.T.N., T.Q.P, H.M.T., N.O., K.M., and T.N. were responsible for methodology, resources, data curation, and original draft preparation, reviewing, and editing. R.K. and T.K. were responsible for conceptualization; methodology; software; validation; formal analysis; investigation; resources; data curation; original draft preparation, reviewing, and editing; visualization; supervision; project administration; and funding support.
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.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
