Abstract
Background
Distant metastasis in thyroid cancer significantly reduces survival in patients with well-differentiated thyroid carcinoma (WDTC). There is limited information available to clinicians regarding pathological features that confer a higher risk of distant metastasis (DM). This study aimed to identify patient and tumor factors that were associated with the development of DM over time in patients with WDTC.
Methods
A retrospective cohort analysis of patients with WDTC (n = 584) at our institution was performed between 2007 and 2017. A total of 39 patients with DM and 529 patients with no DM (NDM) were included. Patient demographics, tumor characteristics and patient survival were compared between the DM and NDM groups using a univariate analysis. Multivariate Cox-proportional hazards model was used to evaluate the risk of developing distant metastasis over time. Kaplan-Meier analysis was used to compare survival between the DM and NDM groups.
Results
Distant metastasis had a substantial impact on disease-specific survival (DSS) at 5 and 10-years in the DM group; 71.0% (SE 8.4%) and 46.9% (SE 11.6%) respectively, compared to 100% survival in the NDM group (p < 0.001). The DM group had significantly higher proportions of males, lymphovascular invasion (LVI), nodal metastasis (NM), large tumor size (TS), extrathyroidal extension (ETE), positive resection margins, multifocality, follicular thyroid cancer (FTC), tall cell variant of papillary thyroid cancer (PTC), and Hurthle cell carcinoma (HCC), when compared to the NDM group (p < 0.05). A TS ≥ 2 cm (Hazard Ratio (HR) 1.370), NM (HR 3.806) and FTC (HR 7.068) were associated with a significantly increased hazard of developing distant metastasis in patients with WDTC.
Conclusions
TS ≥ 2 cm, NM and FTC are associated with a significantly increased propensity for developing DM in our cohort of WDTC patients.
Graphical Abstract
Keywords
Background
The incidence of well-differentiated thyroid cancer (WDTC) has increased substantially in recent years, largely owing to a more widespread use of sensitive diagnostic techniques [1]. Improved multidisciplinary treatment regimens have enhanced patient survival rates, estimated at 95% and higher over 10 years [2]. Although distant metastasis (DM) in WDTC remains uncommon, it is associated with a significantly lower disease-specific survival (DSS), estimated at 30–50% (5–10 years), compared to patients with no DM [3–5].
While several studies have reported the prognosticators of survival [5–10], local recurrence and mortality [6, 11, 12] in WDTC, there remains a lack of information regarding the pathological predictors of DM in this population. Specifically, there is limited information regarding the contributions of aggressive pathological subtypes of WDTC [5, 6] in the development of DM, such as non-classical subtypes of papillary thyroid carcinoma (PTC) (tall cell, oncocytic variants), follicular thyroid carcinoma (FTC) and Hurthle cell carcinoma (HCC). The impact of pathological features such as tumor size, lymphovascular invasion, nodal metastasis, multifocality and positive resection margins also remain controversial [13–16]. While some studies have investigated DM in WDTC, there are limited studies utilizing a statistical model that evaluates the development of DM over time for clinically-relevant pathological features.
In this study, we investigated our cohort of WDTC patients to identify patient and pathological factors that are associated with the development of DM over time. Given that PTC is the most common thyroid malignancy with the most favorable clinical outcomes, a separate subgroup analysis of patients with PTC was also conducted. We also report on one of the largest samples of patients with WDTC and DM in the Canadian literature. The overall aim of this study was to better understand the pathological predictors for developing DM in thyroid carcinoma, which will ultimately better identify, and inform the management of, patients with WDTC.
Methods
The interdisciplinary thyroid cancer database (ITCD) is a computerized registry at our institution which prospectively collects information on all patients referred to the thyroid oncology clinic. A thorough search of the ITCD at our institute from 2007 to 2017 identified 584 patients with a diagnosis of thyroid carcinoma, confirmed on post-operative pathology. Patients with medullary thyroid carcinoma and anaplastic thyroid carcinoma were excluded (n = 16). A total of 39 patients with a diagnosis of thyroid carcinoma and distant metastasis (DM) were identified. Patients were included whether DM was diagnosed on initial presentation or throughout the course of treatment. The remaining 529 patients were included as the non-distant metastasis group (NDM). A diagnosis of DM was confirmed with a combination of biopsy of the site of metastasis, positron emission tomography (PET) scan and/or computed tomography (CT), and Iodine-131 scan.
Variables analyzed included age at diagnosis, sex, pathological type of thyroid cancer, histological subtype, primary tumor size (cm), focality, lymphovascular invasion (LVI), margin status, extrathyroidal extension (ETE) and nodal metastasis (NM). Multifocal disease was defined by the presence of ≥3 tumor foci. This was based on findings by Al Afif et al. where ≥3 foci led to a significantly increased propensity for NM [17]. We also compiled data on therapeutic strategies of patients including surgical resection of distant metastasis, the number of radioactive iodine (RAI) treatments, and whether or not external beam radiation therapy (EBRT) was delivered. Cumulative and disease-specific survival were computed for both patient groups, followed by an analysis of variables contributing to the risk of developing DM. A tumor size of 2 cm was associated with increased mortality in previous articles and was therefore used as a cut-off for our hazard model analysis of DM [16, 18, 19].
Statistical analysis
Univariate analysis was performed using Fisher's exact test or Chi Square test for categorical variables, and Mann-Whitney U test for continuous variables. Cumulative and disease-specific survival plots were generated using the Kaplan-Meier method. Differences in survival were tested using log-rank comparisons. The multivariate Cox proportional hazard model was utilized to evaluate the association of variables with the development of DM over time. All statistical analysis was conducted on SPSS (IBM SPSS Statistics 25). A p-value of < 0.05 was considered significant for all tests.
Results
Demographics and tumor pathology
Clinical and pathological characteristics of DM (n = 39) and NDM (n = 529) groups
DM Distant metastasis, NDM Non-distant metastasis, LVI Lymphovascular invasion, ETE Extrathyroidal extension
Distribution of tumor histopathology in DM and NDM patient groups
DM Distant metastasis, NDM Non-distant metastasis
* p < 0.05
Clinical presentation and therapeutic strategies
Clinical features and therapeutic management of patients in the DM group (n = 39) based on initial presentation with DM (n = 12) or development of DM after initial diagnosis (n = 27)
DM Distant metastasis, RAI Radioactive iodine, EBRT External beam radiation therapy
Univariate analysis of patients in the DM group comparing those with DM on initial presentation (n = 12) and those who presented after their diagnosis of thyroid cancer (n = 27)
LVI Lymphovascular invasion
*p < 0.05
Multivariate analysis
Hazard ratios by Cox proportional hazards model for predicting the risk of distant metastasis in WDTC in DM and NDM patients (n = 568)
LVI Lymphovascular invasion, ETE Extrathyroidal extension, CI Confidence interval
*p < 0.05
Subgroup analysis of papillary thyroid Cancer patients
Hazard ratios by Cox proportional hazards model for all DM and NDM patients with PTC (n = 540) to predict the risk of distant metastasis
LVI Lymphovascular invasion, ETE Extrathyroidal extension
*p < 0.05
Survival
DSS and overall survival (OS) at 10 years are shown in Figs. 1 and 2 respectively. DSS at 5- and 10-years for the DM group were 71% (standard error (SE) 8.4%) and 46.9% (SE 11.6%). DSS for the NDM group at 5- and 10 years were 100 and 100% respectively. OS at 5- and 10-years for the DM group were 68.1% (SE 8.6%) and 44.9% (11.3%). OS for the NDM group at 5- and 10 years were 97.2% (SE 1%) and 87.9% (SE 4.4%) respectively.

