Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Endocrinol., 28 January 2026

Sec. Thyroid Endocrinology

Volume 17 - 2026 | https://doi.org/10.3389/fendo.2026.1709522

Standardizing cervical lymph node evaluation in papillary thyroid carcinoma: diagnostic accuracy of Node-RADS on ultrasound

Hongjun ZhangHongjun Zhang1Danni HeDanni He1Muyi MaoMuyi Mao1Renjie LiRenjie Li1Tinghui YinTinghui Yin2Bo-wen ZhengBo-wen Zheng2Jiayi ZhengJiayi Zheng3Tao WuTao Wu2Jie Ren*Jie Ren2*Zuofeng Xu*Zuofeng Xu1*
  • 1Seventh Affiliated Hospital, Sun Yat-sen University, Shenzhen, China
  • 2Third Affiliated Hospital of Sun Yat-Sen University, Guangzhou, China
  • 3Guangdong Women and Children Hospital, Guangzhou, China

Objectives: To assess the diagnostic performance of the Node Reporting and Data System (Node-RADS) in ultrasonographic evaluation of cervical lymph node (LN) metastasis in papillary thyroid carcinoma (PTC) and compare it with conventional suspicious sonographic features.

Materials and methods: This retrospective study included 227 PTC patients (269 LNs) who underwent preoperative cervical LN ultrasound and biopsy (September 2022–December 2023). Two radiologists independently assigned Node-RADS scores (1–5) based on LN size and morphologic features (shape, border, echogenicity). Diagnostic performance was evaluated using ROC analysis, with histopathology as the reference standard. Node-RADS was compared to five American Thyroid Association (ATA)-defined suspicious features (microcalcifications, cystic changes, etc.). Multivariable analysis assessed an integrated model combining Node-RADS with peripheral vascularity.

Results: Metastatic LNs (63.9%, 172/269) exhibited higher Node-RADS scores (malignancy rates: 35.7% for score 1 vs. 88.2% for score 5). Node-RADS achieved an AUC of 0.720 (95% CI: 0.662–0.773), surpassing individual features (best AUC: 0.703 for peripheral vascularity). Inter-reader agreement was good (κ = 0.67). A combined Node-RADS and peripheral vascularity model significantly improved diagnostic performance (AUC: 0.825, p < 0.01).

Conclusions: Node-RADS offers standardized, reproducible risk stratification for cervical LN metastasis in PTC, with moderate diagnostic accuracy. Integration with peripheral vascularity enhances performance, supporting its potential clinical utility in preoperative assessment.

Introduction

Papillary thyroid carcinoma (PTC) comprises approximately 80–90% of thyroid cancers (1). Its global incidence has been rising steadily, largely due to the increased use of high-resolution imaging and fine-needle aspiration biopsy (FNAB) (2, 3). Although PTC generally exhibits indolent behavior, cervical lymph node (LN) metastasis occurs in up to 80% of patients and plays a critical role in staging, prognosis, and recurrence of PTC, however, accurate preoperative assessment remains challenging (4). Ultrasonography (US) is the primary imaging modality for evaluating thyroid nodules and cervical LNs, favored for its non-invasiveness, accessibility, and high spatial resolution (5, 6). However, its diagnostic accuracy in identifying metastatic lymph nodes in PTC patients remains suboptimal (7, 8). Current US-based evaluations depend critically on operator experience and interpretation of heterogeneous sonographic features, such as shape, margin, echogenicity, vascularity, and calcifications (1, 9). This variability underscores the need for a more objective and reproducible scoring system to improve the diagnostic performance of US in detecting cervical LN metastasis.

Structured imaging reporting systems have demonstrated clinical value in other malignancies. Systems such as the Breast Imaging Reporting and Data System (BI-RADS) and Prostate Imaging Reporting and Data System (PI-RADS) have improved diagnostic consistency by standardizing the interpretation of imaging findings (1013). Despite the high prevalence of LN metastasis in PTC, a comparable standardized system for cervical LN evaluation on US remains unavailable for wide implementation.

To enhance objectivity and standardization in LN assessment, the Node Reporting and Data System (Node-RADS) 1.0 was introduced across CT and MR imaging modalities (14). Preliminary studies suggested its potential utility for LN assessment in various malignancies (1519). However, its performance in US-based assessment of cervical LNs in PTC patients has not been systematically investigated.

This study aimed to investigate the diagnostic performance of Node-RADS for cervical LNs in PTC patients. We further compared its discriminative ability with individual suspicious sonographic features and explored an integrated diagnostic model to improve clinical applicability.

Materials and methods

Ethical approval and study design

This retrospective cohort study was conducted from the third affiliated hospital of Sun Yat-sen University between September 2022 and December 2023. It was approved by the institutional review board of our hospital, and informed consent was waived.

Patient selection criteria

Consecutive patients were included based on the following criteria (1): histologically confirmed PTC (2); presence of cervical LNs suspected malignant on US (3); underwent both diagnostic US and subsequent US-guided FNA or excisional biopsy between September 2022 and December 2023 (4); availability of complete B-mode US and color Doppler flow imaging (CDFI) data.

