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

ORIGINAL RESEARCH article

Front. Med., 08 July 2025

Sec. Ophthalmology

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1623286

This article is part of the Research TopicAdvances in Ocular Autoimmune DiseasesView all 6 articles

Multidimensional predictive model for assessing clinical activity in thyroid eye disease

Yang Li,,&#x;Yang Li1,2,3Guang-Hong Zhang,,&#x;Guang-Hong Zhang1,3,4Man Tian,,Man Tian1,2,3Chuan Hua,,,Chuan Hua1,2,3,5Jian-Ping Zhai,,Jian-Ping Zhai1,3,6Yan-Qiong He,,
Yan-Qiong He1,3,6*Xin-He Zuo,,,
Xin-He Zuo1,2,3,5*
  • 1Thyroid Center of Hubei Provincial Hospital of Traditional Chinese Medicine, Wuhan, China
  • 2Hubei Shizhen Laboratory, Wuhan, China
  • 3Affiliated Hospital of Hubei University of Chinese Medicine, Wuhan, China
  • 4Department of Ophthalmology, Hubei Provincial Hospital of Traditional Chinese Medicine, Wuhan, China
  • 5Hubei Key Laboratory of The Theory and Application Research of Liver and Kidney in Traditional Chinese Medicine, Department of Hepatology, Hubei Provincial Hospital of Traditional Chinese Medicine, Wuhan, China
  • 6Department of Nuclear Medicine, Hubei Provincial Hospital of Traditional Chinese Medicine, Wuhan, China

Objective: Thyroid eye disease (TED) is an autoimmune disorder with complex inflammatory activity that remains challenging to assess accurately. Current method, mainly the Clinical Activity Score (CAS), exhibits limitations in objectivity and comprehensiveness. This study aimed to develop a multidimensional predictive model integrating clinical parameters, SPECT/CT imaging data, and serum biomarkers, to improve TED activity evaluation.

Methods: This retrospective research included 36 TED patients (72 eyes) diagnosed by EUGOGO criteria who underwent SPECT/CT examination. The Clinical Activity Score (CAS) was used to evaluate inflammatory activity. Variables with significant associations with CAS-defined activity were identified using univariate analysis, and Bayesian shrinkage regression (BSR) and the least absolute shrinkage and selection operator (LASSO) were utilized for variable selection in the primary cohort. Predictive models were constructed and evaluated using receiver operating characteristic (ROC) curves (internally validated via five-fold cross-validation), decision curve analysis (DCA), and calibration curves.

Results: Five predictive models were constructed. The comprehensive Model 4, combining clinical, imaging [EX, maximal SPECT/CT uptake ratio (URmax)], and serum biomarkers (TRAb, RBC), achieved superior diagnostic accuracy (AUC: 91.18%; sensitivity: 0.91; specificity: 0.86). Model 5, retaining variables significant in univariate and multivariate analyses, demonstrated robust performance (AUC: 85.97%) with superior stability during cross-validation (ROC mean: 0.8417). Key predictors included male sex (OR = 11.74), TRAb levels, EX, URmax, and RBC count. SPECT/CT-derived URmax correlated strongly with disease activity, while serum biomarkers complemented imaging limitations.

Conclusion: Multidimensional integration of clinical, imaging, and biomarker data significantly enhances TED activity evaluation compared to single-modality approaches. The multidimensional model offers superior diagnostic accuracy, addressing the limitations of conventional methods. These findings advocate for a holistic approach in TED management.

Introduction

Thyroid eye disease (TED), which is also referred to as Graves’ orbitopathy (GO) or thyroid-associated ophthalmopathy (TAO), manifests as an autoimmune disorder involving inflammation and fibrosis of the orbital tissues (1, 2). As a frequent extra-thyroidal manifestation of Graves’ disease, TED exhibits a symptom spectrum spanning mild conjunctival inflammation to severe orbital complications including proptosis, diplopia, and optic nerve impairment (3). The prevalence of TED in Graves’ disease patients ranges from 25 to 50%, and is even higher in uncontrolled hyperthyroidism populations (4).

Accurate assessment of TED activity is critical for guiding treatment decisions. The Clinical Activity Score (CAS), a 7-point scale based on clinical signs (e.g., eyelid swelling, conjunctival redness), is widely used to evaluate inflammatory activity (5). However, CAS exhibits limitations in reproducibility and subclinical inflammation detection, and its results are highly dependent on the doctor’s sensitivity. A cohort study in Asian populations demonstrated that 38% with moderate-to-severe and 10% with vision-threatening TED had CAS <3, suggesting potential underestimation of disease activity (6). These discrepancies highlight the need for objective methods.

Single-photon emission computed tomography/computed tomography (SPECT/CT) has improved the quantification of orbital inflammation (7). By measuring the uptake ratio (UR) of radiotracers 99mTc-DTPA in extraocular muscles relative to the occipital lobe, SPECT/CT provides semi-quantitative data correlating with disease severity (8). A multimodal radiomic model incorporating 99mTc-DTPA SPECT/CT imaging demonstrated statistically significant precision in assessing TED clinical activity (9). Despite these advances, SPECT/CT has its limitations: it primarily reflects localized muscular metabolic activity and fails to capture systemic immune dysregulation, which plays a central role in TED pathogenesis.

