A Clinical Predictive Model of Central Lymph Node Metastases in Papillary Thyroid Carcinoma

Background: Thyroid carcinoma is one of the most common endocrine tumors, and papillary thyroid carcinoma (PTC) is the most common pathological type. Current studies have reported that PTC has a strong propensity for central lymph node metastases (CLNMs). Whether to prophylactically dissect the central lymph nodes in PTC remains controversial. This study aimed to explore the risk factors and develop a predictive model of CLNM in PTC.

Methods: A total of 2,554 patients were enrolled in this study. The basic information, laboratory examination, characteristics of cervical ultrasound, genetic test, and pathological diagnosis were collected. The collected data were analyzed by univariate logistic analysis and multivariate logistic analysis. The risk factors were evaluated, and the predictive model was constructed of CLNM.

Results: The multivariate logistic analysis showed that Age (p < 0.001), Gender (p < 0.001), Multifocality (p < 0.001), BRAF (p = 0.027), and Tumor size (p < 0.001) were associated with CLNM. The receiver operating characteristic curve (ROC curve) showed high efficiency with an area under the ROC (AUC) of 0.781 in the training group. The calibration curve and the calibration of the model were evaluated. The decision curve analysis (DCA) for the nomogram showed that the nomogram can provide benefits in this study.

Conclusion: The predictive model of CLNM constructed and visualized based on the evaluated risk factors was confirmed to be a practical and convenient tool for clinicians to predict the CLNM in PTC.

Introduction

Thyroid carcinoma is one of the most popular endocrine tumors, and papillary thyroid carcinoma (PTC) is the most common pathological type (1, 2). The increased incidence of PTC is attributed to both the truly increased prevalence of thyroid diseases and the advances in imaging technology. Current studies have reported that PTC has a strong propensity for CLNM (3), and it is difficult to effectively detect central lymph node metastases (CLNMs) preoperatively (4–6). The most commonly involved central lymph nodes in thyroid carcinoma are the prelaryngeal (Delphian), pretracheal, and the right and left paratracheal nodes; the paratracheal nodes may be anterior as well as posterior to the recurrent laryngeal nerves (3). Whether to prophylactically dissect the central lymph nodes in PTC remains controversial (7). Dissection of the lymph nodes in the central region of the neck is considered necessary. Central lymph node dissection (CLND) is beneficial in eliminating macroscopic or microscopic metastatic sites. When CLNM is found in patients after surgery, a second operation is often necessary. Extra surgery not only is difficult but also increases the risk of complications (8). Prophylactic lymph node dissection of the central cervical facilitates accurate clinical staging (9). However, it has been argued that prophylactic dissection of the central lymph nodes is not necessary. Routine prophylactic CLND is considered uneconomical (10). Routine prophylactic central compartment dissection is particularly associated with temporary and permanent hypoparathyroidism and recurrent laryngeal nerve injury (4, 11, 12). Although numerous retrospective studies determined the benefits of prophylactic CLND by various investigators, the published results have been inconclusive. Due to the indolent nature of PTC, a very large sample size with extended follow-up for a randomized control trial examination is urgently needed. Therefore, preoperative access is of great clinical significance to accurately assess patients for CLNM. In this study, we analyzed the risk factors of CLNM and constructed a predictive model of CLNM to more accurately assess the risk of CLNM preoperatively and provide a practical and convenient tool for clinicians to predict the CLNM in PTC for clinical decision-making.

Methods

Study Design

We screened all the patients who underwent thyroidectomy for PTC in Thyroid Surgery who were admitted to the First Affiliated Hospital of Zhengzhou University from January 2018 to October 2019. The inclusion criteria were as follows: 1) patients underwent thyroid surgery for the first time; 2) clinically and pathologically diagnosed as PTC. The exclusion criteria included the following: 1) other malignancies combined; 2) preoperative hyperthyroidism or hypothyroidism; 3) other diseases that cause swollen lymph nodes in the neck; 4) two or more thyroidectomies. A total of 2,554 patients including 982 patients with CLNMs (CLNM(+)), and 1,572 patients without CLNMs (CLNM(−)) were enrolled in this study. The enrollment flowchart of the participants is shown in Figure 1. The patients with CLNMs (CLNM(+)) and without CLNMs (CLNM(−)) were divided into the training group (n = 1,787) and the validation group (n = 767).

FIGURE 1

Figure 1 Enrollment flowchart of participants used for model development and validation. PTC, papillary thyroid carcinoma; CLNM, central lymph node metastasis.

