A nomogram based on clinicopathological and ultrasound characteristics to predict central neck lymph node metastases in papillary thyroid cancer

Abstract

Purpose:

Papillary thyroid cancer (PTC) has grown rapidly in prevalence over the past few decades, and central neck lymph node metastasis (CNLNM) is associated with poor prognoses. However, whether to carry out preventive central neck lymph node dissection (CNLND) is still controversial. We aimed to construct a prediction model of CNLNM to facilitate making clinical surgical regimens.

Methods:

A total of 691 patients with PTC between November 2018 and December 2021 were included in our study. Univariate and multivariate analyses were performed on basic information and clinicopathological characteristics, as well as ultrasound characteristics (American College of Radiology (ACR) scores). The prediction model was constructed and performed using a nomogram, and then discriminability, calibrations, and clinical applicability were evaluated.

Results:

Five variables, namely, male, age >55 years, clinical lymph node positivity, tumor size ≥1 cm, and ACR scores ≥6, were independent predictors of CNLNM in the multivariate analysis, which were eventually included to construct a nomogram model. The area under the curve (AUC) of the model was 0.717, demonstrating great discriminability. A calibration curve was developed to validate the calibration of the present model by bootstrap resampling, which indicated that the predicted and actual values were in good agreement and had no differentiation from the ideal model. The decision curve analysis (DCA) indicated that the prediction model has good clinical applicability.

Conclusions:

Our non-invasive prediction model combines ACR scores with clinicopathological features presented through nomogram and has shown good performance and application prospects for the prediction of CNLNM in PTCs.

1 Introduction

The prevalence of thyroid cancer has dramatically increased in the past few decades, among which papillary thyroid cancer (PTC) accounts for the largest proportion of approximately 80%–90% (1, 2). Although PTC has a great prognosis and a low risk, it has a high central neck lymph node metastasis (CNLNM) rate of approximately 20%–90%, which is closely related to the poor postoperative recurrence rate and survival outcome (3, 4). As we all know, central neck lymph node dissection (CNLND) facilitates reducing postoperative recurrence and avoiding a second surgery, but unnecessary prophylactic lymph node dissection may also lead to many other complications, such as permanent hypoparathyroidism and permanent laryngeal nerve damage (5). Accordingly, the prediction of CNLNM is critical for the accurate dissection of CNLNM, the formulation of a preoperative treatment plan, and postoperative personalized follow-up and management plan.

Ultrasound and computed tomography (CT) examination can predict CNLNM of thyroid cancer to a certain extent, but a recent meta-analysis indicated that the specificity and sensitivity were only 0.95 and 0.35, respectively, for ultrasound and 0.46 and 0.88, respectively, for CT in predicting CNLNM (6). Although fine-needle aspiration (FNA) cytology is the gold standard for the diagnosis of PTC and CNLNM (7), the lesion may not be obtained by puncture, and non-invasive examinations are more acceptable to patients. Accordingly, preoperative evaluation of CNLNM is still a crucial clinical challenge, and a new and high-efficiency prediction model is urgently needed to predict the risk of CNLNM.

Thyroid Imaging Reporting and Data System (TI-RADS), a hazard classification of thyroid nodules by scoring their characteristics, presented by the American College of Radiology (ACR) in 2017, provides greater accuracy in the diagnosis of PTC (8). Although the ACR score has been proposed to assess the malignancy of thyroid nodules, several researchers in recent years have suggested that the ACR score or ACR classification is an independent risk factor and can be used as a prediction model for CNLNM in PTC (9–13). For instance, Park et al. proposed that higher ACR TI-RADS scores were a predictor of CNLNM for small PTC (10–20 mm) (10). Nevertheless, all studies did not construct a nomogram when identifying risk factors for CNLNM by single- and multi-factor analyses.

According to our knowledge, this study is the first to construct a prediction model to draw a nomogram based on the ACR score for CNLNM. The purpose of this research is to develop a nomogram based on ACR score, cervical ultrasound features, and clinicopathological characteristics for the prediction of CNLNM in PTCs preoperatively and to guide clinical therapeutic regimen, surgical program, and postoperative follow-up plan.

