- 1Department of Endocrinology, Huzhou First People’s Hospital, The First Affiliated Hospital of Huzhou University, Huzhou, Zhejiang, China
- 2School of Medicine (School of Nursing), Huzhou University, Huzhou, Zhejiang, China
- 3Department of Hepatobiliary and Pancreatic Surgery, Huzhou Central Hospital affiliated to Huzhou University, Huzhou, Zhejiang, China
Background: Nowadays, artificial intelligence (AI) diagnostic models based on ultrasound features have been gradually integrated into the evaluation of thyroid nodules. However, the diagnostic effects of different AI-assisted diagnosis methods vary greatly.
Objective: This study aims to systematically evaluate the performance of the ultrasound-based artificial intelligence diagnostic models in differentiating benign and malignant thyroid nodules and to determine the most effective diagnostic model.
Methods: We conducted a comprehensive literature search in PubMed, Web of Science, and the Cochrane Library using subject-specific keywords to identify studies on AI-assisted thyroid nodule diagnosis. Study quality was assessed using Quality Assessment of Diagnostic Accuracy Studies-2 (QUADAS-2). Meta-analysis was performed using Meta-Disc 1.4, Review Manager 5.4, R 4.4.2, and Stata 17.0. Pooled sensitivity, specificity, diagnostic odds ratio (DOR), and area under the summary receiver operating characteristic curve (SROC-AUC) with 95% confidence intervals (CI) were calculated. Subgroup analyses and clinical applicability assessments were conducted.
Results: Twenty-eight studies involving 134,028 patients, 158,161 thyroid nodules, and 529,479 ultrasound images were included. The AI-assisted diagnostic system demonstrated high diagnostic performance: pooled sensitivity = 0.89 (95% CI: 0.87–0.91), specificity = 0.84 (0.80–0.88), positive likelihood ratio (PLR) = 5.60 (4.40–7.20), negative likelihood ratio (NLR) = 0.13 (0.10–0.16), DOR = 43.94 (30.11–64.14), and SROC-AUC = 0.93 (0.91–0.95). The threshold effect analysis (Spearman correlation = -0.18, P > 0.05) indicated no significant heterogeneity. The diagnostic accuracy is higher in Asian countries, in prospective and multicenter designs, with external validation sets, without cross-validation, with deep learning, and in postoperative patient subgroups. Additionally, improved performance was observed in cohorts with smaller nodule diameters (<20 mm), higher malignancy rates, older patient age (≥50 years), and higher female proportions, though heterogeneity remained significant. Univariate and multivariate meta-regression analyses identified AI type, malignancy rate of nodules as significant sources of heterogeneity. Notably, the EDLC-TN model showed the highest diagnostic accuracy.
Conclusion: AI-assisted diagnostic techniques demonstrate significant potentialin thyroid nodule evaluation, with the EDLC-TN model showing particularly high clinical utility. Optimal diagnostic performance was observed for nodules <20 mm in diameter and in patients aged ≥50 years.
Systematic review registration: https://www.crd.york.ac.uk/PROSPERO/view/CRD42024581421, identifier CRD42024581421.
1 Introduction
Thyroid nodules are localized, abnormal hyperplastic masses in the thyroid tissue, which are common and frequently occurring diseases. They mainly include nodular goiter, neoplastic nodules, cystic lesions, and inflammatory nodules. With the advantages of being economical and noninvasive, ultrasonography has been widely used for the initial diagnosis and follow-up of thyroid nodules and remains the diagnostic tool of choice in clinical practice (1, 2). The most important thing in the diagnosis and treatment of thyroid nodules is to differentiate between their benign and malignant nature. With the wide application of high-resolution ultrasound technology and the increase in people’s health awareness, the detection rate of thyroid nodules is increasing year by year. However, the traditional ultrasound diagnostic mode is highly dependent on the experienced judgment of radiologists, which has inherent defects such as strong subjectivity and limited diagnostic efficiency (3). In the last 5 years, the rapid development of ultrasound-based artificial intelligence diagnostic systems has brought new methods for the diagnosis of thyroid nodules. This artificial intelligence smart diagnostic system has high accuracy and stability. Practical applications have shown that this auxiliary diagnostic system can effectively reduce the overuse of fine-needle aspiration biopsy (FNA) due to subjective judgment errors and unnecessary surgical operations while significantly shortening the waiting time for patients’ diagnosis and treatment (4, 5). In addition, the technique has good diagnostic accuracy in thyroid nodule diagnosis. Currently, intelligent assistive systems for thyroid nodule diagnosis include various methods such as random forest (RF), machine learning (ML), and convolutional neural network (CNN). However, there are significant differences in diagnostic efficacy among different technical routes, and their clinical application value still needs to be further systematically evaluated.
Existing studies have shown that applying this technology to the clinical diagnosis of thyroid nodules can achieve high accuracy. It is worth noting that due to the relatively short application time of this technology, there is still some controversy in the academic community, and the clinical results vary. For example, the results of the diagnostic model used by WeiX’s team showed a high specificity of 94%, while the results of the diagnostic model used by BudaM’s team showed a specificity of only 52%, which is a significant difference (6). Based on the data from the studies we included, we found a large gap in specificity between the different diagnostic models. Therefore, this study intends to use meta-analysis methods to systematically evaluate the overall performance of intelligent auxiliary diagnosis models based on ultrasound images in the diagnosis of thyroid nodules, to screen out the model with the best diagnostic performance, and to provide a scientific basis for the clinical practice of thyroid cancer.
2 Methods
2.1 Search strategy and selection criteria
This study was conducted in strict adherence to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines, and the systematic review was registered on PROSPERO at https://www.crd.york.ac.uk/PROSPERO/view/CRD42024581421. A total of 2,392 articles were identified through a literature search across three databases (PubMed: 1,558; Web of Science: 800; Cochrane Library: 34; database inception to August 31, 2024). After removing 620 duplicate records, an additional screening of titles and abstracts excluded studies that did not meet the inclusion criteria, such as those lacking data (n = 26), being irrelevant (n = 1533), conference proceedings (n = 6), or letters (n = 15). This resulted in 192 articles eligible for full-text evaluation. Following a detailed assessment of the complete texts, further exclusions were made due to the unavailability of full texts (n = 1), incomplete data (n = 97), or irrelevance (n = 66). Ultimately, 28 studies were included in the analysis (6–33). We searched using keywords and medical subject terms (MeSH terms): (Artificial Intelligence OR Machine Learning OR Convolutional Neural Network) AND (Thyroid Nodule OR Thyroid Cancer OR Thyroid Neoplasms). Specific search formulas are provided in the Appendix. Two researchers independently performed the literature selection process (Figure 1), with discrepancies resolved through discussions during a four-person consensus meeting.
Inclusion criteria: (1) All patients were studied as patients with thyroid nodules, and all underwent ultrasonography; (2) Fine-needle aspiration (FNA) or surgical biopsy was used as the “gold standard” for diagnosing benign and malignant thyroid nodules; (3) Data were complete, and four-quadrant contingency tables could be extracted; (4) Articles generally adhered to the STARD (Standards for Reporting of Diagnostic Accuracy) statement for diagnostic accuracy studies in thyroid cancer research, and (5) At least one artificial intelligence (AI) model was included. Exclusion criteria: (1) non-original research articles, including reviews, meta-analyses, conference proceedings, editorials, letters, and case reports; (2) unavailability of full text; (3) non-English-language publications; and (4) studies with overlapping study populations.
2.2 Data extraction
We arranged for two researchers to extract data from each original study and cross-check it independently. Any discrepancies were resolved through discussions within the research team. Specific information extracted from each original study included: first author, year of publication, country, study design, data source, number of cases, study center, malignancy rate of thyroid nodules, number of ultrasound images, reference standard, mean age, proportion of females, mean nodule diameter (in mm), type of artificial intelligence (AI) model, validation set type, depth of study, number of nodules, subgroup status, and diagnostic performance metrics (such us true positives [TP], false positives [FP], false negatives [FN], true negatives [TN]). In this study, for data selection, all external validation sets were used for data extraction and meta-analysis. When multiple external validation sets were available for evaluating AI-assisted diagnostic performance, only the largest cohort was selected for analysis. For nodule attributes, we focused on characteristic data of malignant nodules, including their diameter size, population age, and percentage of females (Tables 1–3).