Data Collection

The basic information, laboratory examination, cervical ultrasound, genetic test, and pathological diagnosis were collected. The basic information included age and gender. The laboratory indices included free triiodothyronine (FT₃), free tetraiodothyronine (FT₄), thyroid-stimulating hormone (TSH), thyroid peroxidase antibodies (TPOAb), thyroglobulin antibodies (TgAb), and thyroglobulin (Tg). The characteristics of cervical ultrasound include multifocality and tumor size. Multifocality was defined as more than one lesion observed in cervical ultrasound and pathologically confirmed as PTC. Tumor size was the maximum diameter of the suspected nodule under ultrasound that was pathologically confirmed as PTC (8). All patients had a review of the cervical ultrasound with the same sonographer preoperatively. Genetic test results included BRAF and TERT. Pathological diagnosis included reports of the paraffin section of the primary lesion and central lymph node.

Statistical Analysis

Multivariate multiple imputations with chained equations were used to deal with a few missing data of several variables to decrease the bias (13). All the statistical analysis processes involved were completed by R software, version 4.1.1. p-Value <0.05 was considered statistically significant. The classification data were expressed as percentages, and means ± SD or medians (quartile 1, quartile 3) were described as continuous variables that satisfy or do not satisfy the normal distribution, respectively. The odds ratio (OR) values were calculated by univariate logistic regression of the variables. After the collinearity among variables was calculated and the colinear factors were eliminated, the potential variables with a p-value <0.05 were selected to perform the multivariate logistic regression. A clinical predictive model of CLNM in PTC was built based on the variables with statistical senses.

Results

Baseline Characteristics of the Variables

The baseline characteristics of CLNM(+) and CLNM(−) are shown in Table 1. All the enrolled patients were randomly divided into the training group (n = 1,787) and the validation group (n = 767). There was no significant difference in the levels of these variables between the two groups (Table 2). Six potential predictors from 12 candidates were considered to have statistical significance (p < 0.05).

TABLE 1

Table 1 Baseline characteristics of enrolled patients.

TABLE 2

Table 2 Baseline characteristics showed there was no statistical difference between the training group and validation group.

Risk Factors and Predictive Model of Central Lymph Node Metastasis

Multivariate logistic analysis showed that Age (OR 0.947–0.966), Gender (OR 1.541–2.589), Multifocality (OR 2.348–3.815), BRAF (OR 1.044–1.926), and Tumor size (OR 4.127–6.416) were associated with CLNM (Table 3). The receiver operating characteristic curve (ROC curve) was drawn to evaluate the diagnostic effectiveness of the model (Figure 2). The ROC showed a high efficiency with an area under the ROC (AUC) of 0.781 (specificity 0.778, sensitivity 0.662, Figure 2A) in the training group. The effectiveness was verified in the validation group with an AUC of 0.736 (specificity 0.631, sensitivity 0.774), and the result is shown in Figure 2B. The model showed a great ability to distinguish the presence or absence of CLNM in the training group with a high value of AUC.

TABLE 3

Table 3 Potential risk factors identified by univariate and multivariate logistic regression analyses.

FIGURE 2

Figure 2 The differential capability of the nomogram. (A) ROC curve based on the potential risk factors identified by multivariate logistic regression analysis showed great ability to distinguish the presence or absence of CLNM in the training group with a high value of AUC. (B) ROC curve based on the potential risk factors identified by multivariate logistic regression analysis showed great ability to distinguish the presence or absence of CLNM in the validation group. ROC, receiver operating characteristic; CLNM, central lymph node metastasis; AUC, area under the curve.

Then the model was evaluated by the calibration curve. The predicted values had good consistency in the training group (mean absolute error = 0.004) and the validation group (mean absolute error = 0.008) with the observed variables. The calibration curve showed that the model had a strong calibration ability (Figure 3).

FIGURE 3

Figure 3 Calibration curve of the predictive nomogram in the (A) training group or (B) validation group.

To visualize the model, we plotted the nomogram of our predictive model based on the five variables: Age, Gender, Focal, BRAF, and Tumor size. Every variable was scored by drawing a straight line upward the “Points” line. The total points were the sum of the points obtained by the five variables. A straight line down to the axis named “CLNM risk” represents the risk of CLNM (Figure 4).

FIGURE 4

Figure 4 Nomogram based on the five variables including Age, Gender, Focal, BRAF, and Tumor size.

In the training group, we constructed a decision curve analysis (DCA) to identify the net benefit of the nomogram (Figure 5). The curve showed that when the threshold probability of patients is between 0.11 and 0.90, the nomogram can provide benefits.