2 Materials and methods

2.1 Patients and study design

This study is a retrospective research that was approved by the Clinical Research Center. All records between November 2018 and December 2021 were derived from Zhujiang Hospital, Southern Medical University, with permission of the research ethics committee and verbal consent of participants. A total of 691 patients enrolled in this research underwent thyroidectomy and CNLND and were diagnosed with PTC by pathology after surgery. Any participant who met the following criteria was excluded: 1) participants with a previous history of thyroid ablation or thyroid surgery; 2) participants with other malignant tumors, such as gastric cancer and medullary thyroid cancer; 3) participants with incomplete clinicopathological and ultrasound data; 4) patients with unclear ultrasound images and foci that cannot be diagnosed and evaluated.

2.2 Sociodemographic information

General information on the following was collected: age (<55 and ≥55 years), sex (female and male), and body mass index (BMI, <25 and ≥25 kg/m²).

2.3 Ultrasonography and ACR TI-RADS scoring procedure

A routine ultrasound examination was performed before surgery to measure the thyroid nodules and neck lymph nodes to determine the subsequent surgical plan. Ultrasonography for thyroid was performed using RESONA 70B (Mindray, Shenzhen, China) or GE Logiq 9, ARIETTA 850 (Hitachi, Tokyo, Japan) equipped with a linear array transducer in 5–20 MHz or 5–13 MHz. Ultrasound images were saved for subsequent data analysis. Two doctors who specialize in the endocrine system with at least 10 years of work experience retrospectively analyzed the images without knowledge of the basic and clinicopathological information, and disputes were resolved by consensus through discussion. The latest ultrasound images were included if the patient has several preoperative ultrasound images. The appearance of ≥2 suspected lesions in the thyroid was considered to be multifocal. The node with the highest ACR TI-RADS score was enrolled in this research in the case of multiple nodes. Clinically node-negative (cN0) PTC refers to the absence of all typical cervical ultrasound features, including microcalcifications, peripheral blood flow features, round lymph nodes, necrosis, unclear lymphatic portal structure, and focal hyperechoic. Diffuse lesions were confirmed by ultrasound reports, which may be chronic lymphocytic thyroiditis.

Ultrasound parameters of thyroid nodules included echo (very hypoechoic, 3 points; hypoechoic, 2 points; isoechoic or hyperechoic, 1 point; anechoic, 0 points), composition (solid or almost completely solid, 2 points; mixed solid and cystic, 1 point; spongiform, cystic, or almost completely cystic, 0 points), shape (taller-than-wide, 3 points; wider-than-tall, 0 points), echogenic foci (punctate echogenic foci, 3 points; peripheral calcifications, 2 points; macrocalcifications, 1 point; large comet-tail artifacts or none, 0 points), and margin (extra-thyroidal extension, 3 points; irregular or lobulated, 2 points; smooth, ill-defined, 0 points). All points were calculated using the 2017 ACR TI-RADS criteria (8).

2.4 Surgery program and pathological examination

All patients were subjected to total thyroidectomy or lobectomy combined with ipsilateral CNLND if the diagnosis of PTC was verified by preoperative FNA or intraoperative frozen pathology. Aside from this, bilateral central and lateral lymph node excision was carried out in cases of positive CNLNM in preoperative ultrasound, preoperative FNA, or intraoperative frozen pathology. FNA was not routinely performed before surgery. Postoperative thyroid and lymph node specimens were paraffin-embedded and used for H&E staining and immunohistochemical staining to clarify the pathological diagnosis. All postoperative specimens obtained were checked retrospectively by two licensed pathologists from Zhujiang Hospital, Southern Medical University. Pathological characteristics included tumor size, tumor pathological type, extrathyroidal extension (including invasion of capsular), multifocality, mutation status of BRAFV600E (mutant, wild type, and not performed), and CNLNM (areas and quantities). Diagnosis of Hashimoto’s thyroiditis was based on the pathology report. The bilateral tumor was described as the existence of cancer in the two thyroid lobes. The diameter of a node with the highest ACR score of multiple lesions was regarded as the maximum size.

