Risk prediction and clinical model building for lymph node metastasis in papillary thyroid microcarcinoma.

The surgical management of papillary thyroid microcarcinoma (PTMC), especially regarding the necessity of central/lateral lymph node dissection, remains controversial. This study investigated the clinicopathologic factors predictive of lymph node metastasis (LNM) in patients diagnosed with PTMC. Multivariate logistic regression analysis was used for PTMC patients identified from the Surveillance, Epidemiology, and End Results database who were treated by surgery between 2002 and 2012, to determine the association of clinicopathologic factors with LNM. According to the results, a total of 31,017 patients met the inclusion criteria of the study. Final histology confirmed 2,135 (6.9%) cases of N1a disease and 1,684 cases (5.4%) of N1b disease. Our multivariate logistic regression analysis identified variables associated with both central LNM and lateral lymph node metastasis (LLNM), including a younger age (<45 years), male sex, non-Hispanic white and other race, classical papillary histology, larger tumor size, multifocality, and extrathyroidal extension; distant metastasis was also significantly associated with LLNM. The significant predictors identified from multivariable logistic regression were integrated into a statistical model that showed that extrathyroidal extension had maximum weight in the predictive role for LNM. LLNM was validated to be a significant risk factor for cancer-specific death in Cox regression analyses, whereas central LNM failed to predict a worse cancer-specific survival according to our data. Therefore, we suggested that central lymph node dissection could be performed in certain patients with risk factors. Given the prevalence of LLNM in PTMC, a thorough inspection of the lateral compartment is recommended in PTMC patients with risk factors for precise staging; from the viewpoint of a radical treatment for tumors, prophylactic lateral lymph node dissection that aims to remove the occult lateral lymph nodes may be an option for PTMC with risk factors. Multicenter studies with long-term follow-up are recommended to better understand the risk factors and surgical management for cervical nodes in PTMC.

Introduction

In many regions of the world, the incidence of papillary thyroid cancer (PTC) is increasing more rapidly than that of any other malignancy.1 This increase may be explained by the increased detection of papillary thyroid microcarcinoma (PTMC), which is defined as a lesion of ≤10 mm in diameter. With the continuously increasing incidence of PTMC worldwide, the clinical care of PTMC has become a management dilemma, and more importantly, a public health issue.2

The management approach to PTMC cannot be overly cavalier because a significant percentage of patients with microcarcinoma have lymph node metastases (LNMs) at the time of presentation, which implies future recurrences, and a very small number of patients may have, or will subsequently develop, distant metastasis.2–13 On the other hand, given the potentially high prevalence of LNM in PTMC, some authors agree that prophylactic lymph node dissection (LND) permits an improved staging of the disease that may guide subsequent treatment;14 however, few studies have demonstrated a meaningful role of prophylactic LND in relation to locoregional control or survival in PTMC. A systematic review that included 8,345 PTMC patients explored the risk factors associated with central lymph node metastasis (CLNM) and proposed that prophylactic central lymph node dissection (CLND) should be considered in PTMC patients with risk factors.15 However, the deficiency of long-term follow-up data for the analysis was a limitation in that study. Finally, for some groups, it appears reasonable to use a selective approach that applies CLND at the time of initial operation only to patients with clinically evident disease based on preoperative evaluation of primary tumor and nodal disease.

The present research aimed to investigate the frequency, pattern, and predictive factors for LNM in PTMC, and evaluate the optimal surgical treatments for regional nodes in PTMC over a long-term follow-up period, using data from the Surveillance, Epidemiology, and End Results (SEER) Program from 2002 to 2012. Finally, we developed a risk calculator based on risk factors to predict LNM risk, and we discuss the role of prophylactic LND in PTMC patients.

Methods

Selection of patients from the SEER database

