Legumes and nuts intake in relation to metabolic health status, serum brain derived neurotrophic factor and adropin levels in adults

There is controversial evidence about the relationship between consumption of legumes and nuts with metabolic disturbances. The present study was undertaken to explore the association of legumes and nuts intake with metabolic health status among Iranian adults. This cross-sectional study was conducted on 527 adults (45.7% female, aged 20–65 years) chosen through a multistage cluster random-sampling approach. Dietary intakes of individuals were assessed using a validated food frequency questionnaire. Fasting blood samples were gathered to evaluate biochemical parameters. Metabolic health status of subjects was determined according to the criteria defined by Wildman. Data of covariates were collected using pre-tested procedures. The overall prevalence of metabolic unhealthy (MU) phenotype was 42.5%. After controlling all confounders, participants with highest intake of legumes and nuts had lower odds of MU status, compared with the lowest intake (OR 0.35; 95% CI 0.18–0.71). This association was stronger in normal-weight rather than overweight/obese adults and also in women rather than men. Higher consumption of legumes and nuts was additionally related to decreased odds of hyperglycemia, hypertriglyceridemia, and hypertension. A marginally inverse association was observed between legumes and nuts intake with low brain-derived neurotrophic factor (BDNF) levels, in fully-adjusted model (ORT3 vs. T1 0.50; 95% CI 0.25–1.01). Each tertile increase in legumes and nuts intake was marginally related to higher adropin levels (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β = 4.06; P = 0.07). In conclusion, this study demonstrated that higher intake of legumes and nuts is associated with lower chance of MU both in normal weight and overweight/obese adults. The association may be facilitated through serum BDNF and adropin.

Obesity, currently affecting about one-third of the world's population, is considered as a major public health threat, due to its effect on risk of chronic diseases such as type 2 diabetes mellitus, cardiovascular diseases (CVDs), and some types of cancers 1-3 .However, recently the concept of metabolic health status has challenged the role of obesity in determining risk of chronic diseases.To be more exact, a subgroup of obese individuals with a favorable metabolic profile who are known as metabolically healthy obese (MHO) is less likely to have chronic diseases compared to metabolically unhealthy obese (MUO) individuals 4 .Additionally, a high risk of cardiometabolic disorders has been reported in normal-weight individuals with a metabolically unhealthy (MU) profile 5,6 .Therefore, it seems that metabolic status could be a better predictor for risk of chronic diseases than obesity, and maintaining a healthy metabolic status could be considered a beneficial strategy for reducing risk of chronic diseases.During recent decades, emerging evidence has pointed to the metabolic role of hormones such as adropin and brain-derived neurotrophic factor (BDNF).Findings of experimental investigations indicated that these secreted molecules could control metabolism of glucose and lipid and regulate energy hemostasis 7,8 .Decreased serum levels of adropin and BDNF have been also reported in people with metabolic disorders such as obesity, diabetes mellitus, and CVDs [9][10][11] .Some reports have also indicated that modifiable factors such as dietary intakes and physical activity could affect adropin and BDNF levels [12][13][14][15][16] .It is postulated that adropin and BDNF play key roles in determining metabolic health status.Therefore, altering their serum levels is probably one of the mechanisms through which behavioral interventions improve metabolic health status.
Previous studies have suggested lifestyle modification as a beneficial approach for improving metabolic parameters.During recent decades, many epidemiological and clinical studies have highlighted the role of diet, among other lifestyle factors, on metabolic status 17 .A great body of evidence from observational and interventional studies has demonstrated the beneficial role of healthy dietary patterns, such as Mediterranean diet, dietary approach to stop hypertension (DASH) diet and plant-based diets on metabolic status [18][19][20] .These dietary patterns are rich in fruits, vegetables, whole grains, legumes, and nuts.Among plant-based foods, legumes, and nuts may affect the metabolic parameters potentially, due to their favorable content of nutrients and bioactive compounds such as vegetable protein, fiber, unsaturated fatty acids, vitamins, minerals, phytosterols, and phytochemicals 21,22 .
The relationship between dietary intake of legumes and nuts with metabolic syndrome (MetS) or its components has been investigated in previous studies.However, the results of these studies are contradictory.Some epidemiological evidences suggest that dietary intake of legumes and nuts is associated with a reduced risk of MetS or its components [23][24][25][26][27] , while other studies have shown conflicting results [28][29][30][31][32] .A recently published investigation has evaluated the link of nuts and legumes intake with metabolic health status in Iranian adolescents 33 , but no previous study has investigated this relationship among adults, especially in Middle Eastern populations.Therefore, the purpose of this study was to investigate the association of nuts and legumes intake with metabolic health status in Iranian men and women adults, considering the potential underlying role of BDNF and adropin.