Emerging evidence suggests that systemic biomarkers could complement imaging modalities (1012). For instance, thyrotrophin receptor antibody (TRAb) levels correlate with TED progression and recurrence (4). Integrating clinical biomarkers with imaging data may enhance diagnostic precision. This study aims to construct a multidimensional model combining clinical data, SPECT/CT data, and serum biomarkers to improve the accuracy of TED activity evaluation.

Methods

Study design and subjects

This retrospective study included 36 TED patients (72 eyes), diagnosed by EUGOGO criteria (13), that were undertook SPECT/CT examination at Hubei Provincial Hospital of Traditional Chinese Medicine between January 2024 and October 2024. The inclusion criteria for TED patient: (1) age >18 years, (2) sufficient image quality, and (3) without systemic glucocorticoid treatment for at least 3 months prior to the study. The exclusion criteria: (1) non-thyroid systemic comorbidities, (2) prior ocular surgery/trauma, and (3) non-TED ophthalmic conditions. Ethical clearance (Approval No. HBZY2024-C04-01) was obtained from Hubei Provincial Hospital of Traditional Chinese Medicine’s Ethics Committee, following the Declaration of Helsinki guidelines. The comprehensive workflow of this research was shown in Figure 1.

Figure 1
Flowchart detailing a retrospective cohort analysis of TED patients using SPECT/CT imaging. Data from 39 records were collected, with three excluded. Final inclusion had 36 subjects, classified into inactive (29 eyes) and active (43 eyes). Lasso regression and BSR screened variables, leading to multidimensional predictive models. Model evaluation used ROC curve analysis, decision curves, and calibration curve.

Figure 1. The comprehensive workflow of this research.

Clinical and systemic measurements

The clinical information, including sex, age, body mass index (BMI), systolic blood pressure (SBP), diastolic blood pressure (DBP), was recorded. Clinical activity assessment was conducted on a per-eye basis using the Clinical Activity Score (CAS) criteria (14). Thyroid function assessment was performed using electrochemiluminescent immunoassays. Concentrations of serum free triiodothyronine (FT3), free thyroxine (FT4), thyroid-stimulating hormone (TSH), thyroglobulin antibodies (TGAb), and thyroid peroxidase antibodies (TPOAb) were analyzed on a SIEMENS Atellica IM analyzer (Siemens Healthineers, Germany) using manufacturer-supplied kits (Ref. Nos. 10995584, 10995588, 10995704, 11201761, and 10995467, respectively). Thyroid-stimulating hormone receptor antibodies (TRAb) were measured on a Roche cobase 411 analyzer (Roche Diagnostics, Germany) using the Immulite TRAb kit (Ref. No. 08496609190). All assays followed standardized protocols per kit instructions: 50 μL serum samples were incubated with paramagnetic particles coated with capture antibodies/antigens, followed by chemiluminescent substrate addition. Signal detection was proportional to analyte concentration, calibrated against traceable reference materials. Serum biomarkers examination included a complete blood count, lipid profile, and the calculation of neutrophil-to-lymphocyte ratio (NLR), platelet-to-lymphocyte ratio (PLR), and systemic immune-inflammatory index (SII) (platelet count × neutrophil count/lymphocyte count).

Ophthalmic examination by 99mTc-DTPA SPECT/CT scan

99mTc-DTPA at a dose of 370–555 MBq was administered through standardized antecubital venous access. Subsequent SPECT/CT scans were obtained 30–45 min later by a dual-head gamma camera system (Discovery NM/CT 670 Pro, GE Healthcare, USA) as described in previous report (15). The uptake ratio (UR) of extraocular muscles of each eye included was calculated by two senior nuclear medicine physicians blinded to the grouping. The degree of the exophthalmos was evaluated on the basis of SPECT/CT. The straight line joining the anterior tips of the zygomatic bones was the line of reference. A distance greater than 20 mm from the reference line to the anterior border of the eye was considered diagnostic of exophthalmos (EX) (8).

Statistical analysis