Exclusion criteria comprised (1): inconclusive pathological results (2); non-diagnostic US image quality (3); prior radioactive iodine (¹³¹I) therapy (4); concurrent malignancies (5); missing clinical or imaging data.

The patient selection process is outlined in Figure 1. Histopathological analysis of US-guided and/or excision biopsies served as the reference standard, with benign LNs undergoing 6-month US follow-up (20, 21).

Figure 1
Flowchart showing the study selection process for cervical lymph nodes (LNs) in PTC patients from September 2022 to December 2023. Out of 337 LNs, 68 were excluded due to previous radioactive iodine treatment, poor ultrasound quality, or no histological result. A total of 269 LNs from 227 patients were included. These were divided into 97 benign LNs, further categorized as reactive hyperplastic lymphadenopathy (95), tuberculous (1), and necrotizing lymphadenitis (1), and 172 malignant LNs, all metastatic.

Figure 1. Flowchart for patient selection. LN, lymph node; PTC, papillary thyroid carcinoma; US, ultrasound.

Sample size calculation

A formal sample size calculation was performed before the analysis. For the development of the multivariable logistic regression model, the rule of thumb of 10 events per variable (EPV) was applied (2224). With 6 predictor variables in the model and an anticipated event (malignancy) rate of approximately 40% based on prior literature in similar populations (1), a minimum of 150 samples (60 events/0.4) was required. Our final cohort of 269 lymph nodes, therefore, provides adequate statistical power for robust model development and validation.

Ultrasound protocol

All US examinations and biopsies were performed using high-resolution linear transducers (7–15 MHz) under standardized protocols, including 5 standardized ultrasonography instruments: ALOKA Prosound Alpha 6 (Hitachi), ALOCA ARIETTA 60 (Hitachi), ALOCA ARIETTA 70 (Hitachi), GE Logiq E9 (GE), and Acuson Sequoia (Siemens). Three fellowship-trained radiologists (10–15 years experience in thyroid imaging) performed all scans and documented LN features, including anatomical location, short-axis diameter, short-to-long axis ratio, cortical thickness, internal echotexture (echogenicity, calcifications, cystic changes), and vascularity pattern (hilar, peripheral, or mixed). Baseline demographics (age, sex) and pathological results were recorded. In accordance with standard clinical practice, all lymph nodes subjected to biopsy were initially marked under ultrasound guidance. During subsequent surgical excision, the nodes were identified and matched based on their anatomical location and correlative skin markings.

Node-RADS assessment

Before the study, two radiologists (with 12 and 4 years of experience, respectively) attended a dedicated training session on the Node-RADS criteria, which included a review of the scoring guidelines and a set of 20 practice cases (not included in the final analysis) to ensure consistent understanding and application of the criteria. In clinical practice, assigning a Node-RADS score requires approximately 1–2 minutes per case. To minimize potential bias, the readers were blinded to all clinical information, subsequent pathological results, and the individual assessments of their fellow readers throughout the Node-RADS scoring process. The Node-RADS score is determined by combining scores from LN size and configuration (including shape, border, and cortical echogenicity) (14). The size category was stratified as normal, enlarged, or bulky, while “configuration” was evaluated based on border regularity, shape, and cortical echogenicity. Each subcategory contributed a partial score, and the cumulative sum determined the final Node-RADS score, ranging from 1 to 5, indicating progressively increasing suspicion of malignancy.

LNs in cervical regions were categorized as “enlarged” if the short-axis diameter was ≥ 10 mm. Vascularity was assessed on CDFI and graded as absent (G0), minimal (G1; one or two discrete color signals < 1 mm), moderate (G2), or marked (G3; ≥4 discrete vessels) (25). Discrepancies in scoring were resolved by choosing the higher of the two scores.

Statistical analyses

LNs were classified as benign or metastatic based on pathological confirmation. Continuous variables were expressed as mean ± standard deviation or median [interquartile range], and categorical variables as counts with percentages. Group comparisons used the independent t-test or Mann–Whitney U test for continuous variables, and the chi-square test for categorical variables.

Receiver operating characteristic (ROC) curves were constructed to evaluate the diagnostic performance of Node-RADS and compared to five American Thyroid Association (ATA)-defined suspicious sonographic features: focal hyperechogenicity, microcalcification, cystic degeneration, round shape, and peripheral vascularity (1). Sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were calculated at the optimal Node-RADS cutoff, determined by the maximum Youden index.

Potential predictors were initially screened using univariable logistic regression. To avoid multicollinearity, the individual components of the Node-RADS score (e.g., size, shape, echogenicity) were not re-evaluated as separate variables. Candidate variables with a p-value of less than 0.10 from the univariable analysis were included in the subsequent multivariable model. A final multivariable logistic regression model was then constructed with the independent predictors. The discriminative performance of this final model was evaluated using ROC curve analysis. To mitigate the risk of overfitting and to substantiate the robustness of the combined diagnostic model, an internal validation was performed using the 10-fold cross-validation method. The average area under the receiver operating characteristic curve (AUC) derived from these validation cycles is reported. Inter-reader agreement for Node-RADS scoring was assessed using Cohen’s kappa statistics (κ), interpreted as follows: excellent (κ > 0.80), good (κ = 0.61–0.80), moderate (κ = 0.41–0.60), fair = 0.21–0.40), and poor (κ < 0.20). In an effort to investigate a potential learning curve, a post-hoc analysis was performed by stratifying the reading process into first and second halves. All statistical analyses were performed using SPSS version 25.0, MedCalc version 20.0, and R software version 44.5.2, with a two-tailed p-value < 0.05 considered statistically significant.

