Study population
Data of the study population came from the NHANES, an ongoing national survey conducted by the Centers for Disease Control and Prevention that focused on Americans' dietary nutrition and general health. Signed informed consent from all participants before participating in the study, and all study protocols were approved by the National Center for Health Statistics' ethical review board. Detailed information about the NHANES database should be visited at https://www.cdc.gov/nchs/nhanes/index.htm. Specifically, data for this study was gathered from the NHANES 2005–2018. Our study only included individuals aged 20 years or above (N = 39,749). Following elimination for missing data on working status or PHQ-9 items (N = 5,639), DII (N = 9,085), and demographics (N = 1,862), finally 23,163 eligible participants left for analysis (Fig. 1).
[Insert Fig. 1 here]
Assessment of depression (outcome)
The PHQ-9 (Patient Health Questionnaire, 9 items) was used by the NHANES to evaluate depression. The PHQ-9 consist of nine signs and symptoms for depression: appetite problems, attention problems, depressed mood, fatigue, interest blank, sleep troubles, somatization disorders, suicidal thoughts, and worthlessness feelings [23]. Each item on a scale from "0" (not at all) to "3" (nearly every day), total score could range from 0 to 27. A PHQ-9 score was further divided into binary categories: depression with a score ≥ 10, while no depression with a score < 10 [24]. The sensitivity and specificity for detecting major depression was 88% at a cut-off of 10 [25].
Assessment of DII (exposure)
We used the modified version for DII calculation developed by Shivappa et al., the further standardized and specific calculation method has been detailed in prior studies [15, 26]. In the present study, 28 of 45 dietary nutrients were incorporated due to the NHANES data limited: alcohol, carbohydrates, caffeine, carotene, cholesterol, energy, fiber, folic acid, iron, magnesium, monounsaturated fatty acids, niacin, n-3 fatty acids, n-6 fatty acids, polyunsaturated fatty acids, protein, saturated fatty acids, selenium, thiamine, total fat, vitamin A, vitamin B2, vitamin B6, vitamin B12, vitamin C, vitamin D, vitamin E, and zinc [17]. Notably, vitamin D data were not contained in the NHANES 2005–2006, it would be ignored. According to DII, participants were divided into pro-inflammatory diet (DII ≥ 0) and anti-inflammatory diet (DII < 0) [27, 28].
Assessment of workforce participation (moderator)
We used four measures of workforce participation at the individual level: work status, hours worked, work types, and shift work [29, 30]. According to prior research [31], work status (working vs not working) in the NHANES was captured by participants reporting whether they were working at a job or business last week. Furthermore, individual disclosed the number of hours they worked last week at all jobs or businesses, and we used hours worked per week as a continuous indicator for working. Private employee, government (including federal, state, and local) employee, self-employed, and unpaid family business or farm were covered as different types of workforce participation. Shift work was assessed with the following response options: regular daytime, evening shifts, night shifts, rotating shifts, or another schedule. Evening and night shifts were combined by referring the study [9].
Assessment of covariates
Based on previous studies [9, 16, 18], demographics, health behaviors, and health status-related covariates that may affect the association between DII, workforce participation, and depression were included in our study. Demographics included age (≤ 39 years, 40–64 years, ≥ 65 years), gender (male, female), race (non-Hispanic White, non-Hispanic Black, Mexican American, other Hispanic, Other/Multiracial), education (less than high school, high school, some college, college or above), marital status (married/living with partner, widowed/divorced/separated, never married), and poverty income ratio (PIR, poor, not poor). According to the studies [28, 32, 33], health behaviors included drinking status (non-drinker, 1–5 drinks/month, 5–10 drinks/month, 10 + drinks/month), smoking status (never smoker, former smoker, current smoker), sleeping time (≤ 6 hours, 7–8 hours, ≥ 9 hours), physical activities (inactive, moderate, vigorous). Health status included body mass index (BMI, underweight, normal, overweight, obese), diabetes (yes, no), hypertension (yes, no), and self-reported perceived health (fair or poor, good, very good, excellent), by referring the research [32, 34]. Generally, a PIR value < 1 was considered as poor, a PIR value ≥ 1 was defined as not poor [23]. For detailed BMI classes: underweight (BMI < 18.5), normal (18.5 ≤ BMI < 25), overweight (25 ≤ BMI < 30), and obese (BMI ≥ 30) [32]. Self-reported diabetes, glycosylated hemoglobin (HbA1c) ≥ 6.5%, or fasting plasma glucose level ≥ 126 mg/dl were considered diabetes [35]. Self-reported high blood pressure told by doctors, the last two times measured average systolic/diastolic blood pressure of at least 140/90 mmHg, or antihypertensive medication being used were considered hypertension [7].
Statistical analysis
In our study, the statistical software R (version 4.4.1) and Storm Statistical Platform (https://www.medsta.cn/software) were used for all statistical analysis. A two-side P-value < 0.05 was considered statistically significant. We adopted the recommendations for accurate reporting in medical research statistics [36], using multiple imputation to fill in missing covariate data. In moderating analysis models, depression (Y, outcome) and DII (X, exposure) were entered as continuous variables, and workforce participation (M, moderator) was entered as a categorical variable [37, 38]. Model "2-way" from the R package "bruceR::PROCESS" was used to examine the moderating effect, with a 95% confidence interval (CI) evaluated by 1000 bootstrap resampling and 12345 seed. We employed restricted cubic spline (RCS) for the binary logistic models to assess the dose–response association of DII and hours worked per week with depression at different levels. In addition, we recruited subgroup analysis and interactions for categorical variables, gender, age, race, education, marital status, PIR, drinking status, smoking status, sleeping time, physical activities, BMI, diabetes, hypertension, and perceived health among covariates. Furthermore, in sensitivity analysis, data of non-multiple imputation and all adult participants were included, also the propensity score matching (PSM) methods were used to test our results stability.