Kaplan-Meier curves for DSS in the DM (red) and NDM groups (blue) for patients with well-differentiated thyroid cancer

Kaplan-Meier curves for OS in the DM (red) and NDM groups (blue) for patients with well-differentiated thyroid cancer
Discussion
Overall, the body of evidence available to clinicians for assessing a patient's risk for developing DM is limited. A recent meta-analysis of articles using univariate comparisons has identified male sex, advanced age, tumor size, multifocality, LVI, ETE and lymph node metastasis as risk factors for DM in WDTC [20]. Our study demonstrates similar outcomes [20]. Our data shows a significantly higher proportion of male patients, larger tumors, LVI, positive resection margins, multifocality (≥3 foci), ETE and nodal metastasis in the DM group when compared with the NDM group [5, 20]. Additionally, PTC was less prevalent in DM patients, with significantly larger proportions of FTC, when compared with the NDM group. We also investigated the impact of PTC subtype and HCC on the development of DM. We found significantly higher proportions of tall cell variant of PTC and HCC in DM patients when compared to the NDM group. A pathological diagnosis of FTC was also a significant predictor of distant metastasis in our Cox proportional hazard model analysis. Moreover, many patients with tall cell variant of PTC or HCC experienced an aggressive progression of their disease with multiple sites of DM. Our findings suggest that the disproportionally higher rates of tall cell variant of PTC, HCC and FTC in DM patients are clinically important risk factors for the development of DM.
Previous reports have mostly focused on survival prognosticators in patients who have already developed DM or performed univariate comparisons between NDM and DM patients [20]. A statistical model assessing the risk or hazard over time for developing DM for selected pathological variables in a cohort of WDTC patients is less common in the literature. In our study, the pathological variables identified in the univariate analysis were used to better assess their impact on the development of DM. Using this model, a tumor size ≥2 cm (HR = 1.4), nodal metastasis (HR = 3.8) and follicular carcinoma (HR = 7) were the only significant variables associated with an increased chance for developing DM over time. Interestingly, a tumor size ≥2 cm has also been correlated with increased mortality and adverse outcomes in patients with WDTC [16]. The age of patients, positive resection margins, ETE and sex were not associated with a risk of developing DM. Although age has been associated with improved survival in DM patients [10], our analysis did not demonstrate age as a significant factor for developing DM. This can partly be a result of wide age ranges in our dataset. However, as our model investigates the development of DM overtime, younger age may only be a prognosticator of survival once a diagnosis of DM is reached as other studies have suggested [10]. Age was a significant factor in the subgroup analysis of PTC patients alone. The literature remains divided on the predictive value of multifocality in regard to tumor recurrence, metastasis and mortality [12, 21]. The presence of ≥3 foci were not a predictor of DM in our full patient cohort despite showing a propensity for nodal metastasis in previous research [17]. However, the presence of ≥3 foci was significant for an increased likelihood of developing DM in our subgroup analysis of PTC patients alone alongside a tumor size ≥2 cm, nodal metastasis and increased age.
The presence of DM had a substantial impact on survival in our cohort of WDTC patients (5-year and 10-year DSS of 71.0 and 46.9% respectively). Despite a multimodal approach to treatment, including RAI, EBRT and surgical resection of distant tumors, the overall survival remained significantly diminished in DM patients (OS 44.9% at 10 years, p < 0.001). These results emphasize the importance of early detection of patients that exhibit a greater risk of developing DM from WDTC.
The present study has several limitations. It is a retrospective chart review from a single institution. The sample size of the DM group is another limitation of our study. This is partly due to the rarity of DM in WDTC. Multi-center cohort studies are needed to mitigate sample size as a limitation considering the rarity of DM in thyroid carcinoma. Sample sizes required to generate more statistical power should roughly be greater than 68 DM patients according to our sample size estimates for a single-variable cox-proportional hazard model [22]. Given the small sample size, it is also difficult to separately analyze patients who developed DM throughout their treatment, compared to those who presented with DM. However, sparing positive resection margins and ETE, there were no significant differences between the groups. Future studies will focus on differentiating between the two groups. Future research is also investigating the influence of tumor driver mutations on the development of distant metastasis.
Conclusion
WDTC patients with a tumor size ≥2 cm, NM and FTC on final pathology demonstrated a significantly increased propensity for developing DM. A subgroup analysis of PTC patients alone revealed increased age, presence of ≥3 foci, NM and tumor size ≥2 cm as significant contributors for the development of DM. Patients with DM have disproportionally higher frequencies of tall cell subtype of PTC, HCC and FTC. Our findings add to the existing literature on pathological risk factors for developing DM in WDTC, which will ultimately better inform the management of this patient group.
Footnotes
Authors’ contributions
UK participated in data collection, data analysis, statistical analysis and writing of manuscript. AA was involved with data analysis and writing of manuscript. AA was involved in data collection and writing of manuscript. CM and MR performed statistical analysis. MR, JT and MT assisted with the overall design of the study methodology and manuscript editing. MR and AI populated the interdisciplinary thyroid cancer database. MB provided pathological reports. RH participated in study conceptualization, study design, manuscript writing and editing. All authors have approved and read the final manuscript.
Funding
There was no funding obtained for this research study.
Availability of data and materials
Datasets are unavailable for this study as it would compromise patient privacy. However, further information regarding the data are available, within limits of patient privacy, upon request.
Abbreviations
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This paper was presented at the Canadian Society of Otolaryngology-Head and Neck Surgery Conference (CSOHNS) in Edmonton, Alberta, Canada on June 3rd, 2019.