Result

Patient characteristics

A total of 227 patients (mean age of 40 ± 14 years; age range, 10–82 years; 71.8% female) with 269 cervical LNs were included in the final analysis. Excluded cases comprised 12 with suboptimal image quality, 6 with inconclusive histology, and 50 with prior 131I therapy (Figure 1).

Pathological confirmation was obtained via FNA alone (n = 85, 31.6%), core-needle biopsy (CNB) alone (n = 37, 13.8%), FNA with excision (n = 115, 42.8%), and CNB with excision (n = 32, 11.9%) (Table 1). The high concordance rate of 92.2% (107/116) between FNA and subsequent surgical excision, which is consistent with the previous literature (26, 27), supports FNA as a robust diagnostic standard in our study. Final pathology revealed 172 metastatic LNs (63.9%) and 97 benign LNs (36.1%), including reactive hyperplasia (n = 95), necrotizing lymphadenitis (n = 1), and tuberculous lymphadenitis (n = 1).

Table 1
www.frontiersin.org

Table 1. Baseline demographic and clinical characteristics of the study.

Node-RADS score distribution and interreader agreement

Among the 269 LNs, the distribution of Node-RADS scores was as follows: score 1 (n = 28), score 2 (n = 55), score 3 (n = 46), score 4 (n = 64), and score 5 (n = 76). The corresponding malignancy rates were 35.7%, 50.9%, 45.7%, 71.9%, and 88.2%, respectively. Representative ultrasound images of each Node-RADS category are presented in Figure 2.

Figure 2
Ultrasound images displaying lymph nodes classified by Node-RADS from one to five. The top row presents 2D ultrasound images, labeled A, C, E, G, and I, showing varying lymph node characteristics. The bottom row shows Color Doppler Flow Imaging (CDFI) images, labeled B, D, F, H, and J, highlighting blood flow differences within each node, using red and blue color mapping for analysis.

Figure 2. Representative cervical LN of different Node-RADS scores. (A, B) Node-RADS score 1: short axis = 4 mm (0 points), homogeneous texture (0 points), smooth borders (0 points), and kidney-bean-like shape with fatty hilum (0 points); (C, D) Node-RADS score 2: short axis = 5 mm (0 points), heterogeneous texture (1 point), smooth border (0 points), and kidney-bean-like shape without fatty hilum (0 points); (E, F) Node-RADS score 3: short axis = 10 mm (1 point), heterogeneous texture (1 point), smooth border (0 points), and oval shape without fatty hilum (0 points); (G, H) Node-RADS score 4: short axis = 9 mm (0 points), heterogeneous (1 point), irregular border (1 point), and spherical without fatty hilum (1 point); (I, J) Node-RADS score 5: short axis = 16 mm (1 point), calcification (3 points), irregular border (1 point), and spherical without fatty hilum (1 point). LN, lymph node; CDFI, Color Doppler flow imaging; Node-RADS, Node Reporting and Data System.

Inter-reader diagnostic consistency for Node-RADS scoring yielded a Cohen’s κ of 0.67 (95% confidence interval [CI]: 0.60–0.73), indicating good concordance between the two radiologists.

The inter-reader agreement was found to be κ = 0.66 in the first half and κ = 0.68 in the second half. While this indicates a trend towards improved concordance, the difference was not statistically significant (P = 0.77). This finding may suggest a rapid familiarization period with the Node-RADS system.

Diagnostic performance of Node-RADS and five suspected US features

Node-RADS demonstrated superior diagnostic performance over individual sonographic features, achieving an area under the curve (AUC) of 0.720 (95% CI: 0.662–0.773, p < 0.001). The optimal cutoff score was > 3, yielding a sensitivity of 65.7% (95% CI: 58.1% - 72.8%) and a specificity of 72.2% (95% CI: 62.1% - 80.8%). The positive predictive value (PPV) and negative predictive value (NPV) were 80.7% (95% CI: 74.9% - 85.4%) and 54.3% (95% CI: 48.3% - 60.2%), respectively (Table 2).

Table 2
www.frontiersin.org

Table 2. Diagnostic performance for detecting cervical lymph node metastasis.

Among individual ultrasound features, peripheral vascularization had the highest AUC (0.703, 95% CI: 0.644 - 0.757, p < 0.001), followed by microcalcification (AUC = 0.663, 95% CI: 0.603 - 0.719, p < 0.001), which demonstrated high specificity (86.6%, 95% CI: 78.2% - 92.7%) but low sensitivity (45.9%, 95% CI: 38.3% - 53.7%). Focal hyperechogenicity and round shape yielded lower diagnostic accuracy (AUC = 0.578 and 0.585, respectively), with the cystic degeneration showing the lowest overall performance (AUC = 0.542, 95% CI: 0.481 - 0.603, p = 0.03) (Table 2, Figure 3).

