Skip to main content

ORIGINAL RESEARCH article

Front. Endocrinol., 30 October 2023
Sec. Thyroid Endocrinology

Risk factor analysis and prediction model for papillary thyroid carcinoma with lymph node metastasis

Juerong LuJuerong Lu1Jintang LiaoJintang Liao1Yunhao ChenYunhao Chen1Jie LiJie Li1Xinyue HuangXinyue Huang1Huajun Zhang,,Huajun Zhang1,2,3Bo Zhang,*Bo Zhang1,4*
  • 1Department of Ultrasonic Imaging, Xiangya Hospital, Central South University, Changsha, Hunan, China
  • 2Department of Oncology, National Health Commission of the People's Republic of China (NHC) Key Laboratory of Cancer Proteomics, Xiangya Hospital, Central South University, Changsha, Hunan, China
  • 3Laboratory of Structural Biology, Xiangya Hospital, Central South University, Changsha, Hunan, China
  • 4Molecular Imaging Research Center of Central South University, Changsha, Hunan, China

Objective: We aimed to identify the clinical factors associated with lymph node metastasis (LNM) based on ultrasound characteristics and clinical data, and develop a nomogram for personalized clinical decision-making.

Methods: A retrospective analysis was performed on 252 patients with papillary thyroid carcinoma (PTC). The patient’s information was subjected to univariate and multivariate logistic regression analyses to identify risk factors. A nomogram to predict LNM was established combining the risk factors. The performance of the nomogram was evaluated using receiver operating characteristic (ROC) curve, calibration curve, cross-validation, decision curve analysis (DCA), and clinical impact curve.

Results: There are significant differences between LNM and non-LNM groups in terms of age, sex, tumor size, hypoechoic halo around the nodule, thyroid capsule invasion, lymph node microcalcification, lymph node hyperechoic area, peak intensity of contrast (PI), and area under the curve (AUC) of the time intensity curve of contrast (P<0.05). Age, sex, thyroid capsule invasion, lymph node microcalcification were independent predictors of LNM and were used to establish the predictive nomogram. The ROC was 0.800, with excellent discrimination and calibration. The predictive accuracy of 0.757 and the Kappa value was 0.508. The calibration curve, DCA and calibration curve demonstrated that the prediction model had excellent net benefits and clinical practicability.

Conclusion: Age, sex, thyroid capsule invasion, and lymph node microcalcification were identified as significant risk factors for predicting LNM in patients with PTC. The visualized nomogram model may assist clinicians in predicting the likelihood of LNM in patients with PTC prior to surgery.

1 Introduction

According to the latest global cancer statistics, thyroid cancer is responsible for 586,000 cases worldwide and ranks in 9th place for incidence in 2020 with an estimated annual mortality of 44,000 (1). PTC is the most common histologic subtype of thyroid cancer and has a rapid increase in incidence globally over recent decades (2). PTC are prone to early metastasizes to lymph nodes with a reported incidence from 30 to 80% at initial diagnosis. LNM has been confirmed to be a significant risk factor affecting the survival rate of PTC patients, and is closely related to distant metastasis, recurrence, and poor prognosis (3).

Currently, surgery remains the primary treatment option for PTC, however, controversy persists regarding lymph node management (4). LNM preoperatively strongly correlates with distant metastasis, high locoregional recurrence, and enhanced death risk (5). Early implementation of appropriate surgical approaches and complete dissection of metastatic lymph nodes can reduce the probability of a second operation, and improve the prognosis and survival rate of the patient (6). Nevertheless, by reducing unnecessary prophylactic lymph node dissection, the occurrence of surgery-related complications can be minimized, such as hypoparathyroidism, recurrent laryngeal nerve damage, nerve injury to the voice, etc (7). Preoperative accurate assessment of whether there is lymph node metastasis in PTC is crucial for clinical decision-making and patient prognosis.

However, there are limited preoperative methods available for assessing LNM in PTC. US is the preferred diagnostic method for LNM with a high accuracy and sensitivity in identifying LNM (8). The additional use of contrast-enhanced ultrasound (CEUS) could provide qualitative and quantitative blood perfusion information and improve the diagnostic accuracy of thyroid nodules (9). Qualitative and quantitative Contrast-Enhanced Ultrasound analysis demonstrated a high diagnostic accuracy for the preoperative diagnosis of metastatic lymph nodes in patients with PTC. However, most sonographic features of metastatic lymph node are atypical, especially in early stages. Even under contrast-enhanced ultrasound, the early metastatic lymph nodes are not easy to detect (10). The assessment of LNM in patients with PTC before surgery presents significant difficulties and challenges for clinical radiologists.