2.5 Statistical analysis

Groups of CNLNM were divided according to the CNLNM status, and the basic and clinical characteristics were listed. Categorical variables are described as the frequency with percentage. Continuous variables are shown as the mean with standard deviation. Odds ratios and 95% confidence intervals (95% CIs) were initially used to describe the potential factors associated with CNLNM in the univariate analysis. Variables with an overall p < 0.05 in single-factor analysis were further enrolled in the multivariable analysis for independent predictor selection using the stepwise method. CNLNM was then predicted using a nomogram derived from the multivariate analysis. Then, this nomogram was evaluated using the receiver operating characteristic (ROC) curve, calibration curves, and decision curve analysis (DCA) developed from two internal cohorts (n = 691) produced by bootstrap resampling.

Statistical analyses were carried out using SAS 9.4 (SAS Institute Inc.) and R software (version 4.2.1, MathSoft Inc.). The p-value was based on two-sided tests. p < 0.05 was considered statistically significant.

3 Results

The study included 691 PTC patients, of which 377 (54.6%) had positive nodes by postoperative pathology. Further, 510 microcarcinomas were included, of which 218 (42.75%) developed CNLNM. In the CNLNM-positive group, 43 (11.4%) were older than 55 years, 117 (31.2%) were male, and 233 (62.3%) had a CN0. In the CNLNM group, 95 (25.2%) had a total ACR score greater than 10, and 266 (70.6%) had a score of 6 to 10, which was higher than that of the CNLNM-negative group (p < 0.001). The basic and clinicopathological features of the PTC are listed in Table 1.

3.1 Univariable and multivariable analysis

Univariable and multivariable analyses were conducted to find out risk factors significantly associated with CNLNM (Table 2), including variables of sex, age, BMI, CN0, Hashimoto’s thyroiditis, tumor size, BRAFV600E mutation, multifocality, tumor position, bilateral tumor, diffuse disease, extrathyroidal extension, and ACR scores. The outcome of the univariable analysis suggested that age, sex, CN0, tumor size, BRAFV600E mutation, multifocality, bilateral tumor, extrathyroidal extension, and ACR scores were significantly correlated with CNLNM (all p < 0.05). A multivariable analysis was further performed, and it showed that ACR scores ≥6 (6 ≤ ACR scores ≤ 10: OR = 3.65, 95% CI: 1.80–7.39, p < 0.001; ACR scores ≥11: OR = 4.37, 95% CI: 2.04–9.36, p < 0.001), clinical lymph node positivity (OR = 2.74, 95% CI: 1.84–4.07, p < 0.001), tumor size ≥1 cm (OR = 3.23, 95% CI: 2.20–4.74, p < 0.001), age > 55 years (OR = 0.52, 95% CI: 0.32–0.84, p = 0.008), and male (OR = 1.64, 95% CI: 1.12–2.41, p = 0.012) were independent risk factors of CNLNM (Table 2).

3.2 Establishment and validation of the nomogram

According to the final results of multivariable analysis, variables of age, sex, tumor size, CN0, and ACR score were eventually included to develop an easy-to-understand nomogram model (Figure 1), which could facilitate the prediction of CNLNM in PTCs. The area under the curve (AUC) of the ROC was 0.717 (Figure 2A). The established nomogram model was further evaluated using two independent cohorts produced by resampling bootstrap analysis, and the AUC values were 0.716 and 0.726 (Figures 2B, C). Three calibration curves were also developed to validate the calibration of the present model by bootstrap resampling (Figure 3). Great concordance was shown in the curves (ideal, apparent, and bias-corrected lines) derived from the training model (mean absolute error (MAE) = 0.010, Figure 3A) and two internal validation cohorts (MAE = 0.018 in Figure 3B; MAE = 0.019 in Figure 3C). The DCA was further conducted to evaluate the clinical practicability of the model in predicting CNLNM. The DCA curves showed that when the threshold probability of three cohorts was between 0.2 and 0.9, our prediction model was better than a no-treat or all-treat regimen.

