Prediction Model of Pathologic Central Lymph Node Negativity in cN0 Papillary Thyroid Carcinoma

Background: Most patients with papillary thyroid carcinoma (PTC) have an excellent prognosis. Although central lymph node invasion is frequent, management via central lymph node dissection (CLND) remains controversial. The present study retrospectively investigated independent predictors of pathologic central lymph node negativity (pCLN-) and established a prediction model for pCLN- in clinical lymph node negativity (cN0) PTC.

Methods: A total of 2,687 patients underwent thyroid surgery for cN0 PTC from 2013 to 2018 at the First Affiliated Hospital of Chongqing Medical University, and lobectomy plus ipsilateral CLND was the basic surgical extent. Clinicopathological characteristics were reviewed and analyzed. Univariate and multivariate analyses were performed to identify factors related to pCLN-. A prediction model was established based on the results of multivariate analyses.

Results: The pCLN- rate was 51.5% (1,383/2,687). Multivariate analysis revealed that sex, age, thyroid stimulating hormone (TSH), size, location, laterality, unifocality and extrathyroidal extension negativity (ETE-) were independent predictors of pCLN-. The nomogram showed good discriminative ability (C-index: 0.784 and 0.787 in derivation and validation groups, respectively) and was well calibrated. We quantified the clinical usefulness of the nomogram by decision curve analysis. The median length of follow-up was 30 (range 12– 83) months, and 190 cases were lost, with a follow-up rate of 92.9% (2,497/2,687). Of the 2,687 patients included, 21 (0.8%) experienced recurrence.

Conclusion: This nomogram, which integrates available preoperative clinicopathological features and intraoperative frozen biopsy outcomes, is a reliable tool with high accuracy to predict pCLN- in cN0 PTC.

Introduction

Thyroid cancer is the most common endocrine malignancy, and papillary thyroid carcinoma (PTC) accounts for more than 90% of cases. The development of ultrasound and fine needle aspiration cytology rapidly improved the diagnosis of PTC in recent years (1–3). Most of PTCs are indolent and carry an excellent prognosis with low mortality. However, early-stage metastasis in cervical lymph nodes may occur, especially in the central compartment (Level VI), which is bound superiorly by the hyoid bone, inferiorly by the innominate artery, and laterally by the common carotid artery. Overall, accurate evaluation of lymph node status is of great significance for PTC treatment, TNM staging and risk stratification, but current methods for evaluating central lymph node (CLN) status have several limitations. Inconsistencies between clinical and pathologic evaluations may also exist. According to the 8th edition of the American Joint Committee on Cancer (AJCC) TNM staging, N0 is subdivided into N0a and N0b. N0a, also known as pN0, refers to negative lymph nodes confirmed by cytology or histology, and N0b, also known as cN0, indicates no metastatic lymph nodes (normal lymph nodes or nonspecific lymphadenitis) found on clinical and imaging examinations (4–7). Ultrasound is the most commonly used tool for the clinical evaluation of lymph nodes. The ultrasonic signs of typical lymph node metastasis primarily include a round shape, disappearance of the lymphatic hilum, mass hyperechoic regions, calcification foci (primarily microcalcification; coarse calcification is rare), cystic changes, and peripheral or mixed blood flow distribution signals. None of these features alone is sufficient to diagnose all metastatic lymph nodes, and the accuracy of the evaluation is limited based on considering the experience of clinicians, the function of instruments and heterogeneity among patients, the accuracy of the evaluation is limited (8–10). Therefore, only some patients with cN0 PTC have occult central lymph node metastasis (CLNM). The 2015 American Thyroid Association (ATA) and 2021 National Comprehensive Cancer Network (NCCN) guidelines recommend that T1 and T2 cN0 PTC should not be treated using prophylactic central lymph node dissection (CLND). However, surgical management of cN0 PTC remains controversial because of the presence of occult metastatic lymph nodes in cN0 PTC and the debate surrounding the effect of these lymph nodes on recurrence and prognosis (5, 11–14). Balancing the risks and benefits of lymph node dissection is crucial for individualized and precise treatment. For cN0 PTC with benign CLNs, it is reasonable to avoid prophylactic CLND, which protects the cervical tissue, reduces the extent of surgery and the corresponding complications, and facilitates the recovery of patients after surgery. This approach essentially promotes the development of scarless endoscopic thyroidectomy (SET). However, there have been few studies on the clinicopathological characteristics of these patients and a practical quantitative tool is urgently needed to identify these patients. We aim to establish a prediction model of pathologic central lymph node negativity (pCLN-) in patients with cN0 PTC.

Materials and Methods

Patient Collection

The local institutional ethics committee approved this retrospective study. Because there was no patient interest or privacy, informed consent was waived after careful review by the ethics committee. A total of 3,916 patients underwent surgery for thyroid diseases at the Department of Endocrine and Breast Surgery, The First Affiliated Hospital of Chongqing Medical University, between January 2013 and December 2018. All patients underwent at least a physical examination and cervical ultrasound examination before surgery. We used the cN0 assessment criteria of Kouvaraki et al. (1). Physical examination did not address enlarged lymph nodes or swollen lymph nodes that had a soft texture (2). Ultrasound examination showed no enlarged lymph nodes or swollen lymph nodes were oval and flat, regularly shaped, with clear boundaries between the cortex and medulla or a clear central fat hilum and no obvious malignant signs (15). Patients with any of the following conditions were excluded: incomplete clinicopathological data (403), CN1 confirmed by preoperative examination (536), non-PTC histology (34), reoperation or other head and neck surgery history (105), drug therapy for hyperthyroidism or hypothyroidism (27), absence of CLND (79), distant metastasis (14), adolescent (22), radiation exposure or family history of PTC (9). A total of 2,687 eligible patients were ultimately included in this study (Figure 1).

