2.1 Study design
This research utilized information from the 2007-2012 National Health and Nutrition Examination Survey (NHANES), a cross-sectional study conducted by the National Centers for Disease Prevention16. Participants in the NHANES provided written informed consent and were approved by the National Center for Health Statistics Ethics Review Board17. Among the 30,442 subjects from 2007-2012, we excluded those who (1) were under the age of 20 (n = 12729), (2) had dietary and lifestyle OBS elemental deficiencies (n = 2033), (3) did not have thyroid function measurements (n = 7648), (4) reported a history of a previous diagnosis of thyroid disease (n = 2033), (5) TSH not in the range of 0.4 to 4.5 mUI/L and FT4 concentrations not in the range of 9 to 25 pmol/L (n = 1291), (6) pregnant at the time the survey was conducted (n = 64), (7) missing education level (n = 7), (8) missing iodine status (n = 186), and (9) missing mortality status (n = 7). 5727 participants remain in the current analysis (Participant Flowchart, Figure 1).
2.2 Oxidative balance score
The OBS was determined using a combination of 16 dietary nutrients and 4 lifestyle factors, encompassing 15 antioxidant elements and 5 pro-oxidant components14,18,19. Physical activity was assessed using the NHANES-recommended metabolic equivalent of task (MET) assignment, which includes one week of “(1) vigorous work-related activity, (2) moderate work-related activity, (3) walking or biking for transportation, (4) intense leisure-time physical activity, and (5) moderate leisure-time physical activity”. The resulting score was determined based on the activity level, with 0 points for low-intensity activity, 1 for moderate-intensity activity, and 2 for high-intensity activity. Participants were categorized into heavy drinkers (≥15 g/d for women and ≥30 g/d for men), non-heavy drinkers (0-15 g/d for women and 0-30 g/d for men), and non-drinkers, receiving 0, 1, and 2 points respectively20. The body mass index (BMI) was categorized based on the World Health Organization's (WHO) classification for obesity (BMI ≥ 30 kg/m2), overweight (BMI = 25-30 kg/m2), and normal weight (BMI < 25 kg/m2), with corresponding scores of 0-2. The remaining OBS antioxidant factors included dietary fiber, carotenoids, riboflavin, niacin, vitamin B6, total folate, vitamin B12, vitamin C, vitamin E, calcium, magnesium, zinc, copper, and selenium. The factors were divided into three categories using a three-quartile approach, and the scores for each category (1-3) ranged from 0 to 2. Higher scores indicated an increase in antioxidant levels. The remaining pro-oxidant factors included total iron and fat. Unlike antioxidants, the top tertile of pro-oxidants received a score of 0, while the bottom tertile received a score of 2. The overall OBS was calculated by adding all the components together. For more information on how each element was scored, please see (Table S1).
2.3 Measurement of serum FT4, TSH, and urinary iodine
FT4 and TSH were measured using immunoenzymatic assays21. Concerning the clinical practice guidelines of the American Association of Clinical Endocrinologists, we defined a participant with normal thyroid function as having a serum TSH concentration of 0.4-4.5 mUI/L and an FT4 concentration of 9-25 pmol/L22,23. The concentration of iodine in the urine was measured using mass spectrometry with inductively coupled plasma and dynamic reaction cells 21.
2.4 Ascertainment of Outcomes
The mortality status and time of death for NHANES participants were obtained through the National Death Index (NDI). The NCHS website provides corresponding methods for accessing this information (https://www.cdc.gov/nchs/data-linkage/mortality-public.htm). The International Classification of Diseases, 10th edition (ICD-10) was used to establish the cause of death. This study examined all-cause mortality as the primary outcome.
2.5 Covariates
We incorporated variables that have been demonstrated to impact thyroid function and the risk of mortality in prior research, such as information on gender, age, race/ethnicity, level of education, status of iodine intake, and health condition. Participants' iodine intake status was defined based on urinary iodine concentration (UIC), with UIC <100 μg/L as iodine deficiency, 100-299 μg/L as usual, and ≥300 μg/L as iodine overdose24; the definition of Diabetes Mellitus (DM) included a physician's report, glycosylated hemoglobin A1c (HbA1c) ≥ 6.5%, fasting blood glucose ≥ 7.0 μmol/L, or an oral glucose tolerance test (OGTT) ≥ 11.1 mmol/L. Hyperlipidemia is defined as total cholesterol ≥ 200 mg/dL, triglycerides ≥ 150 mg/dL, LDL ≥ 130 mg/dL, or HDL < 40 mg/dL25. The definition of hypertension includes self-reported diagnosis by a medical professional or an average of three consecutive measurements indicating systolic blood pressure ≥140 mmHg or diastolic blood pressure ≥90 mmHg. Cardiovascular disease was characterized by physician diagnoses obtained through personal interviews and included conditions such as coronary heart disease, angina, heart failure, heart disease, and stroke.
2.6 Statistical Analyses
In performing the statistical analyses, we followed the CDC recommendations and chose the appropriate weights to apply to the data analyses17. Descriptive analysis was conducted to examine the characteristics of participants, as well as their serum FT4, TSH levels, and overall health status. Continuous variables were presented as median values, while categorical variables were expressed as percentages. Multivariate linear regression models examined the association between OBS and FT4, TSH. Due to the skewed data distribution, natural logarithm conversion is performed for OBS and FT4, TSH. The results are then expressed as percent differences21. We selected the covariates for adjustment based on the findings from the previous NHANES study. Variance Inflation Factor (VIF) less than 10 was considered free of multicollinearity. We used Model 1 (unadjusted), Model 2 (adjusted for age, gender, and race), and Model 3 (adjusted for age, gender, race, education, and iodine status) to explore potential associations between OBS and normal thyroid function. The continuous OBS was converted into categorical variables using quartile methods, and the P value for trends was computed. Sensitivity analyses were performed to exclude participants with serum TPOAb (>9 IU / mL) and TgAb (>4 IU/mL) to minimize the effect of pre-existing immune disorders of thyroid tissue. In addition, we performed a stratified analysis. A multivariate COX proportional hazard model was employed to investigate the association between OBS, FT4, TSH, and all-cause mortality. Restricted cubic spline was employed to investigate the potential for a nonlinear relationship between these three variables and all-cause mortality. When examining the relationship between OBS and all-cause mortality, we also made adjustments for diabetes, hypertension, hyperlipidemia, and cardiovascular disease. In analyzing the association between FT4, TSH, and all-cause mortality, we further refined our model by including adjustments for smoking, alcohol consumption, and BMI. To assess whether the impact of OBS and all-cause mortality were influenced by FT4 or TSH, mediation effect analyses were conducted while controlling for the covariates including age, sex, race, education, diabetes mellitus, hypertension, hyperlipidemia, and Cardiovascular disease26. The R software (version 4.2.2) and EmpowerStats were utilized for all statistical analyses. Statistical significance was determined based on a two-sided P-value below 0.05.