FIGURE 5

Figure 5 Decision curve for the nomogram predicting CLNM in (A) training group or (B) validation group. CLNM, central lymph node metastasis.

Discussion

PTC metastasis occurs most often in the central lymph nodes, with few distant metastases and low mortality, and the rate of cervical lymph node metastasis detection by ultrasound is unsatisfactory, especially for CLNM (5, 14). Information gathered by prospective and randomized clinical studies is the key to determining whether prophylactic dissection of the central lymph nodes in PTC is necessary, but it is probably unavailable currently. Therefore, identification of patients with PTC preoperatively at greater risk of CLNM would be valuable. Prediction models and risk factors analysis based on clinical data have been developed increasingly in a wide variety of diseases in recent years.

To solve these problems, in this study, we collected and analyzed the risk factors and constructed a predictive model of CLNM. The model also showed a high calibration capacity. To the best of our knowledge, it is the first study aiming to construct and visualize a predictive model of CLNM. It is a predictive model that can help clinicians make appropriate treatments and help clinicians assess whether patients need CLND in PTC.

Some studies suggested that CLNM is significantly correlated with age (15, 16), which was similar to our study. Patients with younger age were more likely to be considered at a higher risk of CLNM.

Recent studies have shown that estrogen is a powerful stimulant for benign and malignant thyroid nodules. This explains why thyroid cancer is highly prevalent in women (17). However, in this study, there was a higher rate of CLNM in men. Some scholars confirmed that there were different subtypes of estrogen receptors (ERs) that were considered a protective factor in PTC (17). This may be a potential cause of a higher rate of CLNM in male patients. The detailed knowledge of this regulation in thyroid cancer is still being debated.

Studies have shown that the risk of CLNM increases with the number of foci (11). The multifocal disease was defined as the presence of 2 or more foci of PTC, and each focus was recorded separately. In this study, CLNM rates were high in multifocal PTC with an OR of 2.989 (95% CI 2.348–3.815). With the application of high-resolution ultrasound, the preoperative diagnostic technology of PTC has made a huge breakthrough (18). Therefore, multifocality should be taken seriously during the preoperative ultrasound. When the cervical ultrasound showed that the suspicious lesions have significant multifocality, lobectomy and prophylactic CLNM should be considered.

BRAF genetic test was valuable for the diagnosis, prognosis, and therapy of PTC (19). In addition, some scholars revealed that BRAF mutations were associated with markers of clinical aggressiveness such as larger tumors, lymph node metastases, and poor clinical outcomes (20). In general, BRAF genetic test had strong practicability. Genetic testing in preoperative fine-needle aspiration biopsy (FNAB) was helpful to confirm the diagnosis, patients with a positive result for BRAF genetic test should be further evaluated, and there were stronger recommendations for prophylactic CLND.

Tumor size is an important risk factor for CLNM in PTC in the present study. This finding was similar to other studies that the risk of CLNM in PTC increases with tumor size (21, 22). In clinical practice, suspicious nodules larger than 10 mm should be managed with caution. Clinicians need to evaluate the patient’s cervical lymph nodes to decide whether to perform prophylactic dissection of the lymph nodes. Our research, especially the nomogram, provides a good reference for clinicians.

In the univariate logistic regression model, the titer of FT₃ was statistically significant between the two groups. We searched the relevant literature and found that there is no clear relationship between thyroid hormone and CLNM of PTC. However, one scholar’s research results caught our attention; this study showed a high correlation between PTC microcalcification and thyroid hormones (23). Microcalcification of thyroid nodules indicates that it was more likely to be malignant (24). Therefore, the titer of thyroid hormone may be related to the pathology of thyroid nodules, but it cannot be considered a risk factor for CLNM.

Some scholars have reached similar conclusions (25, 26); younger age, male sex, multifocality, and larger tumor size are the risk factors for cervical lymph node metastases in PTC. The predictive model was also constructed and visualized with a nomogram. However, our study aimed at risk factor analysis and predictive model construction for CLNM with a high AUC. In addition, not only the clinical baseline characteristics but also the genetic test result was included in this study. This was more instructive for patients with preoperative FNAB and BRAF genetic test.

The decision curve and nomogram in this study show the great utility of our model, and the 5 items in the nomogram are routine clinical variables that can easily be obtained by clinicians, indicating that it may be beneficial for clinicians to assess the need for prophylactic CLND. Our study has a large sample size of 2,554 patients and has excellent diagnostic effectiveness of CLNM with an AUC of 0.781 in the training group. The operation of the model is simple and fast, which can provide a reference for timely prophylactic CLND.