Figure 3
ROC curves comparing diagnostic models with sensitivity on the y-axis and one minus specificity on the x-axis. Six curves are shown: Node-RADS (blue, AUC = 0.720), Focal hyperechogenicity (green, AUC = 0.578), Microcalcification (orange, AUC = 0.663), Cystic degeneration (dark green, AUC = 0.542), Round shape (pink, AUC = 0.585), Peripheral vascularisation (brown, AUC = 0.703).

Figure 3. Diagnosis performance of Node-RADS and five suspected US features. Node-RADS, Node Reporting and Data System; US, ultrasound; AUC, Area under the curve.

Diagnostic performance of combining Node-RADS and other features

In multivariate logistic regression, both Node-RADS and peripheral vascularization were identified as independent predictors of cervical LN metastasis. The representative vascular patterns of cervical LNs were shown in Figure 4. The combined model incorporating these two variables significantly improved diagnostic performance (AUC = 0.825, 95% CI: 0.776–0.867, p < 0.0001), compared to Node-RADS alone (AUC = 0.720, 95% CI: 0.662–0.773, p < 0.0001) (p < 0.01) (Table 3, Figure 5). The result of 10-fold cross-validation procedure yielded a mean AUC of 0.814 (standard deviation: ± 0.039), thereby confirming the model’s robust discriminative performance.

Figure 4
Three ultrasound images labeled A, B, and C, each showing color Doppler scans of blood flow in neck tissues. Image A depicts a central vessel with red and blue color streams. Image B shows a similar vessel with concentrated red and blue areas, all indicating blood movement within the tissue. Image C displays a vessel with mixed color flow patterns. Each image has a scale marker of one centimeter.

Figure 4. Representative vascular patterns of cervical LNs, including hilar (A), peripheral (B), and mixed (C). LN, lymph node.

Table 3
www.frontiersin.org

Table 3. Results of univariable and multivariable logistic regression.

Figure 5
ROC curve comparing two models. The solid blue line represents Node-RADS with an ROC of 0.720. The dashed green line represents Node-RADS plus Peripheral Vascularization with an ROC of 0.825. The y-axis indicates sensitivity, and the x-axis indicates 1-specificity.

Figure 5. Diagnostic performance of Node-RADS and Node-RADS + peripheral vascularization. Node-RADS, Node Reporting and Data System; AUC, Area under the curve.

Discussion

While Node-RADS was originally developed for CT/MRI, its systematic adaptation and validation for US in the evaluation of cervical lymph nodes in PTC remains largely unexplored. As a primary study, our study not only confirms the feasibility of applying the Node-RADS framework in US, where its core criteria are fully applicable, but also moves beyond mere translation. We developed and validated a novel combined diagnostic model that integrates Node-RADS with the sonographic feature of peripheral vascularity, which significantly improved diagnostic performance (AUC increased from 0.720 to 0.825). Notably, we found a strong correlation between increasing Node-RADS scores and malignancy rates, with scores of 4 and 5 associated with a 71.9% and 88.2% malignancy risk, respectively. This work provides a critical foundation for and is part of an ongoing research program aimed at establishing a refined, US-specific scoring system for lymph node assessment in PTC.

The diagnostic accuracy of Node-RADS (AUC = 0.720, 95% CI: 0.662–0.773, p < 0.01) exceeded that of commonly used US features such as peripheral vascularity (AUC = 0.703, 95% CI: 0.644 - 0.757, p < 0.001), microcalcification (AUC = 0.663, 95% CI: 0.603 - 0.719, p < 0.001), and round shape (AUC = 0.585, 95% CI: 0.603 - 0.719, p = 0.01). This supported the premise that a composite scoring system integrating multiple morphologic criteria could outperform single-feature evaluations. Among the individual features, peripheral vascularity showed the highest AUC, consistent with previous studies reporting neovascularization as a hallmark of nodal metastasis (26, 28). However, features such as degeneration, although highly specific (92.8%, 95% CI: 85.7% - 97.0%), demonstrated low sensitivity (15.7%, 95% CI: 10.6% - 22.0%), limiting their standalone diagnostic value.

According to previous studies, the Node-RADS system has shown promising results in evaluating lymph node metastasis across various malignancies, including prostate, breast, lung, colon, and gastric cancers, using CT and MRI (2932). However, its performance in distinguishing benign from malignant cervical lymph nodes by US in PTC patients was relatively modest, with an AUC of 0.720 (95% CI: 0.662–0.773, p < 0.01). A similar result was reported by Yu et al. by evaluating the diagnosis performance of Node-RADS using CT images to predict metastatic LNs in patients with PTC, with an AUC of 0.602 (5). This might be attributed to the smaller size of metastatic lymph nodes in PTC, which limits the ability to observe internal structural details. Additionally, the retrospective nature of this study might have introduced inaccuracies in assessing internal lymph node characteristics.