Table 2. Study characteristics and population statistics according to malignant nodule class (test group nodules).
2.3 Quality assessment
The methodological quality of the included studies was assessed using the Quality Assessment of Diagnostic Accuracy Studies-2 (QUADAS-2) tool, which evaluates the risk of bias and applicability concerns across four domains for the 28 articles. Each risk shift associated with a question is evaluated as “Yes,” “No,” or “Uncertain,” and as “High,” “Low,” or “Uncertain” in terms of its applicability. This assessment was independently completed by two researchers. Any differences that emerged were resolved by the research members through group discussions. Based on the evaluation results, the QUADAS-2 scale was completed, and the quality assessment was subsequently performed using Review Manager 5.4 (Figure 2).
2.4 Statistical analysis
We used Meta-Disc 1.4, Review Manager 5.4, R 4.4.2, and Stata 17.0 software for statistical analysis. Diagnostic efficacy was evaluated by constructing and summarizing ROC (SROC) curves by combining sensitivity, specificity, diagnostic advantage ratio (DOR), and AUC at 95% CI. The Deeks’ Funnel Plot Asymmetry Test was then used to assess publication bias. The Spearman correlation coefficient was also used to evaluate the threshold effect between studies. For heterogeneity evaluation, we used Cochran’s Q-test and I-squared (I2) statistic. Univariate and multivariate meta-regression analyses were employed to identify the source of heterogeneity, and subgroup analyses were conducted. Variables that produced heterogeneity were analyzed, and p-values < 0.05 were considered statistically significant.
3 Results
3.1 Literature search and study characteristics
After a comprehensive literature review and checking, 192 records were finally obtained. Unavailability of full text (1), incomplete data (97), and irrelevant studies (66) were excluded. Finally, a total of 28 studies that met the criteria were included in this meta-analysis (Figure 1). Of the 28 included studies, all were in English, 21 studies were from China (75%), 4 studies were from Korea (14%), 2 studies were from the United States (7%), and 1 study was from Spain (4%).25 studies were from Asian countries (89%), and 3 studies were from Western countries (11%).2 studies were prospective studies (7%), and 26 studies were retrospective (93%). In terms of diagnostic models, these included CNN 17 studies (61%), ANN 2 studies (7%), RF 3 studies (11%), and Mixed Model Group 6 studies (21%). Twelve studies employed a single-center design (43%), while 16 studies utilized a multicenter design (57%). 16 studies used an external validation set (57%), and 12 studies used an internal validation set (43%). Seventeen studies performed cross-validation (61%), while 11 studies did not (39%). 21 studies employed deep learning (75%), and 7 studies used non-deep learning (25%). 11 studies selected only cases after thyroid surgery (39%), and 17 studies were FNA cases (61%). Sixteen studies were thyroid nodules with a mean diameter of < 20mm. Nine studies included patients with a mean age of 50 years or older. The malignant rate of thyroid nodules was ≥50% (50%) in 14 studies and < 50% (50%) in the remaining 14 studies (Tables 4, 5). In addition, this meta-analysis included a total of 134,028 patients, 158,161 thyroid nodules, and 529,479 ultrasound images of thyroid nodules (Tables 1–3).
3.2 Study quality assessment
The QUADAS-2 scale was employed to assess the methodological quality of the 28 articles included in this study. Among these, 11 studies exhibited case selection bias because the cases selected by the original authors were exclusively surgical patients, which may have resulted in a significantly higher risk of thyroid cancer in this group. One study demonstrated index test bias due to using the same data set to train and test the same data set to train and test the AI diagnostic model, comprising a total of 172 nodules. However, since the original authors mitigated overfitting through cross-validation, this study can still be rated as low risk in terms of applicability (21). The results showed that the research quality of the studies included in this meta-analysis was satisfactory (Figure 2).
3.3 Diagnostic accuracy and heterogeneity evaluation
We analyzed the threshold effect by Meta-Disc, and the results showed that the Spearman correlation coefficient value was -0.18, P>0.05, no statistical significance, the threshold effect caused no heterogeneity, and we could continue to merge the effect sizes of the diagnostic models and meta-analysis of related data.
Among the 28 studies included in the meta-analysis, the diagnostic performance of AI-assisted diagnostic technology for thyroid nodules was evaluated. The results indicated that the AI-assisted diagnostic system demonstrated good diagnostic efficacy for distinguishing benign and malignant thyroid nodules: the pooled sensitivity was 0.89 (95% CI: 0.87–0.91); the pooled specificity was 0.84 (95% CI: 0.80–0.88); the positive likelihood ratio (PLR) was 5.60 (95% CI: 4.40–7.20); the negative likelihood ratio (NLR) was 0.13 (95% CI: 0.10–0.16); and the diagnostic odds ratio (DOR) was 43.94 (95% CI: 30.11–64.14) (Table 6). Additionally, the area under the curve (AUC) for the summary receiver operating characteristic (SROC) plot was 0.93 (95% CI: 0.91–0.95) (Figure 3).
Although these results indicate that AI-assisted diagnostic techniques exhibit good diagnostic efficacy for thyroid nodules, significant heterogeneity was observed in the pooled analysis of sensitivity and specificity. In the meta-analysis, the heterogeneity of sensitivity for AI-assisted diagnostic techniques was I² = 77.05% (95% CI: 68.84%-85.27%), P = 0, while the heterogeneity of specificity was I² = 94.34% (95% CI: 92.99%-95.69%), P = 0. These findings suggest that although the diagnostic sensitivity and specificity of the thyroid AI-assisted diagnostic system are high, there is substantial heterogeneity among studies (I² > 50%, P = 0) (Figure 4A). The forest plot of the diagnostic odds ratio (DOR) demonstrated a pooled DOR value of 43.94 and a diagnostic score of 3.78, indicating that the AI diagnostic model has strong diagnostic value and can be applied to evaluate individual clinical cases (Figure 4B). In addition, we did forest plots of positive likelihood ratio (PLR) and negative likelihood ratio (NLR), and we found that positive likelihood ratio = 5.61 and negative likelihood ratio = 0.13, which indicates that the AI diagnostic model has a high diagnostic value for thyroid nodules (Figure 4C).

Figure 4. (A) Forest plot for sensitivity and specificity after combination. (B) Forest plot for diagnostic odds ratio and diagnostic score after combination. (C) Forest plot for likelihood ratio after combination (LR+, LR-).
Given the large number of studies included, there may be other factors influencing the overall results. Therefore, we conducted a sensitivity analysis (Figure 5B), sequentially excluding individual studies, and did not find any study with apparent heterogeneity. To identify heterogeneity, we also performed a Bivariate Boxplot analysis, which indicated that four articles exhibited significant heterogeneity: 1 (BudaM2019), 15 (SunC2020), 16 (VelascoPF2024), and 28 (ZhuYC2022) (Figure 5D). Sensitivity analysis results indicated that the four individual studies did not significantly impact the stability of the results. Therefore, we conducted further sensitivity analysis by simultaneously excluding the four articles identified in the Bivariate Boxplot analysis. After exclusion, the sensitivity I² was 73.47% (95% CI: 62.79%–84.15%), P = 0.00. The specificity I² was 93.36% (95% CI: 91.56%–95.16%), P = 0.00, indicating no significant changes in overall sensitivity and specificity heterogeneity. Thus, the results of the 28 studies included in this analysis demonstrate high stability across studies (Supplementary Figure S1). We conducted a publication bias test, with the Deeks’ Funnel Plot of Asymmetry Test yielding a p-value of 0.28, indicating that the studies were generally symmetrical. We also performed an Egger test (p = 0.992), which confirmed the absence of publication bias across the studies (Figure 5C). Given that all included studies in this research were of high quality and stability, the AI type of these four articles does not substantially impact the results, and the findings of this study are highly reliable.

Figure 5. (A) Diagnostic Odds Ratio of different AI-aided diagnostic models after combination. (B) Result of Sensitivity analysis. (C) Results of Deeks’ Funnel Plot of Asymmetry Test for publication bias. (D) Bivariate Boxplot for diagnosis of thyroid nodules by AI.
After excluding the influence of threshold effects on heterogeneity, we conducted a subgroup univariable meta-regression analysis based on the completeness of the collected data, including (country/region, study design, number of centers, validation type, cross-validation, deep learning, reference standard) to determine the sources of heterogeneity in sensitivity and specificity. We found that six subgroups, excluding region, had statistically significant effects on the heterogeneity of sensitivity and specificity. The results of the joint model showed that no subgroup had a statistically significant impact on the heterogeneity of sensitivity and specificity. The diagnostic performance results of the AI diagnostic model indicated that Asian countries, prospective studies, multicenter studies, studies using external validation sets, non-cross-validation, deep learning models, and the subgroup of patients who underwent thyroid surgery had higher diagnostic accuracy, but also exhibited higher heterogeneity. AI models demonstrated high diagnostic performance in Asian populations, with a sensitivity of 0.89 (0.87–0.91) and a specificity of 0.85 (0.82–0.88). In prospective studies, AI models performed notably well, with an AUC of 0.97 (0.95–0.98). In multicenter studies, the AI diagnostic model demonstrated high sensitivity (0.90, 95% CI 0.88–0.93) and specificity (0.85, 95% CI 0.81–0.89). Additionally, when an external validation set was available, the diagnostic model exhibited high diagnostic performance, with an AUC of 0.94 (95% CI 0.91–0.95). The non-cross-validation group had a higher AUC of 0.95 (0.92–0.96). In the deep learning subgroup, the deep learning group had a high AUC of 0.94 (0.91–0.95). Finally, the subgroup of patients who underwent surgery alone had high diagnostic accuracy, with a pooled DOR of 47.32 (27.99–80.02) (Table 4).
We also conducted univariable meta-regression analysis on the remaining eight subgroups (number of patients, overall malignancy rate, AI type, nodule count, nodule malignancy rate, nodule diameter, patient age, and female proportion) using univariable meta-regression analysis. We found that, except for nodule count, which had a statistically significant effect on sensitivity heterogeneity but no statistically significant effect on specificity heterogeneity, the remaining seven subgroups all had statistically significant effects on both sensitivity and specificity heterogeneity (Figure 6, Table 5). The results of the joint model showed that the four subgroups—thyroid nodule malignancy rate, nodule size, patient age, and female proportion — had statistically significant effects on the heterogeneity of sensitivity and specificity. The results of the AI diagnostic model’s diagnostic performance showed that the diagnostic performance of the AI diagnostic model was stronger when the number of patients and nodules was higher, the nodule malignancy rate was higher, the convolutional neural network (CNN) was used, the nodule diameter was smaller, and the patient age and female proportion were higher. The AI model demonstrated high diagnostic performance in studies with a large sample size (≥500 patients), with an AUC of 0.93 (0.91–0.95). When the thyroid nodule malignancy rate was ≥50%, the AI model exhibited a sensitivity of 0.90 (0.88–0.92) and specificity of 0.88 (0.84–0.92). Among AI model subtypes, the convolutional neural network (CNN) model performed best, with a sensitivity of 0.90 (0.88–0.92) and specificity of 0.85 (0.80–0.89). The more nodules included in the study, the better the diagnostic performance of the AI model, with an AUC of 0.94 (0.91–0.96). The higher the malignancy rate of the nodules, the better the diagnostic performance of the model, with an AUC of 0.94 (0.91–0.96). The AI-assisted diagnostic model demonstrated higher sensitivity (0.90 [0.87–0.92]) and specificity (0.86 [0.81–0.89]) for thyroid nodules with an average diameter <20 mm. For patients with thyroid nodules aged 50 years or older, the AI-assisted diagnostic model has higher diagnostic value, although the difference is not statistically significant. The AI-assisted diagnostic model has higher sensitivity (0.90, 95% CI: 0.87–0.92) and specificity (0.85, 95% CI: 0.81–0.89) in populations with a higher proportion of female patients (Table 5). Finally, considering the heterogeneity of the included studies, we conducted a multi-factor meta-regression analysis by including all variables in the model, in addition to the results of the combined model. The results showed that nodule malignancy rate, and AI type had statistically significant heterogeneity in sensitivity and specificity, and were significant sources of heterogeneity (Tables 4, 5).
We divided the AI-assisted diagnostic system into four subgroups: CNN, RF, ANN, and Mixed Model Group. By comparing the DOR and weights of each subgroup, we found significant heterogeneity within the CNN and Mixed Model Group subgroups. The diagnostic odds ratio (DOR) indicates the strength of the association between the diagnostic test results and the disease; a higher value indicates better discriminatory performance. Since the AI diagnostic models differ across subgroups, we compared and analyzed the overall DOR values between different diagnostic model subgroups, finding significant differences between groups. The RF group had the lowest DOR value (22.83), while the CNN group had the highest DOR value (51.66). The AI diagnostic model of the CNN group demonstrated a significantly higher value for distinguishing benign from malignant thyroid nodules compared to the Mixed Model Group, RF, and ANN groups (Figures 7, 5A).
Finally, we evaluated the clinical applicability of thyroid nodules diagnosed using AI-assisted ultrasound diagnosis technology. The results showed that when the pre-test probability was set at 30.00%, the positive likelihood ratio (PLR) was 6.00, and the post-test probability of a positive test result reached 71.00%. The NLR was 0.13, and the post-test probability for negative test results was only 5.00%, indicating that the thyroid AI-assisted diagnostic system has high clinical predictive value (Figure 8). Diagnostic performance was visualized using a likelihood ratio scattergram, where PLR < 10.00 and NLR > 0.10 indicate higher diagnostic accuracy. The combined effect size of all diagnostic models showed that the intelligent auxiliary diagnostic model performed poorly in terms of diagnosing and ruling out malignant thyroid nodules. However, as shown in the figure, the intelligent diagnostic models used in the five studies (11, 13, 14, 18, and 23) demonstrated very high accuracy for thyroid nodules, with the models used being PLS-DA, AIBx V2 (DCNN, ResNet 34), DCAD (FCN, DCNN), EDLC-TN (DCNN, DenseNet), and XGBoost (Figure 9). After comprehensively comparing the sensitivity, specificity, DOR, PLR, and NLR of these five models, we found that the model used in the WeiX2020 (23) study had a sensitivity of 0.93 (95% CI: 0.90–0.96), specificity of 0.94 (95% CI: 0.92–0.96), Diagnostic Score was 5.37 (95% CI: 2.67–5.37), Odds Ratio was 215.25 (95% CI: 127.73–362.73), DLR Positive was 15.44 (95% CI: 11.33–21.03), and DLR Negative was 0.07 (95% CI: 0.05–0.10) (Figures 4A–C). The model used was EDLC-TN (an ensemble deep learning classification model for thyroid nodules). The ROC curve for the AI-based diagnosis of thyroid nodules showed an AUC of 0.93 (95% CI: 0.91–0.95), indicating that the thyroid intelligent assistance model has high accuracy in distinguishing between benign and malignant nodules (Figure 3).