FIGURE 1

Figure 1 Flowchart of this retrospective study.

Surgical Procedure

For all patients enrolled, lobectomy plus ipsilateral CLND was the basic surgical treatment. Total thyroidectomy plus bilateral CLND was further performed for cases with bilateral tumors, tumors with a maximum diameter > 2 cm or intraoperative frozen section biopsy revealing CLNM. We routinely performed CLN intraoperative frozen section biopsy. Three pathologists blindly and independently diagnosed all resected thyroid and lymph node specimens.

Demographic and Clinicopathological Data

Data for sex, age, body mass index (BMI), thyroid stimulating hormone (TSH), size, location, unifocality, laterality, Hashimoto’s thyroiditis (HT) and extrathyroidal extension negativity (ETE-) were collected and analyzed. The diameter of the largest lesion was recorded as the size, and the central location of the largest lesion was recorded as the location. When the tumor occupied almost the entire lobe, its location was recorded as “whole”. Notably, ipsilateral lesions indicated the presence of laterality, and the presence of only one lesion in one lobe indicated unifocality. Tumors in which growth was confined to the thyroid gland without capsule invasion were defined as ETE-, minimal ETE confirmed by intraoperative frozen section biopsy was not considered ETE-.

Follow-up Monitoring

For all patients, the first follow-up was 1 month after surgery, and evaluations included a physical examination and measurement of serum TSH levels. Measurements of serum thyroglobulin (Tg) and Tg antibody levels were only performed for patients who underwent total thyroidectomy. The second follow-up was 6 months postoperatively. Cervical ultrasound examination was also performed for the first time. Regular evaluations were performed every 6-12 months thereafter. Patients with clinical evidence of recurrence, an unstimulated Tg level ≥2 ng/mL, a stimulated Tg level ≥ 20 ng/mL or serial increases in serum Tg (with negative Tg antibody), or an imaging abnormality were suspected of recurrence (16, 17). We performed fine needle aspiration in these patients and reoperations to identify the pathology of the recurrence.

Statistical Analysis

The Mann-Whitney test was performed to compare differences in nonnormally distributed continuous variables (e.g., age, BMI, TSH, and size) between the two groups. Continuous variables with significant differences were converted into categorical variables. The chi-square test or Fisher’s exact test was used to compare categorical variables. The Excel random function was used to assign values to each patient, and the patients corresponding to these values were randomly sorted. Therefore, a total of 2,687 patients were randomly divided into a derivation group (70%) and a validation group (30%). Multivariate analyses were performed using logistic regression analysis to estimate the statistical significance of the relationships between pCLN- and the clinicopathological characteristics in the derivation group. P < 0.05 was considered a statistically significant difference. A nomogram was developed based on the results of these analyses, and the predictive ability of the nomogram was measured using the receiver operating characteristic (ROC) curve, concordance index (C-index) and calibration plot. To evaluate the clinical usefulness of the nomogram, decision curve analysis (DCA) was performed to calculate the net benefits at different threshold probabilities. The Kaplan-Meier method was used to estimate disease-free survival(DFS) rates. All statistical analyses were performed using SPSS version 25.0 (IBM Corp., USA), and graphs were generated using R software.

Results

Patient Characteristics

A total of 2,687 patients were included in the analysis (684 males and 2,003 females). The median age at diagnosis was 43 years (range 18-85 years). The median BMI level was 23.07 (range 15.56-36.98), and the mean serum concentration of TSH was 2.23 μIU/mL (range 0.01-97.18 μIU/mL). The median tumor size was 9 mm (range 3-64 mm). The numbers of tumors located in the upper portion, middle portion, lower portion, isthmus and whole thyroid were 592 (22.0%), 1,133 (42.2%), 806 (30.0%), 85 (3.2%) and 71 (2.6%), respectively. Right, left and bilateral tumors were detected in 1,211 (45.1%), 1,059 (39.4%), and 417 (15.5%) patients, respectively. Unifocality and ETE- were detected in 2,364 (88.0%) and 2,420 (90.1%) cases, respectively. Histological HT was present in 531 patients (19.8%) (Table 1).

TABLE 1

Table 1 Clinicopathological characteristics of all patients (derivation + validation groups).

Lymph Node Yield and Lymph Node Ratio

Overall, the rate of pCLN- was 51.5%. The median central lymph node yield (CLNY), central lymph node ratio (CLNR), and the number of ipsilateral CLNDs in the pCLN- cohorts, were significantly lower than the pCLN+ cohorts (Table 1). There were 1606 (59.8%) ipsilateral CLNDs and 1081 (40.2%) bilateral CLNDs. The median (IQR) CLNY and CLNR in ipsilateral CLNDs were 7 (7) and 0 (0.2), respectively, and medians in bilateral CLNDs were 11 (9.5) and 0.18 (0.43), respectively. (not shown in tables).

