According to the Centers for Disease Control and Prevention (CDC), diabetes constitutes a growing burden on the public’s health that affects 12 % of the US population (over 40 million people with diagnosed or undiagnosed diabetes), with approximately 1.5 million new cases diagnosed each year(1). Diabetes is considered the eighth leading cause of death in the USA, and the risk of premature death among individuals with diabetes is two to three times that of individuals without diabetes(2–Reference Gregg, Gu and Cheng4). Growing evidence suggests that a greater consumption of red meat and processed meat is associated with a higher risk of mortality(Reference Zheng, Li and Satija5,Reference Zhang, Hayden and Jackson6) , chronic diseases, including CHD(Reference Al-Shaar, Satija and Wang7), certain types of cancer(Reference Huang, Wang and Bogumil8,Reference Huang, Cao and Chen9) , diabetes(Reference Papier, Fensom and Knuppel10) and hypertension(Reference Zhang and Zhang11). Several non-pharmacological alternatives have been suggested to prevent diabetes and manage hyperglycaemia and insulin levels, including weight loss, increased physical activity and consumption of healthy diets rich in fruits, vegetables, whole grains and legumes, and low in saturated fats and red meat(Reference Würtz, Jakobsen and Bertoia12–Reference Kasović, Štefan and Kalčik15).
Findings from a previous study using three longitudinal cohorts of the US adult population, confirmed by a meta-analysis, showed that consumption of processed, unprocessed and total red meat was positively associated with the risk of type 2 diabetes(Reference Pan, Sun and Bernstein16). Other previous studies indicated that a greater consumption of red meat was associated with a higher risk of diabetes(Reference van Woudenbergh, Kuijsten and Tigcheler17–Reference Schulze, Manson and Willett21). Moreover, additional studies have suggested that red meat intake, including processed and unprocessed red meat, was associated with several biomarkers of metabolic alterations related to diabetes, including hyperglycaemia, insulin levels, HbA1c levels, and higher BMI(Reference Fretts, Follis and Nettleton22,Reference Hobbs-Grimmer, Givens and Lovegrove23) . The association between total, processed and unprocessed red meat intakes and diabetes remains understudied in a nationally representative sample of the US adult population, a country with among the highest processed and unprocessed red meat consumption, making it directly relevant to the US population and national health policies and guidelines(Reference Frank, Jaacks and Batis24). The present study is among the first to utilise data from the National Health and Nutrition Examination Survey (NHANES) cycles, thereby providing a more generalisable perspective. In addition, our substitution analysis using partition models contributes to the existing literature by offering novel insights into the potential benefits of substituting dietary proteins from foods of plant origin for different types of red meat in the prevention of diabetes. The present study aims to examine the associations of total, unprocessed and processed red meat intake with diabetes and to evaluate the substitution effects of other protein sources for red meat with diabetes. In exploratory analyses, we also assessed whether associations varied by poverty: income ratio (PIR) and food insecurity status.
Methods
Study population
The NHANES programme is designed to collect nationally representative health-related information to assess diet, medical history and other health-related behaviours of the non-institutionalised US population(25). The NHANES evolved from the previous National Health Examination Survey (NHES), which began in the 1960s. However, NHANES, as it is now known, was founded in the early 1970s by combining the NHES with a greater focus on nutrition and health(26). The nutrition data is critical for setting nutritional standards for the population and monitoring trends in healthy and unhealthy eating habits(Reference Ahluwalia, Dwyer and Terry27). The study population consisted of participants of the continuous NHANES 2003 to 2016, who completed a 24-h dietary recall. The NHANES is an ongoing series of nationally representative health surveys conducted by the National Center for Health Statistics of the CDC. The surveys use a multistage probability sampling design to develop a sample representative of the US population(28). The NHANES data are publicly available through the CDC website(29). The NHANES survey protocols are approved annually by the National Center for Health Statistics Research Ethics Review Board, and all participants provide written informed consent(25). Since NHANES data are de-identified and publicly available, the Institutional Review Board at the researchers’ institution does not consider this to be human subject research. Therefore, further human subjects’ approval was unnecessary for the present study.
The present study included participants aged 18 years or older (n 41 907) from the NHANES 2003–2016 data. Participants were excluded if they had borderline diabetes, as defined by the NHANES questionnaire response option ‘Doctor told you have diabetes’, or if they refused to answer or did not know their diabetes status (n 835), reported implausible daily energy intake levels (< 800 kcal or > 4200 kcal for men and < 500 kcal or > 3500 kcal for women) (n 1404), refused or do not know their family history of diabetes (n 670) or had missing dietary sample weights (n 4261). A total of 34 737 participants were included in the present analysis (Figure 1).
Flow chart of the study participants.