Our study aimed to investigate risk factors associated with LNM and develop a user-friendly nomogram incorporating clinical and ultrasound features for quantitative prediction of LNM risk in PTC. The approach will aid clinicians in making informed therapeutic decisions, thereby improving patient prognosis.

2 Materials and methods

2.1 Participants

This prospective study was approved by the Ethics Committee of our Hospital (No. 202206140) and conducted in accordance with guidelines for good clinical practice. The study included 252 consecutive patients who underwent thyroidectomy and lymph node dissection for PTC at our institution from May 2019 to March 2022. Prior to surgery, all patients underwent preoperative serological examination, two-dimensional ultrasound (2D-US) and CEUS of the thyroid. The pathological results of LNM were served as the gold standard. The patients were classified into two groups based on the post-dissection pathological results of the lymph nodes. One group had lymph node metastases (LNM, n=139), while the other group did not (non-LNM, n=113).

The inclusion criteria were as follows: (a) Primary surgical resection was performed with a pathological diagnosis of PTC. (b) The patient underwent partial or total thyroidectomy and lymph node dissection as the surgical procedure. (c) Confirmation of lymph node metastasis was based on postoperative pathology; (d) Preoperative serological examination, 2D-US and CEUS of the thyroid were performed before surgery with complete clinicopathological information. (e) The date between serological examination, ultrasonography and surgery did not exceed one month.

The exclusion criteria were as follows: (a) Other types of thyroid cancer or PTC with other types of thyroid cancer was confirmed through surgical pathology. (b) History of previous neck surgery. (c) Age <18 years old.

2.2 Assessment of variables

The baseline clinical information of patients included age, sex (male, female). Biochemical indicators included free thyroid hormone (FT4), free triiodothyronine (FT3), hypersensitive thyroid stimulating hormone (TSH), thyroglobulin (TG), antithyroglobulin antibody (TGA), thyroid peroxidase antibody (TPOAB), parathyroid hormone (PTH), carcinoembryonic antigen (CEA). The results of FT4, FT3, TSH, TG, TGA, TPOAB, PTH, and CEA were detected within one month before surgery. According to the clinical experience, the thresholds set for FT4, FT3, TSH, TG, TGA, TPOAB, PTH, and CEA were as follows: FT4 12-22 pmol/L, FT3 2.8-7.1 pmol/L, TSH 0.27-4.2 uIU/ml, TG 3.5-77 ng/ml, TGA 0-115.00 IU/mL, TPOAB 0-34.00 IU/mL, PTH 15-65 pg/ml, CEA <5.0 ng/mL. (Table 1 for details).

TABLE 1
www.frontiersin.org

Table 1 Variables of patients’ clinical data in univariate analysis.

2D-US and CEUS indicators for thyroid nodule included the following characteristics: thyroid diffuse lesions, multifocality, tumor size, location, echogenicity, shape, aspect ratio (A/T, wider and tall ratio), boundary, margin, hypoechoic halo (Substantial, not the annular vessel surrounding the nodules), posterior attenuation, calcification, blood flow, thyroid capsule invasion, enhancement level, enhancement mode, enhancement uniformity and perfusion defect (Table 2 for details).

TABLE 2
www.frontiersin.org

Table 2 Variables of thyroid nodule 2D-US and CEUS in univariate analysis.

2D-US indicators of lymph node included the following characteristics: lymph node morphology, microcalcification, hyperechoic change, dermatomedullary boundary and blood flow. (Table 3 for details).

TABLE 3
www.frontiersin.org

Table 3 Variables of lymph node 2D-US in univariate analysis.

Qualitative and quantitative Contrast-Enhanced Ultrasound analysis of thyroid nodule included Base Intensity (BI), Arrival Time (AT), Time to Peak (TTP), Peak Intensity (PI), Ascending Slope (AS), Area Under Curve (AUC). (Table 4 for details).

2.3 US images acquisition and analysis

The US examination was independently performed by one of the two experienced radiologists with more than 20 years of experience in thyroid ultrasound diagnosis. All conventional US and CEUS examinations were performed by Resona R9 & 7S (Mindray Medical, Shenzhen, China) and Acuson Sequoia (Siemens, Erlangen, Germany). The ultrasound contrast agent used in thyroid contrast-enhanced ultrasound was the pure blood-pool contrast agent sulfur hexafluoride microbubble (SonoVue, Bracco, Italy).

Inclusion criteria for quantitative analysis of contrast-enhanced ultrasound of thyroid nodules: (a) CEUS recording lasted longer than one minute. (b) Thyroid nodule appeared in the observation screen completely. (c) Thyroid nodules were not significantly displaced.

