Hostname: page-component-cd9895bd7-gxg78 Total loading time: 0 Render date: 2024-12-25T15:36:53.784Z Has data issue: true hasContentIssue false

Socio-demographic differences in the dietary inflammatory index from National Health and Nutrition Examination Survey 2005–2018: a comparison of multiple imputation versus complete case analysis

Published online by Cambridge University Press:  27 September 2024

Rachel J Meadows*
Affiliation:
Center for Epidemiology & Healthcare Delivery Research, JPS Health Network, 1500 South Main Street, Fort Worth, TX 76104, USA Division of Health Services Management and Policy, College of Public Health, The Ohio State University, Columbus, OH, USA
Electra D Paskett
Affiliation:
Division of Epidemiology, College of Public Health, The Ohio State University, Columbus, OH, USA College of Medicine, Comprehensive Cancer Center, The Ohio State University, Columbus, OH, USA
Julie K Bower
Affiliation:
Division of Epidemiology, College of Public Health, The Ohio State University, Columbus, OH, USA
Gail L Kaye
Affiliation:
Division of Health Behavior and Health Promotion, College of Public Health, The Ohio State University, Columbus, OH, USA
Stanley Lemeshow
Affiliation:
Division of Biostatistics, College of Public Health, The Ohio State University, Columbus, OH, USA
Randall E Harris
Affiliation:
Division of Epidemiology, College of Public Health, The Ohio State University, Columbus, OH, USA
*
*Corresponding author: Email [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Objective:

Studies using the dietary inflammatory index often perform complete case analyses (CCA) to handle missing data, which may reduce the sample size and increase the risk of bias. Furthermore, population-level socio-economic differences in the energy-adjusted dietary inflammatory index (E-DII) have not been recently studied. Therefore, we aimed to describe socio-demographic differences in E-DII scores among American adults and compare the results using two statistical approaches for handling missing data, i.e. CCA and multiple imputation (MI).

Design:

Cross-sectional analysis. E-DII scores were computed using a 24-hour dietary recall. Linear regression was used to compare the E-DII scores by age, sex, race/ethnicity, education and income using both CCA and MI.

Setting:

USA.

Participants:

This study included 34 547 non-Hispanic White, non-Hispanic Black and Hispanic adults aged ≥ 20 years from the 2005–2018 National Health and Nutrition Examination Survey.

Results:

The MI and CCA subpopulations comprised 34 547 and 23 955 participants, respectively. Overall, 57 % of the American adults reported 24-hour dietary intakes associated with inflammation. Both methods showed similar patterns wherein 24-hour dietary intakes associated with high inflammation were commonly reported among males, younger adults, non-Hispanic Black adults and those with lower education or income. Differences in point estimates between CCA and MI were mostly modest at ≤ 20 %.

Conclusions:

The two approaches for handling missing data produced comparable point estimates and 95 % CI. Differences in the E-DII scores by age, sex, race/ethnicity, education and income suggest that socio-economic disparities in health may be partially explained by the inflammatory potential of diet.

Type
Research Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2024. Published by Cambridge University Press on behalf of The Nutrition Society

Chronic inflammation contributes to the initiation and progression of numerous chronic health conditions that differentially burden populations based on demographic and socio-economic factors(Reference Furman, Campisi and Verdin1). Adverse social determinants of health can upregulate systemic inflammation from childhood through adulthood and lead to the early development of chronic diseases(Reference Furman, Campisi and Verdin1,Reference Nusslock and Miller2) . Nonetheless, dietary change is a promising strategy to modify chronic inflammation and thereby reduce chronic disease risks(Reference Calle and Andersen3). Dietary intake is often studied using dietary indices such as the Healthy Eating Index, Dietary Approaches to Stop Hypertension and Dietary Guidelines for Americans(Reference Eicher-Miller, Fulgoni and Keast4,Reference Frank, Jaacks and Avery5) . Substantial socio-economic and racial/ethnic differences in dietary behaviours have been long observed using these traditional dietary indices; however, these indices were not developed to directly measure the pro- or anti-inflammatory effects of diets. Dietary inflammation may be a key contributor to the socio-economic disparities observed in various chronic health conditions.

The dietary inflammatory index (DII) is a novel dietary assessment tool developed to directly measure the inflammatory potential of the overall diet(Reference Shivappa, Steck and Hurley6,Reference Marx, Veronese and Kelly7) . The DII characterises the overall inflammatory potential of a diet using up to forty-five dietary components and has been validated with various markers of inflammation in several populations(Reference Shivappa, Steck and Hurley8,Reference Shivappa, Wirth and Murphy9) . The DII correlates moderately (r = 0·52–0·65) with traditional dietary indices (i.e. the Healthy Eating Index and Dietary Approaches to Stop Hypertension)(Reference Wirth, Hébert and Shivappa10). The DII is sometimes adjusted for overall energy intake, referred to as the energy-adjusted DII (E-DII)(Reference Shivappa, Wirth and Murphy9), which assesses inflammatory potential per 4184 kJ (1000 kcal) of energy intake. The energy-adjusted version is critical for evaluating the effect of the inflammatory potential, independent of the overall energetic intake. To date, over 1200 studies using the DII/E-DII have been published; however, these studies are often limited by not examining missing data mechanisms, resulting in inappropriate handling of missing data, which can potentially result in bias and weaken the validity of the study conclusions(Reference Marx, Veronese and Kelly7,Reference Phillips, Chen and Heude11Reference Pedersen, Mikkelsen and Cronin-Fenton13) .

Missing data mechanisms are often categorised into the following three types: missing completely at random (MCAR), missing at random (MAR) and missing not at random (MNAR)(Reference Rubin14). MCAR occurs when the probability of missing data is the same for each participant; for example, if a weighing scale breaks, some participants have missing BMI data. The missing data are unrelated to the observed or unobserved data. MCAR is the strongest assumption and is unrealistic in typical epidemiologic studies(Reference Perkins, Cole and Harel15). When missing data are MCAR, complete case analysis (CCA) can be valid. Unlike MCAR, MAR occurs when missing data are associated with other collected data. For example, if missing data on alcohol intake were more common in younger age groups than in the older age groups, then the missing alcohol data would be MAR. When missing data are MAR, multiple imputation (MI) methods are recommended(Reference Pedersen, Mikkelsen and Cronin-Fenton13,Reference Harel, Mitchell and Perkins16,Reference Heymans and Twisk17) . Lastly, MNAR occurs when missing data depend on unobserved data, such as the variable itself. For example, if the missing data on sugar intake were more common among those who consume high levels of sugar than among others, regardless of age, sex or other observed variables, then the missing data on sugar intake would be MNAR. Because MNAR data depend on unobserved or unknown information, stronger assumptions are required. A common approach is to conduct sensitivity analyses wherein assumptions are varied and result in several possible conclusions(Reference Perkins, Cole and Harel15,Reference Staudt, Freyer-Adam and Ittermann18) .

Despite strong evidence demonstrating the importance and implementation of statistical methods for the appropriate analysis of missing data, the most basic and commonly used approach, CCA, continues to be used in epidemiologic studies, including dietary research(Reference Cole, Zivich and Edwards12,Reference Eekhout, de Boer and Twisk19,Reference English, Ard and Bailey20) . Using partial data that include only participants with complete data can exclude information that may be vital to studying disparities between subgroups. Because MAR, rather than MCAR or MNAR, is commonly observed in epidemiologic studies(Reference Pedersen, Mikkelsen and Cronin-Fenton13), this study focussed on comparing the basic approach of CCA with a more appropriate approach, i.e. MI(Reference Heymans and Twisk17). This study aimed (1) to investigate differences in E-DII scores by age, sex, race/ethnicity, education and income and (2) to compare the results using two statistical approaches for handling missing data, i.e. CCA and MI.

Methods

Study population

This study was a secondary data analysis of non-Hispanic White (NHW), non-Hispanic Black (NHB) and Hispanic adults aged ≥ 20 years, who participated in the National Health and Nutrition Examination Survey (NHANES) 2005–2018 cycles. The NHANES program includes interviews, physical examinations and laboratory tests. The NHANES uses a clustered design with multistage probability sampling to obtain a study sample representative of the civilian, non-institutionalised U.S. population. Detailed data collection methods for the NHANES are available on the NHANES website(21).

Measures

Questionnaire data

Demographic information

Self-reported information on sex (male or female), age, marital status (married, widowed, divorced, separated, never married or living with a partner), race/ethnicity (NHW, NHB or Hispanic including Mexican Americans and other Hispanics) and education level (less than 9th grade, 9–11th grade, high school graduate or equivalent, some college or associates degree or college graduate or above) was collected. Income was used to calculate the family income-to-poverty ratio (FIPR) based on the U.S. Department of Health and Human Services poverty guidelines for each year. FIPR was categorised based on eligibility levels for insurance subsidies under the Patient Protection and Affordable Care Act as follows: high income (FIPR ≥ 4), middle income (FIPR > 1 and < 4) and at or below the federal poverty level (FIPR ≤ 1)(22).

Dietary intake

The dietary assessment component of NHANES, ‘What We Eat in America’, is a joint programme between the U.S. Department of Health and Human Services and the U.S. Department of Agriculture. Dietary data were collected via a single 24-hour recall at the NHANES mobile examination centres. Interviewers used the U.S. Department of Agriculture’s Automated Multiple-Pass Method, a 5-step interview process, to record the dietary intakes efficiently and accurately(21,Reference Moshfegh, Rhodes and Baer23) . The U.S. Department of Agriculture’s Food Surveys Research Group then used the 24-hour recall to generate nutrient files for each participant(24). The ‘What We Eat in America’ data were used to estimate the prevalence of nutrient adequacy for foods, food groups, nutrients and dietary patterns and to inform nutrition programs and policies for the U.S. population(Reference Ahluwalia, Dwyer and Terry25).

For our study, nutrient files were used to calculate the DII scores; full details of the process have been described previously(Reference Shivappa, Steck and Hurley6). Briefly, the DII comprises forty-five food parameters, including individual nutrients (e.g. n-3 fatty acids), compounds (e.g. flavonoids), foods (e.g. garlic) and drinks (e.g. alcohol), which are associated with an inflammatory effect score. The forty-five food parameters and inflammatory effect scores were based on a review of 1943 published articles, and they were scored on their associations with the following inflammatory biomarkers: C-reactive protein, IL-1β, IL-4, IL-6, IL-10 and TNF-α (Reference Shivappa, Steck and Hurley6). The E-DII scores were calculated per 4184 kJ (1000 kcal) of food consumed to control for the effect of total energy intake, using the energy-standardised version of the global database. The theoretical range of possible values is from −9 to 8, where 0 represents a neutral diet, a score < 0 represents an anti-inflammatory potential of the diet, and a score > 0 represents a pro-inflammatory potential of the diet(Reference Shivappa, Steck and Hurley6). Nonetheless, studies often report ranges between −5 and +5(Reference Phillips, Chen and Heude11,Reference Hébert, Shivappa and Wirth26) . The E-DII has been validated in association with inflammatory biomarkers in several population-based studies(Reference Shivappa, Wirth and Murphy9,Reference Mohammadi, Hosseinikia and Ghaffarian-Bahraman27) .

The following twenty-eight of the forty-five E-DII parameters were available in the NHANES for the calculation of the overall E-DII scores: total energy, carbohydrate, protein, total fat, SFA, MUFA, PUFA, alcohol, fibre, cholesterol, niacin (vitamin B3), riboflavin (vitamin B2), thiamin (vitamin B1), Fe, Mg, Zn, Se, vitamin A, vitamin B6, vitamin B12, vitamin C, vitamin D, vitamin E, folic acid, beta carotene, n-6, n-3 and caffeine(Reference Shivappa, Wirth and Murphy9). Total energy intake was used to standardise the nutrient intake per 4184 kJ (1000 kcal) consumed. The remaining seventeen food parameters were not collected in the NHANES or were not collected consistently across all included years (trans-fat, anthocyanidins, eugenol, flavan-3-ol, flavones, flavonols, flavanones, isoflavones, garlic, ginger, onion, saffron, turmeric, pepper, thyme/oregano, rosemary and green/black tea). Prior studies have shown no significant change in the predictive ability when at least twenty-eight parameters are used compared with the full list of forty-five(Reference Shivappa, Steck and Hurley8). Published studies include twenty-seven food parameters on average, mainly because of the frequently limited representation of dietary information in structured questionnaires(Reference Phillips, Chen and Heude11,Reference Hébert, Shivappa and Wirth28) .

Chronic health conditions

Major chronic health conditions associated with chronic inflammation were self-reported individually by condition(Reference Minihane, Vinoy and Russell29). Indicator variables were created for each of the following conditions: diabetes, cancer, stroke, heart attack, angina, coronary heart disease, congestive heart failure, arthritis, liver disease, asthma, bronchitis and emphysema. Depressive symptoms were assessed using the nine-question Patient Health Questionnaire, employing a cut-off score ≥ 10, which indicates moderate to severe depressive symptoms(Reference Kroenke, Spitzer and Williams30).

Family history of diabetes, heart attack or angina

Participants were asked about their family history (father, mother, sisters or brothers only) for diabetes and heart attack or angina before the age of 50.

Food security

Ten questions from the U.S. Department of Agriculture Adult Food Security Module enquired about adult food insecurity in the last 12 months (e.g. scarcity of food, inability to afford to buy balanced meals and cutting size of meals or skipping meals due to lack of money)(Reference Bickel, Nord and Price31). The NHANES data categorises responses into the following categories: full food security (no affirmative responses), marginal food security (1–2 affirmative responses), low food security (3–5 affirmative responses) and very low food security (6–10 affirmative responses)(32). We merged low food security and very low food security into a single category to maintain sufficient cell size for analyses.

Lifestyle behaviours

Sleep: The participants were asked, ‘How much sleep do you usually get at night on weekdays or workdays?’ Sleep was categorised as short sleep (≤ 6 h/night), normal sleep (7–9 h/night) and long sleep (≥ 10 h).

Smoking: Smoking status was categorised as never smoker (< 100 cigarettes smoked in their lifetime but not currently smoking cigarettes), current smoker (≥ 100 cigarettes smoked in their lifetime and currently smoking cigarettes) or former smoker (≥ 100 cigarettes smoked in their lifetime but not currently smoking cigarettes)(33).

Physical activity: The Global Physical Activity Questionnaire was used to assess daily minutes of moderate- and vigorous-intensity physical activity for transportation (walking or biking), household or yard work and recreation(Reference Armstrong and Bull34). The metabolic equivalents per week were calculated, and the participants were categorised as ‘inactive’ (no reported physical activity), ‘somewhat active’ (> 0 to < 500 metabolic equivalents) and ‘active’ (meeting physical activity guidelines of ≥ 500 metabolic equivalents per week)(35).

Examination data

BMI

Trained health technicians measured participants’ height and weight. BMI was calculated as weight in kilograms divided by height in meters squared and categorised as follows: underweight (< 18·5 kg/m2), healthy weight (18·5–24·9 kg/m2), overweight (25·0–29·9 kg/m2) and obesity (≥ 30·0 kg/m2)(36).

Hypertension

Hypertension was defined as an average systolic blood pressure ≥ 140 mmHg, diastolic blood pressure ≥ 90 mmHg or self-reported use of antihypertensive medications.

Hypercholesterolaemia

Hypercholesterolaemia was defined as total serum cholesterol ≥ 240 mg/dl or self-reported use of cholesterol-lowering medications.

Analyses

Descriptive statistics

We calculated the weighted mean and proportions to describe the characteristics of the overall study population based on non-missing data(Reference Hardt, Herke and Brian37). We used box plots to report medians and interquartile ranges and graphically present the E-DII scores across the NHANES years (Fig. 1).

Fig. 1 Selection of participants for complete case analyses and multiple imputation subpopulations from National Health and Nutrition Examination Survey years 2005–2018. The complete case analyses subpopulation included only participants with no missing data on exposure, outcome or covariates of interest. The multiple imputation analyses subpopulation included all eligible participants, regardless of missing data on exposure, outcome or covariates of interest. NHANES, National Health and Nutrition Examination Survey; NH, non-Hispanic.

Complete case analyses and multiple imputation subpopulations

Complete case analyses subpopulation

This subpopulation was restricted to eligible participants with complete data for all exposure, outcome and covariate variables.

Multiple imputation subpopulation

This subpopulation included all eligible participants, regardless of missing data on exposure, outcome or covariates. Missing data for all the variables were imputed as described below.

Multiple imputation approach

MI is a flexible and commonly used approach that can accommodate variables with up to 50 % missing data(Reference Pedersen, Mikkelsen and Cronin-Fenton13,Reference Heymans and Twisk17) to help minimise bias in national surveys(Reference Reiter, Raghunathan and Kinney38,39) . We used MI by chained equations with the incorporation of NHANES sample weights to impute all missing data for variables stratified by sex (Table 1)(Reference Reiter, Raghunathan and Kinney38). The imputed variables included socio-demographic characteristics, family history, chronic conditions, lifestyle behaviours and E-DII scores. Missing E-DII values indicated adults who did not participate in the 24-hour dietary recall, and an overall E-DII score was imputed (i.e. individual E-DII parameters were not imputed). Among participants who reported a 24-hour dietary recall, variation in the number of food parameters was possible (e.g. Participant A reported a food intake that included 23 E-DII parameters, while Participant B reported a food intake that included twenty-seven E-DII parameters), but individual ‘missing’ E-DII parameters were not imputed. Notably, ‘missing’ individual E-DII parameters indicate that the participant did not report food intake that included that particular E-DII parameter; hence, it would not contribute to their overall E-DII score. Overall, each study variable had < 10 % missing information.

Table 1 Characteristics of American Non-Hispanic White, Non-Hispanic Black and Hispanic adults aged ≥ 20 years from the National Health and Nutrition Examination Survey 2005–2018 (n 34 547)

E-DII, energy-adjusted dietary inflammatory index.

* Includes non-responses and ‘Don’t know’ or ‘Refuse’ responses.

In this study, we generated fifteen imputed datasets, each containing a set of plausible replacement values for missing data points. Multiple imputed datasets were analysed, and overall estimates were obtained by averaging the estimates from each individually imputed dataset. The estimated standard errors and confidence intervals include an imputation variance component(Reference Pedersen, Mikkelsen and Cronin-Fenton13,Reference Harel, Mitchell and Perkins16,Reference Heymans and Twisk17,Reference Reiter, Raghunathan and Kinney38,39) .

Linear regression models

Linear regression models were used to estimate the association between sex and E-DII scores. We used sex-stratified models to assess independent associations between the 5-year age groups, race/ethnicity, education, FIP and E-DII scores. All models were run using the CCA and MI subpopulations to facilitate comparison. The absolute value of the percent change of β was calculated to compare coefficients between the CCA and MI subpopulation estimates for unadjusted and adjusted models. The percent change formula was as follows: | (CCA β – MI β)/CCA β | × 100.

We considered the following potential confounders: NHANES survey year, dietary intake day of the week, age, marital status, BMI, food security, family history of diabetes, family history of heart disease, chronic health conditions, physical activity, sleep and smoking status. In models with age as the main predictor, we categorised age into 5-year groups; in models with age as a covariate, age was modelled as a continuous variable (justified using fractional polynomials). The final models included significant confounders defined by the 10 % change-in-estimate criterion(Reference Rothman, Greenland and Lash40). We also tested whether differences in the unadjusted E-DII scores by sex and race/ethnicity differed by the NHANES year. All analyses were conducted using STATA version 16 (StataCorp LLC.)

We confirmed several assumptions of linear regression for the final models. Scatter plots were used to assess the linearity of the associations between exposure and outcome variables. Histograms were used to assess the distribution of the residual normality. The regression assumptions of homogeneity of variance and collinearity were verified using a visual inspection of residual v. predicted value plots and a variance inflation factor of < 5, respectively. The NHANES survey weights were incorporated into all analyses and were robust to heteroskedasticity, clustering, non-response errors and complex probability sampling designs(21).

Results

Based on non-missing data, Table 1 presents sample characteristics of the 34 547 NHW, NHB and Hispanic adults aged ≥ 20 years, who participated in the NHANES cycles 2005–2018. Among the included participants, 51 % were female and 49 % were male, with a mean age of 48 years. Furthermore, 73 % were NHW, 12 % were NHB and 15 % were Hispanic. In total, 36 % of the participants had an obese BMI, 31 % had an overweight BMI and 49 % reported at least one chronic health condition. Less than half (45 %) met the physical activity recommendations. Table 1 describes the missing data for each variable imputed in the regression models.

The MI and CCA subpopulations comprised 34 547 and 23 955 participants, respectively. Sex and BMI did not differ between the subpopulations (P = 0·373 and P = 0·24, respectively). Age (P < 0·001), race/ethnicity (P < 0·001), education (P < 0·01) and food security (P = 0·004) differed between the subpopulations although absolute differences were minimal. The CCA subpopulation had a slightly lower mean age (47·42 years v. 47·83 years), more NHW participants (74·8 % v. 72·5 %), more college graduates (29·46 % v. 28·02 %) and more individuals with full food security (77·47 % v. 76·79 %) compared with the MI subpopulation.

Overall, 57 % of adults reported E-DII scores indicative of a pro-inflammatory dietary potential, while only 34 % reported scores indicative of an anti-inflammatory dietary potential. The E-DII scores in this sample ranged from –5·96 (anti-inflammatory dietary potential) to a maximum score of 4·90 (pro-inflammatory dietary potential), and the overall mean score was pro-inflammatory (0·44 (95 % CI: 0·39, 0·49)). The average E-DII scores changed slightly over time, decreasing from 2005 to 2012 and increasing from 2012 to 2018 (Fig. 2). Nonetheless, the magnitude of temporal changes in the E-DII score was minimal (β = –0·025 (95 % CI: –0·053, 0·0020).

Fig. 2 Box plots of energy-adjusted dietary inflammatory index scores by two-year National Health and Nutrition Examination Survey cycles. The range of observed E-DII scores was -5·96 (anti-inflammatory dietary potential) to 4·90 (pro-inflammatory dietary potential). (Bars represent the interquartile range; dots represent outliers.) E-DII, energy-adjusted dietary inflammatory index.

Unadjusted associations

Similar patterns were observed in the CCA and MI subpopulations for all unadjusted comparisons of E-DII scores by sex, age, race/ethnicity, education and FIPR (Tables 2 and 3). Mean E-DII scores indicated greater anti-inflammatory potential of diet among women (β = 0·24 (95 % CI: 0·17, 0·30)) than among men (β = 0·66 (95 % CI: 0·60, 0·71)). Regression models using CCA and MI reflected similar patterns in the E-DII scores among women compared with men (CCA: β = –0·43 (95 % CI: –0·49, –0·38); MI: β = –0·42 (95 % CI: –0·47, –0·37); Table 2). The E-DII scores indicated a trend towards greater anti-inflammatory dietary potential with increasing age among both males and females. NHB adults had the highest average E-DII scores (β = 0·84 (95 % CI: 0·78, 0·90)), NHW adults had intermediate average scores (β = 0·39 (95 % CI: 0·33, 0·46)) and Hispanic adults had the lowest average scores (β = 0·35 (95 % CI: 0·30, 0·40)). The E-DII scores indicated a greater anti-inflammatory dietary potential with higher levels of education and FIPR. Differences in unadjusted point estimates between the CCA and MI subpopulations were mostly modest at ≤ 20 %.

Table 2 Unadjusted and adjusted linear regression for the associations between sex and energy-adjusted dietary inflammatory index scores

E-DII, energy-adjusted dietary inflammatory index; CCA, complete case analyses; MI, multiple imputation.

* Absolute value for the percent change in point estimates between the CCA and MI models.

Adjusted for the following covariates: National Health and Nutrition Examination Survey year, dietary intake day of the week, BMI, food insecurity, family history of heart disease, family history of diabetes, comorbidities, smoking, physical activity and sleep.

Table 3 Unadjusted linear regression for the associations between socio-demographic factors and energy-adjusted dietary inflammatory index scores stratified by sex

CCA, complete case analyses; MI, multiple imputation; NHW, non-Hispanic White; NHB, non-Hispanic Black; FIPR, family income-to-poverty ratio.

* Absolute value for the percent change in point estimates between the CCA and MI models.

Adjusted associations

In the adjusted models (Table 4), all patterns of E-DII scores by sex, age, race/ethnicity, education and FIPR remained. Similar patterns were observed in the CCA and MI subpopulations for all adjusted comparisons. Differences in adjusted point estimates between the CCA and MI subpopulations were mostly modest at ≤ 20 %.

Table 4 Adjusted linear regression for the associations between socio-demographic factors and energy-adjusted dietary inflammatory index scores stratified by sex

CCA, complete case analyses; MI, multiple imputation; NHW, non-Hispanic White; NHB, non-Hispanic Black; FIPR, family income-to-poverty ratio.

All models are adjusted for the following covariates: National Health and Nutrition Examination Survey year, dietary intake day of the week, BMI, food security, family history of heart disease, family history of diabetes, comorbidities, smoking, physical activity and sleep.

* Absolute value for the percent change in point estimates between the CCA and MI models.

Adjusted for the additional covariates: race/ethnicity and family income-to-poverty ratio.

Adjusted for the additional covariates: age and family income-to-poverty ratio.

§ Adjusted for the additional covariates: age and race/ethnicity.

Discussion

This study estimated the differences in E-DII scores by sex, age, race/ethnicity, education and FIPR among a U.S. population-based sample of 34 547 NHW, NHB and Hispanic adults, using two different approaches for handling missing data. The CCA and MI approaches resulted in comparable associations between E-DII scores and sex, age, race/ethnicity, education and FIPR. The consistency of point estimates between CCA and MI supports the robustness of MI in handling missing data. The MI approach facilitated the inclusion of an additional 44 % of participants who were ineligible for CCA due to missing data. Overall, 57 % of American adults reported E-DII scores associated with a high potential for supporting pro-inflammatory processes. Anti-inflammatory E-DII scores were associated with increasing age. NHB adults and adults with low education levels and FIPR reported the highest intakes of diets with pro-inflammatory potential.

Our results should be interpreted while considering several limitations. First, diet was assessed using a single 24-hour recall, which may be subject to various biases (e.g. recall and social desirability bias) and does not fully capture long-term dietary intake or patterns. Nonetheless, the use of the automated multiple-pass method reduces reporting bias, and results from dietary studies using the NHANES are often consistent with those of studies that use more robust dietary intake measures(Reference Moshfegh, Rhodes and Baer23,Reference Hébert, Hurley and Steck41) . Additionally, the E-DII as a measure to assess dietary inflammation has some limitations. For instance, many food components may have U-shaped associations with health outcomes that are not fully accounted for in the E-DII, such as vitamins A and E and alcohol. Finally, the crude categorisations of sex, race/ethnicity and socio-economic status in this study may not have fully captured the true complexity of these factors. For example, binary sex and adult socio-economic status capture only a portion of the true range that can influence dietary behaviours. Nonetheless, the results of this study point to methodological improvements and call attention to subgroup differences, which are critical for assessing the impact of dietary behaviours on health outcomes.

Our study demonstrates the advantage of MI over CCA, which is critical for maintaining the national representativeness of the NHANES data. Many prior DII/E-DII studies using national-level data have used CCA, which limited the sample size by at least 20 %(Reference Marx, Veronese and Kelly7,Reference Phillips, Chen and Heude11,Reference Hébert, Shivappa and Wirth26) . Despite the similar associations between CCA and MI in this study, CCA is not guaranteed to be comparable across all associations of interest. The risk of bias from CCA may be higher depending on the outcomes of interest or missing data mechanisms. The reduced sample size also limits the ability to examine potential differences by subgroups or effect measure modifiers, which may be critical when studying dietary behaviour and health outcomes. Others have discussed the advantages of MI in producing unbiased and valid estimates of associations; however, MI has not been commonly applied in the DII/E-DII literature, particularly in studies using NHANES data(Reference Marx, Veronese and Kelly7,Reference Phillips, Chen and Heude11,Reference Pedersen, Mikkelsen and Cronin-Fenton13) . Second, our study described differences in the E-DII on a continuous scale. Many studies categorise the DII/E-DII scores into quartiles or percentiles, which has implications for the loss of information on variability and usefulness in moving the field forward to interventions. For example, the categorisation of the DII/E-DII scores into quartiles does not tell us what level of change in the DII/E-DII score is necessary to reduce the risk of a specific health outcome. This information could help inform the interventions needed to test the impact of the DII/E-DII scores on health outcomes.

The magnitude of socio-demographic differences in the E-DII scores was modest but generally reflected the patterns from other dietary index studies. We found that only 34 % of American adults reported a 24-hour diet intake associated with anti-inflammatory potential. This corroborates the findings of other studies suggesting that most American adults fall short of a healthy diet recommended by the U.S. Dietary Guidelines for Americans based on various dietary assessment measures(Reference Eicher-Miller, Fulgoni and Keast4,Reference Rehm, Peñalvo and Afshin42,Reference Ryu, Shivappa and Veronese43) . Similar to studies assessing diet quality using the healthy eating index, alternate healthy eating index, Mediterranean or dietary approaches to stop hypertension, our study reported the highest anti-inflammatory dietary potential among females, older adults, Hispanics and those with high education levels or income(Reference Eicher-Miller, Fulgoni and Keast4,Reference Ryu, Shivappa and Veronese43) . Nonetheless, the absolute differences in the E-DII scores by the socio-demographic factors were modest. The E-DII scores in this study ranged from −5·96 (anti-inflammatory) to 4·90 (pro-inflammatory), but the average differences between the socio-demographic groups were < 1·0. Previous studies have observed an impact of a 1-point increase in DII score on a 6 % increased risk of depression, 8 % increased risk of CVD, 7 % increased risk of colorectal cancer and 8 % increased risk of overall mortality(Reference Shivappa, Godos and Hebert44Reference Chen, Peng and Lian46). Unfortunately, the DII is not adjusted for overall energy intake; therefore, it is difficult to separate the independent effect of the inflammatory potential. Nevertheless, some prospective studies have found that a 1-point increase in the E-DII score is associated with a 5 % increased risk of all-cause mortality and 16 % increased odds of developing frailty(Reference Chen, Peng and Lian46Reference Liang, Feng and Shivappa48).

Randomised trials are needed to test the magnitude of change required in the E-DII score to reduce the risk of various adverse health outcomes(Reference Hébert, Shivappa and Wirth26,Reference Hariharan, Odjidja and Scott49) . Furthermore, long-term prospective studies are needed to test whether diets with anti-inflammatory E-DII scores confer greater health benefits than other dietary indices. Additionally, potential physiological differences may exist in the effect of an inflammatory diet on health outcomes, considering the presence of other inflammation-inducing factors, such as obesity, financial adversity, racism and discrimination(Reference Nusslock and Miller2,Reference Calle and Andersen3) . Therefore, future studies should investigate potential subgroup differences to clarify the heterogeneity in the associations between the DII/E-DII scores and health outcomes(Reference Marx, Veronese and Kelly7,Reference Phillips, Chen and Heude11,Reference Namazi, Larijani and Azadbakht50) .

Conclusions

Our cross-sectional study of 34 547 American adults compared two different approaches for handling missing data to investigate the socio-demographic differences in the E-DII scores. We found similar patterns of associations with each approach, although MI facilitated the inclusion of an additional 44 % of participants who were ineligible for CCA due to missing data. Based on the E-DII scores, > 50 % of American adults reported 24-hour diet intakes associated with pro-inflammatory potential, with modest differences by sex, age, race/ethnicity, education and income.

Acknowledgements

We would like to thank Matt Jans and Te-Ching Chen from the National Center for Health Statistics for their consultation on the study methods regarding the use of MI for NHANES data.

Financial support

No funding was involved in this work.

Conflicts of interest

Dr. Paskett reports funding to her institution as PI from Pfizer, Merck Foundation, Genentech, Guardant Health and AstraZeneca not related to this work. She is also on Advisory Boards for GSK and Merck, also not related. Other authors have no conflicts of interests to disclose.

Authorship

All authors contributed to study design, methods and interpretation. R. J. M. conducted analyses and wrote the first draft; all coauthors reviewed and commented on subsequent drafts. All authors have approve the final manuscript and agree to be accountable for all aspects of the work.

Ethics of human subject participation

This study was a secondary analysis of publicly available data. The study was conducted according to the guidelines laid down in the Declaration of Helsinki, and all procedures involving research study participants were approved by the National Center for Health Statistics Research Ethics Review Board. Written informed consent was obtained from all subjects/patients.

References

Furman, D, Campisi, J, Verdin, E et al. (2019) Chronic inflammation in the etiology of disease across the life span. Nat Med 25, 18221832.Google Scholar
Nusslock, R & Miller, GE (2016) Early-life adversity and physical and emotional health across the lifespan: a neuroimmune network hypothesis. Biol Psychiatry 80, 2332.Google Scholar
Calle, MC & Andersen, CJ (2019) Assessment of dietary patterns represents a potential, yet variable, measure of inflammatory status: a review and update. Dis Markers 2019, 3102870.Google Scholar
Eicher-Miller, HA, Fulgoni, VL III & Keast, DR (2015) Energy and nutrient intakes from processed foods differ by sex, income status, and race/ethnicity of US adults. J Acad Nutr Diet 115, 18.e6.Google Scholar
Frank, SM, Jaacks, LM, Avery, CL et al. (2024) Dietary quality and cardiometabolic indicators in the USA: a comparison of the Planetary Health Diet Index, Healthy Eating Index-2015, and dietary approaches to stop hypertension. PloS One 19, e0296069.Google Scholar
Shivappa, N, Steck, SE, Hurley, TG et al. (2014) Designing and developing a literature-derived, population-based Dietary Inflammatory Index. Public Health Nutr 17, 16891696.Google Scholar
Marx, W, Veronese, N, Kelly, JT et al. (2021) The Dietary Inflammatory Index and human health: an umbrella review of meta-analyses of observational studies. Adv Nutr 12, 16811690.Google Scholar
Shivappa, N, Steck, SE, Hurley, TG et al. (2014) A population-based Dietary Inflammatory Index predicts levels of C-reactive protein in the Seasonal Variation of Blood Cholesterol Study (SEASONS). Public Health Nutr 17, 18251833.Google Scholar
Shivappa, N, Wirth, MD, Murphy, EA et al. (2019) Association between the Dietary Inflammatory Index (DII) and urinary enterolignans and C-reactive protein from the National Health and Nutrition Examination Survey-2003–2008. Eur J Nutr 58, 797805.Google Scholar
Wirth, MD, Hébert, JR, Shivappa, N et al. (2016) Anti-inflammatory Dietary Inflammatory Index scores are associated with healthier scores on other dietary indices. Nutr Res 36, 214219.Google Scholar
Phillips, CM, Chen, LW, Heude, B et al. (2019) Dietary Inflammatory Index and non-communicable disease risk: a narrative review. Nutrients 11, 1873.Google Scholar
Cole, SR, Zivich, PN, Edwards, JK et al. (2023) Missing outcome data in epidemiologic studies. Am J Epidemiol 192, 610.Google Scholar
Pedersen, AB, Mikkelsen, EM, Cronin-Fenton, D et al. (2017) Missing data and multiple imputation in clinical epidemiological research. Clin Epidemiol 9, 157166.Google Scholar
Rubin, DB (1976) Inference and missing data. Biometrika 63, 581592.Google Scholar
Perkins, NJ, Cole, SR, Harel, O et al. (2018) Principled approaches to missing data in epidemiologic studies. Am J Epidemiol 187, 568575.Google Scholar
Harel, O, Mitchell, EM, Perkins, NJ et al. (2018) Multiple imputation for incomplete data in epidemiologic studies. Am J Epidemiol 187, 576584.Google Scholar
Heymans, MW & Twisk, JWR (2022) Handling missing data in clinical research. J Clin Epidemiol 151, 185188.Google Scholar
Staudt, A, Freyer-Adam, J, Ittermann, T et al. (2022) Sensitivity analyses for data missing at random v. missing not at random using latent growth modelling: a practical guide for randomised controlled trials. BMC Med Res Methodol 22, 250.Google Scholar
Eekhout, I, de Boer, RM, Twisk, JW et al. (2012) Missing data: a systematic review of how they are reported and handled. Epidemiology 23, 729732.Google Scholar
English, LK, Ard, JD, Bailey, RL et al. (2021) Evaluation of dietary patterns and all-cause mortality: a systematic review. JAMA Netw Open 4, e2122277.Google Scholar
National Center for Health Statistics (NCHS) (2023) NHANES Questionnaires, Datasets, and Related Documentation. https://wwwn.cdc.gov/nchs/nhanes/Default.aspx (accessed March 2024).Google Scholar
U.S. Department of Health and Human Services (2023) Prior HHS Poverty Guidelines and Federal Register References. https://aspe.hhs.gov/topics/poverty-economic-mobility/poverty-guidelines/prior-hhs-poverty-guidelines-federal-register-references (accessed March 2024).Google Scholar
Moshfegh, AJ, Rhodes, DG, Baer, DJ et al. (2008) The US Department of Agriculture Automated Multiple-Pass Method reduces bias in the collection of energy intakes. Am J Clin Nutr 88, 324332.Google Scholar
U.S. Department of Agriculture, Agricultural Research Service (2022) Food and Nutrient Database for Dietary Studies. https://www.ars.usda.gov/northeast-area/beltsville-md-bhnrc/beltsville-human-nutrition-research-center/food-surveys-research-group/docs/fndds/ (accessed March 2024).Google Scholar
Ahluwalia, N, Dwyer, J, Terry, A et al. (2016) Update on NHANES Dietary Data: focus on collection, release, analytical considerations, and uses to inform public policy. Adv Nutr 7, 121134.Google Scholar
Hébert, JR, Shivappa, N, Wirth, MD et al. (2019) Perspective: the Dietary Inflammatory Index (DII)-lessons learned, improvements made, and future directions. Adv Nutr 10, 185195.Google Scholar
Mohammadi, S, Hosseinikia, M, Ghaffarian-Bahraman, A et al. (2023) Dietary Inflammatory Index and elevated serum C-reactive protein: a systematic review and meta-analysis. Food Sci Nutr 11, 57865798.Google Scholar
Hébert, JR, Shivappa, N, Wirth, MD et al. (2020) Reply to FJB van Duijnhoven et al. Adv Nutr 11, 179180.Google Scholar
Minihane, AM, Vinoy, S, Russell, WR et al. (2015) Low-grade inflammation, diet composition and health: current research evidence and its translation. Br J Nutr 114, 9991012.Google Scholar
Kroenke, K, Spitzer, RL & Williams, JB (2001) The PHQ-9: validity of a brief depression severity measure. J Gen Intern Med 16, 606613.Google Scholar
Bickel, G, Nord, M, Price, C et al. (2000) Guide to Measuring Household Food Security. https://fns-prod.azureedge.us/sites/default/files/FSGuide.pdf (accessed March 2024).Google Scholar
National Center for Health Statistics (NCHS) (2023) National Health and Nutrition Examination Survey 2017–2018 Data Documentation, Codebook, and Frequencies. https://wwwn.cdc.gov/Nchs/Nhanes/2017-2018/FSQ_J.htm (accessed March 2024).Google Scholar
National Center for Health Statistics (NCHS) (2023) National Health Interview Survey Glossary. https://www.cdc.gov/nchs/nhis/tobacco/tobacco_glossary.htm (accessed March 2024).Google Scholar
Armstrong, T & Bull, F (2006) Development of the World Health Organization Global Physical Activity Questionnaire (GPAQ). J Public Health 14, 6670.Google Scholar
U.S. Department of Health and Human Services (2008) Physical Activity Guidelines for Americans, Publication No. U0036. https://health.gov/our-work/nutrition-physical-activity/physical-activity-guidelines/previous-guidelines/2008-physical-activity-guidelines (accessed March 2024).Google Scholar
Centers for Disease Control and Prevention (2022) About Adult BMI. https://www.cdc.gov/healthyweight/assessing/bmi/adult_bmi/index.html (accessed March 2024).Google Scholar
Hardt, J, Herke, M, Brian, T et al. (2013) Multiple imputation of missing data: a simulation study on a binary response. Open J Statistics 3, 370378.Google Scholar
Reiter, JP, Raghunathan, TE & Kinney, SK (2006) The importance of modeling the sampling design in multiple imputation for missing data. Surv Method 32, 143149.Google Scholar
National Center for Health Statistics (NCHS) (2016) NHANES 1999–2006 DXA Multiple Imputation Data Files. https://wwwn.cdc.gov/nchs/nhanes/dxa/dxa.aspx (accessed March 2024).Google Scholar
Rothman, K, Greenland, S & Lash, TL (2008) Modern Epidemiology, 3rd ed. Philadelphia: Lippincott Williams & Wilkins.Google Scholar
Hébert, JR, Hurley, TG, Steck, SE et al. (2014) Considering the value of dietary assessment data in informing nutrition-related health policy. Adv Nutr 5, 447455.Google Scholar
Rehm, CD, Peñalvo, JL, Afshin, A et al. (2016) Dietary intake among US adults, 1999–2012. JAMA 315, 25422553.Google Scholar
Ryu, S, Shivappa, N, Veronese, N et al. (2019) Secular trends in Dietary Inflammatory Index among adults in the United States, 1999–2014. Eur J Clin Nutr 73, 13431351.Google Scholar
Shivappa, N, Godos, J, Hebert, JR et al. (2018) Dietary Inflammatory Index and cardiovascular risk and mortality – a meta-analysis. Nutrients 10, 200.Google Scholar
Shivappa, N, Godos, J, Heberxt, JR et al. (2017) Dietary Inflammatory Index and colorectal cancer risk – a meta-analysis. Nutrients 9, 1043.Google Scholar
Chen, GQ, Peng, CL, Lian, Y et al. (2021) Association between Dietary Inflammatory Index and mental health: a systematic review and dose-response meta-analysis. Front Nutr 8, 662357.Google Scholar
Millar, CL, Dufour, AB, Shivappa, N et al. (2022) A proinflammatory diet is associated with increased odds of frailty after 12-year follow-up in a cohort of adults. Am J Clin Nutr 115, 334343.Google Scholar
Liang, Z, Feng, Y, Shivappa, N et al. (2022) Dietary Inflammatory Index and mortality from all causes, cardiovascular disease, and cancer: a prospective study. Cancers 14, 4609.Google Scholar
Hariharan, R, Odjidja, EN, Scott, D et al. (2022) The Dietary Inflammatory Index, obesity, type 2 diabetes, and cardiovascular risk factors and diseases. Obes Rev 23, e13349.Google Scholar
Namazi, N, Larijani, B & Azadbakht, L (2018) Dietary Inflammatory Index and its association with the risk of cardiovascular diseases, metabolic syndrome, and mortality: a systematic review and meta-analysis. Horm Metab Res 50, 345358.Google Scholar
Figure 0

Fig. 1 Selection of participants for complete case analyses and multiple imputation subpopulations from National Health and Nutrition Examination Survey years 2005–2018. The complete case analyses subpopulation included only participants with no missing data on exposure, outcome or covariates of interest. The multiple imputation analyses subpopulation included all eligible participants, regardless of missing data on exposure, outcome or covariates of interest. NHANES, National Health and Nutrition Examination Survey; NH, non-Hispanic.

Figure 1

Table 1 Characteristics of American Non-Hispanic White, Non-Hispanic Black and Hispanic adults aged ≥ 20 years from the National Health and Nutrition Examination Survey 2005–2018 (n 34 547)

Figure 2

Fig. 2 Box plots of energy-adjusted dietary inflammatory index scores by two-year National Health and Nutrition Examination Survey cycles. The range of observed E-DII scores was -5·96 (anti-inflammatory dietary potential) to 4·90 (pro-inflammatory dietary potential). (Bars represent the interquartile range; dots represent outliers.) E-DII, energy-adjusted dietary inflammatory index.

Figure 3

Table 2 Unadjusted and adjusted linear regression for the associations between sex and energy-adjusted dietary inflammatory index scores

Figure 4

Table 3 Unadjusted linear regression for the associations between socio-demographic factors and energy-adjusted dietary inflammatory index scores stratified by sex

Figure 5

Table 4 Adjusted linear regression for the associations between socio-demographic factors and energy-adjusted dietary inflammatory index scores stratified by sex