Introduction

Overweight and obesity, defined as abnormal or excessive fat accumulation in the body, in childhood and adolescence has adverse health effects throughout the life1. In 2016, estimates showed that more than 330 million children and adolescents (aged 5–19 years) were overweight or obese2. Global burden of disease (GBD) data reported that by the year 2025, 268 and 124 million children and adolescents will be respectively overweight and obese3. According to a systematic review and meta-analysis published in 2020, 11% of Iranian students were generally obese (based on their body mass index (BMI))4. The existence of these conditions in early life is associated with psychosocial consequences, less educational attainment, and cardiometabolic disorders such as high blood pressure (BP), dyslipidemia and insulin resistance5,6,7. Also, childhood obesity can increase burden on health systems and reduce economic productivity later on8.

Recent attention has focused on two phenotypes for overweight and obesity including metabolically healthy overweight/obese (MHO) and metabolically unhealthy overweight/obese (MUO). As compared to their MUO counterparts, MHO subjects have better lipid profiles, insulin sensitivity, as well as lower cardiovascular disease (CVD) risk factors9,10. It would be possible to shift between metabolic health statuses11. Evaluation of metabolic health in the childhood period may be important to prevent adverse cardiometabolic diseases12. Physical activity, sleep, smoking and dietary intakes are factors that can affect metabolic health status13.

A previous study showed that the Dietary Approaches to Stop Hypertension (DASH) diet may improve insulin sensitivity in overweight children and adolescents14. Also, healthy eating patterns may have a positive effect on metabolically health status15. It is well-known that diet can have a significant effect on acid–base balance. When a diet is rich in animal products, cheese, cereals, and rice and low in fruits and vegetables, can result in endogenous acid production16. Acidosis has adverse health effects such as obesity, overweight, increased BP, and higher triglyceride (TG) levels17,18. Two scores are commonly used to measure dietary acid load (DAL) including the potential renal acid load (PRAL) score and net endogenous acid production (NEAP) score19. Women with higher NEAP scores had higher waist circumference (WC), serum TG, and weight in a cross-sectional study in Iran20. Furthermore, a study on adult Japanese workers found that higher PRAL and NEAP scores were associated with insulin resistance21. Another investigation conducted among young Japanese women demonstrated a link between a more acidic dietary load and the adverse profile of several cardiometabolic risk factors22. In case of children and adolescents, Aslani et al. illustrated a direct association between diet-induced acid load and neck circumference (NC) as well as an inverse association between DAL indices and parental BMI23. However, a recent study could not find a significant association between DAL and metabolic syndrome among children and adolescents24. Since there is no study on the association between DAL and metabolic health status in adolescents, the current study investigated the association of DAL with metabolic health status in Iranian overweight/obese adolescents.

Methods

Study design and population