Figure 1

Figure 2

Figure 3

4 Discussion

The increasing incidence of thyroid cancer, as we all know, is partly due to the advances in technology that allow for more accurate screening techniques now (14). However, along with the increased detection rates that now exist, some unnecessary surgeries have increased, and the resulting postoperative complications, such as hypothyroidism, hypoparathyroidism, and damage to the cervical laryngeal nerve, have been shown to cause lifelong inconvenience to patients (15–19). Hence, a judgment of lymph node status is essential for the decision of CNLND.

After reviewing previous studies that assessed risk factors for predicting CNLNM in PTCs, it was found that only a small number of studies have used total ACR scores. Different from previous guidelines, the 2017 ACR TI-RADS guidelines provide detailed scoring for each node in addition to classifying the nodes (TR1-TR5) only. A study suggested that an increased ACR TI-RADS score was related to CNLNM in both TR5 (diameter of tumor ≥10 mm) PTCs and TR4 PTCs13. Several other articles also believe that the ACR score is an independent risk factor of CNLNM10–12. However, none of them developed an easy-to-understand nomogram for the prediction model. From this, a nomogram based on ACR score, cervical ultrasound characteristics, and clinicopathological characteristics was established in the present study.

Risk factors of CNLNM were evaluated by single-factor and multivariate analyses in PTCs in our research. Sex, age, CN0, tumor size, BRAFV600E mutation, multifocality, bilateral tumor, extrathyroidal extension, and ACR scores were correlated with CNLNM in the univariate analysis. ACR scores ≥6 (6 ≤ ACR scores ≤ 10, p < 0.001; ACR scores ≥11, p < 0.001), clinical lymph node positivity (p < 0.001), age > 55 years (p = 0.008), female (p = 0.012), and tumor size ≥1 cm (p < 0.001) were independent risk factors of CNLNM.

In our study, ACR scores were divided into three groups (ACR scores ≤ 5, 6 ≤ ACR scores ≤ 10, and ACR scores ≥ 11), and in the multifactorial regression analysis, we found that ACR scores of 6–10 and ACR scores ≥11 were independent predictors for CNLNM for PTCs. We can also consider an ACR of 6–10 as 3.65 times and an ACR ≥ 11 as 4.37 times the risk of CNLNM in the group (ACR ≤ 5). Therefore, for PTC patients with an ACR score of ≥6, we should increase the emphasis on preoperative screening of cervical lymph nodes and encourage more preoperative examinations to help surgeons more accurately identify metastatic lymph nodes preoperatively for a more scientific surgical protocol, and more active follow-up is also needed.

Some studies concluded that tumor size ≥1 cm is more likely to develop CNLNM than tumor size <1 cm (20–22), although several studies have not segmented tumors by 1 cm (23–29), and some studies concluded that larger tumor volume is more likely to develop CNLNM (30, 31). In sum, CNLNM is generally more likely to occur with larger tumors. In our research, we discovered that tumor size ≥1 cm was an independent risk factor of CNLNM, which is consistent with previous findings (20–22). It is possible that the larger the tumor diameter, the more advanced the tumor growth, so CNLNM is more likely to occur. However, even if the tumor size is small, CNLNM may occur due to its strong invasiveness (25). We can only confirm the aggressiveness of the tumor by postoperative pathology at present, so we should put importance on the preoperative assessment of thyroid nodules and CNLNM regardless of the size of the tumor.

Many studies suggested that younger-age patients are more likely to develop CNLNM (20–22), and in line with these previous studies, we also found that patients aged <55 years are more likely to develop CNLNM, which may be due to the rapid rate of cell growth and renewal in younger people. However, a study concluded that age < 55 years is not a risk factor for central lymph node metastasis in PTC (32). Therefore, the relationship between CNLNM and age still needs more verification by more studies.