Quantitative analysis of contrast-enhanced ultrasound: the region of interest (ROI) was delineated along the edge of the thyroid nodule, and the ROI was delineated no less than three times. Finally, the highest goodness of fit (GOF) was selected.

2.4 Statistical analysis

Statistical analysis was performed using SPSS 22 software and R software. Univariate and multivariate analyses were conducted for screening risk factors significantly associated with LNM. The statistically significant variables in the univariate analysis were subsequently included in the multivariate logistic regression analysis to construct a risk prediction model-nomogram. The discriminative power and consensus degree of our predictive model were assessed using the receiver operating characteristic (ROC) curve and the calibration curve. The training cohort (70% of patients) and the validation cohort (30% of patients) were randomly sampled for cross-validation. were used to evaluate The consistency of the results was assessed using the accuracy and Kappa values. Decision curve analysis (DCA) and clinical impact curve were applied to assess the feasibility and application value of the model in clinical practice.

3 Results

3.1 Patients’ clinical data univariate analysis

A total of 252 patients with PTC were enrolled in this study. They were divided into LNM and non-LNM groups in accordance with LNM or not. In the LNM group, there were 90 female and 49 male patients. The age of the patients ranged from 18 to 63 years, with a mean age of 41 ± 10 years. In the non-LNM group, there were 93 female and 20 male patients. The age of the patients ranged from 22 to 70 years, with a mean age of 46 ± 12 years.

At univariate analysis, age (OR=0.958, P<0.001), sex (OR=0.395, P=0.002) were found as the potential risk factors associated with LNM in PTC patients. However, the p values of FT4, FT3, TSH, TG, TGA, TPOAB, PTH, and CEA were all higher than 0.05, and there was no significant difference between the two groups (Table 1).

3.2 Thyroid nodule 2D-US and CEUS univariate analysis

At univariate analysis, tumor size (OR=1.043, P=0.031), hypoechoic halo around the nodule (OR=2.364, P=0.009), thyroid capsule invasion (OR=4.040, P<0.001) were found as the potential risk factors associated with LNM in PTC patients (Table 2).

3.3 Lymph node 2D-US univariate analysis

At univariate analysis, lymph node microcalcification (OR=10.966, P<0.001), lymph node hyperechoic change (OR=9.625, P=0.031) were found as the potential risk factors associated with LNM in PTC patients (Table 3).

3.4 Qualitative and quantitative CEUS analysis of thyroid nodule univariate analysis

Quantitative analysis results will be affected by patient’s cough, speech, swallowing or other reasons which will lead to a significant displacement of the lesion. There were some of the videos that were less than a minute tall. Therefore, a total of 151 patients with PTC were included in the final quantitative analysis, including 82 patients with LNM (54%) and 69 patients without LNM (46%).

At univariate analysis, PI (OR=0.954, P=0.029), AUC (OR=0.999, P=0.011) were found as the potential risk factors associated with LNM in PTC patients (Table 4).

TABLE 4
www.frontiersin.org

Table 4 Variables of qualitative and quantitative CEUS analysis in univariate analysis.

3.5 Multivariate logistic regression analyses