Clinicopathological Characteristics of Patients in the Derivation Group

There were no significant differences in CLN pathology, sex, age, BMI, TSH, size, location, unifocality, laterality, ETE-, HT, CLND extent, CLNY or CLNR between the derivation and validation groups (Table 2). We found significant differences in sex, age, BMI, TSH, size, location, unifocality, laterality, and ETE- between the pCLN- and pCLN+ cohorts in the derivation group based on univariate analysis. However, no significant difference was observed for HT (P=0.687). The number of CLNDs was excluded from our analysis due to collinearity.

TABLE 2

Table 2 Difference between the derivation and validation data groups.

Multivariate analysis indicated that age (P < 0.001), sex (P < 0.001), TSH (P = 0.003), size (P < 0.001), location (P < 0.001), laterality (P = 0.005), unifocality (P < 0.001) and ETE- (P < 0.001) were independent predictors of pCLN-. Overall, BMI had little effect on pCLN-, with no statistical significance was found (P= 0.781) (Table 3).

TABLE 3

Table 3 Univariate and multivariate analysis of pCLN- in the derivation group.

Prediction Models

Based on multivariate logistic regression analysis of the derivation group, we developed a prediction model for pCLN-.

$$\text{Model}:\text{sex}+\text{age}+\text{TSH}+\text{size}+\text{location}+\text{laterality}+\text{unifocality}+\text{ETE}-$$

The nomogram shows the score of each variable on each scale. Therefore, the probability of pCLN- is determined by the total score of all variables. Inputting the necessary clinicopathological data concisely estimated the probability of pCLN- (Figure 2).

FIGURE 2

Figure 2 Nomogram for prediction of pCLN-. A line is drawn straight up to the point axis that corresponds with each patient variable to obtain the points. The sum of these points is located on the total score points axis. A line is drawn downwards to the risk axis to determine the possibility of pCLN-.

According to this model, the ROC curves of the derivation and validation groups were drawn, and the area under the curve (AUC) was 0.784 (95% CI 0.763-0.804) in the derivation group and 0.787 (95% CI 0.756-0.818) in the validation group. The C-indexes of the derivation and validation groups were both 0.783, which demonstrated good discriminative ability. The calibration plot revealed excellent agreement between the predictions and actual observations (Figure 3). These results showed that this nomogram had good efficacy in predicting the probability of pCLN-. DCA was performed in the validation group. The results suggested that when the threshold probability was approximately 10% to 90%, compared to a “treat all” or “treat none” strategy, this nomogram yielded a greater net benefit, which indicated that this nomogram had good clinical value (Figure 4).

FIGURE 3

Figure 3 Receiver operating characteristic curves and calibration plot of the model. In the derivation group (A), the sensitivity and specificity were 67.5% and 76.6%, respectively, with an AUC of 0.784 (95% CI 0.763–0.804). In the validation group (B), the sensitivity and specificity were 67.2% and 75.3%, respectively, with an AUC of 0.787 (95% CI 0.756–0.818). The calibration plot depicts the calibration of the model in terms of the agreement between the predicted possibility of pCLN- and the observed outcomes of pCLN-. The derivation group (C) and the validation group (D) indicated good prediction.

FIGURE 4

Figure 4 Decision curve analysis for the nomogram(in the validation group). The possibility of an individual is denoted as Pi. When Pi reaches a certain threshold (denoted as Pt), it is defined as positive, and intervention (not performing CLND) may be taken. Therefore, there are benefits (pros) of intervention for patients with pCLN- and harm (cons) of intervention for patients with pCLN+. There is also the harm (cons) of missed intervention for pCLN+ patients. Pros minus cons is the net benefit. When Pi < Pt, there is no intervention, and the net benefit is 0 (treat-none). When Pi > Pt, all patients receive the intervention, and the net benefit is shown by a gray slanted curve (treat-all). Our decision curve indicates that when the threshold probability is approximately 10% to 90%, the use of this predictive model would accrue greater benefit than a treat-all or treat-none strategy.

Follow-up Status

We followed 2,497 patients (follow-up rate 92.9%) after the initial surgery until Dec. 2019. For patients without regular evaluations, DFS was based on the time of their last evaluation. The median follow-up period was 30 months (range 3–83 months). Among all 2,497 patients, 2,350 (94.1%) and 1,475 (59.1%) patients received TSH suppression and RAI therapy, respectively. Twenty-one (0.8%) patients experienced recurrence, and 3 (0.1%) died (1 due to asphyxia caused by tracheal metastasis and 2 due to other nontumorous causes). The survival curves plotted according to the Kaplan-Meier method showed that the DFS rate of the 2,497 patients was 99.2%. There was no significant difference in the DFS rate between the pCLN- and pCLN+ cohorts (99.4% vs. 98.9%, P=0.126) (Figure 5). We also found significant differences in Tg levels and cervical ultrasound findings between recurrence and nonrecurrence patients (both P < 0.001), which revealed the importance of periodic monitoring of serum Tg levels and cervical ultrasound.(not shown in tables).

FIGURE 5