Study design and population
This cross-sectional study was performed in 2021, on a sample of Iranian adults living in Isfahan, Iran.According to a multistage cluster random-sampling approach, all adults working as teachers, school managers, employees, assistants or crews in several schools of different educational districts of Isfahan were selected.Based on a previous published study 34 , 49.4% of Iranian adults might have metabolic disorders; so, considering type 1 error of 0.05 (confidence interval (CI) of 0.95), precision (d) of 4.5% and power of 80%, 474 subjects were minimally required for this study.Nevertheless, due to high prevalence of covid-19 pandemic and its potential impact on data collection, 600 individuals were invited to participate in our investigation.Response rate was 90.5%.Individuals with the following characteristics were not included in our study: (1) having a history of type 1 diabetes, cardiovascular diseases, stroke and cancer; (2) being pregnant or lactating; (3) following a special diet.Exclusion criteria were: (1) having left more than 70 items on the food frequency questionnaire (FFQ) unfilled; (2) reporting energy intake outside of 800-4200 kcal/day; (3) refusing blood draw.Finally, a total of 527 adults aged 20 to 65 years old were eligible to be included in this analysis.Each participant signed a written informed permission.Local Ethics Committee of Isfahan University of Medical Sciences approved the protocol of this study (no.3402104).

Assessment of dietary intakes
We evaluated dietary intakes of participants using a Willett-format semi-quantitative 168-item FFQ.Validity of this questionnaire has been approved through a prior validation research which found reasonable correlations between the food intakes determined by this questionnaire and those acquired from several 24-h dietary recalls 35 .A comparison of nutrient intakes derived from this FFQ on two occasions, 1 year apart, indicated reliability of this tool as well 35 .All Participants were guided by a registered dietitian to fill out the FFQ by describing the frequency and amount of their eaten foods during the last year.Afterward, the portions of ingested items were changed to g/day using household measurements 36 .Lastly, data of all food items were input into the Nutritionist IV program to calculate daily energy and nutritional intakes.

Assessment of anthropometric indices and blood pressure
Weight, height and waist circumference (WC) of participants were assessed by two skilled dietitians, with little clothing and no shoes.The body composition analyzer (Tanita MC-780MA, Tokyo, Japan) was utilized to determine weight (to the nearest 0.1 kg).Height was measured using a tape measure fixed on the wall (to the nearest 0.1 cm).Then, body mass index (BMI) was computed by dividing weight (kg) by height squared (m 2 ).WC was evaluated after a normal breath with no outside pressure on the body and by applying an unstretched flexible tape measure (to the nearest 0.1 cm).Blood pressure (BP) of each person was taken after sitting for five minutes with an empty bladder and no prior exercise.A digital sphygmomanometer (OMRON, M3, HEM-7154-E, Japan) with an accuracy of 0.5 mmHg was used to assess BP.The measurements were conducted twice, with intervals of 5-10 min, and average of the measurements was reported as the final BP.

Assessment of other variables
Data on sex, age, marital status, education and smoking were collected through a self-reported questionnaire.Furthermore, a validated questionnaire was applied to evaluate socioeconomic status (SES) of participants 39 .This questionnaire assessed individuals in terms of the number of family members, homeownership, type of car, having laptops/computers and traveling in the year.Physical activity (PA) was also measured using a validated International Physical Activity Questionnaire-short form (IPAQ-SF) that evaluates three categories of activity including walking, moderate-intensity activities, and vigorous-intensity activities 40 .The data from this questionnaire were converted to Metabolic Equivalent minutes per week (MET.min/week)and participants were categorized as inactive (< 600 MET.min/week), minimally active (≥ 600 to < 3000 MET.min/week), or active (≥ 3000 MET.min/week).