Although previous studies have suggested that women are more likely to develop thyroid cancer than men, some experts believe that this is due to female estrogen (33). Nevertheless, men are more likely to develop CNLNM once their immune defenses have been breached. The male sex was believed to be an independent risk factor in the previous research (20–23). However, the differences in the extent of thyroid cancer in men versus women may not only be related to biological factors. Several studies pointed out that a higher rate of CNLNM in men can also be attributed to delayed presentation and delayed surgery in men, as body habitus and neck size in men always lead to less palpable lymph nodes than in women (34–39). In our research, we also found that men were more likely to develop CNLNM.

For easier application and visualization of the results of predicting CNLNM risk, an easy-to-use nomogram based on ACR scores and clinicopathological characteristics that we screened out was constructed, which allowed us to quantify the impact of multiple independent variables on the prediction model. Interestingly, the ACR score in our predictive nomogram was the largest contributor to the total score. The nomogram is used as follows: first, draw a line through the values of the variables perpendicular to the axis of points and find the spot of each variable on it. Second, sum the points of all the variables to obtain the total points. Finally, find the spot of the total point axis and draw a line perpendicular to the risk axis, which is the risk of the CNLNM. Based on the nomogram, the AUC and C-index were 71.7% and 0.717, respectively, and the AUC values of two independent cohorts produced by resampling bootstrap analysis were 0.716 and 0.726, which demonstrated the great discriminability of this present prediction model. Calibration curves developed from the training model and two internal validation cohorts showed a good concordance between the predicted value and the actual value and no differentiation from the ideal model. The DCA curves showed great clinical practicability of the model in predicting CNLNM.

Although we obtained a prediction model with good performance, there are still several limitations. First, the present study was single-center retrospective research, and inevitably, confounding factors were present. More multicenter prospective studies are therefore required in the future. Second, since ultrasound features rely heavily on the experience of the doctor, while our ultrasound images are assessed by an experienced thyroid surgeon, subjective bias was inclined to occur, and thus, the objectivity of the findings was easily compromised. Additionally, it is well known that CNLNM is more prone to metastasis and not easily detected by ultrasound due to the interference of the thyroid gland and tracheal cartilage, so the present study only analyzed CNLNM due to the limited number of cases with lateral neck lymph node metastasis. Thus, neck lymph node metastasis in the lateral region demands more analysis urgently in the future. Finally, although our model has demonstrated good diagnostic performance through internal validation by bootstrapping analysis, overfitting of the model may also have occurred, and further external validation is still required.

5 Conclusions

Our model is a non-invasive prediction tool that combines ACR scores with clinicopathological features presented through nomogram, which can predict CNLNM in PTCs and has shown good performance and prospects for application.

References

• 1

  WangJYuFShangYPingZLiuL. Thyroid cancer: incidence and mortality trends in China, 2005-2015. Endocrine (2020) 68:163–73. doi: 10.1007/s12020-020-02207-6

  • CrossRef
  • Google Scholar

• 2

  WangYDuanYLiHYueKLiuJLaiQet al. Detection of thyroglobulin in fine-needle aspiration for diagnosis of metastatic lateral cervical lymph nodes in papillary thyroid carcinoma: A retrospective study. Front Oncol (2022) 12:909723. doi: 10.3389/fonc.2022.909723

  • CrossRef
  • Google Scholar

• 3

  DaviesLHoangJK. Thyroid cancer in the USA: current trends and outstanding questions. Lancet Diabetes Endocrinol (2021) 9:11–2. doi: 10.1016/S2213-8587(20)30372-7

  • CrossRef
  • Google Scholar

• 4

  WangLYGanlyI. Nodal metastases in thyroid cancer: prognostic implications and management. Future Oncol Lond Engl (2016) 12:981–94. doi: 10.2217/fon.16.10

  • CrossRef
  • Google Scholar