Figure 8. Fagan nomogram: Evaluation of clinical applicability of AI-assisted diagnostic techniques in thyroid nodules diagnosis.

Figure 9. Likelihood ratio scattergram: Evaluation of clinical applicability of AI-assisted diagnostic techniques in thyroid nodules diagnosis.
4 Discussion
In recent years, artificial intelligence tools have become increasingly common in various disciplines. Artificial intelligence is beneficial in assessing thyroid nodules, primarily for risk stratification and the diagnosis of benign and malignant thyroid nodules. This meta-analysis encompassed a substantial number of studies that examined various types of AI models, spanning from the past to the present. Our goal is to systematically compare the performance of existing artificial intelligence models in identifying benign and malignant thyroid nodules based on ultrasound image features, identify the optimal model, and provide guidance for developing improved intelligent diagnostic models. We systematically investigated sources of heterogeneity and compared the I2 for subgroup analyses with the pooled effect size, I2, which showed no significant changes in sensitivity or specificity across subgroups. Given the obvious heterogeneity suggested by I2 after the combined effect size, we also performed univariate and multivariate meta-regression analyses for different subgroups. All factors except region were sources of heterogeneity, among which the malignant rate of thyroid nodules, and different AI types were significant sources of heterogeneity. At the same time, the diagnostic accuracy of AI models is constrained by multiple factors, all of which increase the complexity of diagnostic model construction (34, 35).
4.1 The quality of ultrasound images
Ultrasound imaging provides standardized, high-contrast visual features of thyroid nodules, including size, shape, density, and blood flow signal patterns, on which both artificial intelligence algorithms and radiologists rely for assessment. The lack of a uniform ultrasound image acquisition protocol across the original studies included in this study, coupled with the variability in institutional operating procedures and equipment specifications (E.g., Ransducer frequencies, gain settings), greatly limits the generalizability of the current model across different clinical Settings. The frequency of the ultrasound probe will affect the resolution of the image; the higher the frequency, the higher the resolution. Different ultrasound devices and probes were used in each original study, which is also a source of potential heterogeneity. For example, in the study of ChenJH et al., L14-3U transducer of Resona 9 device (Mindray, Shenzhen, China) (frequency: 3–9 MHz) and an L12–5 transducer (frequency: 5–12 MHz) from an iU22 device (Philips, Amsterdam, The Netherlands) (9). However, KimYJ et al. used a higher frequency probe, and the Philips iE33 US instrument with a 12–15 MHz linear array transducer (Philips Medical Systems) provided ultra-high resolution imaging for image acquisition (12). A multicenter study showed that ultrasound equipment from different manufacturers had a significant impact on AI model performance. Studies conducted in Asia may have overrepresented certain regional manufacturers (e.g., Mindray) and underrepresented brands commonly used in the West (e.g., General Electric) (36). In addition, the number of phantom ultrasound devices included in each study varies greatly, which can significantly affect the image quality. YaoJC et al. collected ultrasound images of patients using 26 different devices produced by the General Electric Company, Philips, Esaote, Siemens, and Toshiba (25). Heterogeneity of ultrasound systems and operators from different manufacturers may contribute to variability during training. Of course, different types of ultrasound devices were used in the study, which helped to increase the diversity of the data.
Almost all studies were trained and tested on only a representative single image of each nodule, and the reviewed images were static two-dimensional images, ignoring dynamic features such as hemodynamics and tissue elasticity. Most studies have selected transverse and longitudinal sections, while other sections are rarely included. The selection of representative images and semi-automatic segmentation can theoretically be affected by the operator’s experience, and this method may limit the diagnostic ability of observers and artificial intelligence models. For example, Buda et al. used two orthogonal images of transverse and longitudinal sections as model feature extraction images (6). In clinical practice, physicians often rely on real-time ultrasound information to evaluate thyroid nodules rather than a single representative image. Thus, radiologists were able to assess based on more thorough imaging of each thyroid nodule, which may explain the higher performance of experienced radiologists. At present, due to the complexity of 3D model architecture, the limitation of computational efficiency, and the difference in tools, 3D dynamic image models have not become a research hotspot. At present, researchers have developed artificial intelligence models based on contrast-enhanced ultrasound (CEUS) images. For example, the study by ChenJH et al., included in this study, applied machine learning to extract radiomics features of two-dimensional ultrasound (2D-US) combined with contrast-enhanced ultrasound (CEUS) images of the thyroid to classify and predict benign and malignant thyroid nodules, with an AUC of 0.94 (9). In addition, manually labeled ultrasound images can significantly affect the accuracy of the model. The images mainly extract image Texture Features, Higher-Order Features, and other information, and the labeled images affect the model’s ability to learn features, especially small nodules. For example, the study of ZhuJL et al. included marked thyroid nodule images (33), which is also a potential source of heterogeneity and an important factor affecting the accuracy of the model.
4.2 Size of thyroid nodule
Different nodule sizes can affect the recognition ability of AI diagnostic models. In this study, the diameter sizes of nodules in different original studies were quite different. In general, for image feature extraction, the degree of thyroid nodule margin of large nodules cannot be judged because it exceeds the view boundary, which may lead to poor model performance. Chen et al. showed that nodule size and tumor edge roughness were risk factors for detecting malignant nodules (8, 16, 37, 38). This information will be lost when the nodule is too large. However, for small nodules, especially thyroid micropapillary carcinoma smaller than 10mm, the boundary may not be clear, and the AI system may not have enough identifiable features to make a correct prediction. The study by KoSY et al. showed that the inclusion of thyroid nodules of different sizes would affect the results, and their model findings could not guarantee high accuracy in other studies with large groups of different nodules (14). An international study shows that in Asia, more thyroid cancers are detected through asymptomatic screening (51%), whereas in Western countries, more cancers are detected through symptoms (30%). Nodules detected by screening are usually small and have atypical features, which may lead to inconsistent model criteria for small nodules in different regions (39). In addition, the biopsy rate of small nodules is low in clinical practice. According to ATA guidelines, most thyroid nodules less than 10mm are treated conservatively, which leads to insufficient data in the training set and is also a factor causing bias (37).
4.3 Different subtypes of thyroid nodules
Diagnosis of thyroid nodules is not only about identifying benign and malignant nodules, but also requires accurate identification of benign and malignant nodules of various histologic types of the thyroid gland. Existing training datasets are often inadequate in size, especially for rare histologic subtypes such as follicular carcinoma, which may introduce selection bias and affect model performance. This meta-analysis study included only two (25, 28) that sampled follicular thyroid carcinoma, and there are still fewer studies on rare subtypes, and determining the specific type of follicular thyroid nodule is a difficult challenge in accurately diagnosing thyroid cancer. Fine needle aspiration cytology (FNAC) can only provide a vague Bethesda class IV diagnosis with uncertain malignancy and clinical prognosis, most models only included nodules with malignant or benign clear pathologic findings, and most studies excluded thyroid nodules that were undiagnosable or cytologically or histologically indeterminate, such as Bethesda class IV suspected malignant nodules.YaoJC et al. constructed the swin-transformer (ST) model for follicular thyroid nodules with an AUC of 0.90 (25). In the dataset studied by ParkVY et al., 95.5% of the patients had typical papillary thyroid carcinoma with ultrasound features different from those of follicular carcinoma, medullary carcinoma, and lymphoma. However, other subtypes, such as follicular carcinoma and mesenchymal carcinoma, also had some different ultrasound features compared to benign thyroid lesions. The algorithm may be less efficient in diagnosing medullary carcinoma compared to other subtypes because it has more ultrasound features compared to benign thyroid lesions. Therefore, there is still a need to validate the diagnostic model for nodules with indeterminate cytological findings and therefore characterized as “heterogeneous of uncertain significance” or “follicular lesions of uncertain significance” (19). The proportion of papillary thyroid carcinoma may be higher in Asian populations than in the West, whereas other subtypes, such as follicular carcinoma, are relatively rare (40). Since the majority of malignant nodules in the current studies were papillary carcinomas, it is still not possible to evaluate and compare the performance of machine learning algorithms in each cancer subtype. In addition, Hashimoto nodules and chronic atrophic nodules can exhibit similar ultrasound features to malignant nodules, and there are differences in the ultrasound features of different pathologic types, which may decrease the classification performance when the model is applied to regions with different pathologic spectrums (41).
4.4 Regional differences
About 89% of the subjects in this Meta-analysis were from Asia, and only 11% were from Western countries. Although most studies verified the results in different external test sets, the accuracy and clinical applicability of the results were satisfactory; however, there was inevitably racial bias, and this regional difference may be one of the important sources of heterogeneity in the studies. In a statistical study of 25 countries, the incidence of thyroid cancer in Asian countries was significantly higher than that in Western countries, mainly because the popularization of ultrasound screening in Asian countries led to overdiagnosis, especially in middle-aged and elderly individuals older than 50 years, rather than an increase in true disease (39). The increase in incidence may also be related to economic development, population life pressure, inflammation, and other factors. Asian countries have invested heavily in thyroid disease research and accumulated a large number of high-quality thyroid ultrasound image datasets. After rigorous screening and annotation, these datasets provide a solid foundation for model training, enabling the model to identify the characteristics of thyroid nodules in this population more accurately. A recent Chinese study was trained and tested on ultrasound data from nine hospitals (42). Researchers in Asian countries have optimized and adjusted the diagnostic model according to local medical practice and patient characteristics. By adjusting the model’s parameter settings and refining the feature extraction method, the model becomes more suitable for the medical environment and patient population of Asian countries, thereby improving diagnostic accuracy. It may not perform well in Western populations, particularly in thyroid cancer subtypes with distinct distributions. This is consistent with the results of the subgroup analysis of this study, showing that the diagnostic accuracy of the AI diagnostic model was higher for Asian countries. However, the results of subgroup analysis in this meta-analysis showed that the effect of country and region subgroups on heterogeneity was not statistically significant, possibly due to the small number of studies included in Western countries.
4.5 Feature extraction methods
According to the 2015 American Thyroid Association ATA guidelines, it is divided into six types, including mixed type, intermediate type, and peripheral type, only in terms of calcification. At present, suspicious ultrasound features of thyroid cancer are mainly identified by sonographers, and accurate differentiation of these features is difficult and may be misunderstood or missed (37). Feature extraction is primarily divided into two approaches: experienced radiologists and radiomics. The former is based on human feature extraction, while the latter is based on computer-based image analysis. AI models are typically trained on datasets annotated by experienced radiologists, essentially learning to recognize the same diagnostic patterns as those used by human experts (43). In traditional machine learning models, experts in the relevant domain will first set the most applicable features to reduce the complexity of the data and highlight patterns. However, some hidden relationships may be lost when manually entered, and they may not be able to capture complex, non-intuitive relationships in the data in advance. However, these features may be hidden in a simple image texture, and computer-aided diagnosis (CAD) systems are a novel AI technique developed to achieve automatic analysis of ultrasound images (44). ZhaoCK2021, ZhengYX2024, and other studies included in this study all manually extracted ultrasound features (27, 28), while ChenJH2024, LaiM2023, and other studies automatically extracted features through CNN (9, 15). By comparing the feature extraction methods, we found that ultrasound image features extracted by radiologists were significantly better than image-based methods, such as those used by Chen et al., who extracted empirical features by radiologists for TI-RADS classification, with an AUC of 0.97 (8). However, the model of Lai et al. used the PyRadiomics (45) feature extraction package to extract radiomics features for training machine learning models, and its AUC was only 0.856, which was significantly lower than the empirical features extracted from ultrasound images (15). Features automatically extracted by computer-based image analysis models may contain redundant information (such as artifacts in ultrasound images). Instead of selecting radiomics features with partial thyroid nodule signs in a single image, an experienced radiologist reviews a series of thyroid ultrasound images to make the final decision. It has the best correlation with the nature of the nodule. Traditional methods (such as logistic regression + manual features) can explicitly incorporate clinical prior knowledge (such as Bethesda classification criteria), while AI model extraction of features may ignore these key rules (46). Furthermore, the large number of different radiomics features poses an intractable challenge that hinders the use of machine learning models for diagnosing thyroid nodules in the clinic (47).
A recent study by Chen JH et al. used a machine learning model (DT) to construct feature importance combined with LASSO to select features (9). With different feature screening methods, the final remaining features are different, and the accuracy of the final model is also different. Although the malignant characteristics of different types of TI-RADS, such as C-TIRADS and ATA-TIRADS, are basically the same, there are also differences. There may be differences in the interpretation criteria of ultrasound features by radiologists in different regions. When a standard training model for a single region is used, performance may decline in other regions (48). The different sections of thyroid nodules, such as longitudinal section and transverse section, will affect the extraction of ultrasound image features. The number of representative ultrasound features used in each study varied; for example, the study by Chen et al. (8) included 10 ultrasound features, while the study by zhaoCK et al. (27) included only 6 ultrasound features. Finally, the image acquisition process of KoSY et al. ‘s study, such as selecting input images for thyroid nodules and labeling ROIs, relied on radiologists and was essentially operator dependent due to the perception variation of radiologists and the possibility of ignoring obvious features (14). Aspects of the influencing characteristics are a source of potential heterogeneity and also limit further exploration of potential effects on diagnostic efficacy.
4.6 Malignant rate of thyroid nodules
The results of this study showed that the higher the malignant rate of thyroid nodules, the stronger the diagnostic performance of the AI diagnostic model. However, most of the included studies used retrospective data, which was susceptible to case selection and information bias. Additionally, the cases in both the training set and test set were all patients who underwent thyroid ultrasound screening in hospitals. As can be seen from the extracted data table, the proportion of malignant thyroid nodules in each of the included original studies was significantly higher than that in the normal population (8, 30). The study by ChenYF et al., which included thyroid nodules with surgical resection or FNA, had a high overall malignancy rate, and only thyroid nodules with surgical pathological findings were included, which may have resulted in a higher than usual rate of malignant nodules (10). The real-world data would be a greater proportion of benign nodules and a smaller proportion of malignant nodules, which could lead to overfitting and thus compromise accuracy. The surgical rate of small thyroid nodules (<1cm) in Asian countries (such as China and South Korea) is significantly higher than that in Western countries, and there is a higher proportion of overdiagnosis. This difference leads to a different proportion of “surgically confirmed benign nodules” in the training data, and the proportion of malignant cases in the training data is distorted, which affects the model’s judgment threshold for the true malignant risk (40). Additionally, thyroid nodules have a high risk of malignancy, which may have contributed to selection bias in our sample. Since the cytologic results of fine-needle aspiration biopsy may be inconclusive, only the histopathological results of surgical resection can be used as the reference standard for the final diagnosis of nodules.
In our meta-analysis study, 11 studies were identified with case selection bias, as the original authors had selected only surgical patients. According to the 2015 American ATA guidelines, patients with thyroid nodules, if they are highly suspected of malignancy or have been confirmed to be malignant and the Bethesda reporting system indicates category V or VI, even in low-risk Papillary Thyroid Microcarcinoma(PTMC) patients, surgical treatment is not recommended (37). However, these 11 studies included patients who did not undergo FNA after ultrasound evaluation and chose to proceed directly to surgery. The surgical criteria were outlined as follows: benign nodules exceeding 4.0 cm in size and malignant nodules confirmed by preoperative needle biopsy pathology. For cases where ultrasound strongly suggested malignancy but the standards for fine-needle aspiration biopsy (FNAB) weren’t met, FNAB was advised prior to determining subsequent management (22). Therefore, we considered that the patients were highly suspected of malignancy, such as a history of thyroid radiation in childhood, family history of thyroid cancer, extrathyroidal extension, and thyroid ultrasound assessment of TI-RADS classification 4c or 5. We considered that there was a high risk of patient selection bias in these studies, so only surgical patients were selected as the sample. It was rated as high risk in terms of patient selection. Because the original authors evaluated the models using a validation set, and both achieved high diagnostic accuracy, they were still rated as low risk in terms of Applicability Concerns.
4.7 Age of patients
Additionally, this study found that the AI-assisted diagnosis system demonstrated higher diagnostic performance for patients aged 50 years or older. Results such as those of Chen et al. suggest that age is an important risk factor associated with malignancy. The prevalence of thyroid nodules gradually increases with age, which is consistent with previous research results (8, 49). Age-related tissue degeneration, fibrosis, or calcification may change the ultrasonographic appearance of nodules. Older patients often have multiple coexisting lesions that may interfere with image interpretation, and age-related histological changes such as fibrosis and calcification may affect the ultrasonographic features. As a result, the model’s performance fluctuates when applied across regions, affecting its accuracy (50).
4.8 Gender of patients
The detection rate of thyroid nodules in females is higher than that in males, which is a risk factor for the detection of thyroid nodules. The increasing incidence of thyroid cancer among women has been particularly marked since 2000 (39). The results of this study showed that the larger the proportion of women ≥70%, the higher the intelligent assisted diagnosis system’s sensitivity was, at 0.90 (0.87-0.92), and specificity, at 0.85 (0.81-0.89), in the diagnosis of thyroid nodules, (50). The effects of estrogen may be the primary factor contributing to the incidence of thyroid nodules in women. In this study, the majority of the included original studies showed a significantly higher proportion of female patients than male patients, which made the model training more biased toward females. However, in external validation sets such as those in Chenc et al.'s study, the proportion of male patients was significantly higher than that of female patients, which could also affect accuracy (7, 51). The proportion of female patients was positively correlated with AI accuracy, but due to the limited number of original studies with a higher proportion of male patients, gender-specific conclusions still require further validation.
4.9 Reference standards
The gold standard is the reference standard for the accuracy evaluation of diagnostic tests. It must be the most recognized and accurate diagnostic method at present, and its implementation, testing timing, methods, and interpretation standards should be unified. If the gold standard is interpreted by different people, it may lead to differences in the results and thus introduce bias. The different exclusion and inclusion criteria of each of the included original studies were also a potential source of heterogeneity; the study by ChenJH et al. excluded patients with incomplete FNA pathological results or classified as Bethesda I, III, or IV. However, Koh et al. included only surgically confirmed or cytologically confirmed benign (class II) or malignant (class VI) on the Bethesda system. Different reference standards affect the diagnostic accuracy of AI models (9, 13).
4.10 Model architecture
Notably, compared with traditional machine learning methods, deep learning demonstrated higher accuracy in diagnosing thyroid nodules, with a sensitivity of 0.89 (0.87-0.91) and a specificity of 0.85 (0.80-0.89). Traditional machine learning models often use a variety of classifiers in the development process, such as the study of JinZ et al. (11) using random forest (RF), support vector machine (SVM) and extreme gradient boosting (XGBoost) to build ultrasound and ultrasound and clinical combined machine learning models, respectively. Different classifiers have their characteristics. The decision tree (DT)-based methods achieve classification and regression tasks through conditional probability distributions, but DT, as a non-parametric method, is prone to overfitting. However, DT-based ensemble methods (such as random forest, RF; gradient boosting decision tree, GBDT; and XGBoost) significantly improve performance. By randomly selecting feature subsets and samples for parallel training, RF has the advantages of high training efficiency, low generalization error, strong noise resistance, and simple parameter adjustment. GBDT enhanced feature selection ability by gradient boosting (52). XGBoost provides further support for custom loss functions, regularization terms, and handling of missing values, demonstrating greater flexibility. In contrast, linear models have obvious limitations: Support vector machine (SVM) can deal with nonlinear problems through kernel functions, but the selection of kernel functions is time-consuming and difficult to adapt to large-scale data. Although logistic regression (LR) is suitable for multi-classification tasks with small samples, its function is limited. Overall, RF and XGBoost significantly exceed the performance boundaries of traditional DT and linear models with their stability and flexibility, respectively. These traditional models are built on the paradigm of “artificial features + shallow models”, and their performance highly depends on the quality of feature engineering. Features automatically extracted by deep learning (DL) models may contain redundant information (such as artifacts in ultrasound images). Whereas traditional methods (such as logistic regression + manual features) can explicitly incorporate clinical prior knowledge, DL models may ignore these key rules (46). Although the SVM model can achieve nonlinear classification by kernel function, feature extraction completely relies on manual design, which faces problems such as a high threshold of expertise, limited dimension, and feature interaction. The shallow nature of the model structure also limits the expression ability. Linear models can only learn global linear relationships. Although DT can capture nonlinear features, the tree structure based on the greedy algorithm is easy to overfit and ignore the global pattern. These fundamental defects make traditional machine learning perform poorly in processing high-dimensional and complex data such as images and speech, highlighting its inherent limitations in feature learning and deep pattern capture.
Deep learning algorithms consist of structures called deep neural networks, of which convolutional neural networks are a specific type that is widely used in the field of image processing (53). Its multi-threaded task can be used for more complex radiomics image analysis. For example, Zhou et al. ‘s study adopted a CNN architecture and a transfer learning strategy and achieved a high AUC of 0.97 (29). The study by Yao et al., which constructed a multimodal deep learning model involving 6,032 (Papillary Thyroid Cancer, PTC) cases, was used to develop the DeepThy-Net model, which showed high accuracy (54). Lee, J.H., et al. used deep learning to locate and diagnose metastatic lymph nodes in thyroid cancer and achieved high clinical application value (55). These findings collectively demonstrate the excellent diagnostic performance of artificial intelligence, regardless of imaging type, and the emergence of deep learning as a more powerful approach to diagnosing thyroid cancer. Deep learning has an inherent architectural advantage. The deep learning system trained by pre-simulation can automatically obtain relevant ultrasound features through multi-level nonlinear transformation. The hierarchical feature learning capability of deep learning enables the automatic extraction of complex multi-layered imaging features from raw ultrasound data, whereas machine learning relies on predefined features that may miss subtle details (53). For example, a convolutional neural network (CNN) has a layered architecture that mimics a biological vision system: the bottom layer learns local features such as edges and colors, the middle layer combines textures and parts, and the upper layer forms a complete object representation. This end-to-end hierarchical feature learning has the advantages of distributed encoding and incremental refinement, enabling the model to automatically extract more essential feature representations from the data. Deep networks such as the Transformer architecture can model long-distance dependencies in natural language through a multi-layer self-attention mechanism, and achieve the same expressive ability with fewer parameters. The superior performance of deep learning also stems from its ability to learn complex spatial relationships and identify relevant features, especially in analyzing the complex anatomy of thyroid nodules and surrounding structures (56). The data-dependent nature of deep learning algorithms shows that model performance continues to improve with the increase in the amount of training data, while classical machine learning algorithms tend to be stable. Meanwhile, most well-performing deep learning models are generated from baseline architectures, and the diagnostic performance of AI models will be further improved through innovative modifications of training strategies and algorithm architectures (19, 29). As more and more comprehensive features are included in artificial intelligence models, their prediction results will be closer to pathological results. Deep learning models, such as DCNN, can learn features directly from ultrasound images without the limitations of manually designed features, thereby improving the repeatability of diagnosis.
4.11 External validation set
External validation was designed to explore true differences in characteristics between the development and validation cohorts and to evaluate the performance of diagnostic models. All selected studies were conducted strictly according to the independent external validation set, except for VelascoPF2024, which used the same data set for both internal and external validation. However, some of the studies included in the literature utilize public data sets (such as the Stanford dataset) as internal data sets, which is questionable (18, 21). Although this is a practical option, the overall characteristics of these open-source data sets are often difficult to obtain and limited by inadequately collected data. Diagnostic models generally perform well on training datasets; however, there are differences between the performance of the model on the training set and the performance on the validation set. This difference may be reduced when the training and validation datasets are drawn from the same population, and this imperfect model training and validation result in poor model performance on the external validation set. The results of this study showed that the diagnostic performance of the AI model in the non-cross-validation group was higher than that in the cross-validation group, mainly because the sample size of each study in the non-cross-validation group was much larger than that in the cross-validation group, and more patients’ nodules were included in the model, which eventually had better test power and generalization. Due to the lack of large-scale external validation, most AI applications have not yet been applied to clinical decision-making. Therefore, external validation is necessary to reduce the risk of model overfitting, assess the stability of model performance in different populations, and potentially enhance clinicians’ confidence in AI-assisted diagnostic tools.
4.12 Prospective study
In prospective studies, the performance of AI-assisted diagnostic techniques, particularly in terms of sensitivity, specificity, and the area under the curve (AUC), is outstanding. Prospective studies can be enrolled at the onset of symptoms and performed according to prespecified diagnostic criteria and trial algorithms, reducing selection bias. However, retrospective studies may have problems, such as selectivity, including the retrospective selection of patients who have been diagnosed, which can miss some undetected cases and lead to bias. However, only two prospective studies were included in this meta-analysis, accounting for only 7% of the total studies, which may affect the accuracy of the results (19, 30).
In addition, this meta-analysis included various artificial intelligence models, including some commercial computer-aided diagnosis (CAD) and machine learning (ML) models. Although several commercial AI-assisted diagnostic systems have been introduced into clinical practice in recent years, such as the AI-SONIC™ Thyroid model mentioned in the ZhouTH2024 study, which achieved a sensitivity of 95% for identifying malignant thyroid nodules, some manufacturers may adjust the model thresholds to improve sensitivity, which could affect the accuracy of the results in this study. The heterogeneity of imaging systems and operators across different manufacturers may introduce variability during the training process, potentially limiting the interoperability of devices manufactured by other suppliers (30, 57). However, due to differences in algorithm architecture, training dataset quality, and validation methods, significant disparities remain in the diagnostic performance of these systems, and their clinical application value requires further validation and optimization (30).
The bivariate boxplot results of this study suggest significant heterogeneity among the four articles, specifically studies 1, 15, 16, and 28 (Figure 5D). Although this does not affect the reliability of the final pooled effect size, we carefully analyzed these four articles from various aspects, including population characteristics, diagnostic thresholds, reference standards, study design, and sample size, to identify the specific reasons for the outlier studies. We found that studies 1 and 16, BudaM2019 and VelascoPF2024, used cytopathology from fine-needle aspiration (FNA) as the sole reference standard and had lower malignant rates of included nodules compared to other studies, at 10.3% and 11.1%, respectively. Consequently, the specificity for diagnosing thyroid nodules was only 52.0% (6, 21). We analyzed studies 15 and 28, SunC et al.’s DCNN (VGG-F) model, which actually combines multiple features, using SVM for feature extraction, and combines the deep features prioritized by the convolutional neural network (CNN) with hand-crafted features (20). ZhuYC2022 et al.’s ANN (TDUS-Net) model employs color Doppler ultrasound (CDUS) features extracted via deep learning (whole ratio, intranodular ratio, peripheral ratio, and number of vessels) and gray-scale ultrasound (US) features, which differ significantly from those of other models (33). These factors may be a significant source of heterogeneity.
Ultimately, we selected the WeiX2020 study as having the optimal model. The core advantage of the EDLC-TN model lies in its innovative ensemble learning framework and multi-stage design (23). The model first extracts the region of interest (ROI) through precise nodule segmentation, eliminating interference from irrelevant background in ultrasound images, so that subsequent classification can focus on key nodule features (such as edge morphology and internal echoes). In the classification stage, EDLC-TN employs a multi-model ensemble strategy based on DenseNet (58): by training three structurally distinct weak classifiers (based on ROI, mask, and fused features, respectively) and integrating their outputs via voting and averaging methods, the model retains DenseNet’s efficient capture of subtle features through dense connections while addressing the limitations of single algorithms through model diversity. This design enables the model to achieve an accuracy rate of 98.51% (AUC 0.941) on the test set, far surpassing other comparison models. Additionally, the superiority of EDLC-TN is evident in its data and training strategies. The study utilized a multi-center dataset (26,541 images spanning four hospitals and various devices), implemented dynamic learning rate adjustment (Adam optimizer gradually decaying from 0.1), and enforced strict data standardization. In external validation (unseen GE ultrasound device-acquired image data), the model maintained an accuracy rate of 95.76%, demonstrating its cross-device adaptability. The constructed EDLC-TN model is a universal network platform that can be applied to ultrasound images from different medical centers. Whether applied to ultrasound images from hospitals with completely different types of ultrasound devices or compared with the performance of radiologists, the model achieved excellent accuracy, sensitivity, and specificity. This indicates that the EDLC-TN model has the potential to effectively learn from different types of medical images and possesses high generalizability. This will be beneficial for screening programs and lead to more effective referral systems across all medical fields, thereby yielding widespread beneficial clinical and public health outcomes. However, due to the limited number of studies on which the model is based, its generalizability and applicability across different medical centers require further validation through large-scale, multi-center prospective studies.
Additionally, interpretability models remain a significant research challenge. At present, many scholars have doubts about the reliability of AI models. As a “black box” system, these models lack the interpretability of the diagnosis and decision-making process, resulting in a cognitive gap between doctors and models and weakening clinical trust. Interpretable AI is a set of tools and methods that help people understand and interpret the predictions made by machine learning algorithms. This includes an interpretable model and an interpretable interface. This contributes to the accuracy, fairness, and transparency of diagnostic models and understanding of the results of AI-driven decision making (59). The reliability and interpretability of diagnostic AI tools can be improved through collaborative engagement between AI developers and clinical practitioners. In addition, most models only output static binary classification probability results (such as malignant probability) and lack an interactive feedback mechanism with clinicians. This “silent box” feature makes doctors more inclined to use traditional transparent diagnostic methods, which leads to doubts about the actual value of AI and hinders the deep application of big data technology in the medical field (8, 26). In the study of Zheng et al., the model was constructed, and the interpretability analysis was performed through the Shapley Additive Explanations (SHAP) to achieve the transparent effect of the model (28). The study by Yao et al. developed a multimodal generative pre-trained transducer model, ThyGPT, which provides a transparent, interpretable, publicly available, and efficient AI-aided tool for the diagnosis and management of thyroid nodules (42). Model transparency requires further optimization to explain the decision-making process of algorithms, and the accurate conclusions of deep learning models still need to be verified on external test sets with large samples (60–62).
The artificial intelligence system offers a new option for diagnosing thyroid nodules and has significant clinical application value. This meta-analysis presents a comprehensive subgroup analysis of AI models, addressing emerging trends in imaging technology. These methodological advances contribute to a more complete understanding of the role of ultrasonographic AI in the diagnosis of thyroid nodules, although further validation in different populations is still needed. AI diagnostic models (especially deep learning algorithms) can significantly improve the diagnostic accuracy and sensitivity of radiologists at all levels (63), which can effectively improve the diagnostic accuracy of young radiologists and radiologists in primary hospitals, and shorten the growth cycle of radiologists. Urban-rural medical resources are not balanced in China and many countries worldwide. The research and development of artificial intelligence systems can help reduce barriers and provide convenient ways for community hospitals to improve the diagnosis of thyroid cancer. The implementation of ultrasound-based AI in the primary healthcare system could enable the early detection and management of thyroid nodules, particularly thyroid cancer, in remote areas with limited diagnostic resources. AI models can also reduce the time required for image interpretation by doctors, decrease clinical workload, decrease labor intensity, and mitigate the influence of subjective factors on diagnosis. Additionally, they can serve as a reference for clinical experts to enhance the quality of clinical diagnosis of thyroid cancer, particularly for smaller thyroid nodules (64). AI technology can efficiently extract image features, transform the traditional image diagnosis from subjective qualitative analysis to objective quantitative assessment, and realize the whole process of auxiliary diagnosis from nodule detection, benign and malignant judgment to pathological classification and prognosis prediction, which is helpful to promote the transformation from population to individual diagnosis and treatment of thyroid nodules (65). If widely and continuously implemented, it is expected to reduce the relatively invasive FNA biopsy in the diagnosis of thyroid nodules, avoid the discomfort of patients in the process of FNA biopsy and the risk of biopsy-related side effects or complications, realize the early and accurate diagnosis of thyroid nodules, improve the influence on treatment options and prognosis, reduce the risk of delayed diagnosis caused by false negative results, and reduce the anxiety of patients. Reduce healthcare spending and improve patient satisfaction and patient experience. The use of AI models did not lead to a change in the mortality rate despite an increase in the detection rate of thyroid cancer. The medical community also needs to pay attention to the problem of overdiagnosis. For example, it may be more appropriate to use AI models to adopt active surveillance strategies for small papillary carcinomas. In the future, the workflow optimization of AI initial screening and physician confirmation, as well as dynamic monitoring integrated into portable devices, can help doctors dynamically monitor disease progression and optimize the management of high-risk patients, which is expected to further improve the accuracy and accessibility of thyroid nodule management.
In the actual clinical scenario, when physicians need to perform a comprehensive evaluation of patients with thyroid nodules, they may encounter a large amount of information from case data, clinical features, radiomics, and genomics. It is important to study the ability to build a unified model from both macro and micro levels, integrate and analyze information such as pathology, genomics, and clinical data, so as to achieve truly multimodal and cross-body assisted diagnosis. At present, diagnostic models based on thyroid ultrasound image features are gradually favoring hybrid models and multimodal large language models. For example, the weights of the hybrid assisted multi-classification model consisting of Deep Maxout and CNN built by Gulame et al. are automatically learned through transfer learning, which is significantly better than the earlier method of a single model (66). Yao et al. proposed an AI-assisted diagnosis system of a large language model for thyroid nodule risk assessment. The ThyGPT model was trained and verified on large sample data sets, and the LlaMA2-13B model framework was used to segment and store thyroid-related knowledge and cases combined with the Lang-Chain framework. The hybrid architecture of Swin-Transformer and DCNN is used for image analysis, which integrates multi-source information to generate feature markers, and intuitively displays the decision-making basis and the contribution of ultrasound features through human-computer interaction. This model provides a new direction for the development of next-generation CAD systems (67).
Despite the advantages of AI models, enormous challenges remain. First, the use of artificial intelligence models in clinical practice raises ethical and legal issues. Regulators must develop clear guidelines that address issues such as patient privacy, data security, and diagnostic decision-making responsibility to ensure the safe, effective, responsible, and ethical use of AI technologies in clinical Settings. Addressing these challenges is critical for the successful implementation and widespread adoption of auxiliary diagnostic models for thyroid nodules in clinical practice. Second, future AI developments should prioritize adaptability, algorithmic rhythm transparency, and interpretability to facilitate professional acceptance and integration of AI-assisted diagnostic technologies. At the same time, unresolved issues regarding data security protocols and ambiguous accountability frameworks in smart healthcare technologies continue to impede clinical adoption. Third, at present, AI is more suitable for use as an adjunct rather than a stand-alone diagnostic method, and its results must be reviewed by certified physicians. Further research and method standardization are needed to address these discrepancies in the future, to integrate AI software seamlessly into physicians’ workflows, and external validation studies with large samples should be conducted (68). Fourth, clinical ultrasound is a highly dynamic diagnosis and treatment process, and doctors rely on real-time images for comprehensive judgment. Dynamic data processing will face many technical challenges. The amount of data to be processed by real-time video analysis is hundreds of times that of static images, which requires the development of new neural network architectures (such as 3D CNN+Transformer). Instead of the current mainstream 2D CNN model, the inference delay needs to be controlled to ensure clinical usability, which puts forward extremely high requirements for algorithm optimization. A lightweight model (such as deformable convolution) specifically for the characteristics of ultrasound video is developed, which integrates color Doppler, pulse wave, and other multimodal timing data, and the chip inside the ultrasound probe realizes end-to-end real-time processing. In the future, dynamic ultrasound database standards (such as probe motion trajectory and time synchronization markers) will be established. Evaluation frameworks such as real-time human-computer interaction tests will be developed through databases to simulate clinical environments. Progressive verification processes will be designed from static frames to video clips, and then stepwise verification of complete clinical scenarios will be formed through real-time streaming. Finally, international multicenter data sets should be established in the future to ensure that the samples of each subgroup are balanced, include data from different ultrasound devices and different acquisition protocols, and adopt unified image annotation standards and pathological confirmation procedures to construct diagnostic models for different subtypes of thyroid nodules. The model needs to create a well-defined set of features that can be efficiently captured by human visual assessment and computational analysis (69). The model can incorporate the characteristics of cervical lymph nodes and also provide value for surgical methods used in treating thyroid cancer patients, such as lymph node dissection. AI models offer promising prospects for clinical applications, and factors such as human selection bias and reliance on human labor can be improved.
This meta-analysis should acknowledge several important limitations. First, the heterogeneity among the included studies is significant. In addition to the various factors discussed in the analysis section, there may be some unrecognized factors that have been overlooked or that could not be further analyzed due to the limited number of original studies. Second, over half of the studies employed retrospective designs, which may introduce potential selection and bias. Since researchers developed AI models, the training phase required the construction of models with accurate labels based on reference standards, a step that necessitated a retrospective design, whereas the predictive phase, which involved a validation set (or test set), did not. Third, the generalizability of these findings may be limited by the fact that most studies were conducted in Asian regions, so our findings may have limited generalizability to populations in other regions. Fourth, since some basic characteristics of thyroid nodules, such as nodule echogenicity, internal structure, and calcification, were not included, it is not possible to further explore the impact of these basic characteristics on diagnostic performance. Finally, we only included studies written in English, which may make this meta-analysis more susceptible to publication bias.
5 Conclusion
This meta-analysis examined the clinical value of an artificial intelligence-assisted diagnosis system based on ultrasound images in the diagnosis of benign and malignant thyroid nodules, and the results demonstrated high clinical potential. Among them, the EDLC-TN model shows higher diagnostic accuracy and clinical effectiveness in the diagnosis of thyroid nodules. For thyroid nodules in female patients with an average diameter of <20mm and an age of ≥50 years, artificial intelligence-assisted diagnostic models are more effective, especially deep learning models.
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.
Author contributions
JFZ: Writing – original draft, Investigation, Resources, Visualization, Software, Data curation, Project administration, Conceptualization, Validation, Writing – review & editing, Methodology, Formal analysis. JZ: Conceptualization, Investigation, Validation, Formal analysis, Methodology, Data curation, Writing – review & editing, Writing – original draft. SZ: Data curation, Conceptualization, Writing – original draft. LN: Writing – original draft, Investigation, Conceptualization, Data curation. CZ: Investigation, Conceptualization, Writing – original draft, Data curation. JH: Project administration, Supervision, Writing – review & editing, Investigation, Funding acquisition, Resources.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. The program of the medicine project in the Huzhou Science and Technology Bureau.
Acknowledgments
The authors would like to thank any researcher who contributed to this study.
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.
Generative AI statement
The author(s) declare that no Generative AI was 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2025.1614603/full#supplementary-material
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:209–49. doi: 10.3322/caac.21660
2. Kong W, Yue X, Ren J, and Tao X. A comparative analysis of diffusion-weighted imaging and ultrasound in thyroid nodules. BMC Med Imaging. (2019) 19:92. doi: 10.1186/s12880-019-0381-x
3. Mai W, Zhou M, Li J, Yi W, Li S, Hu Y, et al. The value of the Demetics ultrasound-assisted diagnosis system in the differential diagnosis of benign from Malignant thyroid nodules and analysis of the influencing factors. Eur Radiol. (2021) 31:7936–44. doi: 10.1007/s00330-021-07884-z
4. Eidt LB, Nunes de Oliveira C, Lagos YBB, Solera GLM, Izquierdo R, Meyer ELS, et al. A prospective comparison of ACR-TIRADS and EU-TIRADS in thyroid nodule assessment for FNA-US. Clin Endocrinol (Oxf). (2023) 98:415–25. doi: 10.1111/cen.14799
5. Cappelli C, Pirola I, Agosti B, Tironi A, Gandossi E, Incardona P, et al. Complications after fine-needle aspiration cytology: a retrospective study of 7449 consecutive thyroid nodules. Br J Oral Maxillofac Surg. (2017) 55:266–9. doi: 10.1016/j.bjoms.2016.11.321
6. Buda M, Wildman-Tobriner B, Hoang JK, Thayer D, Tessler FN, Middleton WD, et al. Management of thyroid nodules seen on US images: deep learning may match performance of radiologists. Radiology. (2019) 292:695–701. doi: 10.1148/radiol.2019181343
7. Chen C, Jiang Y, Yao J, Lai M, Liu Y, Jiang X, et al. Deep learning to assist composition classification and thyroid solid nodule diagnosis: a multicenter diagnostic study. Eur Radiol. (2024) 34:2323–33. doi: 10.1007/s00330-023-10269-z
8. Chen D, Hu J, Zhu M, Tang N, Yang Y, and Feng Y. Diagnosis of thyroid nodules for ultrasonographic characteristics indicative of Malignancy using random forest. BioData Min. (2020) 13:14. doi: 10.1186/s13040-020-00223-w
9. Chen JH, Zhang YQ, Zhu TT, Zhang Q, Zhao AX, and Huang Y. Applying machine-learning models to differentiate benign and Malignant thyroid nodules classified as C-TIRADS 4 based on 2D-ultrasound combined with five contrast-enhanced ultrasound key frames. Front Endocrinol (Lausanne). (2024) 15:1299686. doi: 10.3389/fendo.2024.1299686
10. Chen Y, Gao Z, He Y, Mai W, Li J, Zhou M, et al. An artificial intelligence model based on ACR TI-RADS characteristics for US diagnosis of thyroid nodules. Radiology. (2022) 303:613–9. doi: 10.1148/radiol.211455
11. Jin Z, Pei S, Ouyang L, Zhang L, Mo X, Chen Q, et al. Thy-Wise: An interpretable machine learning model for the evaluation of thyroid nodules. Int J Cancer. (2022) 151:2229–43. doi: 10.1002/ijc.34248
12. Kim YJ, Choi Y, Hur SJ, Park KS, Kim HJ, Seo M, et al. Deep convolutional neural network for classification of thyroid nodules on ultrasound: Comparison of the diagnostic performance with that of radiologists. Eur J Radiol. (2022) 152:110335. doi: 10.1016/j.ejrad.2022.110335
13. Koh J, Lee E, Han K, Kim EK, Son EJ, Sohn YM, et al. Diagnosis of thyroid nodules on ultrasonography by a deep convolutional neural network. Sci Rep. (2020) 10:15245. doi: 10.1038/s41598-020-72270-6
14. Ko SY, Lee JH, Yoon JH, Na H, Hong E, Han K, et al. Deep convolutional neural network for the diagnosis of thyroid nodules on ultrasound. Head Neck. (2019) 41:885–91. doi: 10.1002/hed.25415
15. Lai M, Feng B, Yao J, Wang Y, Pan Q, Chen Y, et al. Value of artificial intelligence in improving the accuracy of diagnosing TI-RADS category 4 nodules. Ultrasound Med Biol. (2023) 49:2413–21. doi: 10.1016/j.ultrasmedbio.2023.08.008
16. Li L, Deng H, Chen W, Wu L, Li Y, Wang J, et al. Comparison of the diagnostic effectiveness of ultrasound imaging coupled with three mathematical models for discriminating thyroid nodules. Acta Radiol. (2024) 65:441–8. doi: 10.1177/02841851231221912
17. Li X, Zhang S, Zhang Q, Wei X, Pan Y, Zhao J, et al. Diagnosis of thyroid cancer using deep convolutional neural network models applied to sonographic images: a retrospective, multicohort, diagnostic study. Lancet Oncol. (2019) 20:193–201. doi: 10.1016/S1470-2045(18)30762-9
18. Nair G, Vedula A, Johnson ET, Thomas J, Patel R, Cheng J, et al. Combining image similarity and predictive artificial intelligence models to decrease subjectivity in thyroid nodule diagnosis and improve Malignancy prediction. Endocr Pract. (2024) 30:1031–7. doi: 10.1016/j.eprac.2024.08.001
19. Park VY, Han K, Seong YK, Park MH, Kim EK, Moon HJ, et al. Diagnosis of thyroid nodules: performance of a deep learning convolutional neural network model vs. Radiologists. Sci Rep. (2019) 9:17843. doi: 10.1038/s41598-019-54434-1
20. Sun C, Zhang Y, Chang Q, Liu T, Zhang S, Wang X, et al. Evaluation of a deep learning-based computer-aided diagnosis system for distinguishing benign from Malignant thyroid nodules in ultrasound images. Med Phys. (2020) 47:3952–60. doi: 10.1002/mp.14301
21. Fernandez Velasco P, Perez Lopez P, Torres Torres B, Delgado E, de Luis D, and Diaz Soto G. Clinical evaluation of an artificial intelligence-based decision support system for the diagnosis and american college of radiology thyroid imaging reporting and data system classification of thyroid nodules. Thyroid. (2024) 34:510–8. doi: 10.1089/thy.2023.0603
22. Wang L, Yang S, Yang S, Zhao C, Tian G, Gao Y, et al. Automatic thyroid nodule recognition and diagnosis in ultrasound imaging with the YOLOv2 neural network. World J Surg Oncol. (2019) 17:12. doi: 10.1186/s12957-019-1558-z
23. Wei X, Gao M, Yu R, Liu Z, Gu Q, Liu X, et al. Ensemble deep learning model for multicenter classification of thyroid nodules on ultrasound images. Med Sci Monit. (2020) 26:e926096. doi: 10.12659/MSM.926096
24. Wu GG, Lv WZ, Yin R, Xu JW, Yan YJ, Chen RX, et al. Deep learning based on ACR TI-RADS can improve the differential diagnosis of thyroid nodules. Front Oncol. (2021) 11:575166. doi: 10.3389/fonc.2021.575166
25. Yao J, Zhang Y, Shen J, Lei Z, Xiong J, Feng B, et al. AI diagnosis of Bethesda category IV thyroid nodules. iScience. (2023) 26:108114. doi: 10.1016/j.isci.2023.108114
26. Zhang B, Tian J, Pei S, Chen Y, He X, Dong Y, et al. Machine learning-assisted system for thyroid nodule diagnosis. Thyroid. (2019) 29:858–67. doi: 10.1089/thy.2018.0380
27. Zhao CK, Ren TT, Yin YF, Shi H, Wang HX, Zhou BY, et al. A comparative analysis of two machine learning-based diagnostic patterns with thyroid imaging reporting and data system for thyroid nodules: diagnostic performance and unnecessary biopsy rate. Thyroid. (2021) 31:470–81. doi: 10.1089/thy.2020.0305
28. Zheng Y, Zhang Y, Lu K, Wang J, Li L, Xu D, et al. Diagnostic value of an interpretable machine learning model based on clinical ultrasound features for follicular thyroid carcinoma. Quant Imaging Med Surg. (2024) 14:6311–24. doi: 10.21037/qims-24-601
29. Zhou H, Jin Y, Dai L, Zhang M, Qiu Y, Wang K, et al. Differential diagnosis of benign and Malignant thyroid nodules using deep learning radiomics of thyroid ultrasound images. Eur J Radiol. (2020) 127:108992. doi: 10.1016/j.ejrad.2020.108992
30. Zhou T, Xu L, Shi J, Zhang Y, Lin X, Wang Y, et al. US of thyroid nodules: can AI-assisted diagnostic system compete with fine needle aspiration? Eur Radiol. (2024) 34:1324–33. doi: 10.1007/s00330-023-10132-1
31. Zhu J, Zhang S, Yu R, Liu Z, Gao H, Yue B, et al. An efficient deep convolutional neural network model for visual localization and automatic diagnosis of thyroid nodules on ultrasound images. Quant Imaging Med Surg. (2021) 11:1368–80. doi: 10.21037/qims-20-538
32. Zhu LC, Ye YL, Luo WH, Su M, Wei HP, Zhang XB, et al. A model to discriminate Malignant from benign thyroid nodules using artificial neural network. PloS One. (2013) 8:e82211. doi: 10.1371/journal.pone.0082211
33. Zhu YC, Du H, Jiang Q, Zhang T, Huang XJ, Zhang Y, et al. Machine learning assisted doppler features for enhancing thyroid cancer diagnosis: A multi-cohort study. J Ultrasound Med. (2022) 41:1961–74. doi: 10.1002/jum.15873
34. Mauri G, Nicosia L, Della Vigna P, Varano GM, Maiettini D, Bonomo G, et al. Percutaneous laser ablation for benign and Malignant thyroid diseases. Ultrasonography. (2019) 38:25–36. doi: 10.14366/usg.18034
35. Pacella CM. Image-guided thermal ablation of benign thyroid nodules. J Ultrasound. (2017) 20:347–9. doi: 10.1007/s40477-017-0269-y
36. Xu W, Jia X, Mei Z, Gu X, Lu Y, Fu CC, et al. Generalizability and diagnostic performance of AI models for thyroid US. Radiology. (2023) 307:e221157. doi: 10.1148/radiol.221157
37. 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
38. Papini E, Guglielmi R, Bianchini A, Crescenzi A, Taccogna S, Nardi F, et al. Risk of Malignancy in nonpalpable thyroid nodules: predictive value of ultrasound and color-Doppler features. J Clin Endocrinol Metab. (2002) 87:1941–6. doi: 10.1210/jcem.87.5.8504
39. Miranda-Filho A, Lortet-Tieulent J, Bray F, Cao B, Franceschi S, Vaccarella S, et al. Thyroid cancer incidence trends by histology in 25 countries: a population-based study. Lancet Diabetes Endocrinol. (2021) 9:225–34. doi: 10.1016/S2213-8587(21)00027-9
40. Vuong HG, Ngo HTT, Bychkov A, Jung CK, Vu TH, Lu KB, et al. Differences in surgical resection rate and risk of Malignancy in thyroid cytopathology practice between Western and Asian countries: A systematic review and meta-analysis. Cancer Cytopathol. (2020) 128:238–49. doi: 10.1002/cncy.22228
41. Han Y, Wu JQ, Hou XJ, Sun JW, Piao ZY, Teng F, et al. Strain imaging in the evaluation of thyroid nodules: the associated factors leading to misdiagnosis. Ultrasound Med Biol. (2021) 47:3372–83. doi: 10.1016/j.ultrasmedbio.2021.08.016
42. Yao J, Wang Y, Lei Z, Wang K, Feng N, Dong F, et al. Multimodal GPT model for assisting thyroid nodule diagnosis and management. NPJ Digit Med. (2025) 8:245. doi: 10.1038/s41746-025-01652-9
43. Mei X, Liu Z, Robson PM, Marinelli B, Huang M, Doshi A, et al. RadImageNet: an open radiologic deep learning research dataset for effective transfer learning. Radiol Artif Intell. (2022) 4:e210315. doi: 10.1148/ryai.210315
44. Chang Y, Paul AK, Kim N, Baek JH, Choi YJ, Ha EJ, et al. Computer-aided diagnosis for classifying benign versus Malignant thyroid nodules based on ultrasound images: A comparison with radiologist-based assessments. Med Phys. (2016) 43:554. doi: 10.1118/1.4939060
45. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational radiomics system to decode the radiographic phenotype. Cancer Res. (2017) 77:e104–e7. doi: 10.1158/0008-5472.CAN-17-0339
46. Suzuki K. Overview of deep learning in medical imaging. Radiol Phys Technol. (2017) 10:257–73. doi: 10.1007/s12194-017-0406-5
47. Gillies R, Kinahan P, and Hricak H. Radiomics: images are more than pictures, they are data. Radiology. (2016) 278:563–77. doi: 10.1148/radiol.2015151169
48. Ha EJ, Lee JH, Lee DH, Moon J, Lee H, Kim YN, et al. Artificial intelligence model assisting thyroid nodule diagnosis and management: A multicenter diagnostic study. J Clin Endocrinol Metab. (2024) 109:527–35. doi: 10.1210/clinem/dgad503
49. Mu C, Ming X, Tian Y, Liu Y, Yao M, Ni Y, et al. Mapping global epidemiology of thyroid nodules among general population: A systematic review and meta-analysis. Front Oncol. (2022) 12:1029926. doi: 10.3389/fonc.2022.1029926
50. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, and Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2018) 68:394–424. doi: 10.3322/caac.21492
51. Xu S, Chen G, Peng W, Renko K, and Derwahl M. Oestrogen action on thyroid progenitor cells: relevant for the pathogenesis of thyroid nodules? J Endocrinol. (2013) 218:125–33. doi: 10.1530/JOE-13-0029
52. Zhang Z and Jung C. GBDT-MO: gradient-boosted decision trees for multiple outputs. IEEE Trans Neural Netw Learn Syst. (2021) 32:3156–67. doi: 10.1109/TNNLS.2020.3009776
53. Shen D, Wu G, and Suk HI. Deep learning in medical image analysis. Annu Rev BioMed Eng. (2017) 19:221–48. doi: 10.1146/annurev-bioeng-071516-044442
54. Yao J, Lei Z, Yue W, Feng B, Li W, Ou D, et al. DeepThy-net: A multimodal deep learning method for predicting cervical lymph node metastasis in papillary thyroid cancer. Advanced Intelligent Syst. (2022) 4:2200100. doi: 10.1002/aisy.202200100
55. Lee JH, Baek JH, Kim JH, Shim WH, Chung SR, Choi YJ, et al. Deep learning-based computer-aided diagnosis system for localization and diagnosis of metastatic lymph nodes on ultrasound: A pilot study. Thyroid. (2018) 28:1332–8. doi: 10.1089/thy.2018.0082
56. Svecic A, Roberge D, and Kadoury S. Prediction of inter-fractional radiotherapy dose plans with domain translation in spatiotemporal embeddings. Med Image Anal. (2020) 64:101728. doi: 10.1016/j.media.2020.101728
57. Kim PH, Suh CH, Baek JH, Chung SR, Choi YJ, and Lee JH. Unnecessary thyroid nodule biopsy rates under four ultrasound risk stratification systems: a systematic review and meta-analysis. Eur Radiol. (2021) 31:2877–85. doi: 10.1007/s00330-020-07384-6
58. Huang. G, Liu. Z, Maaten. LVD, and Weinberger. KQ. Densely connected convolutional networks. In: Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Honolulu, HI, USA: IEEE (2017), 2261–9. doi: 10.1109/CVPR.2017.243
59. Chen H, Gomez C, Huang CM, and Unberath M. Explainable medical imaging AI needs human-centered design: guidelines and evidence from a systematic review. NPJ Digit Med. (2022) 5:156. doi: 10.1038/s41746-022-00699-2
60. Grimes DA and Schulz KF. Bias and causal associations in observational research. Lancet. (2002) 359:248–52. doi: 10.1016/S0140-6736(02)07451-2
61. Liu Y, Feng Y, Qian L, Wang Z, and Hu X. Deep learning diagnostic performance and visual insights in differentiating benign and Malignant thyroid nodules on ultrasound images. Exp Biol Med (Maywood). (2023) 248:2538–46. doi: 10.1177/15353702231220664
62. Yu J, Deng Y, Liu T, Zhou J, Jia X, Xiao T, et al. Lymph node metastasis prediction of papillary thyroid carcinoma based on transfer learning radiomics. Nat Commun. (2020) 11:4807. doi: 10.1038/s41467-020-18497-3
63. Zhang MB, Meng ZL, Mao Y, Jiang X, Xu N, Xu QH, et al. Cervical lymph node metastasis prediction from papillary thyroid carcinoma US videos: a prospective multicenter study. BMC Med. (2024) 22:153. doi: 10.1186/s12916-024-03367-2
64. Zhou LQ, Zeng SE, Xu JW, Lv WZ, Mei D, Tu JJ, et al. Deep learning predicts cervical lymph node metastasis in clinically node-negative papillary thyroid carcinoma. Insights Imaging. (2023) 14:222. doi: 10.1186/s13244-023-01550-2
65. Li LR, Du B, Liu HQ, and Chen C. Artificial intelligence for personalized medicine in thyroid cancer: current status and future perspectives. Front Oncol. (2020) 10:604051. doi: 10.3389/fonc.2020.604051
66. Gulame MB and Dixit VV. Hybrid deep learning assisted multi classification: Grading of Malignant thyroid nodules. Int J Numer Method BioMed Eng. (2024) 40:e3824. doi: 10.1002/cnm.3824
67. Jincao Y, Yunpeng W, Zhikai L, Kai W, Xiaoxian L, Jianhua Z, et al. AI-generated content enhanced computer-aided diagnosis model for thyroid nodules :A chatGPT-style assistant. arXiv. (2024). doi: 10.48550/arXiv.2402.02401
68. Tessler FN and Thomas J. Artificial intelligence for evaluation of thyroid nodules: A primer. Thyroid. (2023) 33:150–8. doi: 10.1089/thy.2022.0560
Keywords: artificial intelligence, diagnosis, thyroid nodule, thyroid neoplasms, ultrasound
Citation: Zhan J, Zhang J, Zhu S, Ni L, Zhang C and Hu J (2025) Diagnostic performance of ultrasound characteristics-based artificial intelligence models for thyroid nodules: a systematic review and meta-analysis. Front. Oncol. 15:1614603. doi: 10.3389/fonc.2025.1614603
Received: 22 April 2025; Accepted: 11 August 2025;
Published: 03 September 2025.
Edited by:
Marta Zerunian, Sapienza University of Rome, ItalyReviewed by:
Jincao Yao, University of Chinese Academy of Sciences, ChinaGe Zhang, Wuhan University of Science and Technology, China
Copyright © 2025 Zhan, Zhang, Zhu, Ni, Zhang and Hu. 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: Jia Hu, NTA3MjBAempodS5lZHUuY24=