Importantly, our multivariate analysis revealed that both Node-RADS and peripheral vascularity were independent predictors of nodal metastasis. The combined diagnostic model (Node-RADS + peripheral vascularity) significantly improved diagnostic performance (AUC = 0.825, 95% CI: 0.776–0.867, p < 0.001), underscoring the complementary role of vascular features in enhancing malignancy prediction. The improvement in AUC is not only statistically significant but also suggests a clinically meaningful enhancement in risk stratification. This could potentially refine clinical pathways by more accurately identifying patients who would benefit directly from surgical intervention versus those for whom a confirmatory FNA remains the most prudent step. This performance was comparable to the risk stratification model for metastatic LN of PTC using US images developed by Ni et al. (AUC = 0.827) (33). These findings aligned with prior work suggesting that vascular flow patterns, particularly peripheral or chaotic vascularization, add incremental value to grayscale-based assessments (34, 35). The biological plausibility of this observation is supported by histopathologic studies. Specifically, peripheral vascularity on Doppler US has been shown to correlate with significantly increased microvessel density in the subcapsular sinus of metastatic lymph nodes (36). The underlying mechanism can be explained by the process of nodal metastasis: in the early stages, tumor cells arrive via lymphatic vessels and initially reside in the subcapsular sinus. These cells secrete pro-angiogenic factors such as VEGF, which is significantly overexpressed in metastatic lymph nodes of PTC (37), thereby inducing neoangiogenesis and resulting in the observed peripheral vascularization. As the tumor progresses, cells infiltrate deeper into the cortex and medulla, disrupting the normal nodal architecture and leading to more diffuse vascularization (3841). This established pathologic sequence underscores the importance of peripheral vascularity as a critical indicator for assessing lymph node metastasis in PTC patients.

Our findings demonstrated that Node-RADS was a reproducible tool, with a substantial interreader agreement (κ = 0.67, 95% CI: 0.60–0.73), which suggested Node-RADS could facilitate more consistent interpretation among radiologists with varying experience levels. Previous studies have highlighted variability in LN characterization using traditional ultrasound criteria (7, 9); thus, standardized reporting systems such as Node-RADS might bridge the gap between subjectivity and reproducibility in clinical practice.

Our study had several limitations. First, this study is a single-center retrospective analysis, and the inclusion of only sonographically suspicious lymph nodes may introduce spectrum bias, potentially overestimating diagnostic performance. Additionally, the duration of clinical follow-up for benign lymph nodes was relatively limited. Also, the inter-scanner variability could affect the echogenicity and vascularity assessment. Therefore, the generalizability of our findings to unselected populations requires further validation through our ongoing prospective multi-center studies utilizing the same equipment and a standardized image archiving protocol with long-term follow-up. Second, despite rigorous reader training and blinding, subjective interpretation of ultrasound features might still influence Node-RADS scoring. Third, pathologic confirmation was not uniform, with some cases relying solely on cytology, which might carry limitations in diagnostic certainty. However, the high concordance rate of 92.2% (107/116) between FNA and subsequent surgical excision, in accordance with the literature (26, 27), supports the reliability of FNA as a precise technique, and the limitation can be overlooked.

In conclusion, our study is an exploratory, preliminary validation of Node-RADS as an effective and reproducible scoring system for US evaluation of cervical LNs in PTC patients. Its diagnostic performance surpassed that of individual ultrasound features, and combination with peripheral vascularity further enhanced accuracy. Future studies with larger and more diverse populations are warranted to refine the scoring criteria and validate its application across different clinical settings and disease types.

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 humans were approved by the Third Affiliated Hospital of Sun Yat-sen University. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent from the participants or the participants’ legal guardians/next of kin because of the retrospective nature of the study.

Author contributions

HZ: Methodology, Writing – original draft. DH: Methodology, Writing – original draft. MM: Methodology, Writing – original draft. RL: Investigation, Software, Writing – review & editing. TY: Investigation, Project administration, Writing – review & editing. B-wZ: Project administration, Writing – review & editing. JZ: Data curation, Writing – review & editing. TW: Formal analysis, Writing – review & editing. JR: Conceptualization, Resources, Supervision, Validation, Visualization, Writing – review & editing. ZX: Conceptualization, Resources, Supervision, Validation, Visualization, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the Seventh Affiliated Hospital, Sun Yat-Sen University Clinical Research 735 program (No. ZSQY735202210).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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

ATA, American Thyroid Association; AUC, Area under the curve; BI-RADS, Breast Imaging Reporting and Data System; CDFI, Color Doppler flow imaging; CI, Confidence interval; CNB, Core-needle biopsy; FNAB, Fine-needle aspiration biopsy; LN, Lymph node; Node-RADS, Node Reporting and Data System; NPV, Negative predictive value; PI-RADS, Prostate Imaging Reporting and Data System; PPV, Papillary thyroid carcinoma; PTC, Positive predictive value; ROC, Receiver operating characteristic; US, Ultrasonography.

References

1. Haugen BR, Alexander EK, Bible KC, Doherty GM, Mandel SJ, Nikiforov YE, et al. 2015 american thyroid association management guidelines for adult patients with thyroid nodules and differentiated thyroid cancer: the american thyroid association guidelines task force on thyroid nodules and differentiated thyroid cancer. Thyroid. (2016) 26:1–133. doi: 10.1089/thy.2015.0020

PubMed Abstract | Crossref Full Text | Google Scholar

2. Li M, Dal Maso L, and Vaccarella S. Global trends in thyroid cancer incidence and the impact of overdiagnosis. Lancet Diabetes Endocrinol. (2020) 8:468–70. doi: 10.1016/S2213-8587(20)30115-7