Statistical analysis
Energy-adjusted intake of legumes and nuts was computed based on residual method.Then, participants were classified according to tertiles of energy-adjusted legumes and nuts intake.Continues and categorical variables were reported as mean ± SD/SE and number (percentage), respectively.We used one-way ANOVA and chi-square test to compare characteristics of subjects across tertiles of energy-adjusted legumes and nuts intake.Dietary intakes of individuals were also assessed through ANCOVA by adjusting sex, age and energy intake.Possible associations between intake of legumes and nuts with metabolic unhealthy (MU) status and its components were determined using binary logistic regression.Confounding effects of potential variables (including age, sex and energy intake in the first model, education, marital status, smoking status, SES and PA in the second model, intake of fruits, vegetables, dairy, whole and refined grains in the third model, and BMI in the last model) were controlled in the analyses 41,42 .All odds ratios (ORs) were calculated by considering the first tertile of legumes and nuts as the reference.Tertiles of energy-adjusted legumes and nuts intake were regarded as an ordinal variable in logistic regression models to evaluate trend of ORs across tertiles of legumes and nuts intake.Stratified analyses were also conducted by BMI categories (normal-weight vs. overweight/obese) and sex (women vs. men).
Multivariable-adjusted odds of low BDNF values (< 0.47 ng/mL) in tertiles of energy-adjusted legumes and nuts intake were determined by considering the effects of age, sex, PA, high BP, high TG and high FBG.In addition, linear regression analysis was applied to assess adropin levels in tertiles of energy-adjusted legumes and nuts intake with adjusting age, sex, energy intake, PA and BMI as covariates.SPSS software version 26 (IBM, Chicago, IL) was utilized for all analyses and P-values < 0.05 were regarded to be statistically significant.

Ethical approval and consent to participate
The study procedure was performed according to declaration of Helsinki and STROBE checklist.All participants provided informed written consent.The study protocol was approved by the local Ethics Committee of Isfahan University of Medical Sciences.

Consent to participate
Informed consent was obtained from all participants involved in the study.

