---
title: "Link Between Metabolic Score, Insulin Resistance, Uric Acid, and Thyroid Nodules"
id: "plos-one-6-associations-of-the-metabolic-score-for-insulin-resistance-and-uric-acid-with"
canonical_url: "https://medichelpline.com/clinical-feed/plos-one-6-associations-of-the-metabolic-score-for-insulin-resistance-and-uric-acid-with"
content_type: "clinical_feed_article"
specialty: "Endocrinology"
source_name: "PLOS ONE (Medicine)"
source_url: "https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0358879"
published_at: "2026-09-21T14:00:00.000Z"
evidence_level: "Journal Feed"
license: "CC-BY-NC-4.0 / Informational Use"
---
# Link Between Metabolic Score, Insulin Resistance, Uric Acid, and Thyroid Nodules
## Provenance & Clinical Metadata
- **Canonical URL:** https://medichelpline.com/clinical-feed/plos-one-6-associations-of-the-metabolic-score-for-insulin-resistance-and-uric-acid-with
- **Specialty:** [Endocrinology](https://medichelpline.com/clinical-feed/endocrinology.md)
- **Primary Source:** PLOS ONE (Medicine)
- **Source URL:** [Original Journal Publication](https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0358879)
- **Published At:** 2026-09-21T14:00:00.000Z
- **Evidence Rating:** Journal Feed
## Executive GIST (TL;DR)
- The study investigates the associations of the **Metabolic Score for Insulin Resistance** (METS-IR) and **uric acid** (UA) with prevalent thyroid nodules. - A total of 5,232 participants from health screenings were analyzed, with significant findings for those in the highest quartiles of METS-IR and UA. - Both METS-IR and UA were linked to increased odds of thyroid nodules, with odds ratios of 3.928 and 3.627, respectively. - Nonlinear associations were observed across C-TIRADS categories, particularly for categories 2-3 and 4. - Systolic blood pressure (SBP) was evaluated as a potential mediator, showing proportional mediation effects for METS-IR and UA but limited effects in specific C-TIRADS subgroups. - Joint elevated levels of METS-IR and UA correlated with the highest odds of having thyroid nodules, although interaction analyses suggested limited additive synergy, indicating complexity in the associations.
## Clinical Analysis & Structured Key Points
Associations of the metabolic score for insulin resistance and uric acid with prevalent thyroid nodules: A cross-sectional study with systolic blood pressure mediation analysis and C-TIRADS category-specific analyses | PLOS One Browse Subject Areas ? Click through the PLOS taxonomy to find articles in your field. For more information about PLOS Subject Areas, click here . Article Authors Metrics Comments Media Coverage Reader Comments Figures Figures Abstract Objective To investigate the associations of metabolic score for insulin resistance (METS-IR) and uric acid (UA) with the odds of prevalent thyroid nodules, including nonlinear associations across C-TIRADS categories; to explore SBP-related indirect associations; and to evaluate the joint association and statistical interaction of METS-IR and UA. Methods Among 12,256 individuals initially screened at the Health Management Center of the First Hospital of Shanxi Medical University between January 2018 and December 2021, 5,232 participants were included in the final complete-case analysis. Data were retrieved from the hospital health examination database on January 25, 2026. METS-IR and UA were categorized into quartiles. Multivariable logistic regression, restricted cubic spline (RCS) regression, exploratory mediation analysis, joint classification analysis, and interaction analyses were performed. Results Compared with the lowest quartile, the highest quartiles of METS-IR and UA were associated with higher odds of prevalent thyroid nodules in the fully adjusted models (METS-IR: OR =3.928, 95% CI : 3.120 to 4.946; UA: OR =3.627, 95% CI : 2.839 to 4.634; both P < 0.001). RCS analyses showed nonlinear associations of METS-IR and UA with C-TIRADS 2–3 and C-TIRADS 4 thyroid nodules (all P for overall association <0.001; all P for nonlinearity ≤0.049). Exploratory mediation analyses yielded outcome-dependent findings: in the overall analysis, the estimated proportions mediated through SBP were 6.62% for METS-IR and 1.44% for UA; no significant indirect associations were observed for C-TIRADS 2–3 nodules, whereas only the indirect association for METS-IR was significant for C-TIRADS 4 nodules. Participants with concurrently elevated METS-IR and UA had the highest adjusted odds of prevalent thyroid nodules ( OR , 4.562; 95% CI , 3.753–5.544; P < 0.001); however, interaction analyses indicated a submultiplicative joint association and no evidence of positive additive interaction. Conclusion METS-IR and UA were associated with higher odds of prevalent thyroid nodules, with nonlinear associations observed across the exposure distributions. Concurrent elevations in METS-IR and UA identified a subgroup with the highest observed odds of thyroid nodules, although interaction analyses did not support positive synergy. Exploratory mediation analyses provided limited and outcome-dependent evidence of SBP-related indirect associations. The cross-sectional design precludes conclusions regarding temporal ordering or causal mediation. Citation: Ding J, Hu X, Sun M, Ge S (2026) Associations of the metabolic score for insulin resistance and uric acid with prevalent thyroid nodules: A cross-sectional study with systolic blood pressure mediation analysis and C-TIRADS category-specific analyses. PLoS One 21(9): e0358879. https://doi.org/10.1371/journal.pone.0358879 Editor: Guang Chen, Hong Kong Baptist University, HONG KONG Received: May 15, 2026; Accepted: September 6, 2026; Published: September 21, 2026 Copyright: © 2026 Ding et al. This is an open access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: All relevant data are within the paper and its Supporting information files. Funding: This work was supported by the 2025 Shanxi Provincial Public Management Major Special Research Project (Grant No. SXSGGGLYB2517). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist. 1. Introduction Thyroid nodules are among the most common endocrine disorders, and their detection rate has increased substantially in recent years [ 1 ]. According to the latest Chinese guidelines for the diagnosis and treatment of thyroid nodules, the prevalence of thyroid nodules in the adult is approximately 20.4%, and nearly 8%−16% of nodules are malignant [ 2 ]. Thyroid nodule occurrence has been associated with sex, age, genetic susceptibility, and iodine nutritional status [ 3 ]. Increasing evidence also suggests that metabolic abnormalities are associated with the presence and burden of thyroid nodules [ 4 ]. Insulin resistance (IR)is a central feature of metabolic syndrome [ 5 ] and has been reported to be associated with thyroid nodules after adjustment for conventional risk factors [ 6 ]. Homeostasis Model Assessment of Insulin Resistance (HOMA-IR) and Metabolic Score for Insulin Resistance (METS-IR) are commonly used surrogate markers of IR [ 7 , 8 ]. HOMA-IR is calculated from fasting glucose and insulin concentrations [ 9 ], whereas METS-IR incorporates fasting blood glucose (FBG), triglycerides (TG), body mass index (BMI), and high-density lipoprotein cholesterol (HDL-C) [ 10 ]. By incorporating both glycemic and lipid-related measures, METS-IR may capture metabolic abnormalities relevant to IR beyond those represented by HOMA-IR [ 11 , 12 ]. Previous studies have also reported that METS-IR has useful predictive performance for several metabolic outcome [ 13 ]. Elevated uric acid (UA), another common metabolic abnormality, has also been associated with a higher prevalence of thyroid nodules [ 14 , 15 ]. However, evidence regarding the joint association of IR-related metabolic dysfunction and UA with thyroid structural abnormalities remains limited. Both IR and elevated UA have been linked to higher blood pressure through mechanisms that may include sodium retention, activation of the renin-angiotensin-aldosterone system, endothelial dysfunction, and chronic low-grade inflammation [ 16 – 19 ]. Systolic blood pressure (SBP) may be particularly relevant to these associations [ 20 ]. In turn, elevated SBP has been associated with thyroid structural changes, potentially through impaired microcirculation, oxidative stress, and endothelial dysfunction [ 21 ]. Together, these observations provide a rationale for examining SBP as a potential intermediate factor in the associations of METS-IR and UA with thyroid nodules. Therefore, this study aimed to examine the separate and joint associations of METS-IR and UA with prevalent thyroid nodules, including nonlinear associations across C-TIRADS categories. We also explored SBP-related indirect associations and assessed whether the joint association of METS-IR and UA departed from multiplicativity or additivity. By evaluating these metabolic indicators together, this study may provide additional evidence on their associations with prevalent thyroid nodules. 2. Materials and methods 2.1. Study population A total of 5,232 individuals who underwent routine health examinations at the Health Management Center of the First Hospital of Shanxi Medical University between January 2018 and December 2021 were enrolled in this study, including 38.0% males and 62.0% females. All participants completed demographic data collection, biochemical examinations, and thyroid ultrasonography during the same health examination visit. The data used in this study were retrieved from the hospital health examination database on 25/01/2026. Based on these data, the associations between metabolic-related indicators and the prevalence of thyroid nodules were evaluated. The inclusion criteria were as follows: (1) age ≥18 years; (2) complete demographic information; and (3) availability of complete biochemical test results and thyroid ultrasonography data. The exclusion criteria were as follows: (1) presence of thyroid-related diseases, including hyperthyroidism, hypothyroidism, or Hashimoto’s thyroiditis; (2) use of medications for thyroid diseases, such as thyroid hormone preparations or iodine-containing agents, during the study period; (3) presence of severe organic diseases, including hepatic or renal diseases, cardiovascular or cerebrovascular diseases, and malignancies; and (4) pregnancy or lactation. During the study period, 12,256 individuals were initially screened. We sequentially excluded 152 individuals aged <18 years and 2,495 records because they were duplicate, incomplete, or otherwise invalid. Among the 2,495 excluded records, 43 were duplicate records, 140 had missing or invalid demographic data, 1,267 had missing or invalid biochemical data, and 1,045 had missing or invalid thyroid ultrasonography data. These exclusion categories were applied sequentially and were mutually exclusive. We subsequently excluded 1,586 individuals with thyroid-related diseases, 658 individuals using thyroid-related medications, 1,982 individuals with severe hepatic, renal, cardiovascular, or cerebrovascular diseases or malignancies, and 151 individuals who were pregnant or lactating. Finally, 5,232 participants with complete and valid records were included in the final analysis ( Fig 1 ). Download: PNG larger image TIFF original image Fig 1. Flow diagram of participant selection. https://doi.org/10.1371/journal.pone.0358879.g001 Among the 2,495 excluded records, 43 were duplicate records, 140 had missing or invalid demographic data, 1,267 had missing or invalid biochemical data, and 1,045 had missing or invalid thyroid ultrasonography data. Exclusion categories were applied sequentially and were mutually exclusive. This study was approved by the Ethics Committee of the First Hospital of Shanxi Medical University (Approval No. KYLL-2024-045). Written informed consent was obtained from all participants prior to enrollment. No formal sample size calculation was performed because this study included all eligible individuals with complete and valid records available during the predefined study period. 2.2. Data collection The study was initiated in 2018. Basic information for all participants was collected by trained nurses using a standardized procedure, including demographic characteristics, current medical history, past medical history, medication history, and family history. Anthropometric measurements were performed using an SK-CK ultrasonic physical examination device. Participants were instructed to remove their shoes and hats and wear light clothing while standing upright on the measurement platform. Height, weight, and body mass index (BMI) were subsequently recorded. Blood pressure was measured by trained nurses using an automated sphygmomanometer (HBP-9021; Omron Healthcare). Participants were required to rest in a seated position for at least 5 minutes before measurement. Blood pressure was measured twice, and the average value was used for analysis. For laboratory examinations, participants were instructed to fast for at least 8 hours overnight before venous blood samples were collected in the morning. Biochemical parameters, including total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), triglycerides (TG), fasting plasma glucose (FPG), and uric acid (UA), were analyzed using an AU5800 series automatic biochemical analyzer manufactured by Beckman Coulter. Thyroid ultrasonography was performed by qualified sonographers. Participants were examined in the supine position with the neck fully exposed. Thyroid nodules were assessed and recorded using a color Doppler ultrasound system (EPIQ 7C; Philips Healthcare). 2.3. Calculation of METS-IR The Metabolic Score for Insulin Resistance (METS-IR) was calculated using the following formula [ 22 ] where fasting plasma glucose (FPG), triglycerides (TG), and high-density lipoprotein cholesterol (HDL-C) were expressed in mg/dL, and body mass index (BMI) was expressed in kg/m 2 . 2.4. Diagnostic criteria According to the 2020 Chinese Guidelines for Ultrasound Malignancy Risk Stratification of Thyroid Nodules: C-TIRADS [ 23 ], thyroid nodules were evaluated using the Chinese Thyroid Imaging Reporting and Data System (C-TIRADS). Participants classified as C-TIRADS categories 2–4 were defined as the thyroid nodule group, whereas those classified as C-TIRADS category 1 or without ultrasonographic abnormalities were assigned to the non-nodule group. The number, size, and morphological characteristics of thyroid nodules were independently assessed and recorded by experienced radiologists. For further analysis, participants were additionally stratified according to clinically relevant risk categories. Specifically, C-TIRADS categories 2–3 were defined as the low-risk nodule subgroup, whereas category 4 was defined as the high-risk nodule subgroup. Because of the limited number of participants classified as C-TIRADS category 5, these cases were not included in the subgroup analysis. 2.5. Group definitions Participants were stratified according to quartiles of METS-IR and uric acid (UA). The quartiles of METS-IR were defined as follows: Q1 (<30.26), Q2 (≥30.26 to <34.84), Q3 (≥34.84 to <40.00), and Q4 (≥40.00). The quartiles of UA were defined as follows: Q1 (<254 μmol/L), Q2 (≥254 to <311 μmol/L), Q3 (≥311 to <384 μmol/L), and Q4 (≥384 μmol/L). In addition, participants were categorized into four groups according to the median levels of METS-IR (34.84) and UA (311 μmol/L): (1) METS-IR(-)UA(-) group, defined as METS-IR < 34.84 and UA < 311 μmol/L; (2) METS-IR(-)UA(+) group, defined as METS-IR < 34.84 and UA ≥ 311 μmol/L; (3) METS-IR(+)UA(-) group, defined as METS-IR ≥ 34.84 and UA < 311 μmol/L; and (4) METS-IR(+)UA(+) group, defined as METS-IR ≥ 34.84 and UA ≥ 311 μmol/L. 2.6 Quality control Strict quality control procedures were implemented throughout the study. Participants were enrolled in accordance with predefined inclusion and exclusion criteria. Prior to participation, all individuals were fully informed of the study objectives, procedures, and contents, and written informed consent was obtained. All medical staff involved in questionnaire administration, physical examinations, and laboratory testing received standardized training before the study. Standard operating procedures for health examinations were strictly followed during data collection. The same model of examination equipment was used for all participants, and all devices were uniformly calibrated before use. In addition, questionnaire data and examination results were independently verified and entered by two trained researchers. 3. Statistical analysis All statistical analyses were performed using R software (version 4.5.1) and IBM SPSS Statistics (version 27.0). Categorical variables were summarized as frequencies and percentages and compared using the chi-square test. Non-normally distributed continuous variables were summarized as medians and interquartile ranges and compared using the Kruskal-Wallis test. Analyses were restricted to participants with complete data; no statistical imputation was performed. All tests were two-sided, and P < 0.05 was considered statistically significant. Multivariable logistic regression was used to estimate the associations of METS-IR and UA with the odds of prevalent thyroid nodules. Model 1 was unadjusted, and Model 2 was adjusted for age and sex. For the METS-IR analysis, Model 3 was adjusted for age, sex, smoking history, alcohol consumption history, BUN, and UA. For the UA analysis, Model 3 was adjusted for age, sex, smoking history, alcohol consumption history, BUN, HDL-C, LDL-C, TG, TC, and FPG. SBP and DBP were not included in the primary association models because blood pressure variables may lie on the hypothesized pathways linking METS-IR and UA with thyroid nodules; SBP was examined separately as a potential mediator. Results are presented as odds ratios (ORs) with 95% confidence intervals (CIs). For the joint classification analysis, the four groups defined in Section 2.5 were examined, with the METS-IR(-)UA(-) group as the reference. Model 1 was unadjusted, Model 2 was adjusted for age and sex, and Model 3 was additionally adjusted for smoking history, alcohol consumption history, and BUN. Multiplicative interaction was assessed by including a product term between dichotomized METS-IR and UA in the fully adjusted logistic regression model. Additive interaction was explored using the relative excess risk due to interaction (RERI), attributable proportion (AP), and synergy index (S), calculated on the OR scale. Because thyroid nodules were relatively common in this study population, the OR-based additive interaction measures were interpreted as exploratory rather than as formal risk-scale measures. Multivariable restricted cubic spline (RCS) logistic regression was used to examine the potential nonlinear associations of METS-IR and UA with C-TIRADS 2–3 and C-TIRADS 4 nodules separately. Each RCS model used the same covariate set as the corresponding exposure-specific Model 3; that is, the METS-IR RCS models were adjusted for age, sex, smoking history, alcohol consumption history, BUN, and UA, whereas the UA RCS models were adjusted for age, sex, smoking history, alcohol consumption history, BUN, HDL-C, LDL-C, TG, TC, and FPG. Models with three to five knots were compared using the Akaike information criterion (AIC). Five knots were selected for both METS-IR models and for the UA model for C-TIRADS 2–3 nodules, whereas four knots were selected for the UA model for C-TIRADS 4 nodules. Category-specific median exposure values were used as the reference values (OR = 1): 34.58 and 33.75 for METS-IR and 309.00 and 297.00 μmol/L for UA in the C-TIRADS 2–3 and C-TIRADS 4 analyses, respectively. Overall and nonlinear associations were evaluated using Wald chi-square tests. Mediation analyses were performed using the mediation package in R and were considered exploratory. METS-IR and UA were standardized as z scores, with the exposure contrast defined as an increase from the mean to one standard deviation above the mean. A linear regression model was fitted for SBP as the mediator, and a logistic regression model was fitted for the outcome, including the standardized exposure, SBP, and covariates from the corresponding Model 3. In the overall analysis, participants with any C-TIRADS 2–4 thyroid nodule were compared with those without thyroid nodules (n = 5,232). In category-specific analyses, participants with C-TIRADS 2–3 nodules were compared with those without thyroid nodules after excluding participants with C-TIRADS 4 nodules (n = 5,013); participants with C-TIRADS 4 nodules were similarly compared after excluding those with C-TIRADS 2–3 nodules (n = 3,950). The average causal mediation effect (ACME), average direct effect (ADE), total effect, and proportion mediated were estimated on the probability scale, with percentile-based 95% CIs derived from 5,000 nonparametric bootstrap simulations. 4. Results 4.1. Comparison of clinical characteristics between the non-thyroid nodule and thyroid nodule groups A total of 5,232 participants were included in this study, comprising 3,731 individuals in the non-thyroid nodule group, 1,282 individuals in the C-TIRADS 2-3 group, and 219 individuals in the C-TIRADS 4 group. Among participants in the C-TIRADS 2-3 group, there were 639 males and 643 females, with females accounting for a slightly higher proportion. In the C-TIRADS 4 group, there were 123 males and 96 females, with males predominating. Age, body mass index (BMI), systolic blood pressure (SBP), diastolic blood pressure (DBP), fasting plasma glucose (FPG), total cholesterol (TC), triglycerides (TG), uric acid (UA), and metabo
## Related Clinical Research

- [Diagnostic accuracy of GH stimulation tests after hypothalamic–pituitary radiotherapy: systematic](https://medichelpline.com/clinical-feed/bmj-open-9-diagnostic-accuracy-of-growth-hormone-stimulation-tests-for-growth-hormone.md)
- [Scalable encapsulation of pluripotent stem cell-derived islets in high-concentration alginate bead](https://medichelpline.com/clinical-feed/biorxiv-13-bridging-scale-up-to-transplantation-pluripotent-stem-cell-derived-pancreatic.md)
- [INSR–BECN1–KIAA0825 complex prevents hepatocytic insulin receptor secretion and guards against mon](https://medichelpline.com/clinical-feed/biorxiv-17-an-insr-becn1-kiaa0825-complex-prevents-hepatocytic-insulin-receptor-secretion.md)
- [Premature adrenarche and long-term metabolic risk: association with insulin resistance and hyperan](https://medichelpline.com/clinical-feed/pubmed-42375004.md) (DOI: 10.1210/clinem/dgag258)
- [INSULIN Act Advances: Senate HELP Committee Passes Bill to Expand $35 Insulin Cap](https://medichelpline.com/clinical-feed/endocrine-news-0-endocrine-society-advocacy-win-senate-help-committee-passes-insulin-act.md)

## Navigation
- [← Back to Endocrinology Feed](https://medichelpline.com/clinical-feed/endocrinology.md)
- [← All Clinical Specialties](https://medichelpline.com/clinical-feed.md)
## Medical & Regulatory Disclaimer

> [!CAUTION]
> MedicHelpline content is structured for research, educational, and professional discovery purposes. It does not constitute individual medical advice, clinical diagnosis, or treatment recommendations.
> Always verify dosing, contraindications, and regulatory alerts against official product labeling and primary regulatory sources before clinical decision-making.