• 5

  ZhaoY-ZHeN-AYeX-JJinFLiM-XJiangX. Analysis of risk factors associated with central lymph node metastasis in papillary thyroid carcinoma with cT1N0 stage. Front Endocrinol (2022) 13:880911. doi: 10.3389/fendo.2022.880911

  • CrossRef
  • Google Scholar

• 6

  YangJZhangFQiaoY. Diagnostic accuracy of ultrasound, CT and their combination in detecting cervical lymph node metastasis in patients with papillary thyroid cancer: a systematic review and meta-analysis. BMJ Open (2022) 12:e051568. doi: 10.1136/bmjopen-2021-051568

  • CrossRef
  • Google Scholar

• 7

  QiaoD-HHeX-MDengXJiY-CYangHChengLet al. Aberrant expression of five miRNAs in papillary thyroid carcinomas. J Clin Lab Anal (2021) 35:e23907. doi: 10.1002/jcla.23907

  • CrossRef
  • Google Scholar

• 8

  TesslerFNMiddletonWDGrantEGHoangJKBerlandLLTeefeySAet al. ACR thyroid imaging, reporting and data system (TI-RADS): white paper of the ACR TI-RADS committee. J Am Coll Radiol JACR (2017) 14:587–95. doi: 10.1016/j.jacr.2017.01.046

  • CrossRef
  • Google Scholar

• 9

  HuTLiZPengCHuangLLiHHanXet al. Nomogram to differentiate benign and Malignant thyroid nodules in the American College of Radiology Thyroid Imaging Reporting and Data System level 5. Clin Endocrinol (Oxf) (2023) 98:249–58. doi: 10.1111/cen.14824

  • CrossRef
  • Google Scholar

• 10

  ParkHMLeeJKwakJYParkVYRhoMLeeMet al. Using ultrasonographic features to predict the outcomes of patients with small papillary thyroid carcinomas: a retrospective study implementing the 2015 ATA patterns and ACR TI-RADS categories. Ultrason Seoul Korea (2022) 41:298–306. doi: 10.14366/usg.21097

  • CrossRef
  • Google Scholar

• 11

  ZhongMZhangZXiaoYHeYChenYHuangWet al. The predictive value of ACR TI-RADS classification for central lymph node metastasis of papillary thyroid carcinoma: A retrospective study. Int J Endocrinol (2022) 2022:4412725. doi: 10.1155/2022/4412725

  • CrossRef
  • Google Scholar

• 12

  YangJ-RSongYChangS-JShiL-L. Prediction of central compartment nodal metastases in papillary thyroid cancer using TI-RADS score, blood flow, and multifocality. Acta Radiol Stockh Swed 1987 (2022) 63:1374–80. doi: 10.1177/02841851211041811

  • CrossRef
  • Google Scholar

• 13

  ShangguanRKanGYangJBaoL. Predicting papillary thyroid carcinoma cervical lymph node metastases: an algorithm using the American College of Radiology Thyroid Imaging, Reporting and Data System. Acta Radiol Stockh Swed 1987 (2023) 64:101–7. doi: 10.1177/02841851211069779

  • CrossRef
  • Google Scholar

• 14

  XuYChenJYangZXuL. Identification of RNA Expression Profiles in Thyroid Cancer to Construct a Competing Endogenous RNA (ceRNA) Network of mRNAs, Long Noncoding RNAs (lncRNAs), and microRNAs (miRNAs). Med Sci Monit Int Med J Exp Clin Res (2019) 25:1140–54. doi: 10.12659/MSM.912450

  • CrossRef
  • Google Scholar

• 15

  KannanSRajuNKekatpureVChandrasekharNHPillaiVKeshavamurthyARet al. Improving bethesda reporting in thyroid cytology: A team effort goes a long way and still miles to go…. Indian J Endocrinol Metab (2017) 21:277–81. doi: 10.4103/ijem.IJEM_350_16

  • CrossRef
  • Google Scholar