Figure 5 The Kaplan–Meier curves for the 2,497 patients. The curves show that the DFS rate was 99.2% in 2,497 patients (A), and there was no significant difference in the DFS rate between the pCLN- and pCLN+ cohorts (B).

Discussion

The present study indicated that for cN0 PTC, sex, age, TSH, size, location, unifocality, laterality, and ETE- were independent predictors of pCLN- but BMI and HT were not predictive.

The pCLN- cohorts had significantly lower CLNY and CLNR. This result was primarily because once the presence of CLNM was confirmed on intraoperative frozen section biopsy, prophylactic CLND was changed to therapeutic CLND, and the extent of CLND was further expanded.

The morbidity of PTC is higher in females than males, but the incidence of CLNM is lower (2, 18, 19). The present study supports these observations, and there was a clear tendency for females to present pCLN- (OR 2.069, P < 0.001). Age is an important basis for TNM staging, and it is one of the factors of many classic prognostic systems: AMES [age, metastasis (distance), extent, size], AGES (age, sex, extent, size), and MACIS [metastasis (distance), age, completeness of surgery, invasiveness, size], which suggests that age is closely related to the biological behavior of PTC. Prognosis and staging may be poorer with increasing age (6, 20, 21). We found that older patients (age > 45 years) generally had a significantly higher incidence of pCLN- (22, 23), especially patients over 65 years old, due to the susceptibility of older patients to lateral or distant metastasis (24, 25). Although no significant difference in pCLN- possibility was found among the three groups in young patients (age ≤ 45 years), patients aged 35 to 45 scored higher on the nomogram scale.

TSH regulates the proliferation and function of thyrocytes, and TSH suppression is an important aspect of the comprehensive treatment for PTC. Many studies showed that TSH was associated with tumorigenesis in PTC (26, 27). An increase in TSH may promote tumor progression, and a higher TSH level is a risk factor for PTC and lymph node metastasis (28, 29). Patients on levothyroxine or antithyroid drugs preoperatively were excluded from the present study, and we further divided the TSH plasma concentrations into 6 ranges. We found that in the 3 ranges of TSH ≤ 3 μIU/mL, the probability of PCLN- increased with decreasing TSH levels, and the opposite was true in the other 3 ranges of TSH > 3 μIU/mL. This result suggests that the level of TSH and the probability of PCLN- are not simple linear correlations.

Size reflects solid tumor staging, and it is related to many aspects, such as treatment and prognosis. PTC with a size ≤10 mm is defined as papillary thyroid microcarcinoma (PTMC), and most of these cases tend to exhibit less CLNM and a better prognosis, which suggests indolent behavior. Yu et al. included 11,845 PTMC patients in the SEER database and found that the LNM incidence of PTMC was 12%. Similarly, Ito et al. revealed an LNM incidence of <5% based on 5-10 years of PTMC data, and all new metastases appeared in the lateral compartment instead of the central compartment (5, 30, 31). Compared to PTC, which has an incidence range of 30-90%, PTMC is less prone to LN involvement, which suggests that PTC of a small size exhibits less nodal invasion (16, 32). We found that even for PTMC, compared to size > 5 but ≤ 10 mm, size ≤ 5 mm significantly increased the possibility of pCLN-, which obtained the highest score on the nomogram. Notably, the tumor size almost linearly negatively correlated with the probability of pCLN- except for the size > 25 but ≤ 30 mm group. The number of patients with tumor sizes between 25 and 30 mm in this study was the lowest(71 in total and 47 in the derivation group). More detailed studies are needed to determine whether this result was due to the decrease in CLNM when the tumor developed to a certain critical size or to the bias caused by an insufficient sample size.

Lymphatic vessels originating at the thyroid form a subcapsular lymphatic network on the surface of the gland. In addition to the superior and inferior thyroid arteries and veins, some lymphatic vessels run along the middle thyroid vein, and some do not accompany the blood vessels. This network is how thyroid-centered lymphatic vessels connect to the surrounding cervical lymph nodes. Different thyroid gland locations have different lymphatic drainage paths. Lymphatic drainage of the upper thyroid is primarily collected by the lymphatic vessels accompanying the superior thyroid artery, with some drainage flowing from the central to the lateral compartment, and some drainage flowing into the venous system through the lymph nodes of the lateral compartment. Lymphatic drainage of the lower thyroid is primarily collected by the lymphatic vessels accompanying the inferior thyroid artery and flowing into the lateral lymph nodes through the CLNs (especially the paratracheal lymph nodes) and ultimately into the venous system. Lymphatic drainage from the isthmus of the thyroid primarily descends to the mediastinal lymph nodes (33, 34). We found that the middle portion was the most common location for PTC (43.7%) and the location with the highest likelihood of pCLN-. This result may be explained by the more complex and crisscrossing lymphatic drainage pathways in the middle portion, and drainage of this region may be carried up or down before proceeding through the upper or lower lymphatic vessels to the corresponding lymph nodes. Therefore, middle tumors progress more slowly than tumors in other portions. There is also a lateral drainage pathway in the middle of the thyroid that passes through the deep surface of the cervical sheath. However, these hypotheses emphasize the need for more accurate anatomical studies to understand the patterns of metastatic spread in PTC. Previous studies reported that an upper location was a risk factor for LLNM, and tumors located in the isthmus and lower portion were more likely to result in CLNM (35, 36). Consistent with these studies, we found that the probability of pCLN- in upper tumors was second only to middle tumors. The probability of pCLN- in isthmic and lower tumors was reduced even further.