However, there are also several limitations in our study. First, all the enrolled patients came from the same hospital without external validation. Moreover, we still need to expand the sample size to reduce the heterogeneity. In addition, our predictive model is only suitable for PTC; there is still a lack of predictive ability of our model for other types of thyroid carcinoma.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

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

Author Contributions

DY conceived of the idea and provided guidance. ZW wrote the manuscript and completed the figures. QC and HZ contributed to organizing the database. SL, YL, and HS carefully reviewed the manuscript. GD made critical revisions to the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This study was also funded by the Key Medical Science and Technology Project of Henan Province (SBGJ202101014), which is also from the corresponding author DY, should be placed before funding Major Scientific Research Projects of Traditional Chinese Medicine in Henan Province (No. 20-21ZYZD14).

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

We acknowledge Qungang Chang for the data collection.

References

1. Xu S, Huang H, Qian J, Liu Y, Huang Y, Wang X, et al. Prevalence of Hashimoto Thyroiditis in Adults With Papillary Thyroid Cancer and Its Association With Cancer Recurrence and Outcomes. JAMA Netw Open (2021) 4(7):e2118526. doi: 10.1001/jamanetworkopen.2021.18526

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Yasuhiro Ito MK, Takamura Y, Kobayashi K, Miya A, Miyauchi A. Prognosis and Prognostic Factors of Patients With Papillary Thyroid Carcinoma Requiring Resection of Recurrent Laryngeal Nerve Due to Carcinoma Extension. Endocrine J (2012) 59(3):247–52. doi: 10.1507/endocrj.ej11-0355

CrossRef Full Text | Google Scholar

3. Carty SE, Doherty GM, Duh Q-Y, Kloos RT, Mandel SJ, Randolph GW, et al. Consensus Statement on the Terminology and Classification of Central Neck Dissection for Thyroid Cancer. Thyroid (2009) 11:2009. doi: 10.1089/thy.2009.0159

CrossRef Full Text | Google Scholar

4. Sturgeon C, Yang A, Elaraj D. Surgical Management of Lymph Node Compartments in Papillary Thyroid Cancer. Surg Oncol Clin N Am (2016) 25(1):17–40. doi: 10.1016/j.soc.2015.08.013

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Kim KE, Kim EK, Yoon JH, Han KH, Moon HJ, Kwak JY. Preoperative Prediction of Central Lymph Node Metastasis in Thyroid Papillary Microcarcinoma Using Clinicopathologic and Sonographic Features. World J Surg (2013) 37(2):385–91. doi: 10.1007/s00268-012-1826-3

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Mazzaferri ELMD, Doherty GM, Steward DL. The Pros and Cons of Prophylactic Central Compartment Lymph Node Dissection for Papillary Thyroid Carcinoma. Thyroid (2009) 19(7) 2009. doi: 10.1089/thy.2009.1578

CrossRef Full Text | Google Scholar

7. Fritze D, Doherty GM. Surgical Management of Cervical Lymph Nodes in Differentiated Thyroid Cancer. Otolaryngol Clin North Am (2010) 43(2):285–300. doi: 10.1016/j.otc.2010.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Guang Y, He W, Zhang W, Zhang H, Zhang Y, Wan F. Clinical Study of Ultrasonographic Risk Factors for Central Lymph Node Metastasis of Papillary Thyroid Carcinoma. Front Endocrinol (Lausanne) (2021) 12:791970. doi: 10.3389/fendo.2021.791970

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Hughes DT, White ML, Miller BS, Gauger PG, Burney RE, Doherty GM. Influence of Prophylactic Central Lymph Node Dissection on Postoperative Thyroglobulin Levels and Radioiodine Treatment in Papillary Thyroid Cancer. Surgery (2010) 148(6):1100–6; discussion 006-7. doi: 10.1016/j.surg.2010.09.019

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Garcia A, Palmer BJ, Parks NA, Liu TH. Routine Prophylactic Central Neck Dissection for Low-Risk Papillary Thyroid Cancer Is Not Cost-Effective. Clin Endocrinol (Oxf) (2014) 81(5):754–61. doi: 10.1111/cen.12506

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Al Afif A, Williams BA, Rigby MH, Bullock MJ, Taylor SM, Trites J, et al. Multifocal Papillary Thyroid Cancer Increases the Risk of Central Lymph Node Metastasis. Thyroid (2015) 25(9):1008–12. doi: 10.1089/thy.2015.0130

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Sancho JJ, Lennard TW, Paunovic I, Triponez F, Sitges-Serra A. Prophylactic Central Neck Disection in Papillary Thyroid Cancer: A Consensus Report of the European Society of Endocrine Surgeons (ESES). Langenbecks Arch Surg (2014) 399(2):155–63. doi: 10.1007/s00423-013-1152-8

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Buuren SG-O. Mice: Multivariate Imputation by Chained Equations in R. J Stat Software (2011) 45(3):1–67. doi: 10.18637/jss.v045.i03