• 16

  PolistenaAMonacelliMLucchiniRTriolaRContiCAveniaSet al. Surgical morbidity of cervical lymphadenectomy for thyroid cancer: A retrospective cohort study over 25 years. Int J Surg Lond Engl (2015) 21:128–34. doi: 10.1016/j.ijsu.2015.07.698

  • CrossRef
  • Google Scholar

• 17

  ArditoGRevelliLPolistenaALucchiniRGiustozziEGuidiMLet al. Complications of neck dissections in papillary thyroid carcinoma: A modified procedure to reduce parathyroid morbidity. Vivo Athens Greece (2016) 30:303–8.

  • Google Scholar

• 18

  McMullenCRockeDFreemanJ. Complications of bilateral neck dissection in thyroid cancer from a single high-volume center. JAMA Otolaryngol– Head Neck Surg (2017) 143:376–81. doi: 10.1001/jamaoto.2016.3670

  • CrossRef
  • Google Scholar

• 19

  De NapoliLMatroneAFavillaKPiaggiPGalleriDAmbrosiniCEet al. Role of prophylactic central compartment lymph node dissection on the outcome of patients with papillary thyroid carcinoma and synchronous ipsilateral cervical lymph node metastases. Endocr Pract Off J Am Coll Endocrinol Am Assoc Clin Endocrinol (2020) 26:807–17. doi: 10.4158/EP-2019-0532

  • CrossRef
  • Google Scholar

• 20

  FengYMinYChenHXiangKWangXYinG. Construction and validation of a nomogram for predicting cervical lymph node metastasis in classic papillary thyroid carcinoma. J Endocrinol Invest (2021) 44:2203–11. doi: 10.1007/s40618-021-01524-5

  • CrossRef
  • Google Scholar

• 21

  HuangYHuangZCaiHZhugeLWangSYanDet al. Evaluation of serum B7-H3 expression, ultrasound and clinical characteristics to predict the risk of cervical lymph node metastases in papillary thyroid carcinoma by nomogram. J Clin Lab Anal (2023) 37:e24811. doi: 10.1002/jcla.24811

  • CrossRef
  • Google Scholar

• 22

  ChangQZhangJWangYLiHDuXZuoDet al. Nomogram model based on preoperative serum thyroglobulin and clinical characteristics of papillary thyroid carcinoma to predict cervical lymph node metastasis. Front Endocrinol (2022) 13:937049. doi: 10.3389/fendo.2022.937049

  • CrossRef
  • Google Scholar

• 23

  YeJFengJ-WWuW-XHuJHongL-ZQinA-Cet al. Papillary thyroid microcarcinoma: A nomogram based on clinical and ultrasound features to improve the prediction of lymph node metastases in the central compartment. Front Endocrinol (2021) 12:770824. doi: 10.3389/fendo.2021.770824

  • CrossRef
  • Google Scholar

• 24

  FengJ-WWuW-XQiG-FHongL-ZHuJLiuS-Yet al. Nomograms based on sonographic and clinicopathological characteristics to predict lateral lymph node metastasis in classic papillary thyroid carcinoma. J Endocrinol Invest (2022) 45:2043–57. doi: 10.1007/s40618-022-01825-3

  • CrossRef
  • Google Scholar

• 25

  NieXTanZGeMJiangLWangJZhengC. Risk factors analyses for lateral lymph node metastases in papillary thyroid carcinomas: a retrospective study of 356 patients. Arch Endocrinol Metab (2016) 60:492–9. doi: 10.1590/2359-3997000000218

  • CrossRef
  • Google Scholar

• 26

  TaşkınOÇArmutluAAğcaoğluOPekerÖTerzioğluTDemirkolMOet al. Tumor border pattern and size help predict lymph node status in papillary microcarcinoma: A clinicopathologic study. Ann Diagn Pathol (2020) 48:151592. doi: 10.1016/j.anndiagpath.2020.151592

  • CrossRef
  • Google Scholar