Bilaterality or multifocality may develop via intraglandular spread from a single primary tumor or multiple synchronous primary tumors (37). The invasiveness of unilateral and unifocal PTC is weaker than bilateral and multifocal PTC (38, 39). The present study also suggested that CLNs were less prone to invasion in unilateral or unifocal PTC. We further found no significant difference in the probability of pCLN- between the left and right lobes. The left and right lobes of the thyroid gland were not completely symmetrical, but the right lobe was generally larger (40).

ETE- refers to tumors confined to the thyroid capsule. It is easier for the tumor to invade lymph nodes through lymphatics on the surface of the capsule after the thyroid capsule is breached (4, 41). We also found that ETE- was separately and independently associated with pCLN- in cN0 PTC.

The incidence of being overweight increased synchronously with the morbidity of PTC (42), but the relationship between BMI and PTC is not clear. A study of 11,323,006 adults showed that BMI and waist circumference were associated with PTC, and weight loss was associated with lower PTC morbidity (43). Kang et al. and Ray et al. found that fibroblast growth factor 21 (FGF21) positively correlated with BMI, and leptin played a role as a proinflammatory factor in overweight people. FGF21 and leptin were also involved in the regulation of PTC development, which indicated that BMI was also related to invasion and metastasis in PTC (44, 45).The BMI level of the pCLN- cohort in the present study was slightly lower, but there was no significant correlation between BMI and pCLN-, which is consistent with the results obtained by Grani et al. and Paes et al. (46, 47). Because being overweight may involve aromatase activation, chronic inflammation and oxidative stress, more molecular biological mechanisms that affect the relationship between BMI, PTC and invasiveness must be explored.

HT is a type of autoimmune thyroiditis that is characterized by an abundance of lymphocytes, and it also known as chronic lymphocyte thyroiditis. The relationship between HT and CLN status in PTC is inconclusive (48, 49), and the differences in diagnostic criteria, severity, antibody levels, and the heterogeneity of HT may underlie these controversies. HT was confirmed via postoperative pathology in this study, and we found that HT barely affected CLN status. However, this study did not determine the level of autoantibodies (primarily TPO-Ab and Tg-Ab), and quantitative analysis of the degree of HT before surgery needs further study.

The follow-up results showed that very few patients experienced lymph node recurrence, which confirmed that the lymph node pathologic status, especially in the pCLN- cohort was reliable. Specifically, the derivation basis of the nomogram was reliable. Because the present study was a single-center retrospective study, selection deviation is inevitable. Notably, the lack of external validation is a limitation, and relevant prospective multicenter clinical studies should be performed in the future to evaluate the accuracy of the proposed model. Analysis of potential confounders of future metastasis was not possible due to the relatively short follow-up, which suggests that prolonged follow-up and more data are required in our further work.

There have been many prediction models of CLNM in cN0 PTC (50, 51). However, available prediction models for other patients who present pCLN- and constitute the majority of the study population (the pCLN- rate was 51.5% in this study) are lacking. In conclusion, we established the first quantified prediction model for pCLN- in cN0 PTC. Clinicians can use this nomogram to evaluate patients with PTC based on clinicopathological characteristics. For cN0 PTC, which has a high probability of pCLN-, unnecessary prophylactic CLND can may be avoided.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee of the First Affiliated Hospital of Chongqing Medical University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

XJS and XLS designed this study. XJS, LT, DH, YW, PY, ZY, CD, and DW collected the data. XJS analyzed the data. All authors contributed to the article and approved the submitted version.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors appreciate all of the staff involved in the preparation of the study.

Abbreviations

cN0, clinical lymph node negativity; pN0, pathologic lymph node negativity; pCLN-, pathologic central lymph node negativity; PTC, papillary thyroid carcinoma; PTMC, papillary thyroid microcarcinoma; CLN, central lymph node; CLNM, central lymph node metastasis; CLND, central lymph node dissection; SET, scarless endoscopic thyroidectomy; BMI, body mass index; TSH, thyroid stimulating hormone; Tg, thyroglobulin; HT, Hashimoto’s thyroiditis; ETE-, extrathyroidal extension negativity (tumor confined to the thyroid); IQR, interquartile range; CLNY,central lymph node yield; CLNR,central lymph node ratio; AUC, area under curve; DFS, disease-free survival; OR, odd ratio; CI, confidence index.

References

1. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2020. CA Cancer J Clin (2020) 70(1):7–30. doi: 10.3322/caac.21590

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Cabanillas ME, McFadden DG, Durante C. Thyroid Cancer. Lancet (2016) 388(10061):2783–95. doi: 10.1016/S0140-6736(16)30172-6

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Burgess JR, Tucker P. Incidence Trends for Papillary Thyroid Carcinoma and Their Correlation With Thyroid Surgery and Thyroid Fine-Needle Aspirate Cytology. Thyroid (2006) 16(1):47–53. doi: 10.1089/thy.2006.16.47

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Goyal N, Pakdaman M, Kamani D, Caragacianu D, Goldenberg D, Randolph GW. Mapping the Distribution of Nodal Metastases in Papillary Thyroid Carcinoma: Where Exactly are the Nodes? Laryngoscope (2017) 127(8):1959–64. doi: 10.1002/lary.26495

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Haugen BR, Alexander EK, Bible KC, Doherty GM, Mandel SJ, Nikiforov YE, et al. 2015 American Thyroid Association Management Guidelines for Adult Patients With Thyroid Nodules and Differentiated Thyroid Cancer: The American Thyroid Association Guidelines Task Force on Thyroid Nodules and Differentiated Thyroid Cancer. Thyroid (2016) 26(1):1–133. doi: 10.1089/thy.2015.0020

PubMed Abstract | CrossRef Full Text | Google Scholar

6. American Thyroid Association Surgery Working G, American Association of Endocrine S, American Academy of O-H, Neck S, American H, Neck S, et al. Consensus Statement on the Terminology and Classification of Central Neck Dissection for Thyroid Cancer. Thyroid (2009) 19(11):1153–8. doi: 10.1089/thy.2009.0159

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Amin MB, Gress DM. AJCC Cancer Staging Manual. Eight Edition. New York: Springer (2017). doi: 10.1007/s00106-017-0391-3

CrossRef Full Text | Google Scholar

8. Eun NL, Son EJ, Kim JA, Gweon HM, Kang JH, Youk JH. Comparison of the Diagnostic Performances of Ultrasonography, CT and Fine Needle Aspiration Cytology for the Prediction of Lymph Node Metastasis in Patients With Lymph Node Dissection of Papillary Thyroid Carcinoma: A Retrospective Cohort Study. Int J Surg (2018) 51:145–50. doi: 10.1016/j.ijsu.2017.12.036

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Rosario PW, de Faria S, Bicalho L, Alves MF, Borges MA, Purisch S, et al. Ultrasonographic Differentiation Between Metastatic and Benign Lymph Nodes in Patients With Papillary Thyroid Carcinoma. J Ultrasound Med (2005) 24(10):1385–9. doi: 10.7863/jum.2005.24.10.1385

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Hoang JK, Lee WK, Lee M, Johnson D, Farrell S. US Features of Thyroid Malignancy: Pearls and Pitfalls. Radiographics (2007) 27(3):847–60. doi: 10.1148/rg.273065038 discussion 861-845.

PubMed Abstract | CrossRef Full Text | Google Scholar

11. National Comprehensive Cancer Network. (NCCN) Clinical Practice Guidelines in Oncology, in: Thyroid Carcinoma Version1 (2021). Available at: https://www.nccn.org (Accessed 25 Aug 2021).

Google Scholar

12. Roh JL, Park JY, Park CI. Total Thyroidectomy Plus Neck Dissection in Differentiated Papillary Thyroid Carcinoma Patients: Pattern of Nodal Metastasis, Morbidity, Recurrence, and Postoperative Levels of Serum Parathyroid Hormone. Ann Surg (2007) 245(4):604–10. doi: 10.1097/01.sla.0000250451.59685.67

PubMed Abstract | CrossRef Full Text | Google Scholar