The nine variables selected by univariate analysis (age, sex, size, hypoechoic halo around the nodule, thyroid capsule invasion, lymph node microcalcification, lymph node hyperechoic change, PI and AUC were entered into the regression model for multivariate logistic regression analysis. The results showed that age (P=0.024), sex (P=0.021), thyroid capsule invasion (P<0.001), and lymph node microcalcification (P=0.005) were found as risk factors highly associated with lymph node metastasis in PTC (Table 5).

TABLE 5
www.frontiersin.org

Table 5 Nine variables selected by univariate analysis in multivariate analysis.

3.6 Development of prediction model and nomogram for predicting LNM in PTC patients

Based on the results of the multivariate logistic regression analysis, the logistic prediction model was developed as follows: logit (P) = 2.380-1.037* sex-0.044 * age +1.627* thyroid capsule invasion +2.009* lymph node microcalcification (0 for male, 1 for female; 0 for no thyroid capsule invasion, 1 for thyroid capsule invasion; 0 for no lymph node microcalcification, 1 for lymph node microcalcification).

In order to facilitate the clinical application and make the model become convenient and visual, we use R software and “rms” software package to construct the nomogram and make various forms of nomogram. Finally, we apply and register the web nomogram: https://juenomogram.shinyapps.io/DynNomapp/ (Figures 1, 2).

FIGURE 1
www.frontiersin.org

Figure 1 (A) Concise nomogram for predicting LNM in PTC patients. (B) Dynamic nomogram for predicting LNM in PTC patients.

FIGURE 2
www.frontiersin.org

Figure 2 (A) Dynamic web-based nomogram for predicting LNM in PTC patients (screen shot 1). (B) Dynamic web-based nomogram for predicting LNM in PTC patients (screen shot 2).

3.7 Evaluation of clinical prediction model

The ROC curve was drawn by R language. The nomogram yielded the AUC (0.800; 95% CI: 0.747-0.853, p < 0.001) which indicated that the prediction model was good and had good clinical predictive ability.The calibration curve was drawn by bootstrap 1000 resampling method. The calibration curve showed that a mostly agreement between the predicted and actual results of the nomogram model, indicating that the calibration degree of the model was good (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3 (A) ROC curve analysis to predict LNM in PTC patients. (B) Calibration curve of the nomogram. The calibration curve illustrates the calibration of the nomogram in terms of the agreement between the predicted risk of LNM and the observed outcomes of LNM. The 45° dotted gray line represents a perfect prediction, and the solid black line represents the predictive performance of the nomogram. The solid line has a closer fit to the dotted line, which indicates better predictive accuracy of the nomogram.

3.8 Validation of clinical prediction model

The R language was used for operation. 176 cases (70%) were selected as the training set, 76 cases (30%) were selected as the validation set. The number of iterations was 1. The results were: accuracy =0.757, Kappa=0.508.

The accuracy is greater than 0.7, indicating that the prediction result is good. Kappa is used to measure the stability of the model. The higher the Kappa value, the better the consistency. When the Kappa value is between 0.4 and 0.6, the consistency is acceptable.

3.9 Evaluation of clinical decision-making ability and practical application

The DCA curve showed that the nomogram curve had a higher net benefit than the extreme curve and other ultrasound indexes, suggesting that the model had a good clinical application value for the occurrence of LNM in PTC. The clinical impact curve of this study suggested that the nomogram model had good clinical practicability (Figure 4). Representative examples of dynamic web-based nomogram prediction for LNM in PTC patients are shown in Figures 5, 6.

FIGURE 4
www.frontiersin.org

Figure 4 (A) Decision curve analysis for the nomogram. The y-axis measures the net benefit. The x-axis measures the threshold probability. (B) Clinical impact curve of prediction model.

FIGURE 5
www.frontiersin.org

Figure 5 (A–D) Title: Representative example of dynamic web-based nomogram 520 prediction for a 36-year-old female PTC patient with LNM. (A) US images showed the maximum diameter of hypoechoic thyroid nodule was 23.8mm in left lobe with thyroid capsule invasion. (B) US images showed there was microcalcification in the left lymph node. (C) Dynamic web-based nomogram predicted that the risk of LNM in this patient was 0.965. (D) Pathology showed PTC in left lobe with thyroid capsule invasion and LNM.

FIGURE 6
www.frontiersin.org

Figure 6 (A–D) Title: Representative example of dynamic web-based nomogram prediction for a 48-year-old female PTC patient without LNM. (A) US images showed the maximum diameter of hypoechoic thyroid nodule was 5.4mm in left lobe with no thyroid capsule invasion. (B) US images showed there was no microcalcification in the left lymph node. (C) Dynamic web-based nomogram predicted that the risk of LNM in this patient was 0.217. (D) Pathology showed PTC in left lobe with no thyroid capsule invasion and no LNM.

4 Discussion

PTC is among the most frequently occurring types of cancer (11). While the mortality rate of PTC is not high, LNM significantly impacts the patient prognosis. Early detection of lymph node metastasis allows for lymph node dissection, which improves patient survival rates (12). If there is no LNM, unnecessary operations can be reduced. In more extensive surgery, interventions are influenced by more severe complications. Complications may arise in a significant number of patients, even in the most skilled hands and in high-volume centers (13). Postoperative hypocalcemia is observed in up to one third of thyroidectomy patients (14). There are other complications, such as hypoparathyroidism, injury to the recurrent laryngeal nerve, cervical hematoma, and wound infection (13, 15). These complications result in a decrease in the patients’ quality of life and increased costs for healthcare systems (13).

The present study identified nine variables, including gender, age, size, hypoechoic halo around the nodule, invasion of the thyroid capsule, microcalcification in lymph nodes, hyperechoic change in lymph nodes, PI, and AUC, which are associated with a high volume of CLNM in PTMC patients during the binary univariate logistic regression analysis.

While the incidence of papillary thyroid cancer is higher in women compared to men, men have a great risk of LNM than women in this type of cancer (16). Currently, the specific mechanism is unknown, but relevant studies suggested that it may be related to gender genes and hormone (17). Although some studies have shown that BRAF mutation and estrogen were not equally distributed between the sexes, whether they affect the occurrence of LNM in PTC needs to be further confirmed (18). Men have a higher basal metabolic rate compared to women, leading to accelerated tumor cell growth. Consequently, men face a higher risk for thyroid capsule invasion and LNM.

Previous studies used the TNM stage of AJCC thyroid cancer was used as a reference, with a cut-off age of 45 or 55 years (19). However, some scholars argue that age is negatively correlated with LNM in PTC and that there is no optimal cut-off value for age (20). While the incidence of PTC increases with age, the rate of LNM decreases (21). Our hypothesis is that the younger individuals have higher the metabolic rate and provide a more favorable environment for tumor cells. Therefore, younger patients with PTC are more likely to have LNM than older patients with PTC. In this study, age was included in the quantitative measures. The prediction model could give higher accurate individual LNM risk probability for each patient at different ages. Chan et al. (22). found that a distance from the capsule < 1.9 mm was associated with CLNM in cN0 PTC patients. A Meta-Analysis of clinicopathological significance of PTC located in the isthmus showed that isthmus PTC was associated with an increased risk of central LN metastasis (23). We conclude this for two reasons: first, the closer the lesion is to the thyroid capsule, the shorter the path for cancer cells to metastasize through the lymph nodes; second, the faster the cancer cells grow, the easier it is for them to invade the thyroid capsule, indicating that the malignant biological behavior of cancer cells is active.

Microcalcification is usually defined as calcification less than 1mm in diameter, which is a characteristic feature of PTC (24). This may be attributed to tumor thrombus from lymphatic or blood vessels. Or the glycoproteins and mucopolysaccharides produced by cancer cells leading to the formation of microcalcifications (25). The metastatic lymph nodes may have the characteristic features of the primary lesion, so it is inferred that the microcalcification in lymph nodes may originate from PTC. It has been demonstrated that hyperechoic changes in lymph nodes are caused by aggregates of thyroglobulin. Metastatic cancer cells are able to synthesize thyroglobulin but they lack an intact follicular structure. The synthesized thyroglobulin cannot be utilized and clumps form, presenting as hyperechoic change (26).

Zhang et al. (27) found that the maximum diameter of the nodule was the most important variable based on the study of 337 patients who had undergone the surgery on lymph node metastasis of PTC. Xue et al. (28) showed in their study that the risk of CLNM is greater when tumor size ≥1 cm. The maximum diameter of the nodule is related to tumor growth and proliferation, indicating tumors are more aggressive and prone to lymph node metastasis. The larger the nodule, the closer it is to the thyroid capsule. Larger nodules are prone to lymph node metastasis.

The hypoechoic halo may result from a progressing desmoplastic reaction or fibrosis of the capsule after the presence of microscopic capsular invasion (29). Some previous studies showed that the halo around thyroid nodules was caused by the rapid growth of the lesions and the extrusion of surrounding tissues (30). Our study showed that there was significant differences between LNM and non-LNM groups in terms of hypoechoic halo around the nodule. Hypoechoic halo may be associated with LNM, but it’s not a risk factor highly associated with LNM in PTC. This result corresponds with the relevant literature (3133). We speculate that the appearance of hypoechoic halo indicates the malignant biological behavior of PTC cancer cells, which means that the tumor has a relatively strong ability to grow and infiltrate, which may be an important reason for its association with lymph node metastasis.

PI means maximum dose of microbubbles filling the ROI within a certain period of time, and AUC means total dose of microbubbles filling the ROI within a certain period of time, compared with those of adjacent tissue (34). Fast growing tumors are easy to ischemic necrosis and lack of nourishing blood vessels. The blood vessels of the high malignant tumors are messy and tortuous, which are easily blocked by tiny emboli and lead to stenosis (35). Therefore, the low values of PI and AUC have reference significance in predicting LNM of PTC.

The prediction model was developed through univariate and multivariate logistic regression analyses. And the model was evaluated and validated by ROC, calibration curve, cross-validation, DCA, and clinical impact curve. The results showed the prediction model was good and had favorable clinical practicability. We applied and registered our nomogram prediction model online. This model can be used to predict LNM through mobile phone, computer and other network login, which is convenient for clinical application.

Despite the predictive ability of our nomogram, there are still several limitations. Firstly, the number of cases for qualitative and quantitative CEUS analysis was insufficient. Therefore, the sample size can be expanded in the future. Secondly, indicators of gene mutations associated with PTC, such as BRAFV600E mutation (36) should be also involved in future analyses. The model can be further improved by incorporating ultrasound elastography, CT and other imaging techniques, as well as computer-aided diagnosis, to enhance its diagnostic ability.

5 Conclusion

Our study identified that male sex, younger age, hypoechoic halo around the nodule, thyroid capsule invasion, lymph node microcalcification, lymph node hyperechoic change, PI and AUC are associated with an increased risk of LNM in patients with PTC. Nomogram based on sex, age, thyroid capsule invasion, lymph node microcalcification has the potential for preoperative LNM risk assessment. The nomogram serves as a valuable clinical tool for aiding in clinical decision-making. It conveniently provides for individual preoperative prediction and might potentially improve the survival outcome in PTC patients with LNM.

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 humans were approved by Ethics Committee of Xiangya Hospital (No. 202206140). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

JuL: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Resources, Software, Visualization, Writing – original draft, Writing – review & editing. JinL: Conceptualization, Data curation, Investigation, Project administration, Resources, Supervision, Validation, Writing – review & editing. YC: Data curation, Formal Analysis, Resources, Writing – original draft. JiL: Data curation, Resources, Software, Writing – original draft. XH: Data curation, Resources, Writing – original draft. HZ: Data curation, Formal Analysis, Supervision, Writing – review & editing. BZ: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing, Writing – original draft.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study has received funding by Natural Science Foundation of Hunan Province (2022JJ30982).

Acknowledgments

I would like to thank my supervisor, BZ, whose expertise was invaluable in creating the study design and experimentation.

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.

Abbreviations

PTC, papillary thyroid carcinoma; LNM, lymph node metastasis; ROC, receiver operating characteristic; DCA, decision curve analysis; PI, peak intensity of contrast; AUC, area under the curve; CEUS, contrast-enhanced ultrasound; 2D-US, two-dimensional ultrasound; FT4, free thyroid hormone; FT3, free triiodothyronine; TSH, hypersensitive thyroid stimulating hormone; TG, thyroglobulin; TGA, antithyroglobulin antibody; TPOAB, thyroid peroxidase antibody; PTH, parathyroid hormone; CEA, carcinoembryonic antigen; BI, Base Intensity; AT, Arrival Time; TTP, Time to Peak; PI, Peak Intensity; AS, Ascending Slope; ROI, region of interest; GOF, highest goodness of fit.

References

1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Baloch Z, LiVolsi VA. Fifty years of thyroid pathology: concepts and developments. Hum Pathol (2020) 95:46–54. doi: 10.1016/j.humpath.2019.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Wang W, Yang Z, Ouyang Q. A nomogram to predict skip metastasis in papillary thyroid cancer. World J Surg Oncol (2020) 18(1):167. doi: 10.1186/s12957-020-01948-y

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Patel KN, Yip L, Carty SE. Response to the comment on “The american association of endocrine surgeons guidelines for the definitive surgical management of thyroid disease in adults”. Ann Surg (2021) 274(6):e746–7. doi: 10.1097/SLA.0000000000004274

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Sun J, Jiang Q, Wang X, Liu W, Wang X. Nomogram for preoperative estimation of cervical lymph node metastasis risk in papillary thyroid microcarcinoma. Front Endocrinol (Lausanne) (2021) 12:613974. doi: 10.3389/fendo.2021.613974

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Ryu YJ, Kwon SY, Lim SY, Na YM, Park MH. Predictive factors for skip lymph node metastasis and their implication on recurrence in papillary thyroid carcinoma. Biomedicines (2022) 10(1):179. doi: 10.3390/biomedicines10010179

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Yu P, Wu X, Li J, Mao N, Zhang H, Zheng G, et al. Extrathyroidal extension prediction of papillary thyroid cancer with computed tomography based radiomics nomogram: A multicenter study. Front Endocrinol (Lausanne) (2022) 13:874396. doi: 10.3389/fendo.2022.874396

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Tsou P, Wu CJ. Mapping driver mutations to histopathological subtypes in papillary thyroid carcinoma: applying a deep convolutional neural network. J Clin Med (2019) 8(10):1675. doi: 10.3390/jcm8101675

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Bojunga J. Ultrasound of thyroid nodules. Sonografie von schilddrüsenknoten. Ultraschall Med (2018) 39(5):488–511. doi: 10.1055/a-0659-2350

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Chen L, Chen L, Liu J, Wang B, Zhang H. Value of qualitative and quantitative contrast-enhanced ultrasound analysis in preoperative diagnosis of cervical lymph node metastasis from papillary thyroid carcinoma. J Ultrasound Med (2020) 39(1):73–81. doi: 10.1002/jum.15074

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Guo Q, Sun C, Chang Q, Wang Y, Chen Y, Wang Q, et al. Contrast-enhanced ultrasound-based nomogram for predicting Malignant involvements among sonographically indeterminate/suspicious cervical lymph nodes in patients with differentiated thyroid carcinoma. Ultrasound Med Biol (2022) 48(8):1579–89. doi: 10.1016/j.ultrasmedbio.2022.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Wei Y, Yu MA, Niu Y, Hao Y, Di JX, Zhao ZL, et al. Combination of lymphatic and intravenous contrast-enhanced ultrasound for evaluation of cervical lymph node metastasis from papillary thyroid carcinoma: A preliminary study. Ultrasound Med Biol (2021) 47(2):252–60. doi: 10.1016/j.ultrasmedbio.2020.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Canu GL, Medas F, Cappellacci F, Giordano ABF, Gurrado A, Gambardella C, et al. Risk of complications in patients undergoing completion thyroidectomy after hemithyroidectomy for thyroid nodule with indeterminate cytology: an italian multicentre retrospective study. Cancers (Basel) (2022) 14(10):2472. doi: 10.3390/cancers14102472

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Tolone S, Roberto R, del Genio G, Brusciano L, Parmeggiani D, Amoroso V, et al. The impact of age and oral calcium and vitamin D supplements on postoperative hypocalcemia after total thyroidectomy. A prospective study. BMC Surg (2013) 13(2):S11. doi: 10.1186/1471-2482-13-S2-S11

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Gambardella C, Offi C, Romano RM, De Palma M, Ruggiero R, Candela G, et al. Transcutaneous laryngeal ultrasonography: a reliable, non-invasive and inexpensive preoperative method in the evaluation of vocal cords motility-a prospective multicentric analysis on a large series and a literature review. Updates Surg (2020) 72(3):885–92. doi: 10.1007/s13304-020-00728-3

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Filetti S, Durante C, Hartl D, Leboulleux S, Locati LD, Newbold K, et al. Thyroid cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up†. Ann Oncol (2019) 30(12):1856–83. doi: 10.1093/annonc/mdz400

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Chen GG, Vlantis AC, Zeng Q, van Hasselt CA. Regulation of cell growth by estrogen signaling and potential targets in thyroid cancer. Curr Cancer Drug Targets (2008) 8(5):367–77. doi: 10.2174/156800908785133150

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Lorenz K, Schneider R, Elwerr M. Thyroid carcinoma: do we need to treat men and women differently? Visc Med (2020) 36(1):10–4. doi: 10.1159/000505496

PubMed Abstract | CrossRef Full Text | Google Scholar

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

20. Liu Y, Wang Y, Zhao K, Li D, Chen Z, Jiang R, et al. Lymph node metastasis in young and middle-aged papillary thyroid carcinoma patients: a SEER-based cohort study. BMC Cancer (2020) 20(1):181. doi: 10.1186/s12885-020-6675-0

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Lang BH, Ng CP, Au KB, Wong KP, Wong KK, Wan KY. Does preoperative neutrophil lymphocyte ratio predict risk of recurrence and occult central nodal metastasis in papillary thyroid carcinoma? World J Surg (2014) 38(10):2605–12. doi: 10.1007/s00268-014-2630-z

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Seong CY, Chai YJ, Lee SM, Kim SJ, Choi JY, Lee KE, et al. Significance of distance between tumor and thyroid capsule as an indicator for central lymph node metastasis in clinically node negative papillary thyroid carcinoma patients. PloS One (2018) 13(7):e0200166. doi: 10.1371/journal.pone.0200166

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Lyu YS, Pyo JS, Cho WJ, Kim SY, Kim JH. Clinicopathological significance of papillary thyroid carcinoma located in the isthmus: A meta-analysis. World J Surg (2021) 45(9):2759–68. doi: 10.1007/s00268-021-06178-1

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Yin L, Zhang W, Bai W, He W. Relationship between morphologic characteristics of ultrasonic calcification in thyroid nodules and thyroid carcinoma. Ultrasound Med Biol (2020) 46(1):20–5. doi: 10.1016/j.ultrasmedbio.2019.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Ferreira LB, Gimba E, Vinagre J, Sobrinho-Simões M, Soares P. Molecular aspects of thyroid calcification. Int J Mol Sci (2020) 21(20):7718. doi: 10.3390/ijms21207718

PubMed Abstract | CrossRef Full Text | Google Scholar

26. d'Herbomez M, Lion G, Béron A, Wémeau JL, DoCao C. Advances in thyroglobulin assays and their impact on the management of differentiated thyroid cancers. Ann Biol Clin (Paris) (2016) 74(1):21–7. doi: 10.1684/abc.2015.1106

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Zhang C, Cheng L, Zhu W, Zhuang J, Zhao T, Li X, et al. Construction of a diagnostic model for lymph node metastasis of the papillary thyroid carcinoma using preoperative ultrasound features and imaging omics. J Healthc Eng (2022) 2022:1872412. doi: 10.1155/2022/1872412

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Xue T, Liu C, Liu JJ, Hao YH, Shi YP, Zhang XX, et al. Analysis of the relevance of the ultrasonographic features of papillary thyroid carcinoma and cervical lymph node metastasis on conventional and contrast-enhanced ultrasonography. Front Oncol (2021) 11:794399. doi: 10.3389/fonc.2021.794399

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Zhang F, Chen W. Sonographic features of follicular variant of papillary thyroid carcinoma (FV-PTC) and diagnostic performance of the 2017 ACR TI-RADS in FV-PTC. Endocrine (2020) 67(2):379–86. doi: 10.1007/s12020-019-02184-5

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Ng SC, Kuo SF, Hua CC, Huang BY, Chiang KC, Chu YY, et al. Differentiation of the follicular variant of papillary thyroid carcinoma from classic papillary thyroid carcinoma: an ultrasound analysis and complement to fine-needle aspiration cytology. J Ultrasound Med (2018) 37(3):667–74. doi: 10.1002/jum.14377

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Sun G, Xie YD, Xu CL, Yang B. Clinical value of a nomogram based on ultrasonic and clinical features for predicting central and lateral cervical lymph node metastases of thyroid papillary carcinoma. Chin J Med Ultrasound (Electronic Edition) (2023) 20:734–42. doi: 10.3877/cma.j.issn.1672-6448.2023.07.012

CrossRef Full Text | Google Scholar

32. Zhou J, Zhou SC, Li JW, Wang Y, Chen YL, Wang F, et al. Risk factors of predicting lateral neck lymph node metastasis following solitary papillary thyroid carcinoma. Chin J Ultrasonography (2019) 28:971–5. doi: 10.3760/cma.j.issn.1004-4477.2019.11.009

CrossRef Full Text | Google Scholar

33. Lan Y, Song Q, Jin Z, Zhang Y, Lin L, Luo YK. Correlation of routine ultrasound features and BRAFV600E gene to the cervical lymph node metastasis of thyroid papillary carcinoma. Med J Chin People's Liberation Army (2019) 44:747–52. doi: 10.11855/j.issn.0577-7402.2019.09.06

CrossRef Full Text | Google Scholar

34. Fang F, Gong Y, Liao L, Ye F, Zuo Z, Li X, et al. Value of contrast-enhanced ultrasound for evaluation of cervical lymph node metastasis in papillary thyroid carcinoma. Front Endocrinol (Lausanne) (2022) 13:812475. doi: 10.3389/fendo.2022.812475

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Jiang W, Wei HY, Zhang HY, Zhuo QL. Value of contrast-enhanced ultrasound combined with elastography in evaluating cervical lymph node metastasis in papillary thyroid carcinoma. World J Clin cases (2019) 7(1):49–57. doi: 10.12998/wjcc.v7.i1.49

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Siraj AK, Parvathareddy SK, Pratheeshkumar P, Divya SP, Al-Sobhi SS, Al-Dayel F, et al. PD-L1 is an independent prognostic marker in middle eastern PTC and its expression is upregulated by BRAFV600E mutation. Cancers (2021) 13(3):555. doi: 10.3390/cancers13030555

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: papillary thyroid carcinoma, lymph node metastasis, ultrasound, risk factor, nomogram

Citation: Lu J, Liao J, Chen Y, Li J, Huang X, Zhang H and Zhang B (2023) Risk factor analysis and prediction model for papillary thyroid carcinoma with lymph node metastasis. Front. Endocrinol. 14:1287593. doi: 10.3389/fendo.2023.1287593

Received: 02 September 2023; Accepted: 16 October 2023;
Published: 30 October 2023.

Edited by:

Mehmet Uludag, Şişli Hamidiye Etfal Education and Research Hospital, Türkiye

Reviewed by:

Gökhan Içöz, Ege University, Türkiye
Ludovico Docimo, University of Campania Luigi Vanvitelli, Italy
Peisong Wang, The 1st Hospital of Jilin University, China

Copyright © 2023 Lu, Liao, Chen, Li, Huang, Zhang and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Bo Zhang, zhangbo8095@csu.edu.cn

Disclaimer: 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.