For this study, data were extracted from the SEER cancer registry, which is a population-based registry sponsored by the National Cancer Institute which collects information on cancer incidence and survival from 17 population-based cancer registries, including ~28% of the US population.16 The SEER data contain no identifiers and are publicly available for studies on cancer-based epidemiology and health policy. The National Cancer Institute’s SEER*Stat software (Version 8.1.2; Surveillance Research Program, National Cancer Institute, Bethesda, MD, USA; www.seer.cancer.gov/seerstat) was used to identify patients who were diagnosed with single primary PTMC from 2004 to 2012. Patients who had undergone surgery for PTMC as the only malignancy or the first one of more malignancies were included. Histology types were limited to papillary cancer (8,260) and papillary cancer follicular variant (8,340). The patients with insufficient data or unknown clinicopathologic profile, undetermined histology, or other types of thyroid cancers (follicular thyroid cancer, medullary thyroid cancer, anaplastic thyroid cancer, etc) were excluded.

Ethics statement

This study was performed in compliance with the Declaration of Helsinki, approved by an independent ethics committee/institutional review board at Fudan University Shanghai Cancer Center, and conducted in accordance with the approved guidelines. Data released from the SEER database do not require informed patient consent because they contain no identifiers and are publicly available. We obtained permission to access the research data file in the SEER Program from the National Cancer Institute (reference number 13579-Nov 2014).

Assessment of clinicopathologic variables

The following variables were used to analyze the risk factors for CLNM and lateral lymph node metastasis (LLNM): sex, age at diagnosis, race, maximal tumor size, multifocality, and extrathyroidal extension (ETE), surgical procedures, adjuvant therapies, and survival time were extracted from the SEER database. Race was categorized as African American, non-Hispanic white, and others (American Indian/Alaska [AK] Native, Asian/Pacific Islander) as provided by the SEER database.17 According to the current staging system, an age of 45 years was used as the cutoff point to divide all patients into two groups: younger patients (<45 years) and older patients (≥45 years).18 Multifocal primary lesions were defined as two or more cancer sites within the thyroid; a solitary intrathyroidal primary lesion represented only a single cancer site within the thyroid. An ETE of the primary tumor consisted of thyroid capsule, minimal ETE (strap muscle), and gross ETE (nerves, esophagus, larynx, sternocleidomastoid muscle, etc). The multifocal lesion and ETE were both assessed in the final pathology. The end point of the present study was cancer-specific survival (CSS) of PTMC which was calculated from the day of diagnosis to the day of cancer-specific death and was shown as “SEER cause-specific survival” in the SEER database. We followed the guidance of the 2010 TNM classification of American Joint Committee on Cancer.17,18

Statistical analysis