13. So YK, Son YI, Hong SD, Seo MY, Beak CH, Jeong HS, et al. Subclinical Lymph Node Metastasis in Papillary Thyroid Microcarcinoma: A Study of 551 Resections. Surgery (2010) 148(3):526–31. doi: 10.1016/j.surg.2010.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Popadich A, Levin O, Lee JC, Smooke-Praw S, Ro K, Fazel M, et al. A Multicenter Cohort Study of Total Thyroidectomy and Routine Central Lymph Node Dissection for Cn0 Papillary Thyroid Cancer. Surgery (2011) 150(6):1048–57. doi: 10.1016/j.surg.2011.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Kouvaraki MA, Shapiro SE, Fornage BD, Edeiken-Monro BS, Sherman SI, Vassilopoulou-Sellin R, et al. Role of Preoperative Ultrasonography in the Surgical Management of Patients With Thyroid Cancer. Surgery (2003) 134(6):946–54. doi: 10.1016/s0039-6060(03)00424-0 discussion 954-945.

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Hall FT, Beasley NJ, Eski SJ, Witterick IJ, Walfish PG, Freeman JL. Predictive Value of Serum Thyroglobulin After Surgery for Thyroid Carcinoma. Laryngoscope (2003) 113(1):77–81. doi: 10.1097/00005537-200301000-00014

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Nascimento C, Borget I, Troalen F, Al Ghuzlan A, Deandreis D, Hartl D, et al. Ultrasensitive Serum Thyroglobulin Measurement is Useful for the Follow-Up of Patients Treated With Total Thyroidectomy Without Radioactive Iodine Ablation. Eur J Endocrinol (2013) 169(5):689–93. doi: 10.1530/EJE-13-0386

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Smith VA, Sessions RB, Lentsch EJ. Cervical Lymph Node Metastasis and Papillary Thyroid Carcinoma: Does the Compartment Involved Affect Survival? Experience From the SEER Database. J Surg Oncol (2012) 106(4):357–62. doi: 10.1002/jso.23090

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Zheng X, Peng C, Gao M, Zhi J, Hou X, Zhao J, et al. Risk Factors for Cervical Lymph Node Metastasis in Papillary Thyroid Microcarcinoma: A Study of 1,587 Patients. Cancer Biol Med (2019) 16(1):121–30. doi: 10.20892/j.issn.2095-3941.2018.0125

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Ito Y, Kudo T, Kobayashi K, Miya A, Ichihara K, Miyauchi A. Prognostic Factors for Recurrence of Papillary Thyroid Carcinoma in the Lymph Nodes, Lung, and Bone: Analysis of 5,768 Patients With Average 10-Year Follow-Up. World J Surg (2012) 36(6):1274–8. doi: 10.1007/s00268-012-1423-5

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Kauffmann RM, Hamner JB, Ituarte PHG, Yim JH. Age Greater Than 60 Years Portends a Worse Prognosis in Patients With Papillary Thyroid Cancer: Should There be Three Age Categories for Staging? BMC Cancer (2018) 18(1):316. doi: 10.1186/s12885-018-4181-4

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Vriens MR, Moses W, Weng J, Peng M, Griffin A, Bleyer A, et al. Clinical and Molecular Features of Papillary Thyroid Cancer in Adolescents and Young Adults. Cancer (2011) 117(2):259–67. doi: 10.1002/cncr.25369

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Jianming L, Jibin L, Linxue Q. Suspicious Ultrasound Characteristics Correlate With Multiple Factors That Predict Central Lymph Node Metastasis of Papillary Thyroid Carcinoma: Significant Role of HBME-1. Eur J Radiol (2020) 123:108801. doi: 10.1016/j.ejrad.2019.108801

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Shaha AR, Shah JP, Loree TR. Differentiated Thyroid Cancer Presenting Initially With Distant Metastasis. Am J Surg (1997) 174(5):474–6. doi: 10.1016/s0002-9610(97)00158-x

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Back K, Kim JS, Kim JH, Choe JH. Superior Located Papillary Thyroid Microcarcinoma is a Risk Factor for Lateral Lymph Node Metastasis. Ann Surg Oncol (2019) 26(12):3992–4001. doi: 10.1245/s10434-019-07587-2

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Fiore E, Rago T, Latrofa F, Provenzale MA, Piaggi P, Delitala A, et al. Hashimoto’s Thyroiditis is Associated With Papillary Thyroid Carcinoma: Role of TSH and of Treatment With L-Thyroxine. Endocr Relat Cancer (2011) 18(4):429–37. doi: 10.1530/ERC-11-0028

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Kim SS, Lee BJ, Lee JC, Song SH, Kim BH, Son SM, et al. Preoperative Serum Thyroid Stimulating Hormone Levels in Well-Differentiated Thyroid Carcinoma is a Predictive Factor for Lateral Lymph Node Metastasis as Well as Extrathyroidal Extension in Korean Patients: A Single-Center Experience. Endocrine (2011) 39(3):259–65. doi: 10.1007/s12020-010-9430-5

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Kohler H, Latteyer S, Hones GS, Theurer S, Liao XH, Christoph S, et al. Increased Anaplastic Lymphoma Kinase Activity Induces a Poorly Differentiated Thyroid Carcinoma in Mice. Thyroid (2019) 29(10):1438–46. doi: 10.1089/thy.2018.0526

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Fiore E, Rago T, Provenzale MA, Scutari M, Ugolini C, Basolo F, et al. Lower Levels of TSH are Associated With a Lower Risk of Papillary Thyroid Cancer in Patients With Thyroid Nodular Disease: Thyroid Autonomy may Play a Protective Role. Endocr Relat Cancer (2009) 16(4):1251–60. doi: 10.1677/ERC-09-0036

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Ito Y, Miyauchi A, Inoue H, Fukushima M, Kihara M, Higashiyama T, et al. An Observational Trial for Papillary Thyroid Microcarcinoma in Japanese Patients. World J Surg (2010) 34(1):28–35. doi: 10.1007/s00268-009-0303-0

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Yu XM, Wan Y, Sippel RS, Chen H. Should All Papillary Thyroid Microcarcinomas be Aggressively Treated? An Analysis of 18,445 Cases. Ann Surg (2011) 254(4):653–60. doi: 10.1097/SLA.0b013e318230036d

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Maciel J, Donato S, Simões H, Leite V. Clinical Outcomes of a Conservative Approach in Cervical Lymph Node Metastases of Thyroid Cancer. Clin Endocrinol (Oxf) (2020) 94(3):460–5. doi: 10.1111/cen.14306

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Qubain SW, Nakano S, Baba M, Takao S, Aikou T. Distribution of Lymph Node Micrometastasis in Pn0 Well-Differentiated Thyroid Carcinoma. Surgery (2002) 131(3):249–56. doi: 10.1111/cen.14306

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Likhterov I, Reis LL, Urken ML. Central Compartment Management in Patients With Papillary Thyroid Cancer Presenting With Metastatic Disease to the Lateral Neck: Anatomic Pathways of Lymphatic Spread. Head Neck (2017) 39(5):853–9. doi: 10.1002/hed.24568

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Hu D, Zhou J, He W, Peng J, Cao Y, Ren H, et al. Risk Factors of Lateral Lymph Node Metastasis in Cn0 Papillary Thyroid Carcinoma. World J Surg Oncol (2018) 16(1):30. doi: 10.1186/s12957-018-1336-3

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Zhang TT, Qi XZ, Chen JP, Shi RL, Wen SS, Wang YL, et al. The Association Between Tumor’s Location and Cervical Lymph Nodes Metastasis in Papillary Thyroid Cancer. Gland Surg (2019) 8(5):557–68. doi: 10.21037/gs.2019.10.02

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Bansal M, Gandhi M, Ferris RL, Nikiforova MN, Yip L, Carty SE, et al. Molecular and Histopathologic Characteristics of Multifocal Papillary Thyroid Carcinoma. Am J Surg Pathol (2013) 37(10):1586–91. doi: 10.1097/PAS.0b013e318292b780

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Wang W, Su X, He K, Wang Y, Wang H, Wang H, et al. Comparison of the Clinicopathologic Features and Prognosis of Bilateral Versus Unilateral Multifocal Papillary Thyroid Cancer: An Updated Study With More Than 2000 Consecutive Patients. Cancer (2016) 122(2):198–206. doi: 10.1002/cncr.29689

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Dong S, Xie XJ, Xia Q, Wu YJ. Indicators of Multifocality in Papillary Thyroid Carcinoma Concurrent With Hashimoto’s Thyroiditis. Am J Cancer Res (2019) 9(8):1786–95.