Spearman’s correlation was used to assess associations between variables and CAS-defined activity. All variables were subjected to univariate analysis, and those with statistically significant associations (p < 0.05) were included in the subsequent multivariate model. To select variables in the primary cohort and avoid over-fitting or under-fitting, two multivariate regression methods—Bayesian shrinkage regression (BSR) and least absolute shrinkage and selection operator (LASSO)—were employed. For BSR, variable selection was guided by the Bayesian information criterion (BIC). Predictive models were constructed and evaluated using receiver operating characteristic (ROC) curves, assessed through internal five-fold cross-validation for stability, followed by decision curve analysis (DCA), and calibration curves. All the analyses were conducted using software packages R (version 4.2.2, The R Foundation, http://www.R-project.org) and Free Statistic software (version 2.1, Beijing, China).

Results

Basic demographic data and clinical characteristics

A total of 29 inactive eyes and 43 active eyes of TED patients were included in this study. The demographic data and clinical characteristics were depicted in Table 1. Significant differences between the two groups were observed in sex, SBP, EX, URmax, TRAb, RBC, and LDL-C (all p < 0.05). Representative unilateral radiographic manifestations of active and inactive TED defined by CAS were shown in Figure 2. The active TED exhibited increased radionuclide uptake in orbital tissues, including the extraocular muscles, lacrimal glands, and paranasal sinuses.

Table 1
www.frontiersin.org

Table 1. Comparison of basic demographic data and clinical characteristics between active and inactive TED patients defined by CAS.

Figure 2
This figure demonstrates comparative SPECT/CT imaging (CT, SPECT, and fused images) of active and inactive TED. The top row shows inactive states with less defined features, while the bottom row shows active states with more distinct features, especially in the SPECT/CT image where there is visible yellow highlighting.

Figure 2. Representative unilateral radiographic manifestations of active and inactive TED defined by CAS.

Univariate analysis of variables associated with CAS-defined activity

Univariate analysis revealed that sex, SBP, EX, URmax, TRAb, and RBC were significantly associated with CAS-defined disease activity (Table 2). Males showed a higher odds ratio (OR) for clinical activity compared to females (OR = 11.74, p = 0.002). SBP (OR = 1.05, p = 0.015), EX (OR = 1.49, p = 0.002), URmax (OR = 1.55, p = 0.028), TRAb (OR = 1.04, p = 0.038), and RBC (OR = 3.28, p = 0.013) were positively correlated with clinical activity. In contrast, no significant associations were found for age, BMI, DBP, FT3, FT4, TSH, CRP, WBC, PLT, NEU, LYM, NLR, PLR, SII, TGAb, and TPOAb (all p > 0.05) (Table 2).

Table 2
www.frontiersin.org

Table 2. Univariate analysis of variables associated with CAS-defined activity.

Variables selection using the BSR and the LASSO

The combined application of BSR and LASSO identified sex, age, DBP, EX, URmax, FT3, FT4, TSH, TRAb, RBC, and PLR as significant predictors (Figures 3AC). Five predictive models were constructed:

(i) Model 1: exclusively incorporated imaging data URmax and EX.

(ii) Model 2: integrated clinical characteristics (sex, age, and DBP) with imaging data (URmax and EX).

(iii) Model 3: integrated clinical characteristics (sex, age, and DBP) with blood markers (FT3, FT4, TSH, TRAb, RBC, and PLR).

(iv) Model 4: comprehensively combined clinical characteristic, imaging data, and blood markers.

(v) Model 5: selected variables with statistical significance from both univariate and regression analyses (sex, EX, URmax, TRAb, and RBC) (Figure 3D).

Figure 3
A four-panel figure displays different statistical visualizations. Panel A shows a line plot of coefficients versus log lambda for various variables. Panel B presents a plot of binomial deviance against log lambda, highlighting a curve with red dots. Panel C is a bar chart illustrating the distribution of different variables on the x-axis with corresponding values on the y-axis. Panel D exhibits a network diagram with interconnected nodes labeled with variables, representing relationships identified through Lasso and univariate analysis.

Figure 3. Variables for predicting model screening results. (A,B) Variables screening based on the LASSO regression. (C) Variables screening based on the BSR. (D) Variables selected from both univariate and regression analyses with statistical significance.

ROC curve analysis for multidimensional predictive models

ROC analysis demonstrated varying diagnostic performance across models (Table 3 and Figure 4A). Model 4 achieved the highest discriminative power, with an AUC of 91.18% (95% CI: 83.82–98.54%) and an optimal threshold of 0.55, achieving balanced sensitivity (0.91), specificity (0.86), and accuracy (0.89), alongside strong predictive values (PPV: 0.91, NPV: 0.86). Model 3 ranked second in AUC (87.61, 95% CI: 77.49–97.73%) with the highest sensitivity (0.95) at a threshold of 0.48, moderate specificity (0.83) and accuracy (0.90). In contrast, Model 2 prioritized specificity (0.90) over sensitivity (0.65) at a high threshold (0.71), yielding the highest PPV (0.90) but lower accuracy (0.75). Model 5 exhibited an AUC of 85.97% (95% CI: 76.35–95.59%) at a threshold of 0.51, with balanced sensitivity (0.86) and specificity (0.79). Model 1 showed the weakest performance (AUC: 74.62, 95% CI: 62.55–86.69%), characterized by high sensitivity (0.93) but low specificity (0.52).

Table 3
www.frontiersin.org

Table 3. ROC parameters of multidimensional predictive models.

Figure 4
Four-panel figure evaluating model performance. Panel A: ROC curves with AUC values for five models, showing sensitivity versus 1-specificity. Panel B: Decision curve analysis plotting net benefit versus risk threshold for the models and reference lines for

Figure 4. (A) ROC curve of multidimensional predictive models. AUC, area under the curve. (B) DCA of multidimensional predictive models. (C,D) Calibration curve of multidimensional predictive models.

Using five-fold cross-validation to evaluate model stability and generalizability (Figure 5), we observed minimal divergence between the primary ROC (0.746) and cross-validated ROC mean (0.7496) for Model 1. Model 2 exhibited a moderate reduction (primary ROC: 0.832; ROC mean: 0.7931). Model 3 showed greater performance reductions (primary ROC: 0.876; ROC mean: 0.8196). Model 4, which incorporated the broadest variable set, yielded the largest discrepancy between primary and cross-validated performance (primary ROC: 0.912; ROC mean: 0.8196). Conversely, Model 5 maintained high consistency (primary ROC: 0.860; ROC mean: 0.8417), indicating superior stability. The mean Brier scores for Model 3, Model 4, and Model 5 were consistently below 0.2.

Figure 5
Horizontal bar chart comparing five models' performance using five-fold cross-validated ROC and Brier scores. For each model, both scores and their cross-validated means are displayed. Model 1 has ROC values of 0.746 and 0.7496, and Brier scores of 0.198 and 0.2203. Model 2 shows ROC values of 0.832 and 0.7931, with Brier scores of 0.165 and 0.2013. Model 3 has ROC values of 0.876 and 0.8196, and Brier scores of 0.113 and 0.1661. Model 4 shows ROC values of 0.912 and 0.8196, and Brier scores of 0.109 and 0.1859. Model 5 has ROC values of 0.860 and 0.8417, and Brier scores of 0.142 and 0.1779.

Figure 5. The five-fold cross-validation performance.

DCA analysis for multidimensional predictive models

The DCA analysis of multidimensional predictive models was presented in Figure 4B. Model 4 demonstrated the highest net benefit across a broad range of probability thresholds, significantly outperforming the other models. Model 3 also showed a relatively high net benefit, but it was slightly lower than Model 4 in most cases. Model 2 had a moderate net benefit compared to other models. Model 5 exhibited a net benefit that was lower than Model 4 but higher than Model 1 across most probability thresholds. In contrast, Model 1 consistently exhibited the lowest net benefit across nearly all probability thresholds, indicating its limited effectiveness in this analysis.

Calibration curve analysis for multidimensional predictive models

The calibration curve analysis for the multidimensional predictive models, as illustrated in Figures 4C,D, evaluated the alignment between predicted probabilities and actual outcomes across the five models. Model 4 exhibited the best calibration, with its curve closely following the line of identity (ideal calibration), indicating that its predicted probabilities align well with observed outcomes across different probability thresholds. Model 3 also revealed relatively good calibration, though with slight deviations from the ideal line in certain regions. Model 2 exhibited moderate calibration performance, with larger deviations from the ideal line compared to Models 3 and 4. Model 5 achieved moderate calibration, with minor deviations in both low and high probability extremes. In contrast, Model 1 had the poorest calibration, as its curve deviated significantly from the line of identity across most probability thresholds.

Discussion

This study developed a multidimensional predictive model integrating clinical characteristics, SPECT/CT imaging data, and serum biomarkers to assess clinical activity in thyroid eye disease (TED). The results demonstrated that the comprehensive model (Model 4), which combined all three data dimensions, outperformed models relying on single or partial parameters, achieving an AUC of 91.18% with balanced sensitivity (0.91) and specificity (0.86). This highlights its potential utility in guiding clinical decision-making.

CAS, as an internationally recognized assessment of thyroid associated eye disease activity, has irreplaceable advantages such as simplicity, non-invasiveness, and low cost. However, due to its reliance on inspectors, CAS has been shown to have a certain degree of subjectivity and unreliable inter-observer variability (16). In our dataset, 8 of 72 eyes with a CAS score ≤2 exhibited orbital inflammatory edema and extraocular muscle involvement, with SPECT/CT-derived uptake ratio (UR) values exceeding 5.0. These findings suggest that exclusive reliance on CAS scores may fail to fully capture the underlying pathological processes of inflammatory activity. Thus, establishing a comprehensive assessment method that incorporates multiple factors is essential to enhance diagnostic accuracy.

Orbital SPECT/CT has incomparable advantages over orbital CT plain scan (9). SPECT/CT technology offers metabolic and functional insights into orbital tissues by measuring radioactive tracer uptake, thereby providing an intuitive reflection of orbital inflammatory activity (17). In accordance with previous studies (15, 18), our results showed that the orbital SPECT/CT UR in the active phase group were higher than those in the stable phase group, indicating that orbital SPECT/CT UR are meaningful for assessing the activity of TED patients. While SPECT/CT is effective for TED evaluation, it also has limitations compared to other imaging techniques like MRI, which excels in soft tissue resolution and multi-parametric imaging (18). Future research should explore combining SPECT/CT with other imaging modalities to leverage their complementary strengths and further enhance the accuracy of TED assessment.

SPECT/CT imaging effectively reflects the morphological and metabolic characteristics of local orbital tissues 22 (19). However, as it primarily captures cellular metabolic activity within the orbit, it may not fully represent systemic inflammatory processes or immune reactivity. Moreover, radionuclide physiological accumulation in the nasal sinus mucosa near the orbit and in tears can cause false-positive results in SPECT/CT imaging (20, 21). In our data, 18 of 72 eyes with UR values <5.0 showed CAS levels exceeding 3.0, highlighting the limitations of using UR as a sole indicator. TED disease activity is closely associated with clinical markers and blood biomarkers, such as TRAb, NLR, PLR, and SII, which play critical roles in disease progression and therapeutic response (2225). This understanding has driven the evolution from traditional imaging-based approaches toward more integrative models that incorporate systemic health indicators. By combining SPECT/CT findings with serum markers such as TRAb, we can develop a more comprehensive assessment framework. This multimodal approach not only enhances the characterization of disease activity but also improves diagnostic accuracy in distinguishing active from inactive phases of TED, thereby optimizing clinical decision-making.

In the present study, univariate analysis identified sex, SBP, EX, URmax on SPECT/CT, TRAb levels, and RBC count as significant predictors of TED activity. Regarding sex differences, our results corroborate previous findings. A meta-analysis identified a significant positive correlation between male sex and active TED (26), and a large Chinese study reported a substantially higher prevalence of active-stage TED among males compared to females (27). Notably, our investigation demonstrated a markedly elevated risk of active TED in males (OR = 11.74), reflecting sex-related differences in immune regulation or genetic susceptibility, which needs further mechanistic exploration. Specifically regarding SBP, our findings confirm its association with TED activity. This observation aligns with previous research suggesting that abnormal blood pressure, including elevated SBP in TED patients, may impact retinal vascular perfusion. Furthermore, elevated pulse pressure has been proposed to affect retinal vasculature prior to the onset of mechanical compression factors (28). The precise relationship between blood pressure dysregulation and TED pathogenesis warrants further investigation. EX, a key clinical feature of TED and a critical indicator for evaluating treatment efficacy (29), was also significantly associated with clinical activity in our cohort, consistent with established literature. TRAb, a core biomarker in TED pathogenesis, consistently correlated with disease activity, reinforcing its role in autoimmune dysregulation. Interestingly, RBC count was demonstrated significant differences between patients in the active and inactive phases of TED. Subsequent Lasso regression and BSR regression analyses consistently identified RBC count as a robust predictor, demonstrating its significant association with TED clinical activity. RBCs are known to interact with various immune cells and participate in immune responses (30). They may affect immune homeostasis through triggering vascular endothelial dysregulation and altering complement-mediated immunomodulation (31). Another potential mechanism involves RBCs in oxygen transport and tissue oxygenation, autoregulating blood flow via S-nitroso-hemoglobin (SNO-Hb) -mediated nitric oxide (NO) bioactivity export (32). TED is characterized by orbital inflammation, edema, and increased fibrosis, which can lead to reduced blood flow and hypoxia in the affected tissues. A recent study revealed that hypoxia enhances the differentiation of Th1 cells, Th17 cells, and antigen-specific CD4+ cytotoxic T lymphocytes (CTLs), representing a key mechanistic pathway in T-cell-mediated TED pathogenesis (33). RBCs are responsible for delivering oxygen to tissues, and any alteration in RBC function or number may impact tissue oxygenation in the orbit. Therefore, RBCs may indirectly contribute to the development of TED by influencing the oxygenation status of orbital tissues.

By employing BSR and LASSO regressions for variable selection, this study constructed five predictive models. Among the five models, Model 1 exclusively incorporated imaging data, yielding a diagnostic efficacy of merely 0.7462. This suggests that relying solely on SPECT/CT imaging indicators is insufficient for evaluating active and inactive TED phases. Model 2, which integrates clinical characteristics with imaging data, demonstrated a significant enhancement in diagnostic performance (primary ROC: 0.832). However, five-fold cross-validation revealed moderate performance reduction (ROC mean: 0.7931), indicating potential sensitivity to data variability. Model 3, combining only clinical and blood markers, achieved a diagnostic efficiency exceeding 0.8 (primary ROC: 0.876). This model is particularly advantageous in primary hospitals lacking advanced detection equipment like SPECT/CT, as it still enables effective determination of TED phases. Notably, its cross-validated ROC mean (0.8196) showed greater divergence from the primary value, suggesting limited generalizability despite clinical utility. Model 4, which comprehensively combines clinical, imaging, and blood markers, achieved the highest diagnostic efficiency of 0.9118, indicating its superiority as a multidimensional model for assessing TED activity. Additionally, Model 4 demonstrated high calibration accuracy and the highest net benefit across a broad threshold range in DCA, suggesting its clinical utility in balancing sensitivity and specificity. Nevertheless, it exhibited the largest cross-validation discrepancy (ROC mean: 0.8196), likely reflecting overfitting from its broad variable set. This underscores the trade-off between model comprehensiveness and generalizability. In contrast, Model 5, constructed by intersecting univariate and multivariate significant variables (p < 0.05), demonstrated robust primary accuracy (AUC: 0.860) with superior stability during cross-validation (ROC mean: 0.8417). This optimized variable selection strategy effectively balanced discrimination power with generalizability, positioning it as a reliable pragmatic tool, particularly for rapid clinical assessments where resource efficiency is critical.

Limitations

This study has some limitations. First, the limited sample size of 36 patients and the single-center retrospective nature might lead to selection bias and restrict the generalizability of the findings. Second, external validation in multicenter prospective cohorts is essential to confirm the model’s robustness and refine its predictive thresholds. Third, the biological mechanisms linking certain variables (e.g., RBC count, SBP) to TED activity remain unclear, necessitating mechanistic studies.

Conclusion

This study developed a multidimensional predictive model for assessing clinical activity in TED by integrating clinical data, SPECT/CT imaging, and serum biomarkers. The concise and comprehensive Model 5, incorporating statistically screened predictors (sex, TRAb, EX, URmax, RBC), achieved balanced diagnostic accuracy with enhanced stability. This suggests that a multidimensional approach can provide a more holistic and precise assessment of TED activity than any single modality alone, offering a promising direction for improving clinical decision-making.

Data availability statement

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

Ethics statement

The studies involving humans were approved by the Ethics Committee of the Hubei Provincial Hospital of Traditional Chinese Medicine. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

YL: Writing – original draft, Data curation, Validation, Conceptualization. G-HZ: Conceptualization, Writing – review & editing, Data curation. MT: Writing – original draft, Data curation. CH: Formal analysis, Writing – original draft. J-PZ: Supervision, Writing – original draft. Y-QH: Validation, Data curation, Writing – original draft, Methodology, Conceptualization, Software, Formal analysis. X-HZ: Methodology, Writing – review & editing, Conceptualization.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by the Natural Science Foundation of Hubei Province (Grant No. 2023AFD110), Hubei Provincial Administration of Traditional Chinese Medicine Research Project (Grant No. ZY2025L278), and the Shizhen Talent Program of Hubei Province for Scientific Research [Grant No. Hubei Health Document (2024) No. 256].

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 authors declare that no Gen AI was used in the creation of this manuscript.

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.

References

1. Kulbay, M, Tanya, SM, Tuli, N, Dahoud, J, Dahoud, A, Alsaleh, F, et al. A comprehensive review of thyroid eye disease pathogenesis: from immune dysregulations to novel diagnostic and therapeutic approaches. Int J Mol Sci. (2024) 25:11628. doi: 10.3390/ijms252111628

PubMed Abstract | Crossref Full Text | Google Scholar

2. Spadaro, JZ, and Kohli, AA. Pathogenesis of thyroid eye disease. Int Ophthalmol Clin. (2023) 63:65–80. doi: 10.1097/IIO.0000000000000464

PubMed Abstract | Crossref Full Text | Google Scholar

3. Weiler, DL. Thyroid eye disease: a review. Clin Exp Optom. (2017) 100:20–5. doi: 10.1111/cxo.12472

PubMed Abstract | Crossref Full Text | Google Scholar

4. Moledina, M, Damato, EM, and Lee, V. The changing landscape of thyroid eye disease: current clinical advances and future outlook. Eye. (2024) 38:1425–37. doi: 10.1038/s41433-024-02967-9

PubMed Abstract | Crossref Full Text | Google Scholar

5. Burch, HB, Perros, P, Bednarczuk, T, Cooper, DS, Dolman, PJ, Leung, AM, et al. Management of thyroid eye disease: a consensus statement by the American Thyroid Association and the European Thyroid Association. Thyroid. (2022) 32:1439–70. doi: 10.1089/thy.2022.0251

PubMed Abstract | Crossref Full Text | Google Scholar

6. Lai, KKH, Aljufairi, FMAA, Sebastian, JU, Lai, KH, Chan, KKW, Chin, JKY, et al. Low clinical activity score, “progressive” thyroid eye disease: presentations of 1439 patients from a tertiary centre in Hong Kong. Br J Ophthalmol. (2025) 109:bjo-2024-325346. doi: 10.1136/bjo-2024-325346

PubMed Abstract | Crossref Full Text | Google Scholar

7. Li, S, and Liu, Y-J. Application of quantitative orbital analysis to assess the activity of Graves’ ophthalmopathy. Am J Nucl Med Mol Imaging (2023) 13:259–268. Available online at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10774606/. (Accessed April 16, 2025)

Google Scholar

8. Szumowski, P, Abdelrazek, S, Żukowski, Ł, Mojsak, M, Sykała, M, Siewko, K, et al. Efficacy of 99mTc-DTPA SPECT/CT in diagnosing orbitopathy in Graves’ disease. BMC Endocr Disord. (2019) 19:10. doi: 10.1186/s12902-019-0340-0

PubMed Abstract | Crossref Full Text | Google Scholar

9. Liu, Y, Li, M, Chen, H, Liu, W, Hu, M, Hu, F, et al. Diagnostic precision in thyroid-associated ophthalmopathy using multi-center radiomics with 99mTc-DTPA SPECT/CT. Sci Rep. (2024) 14:25810. doi: 10.1038/s41598-024-76018-4

PubMed Abstract | Crossref Full Text | Google Scholar

10. Gong, J, Liu, J-Y, Jiang, Y-J, Sun, X-W, Zheng, B, and Nie, S-D. Fusion of quantitative imaging features and serum biomarkers to improve performance of computer-aided diagnosis scheme for lung cancer: a preliminary study. Med Phys. (2018) 45:5472–81. doi: 10.1002/mp.13237

PubMed Abstract | Crossref Full Text | Google Scholar

11. Deng, Z, Chen, L, Tan, J, Wang, S, Liu, D, Wang, J, et al. Combination model of thyrotrophin receptor antibody and volumetric orbital apex crowding index as an indicator of dysthyroid optic neuropathy. Dis Markers. (2021) 2021:1–7. doi: 10.1155/2021/9964232

PubMed Abstract | Crossref Full Text | Google Scholar

12. Yan, Q, Wu, M, Zhang, J, Yang, J, Lv, G, Qu, B, et al. MRI radiomics and nutritional-inflammatory biomarkers: a powerful combination for predicting progression-free survival in cervical cancer patients undergoing concurrent chemoradiotherapy. Cancer Imaging. (2024) 24:144. doi: 10.1186/s40644-024-00789-2

PubMed Abstract | Crossref Full Text | Google Scholar

13. Bartalena, L, Kahaly, GJ, Baldeschi, L, Dayan, CM, Eckstein, A, Marcocci, C, et al. The 2021 European group on Graves’ orbitopathy (EUGOGO) clinical practice guidelines for the medical management of Graves’ orbitopathy. Eur J Endocrinol. (2021) 185:G43–67. doi: 10.1530/EJE-21-0479

PubMed Abstract | Crossref Full Text | Google Scholar

14. Dolman, PJ. Grading severity and activity in thyroid eye disease. Ophthalmic Plast Reconstr Surg. (2018) 34:S34–40. doi: 10.1097/IOP.0000000000001150

PubMed Abstract | Crossref Full Text | Google Scholar

15. Jiang, C, Deng, Z, Huang, J, Deng, H, Tan, J, Li, X, et al. Monitoring and predicting treatment response of extraocular muscles in Grave’s orbitopathy by 99mTc-DTPA SPECT/CT. Front Med. (2021) 8:791131. doi: 10.3389/fmed.2021.791131

PubMed Abstract | Crossref Full Text | Google Scholar

16. Perros, P, Žarković, M, Pearce, SH, Razvi, S, Kolli, H, and Dickinson, AJ. Inter-observer variability of clinical activity score: assessments in patients with thyroid eye disease. Am J Ophthalmol. (2023) 252:94–100. doi: 10.1016/j.ajo.2023.03.027

PubMed Abstract | Crossref Full Text | Google Scholar

17. Yao, N, Li, L, Gao, Z, Zhao, C, Li, Y, Han, C, et al. Deep learning-based diagnosis of disease activity in patients with Graves’ orbitopathy using orbital SPECT/CT. Eur J Nucl Med Mol Imaging. (2023) 50:3666–74. doi: 10.1007/s00259-023-06312-2

PubMed Abstract | Crossref Full Text | Google Scholar

18. Hu, Z, Liu, J, Deng, H, Chen, N, Chen, L, Wang, S, et al. Evaluation of inflammatory activity of extraocular muscles in thyroid associated orbitopathy by [68Ga]DOTATATE PET/CT. Mol Imaging Biol. 27:120–30. doi: 10.1007/s11307-024-01970-6

Crossref Full Text | Google Scholar

19. Li, H, Wang, Y, Zhou, J, You, D, Song, L, Wang, M, et al. Evaluation of fibroinflammatory activity in thyroid eye disease using [18F]AlF-NOTA-FAPI-04 PET/CT: a prospective study. Eur J Nucl Med Mol Imaging. (2025) 52:1890–900. doi: 10.1007/s00259-024-07015-y

PubMed Abstract | Crossref Full Text | Google Scholar

20. Mulazimoglu, M, Koca, S, Tamam, MO, Uyanik, E, and Ozpacaci, T. False-positive findings in post-treatment iodine-131 whole-body scintigraphy in a nasolacrimal sac cyst, confirmed with SPECT/CT and MRI. Clin Nucl Med. (2011) 36:805–7. doi: 10.1097/RLU.0b013e318219b4b3

PubMed Abstract | Crossref Full Text | Google Scholar

21. Lee, JY, Song, HS, and Kang, M. Physiological uptake of radioactive iodine around an artificial eyeball observed with single-photon emission computed tomography/computed tomography after radioactive iodine treatment. Nucl Med Mol Imaging. (2020) 54:204–6. doi: 10.1007/s13139-020-00651-w

Crossref Full Text | Google Scholar

22. Zhang, H, Fan, J, Qu, J, Han, Q, Zhou, H, and Song, X. Predictive markers for anti-inflammatory treatment response in thyroid eye disease. Front Endocrinol. (2023) 14:1292519. doi: 10.3389/fendo.2023.1292519

PubMed Abstract | Crossref Full Text | Google Scholar

23. Abounoori, M, Pourazizi, M, Bahmani Kashkouli, M, Akha, O, Jafari, R, and Movahedirad, M. Novel immunoinflammatory blood markers in Graves’ orbitopathy: insights into activity and severity. BMJ Open Ophthalmol. (2024) 9:e001744. doi: 10.1136/bmjophth-2024-001744

PubMed Abstract | Crossref Full Text | Google Scholar

24. Yuksel, N, Saritas, O, and Yuksel, E. Effect of thyroid hormone status on complete blood cell count-derived inflammatory biomarkers in patients with moderate-to-severe Graves’ ophthalmopathy. Int Ophthalmol. (2023) 43:3355–62. doi: 10.1007/s10792-023-02742-x

PubMed Abstract | Crossref Full Text | Google Scholar

25. Yeter, V, Koçak, N, Arslan, MT, and Kan, EK. Blood count-derived immunoinflammatory markers in thyroidassociated ophthalmopathy. Korean J Ophthalmol. (2021) 35:198–206. doi: 10.3341/kjo.2021.0012

PubMed Abstract | Crossref Full Text | Google Scholar

26. Lee, MH, Chin, YH, Ng, CH, Nistala, KRY, Ow, ZGW, Sundar, G, et al. Risk factors of thyroid eye disease. Endocr Pract. (2021) 27:245–53. doi: 10.1016/j.eprac.2020.11.011

PubMed Abstract | Crossref Full Text | Google Scholar

27. Yang, M, and He, W. Age and gender influence on clinical manifestations of thyroid-associated ophthalmopathy: a case series of 2479 Chinese patients. Front Endocrinol. (2024) 15:1434155. doi: 10.3389/fendo.2024.1434155

PubMed Abstract | Crossref Full Text | Google Scholar

28. Ye, J, Liu, W, Hu, X, Jiang, H, Xu, M, Jin, H, et al. Elevated pulse pressure correlated with reduced retinal peripapillary capillary in thyroid-associated ophthalmology with visual field defect. Front Endocrinol. (2022) 13:941051. doi: 10.3389/fendo.2022.941051

PubMed Abstract | Crossref Full Text | Google Scholar

29. Wiersinga, WM, Eckstein, AK, and Žarković, M. Thyroid eye disease (Graves’ orbitopathy): clinical presentation, epidemiology, pathogenesis, and management. Lancet Diabetes Endocrinol. (2025) 13:600–14. doi: 10.1016/S2213-8587(25)00066-X

PubMed Abstract | Crossref Full Text | Google Scholar

30. Buttari, B, Profumo, E, and Riganò, R. Crosstalk between red blood cells and the immune system and its impact on atherosclerosis. Biomed Res Int. (2015) 2015:616834. doi: 10.1155/2015/616834

PubMed Abstract | Crossref Full Text | Google Scholar

31. Dobkin, J, and Mangalmurti, NS. Immunomodulatory roles of red blood cells. Curr Opin Hematol. (2022) 29:306–9. doi: 10.1097/MOH.0000000000000734

PubMed Abstract | Crossref Full Text | Google Scholar

32. Reynolds, JD, Jenkins, T, Matto, F, Nazemian, R, Farhan, O, Morris, N, et al. Pharmacologic targeting of red blood cells to improve tissue oxygenation. Clin Pharmacol Ther. (2018) 104:553–63. doi: 10.1002/cpt.979

PubMed Abstract | Crossref Full Text | Google Scholar

33. Zhang, M, Qi, X, Zhou, X, Liu, Y, He, M, Wang, J, et al. A preclinical state of Graves’ ophthalmopathy characterized by hypoxia of T-cells identified via multi-omics analysis. J Clin Endocrinol Metab. (2025):dgaf230. doi: 10.1210/clinem/dgaf230

Crossref Full Text | Google Scholar

Keywords: thyroid eye disease, SPECT/CT, clinical activity score, serum biomarkers, multidimensional model

Citation: Li Y, Zhang G-H, Tian M, Hua C, Zhai J-P, He Y-Q and Zuo X-H (2025) Multidimensional predictive model for assessing clinical activity in thyroid eye disease. Front. Med. 12:1623286. doi: 10.3389/fmed.2025.1623286

Received: 05 May 2025; Accepted: 25 June 2025;
Published: 08 July 2025.

Edited by:

Ban Luo, Huazhong University of Science and Technology, China

Reviewed by:

Han Wu, Zhejiang University, China
Benping Zhang, Huazhong University of Science and Technology, China

Copyright © 2025 Li, Zhang, Tian, Hua, Zhai, He and Zuo. 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: Yan-Qiong He, aGV5YW5xaW9uZ0BoYmh0Y20uY29t; Xin-He Zuo, WnVveGluaGVAaGJodGNtLmNvbQ==

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.