Assessment of diet
Participants who were examined during NHANES were also eligible to complete 24-h dietary recall interviews in which respondents reported all foods and beverages consumed in the previous 24 h. During NHANES 2003–2016, among adults aged 18 years and older, over 85 % completed both 24-h dietary recalls(Reference Herrick, Rossen and Parsons30). While not all eligible participants completed both recalls, more than 95 % completed at least one recall. On the first day (Day 1), dietary recall surveys were collected in person at the Mobile Examination Center by trained interviewers. The second day (Day 2) of dietary recall interviews was conducted by telephone, 3–10 days after the Mobile Examination Center interviews(31). The NHANES 24-h dietary recalls were collected using the computerised method of the US Department of Agriculture (USDA) Automated Multiple-Pass Method, which accounts for day-to-day variation(Reference Moshfegh, Rhodes and Baer32). Unprocessed red meat includes beef, veal, pork, lamb and game meat. We used the MyPyramid Equivalents Database 2.0 for 2003–2004 and the USDA Food Patterns Equivalents Database for 2005–2016 to estimate food patterns, including red and processed meat intake. The USDA Food Patterns Equivalents Database is the updated version of the MyPyramid Equivalents Database. Red and processed meat intakes were directly obtained from these NHANES-specific databases. Food pattern components are reported in cup equivalents for fruits, vegetables and dairy products and in ounce equivalents for grains and protein foods(33). Processed meat included frankfurters, sausage and luncheon meats. For this study, total red meat was defined as the combination of unprocessed red meat and processed meat. For participants with two 24-h dietary recalls (Day 1 and Day 2), dietary data were averaged; for those with only one recall, the single-day data were used.
Ascertainment of diabetes
For the present study, diabetes was defined using any of the following criteria: self-report of diabetes diagnosis by a physician or other health professional, having a fasting plasma glucose level of 126 mg/dl or more, an HbA1c level of 6·5 % or more, or use of antidiabetic drugs in the past 30 d prior to the survey interview(Reference ElSayed, Aleppo and Aroda34).
Assessment of covariates
The following covariates were directly extracted from the NAHNES databases: sociodemographic factors including age (years), sex (men/women), ethnicity-race (Mexican American, other Hispanic, Non-Hispanic White, other race-multi-racial), education (less than high school, high school degree, more than high school), marital status (never married/married), BMI (kg/m2), food insecurity status (yes/no), PIR (≤ 1·30, > 1·30), psychological distress (yes/no), lifestyle factors including smoking status (non-smokers, former, current smokers), physical activity expressed as metabolic equivalent of task (MET) minutes per week, alcohol intake (g/d), high cholesterol, which was defined based on self-reported physician diagnosis of high cholesterol (yes/no), family history of diabetes (yes/no) and CVD (yes/no) were recorded as separate binary variables, and total calories intake (kcal/d) were included. Other dietary-related variables, including but not limited to the consumption of fish and shellfish (serving/d), poultry (serving/d), vegetables (serving/d), fruit (serving/d) and sugar-sweetened beverages (kcal/d) were also included in this analysis. During the Mobile Examination Center, weight and height were measured by trained technicians. BMI was calculated as weight (in kilograms) divided by height (in metres) squared. The PIR is the ratio of family income to the poverty threshold. It is calculated by dividing the family income by poverty guidelines specific to family size, year and state. Participants with a PIR ≤ 1·30 are considered to have low income (below 130 % of the federal poverty line), while those with a PIR > 1·30 are considered to have higher income (above 130 % of the federal poverty line)(35). More than high school’ refers to individuals with some college or an associate degree, or those with a college degree or higher (bachelor’s degree or above). Psychological distress was defined as the use of anti-anxiety or antidepressant medications. NHANES measures of physical activity were self-reported using the Global Physical Activity Questionnaire. The total physical activity score (MET-minutes per week) was calculated by summing MET-minutes per week across all activity domains. This score was then categorised into three groups: < 600, 600–1199 and ≥ 1200 MET-minutes per week. Food security status was assessed using the Household Food Security Module developed by the USDA(Reference Nord, Andrews and Carlson36). The module included eighteen items for households with children younger than 18 years and ten items for households without children. Food insecurity status was treated as a binary variable. Consistent with previous studies, participants who reported three or more affirmative responses were classified as food insecure (coded as ‘yes’), corresponding to low or very low food security; otherwise, responses were coded as ‘no’ (food secure)(Reference Jun, Cowan and Dodd37).
Statistical analysis
All analyses followed the National Center for Health Statistics guidelines, using appropriate weights to account for the complex sampling design of the NHANES survey, including strata, cluster and mobile examination centre weights(38). To describe the characteristics of the study population, we performed univariable analyses. Imputation was performed with the hot-deck technique for participants with missing demographic and lifestyle variables using PROC SURVEYIMPUTE(Reference Zhang, Muscat and Kris-Etherton39). Data missingness ranged from 0·08 % for education level to 29·6 % for the physical activity variable. Distributions of all variables were examined before and after imputation, and the overall data structure was well preserved. Multivariable logistic regression models (PROC SURVEYLOGISTIC; SAS Institute) were used to examine the independent associations of total red meat, unprocessed and processed meat intake with diabetes, adjusting for potential confounders, including total energy intake, demographics, lifestyle factors and other dietary factors. We conducted a nutritional substitution analysis to compare the health impacts of substituting 1 serving/d of other protein sources, such as poultry, fish, eggs, dairy products and plant-based protein sources (nuts, seeds, legumes and soya) for total, processed and unprocessed red meat. To facilitate interpretation of the results and maintain consistency with previous studies, we combined all plant-based protein sources, as defined and reported in NHANES, into a single composite variable. Food substitution analyses complement overall dietary pattern analyses but do not replace them. However, substitution results are generally easier to interpret than analyses of single foods that do not specify a comparison. This distinction is important because increasing the intake of one food typically decreases the intake of another; therefore, the overall health impact depends on both foods(Reference Gardner, Mehta and Bernstein40). The dietary components were originally expressed using USDA food patterns (cup equivalents for fruits, vegetables and dairy products; ounce equivalents for grains and protein foods)(33). To improve interpretability, we converted these into servings/d using standardised factors (e.g. 1 serving of poultry, fish or soya = 3·5 oz; legumes = 2 oz; others, such as fruits and vegetables = 1 oz or 1 cup equivalent).
As in a previous study, one serving of unprocessed red meat equals 3 oz of beef, veal, pork, lamb or game; one serving of processed red meat equals 1·6 oz of sausages, frankfurters and luncheon meats made from beef, pork or poultry(Reference Gu, Drouin-Chartier and Sacks41).
The association of substituting 1 serving/d of protein sources for 1 serving/d of total, processed and unprocessed red meat with diabetes was assessed by including the alternative foods as continuous variables in the same multivariable logistic regression model after adjusting for total energy intake in addition to dietary and non-dietary factors(Reference Song and Giovannucci42). The difference in the regression coefficients of the two foods being compared, and their variances and covariances, was used to estimate the β coefficient and variance for the substitution effect, which was used to calculate OR and 95 % CI for the substitution effect(Reference Pan, Sun and Bernstein16,Reference Song and Giovannucci42–Reference Halton, Willett and Liu44) . Linear trends were tested for significance using the median intake for each quintile of total, processed and unprocessed red meat, modelled as a single continuous variable. Effect modification was tested by including the multiplicative interaction term between total red meat intake and sex, age, PIR, psychological distress and food insecurity. We conducted a sensitivity analysis by excluding individuals classified as experiencing psychological distress. The multivariable logistic regression results were presented as adjusted OR and the 95 % CI. Data were analysed in SAS software (version 9.4; SAS Institute) at a two-tailed α level of 0·05.
Results
We included a total of 34 737 participants (mean (sd) age of 45·8 (17·5) years) in the analysis, of whom 10·5 % had diabetes. Participants with the highest category of total red meat intake were more likely to be men (72·6 %), non-Hispanic White (69·4 %), younger, non-smokers (48·7 %), married (63·8 %), more educated (55·7 %) and physically more active; they also had higher alcohol intake, higher BMI and higher total calorie intake compared with participants in lower quintiles (Table 1). Approximately 14 % of the overall sample were food insecure, with food insecurity varying slightly across quintiles of total red meat intake but remaining generally similar overall (Table 1). Compared with the lower fifths of the distribution, participants in the highest fifth of total red meat intake reported a slightly higher prevalence of family history of diabetes and a lower prevalence of CVD history, except when compared with the lowest fifth (Table 1). Moreover, dairy products intake increased slightly across quintiles of total red meat intake. Although whole-grain intake was lower only relative to the bottom fifth, participants with higher total red meat consumption had lower intakes of fruit, fish, and poultry, plant-based proteins, and higher coffee consumption (Table 1).
Baseline characteristics of participants (n 34 737) according to quintiles of total red meat intake based on NHANES 2003–2016 data

