Red meat consumption (RMC) has been associated with a pro-inflammatory status, which in turn has been related to a higher risk of type 2 diabetes( Reference van Woudenbergh, Kuijsten and Tigcheler1 ), the metabolic syndrome( Reference Azadbakht and Esmaillzadeh2 ) and CHD( Reference Azadbakht and Esmaillzadeh2 ). The Multi-Ethnic Cohort Study and Nurses’ Health Study reported that a high RMC was linked to approximately 40 % greater risk of type 2 diabetes over a follow-up time of 14 and 4 years( Reference Steinbrecher, Erber and Grandinetti3 , Reference Pan, Sun and Bernstein4 ). However, findings are not consistent across studies as a non-significant association was also found between type 2 diabetes risk and red meat intake, especially between unprocessed red meat and diabetes risk( Reference Pan, Sun and Bernstein5 , Reference Ley, Sun and Willett6 ).
Despite the conflicting literature surrounding RMC and risk of type 2 diabetes, there is evidence to suggest that increased RMC is associated with increased whole-body and central adiposity( Reference Wang and Beydoun7 ). These findings may contribute to explaining the significant association between RMC and type 2 diabetes as central adiposity is a more sensitive predictor of type 2 diabetes and other obesity-related chronic diseases compared with BMI( Reference Wang, Rimm and Stampfer8 ).
The role of visceral fat as a causal factor connecting obesity and weight gain to the pathogenesis of insulin resistance (IR) and atherosclerosis is established( Reference Grundy9 ). Waist circumference (WC) has been proposed as a rapid and simple measurement for the assessment of abdominal adiposity but, like BMI, is not able to discriminate between subcutaneous and visceral abdominal fat depots( Reference Romero-Corral, Somers and Sierra-Johnson10 ). Hence, additional simple and integrated indexes have been recently proposed by combining physical (i.e. WC, BMI, thigh circumference, age) and biochemical measures (i.e. TAG, blood glucose or HDL). The lipid accumulation product (LAP) index is a marker of central fat accumulation derived from the measurements of WC and circulating TAG( Reference Onat, Avci and Barlan11 ), which has been proposed as a predictor of IR, the metabolic syndrome (MetS), type 2 diabetes mellitus and CVD( Reference Mirmiran, Bahadoran and Azizi12 , Reference Wakabayashi and Daimon13 ). The visceral adiposity index (VAI) is another indicator of adipose tissue distribution which has been used in the stratification of adult obesity phenotypes( Reference Du, Yu and Zhang14 ) and to improve the prediction of cardio-metabolic risk( Reference Amato, Giordano and Galia15 ).
Regression analysis is frequently used to evaluate the association between dietary factors and disease risk, but it may be characterised by a limited capacity to identify putative biological mechanisms that could possibly explain the association between red meat intake and the risk of high inflammation and impaired glucose control( Reference Fairchild and McDaniel16 ). Mediation analysis is a more sensitive statistical approach that can be used to explore and quantify the extent to which the relationship between an exposure and an outcome of interest is established through the effect of a third variable( Reference Fairchild and McDaniel16 , Reference Baron and Kenny17 ). The traditional approach to mediation analysis tends to produce a bias when the interaction between exposure and mediator is undefined( Reference Richiardi, Bellocco and Zugna18 , Reference Robins and Greenland19 ). In addition, unbiased valid estimates of direct and indirect effects can be obtained using the counterfactual framework in causal mediation analysis( Reference Richiardi, Bellocco and Zugna18 , Reference Robins and Greenland19 ). It is unclear to what extent the adjustment for adiposity modifies or attenuates the association between meat consumption and cardio-protective parameters. Mediation analysis could clarify the role of adiposity underlying the relation between meat consumption and cardio-protective factors.
Previous studies have investigated the association between red meat intake with biomarkers of inflammation and glucose/insulin metabolism with mixed results( Reference van Woudenbergh, Kuijsten and Tigcheler1 , Reference Azadbakht and Esmaillzadeh2 , Reference Babio, Sorlí and Bulló20 , Reference Montonen, Boeing and Fritsche21 ); however, none of these studies has attempted to identify the intermediate factors that connect the exposure to red meat with the selected health outcomes. The present analysis aims to specifically investigate the link between red meat intake with C-reactive protein (CRP) and glucose/insulin haemostasis and identify adiposity factors that may mediate these associations. These factors include markers of adiposity (WC, BMI, visceral adipose tissue (VAT)), LAP and VAI in a representative adult population of USA using National Health and Nutrition Examination Survey (NHANES) database. We hypothesised that a higher red meat intake would be associated with unfavourable concentrations of inflammatory and glucose/insulin haemostasis biomarkers among adults and that these associations would be partly or fully mediated by adiposity markers.
Methods
Population characteristics
The NHANES programme is implemented by the US National Center for Health Statistics (NCHS)( Reference Kalk and Joffe22 ). NHANES uses a complex, multistage and stratified sampling design to select a representative sample of the civilian and non-institutionalised resident population of the USA. The NCHS Research Ethics Review Board approved the NHANES protocol, and consent was obtained from all participants( Reference Kalk and Joffe22 ). The present study was based on analysis of data collected from 2005 to 2010. Data collection on demographics occurs through in-home-administered questionnaires, while anthropometric and biochemistry data are collected by trained personnel using mobile examination centres (MEC). More detailed information is available elsewhere( Reference Kalk and Joffe22 , 23 ).
For the assessment of height and weight during the physical examination, participants were dressed in underwear, disposable paper gowns and foam slippers. A digital scale (‘Mettler Toledo, Panther’) was used to measure weight to the nearest 100 g, a fixed stadiometer was used to measure height to the nearest millimetre. BMI was calculated as weight in kilograms divided by the square of height in metres. WC was measured at the iliac crest to the nearest millimetre( 23 ).
A blood specimen was drawn from the participant’s antecubital vein by a trained phlebotomist. Glycated Hb (HbA1c) was measured using a Tosoh A1C 2.2 Plus Glycohemoglobin Analyzer. Fasting blood glucose (FBG) was measured by a hexokinase method using a Roche/Hitachi 911 Analyzer and Roche Modular P Chemistry Analyzer. Insulin was measured using an ELISA immunoassay (Merocodia)( Reference Mazidi, Michos and Banach24 ). Other laboratory test details are available in the NHANES Laboratory/Medical Technologists Procedures Manual( 25 ). Apo B was measured by radial immunodiffusion( 23 ). Details on the measurement of CRP concentrations are available elsewhere( 23 ). Homoeostatic model assessment of IR (HOMA-IR), β-cell function (HOMA-B) and insulin sensitivity (HOMA-IS) were calculated as follows: HOMA-IR = (FBG (nmol/l) × insulin (mU/ml)/22·5), and HOMA-β = (20 × insulin (μU/ml))/(FBG (mmol/l) − 3·5)( Reference Musso, Gambino and Bo26 ). The TAG glucose (TyG) index was calculated as the ln(TAG (TAG, mg/dl) × FBG (mg/dl)/2)( Reference Simental-Mendia, Rodriguez-Moran and Guerrero-Romero27 ). Anthropometrically predicted VAT (apVAT) was predicted with sex-specific validated equations that included age, BMI, WC and thigh circumference ( Reference Samouda, Dutour and Chaumoitre28 ). The equation for men was 6 × WC – 4·41 × proximal thigh circumference + 1·19 × age – 213·65; and the equation for women was 2·15 × WC – 3·63 × proximal thigh + 1·46 × age + 6·22 × BMI – 92·713( Reference Samouda, Dutour and Chaumoitre28 ). VAI was calculated using sex-specific formulas: males (WC/39·68 + (1·88 × BMI)) × (TAG/1·03) × (1·31/HDL); females: (WC/36·58 + (1·89 × BMI)) × (TAG/0·81) × (1·52/HDL), where both TAG and HDL levels are expressed in mmol/l( Reference Amato, Giordano and Galia15 ). LAP was calculated as (WC − 65) × TAG in men, and (WC − 58) × TAG in women( Reference Onat, Avci and Barlan11 ). Smoking status was self-reported and participants classified as current smoker or not. Metabolic equivalent of task (MET) is used to measure the intensity level of physical activity and indicated the rate of energy consumption for a specific activity. MET is defined as 1 kcal/kg per h (4·184 kJ/kg per h) that is roughly equal to the energy cost of being at rest. Physical activity was categorised into three intensity levels based on MET score: light, moderate and vigorous( Reference Liu29 ). Subjects with diabetes were excluded from the study.
Red meat consumption
Dietary intake was assessed via a 24-h recall obtained by a trained interviewer during the MEC visit, using a computer-assisted dietary interview system with standardised probes, that is, the US Department of Agriculture Automated Multiple-Pass Method (AMPM)( Reference Ahluwalia, Andreeva and Kesse-Guyot30 , Reference Ahluwalia, Dwyer and Terry31 ). Briefly, information on the type and quantity of all food and beverages consumed in a 24-h period before the dietary interview (from midnight to midnight) was collected using the AMPM. The AMPM is designed to enhance complete and accurate data collection while reducing respondent burden( Reference Ahluwalia, Dwyer and Terry31 , Reference Moshfegh, Rhodes and Baer32 ). Detailed descriptions of the dietary interview methods are provided in the NHANES Dietary Interviewer’s Training Manual( Reference Bowman, Friday and Moshfegh33 ). The MyPyramid Equivalents Database for USDA Survey Food Codes was used to calculate RMC( Reference Bowman, Friday and Moshfegh33 ). In the present study, red meat intake was calculated as the sum of beef, pork, lamb, veal and game consumption and expressed as g/d.
Statistical analysis
Analyses were conducted using the SPSS software (version 22) according to the guidelines set forth by the Centers for Disease Control and Prevention for analysis of complex NHANES data sets, accounting for the masked variance and using the proposed weighting methodology( 34 ). We used means and standard errors of the mean for continuous measures (ANOVA) and percentages for categorical variables (χ 2). ANCOVA was used to compute age, race, energy intake and sex-adjusted means of IR markers or inflammation across quartiles of meat consumption.
The counterfactual framework assessed the total, direct and indirect effects of RMC on markers of IR or inflammation with BMI, WC, apVAT, VAI and LAP as mediators( Reference VanderWeele35 , Reference VanderWeele36 ). In this approach, the ‘total effect’ can be decomposed into a ‘direct effect’ (not mediated by BMI, WC, apVAT, VAI and LAP; online Supplementary Fig. S1) and an ‘indirect effect’ (mediated by BMI, WC, apVAT, VAI and LAP; online Supplementary Fig. S1). The analysis was conducted using the SPSS Macro developed by Preacher and Hayes( Reference Preacher and Hayes37 ). A product-of-coefficients test was used as it has the potential to detect significant mediation effects in the absence of a significant intervention effect( Reference VanderWeele35 , Reference VanderWeele36 ). In brief, the macro generates outputs that include the following steps. First, the total effect (γ-coefficient) of the exposure on the outcome variable (i.e. markers of IR or inflammation) is estimated by regressing the markers of IR or inflammation (outcomes) on RMC (independent variable) while adjusting for the covariates used in the first step but without adjusting for mediators. The ‘action theory’ test is then used to examine the effect of the exposure (meat consumption) on the hypothesised mediators (α-coefficient, BMI, WC, apVAT, VAI and LAP; online Supplementary Fig. S1). The ‘conceptual theory’ test examines the association between changes in the hypothesised mediators and changes in outcome variables (i.e. markers of IR or inflammation; β-coefficient, online Supplementary Fig. S1). The programme also estimates the direct (γ’-coefficient) and indirect (α × β product of coefficients) effects. The proportion of the mediation effect was calculated using the following equation (α × β/ (α × β + γ)). Full or complete mediation is present when the total effect (the γ’ path) is significant, the direct effect (the γ’ path) is NS and α × β is significant, whereas partly or incomplete mediation is present when the direct effect (the γ’ path) is also significant. Inconsistent mediation is when neither total nor direct effect is significant and α × β is significant( Reference MacKinnon, Fairchild and Fritz38 ). All estimates were adjusted for age, sex, race/ethnicity, educational, smoking and level of physical activity.
Results
General characteristics
A total of 16 621 subjects met the criteria for inclusion in the present analyses. Overall 8607 (48·3 %) participants were men and the mean age was 47·1 years. Non-Hispanic White (69·4 %) was the largest racial group and other Hispanic (4·5 %) the smallest racial group. Furthermore, 56·1 % of the participants were married, while 56·4 % had achieved more than high school education. Means for BMI, WC and apVAT were 28·7 (se 0·05) kg/m2, 98·2 (se 0·12) cm and 179·2 (se 1·18), respectively. Overall, 20·1 % were current smokers including 24·7 % of men and 15·7 % of women. Participants engaging in vigorous physical activity represented 5·2 % of the participants, and those engaging in little/no physical activity represented 24·3 %. Age, sex and race-adjusted mean of markers of IR and inflammation (high-sensitivity CRP, apo B, FBG, insulin, HOMA-IR, HOMA-B and TyG index) significantly increased across quartiles of RMC (all Ps <0·001; Table 1); HbA1c was not associated with RMC (Table 1).
hs-CRP, high-sensitivity C-reactive protein; HOMA-IR, homoeostatic model assessment of insulin resistance; HOMA-β, homoeostatic model assessment of β-cell function; HbA1c, Hb A1c; TyG index, TAG glucose index.
* P values for linear trend across quartiles of hs-CRP. Variables were compared across quartiles of red meat consumption using ANCOVA test.
Red meat intake, anthropometry, insulin resistance and inflammation
Action theory. After adjusting for covariates, a significant association was found between red meat intake and BMI (β: 0·345, P < 0·001), WC (β: 0·912, P < 0·001), apVAT (β: 3·27, P < 0·001), VAI (β: 0·054, P < 0·001) and LAP (β: 0·066, P < 0·001) (Table 2).
hs-CRP, high-sensitivity C-reactive protein; HOMA-IR, homoeostatic model assessment of insulin resistance ; HOMA-β, homoeostatic model assessment of β-cell function; HbA1c, glycated Hb; TyG index, TAG glucose index.
* All estimates were adjusted for age, sex, race/ethnicity, education, smoking and level of physical activity. Estimates for mediator and outcomes correspond to the regression coefficients α and γ, respectively, in online Supplementary Fig. S1.
Total effect. This was calculated by examining the association between red meat intake and markers of IR or inflammation in multivariate models without adjusting for potential mediators. Results showed that, with the exception of HbA1c and HOMA-β, all the markers of IR or inflammation were positively and significantly associated with red meat intake (all Ps < 0·04; Table 2).
Conceptual theory. This analysis tested the association between mediators (BMI, WC, apVAT, VAI and LAP) and markers of IR or inflammation; all potential mediators had significant and positive association with markers of IR or inflammation (all Ps < 0·001; Table 3).
hs-CRP, high-sensitivity C-reactive protein; HOMA-IR, homoeostatic model assessment of insulin resistance; HOMA-β, homoeostatic model assessment of β-cell function; HbA1c, glycated Hb; TyG index, TAG glucose index.
* All estimates were adjusted for age, sex, race/ethnicity, educational, smoking and level of physical activity. Regression coefficient β is shown in online Supplementary Fig. S1.
Direct and indirect effects of red meat consumption on insulin resistance and inflammation
Table 4 shows the direct effect, indirect effect, proportion of mediation effect and Sobel statistics for testing indirect effects. Both BMI and WC significantly mediated the association between markers of IR and inflammation and red meat intake (all Ps < 0·001); BMI and WC showed the greatest effect on FBG (β = 0·312 and β = 0·371, respectively). apVAT was a significant mediator for the association between red meat intake and CRP, FBG, HbA1c and TyG index (all Ps < 0·001); similarly FBG was the variable with the strongest association with apVAT (β = 0·293). Both VAI and LAP mediated the association between red meat intake and markers of IR and inflammation (all Ps < 0·001); serum apo B was the variable with the strongest association with VAI and LAP (β = 0·682 and β = 0·808, respectively).
hs-CRP, high-sensitivity C-reactive protein; HOMA-IR, homoeostatic model assessment of insulin resistance; HOMA-β, homoeostatic model assessment of β-cell function; HbA1c, glycated Hb; TyG index, TAG glucose index.
* All estimates were adjusted for age, sex, race/ethnicity, educational, smoking and level of physical activity. Regression coefficients α, β and γ , are shown in online Supplementary Fig. S1.
Discussion
In the present study, we demonstrated that red meat intake was significantly associated with all anthropometric outcomes (BMI, WC, apVAT, VAI and LAP) in fully adjusted models. Red meat intake was also significantly associated with markers of IR (except for HbA1c and HOMA-B) and with inflammation. In addition, mediation analyses suggested that these significant associations were partly or fully mediated by central adiposity.
A systematic review reported a significant link between red meat intake, especially processed varieties, with risk of breast cancer( Reference Farvid, Stern and Norat39 ). Existing observational and intervention studies testing the association between RMC and CRP levels have reported mixed results. Findings were significant in some studies( Reference Azadbakht and Esmaillzadeh2 , Reference Ley, Sun and Willett6 , Reference Montonen, Boeing and Fritsche21 ), which is in line with that found in the present study. However, other studies have findings contradictory to our study findings, which demonstrated that while processed meat was positively associated with CRP, RMC alone was not( Reference van Woudenbergh, Kuijsten and Tigcheler1 ) and in another study where lean red meat was not associated with CRP levels( Reference Hodgson, Ward and Burke40 ). The content of cholesterol( Reference Damião, Castro and Cardoso41 ), Fe( Reference Wagener, Volk and Willis42 ) and SFA( Reference Damião, Castro and Cardoso41 , Reference Esmaillzadeh, Kimiagar and Mehrabi43 ) in red meat to some extent explains the association between red meat and adverse health outcomes( Reference Krafft, Perewusnyk and Hänseler44 , Reference Schümann, Kroll and Weiss45 ). Previous studies on the effect of meat consumption on glucose/insulin homoeostasis have been inconsistent, with some finding an association, in line with our study( Reference Pan, Sun and Bernstein5 , Reference Feskens, Sluik and van Woudenbergh46 – Reference van Dam, Rimm and Willett48 ), and others failing to show such an association( Reference Babio, Sorlí and Bulló20 , Reference Montonen, Boeing and Fritsche21 ). It has been reported that RMC may have an impact on glucose/insulin metabolism through Fe-related metabolic pathways( Reference Fernández-Real, López-Bermejo and Ricart49 ). Fe is a strong pro-oxidant that catalyses several cellular reactions involved in the production of reactive oxygen species and hence increases the oxidative stress level( Reference Rajpathak, Crandall and Wylie-Rosett50 ). This can cause damage to cellular structures, including pancreatic beta cells, and high body Fe stores have been shown to be associated with an elevated risk of diabetes( Reference Rajpathak, Crandall and Wylie-Rosett50 ). Once Fe is accumulated in the liver, it could interrupt with role of insulin and also constrain glucose production( Reference Niederau, Berger and Stremmel51 ). Increased Fe accumulation might lead to IR by constraining glucose uptake in different tissues( Reference Fernández-Real, López-Bermejo and Ricart49 ). Clinical studies have shown no significant effect of Fe supplementation on CRP levels( Reference Wagener, Volk and Willis42 ).
Additionally, the effect of red meat on uric acid levels could constitute another pathway linking this dietary component with glucose/insulin homoeostasis dysregulation( Reference Dehghan, Van Hoek and Sijbrands52 ). It has been reported that uric acid could play a role in oxidative stress( Reference Butler, Morris and Belch53 ) and inflammatory factors( Reference Butler, Morris and Belch53 ), which are both linked to the progress of unfavourable glucose/insulin homoeostasis( Reference Dehghan, Van Hoek and Sijbrands52 , Reference Butler, Morris and Belch53 ). Further, an experiment in animal models reported that fructose-induced hyperuricaemia plays a pathogenic role in the development of cardio-metabolic risk factors( Reference Nakagawa, Hu and Zharikov54 ).
A prolonged intake of SFA is correlated with the MetS and is known to contribute to weight gain and inflammation if consumed in excess( Reference Lundman, Boquist and Samnegard55 , Reference Johnson, Mander and Jones56 ). In particular, SFA are known to contribute to increases in influence of white adipose tissue increasing inflammatory response( Reference Kien, Bunn and Ugrasbul57 – Reference Vergnaud, Norat and Romaguera59 ). Hence, we hypothesised that a high energy diet which contains excess red meat and is high in SFA may be associated with weight gain and thus increased adiposity and subsequently contribute to developing unfavourable glucose/insulin homoeostasis enhancing low-grade inflammation, which is strongly linked to the pathogenesis of CVD and other non-communicable diseases. The role of dietary factors, such as excess refined sugar or saturated fat intake, in triggering low-grade chronic inflammatory response has received further scientific support recently, reiterating the link with age-related chronic conditions. These inflammatory responses are thought to interact with the ageing process and, if persisting, could play a key role in the pathogenic mechanisms, leading to the onset of chronic metabolic and CVD( Reference Hotamisligil60 ).
The main strength of the present study is the investigation of the mediation effects using various markers of adiposity including not only BMI and WC, which are the markers of general and abdominal obesity, but also apVAT, VAI and LAP. Moreover, we have used a randomly selected, large and representative sample, and our results can be extrapolated to the general population. Lastly, the analyses included extensive adjustment for potential confounders which reduces the chance of the residual confounders.
The study has limitations. The study is cross-sectional and focused on adults only. Cohort studies may better address the causal relation between red meat intake and relevant health outcomes but may not be feasible because of the nature of the exposure and ethical issues. Consumption reported associations do not necessarily mean causation. Consumption of red meat in different life stages (childhood, adolescents adulthood and during ageing) has been previously shown to affect risk estimates( Reference Farvid, Cho and Chen61 ). Furthermore, different patterns of exposure over time could also affect the results. Cooking methods and variations in animal farming and meat preparation can alter the quality as well as the health effects of red meats( Reference Binková, Šmerhovský and Strejc62 , Reference Lakshmi, Schut and Zenser63 ). This introduces additional residual confounding into the analyses( Reference Bogen and Keating64 – Reference Kouvari, Tyrovolas and Panagiotakos66 ). This information was not available and therefore could not be controlled. The mediating effect of WC may be affected by BMI, or vice versa, because of the high colinearity between these two variables. This issue could be resolved by adding BMI and WC simultaneously to the mediation model( Reference Wagener, Volk and Willis42 ). However, this approach was not possible in our analyses because of the complex survey design of the present study. Hence, in an attempt to overcome this limitation, we have added other validated adiposity factors to the models such as apVAT, LAP and VAI, which provide an independent evaluation of the mediating role of central adiposity in red meat intake and markers of diabetes risk. Lastly, although BMI and WC are regularly applied to determine adiposity, these indicators are still imprecise and can lead to bias in determining obesity. For instance, BMI, is usually limited compared with direct measures of obesity due to the fact that it does not consider age, sex, bone structure, fat distribution or muscle mass info( Reference Rothman67 ). However, we also used other markers of adiposity (apVAT and LAP) which are sensitive to the age and sex.
Conclusion
Our study finding suggests that high RMC could negatively affect glucose/insulin homoeostasis and inflammatory profile, via mechanisms involving central fat accumulation. Future research is warranted to explore the effect of reducing red meat intake on glucose/insulin homoeostasis, which in turn could inform dietary strategies to reduce the risk of type 2 diabetes mellitus.
Acknowledgements
The corresponding author is responsible for submitting a competing financial interest statement on behalf of all authors of the paper. The authors declare no competing financial interests.
M. M. was involved in study conception and design, analysis and interpretation of data and drafting of manuscript. A. P. K., E. S. G. and M. S. were involved in critical revision and interpretation of data.
The authors have no conflicts of interest to declare.
For the data collection and physical examination of the NHANES, informed consent was obtained from all adult participants, and the National Centre for Health Statistics Research Ethics Review Board approved the protocol.
For the data collection and physical examination of the NHANES, informed consent (publication) was obtained from all adult participants, and the National Centre for Health Statistics Research Ethics Review Board approved the protocol.
All the data are from a public access database.
Supplementary material
For supplementary material/s referred to in this article, please visit https://doi.org/10.1017/S0007114519002149