The results are expressed as n (%) or mean ± standard deviation. The statistical analyses were performed using Student’s t-test, the χ2 test, or the Mann–Whitney test as appropriate. The odds ratio (OR) and the 95% confidence interval (CI) for the relationships between each variable and LNM (yes or no) were calculated using binary logistic regression. A final model selection was performed by a backward stepdown selection process with the Akaike information criterion.19 The discrimination and calibration of the simplified scoring system were then assessed and internally validated using bootstrapping techniques. A Cox proportional hazards regression model was built to evaluate the effects of variables on cancer mortality in PTMC patients and then adjusted by multivariate analysis for patient characteristics and tumor characteristics. The hazard ratios for relationships between each variable and CSS were calculated using a binary Cox regression model. All statistical analyses were performed using package of rms20 in R version 3.2.2 (The R Foundation, Vienna, Austria; http://www.r-project.org/). Two-sided P<0.05 was considered statistically significant.

Results

Baseline characteristics

We identified 75,791 eligible patients with PTC from the SEER database treated between 2002 and 2012; of these patients, 31,017 (40.9%) were diagnosed with PTMC histologically. The age of the enrolled patients ranged from 3 to 99 years (mean 50.9±13.9). Regarding surgical procedures, 23,572 (76.0%) patients underwent total thyroidectomy (TT). Of the patients who underwent TT, radioisotope therapy was performed in 8,886 cases (28.6%) as an adjuvant therapy. Among the patients who underwent LND, the final pathological results revealed that 2,135 patients (6.9%) had positive malignant LNM in central neck, while 1,684 (5.4%) had it in lateral neck. The clinicopathologic characteristics of all patients are summarized in Table 1.

Table 1

Characteristics of patients with PTMC with LNM identified from the SEER database

Variable	Total	CLNM-negative	CLNM-positive	P-value	LLNM-negative	LLNM-positive	P-value
Total	31,017 (100.0)	28,882 (93.1)	2,135 (6.9)		29,333 (94.6)	1,684 (5.4)
Sex
Male	5,725 (18.5)	5,201 (18.0)	524 (24.5)	0.001	5,117 (17.4)	608 (36.1)	0.001
Female	25,292 (81.5)	23,681 (82.0)	1,611 (75.5)		24,216 (82.6)	1,076 (63.9)
Age (years)
<45	10,272 (33.1)	9,362 (32.4)	910 (42.6)	0.001	9,485 (32.3)	787 (46.7)	0.001
≥45	20,745 (66.9)	19,520 (67.6)	1,225 (57.4)		19,848 (67.7)	897 (53.3)
Race
Black	1,907 (6.1)	1,837 (6.4)	70 (3.3)	0.001	1,872 (6.4)	35 (2.1)	0.001
White	25,958 (83.7)	24,139 (83.6)	1,819 (85.2)		24,518 (83.6)	1,440 (85.5)
Othera	2,784 (9.0)	2,573 (8.9)	211 (9.9)		2,597 (8.9)	187 (11.1)
Size (mm)	5.4±3.0	5.3±3.0	6.9±2.7	0.001	5.4±3.0	6.7±2.7	0.001
Histology
FV-PTMC	22,482 (72.5)	20,771 (71.9)	1,711 (80.1)	0.001	21,113 (72.0)	1,369 (81.3)	0.001
Papillary	8,535 (27.5)	8,111 (28.1)	424 (19.9)		8,220 (28.0)	315 (18.7)
Multifocality	10,751 (34.7)	9,664 (35.5)	1,087 (50.9)	0.001	9,763 (33.3)	988 (58.7)	0.001
ETE
Capsule	1,000 (3.2)	910 (3.2)	90 (4.2)	0.001	904 (3.1)	96 (5.7)	0.001
Minimal	1,643 (5.3)	1,330 (4.6)	313 (14.7)		1,297 (4.4)	346 (20.5)
Gross	161 (0.5)	135 (0.5)	26 (1.2)		112 (0.4)	49 (2.9)
Metastasisb	92 (0.3)	81 (0.3)	11 (0.5)	0.001	46 (0.2)	46 (2.7)	0.001
TT	23,572 (76.0)	21,661 (75.0)	1,911 (89.5)		21,979 (74.9)	1,593 (94.6)	0.001
Radioisotopes	8,886 (28.6)	7,585 (26.3)	1,301 (60.9)	0.001	7,688 (26.2)	1,198 (71.1)	0.001

Notes: Data are presented as n (%) or mean ± standard deviation.

Including American Indian/Alaska Native and Asian/Pacific Islander.

Including the lungs, bones, and brains.

Abbreviations: PTMC, papillary thyroid microcarcinoma; SEER, Surveillance, Epidemiology, and End Results; LNM, lymph node metastasis; CLNM, central lymph node metastasis; LLNM, lateral lymph node metastasis; FV-PTMC, follicular variant of papillary thyroid microcarcinoma; ETE, extrathyroidal extension; TT, total thyroidectomy.

Risk factors for CLNM

The primary tumor size as a continuous variable was significantly different between the CLNM-negative and CLNM-positive groups (P=0.001, Table 1). Tumor size was associated with CLNM according to a univariate logistic regression (OR 1.199, 95% CI 1.181–1.219, P=0.001). Other variables, such as male sex (OR 0.675, 95% CI 0.609–0.749, P=0.001), younger age (<45 years, OR 0.646, 95% CI 0.591–0.706, P=0.001), non-Hispanic white (OR 1.978, 95% CI 1.550–2.522, P=0.001) and other (American Indian/AK Native, Asian/Pacific Islander) race (OR 2.152, 95% CI 1.632–2.839, P=0.001), classical papillary histology (OR 0.635, 95% CI 0.569–0.708, P=0.001), multifocality (OR 2.063, 95% CI 1.888–2.253, P=0.001), and ETE (capsule, OR 1.537, 95% CI 1.232–1.919, P=0.001; minimal ETE, OR 2.992, 95% CI 1.961–4.566, P=0.001; gross ETE, OR 3.657, 95% CI 3.203–4.175, P=0.001) were also found to be significantly associated with CLNM via the univariate logistic regression. Distant metastasis (OR 1.841, 95% CI 0.979–3.462, P=0.058) was not associated with CLNM. When these variables were included in the multivariate logistic regression models, they were all found to be significant predictors of CLNM except distant metastasis (Table 2).

Table 2

Univariate and multivariate logistic regression analysis for CLNM in patients with PTMC

Independent variable	Univariate		Multivariate
	OR (95% CI)	P-value	OR (95% CI)	P-value
Sex (female vs male)	0.675 (0.609–0.749)	0.001	0.673 (0.605–0.748)	0.001
Age (≥45 vs <45 years)	0.646 (0.591–0.706)	0.001	0.683 (0.623–0.749)	0.001
Race
Black	1 (reference)		1 (reference)
White	1.978 (1.550–2.522)	0.001	1.656 (1.295–2.117)	0.001
Othera	2.152 (1.632–2.839)	0.001	1.578 (1.191–2.089)	0.001
Size	1.199 (1.181–1.219)	0.001	1.172 (1.153–1.191)	0.001
Histology (PTMC vs FV-PTMC)	0.635 (0.569–0.708)	0.001	0.641 (0.574–0.716)	0.001
Multifocality	2.063 (1.888–2.253)	0.001	1.662 (1.516–1.821)	0.001
ETE
Negative	1 (reference)		1 (reference)
Capsule	1.537 (1.232–1.919)	0.001	1.232 (1.183–1.543)	0.007
Minimal	2.992 (1.961–4.566)	0.001	2.119 (1.377–3.263)	0.001
Gross	3.657 (3.203–4.175)	0.001	2.546 (2.218–2.921)	0.001
Metastasisb	1.841 (0.979–3.462)	0.058	1.152 (0.598–2.222)	0.672

Notes:

Including American Indian/Alaska Native and Asian/Pacific Islander.

Including the lungs, bones, and brains.

Abbreviations: CLNM, central lymph node metastasis; PTMC, papillary thyroid microcarcinoma; OR, odds ratio; CI, confidence interval; FV-PTMC, follicular variant of papillary thyroid microcarcinoma; ETE, extrathyroidal extension.

Risk factors for LLNM

The risk factors for LLNM were evaluated using univariate and multivariate logistic regression analyses for the patients who had positive LLNM according to the final pathological results. An increase in the primary tumor size was associated with a corresponding increase in the risk of LLNM (OR 1.112, 95% CI 1.091–1.133, P=0.001). With the exception of the primary tumor size, the univariate analyses also indicated that male sex (OR 0.374, 95% CI 0.337–0.415, P=0.001), younger age (<45 years, OR 0.545, 95% CI 0.493–0.601, P=0.001), non-Hispanic white (OR 3.141, 95% CI 2.239–4.407, P=0.001) and other (American Indian/AK Native, Asian/Pacific Islander) race (OR 3.851, 95% CI 1.971–5.868, P=0.001), classical papillary histology (OR 0.591, 95% CI 0.521–0.670, P=0.001), multifocality (OR 2.845, 95% CI 2.575–3.145, P=0.001), and ETE (capsule, OR 2.405, 95% CI 1.934–2.992, P=0.001; minimal ETE, OR 6.042, 95% CI 5.295–6.895, P=0.001; gross ETE, OR 9.909, 95% CI 7.048–13.931, P=0.001) were also found to be significantly associated with LLNM via the univariate logistic regression. In particular, distant metastasis was a strong predictor of LLNM (OR 17.880, 95% CI 11.846–26.987, P=0.001). When these variables were included in the multivariate logistic regression models, they were all found to be associated with LLNM (Table 3).

Table 3

Univariate and multivariate logistic regression analysis for LLNM in patients with PTMC

Independent variable	Univariate		Multivariate
	OR (95% CI)	P-value	OR (95% CI)	P-value
Sex (female vs male)	0.374 (0.337–0.415)	0.001	0.360 (0.322–0.401)	0.001
Age (≥45 vs <45 years)	0.545 (0.493–0.601)	0.001	0.503 (0.453–0.559)	0.001
Race
Black	1 (reference)		1 (reference)
White	3.141 (2.239–4.407)	0.001	2.441 (1.730–3.443)	0.001
Othera	3.851 (1.971–5.868)	0.001	2.620 (1.802–3.808)	0.001
Size	1.170 (1.150–1.191)	0.001	1.112 (1.091–1.133)	0.001
Histology (PTMC vs FV-PTMC)	0.591 (0.521–0.670)	0.001	0.625 (0.549–0.711)	0.001
Multifocality	2.845 (2.575–3.145)	0.001	2.265 (2.039–2.517)	0.001
ETE
Negative	1 (reference)		1 (reference)
Capsule	2.405 (1.934–2.992)	0.001	2.105 (1.681–2.637)	0.001
Minimal	6.042 (5.295–6.895)	0.001	4.601 (3.994–5.300)	0.001
Gross	9.909 (7.048–13.931)	0.001	7.093 (4.910–10.246)	0.001
Metastasisb	17.880 (11.846–26.987)	0.001	11.948 (7.523–18.975)	0.001

Notes:

Including American Indian/Alaska Native and Asian/Pacific Islander.

Including the lungs, bones, and brains.

Abbreviations: LLNM, lateral lymph node metastasis; PTMC, papillary thyroid microcarcinoma; OR, odds ratio; CI, confidence interval; FV-PTMC, follicular variant of papillary thyroid microcarcinoma; ETE, extrathyroidal extension.

Prognostic model for cervical LNM

The final method of displaying the relationships among tumor size, other predictors, and probability of LNM is to construct a model. A predictive model not only sheds light on the effect of each predictor on the probability of response by levels of other factors but also allows quick estimation of probability of response for individual subjects.21 For convenience, we define the response to be Y=0 or 1, with Y=1 denoting positive development of LNM. The significant predictors identified from multivariable logistic regression were integrated into a statistical model that shows the relationship between tumor size and LNM occurrence in PTMC patients (Figure S1).21 The values of the predictors derived from the coefficients of logistic regression, and the regression equation for LNM risk are given by the following equations: where

=−2.58+0.14 race −0.39 sex −0.47 histology +0.15 size +0.18 capsule +0.86 minimal ETE +0.66 gross ETE +0.51 multifocality −0.38 age where

= −1.58+0.21 race, −1.04 sex, −0.51 histology +0.09 size, 0.13+ distant metastasis, +0.70 capsule, +1.42 minimal ETE, +1.96 gross ETE, +0.81 multifocality, −0.67 age distant metastasis, ETE, multifocality, >45 yr =1, 0 otherwise; race[black] =1, [white] =2, [other] =3; histology[PTC] =1, [FV-PTC] =2.

The output of validation indicates minor overfitting. Overfitting would have been worse had the risk factors not been so strong. The closeness of the calibration curve to the 45° line in Figure S2 demonstrates excellent validation on an absolute probability scale. However, the extent of missing data casts some doubt on the validity of this model, and on the efficiency of its parameter estimates.21

Predictors of CSS in patients with PTMC

The mean length of follow-up was 43.5±29.9 months and ranged from 0 to 107 months. During the follow-up period, only 63 (0.2%) died as a result of a cancer-specific cause. Male sex, elder age (≥45 years), gross ETE, LLNM, and distant metastasis were validated to be significant risk factors for cancer-specific death in both the univariate and multivariate Cox regression analyses. The patients with minimal ETE exhibited a significantly worse CSS compared with the patients without ETE; however, minimal ETE was not a significant predictor of cancer-specific death in the multivariate Cox regression analysis. Other factors, such as race, tumor size, histology, and multifocality, were not associated with CSS in univariate or multivariate analyses. In particular, TT or postoperative radioisotopes therapy did not improve the CSS significantly in patients with PTMC (Table 4).

Table 4

Univariate and multivariate Cox regression analysis of CSS in patients with PTMC

Independent variable	Univariate		Multivariate
	HR (95% CI)	P-value	HR (95% CI)	P-value
Sex (female vs male)	0.281 (0.169–0.465)	0.001	0.544 (0.321–0.923)	0.024
Age (≥45 vs <45 years)	5.709 (2.287–14.253)	0.001	7.419 (2.898–18.990)	0.001
Race
Black	1 (reference)		1 (reference)
White	3.749 (0.519–27.101)	0.190	2.052 (0.281–15.007)	0.479
Othera	4.117 (0.496–34.199)	0.190	2.319 (0.275–19.548)	0.439
Size	1.054 (0.969–1.147)	0.222	1.001 (0.914–1.096)	0.989
Histology (FV-PTMC vs PTMC)	0.721 (0.391–1.331)	0.295	0.845 (0.452–1.577)	0.596
Multifocality	1.608 (0.972–2.661)	0.065	0.938 (0.543–1.621)	0.819
ETE
Negative	1 (reference)		1 (reference)
Capsule	0.759 (0.104–5.513)	0.785	0.625 (0.085–4.564)	0.643
Minimal	3.262 (1.531–6.948)	0.002	1.358 (0.604–3.053)	0.458
Gross	38.880 (19.504–77.508)	0.001	5.844 (2.617–13.051)	0.001
CLNM	1.840 (0.837–4.043)	0.129	2.345 (0.989–5.563)	0.053
LLNM	10.004 (5.960–16.792)	0.001	4.881 (2.496–9.544)	0.001
Metastasisb	85.311 (45.373–160.404)	0.001	18.051 (8.732–37.314)	0.001
TT	1.318 (0.701–2.479)	0.391	0.794 (0.396–1.595)	0.518
Radioisotopes	2.637 (0.535–4.529)	0.068	1.563 (0.841–2.904)	0.158

Notes:

Including American Indian/Alaska Native and Asian/Pacific Islander.

Including the lungs, bones, and brains.

Abbreviations: CSS, cancer-specific survival; PTMC, papillary thyroid microcarcinoma; HR, hazard ratio; CI, confidence interval; FV-PTMC, follicular variant of papillary thyroid microcarcinoma; ETE, extrathyroidal extension; CLNM, central lymph node metastasis; LLNM, lateral lymph node metastasis; TT, total thyroidectomy.

Discussion

The incidence of PTC has increased rapidly in the past years; small PTCs account for most new diagnoses, which may result from the increased detection of small tumors that are incidentally noted on diagnostic imaging studies.1,22 Recently, the increased use of imaging procedures, such as ultrasonography of the neck, has identified a large number of small thyroid nodules, with a similar increasing trend in the proportion of PTMC in all PTC patients.

There were 31,017 patients (40.9%) histologically diagnosed with PTMC based on the SEER data from 2002 to 2012. The prevalence rates of CLNM and LLNM in present study were 6.9% and 5.4%, respectively, values that were lower compared to previous studies that reported a 17.8%–64.1% rate for CLNM and a 3.7%–44.5% rate for LLNM.4,8,23–30 The variation in the prevalence of cervical LNM across different studies may result from the differences in the number of patients, the different rationales (therapeutic or prophylactic) underlying the use of LND in PTMC, and the subtle variations in the surgical completeness among surgeons and the number of lymph nodes in histopathological findings detected among pathologists.15

In our study, multivariate logistic analyses were used to evaluate the risk factors for LNM. We determined that male sex, younger age (<45 years), non-African American race, increased tumor size, classical papillary histology, multifocality, and ETE were significantly associated with LNM; distant metastasis was also significantly associated with LLNM.

Male sex was a risk factor for LNM, which suggests that males require specialized neck checkups of the thyroid to enable the early detection of thyroid tumors.15 Patient age is known to be an important risk factor for patients with PTC >1 cm;31,32 our results indicated that younger age (<45 years) was associated with an increased risk of LNM in PTMC patients. In particular, children and adolescents (<20 years of age) are more likely to have more advanced nodal stage at the time of diagnosis. One large study found that 64% of children had cervical LNMs and 23% had distant metastases at the time of diagnosis. In contrast, studies of adults have found that up to 40% have LNMs and only 5% present with distant metastases.33 Despite the extent of disease at the time of diagnosis, younger patients generally have excellent survival rates. One large study found a 2% cause-specific mortality rate after 40 years of follow-up.33,34 This finding may be explained by the results of our survival analysis or current staging and risk-stratification systems. First, younger age was not associated with CSS according to our data; moreover, even when younger patients showed nodal involvement, they were still staged as I, and their risk of mortality did not increase. On the other hand, most authorities suggest that younger patients have been treated with TT and radioiodine, and rigorous initial surgical and radioiodine treatment followed by thyrotropin suppression has resulted in favorable outcomes. Therefore, younger age (<45 years) had a greater impact on LNM than on prognosis in PTMC patients.

A greater tumor size was significantly associated with LNM, which was consistent with a previous study.31 As noted, the evaluation of LLNM was not simply dependent on the primary tumor size but, more importantly, on the results of the preoperative imaging examination, especially in the PTMC patients with increased tumor size.

Interestingly, our data showed that the presence of LNM in papillary cancer follicular variant at diagnosis was unlikely to be PTC. This sub-variant histology has been consistently shown in many studies to be associated with increased mortality and poorer outcome; its wide invasiveness is probably owing to the distinct biological features of Hurthle cell carcinomas from PTC, such as the demonstration of vascular or capsular invasion and the metastatic spread to the blood.

Multifocal lesion and ETEs were also significant factors for LNM. A recent study suggested that multifocal lesion represented multiple lesions with individual characteristics;35 thus, the characteristics of the individual tumor, rather than multifocal lesion itself, tend to correlate more with the risk of LNM. Although ETE was not a common finding in PTMC, it could be well predicted with a preoperative, high-resolution ultrasound examination.36–38 A preoperative imaging examination helps to predict cervical LNM in PTMC patients. On the other hand, it is also important to differentiate the clinical significance of minimal ETE (American Thyroid Association [ATA] intermediate risk) from gross ETE of surrounding structures (ATA high risk) by pathology.28 The risk of recurrence associated with minimal ETE (pT3 disease manifest by ETE) ranges from 3% to 9%, while the risk of recurrence in patients with gross ETE (pT4a disease involving the subcutaneous soft tissues, larynx, trachea, esophagus, or recurrent laryngeal nerve) ranges from 23% to 40%.39–46

Based on our finding of the association between risk factors and LNM, we further designed a model that systematically considers multiple variables to estimate an individual risk of LNM in patients with PTMC. Our model employs easily accessible clinical information, with high concordance on an absolute probability scale. Coefficients from the model indicate a change in the log(OR) per unit increment in tumor size (1 mm, as continuous variable) or presenting with risk factors (as binary variables), when they represent a single factor that is linear and does not interact with other factors and if all other factors are held constant. For either CLNM or LLNM, ETE had maximum weight in the predictive model. As a form of tumor invasion, ETE was consistent with the previously described risk factors in the criteria for the assessment of progression. On the other hand, the ATA developed consensus-based criteria for the assessment of recurrence risk, and this system has been validated retrospectively with clinical data.47 PTMC was divided into low-, intermediate-, and high-risk groups, which represented different risks of poor prognosis. The basis of the classification included risk factors such as local metastases (LNM), residual tumor, tumor invasion, and aggressive histology. Therefore, ETE and other risk factors listed by the ATA criteria might provide guidance in the initial management of PTMC patients.

The optimal management of PTMC is a controversial issue. Given the high incidence of multifocal lesion, extensive TT surgery was considered to be more likely to remove all disease foci and improve the prognosis. However, TT did not improve the prognosis in the present study. Namely, the evidence from the present study may favor lobectomy for most PTMC patients without risk factors because of the comparative impact of the prognosis and the lower rate of surgical morbidity reported compared with TT.48 As noted, the presence of distant metastases was a strong risk factor for cancer-specific death. TT and postoperative radioisotopes therapy were an important approach to detect positive distant metastases at diagnosis or recurrence during follow-up by imaging or whole-body scan. However, given the small number of patients who were presenting with radioactive iodine at diagnosis (0.3%), it might be limited to perform postoperative radioisotopes therapy in all PTMC patients. Further studies are needed to explore the indication for RAI therapy in PTMC patients.

In the current guidelines,47,49,50 prophylactic CLND is not recommended for patients with clinically node-negative (cN0) PTMC. However, some authors have argued its application in cN0 PTMC because of the high prevalence of CLNM.25,27,35,51 The removal of the occult nodes in the central compartment may confirm the status of the central lymph node for predicting LLNM and guiding treatment, with the aim of reducing the disease persistence.52–55 Although CLNM was not significantly associated with prognosis according to survival analysis, from the viewpoint of the radical treatment of tumors, prophylactic CLND must be considered in PTMC patients who present with high-risk factors.

The rationale for prophylactic LLND in PTC is controversial. The relation between prophylactic LLND and locoregional control or survival has not been demonstrated in previous studies. In the present study, cervical LNM was a prognostic indicator in the PTMC patients, and similar results have also been demonstrated in several previous studies.56 However, the associations between the node status and the recurrence prognosis with a long-term follow-up period and in a larger PTMC population remain unknown. The majority of reports on the predictive factors associated with a high risk of node metastasis have included clinically lymph node-positive patients and have not addressed the issue of occult lateral lymph nodes. Lim et al57 conducted the first prospective quantitative evaluation of LLNM in PTMC without clinical LLNM. In their study, the incidence of occult LNM was as high as 55%, and the primary tumor size and the CLNM were two important risk factors for LLNM. Therefore, it appears that prophylactic LLND may be recommended in selected patients with high-risk factors.

Our study has some limitations. First, the confirmation for LNM was mainly developed from the collection of therapeutic LNDs. The low incidence of LNM in this cohort suggests the possibility of underestimation of LNM in patients with PTMC. However, previous studies have shown that ultrasound is an important tool in the detection of metastatic nodes, with a high specificity and a positive predictive value. Therefore, the current ultrasound system is relatively diagnostically accurate for clinical LNM, and seems to be improved by our findings or the model from this study, which integrates a series of clinicopathologic variables. Second, in patients with LLNM from the SEER database, we could not differentiate the status of the central lymph node. Although CLNM is reported to predict LLNM, our risk evaluation for LLNM did not consider CLNM as a risk factor in the multivariate analysis. Diagnostic algorithms for LNM suggested by this analysis have not been tested and will require prospective study before recommendations can be made. We also recommend prospective multicenter studies with long-term follow-ups to ascertain the optimal surgical approach for patients with PTMC.

In conclusion, the results of our study demonstrate that male sex, younger age (<45 years), white and other race (including American Indian/AK Native, Asian/Pacific Islander), a greater tumor size, ETE, and multifocality are risk factors for cervical LNM. Distant metastasis is also a predictor of LLNM. CLND is recommended in select patients with risk factors. Surgeons could identify some of the high-risk factors or features which are associated with LLNM preoperatively by ultrasound, imaging exams, fine-needle biopsy, or molecular analysis. Given the prevalence of LLNM in PTMC, a thorough inspection of the lateral compartment is recommended in PTMC patients with risk factors; from the viewpoint of radical treatment for tumors, prophylactic LLND that aims to remove occult lateral lymph nodes may be an option for PTMC with high-risk factors. We recommend prospective multicenter studies with long-term follow-ups to ascertain the optimal surgical approach for patients with PTMC.

Lymph node metastasis predictive model.

Notes: A model selection was performed by a backward stepdown selection process with the Akaike information criterion. (A) Central lymph node metastasis. (B) Lateral lymph node metastasis.

Abbreviation: ETE, extrathyroidal extension.

The validation for the model by the calibration curve.

Notes: (A) Central lymph node metastasis. (B) Lateral lymph node metastasis.

Abbreviations: CLNM, central lymph node metastasis; LLNM, lateral lymph node metastasis.