CrossRef Full Text | Google Scholar

14. Zhao W, He L, Zhu J, Su A. A Nomogram Model Based on the Preoperative Clinical Characteristics of Papillary Thyroid Carcinoma With Hashimoto’s Thyroiditis to Predict Central Lymph Node Metastasis. Clin Endocrinol (Oxf) (2021) 94(2):310–21. doi: 10.1111/cen.14302

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Ito Y, Miyauchi A, Kihara M, Higashiyama T, Kobayashi K, Miya A. Patient Age is Significantly Related to the Progression of Papillary Microcarcinoma of the Thyroid Under Observation. Thyroid (2014) 24(1):27–34. doi: 10.1089/thy.2013.0367

PubMed Abstract | CrossRef Full Text | Google Scholar

16. 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

17. Derwahl M, Nicula D. Estrogen and Its Role in Thyroid Cancer. Endocr Relat Cancer (2014),21(15) T273–83. doi: 10.1530/ERC-14-0053

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Khokhar MT, Day KM, Sangal RB, Ahmedli NN, Pisharodi LR, Beland MD, et al. Preoperative High-Resolution Ultrasound for the Assessment of Malignant Central Compartment Lymph Nodes in Papillary Thyroid Cancer. Thyroid (2015) 25(12):1351–4. doi: 10.1089/thy.2015.0176

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Chen H, Song A, Wang Y, He Y, Tong J, Di J, et al. BRAF(V600E) Mutation Test on Fine-Needle Aspiration Specimens of Thyroid Nodules: Clinical Correlations for 4600 Patients. Cancer Med (2021) 11(1):40–9. doi: 10.1002/cam4.4419

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Melo M, Gaspar da Rocha A, Batista R, Vinagre J, Martins MJ, Costa G, et al. TERT, BRAF, and NRAS in Primary Thyroid Cancer and Metastatic Disease. J Clin Endocrinol Metab (2017) 102(6):1898–907. doi: 10.1210/jc.2016-2785

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Machens A, Holzhausen HJ, Dralle H. The Prognostic Value of Primary Tumor Size in Papillary and Follicular Thyroid Carcinoma. Cancer (2005) 103(11):2269–73. doi: 10.1002/cncr.21055

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Machens A, Hinze R, Thomusch O, Dralle H. Pattern of Nodal Metastasis for Primary and Reoperative Thyroid Cancer. World J Surg (2002) 26(1):22–8. doi: 10.1007/s00268-001-0176-3

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Ha J, Lee J, Jo K, Han JS, Kim MH, Jung CK, et al. Calcification Patterns in Papillary Thyroid Carcinoma Are Associated With Changes in Thyroid Hormones and Coronary Artery Calcification. J Clin Med (2018) 7(8):183. doi: 10.3390/jcm7080183

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Bu Kyung Kim YSC, Kwon HJ, Lee JS, Heo JJ, Han YJ, Park Y-H, et al. Relationship Between Patterns of Calcification in Thyroid Nodules and Histopathologic Findings. Endocrine J (2012) 60(2):150–60. doi: 10.1507/endocrj.ej12-0294

CrossRef Full Text | Google Scholar

25. Feng Y, Min Y, Chen H, Xiang K, Wang X, Yin G. Construction and Validation of a Nomogram for Predicting Cervical Lymph Node Metastasis in Classic Papillary Thyroid Carcinoma. J Endocrinol Invest (2021) 44(10):2203–11. doi: 10.1007/s40618-021-01524-5

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Tian X, Song Q, Xie F, Ren L, Zhang Y, Tang J, et al. Papillary Thyroid Carcinoma: An Ultrasound-Based Nomogram Improves the Prediction of Lymph Node Metastases in the Central Compartment. Eur Radiol (2020) 30(11):5881–93. doi: 10.1007/s00330-020-06906-6

PubMed Abstract | CrossRef Full Text | Google Scholar