PubMed Abstract | Crossref Full Text | Google Scholar

3. Bray F, Laversanne M, Sung H, Ferlay J, Siegel RL, Soerjomataram I, et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2024) 74:229–63. doi: 10.3322/caac.21834

PubMed Abstract | Crossref Full Text | Google Scholar

4. Sawka AM, Carty SE, Haugen BR, Hennessey JV, Kopp PA, Pearce EN, et al. American thyroid association guidelines and statements: past, present, and future. Thyroid. (2018) 28:692–706. doi: 10.1089/thy.2018.0070

PubMed Abstract | Crossref Full Text | Google Scholar

5. Yu P, Wang C, Zhang H, Zheng G, Jia C, Liu Z, et al. Deep learning-based automatic pipeline system for predicting lateral cervical lymph node metastasis in patients with papillary thyroid carcinoma using computed tomography: A multi-center study. Chin J Cancer Res. (2024) 36:545–61. doi: 10.21147/j.issn.1000-9604.2024.05.07

PubMed Abstract | Crossref Full Text | Google Scholar

6. Li L, Cheng SN, Zhao YF, Wang XY, Luo DH, and Wang Y. Diagnostic accuracy of single-source dual-energy computed tomography and ultrasonography for detection of lateral cervical lymph node metastases of papillary thyroid carcinoma. J Thorac Dis. (2019) 11:5032–41. doi: 10.21037/jtd.2019.12.45

PubMed Abstract | Crossref Full Text | Google Scholar

7. Zhao H and Li H. Meta-analysis of ultrasound for cervical lymph nodes in papillary thyroid cancer: Diagnosis of central and lateral compartment nodal metastases. Eur J Radiol. (2019) 112:14–21. doi: 10.1016/j.ejrad.2019.01.006

PubMed Abstract | Crossref Full Text | Google Scholar

8. Yang J, Zhang F, and Qiao Y. 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

PubMed Abstract | Crossref Full Text | Google Scholar

9. Wei B, Yao J, Peng C, Zhao S, Wang H, Wang L, et al. Clinical features and imaging examination assessment of cervical lymph nodes for thyroid carcinoma. BMC Cancer. (2023) 23:1225. doi: 10.1186/s12885-023-11721-5

PubMed Abstract | Crossref Full Text | Google Scholar

10. Ekpo EU, Ujong UP, Mello-Thoms C, and McEntee MF. Assessment of interradiologist agreement regarding mammographic breast density classification using the fifth edition of the BI-RADS Atlas. AJR Am J Roentgenol. (2016) 206:1119–23. doi: 10.2214/AJR.15.15049

PubMed Abstract | Crossref Full Text | Google Scholar

11. An JY, Unsdorfer KML, and Weinreb JC. BI-RADS, C-RADS, CAD-RADS, LI-RADS, lung-RADS, NI-RADS, O-RADS, PI-RADS, TI-RADS: reporting and data systems. Radiographics. (2019) 39:1435–6. doi: 10.1148/rg.2019190087

PubMed Abstract | Crossref Full Text | Google Scholar

12. Weinreb JC, Barentsz JO, Choyke PL, Cornud F, Haider MA, Macura KJ, et al. PI-RADS prostate imaging - reporting and data system: 2015, version 2. Eur Urol. (2016) 69:16–40. doi: 10.1016/j.eururo.2015.08.052

PubMed Abstract | Crossref Full Text | Google Scholar

13. Werner RA, Hartrampf PE, Fendler WP, Serfling SE, Derlin T, Higuchi T, et al. Prostate-specific membrane antigen reporting and data system version 2.0. Eur Urol. (2023) 84:491–502. doi: 10.1016/j.eururo.2023.06.008

PubMed Abstract | Crossref Full Text | Google Scholar

14. Elsholtz FHJ, Asbach P, Haas M, Becker M, Beets-Tan RGH, Thoeny HC, et al. Introducing the Node Reporting and Data System 1.0 (Node-RADS): a concept for standardized assessment of lymph nodes in cancer. Eur Radiol. (2021) 31:6116–24. doi: 10.1007/s00330-020-07572-4

PubMed Abstract | Crossref Full Text | Google Scholar

15. Lu F, Zhao Y, Wang Z, and Feng N. Diagnostic accuracy of node-RADS for the detection of lymph node invasion: a systematic review and meta-analysis. Eur Radiol. (2025) 35:4429–37. doi: 10.1007/s00330-024-10867-5

PubMed Abstract | Crossref Full Text | Google Scholar

16. Bai X, Peng C, Liu B, Zhou S, Liu H, Chen Y, et al. Diagnostic and prognostic value of MRI-based Node-RADS for the assessment of regional lymph node metastasis in renal cell carcinoma. Diagn Interv Imaging. (2025) 106:117–25. doi: 10.1016/j.diii.2024.10.005

PubMed Abstract | Crossref Full Text | Google Scholar

17. Yang X, Yang J, Li J, Leng J, Qiu Y, and Ma X. Diagnostic performance of Node Reporting and Data System magnetic resonance imaging score in detecting metastatic cervical lymph nodes of nasopharyngeal carcinoma. Clin Med Insights Oncol. (2024) 18:11795549241231564. doi: 10.1177/11795549241231564