This cross-sectional study included 102 girls and 101 boys (12–< 18 years old) who were selected by a multistage cluster random sampling method. Participants were randomly selected from 6 different districts of the city of Isfahan, Iran, using a stratified, multi-stage cluster sampling design. In total, 16 schools were chosen. Body weight and height of all students were measured and BMI was calculated based on weight (kg)/height2 (m). Then, adolescents were classified into 3 groups: normal-weight, overweight, and obese25. Overweight and obese students with different socioeconomic status (SES) were considered to participate in the study. The students who had the following conditions were not included in this study: having genetic or endocrine disorders (such as type 1 diabetes mellitus, hypothyroidism, and Cushing's syndrome), being on a weight-loss diet, taking vitamin and mineral supplements or medications that might affect body weight, hypertension, blood glucose or lipid profile. All participants and their parents signed the written informed consent. The study protocol was approved by the local Ethics Committee of Isfahan University of Medical Sciences.

Assessment of dietary intakes

Dietary intake was assessed using a validated 147-item semi-quantitative food frequency questionnaire (FFQ)26. According to this questionnaire, dietary intake in the previous year was calculated by determining the frequency of consumption of each food item on a daily, weekly, or monthly basis. The amount of food consumed was evaluated by asking a standard common portion size. We converted the portion sizes of consumed foods to grams per day using household measures. Finally, the grams of food intakes were entered into Nutritionist IV software to derive nutrient intake data. Nutritionist IV software was based on USDA food composition data; some Iranian foods were also incorporated into it.

Calculation of DAL indices

Urinary net acid excretion depends on dietary intake, and direct measurement is difficult. Therefore, two dietary indices (PRAL and NEAP) have been defined to calculate body acid production or DAL. In a previous study by Remer et al., the validity of both PRAL and NEAP scores compared to a 24-h urinary acid load was confirmed27. The PRAL score was calculated based on the following formula27:

$$ {\text{PRAL}}\,\left( {\text{mEq/day}} \right)\, = \,\left( {0.{4888}\, \times \,{\text{protein}}\,{\text{intake}}\,\left( {\text{g/day}} \right)} \right)\, + \,\left( {0.0{366}\, \times \,{\text{phosphorus}}\,\left( {\text{mg/day}} \right)} \right)\,{-}\,\left( {0.0{2}0{5}\, \times \,{\text{potassium}}\,\left( {\text{mg/day}} \right)} \right)\,{-}\,\left( {0.0{125}\, \times \,{\text{calcium}}\,\left( {\text{mg/day}} \right)} \right){-}\left( {0.0{263}\, \times \,{\text{magnesium}}\,\left( {\text{mg/day}} \right)} \right). $$

In addition, we estimated the NEAP score as follows28:

$$ {\text{NEAP}}\,\left( {\text{mEq/day}} \right) = \left[ {{54}.{5}\, \times \,{\text{protein}}\,{\text{intake}}\,\left( {\text{g/day}} \right)\, \div \,{\text{potassium}}\,{\text{intake}}\,\left( {\text{mEq/day}} \right)} \right]\, - \,{1}0.{2}. $$

Assessment of anthropometric indices and cardiometabolic risk factors

Two trained nutritionists measured weight using a calibrated electronic scale (Seca Instruments, Germany) in minimal clothing and without shoes to the nearest 0.1 kg. Also, standing height was measured without shoes with a stadiometer to the nearest 0.1 cm. Then, body mass index was calculated as kg/m2. Based on the sex- and age-specific World Health Organization (WHO) cut-off points, participants were classified as overweight or obese adolescents25. Waist circumference was recorded twice at the midway between the lowest rib and the superior border of the iliac crest by the use of an unstretched flexible anthropometric tape at 0.1 cm precision, after a normal expiration and without any pressure on the body surface. The mean of the two measurements for each student was considered as WC. Systolic blood pressure (SBP) and diastolic blood pressure (DBP) were measured twice at the right arm using a mercury sphygmomanometer with a suitable cuff size after a rest period of 15 min. We considered the average of the two measurements in the analysis. Venous blood samples after 12 h of fasting were collected to determine biochemical values. Fasting blood glucose (FBG), high-density lipoprotein cholesterol (HDL-c), triglycerides (TG), and insulin concentrations were measured with standard methods. Homeostasis Model Assessment Insulin Resistance (HOMA-IR) was calculated to estimate Insulin resistance (IR) by using the following formula: HOMA-IR = (insulin (µUI/mL) × glucose (mg/dL))/405.

Assessment of metabolic status

Two methods were used to classify participants into MHO or MUO individuals. Students with two or more of the four following risk factors were defined as MUO subjects, based on the modified International Diabetes Federation (IDF) criteria29: (1) increased fasting blood glucose (≥ 100 mg/dL), (2) increased blood pressure (≥ 130/85 mmHg), (3) increased triglycerides (≥ 150 mg/dL) and (4) decreased HDL-c (< 40 mg/dL for the age of < 16 y, and < 50 mg/dL in girls/ < 40 mg/dL in boys for the age of ≥ 16 y). MHO adolescents were those with one or no risk factors. In the second method, presence of insulin resistance (defined by HOMA-IR score) was added to the IDF criteria that were used in the first definition30. Students with HOMA-IR score ≥ 3.16 that had two or more above-mentioned metabolic risk factors were considered as MUO individuals and those with HOMA-IR < 3.16 were classified considered as MHO adolescents. Based on some previous studies on children and adolescents, the cut-off value for HOMA-IR was set at 3.1631,32,33.

Assessment of other variables

The Physical Activity Questionnaire for Adolescents (PAQ-A), which contains 9 questions on various aspects of physical activity over the past week, was used to determine the physical activity level of participants34. Items 1 to 8 of this questionnaire were scored from 1 to 5, with 1 representing the least physical activity and 5 representing the most. The ninth question examined the unusual activity of adolescents during the previous 7 days. Based on their total scores, adolescents were categorized as highly active (score ≥ 4), moderately active (4 < score ≤ 3), lowly active (3 < score ≤ 2), sedentary, or not having an orderly week of physical activity (score < 2). A validated demographic questionnaire was used by trained investigators to ascertain SES of students (29), based on the following items: parental education level, parental job, number of family members, having cars in the family, having personal room, having computers/laptops and taking trips in the year. Finally, a total score of SES (ranged between 4 and 27) was calculated. Information of sex, age, history of diseases, and use of medications and dietary supplements of participants was gathered through a questionnaire.

Statistical analysis

We presented continuous variables as mean ± SD/SE and qualitative variables as frequency (percentage). Participants were categorized into tertiles of PRAL and NEAP (based on their scores). We used one-way analysis of variance (ANOVA) and the χ2 test to compare quantitative and categorical variables between tertiles of PRAL and NEAP. Age-, sex- and energy-adjusted dietary macro and micronutrient intakes of participants were compared across tertiles of PRAL and NEAP using analysis of covariance (ANCOVA). Multivariable logistic regression was applied to identify the association between PRAL and NEAP with MUO status. The odds ratio (OR) and 95% confidence interval (CI) for MUO status were computed in crude and adjusted models. In the first model, we adjusted for sex, age, and energy intake. In the second model, further adjustments were done for physical activity levels and SES. In the last model, more adjustment for BMI was done. The first tertile of PRAL and NEAP was considered as the reference category, in all models. The overall trend of OR across increasing PRAL and NEAP tertiles was examined by considering tertiles of PRAL and NEAP as a continuous variable. Stratified analysis was done to obtain OR for MUO in different categories of BMI (overweight vs. obese individuals). SPSS software version 23 (IBM, Chicago, IL) was used for all analyses. P values < 0.05 (two-sided) were considered as 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.

Results

The present study included 101 (49.8%) boys and 102 (50.2%) girls with a mean age of 13.34 ± 1.38 and 14.61 ± 1.57 and BMI of 26.67 ± 2.82 and 28.03 ± 3.48, respectively. General characteristics and cardiometabolic risk factors of study participants across tertiles of PRAL and NEAP are presented in Table 1. Adolescents in the highest tertile of PRAL in comparison to those in the lowest tertile had higher FBG and TG (P < 0.05). Furthermore, subjects in the top category of PRAL had lower physical activity level and different SES status (P < 0.05). Adolescents in the highest tertile of NEAP had higher FBG than those in the lowest tertile and were different in terms of physical activity from those in the first category (P < 0.05). There were no significant differences in other general characteristics and cardiometabolic risk factors among tertiles of PRAL and NEAP.

Table 1 General characteristics and cardiometabolic factors of study participants across tertiles of PRAL and NEAP1.

Dietary intakes of participants across tertiles of PRAL and NEAP are presented in Table 2. Participants in the top tertile of PRAL, compared to the bottom tertile, had greater intake of energy, protein, thiamin, niacin, selenium and sodium and lower intake of fat, MUFA, vitamin C, vitamin A, riboflavin, vitamin B6, folate, vitamin B12, magnesium, total fiber, potassium, and calcium. Also, highest vs. lowest tertile of NEAP was associated with increased intake of carbohydrate, thiamin and niacin. However, adolescents in the first tertile of NEAP, compared to those in the third tertile, had higher intake of fat, cholesterol, SFA, MUFA, vitamin C, vitamin A, riboflavin, vitamin B6, folate, vitamin B12, magnesium, zinc, total fiber, potassium, calcium and phosphorus.

Table 2 Dietary intakes (energy and macro/micro nutrients) of study participants across tertiles of PRAL and NEAP 1.

The prevalence of MUO (based on IDF and IDF/HOMA-IR criteria) across tertiles of PRAL and NEAP are reported in Fig. 1. Based on the IDF definition for metabolic health status, the prevalence of MUO in the highest tertile of PRAL, compared to the bottom tertile, was significantly higher (50.0 vs. 26.9%, P = 0.02). Also, the higher scores of NEAP were associated with more prevalence of MUO (third tertile vs. first tertile of NEAP: 47.1 vs. 25.4%, P = 0.02). According to the second definition of MUO (IDF/HOMA-IR criteria), the prevalence of MUO in the top tertile of PRAL and NEAP vs. the bottom tertile was marginally significant (42.6 vs. 23.9%, P = 0.06 and 38.2 vs. 22.4%, P = 0.07, respectively).

Figure 1
figure 1

Prevalence of MUO across tertiles of PRAL and NEAP. (A) Based on IDF definition, (B) Based on IDF/HOMA-IR definition.

Crude and multivariate adjusted odds ratio and 95% CI for MUO phenotype across tertiles of PRAL and NEAP are shown in Table 3. In the crude model, adolescents in the top category of PRAL compared with the bottom level had 172% higher odds of MUO status, based on IDF criteria (OR 2.72; 95% CI 1.32–5.59). This association was significant after adjustment for age, sex, and energy intake (OR 2.42; 95% CI 1.13–5.15). But further adjustment for physical activity levels, SES, and BMI removed the significant association between PRAL categories and odds of MUO phenotype (fully-adjusted model: OR 1.55; 95% CI 0.65–3.66). The same findings were obtained for the association between NEAP tertiles with odds of MUO phenotype (crude model: OR 2.61; 95% CI 1.26–5.41; model 1: OR 2.52; 95% CI 1.17–5.41; and fully-adjusted model: OR 1.47; 95% CI 0.61–3.49). According to the IDF/HOMA-IR definition, the significant direct association between the highest category of PRAL and odds of MUO phenotype was observed only in the crude model (OR 2.37; 95% CI 1.13–4.96). This association was not significant in the crude model for NEAP. After adjustments for potential confounders, no significant relation was found between PRAL or NEAP and MUO phenotype (fully-adjusted model: OR 1.13; 95% CI 0.46–2.79; OR 1.21; 95% CI 0.48–3.07, respectively).

Table 3 Multivariate adjusted odds ratio (OR) and 95% confidence interval (CI) for MUO phenotype across tertiles of PRAL and NEAP 1.

Table 4 describes crude and multivariate adjusted odds ratio and 95% CI for MUO phenotype across tertiles of PRAL and NEAP, stratified by BMI categories. PRAL categories were not significantly associated with MUO status in overweight or in obese adolescents, based on both metabolic health status definitions. Among overweight adolescents, the highest category of NEAP, compared to the reference category, had nearly four times higher odds of MUO status, based on IDF criteria, in the crude model and model 1 (crude model: OR 4.04; 95% CI 1.35–12.06; model 1: OR 4.82; 95% CI 1.47–15.76). However, after adjustment for physical activity and SES, this association disappeared (fully-adjusted model: OR 2.03; 95% CI 0.47–8.75). When we examined these associations based on the IDF/HOMA-IR criteria, the same findings were found (crude model: OR 3.96; 95% CI 1.13–13.86; model 1: OR 4.96; 95% CI 1.23–19.89; and fully-adjusted model: OR 1.44; 95% CI 0.26–7.71).

Table 4 Multivariate adjusted odds ratio (OR) and 95% confidence interval (CI) for MUO phenotype across tertiles of PRAL and NEAP, stratified by BMI categories1.

Discussion

We found significant direct associations between PRAL and NEAP with odds of MUO status in Iranian adolescents, based on IDF criteria. However, these associations were dependant to covariates. The same finding was observed between PRAL and odds of MUO based on IDF/HOMA-IR definition, before adjustments for potential confounders. Stratified analysis indicated that greater NEAP categories were associated with higher odds for MUO among overweight adolescents, based on both definitions of metabolic health status. The present study is one of the first investigations on the association between DAL and metabolic health status in adolescents. The attainment of MHO may be a worthwhile first goal in the management of obesity and its complications, although MHO is an intermediate state before the development of MUO. Diet is one of the most important factors in determining obesity phenotype that should be considered to improve general health status13.

Our study showed possible beneficial relations between lower acid load of diet and metabolic health. Some prior studies have also investigated the association between DAL and obesity or metabolic syndrome. It is also worth noting that although some components of MHO/MUO and metabolic syndrome were the same, these are two different outcomes. Arisawa et al. have reported significant positive trends between NEAP scores with metabolic syndrome and some components of MetS in adults35. A systematic review and meta-analysis on observational studies additionally revealed that high DAL content was associated with higher obesity prevalence and serum TG concentrations17. A recent study demonstrated that the odds of obesity and abdominal adiposity significantly increased in tertile of DAL and also DAL could be an indicator of diet quality36. However, a previous cross-sectional study by Rahbarinejad et al. indicated no significant association between PRAL and NEAP and MetS in children before or after adjustment for potential confounders24. Also, another study found the same finding in adults, but higher scores of NEAP were associated with increased risk of impaired fasting blood glucose37. A national multi-center cross-sectional survey on 5326 Iranian students aged 6–18 years after evaluating anthropometric indices, showed a direct association between diet-induced acid load with neck circumference23. Furthermore, an inverse relationship was found between DAL indices and parental BMI23.

We observed that greater scores of PRAL and NEAP were related to higher FBG. Also, individuals with higher score of PRAL had higher TG. Some other investigations have also evaluated the relationship between DAL and cardiometabolic risk factors. A cross-sectional study on 371 Iranian adult women pointed out that women with higher NEAP scores had higher weight, WC, and TG concentrations. In case of PRAL, only a significant relation with TG was found20, which was similar to our finding. Lee et al. reported that diets with high acid-forming potential (high PRAL and NEAP scores) were positively associated with higher risk of IR during 7.3 years follow-up in a large population-based cohort study38. In another study on Japanese workers, high DAL was associated with IR, but not with FBG or HbA1c levels21. Another study by Haghighatdoost et al. demonstrated that PRAL and protein to potassium ratio were directly related to HbA1c among diabetic nephropathy patients39. Dehghan et al. revealed a potent negative relation between high DAL, particularly with PRAL scores, and cardiometabolic risk factors (including SBP, DBP, insulin, and odds of diabetes) in a systematic review and meta-analysis40. Differences in the findings of previous studies might be associated to the study populations, dietary intake assessment tools, sample size of the studies, applied methods to calculate dietary acid load, definition of dietary patterns, and various controlled confounders.

Several possible mechanisms might account for the relationship between DAL and MUO. Our study showed that adolescents in the top tertiles of PRAL and NEAP had greater intake of protein and carbohydrate, respectively. Animal food sources such as meat, fish, egg, chicken, cheese, and also cereals are rich in sulfur-containing amino acids, phosphorous and chloride. Consumption of these foods can contribute to higher acid production in the body; while high content of malate, citrate and glutamate in fruits and vegetables can contribute to base production. Therefore, animal-based foods, which are high in western diets, can be considered the most acid-producing diets and are associated with metabolic diseases41. These foods are pro-inflammatory and can increase inflammatory cytokines42, which are involved in lipid disorders and IR, particularly in obese subjects43. Higher DAL stimulates the production of glucocorticoids44. High cortisol levels induce lipase activity (lipoprotein lipase and hormone-sensitive lipase), which can lead to the increased efflux of fatty acids to the bloodstream and augment the production of very-low-density lipoproteins in the liver45. Furthermore, increased acidity of the body is related to the renal excretion of magnesium, calcium, and potassium46, although we did not measure serum concentration of these nutrients in our study. Increased BP and IR and obesity may result from these processes47,48. Also, diet-induced acidosis causes a decrease in adipokines (leptin and adiponectin), which can impair appetite suppression49,50. In addition, it has been demonstrated that high DAL reduces lean body mass and ultimately leads to increased body fat synthesis51.

The present study has some strength. First, the novelty of this study was the evaluation of DAL-metabolic health relation in a sample of adolescents from all socioeconomic districts of a large central city in Iran. Second, two different definitions of metabolic health status were used to determine MUO/MHO status. Third, several potential confounders were adjusted in the analysis. However, we acknowledge some limitations. The cross-sectional design of the study did not allow us to examine causality; so, conducting prospective studies is necessary. In addition, recall bias and other potential reporting biases may have affected the results. Also, the interview setting was used to collect dietary intakes data, which may be subject to social desirability biases. Despite controlling many potential confounders, residual confounders (such as the pattern of dietary habits, puberty, and sleep health) may affect the association between DAL scores and MUO. Finally, although BMI and WC were measured to define obesity and abdominal obesity, we could not measure body composition and fat distribution that could be involved in metabolic health status.

In conclusion, our findings indicated that diet-induced acid load might be associated with higher odds of MUO phenotype in overweight/obese adolescents. Thus, lower acidogenic food ingredients in the diets and high consumption of fruits and vegetables may improve the metabolic health status of adolescents. Further studies, particularly prospective nature, are needed to confirm these findings.