SSB, sugar-sweetened beverages.
* Unweighted sample size (n); weighted mean (standard deviation) for continuous variable, and weighted proportion for categorical variable.
† Psychological distress was defined as the use of anti-anxiety or antidepressant medications; plant-based proteins (nuts, seeds, legumes and soya).
In this study, and after adjustment for major confounders including demographics, lifestyle factors, psychological distress, diet quality including intake of fruits and vegetables, sugar-sweetened beverages, coffee, whole grains, and major protein food sources such as poultry, eggs, dairy products, legumes, nuts, and seeds, and BMI (Model 4), intakes of total red meat, processed red meat and unprocessed red meat were all associated with higher odds of diabetes (Table 2). Participants in the highest quintile of total red meat, compared with those in the lowest quintile, had 49 % higher odds of diabetes (ORQ5 v. Q1 = 1·49, 95 % CI 1·22, 1·81; P-trend < 0·001). Comparing the highest quintiles, processed red meat intake was also associated with 47 % higher odds of diabetes (ORQ5 v. Q1 = 1·47, 95 % CI 1·17, 1·84; P-trend = 0·001), whereas unprocessed red meat intake was associated with 24 % higher odds of diabetes (ORQ5 v. Q1 = 1·24, 95 % CI 1·06, 1·44; P-trend = 0·006) (Table 2).
Associations between total, processed and unprocessed red meat intake and diabetes in NHANES 2003–2016, showing the odds ratios (ORs) and 95% confidence intervals (CIs)