Results
This study was conducted on 527 adults (  1. Individuals in the top tertile of legumes and nuts intake were more likely to be male, had a lower body weight and also lower prevalence of hypertriglyceridemia, compared to those at the bottom tertile.However, no substantial differences were observed in other demographic and cardiometabolic variables including age, BMI, WC, education, marital status, smoking, SES, PA, high BP, high FBG, IR (high HOMA-IR score), high hs-CRP, low HDL-c, low BDNF values and adropin levels.
Dietary intakes of subjects across energy-adjusted tertiles of legumes and nuts intake are illustrated in Table 2. Participants with higher consumption of legumes and nuts had also higher intakes of proteins, dietary fiber, folate, magnesium and potassium, and lower intakes of total energy and refined grains.While, consumptions of fruits, vegetables, meats, fish, whole grains, dairy, carbohydrates, fats, vitamin C and calcium did not differ significantly across tertiles of legumes and nuts intake.
Figure 1 represents the prevalence of MU status across energy-adjusted tertiles of legumes and nuts intake.Prevalence of MU phenotype was 50.9% in first tertile, 40.3% in the second tertile and 36.4% in the third tertile of legumes and nuts intake.This prevalence was significantly different across tertiles of legumes and nuts intake (P value = 0.02).
Crude and multivariable-adjusted ORs for MU phenotype across energy-adjusted tertiles of legumes and nuts intake are reported in Table 3.In comparison with the first tertile of legumes and nuts intake, individuals at the last tertile had a decreased odds of MU status in crude model (OR T3 vs.T1 0.55; 95% CI 0.36-0.85).After controlling potential variables, this association was strengthened (OR T3 vs.T1 0.35; 95% CI 0.18-0.71).A significant decreasing trend was also observed for MU phenotype across tertiles of legumes and nuts intake in all models (P trend = 0.01).Each tertile increase in intake of legumes and nuts was significantly associated with reduced odds of MU status in crude and maximally-adjusted models as well.Stratified analysis by BMI revealed an inverse relationship between legumes and nuts intake and likelihood of MUNW profile among normal-weight participants, in both crude (OR T3 vs.T1 0.35; 95% CI 0.14-0.87)and fully-adjusted (OR T3 vs.T1 0.15; 95% CI 0.03-0.85)models.Among participants with overweight/obesity, legumes and nuts intake was marginally associated with reduced odds of MUOW profile in crude model (OR T3 vs.T1 0.60; 95% CI 0.36-1.00).This association was strengthened after considering confounders (OR T3 vs.T1 0.29; 95% CI 0.12-0.69).
Table 4 presents the association between MU status and intake of legumes and nuts, stratified by sex.Among women, an inverse association was observed between dietary intake of legumes and nuts with MU status in both crude (OR T3 vs.T1 0.46; 95% CI 0.23-0.94)and fully-adjusted (OR T3 vs.T1 0.28; 95% CI 0.09-0.091)models.Legumes and nuts intake of men was not substantially related to MU phenotype in crude model (OR T3 vs.T1 0.68; 95% CI 0.40-1.19).However, by adjusting potential confounders this association became marginally significant (OR T3 vs.T1 0.40; 95% CI 0.16-1.00).
The average serum levels of BDNF and adropin among participants were 1.25 ng/mL and 56.59 pg/mL, respectively.Multivariate-adjusted ORs for low BDNF levels across tertiles of legumes and nuts intake are depicted in Fig. 2. Compared to the reference group, higher intake of legumes and nuts was marginally related to decreased likelihood of low BDNF levels, in both crude (OR T3 vs.T1 0.56; 95% CI 0.28-1.10)and multivariable-adjusted (OR T3 vs.T1 0.50; 95% CI 0.25-1.01)models.As shown in Fig. 3, linear regression analysis revealed that each tertile increment in legumes and nuts intake was substantially associated with an increase of 4.78 pg/mL in levels of adropin in crude model (95% CI 0.47-9.09;P value = 0.03).After adjusting all covariates, a marginally significant direct association was observed between each tertile increase in legumes and nuts intake with adropin levels ( β = 4.06; 95% CI − 0.32, 8.44; P value = 0.07).No significant differences were found in circulating BDNF and adropin levels between individuals with metabolically healthy and unhealthy phenotypes, even after stratifying analysis by sex or age categories (Supplemental Figs. 1, 2).Table 1.Demographic and cardiometabolic features of participants across energy-adjusted tertiles of legumes and nuts intake (n = 527).Quantitative variables: means ± standard deviation (SD).Qualitative variables: frequency (percentage). 1Legumes and nuts intake was adjusted for energy intake based on residual method. 2 Resulted from ANOVA for quantitative variables and chi-square test for categorical variables.

Discussion
The findings of the present cross-sectional study showed that more than 40% of the Iranian population had an MU phenotype with a significantly higher prevalence among obese/overweight than normal-weight individuals (79.5% vs. 20.5%).Higher consumption of legumes and nuts was associated with lower odds of MU phenotype both in normal-weight and obese/overweight individuals.This association was more prominent among females.A significant inverse relationship was also found between legumes and nuts intake with some components  Table 3. Multivariable-adjusted odds ratio for metabolic unhealthy status across energy-adjusted tertiles of legumes and nuts intake.All values are odds ratios and 95% confidence intervals.Model 1: Adjusted for age, sex, energy intake.Model 2: More adjustments for physical activity, socioeconomic status, education, marital status, smoking status.Model 3: Further adjustments for dietary intake of fruits, vegetables, dairy, whole and refined grains.Model 4: More adjustment for BMI.P trend was obtained by the use of tertiles of legumes and nuts intake as an ordinal variable in the model. 1 Legumes and nuts intake was adjusted for total energy intake based on residual method.Table 4. Multivariable-adjusted odds ratio for metabolic unhealthy status across energy-adjusted tertiles of legumes and nuts intake stratified by sex.All values are odds ratios and 95% confidence intervals.Model 1: Adjusted for age, and energy intake.Model 2: More adjustments for physical activity, socioeconomic status, education, marital status, smoking status.Model 3: Further adjustments for dietary intake of fruits, vegetables, dairy, whole and refined grains.Model 4: More adjustment for BMI.P trend was obtained by the use of tertiles of legumes and nuts intake as an ordinal variable in the model. 1 Legumes and nuts intake was adjusted for total energy intake based on residual method. of metabolic health status including hyperglycemia, hypertriglyceridemia, and hypertension.These beneficial associations might be facilitated through higher BDNF and adropin levels.Non-communicable diseases (NCDs), mainly CVDs, cancers, respiratory diseases, and diabetes, account for over 70% of deaths worldwide 43 .Among the modifiable lifestyle behaviors increasing the risk of NCDs, diet plays a key role through its contribution in metabolic disturbances 44 .Therefore, dietary intake assessments and implementing early clinical interventions could be beneficial approaches for reducing risk of metabolic disturbances and related NCDs.Findings of this study revealed that daily consumption of more than 54 g of legumes and nuts could reduce MU odds.Therefore, providing nutritional education and interventions to consume a Table 5. Multivariable-adjusted odds ratio for metabolic components across tertiles of legumes and nuts intake.All values are odds ratios and 95% confidence intervals.P trend was obtained by the use of tertiles of legumes and nuts intake as an ordinal variable in the model.FBG fasting blood glucose, TG triglycerides, HDL high density lipoprotein, BP blood pressure, HOMA-IR homeostasis model assessment insulin resistance, hs-CRP high sensitive C-reactive protein. 1 Legumes and nuts intake was adjusted for total energy intake based on residual method. 2 Adjusted for age, sex, energy intake, physical activity, socioeconomic status, education, marital status, smoking status, dietary intake of fruits, vegetables, dairy, whole and refined grains and BMI. 3 HDL-c < 40 mg/dL in men, and < 50 mg/dL in women.healthy diet containing appropriate amounts of legumes and nuts is recommended as a beneficial approach to reduce risk of MU status.
To the best of our knowledge, this report is the first observational study investigating the relationship between the intake of legumes and nuts with metabolic health status.However, legume and nut intake in relation to MetS or its components has been investigated previously and revealed controversial results.Contrary to our results, a cross-sectional study involving 420 Iranian female nurses has failed to find any association between the consumption of legumes and nuts with MetS or its components 28 .Additionally, a systematic review and meta-analysis of 7 observational studies containing 56,028 participants demonstrated no significant association between legume intake and odds of MetS 41 .However, the beneficial effect of legume intake on total cholesterol and low-density lipoprotein cholesterol (LDL-c) as well as BP in obese and overweight individuals has been reported by two other studies 45,46 .A meta-analysis of 11 observational studies showed an inverse significant association between nut intake and risk of MetS 42 .There are also other reports on the protective role of nut consumption on cardiometabolic disorders.Nut intake has been shown to have a favorable effect on HOMA-IR and fasting insulin in a meta-analysis of 40 randomized controlled trials 47 .In addition, the beneficial effect of nut consumption on hypertension has been reported in other studies 48 .However, in a meta-analysis of 6 randomized controlled trials, no significant association was found between total nut intake and lipid profile 49 .These discrepant findings can be explained by different genetic backgrounds and dietary habits of the studied population, tools used to assess dietary intakes, and controlling for various confounders.
Many NCDs in adulthood, especially endocrine dysfunctions, have roots in health status, lifestyle and dietary intakes in childhood.According to a cohort survey by Liang et al., children with obesity had higher risk of diabetes in adulthood 50 .Another study also revealed that childhood risk factors including unhealthy diets had long-term effects on metabolic health and CVD risk in adulthood 51 .Therefore, maintaining the health status in early life would be of great importance in every society.In our previous study, we discovered that higher intake of legumes and nuts was associated with lower odds of metabolic unhealthy among Iranian adolescents aged 12 to 18 years 33 .The same findings were obtained in the current investigation in adults.Although there were no sufficient data regarding the metabolic status of the present study participants in their childhood, it could be mentioned that the current metabolic status of these individuals might be affected by their childhood lifestyle.Therefore, it would be more valuable to modify food habits and intakes from childhood and adolescence in order to prevent the occurrence of many diseases in adulthood.
Our results showed that the inverse association between legumes and nuts intake with MU phenotype was more prominent in females.This finding might be explained by healthier dietary knowledge and behavior among women than men 52 .Also, better metabolic profiles among women due to the differences in substrate use, accumulation and metabolism in important metabolic organs including liver, skeletal muscle and fat tissue may explain the observed finding 53 .Furthermore, premenopausal women have reduced odds of metabolic disturbances which highlights the key role of sex-hormones in metabolic health 53 .In addition, the observed association was stronger among normal-weight rather than over weight/obese subjects.Different physiological responses to environmental factors, such as diet, and differences in accuracy of reported food intakes might explain the observed findings.However, further studies are required to confirm these hypotheses.
We found a marginally inverse association between low BDNF levels with legumes and nuts intake.A positive linear association was also observed between levels of adropin and legumes and nuts intake.So far, several studies have examined the association between some dietary factors with BDNF and adropin levels in both animals and humans 12,13,16,54 .A clinical trial conducted on Spanish adults found that the Mediterranean diet supplemented by nuts was related to reduced odds of low BDNF levels 55 .There were no other related investigations that examined the association of legumes or nuts intake with BDNF and adropin levels.The current epidemiologic study www.nature.com/scientificreports/revealed no substantial differences in BDNF or adropin levels between metabolically healthy and unhealthy groups.Although some previous studies have documented that higher concentrations of BDNF and adropin were related to decreased odds of metabolic disorders through regulating energy hemostasis and its related signaling pathways [7][8][9] , there are insufficient data regarding this association and the possible mechanisms.Our findings could broaden insights into further related studies.The inverse relationship between the intake of legumes and nuts with MU phenotype can be explained by several feasible mechanisms.Legumes and nuts are rich in fiber, minerals, and bioactive compounds which might have beneficial effects on metabolic health status.Magnesium content of legumes and nuts plays an important role in decreasing inflammation and insulin resistance 56,57 .Additionally, a large body of evidence indicates that dietary fiber contributes in improving metabolic health status components such as hypertension, dyslipidemia, insulin sensitivity, and inflammatory biomarkers levels 58 .Another explanation for the beneficial influence on metabolic health status may be the antioxidant and anti-inflammatory properties of bioactive phytochemicals such as polyphenols which have modulatory effects on metabolic processes 59 .Moreover, low glycemic index of legumes and nuts reduces the risk of insulin resistance and improves metabolic health status 60 .
The current study has several limitations that should be acknowledged.The study does not provide evidence for a causal relationship between the consumption of legumes and nuts with MU, due to its cross-sectional design.Self-reported dietary intakes might be subject to recall bias and misclassification, despite using a validated FFQ for assessment of dietary intakes.In addition, due to low amounts of nuts intake in our study population, no separate analysis was performed for nuts and legumes intake.On the other hand, to the best of our knowledge, this is the first study that investigates the links between the intake of legumes and nuts with metabolic health status in a somehow representative sample of Iranian adults.Furthermore, a relatively comprehensive definition of metabolic health status (Wildman et al. method) was used that includes an inflammatory index (hs-CRP) in addition to HOMA-IR and traditional cardiometabolic risk factors.The confounding role of several variables was also taken into account in the statistical analyses.The final strength of the study was the assessment of serum levels of adropin and BDNF which rarely have been investigated in epidemiological studies of nutrition.
In summary, a significant relationship was found between consumption of legumes and nuts with decreased odds of MUNW and MUOW phenotypes, particularly among females.Higher consumption of legumes and nuts was also inversely associated with hyperglycemia, hypertriglyceridemia, and hypertension.The association might be facilitated through BDNF and adropin.Further population-based prospective studies are warranted to confirm these findings. https://doi.org/10.1038/s41598-023-43855-8 https://doi.org/10.1038/s41598-023-43855-8

Figure 1 .
Figure 1.Prevalence of metabolically unhealthy (MU) status across energy-adjusted tertiles of legumes and nuts intake in the study population.

Figure 2 .
Figure 2.Multivariable-adjusted odds ratio and 95% confidence intervals for low BDNF levels across tertiles of legumes and nuts intake.Model 1: Adjusted for age and sex; Model 2: More adjustments for physical activity, history of high blood pressure, high triglyceride and high fasting blood glucose.P trend was obtained by the use of tertiles of legumes and nuts intake as an ordinal variable in the model.

Figure 3 .
Figure 3. Linear association between tertiles of legumes and nuts intake with adropin levels.All values are regression coefficients and 95% confidence intervals.Tertiles of legumes and nuts intake were considered as an ordinal variable in linear regression analysis.Model 1: Adjusted for age, sex, and energy intake; Model 2: More adjustments for physical activity and BMI.

Table 2 .
Multivariable -adjusted intakes of selected food groups and nutrients of study participants across energy-adjusted tertiles of legumes and nuts intake (n = 527).Values are means ± standard error (SE).Energy intake and macronutrients were adjusted for age and sex; all other values were adjusted for age, sex and energy intake. 1 Legumes and nuts intake was adjusted for energy intake based on residual method. 2 P-value obtained from ANCOVA test for adjustment of energy intake.