PubMed Abstract | Google Scholar

40. Dou Y, Chen Y, Hu D, Su X. The Recovery of Thyroid Function in Low-Risk Papillary Thyroid Cancer After Lobectomy: A 3-Year Follow-Up Study. Front Endocrinol (Lausanne) (2020) 11:619841. doi: 10.3389/fendo.2020.619841

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Kim SK, Park I, Woo JW, Lee JH, Choe JH, Kim JH, et al. Predictive Factors for Lymph Node Metastasis in Papillary Thyroid Microcarcinoma. Ann Surg Oncol (2016) 23(9):2866–73. doi: 10.1245/s10434-016-5225-0

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Kitahara CM, Pfeiffer RM, Sosa JA, Shiels MS. Overweight and Obesity on US Papillary Thyroid Cancer Incidence Trends (1995-2015). J Natl Cancer Inst (2020) 112(8):810–7. doi: 10.1093/jnci/djz202

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Kwon H, Han KD, Park CY. Weight Change is Significantly Associated With Risk of Thyroid Cancer: A Nationwide Population-Based Cohort Study. Sci Rep (2019) 9(1):1546. doi: 10.1038/s41598-018-38203-0

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Kang YE, Kim JT, Lim MA, Oh C, Liu L, Jung SN, et al. Association Between Circulating Fibroblast Growth Factor 21 and Aggressiveness in Thyroid Cancer. Cancers (Basel) (2019) 11(8). doi: 10.3390/cancers11081154

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Ray A, Cleary MP. The Potential Role of Leptin in Tumor Invasion and Metastasis. Cytokine Growth Factor Rev (2017) 38:80–97. doi: 10.1016/j.cytogfr.2017.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Grani G, Lamartina L, Montesano T, Ronga G, Maggisano V, Falcone R, et al. Lack of Association Between Obesity and Aggressiveness of Differentiated Thyroid Cancer. J Endocrinol Invest. (2019) 42(1):85–90. doi: 10.1007/s40618-018-0889-x

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Paes JE, Hua K, Nagy R, Kloos RT, Jarjoura D, Ringel MD. The Relationship Between Body Mass Index and Thyroid Cancer Pathology Features and Outcomes: A Clinicopathological Cohort Study. J Clin Endocrinol Metab (2010) 95(9):4244–50. doi: 10.1210/jc.2010-0440

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Jara SM, Carson KA, Pai SI, Agrawal N, Richmon JD, Prescott JD, et al. The Relationship Between Chronic Lymphocytic Thyroiditis and Central Neck Lymph Node Metastasis in North American Patients With Papillary Thyroid Carcinoma. Surgery (2013) 154(6):1272–80. doi: 10.1016/j.surg.2013.07.021 discussion 1280-1282.

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Wen X, Wang B, Jin Q, Zhang W, Qiu M. Thyroid Antibody Status is Associated With Central Lymph Node Metastases in Papillary Thyroid Carcinoma Patients With Hashimoto’s Thyroiditis. Ann Surg Oncol (2019) 26(6):1751–8. doi: 10.1245/s10434-019-07256-4

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Thompson AM, Turner RM, Hayen A, Aniss A, Jalaty S, Learoyd DL, et al. A Preoperative Nomogram for the Prediction of Ipsilateral Central Compartment Lymph Node Metastases in Papillary Thyroid Cancer. Thyroid (2014) 24(4):675–82. doi: 10.1089/thy.2013.0224

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Wang Y, Guan Q, Xiang J. Nomogram for Predicting Central Lymph Node Metastasis in Papillary Thyroid Microcarcinoma: A Retrospective Cohort Study of 8668 Patients. Int J Surg (2018) 55:98–102. doi: 10.1016/j.ijsu.2018.05.023

PubMed Abstract | CrossRef Full Text | Google Scholar