* Dietary intake was collected using 24-h dietary recalls. One serving of unprocessed red meat equals 3 oz of beef, veal, pork, lamb or game meat; one serving of processed red meat equals 1·6 oz of sausages, frankfurters or luncheon meats made from beef, pork or poultry.
Model 1: Adjusted for age and sex.
Model 2: Further adjusted for model 1 covariates plus smoking status, alcohol intake, race/ethnicity, education, family poverty:income ratio, marital status, physical activity, family history of diabetes, high cholesterol, history of CVD, total energy intake, food insecurity status and survey year.
Model 3: Further adjusted for model 2 covariates plus dietary intakes of fruits and vegetables, poultry, fish, eggs, dairy products, plant-based protein sources (nuts, seeds, legumes and soya), whole grains and coffee.
Model 4: Further adjusted for model 3 covariates plus sugar-sweetened beverage (SSB) intake, psychological distress and BMI.
In continuous analyses (model 4), every 1 serving/d increment in total red meat was associated with 1·16 times higher odds of diabetes (95 % CI 1·07, 1·25) (Table 2). Moreover, every 1 serving/d increment of processed red meat was associated with 1·10 times higher odds of diabetes (95 % CI 1·03, 1·19), and that increment of unprocessed red meat was associated with 1·10 times higher odds of diabetes (95 % CI 1·02, 1·19) (Table 2).
In our exploratory analysis, we did not observe significant interactions between continuous total red meat, unprocessed red meat or processed red meat intake and PIR, or between these red meat types and food insecurity status (P-interaction > 0·05 for all).
Furthermore, we did not observe any significant interaction between total red meat intake and sex, age, psychological distress or food insecurity status (P-interaction > 0·05, for all).
In substitution analyses using 1 serving/d of plant-based protein sources as the replacement, replacing 1 serving/d of total red meat was associated with 14 % lower odds of diabetes (OR = 0·86, 95 % CI 0·79, 0·94, Figure 2, Panel A). Similarly, substitutions of 1 serving/d of poultry, dairy products and whole grains for total red meat were associated with 11–12 % lower odds of diabetes (poultry: OR = 0·88, 95 % CI 0·79, 0·99; dairy products: OR = 0·89, 95 % CI 0·82, 0·97; whole grains: OR = 0·89, 95 % CI 0·82, 0·98, Figure 2, Panel A). For processed red meat, substituting 1 serving/d of plant-based protein sources for 1 serving/d of processed red meat was associated with 9 % lower odds of diabetes (OR = 0·91, 95 % CI 0·84, 0·99; Figure 2, Panel B). Similarly, replacing 1 serving/d of unprocessed red meat with 1 serving/d of plant-based protein sources was associated with 9 % lower odds of diabetes (OR = 0·91, 95 % CI 0·84, 0·98, Figure 2, Panel C).
Associations between substituting one serving of other protein sources for red meat intake and the odds of diabetes in the NHANES (2003–2016). The models were adjusted for age, sex, smoking status, alcohol intake, race/ethnicity, education level, poverty: income ratio, marital status, physical activity, family history of diabetes, high cholesterol, history of CVD, total energy intake, food insecurity status, survey year, dietary intakes of fruits and vegetables, sugar-sweetened beverages (SSB), psychological distress and BMI, excluding the foods being substituted.