PubMed Abstract | Crossref Full Text | Google Scholar

18. Wu Q, Lou J, Liu J, Dong L, Wu Q, Wu Y, et al. Performance of node reporting and data system (node-RADS): a preliminary study in cervical cancer. BMC Med Imaging. (2024) 24:28. doi: 10.1186/s12880-024-01205-8

PubMed Abstract | Crossref Full Text | Google Scholar

19. Ninkova RV, Calabrese A, Curti F, Riccardi S, Gennarini M, Miceli V, et al. The performance of the node reporting and data system 1.0 (Node-RADS) and DWI-MRI in staging patients with cervical carcinoma according to the new FIGO classification (2018). Radiol Med. (2024) 129:1062–75. doi: 10.1007/s11547-024-01824-9

PubMed Abstract | Crossref Full Text | Google Scholar

20. Fu Y, Cui LG, Ma JY, Fang M, Lin YX, and Li N. Development of a novel contrast-enhanced ultrasound-based nomogram for superficial lymphadenopathy differentiation: postvascular phase value. Ultrasound Med Biol. (2024) 50:852–9. doi: 10.1016/j.ultrasmedbio.2024.02.009

PubMed Abstract | Crossref Full Text | Google Scholar

21. Zhu Y, Meng Z, Fan X, Duan Y, Jia Y, Dong T, et al. Deep learning radiomics of dual-modality ultrasound images for hierarchical diagnosis of unexplained cervical lymphadenopathy. BMC Med. (2022) 20:269. doi: 10.1186/s12916-022-02469-z

PubMed Abstract | Crossref Full Text | Google Scholar

22. Peduzzi P, Concato J, Feinstein AR, and Holford TR. Importance of events per independent variable in proportional hazards regression analysis. II. Accuracy and precision of regression estimates. . J Clin Epidemiol. (1995) 48:1503–10. doi: 10.1016/0895-4356(95)00513-7

PubMed Abstract | Crossref Full Text | Google Scholar

23. Peduzzi P, Concato J, Kemper E, Holford TR, and Feinstein AR. A simulation study of the number of events per variable in logistic regression analysis. J Clin Epidemiol. (1996) 49:1373–9. doi: 10.1016/S0895-4356(96)00236-3

PubMed Abstract | Crossref Full Text | Google Scholar

24. Concato J, Peduzzi P, Holford TR, and Feinstein AR. Importance of events per independent variable in proportional hazards analysis. I. Background, goals, and general strategy. J Clin Epidemiol. (1995) 48:1495–501. doi: 10.1016/0895-4356(95)00510-2

PubMed Abstract | Crossref Full Text | Google Scholar

25. Zu DM, Feng LL, Zhang L, Ma SL, and Zhu YC. Evaluation of mesenteric lymph nodes in a pediatric population with mesenteric lymphadenitis using superb microvascular imaging. Med Sci Monit. (2019) 25:5336–42. doi: 10.12659/MSM.914903

PubMed Abstract | Crossref Full Text | Google Scholar

26. Aladag Kurt S, Kayadibi Y, Onur I, Uslu Besli L, Necati Sanli A, and Velidedeoglu M. Predicting axillary nodal metastasis based on the side of asymmetrical cortical thickening in breast cancer: Evaluation with grayscale and microvascular imaging findings. Eur J Radiol. (2023) 158:110643. doi: 10.1016/j.ejrad.2022.110643

PubMed Abstract | Crossref Full Text | Google Scholar

27. Usmani S, Ahmed N, Al Saleh N, abu Huda F, Amanguno HG, Amir T, et al. The clinical utility of combining pre-operative axillary ultrasonography and fine needle aspiration cytology with radionuclide guided sentinel lymph node biopsy in breast cancer patients with palpable axillary lymph nodes. Eur J Radiol. (2015) 84:2515–20. doi: 10.1016/j.ejrad.2015.10.003

PubMed Abstract | Crossref Full Text | Google Scholar

28. Ferroni G, Sabeti S, Abdus-Shakur T, Scalise L, Carter JM, Fazzio RT, et al. Noninvasive prediction of axillary lymph node breast cancer metastasis using morphometric analysis of nodal tumor microvessels in a contrast-free ultrasound approach. Breast Cancer Res. (2023) 25:65. doi: 10.1186/s13058-023-01670-z

PubMed Abstract | Crossref Full Text | Google Scholar

29. Lucciola S, Pisciotti ML, Frisenda M, Magliocca F, Gentilucci A, Del Giudice F, et al. Predictive role of node-rads score in patients with prostate cancer candidates for radical prostatectomy with extended lymph node dissection: comparative analysis with validated nomograms. Prostate Cancer Prostatic Dis. (2023) 26:379–87. doi: 10.1038/s41391-022-00564-z

PubMed Abstract | Crossref Full Text | Google Scholar

30. Pediconi F, Maroncelli R, Pasculli M, Galati F, Moffa G, Marra A, et al. Performance of MRI for standardized lymph nodes assessment in breast cancer: are we ready for Node-RADS? Eur Radiol. (2024) 34:7734–45. doi: 10.1007/s00330-024-10828-y