• 27

  CelikHAkgulOYıldızBDSaylamBTezM. Predictors of central lymph node metastasis in papillary thyroid cancer. Ann Ital Chir (2017) 88:193–7.

  • Google Scholar

• 28

  KimSKParkIWooJ-WLeeJHChoeJ-HKimJ-Het al. Predictive factors for lymph node metastasis in papillary thyroid microcarcinoma. Ann Surg Oncol (2016) 23:2866–73. doi: 10.1245/s10434-016-5225-0

  • CrossRef
  • Google Scholar

• 29

  JiYBYooHSSongCMParkCWLeeCBTaeK. Predictive factors and pattern of central lymph node metastasis in unilateral papillary thyroid carcinoma. Auris Nasus Larynx (2016) 43:79–83. doi: 10.1016/j.anl.2015.09.005

  • CrossRef
  • Google Scholar

• 30

  ZhouLYaoJOuDLiMLeiZWangLet al. A multi-institutional study of association of sonographic characteristics with cervical lymph node metastasis in unifocal papillary thyroid carcinoma. Front Endocrinol (2022) 13:965241. doi: 10.3389/fendo.2022.965241

  • CrossRef
  • Google Scholar

• 31

  ParkKNKangKYHongHSJeongH-SLeeSW. Predictive value of estimated tumor volume measured by ultrasonography for occult central lymph node metastasis in papillary thyroid carcinoma. Ultrasound Med Biol (2015) 41:2849–54. doi: 10.1016/j.ultrasmedbio.2015.02.018

  • CrossRef
  • Google Scholar

• 32

  WangZChangQZhangHDuGLiSLiuYet al. A clinical predictive model of central lymph node metastases in papillary thyroid carcinoma. Front Endocrinol (2022) 13:856278. doi: 10.3389/fendo.2022.856278

  • CrossRef
  • Google Scholar

• 33

  DerwahlMNiculaD. Estrogen and its role in thyroid cancer. Endocr Relat Cancer (2014) 21:T273–283. doi: 10.1530/ERC-14-0053

  • CrossRef
  • Google Scholar

• 34

  WangFZhaoSShenXZhuGLiuRViolaDet al. BRAF V600E confers male sex disease-specific mortality risk in patients with papillary thyroid cancer. J Clin Oncol Off J Am Soc Clin Oncol (2018) 36:2787–95. doi: 10.1200/JCO.2018.78.5097

  • CrossRef
  • Google Scholar

• 35

  SuteauVMunierMBrietCRodienP. Sex bias in differentiated thyroid cancer. Int J Mol Sci (2021) 22:12992. doi: 10.3390/ijms222312992

  • CrossRef
  • Google Scholar

• 36

  GermanoASchmittWAlmeidaPMateus-MarquesRLeiteV. Ultrasound requested by general practitioners or for symptoms unrelated to the thyroid gland may explain higher prevalence of thyroid nodules in females. Clin Imaging (2018) 50:289–93. doi: 10.1016/j.clinimag.2018.05.003

  • CrossRef
  • Google Scholar

• 37

  BertakisKDAzariRHelmsLJCallahanEJRobbinsJA. Gender differences in the utilization of health care services. J Fam Pract (2000) 49:147–52.

  • Google Scholar

• 38

  VaccarellaSDal MasoLLaversanneMBrayFPlummerMFranceschiS. The impact of diagnostic changes on the rise in thyroid cancer incidence: A population-based study in selected high-resource countries. Thyroid Off J Am Thyroid Assoc (2015) 25:1127–36. doi: 10.1089/thy.2015.0116

  • CrossRef
  • Google Scholar

• 39

  YooJAhnHSKimSJParkSHSeoMChongS. Evaluation of diagnostic performance of screening thyroid ultrasonography and imaging findings of screening-detected thyroid cancer. Cancer Res Treat (2018) 50:11–8. doi: 10.4143/crt.2016.600

  • CrossRef
  • Google Scholar