Discussion
In this large, nationally representative, unweighted sample of 34 747 US adults using NHANES data from 2003 to 2016, the overall weighted prevalence of diabetes was 10·5 % (95 % CI 9·97, 11·09). After adjustment for sociodemographic, lifestyle and dietary factors, higher intakes of total red meat, processed red meat and unprocessed red meat were all positively associated with higher odds of diabetes. More importantly, the findings were consistent after adjusting for BMI, suggesting that the observed associations were independent of BMI. Consistent with previous prospective cohort studies, all observed associations were significant, but the association was stronger for processed meat than for unprocessed red meat(Reference Würtz, Jakobsen and Bertoia12,Reference Pan, Sun and Bernstein16,Reference Gu, Drouin-Chartier and Sacks41) . In addition, excluding participants experiencing psychological distress did not alter the observed associations for different types of red meat intake. Our substitution analyses indicate that substituting red meat with alternative protein sources is associated with lower odds of diabetes. Specifically, substituting total red meat with poultry, dairy products, whole grains or plant-based protein sources was associated with 11–14 % lower odds, with plant-based protein sources showing the strongest association. For both processed and unprocessed red meat, replacing 1 serving/d with plant-based proteins was consistently associated with 9 % lower odds of diabetes, underscoring the benefit of limiting red meat intake and reinforcing the importance of alternative protein sources for diabetes prevention.
To the best of our knowledge, this study is among the first studies to examine the association between red meat intake and diabetes using a nationally representative sample of the US adult population. Furthermore, our food substitution analyses serve to complement, rather than replace, overall dietary pattern analyses. Substitution models are generally easier to interpret than analyses of single foods that do not specify a comparison. This distinction is important, as increasing the intake of one food necessarily decreases the intake of another; thus, the overall health impact depends on both foods(Reference Gardner, Mehta and Bernstein40).
Our findings are directly generalisable to the US population in terms of national policies and public health relevance. Findings from a previous study using three longitudinal cohorts with a comprehensive dietary assessment using a validated FFQ, confirmed by a meta-analysis, as well as a recent study using Nurses’ Health Study, II, and Health Professionals Follow-up Study showed that greater intakes of total processed and unprocessed red meat were positively associated with the risk of type 2 diabetes(Reference Pan, Sun and Bernstein16,Reference Gu, Drouin-Chartier and Sacks41) . Other previous studies indicated that a greater consumption of red meat is associated with a higher risk of diabetes(Reference van Woudenbergh, Kuijsten and Tigcheler17–Reference Schulze, Manson and Willett21). A recent study, the Nurses’ Health Study II and the Health Professionals Follow-up Study, also indicated that intake of red meat, including both processed and unprocessed red meat, was strongly associated with a higher risk of type 2 diabetes(Reference Gu, Drouin-Chartier and Sacks41). A previous study using the European Prospective Investigation into Cancer and Nutrition data indicated a positive association between high consumption of total red meat and incident type 2 diabetes(Reference Bendinelli, Palli and Masala45). Moreover, other studies have suggested that red meat consumption is associated with higher fasting blood glucose, insulin concentrations and HbA1c levels(Reference Fretts, Follis and Nettleton22,Reference Hobbs-Grimmer, Givens and Lovegrove23) .
The higher content of saturated fats and heme iron(Reference Rajpathak, Crandall and Wylie-Rosett46,Reference White and Collinson47) in red meat, in addition to the high levels of Na(Reference Männistö, Kontto and Kataja-Tuomola48) in processed meat, may contribute to adverse health effects and increase the risk of diabetes. These dietary components are known to be positively associated with endothelial dysfunction, inflammation, insulin resistance and oxidative stress, which can cause tissue damage, especially reduced pancreatic β-cell function, thereby providing a plausible biological mechanism for the observed associations(Reference Rajpathak, Crandall and Wylie-Rosett46,Reference White and Collinson47,Reference Peppa, Goldberg and Cai49–Reference Risérus51) .
The higher content of nitrates and Na in processed meat may contribute to adverse health effects, such as insulin resistance, and help explain the strongest relationship between red meat consumption and diabetes observed in the present study(Reference Micha, Michas and Mozaffarian19). Findings from our substitution analyses are consistent with previous studies and suggest that alternative major plant-based protein sources, such as nuts and legumes, fish and eggs, were associated with the most substantial lower odds of diabetes(Reference Würtz, Jakobsen and Bertoia12,Reference Gu, Drouin-Chartier and Sacks41) .
Limitations
Our study, however, has some limitations that need to be addressed. First, diet was assessed using up to two 24-h dietary recalls, which may not have adequately captured within-person variation in total, processed and unprocessed red meat intake. In this context, random within-person error likely attenuated the point estimates and reduced the precision of the observed associations in the present study. In addition, the FPED does not disaggregate processed meat into processed red and processed poultry, which may result in non-differential misclassification of exposure. Second, because of its cross-sectional design, reverse causality emerging from changes in diet due to recent symptoms is possible. Moreover, the cross-sectional nature of the study limits our ability to assess causality beyond observed associations. Third, as in any observational study, residual confounding may be present. However, we carefully controlled for a set of major potential confounders, including sociodemographics, major lifestyle factors and other dietary risk factors, including total energy intake. As mentioned earlier, this study supports previous findings, highlighting that the focus is on higher v. lower red meat consumption rather than red meat intake per se. Therefore, our findings should be interpreted with caution. Lastly, because of the NHANES diabetes questionnaire’s structure, we were unable to differentiate between type 1 diabetes and type 2 diabetes in this study.
Conclusions
In this study, greater intakes of total, processed and unprocessed red meat were associated with higher odds of diabetes, adding to the evidence of potential health concerns associated with higher red meat consumption compared with lower intake. Consistent with previous studies, substitution of alternative dietary components, particularly plant-based protein sources, poultry, dairy products, eggs or whole grains for total, processed or unprocessed red meat was associated with lower odds of diabetes in this study. These findings underscore the potential benefits of replacing red meat with alternative protein sources for diabetes prevention.
Acknowledgements
Not applicable.
D. M. B., PhD, MPH, was supported by the NIDDK training grant 1K01DK140625-01. The content is solely the authors’ responsibility and does not necessarily represent the official views of the National Institutes of Health.
Research design (project conception and development of the overall research plan): D. M. B. Data extraction and study oversight: D. M. B. Data analysis: D. M. B. Statistical analysis: D. M. B. Drafting the manuscript: D. M. B. Review and editing: Y. Z., T. Q., A. R. P-L., N. R-K., X. G., V. M. C. and D. M. B. All authors have read and approved the final manuscript.
The authors have no conflicts of interest to disclose.
The NHANES datasets used and/or analysed during the current study are publicly available at https://wwwn.cdc.gov/nchs/nhanes/.