PubMed Abstract | Crossref Full Text | Google Scholar

31. Meyer HJ, Schnarkowski B, Pappisch J, Kerkhoff T, Wirtz H, Hohn AK, et al. CT texture analysis and node-RADS CT score of mediastinal lymph nodes - diagnostic performance in lung cancer patients. Cancer Imaging. (2022) 22:75. doi: 10.1186/s40644-022-00506-x

PubMed Abstract | Crossref Full Text | Google Scholar

32. Maggialetti N, Greco CN, Lucarelli NM, Morelli C, Cianci V, Sasso S, et al. Applications of new radiological scores: the Node-rads in colon cancer staging. Radiol Med. (2023) 128:1287–95. doi: 10.1007/s11547-023-01703-9

PubMed Abstract | Crossref Full Text | Google Scholar

33. Ni X, Xu S, Zhan W, and Zhou W. A risk stratification model for metastatic lymph nodes of papillary thyroid cancer: a retrospective study based on sonographic features. Front Endocrinol (Lausanne). (2022) 13:942569. doi: 10.3389/fendo.2022.942569

PubMed Abstract | Crossref Full Text | Google Scholar

34. Lee S, Lee JY, Yoon RG, Kim JH, and Hong HS. The value of microvascular imaging for triaging indeterminate cervical lymph nodes in patients with papillary thyroid carcinoma. Cancers (Basel). (2020) 12:2839. doi: 10.3390/cancers12102839

PubMed Abstract | Crossref Full Text | Google Scholar

35. Huang T, Huang PT, Luo ZY, Lv JF, Jin PL, Zhang T, et al. Use superb microvascular imaging to diagnose and predict metastatic cervical lymph nodes in patients with papillary thyroid carcinoma. J Cancer Res Clin Oncol. (2024) 150:268. doi: 10.1007/s00432-024-05770-x

PubMed Abstract | Crossref Full Text | Google Scholar

36. Yasuoka H, Nakamura Y, Zuo H, Tang W, Takamura Y, Miyauchi A, et al. VEGF-D expression and lymph vessels play an important role for lymph node metastasis in papillary thyroid carcinoma. Mod Pathol. (2005) 18:1127–33. doi: 10.1038/modpathol.3800402

PubMed Abstract | Crossref Full Text | Google Scholar

37. Lee JH, Shin HJ, Yoon JH, Kim EK, Moon HJ, Lee HS, et al. Predicting lymph node metastasis in patients with papillary thyroid carcinoma by vascular index on power Doppler ultrasound. Head Neck. (2017) 39:334–40. doi: 10.1002/hed.24592

PubMed Abstract | Crossref Full Text | Google Scholar

38. García-Silva S and Peinado H. Mechanisms of lymph node metastasis: An extracellular vesicle perspective. Eur J Cell Biol. (2024) 103:151447. doi: 10.1016/j.ejcb.2024.151447

PubMed Abstract | Crossref Full Text | Google Scholar

39. Yu XM, Lo CY, Chan WF, Lam KY, Leung P, and Luk JM. Increased expression of vascular endothelial growth factor C in papillary thyroid carcinoma correlates with cervical lymph node metastases. Clin Cancer Res. (2005) 11:8063–9. doi: 10.1158/1078-0432.CCR-05-0646

PubMed Abstract | Crossref Full Text | Google Scholar

40. Liu Z and Li C. Correlation of lymph node metastasis with contrast-enhanced ultrasound features, microvessel density and microvessel area in patients with papillary thyroid carcinoma. Clin Hemorheol Microcirc. (2022) 82:361–70. doi: 10.3233/CH-221545

PubMed Abstract | Crossref Full Text | Google Scholar

41. Loch FN, Beyer K, Kreis ME, Kamphues C, Rayya W, Schineis C, et al. Diagnostic performance of Node Reporting and Data System (Node-RADS) for regional lymph node staging of gastric cancer by CT. Eur Radiol. (2024) 34:3183–93. doi: 10.1007/s00330-023-10352-5

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: papillary thyroid carcinoma, cervical lymph node, ultrasonography, Node reporting and data system, risk stratification

Citation: Zhang H, He D, Mao M, Li R, Yin T, Zheng B-w, Zheng J, Wu T, Ren J and Xu Z (2026) Standardizing cervical lymph node evaluation in papillary thyroid carcinoma: diagnostic accuracy of Node-RADS on ultrasound. Front. Endocrinol. 17:1709522. doi: 10.3389/fendo.2026.1709522

Received: 20 September 2025; Accepted: 09 January 2026; Revised: 25 November 2025;
Published: 28 January 2026.

Edited by:

Terry Francis Davies, Icahn School of Medicine at Mount Sinai, United States

Reviewed by:

İlhan Hekimsoy, Ege University, Türkiye
Yanjun Gao, George Washington University, United States
Jianxue Liu, Baoji Central Hospital, China

Copyright © 2026 Zhang, He, Mao, Li, Yin, Zheng, Zheng, Wu, Ren and Xu. 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: Jie Ren, cmVuakBtYWlsLnN5c3UuZWR1LmNu; Zuofeng Xu, eHV6ZkBtYWlsLnN5c3UuZWR1LmNu

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.