Abstracts

Introduction

Though several studies reported that dairy and dietary calcium could decrease the risk of obesity in children and adolescents ¹1. Abreu S, Santos R, Moreira CS, Vale S, Santos PC, Soares-Miranda L, et al. Association between dairy product intake and abdominal obesity in Azorean adolescents. Eur J Clin Nutr 2012; 66:830-5.^,22. Lee H-J, Cho J-i, Lee H-SH, Kim C-i, Cho E. Intakes of dairy products and calcium and obesity in Korean adults: Korean National Health and Nutrition Examination Surveys (KNHANES) 2007-2009. PLoS One 2014; 9:e99085.^,33 Jürimäe J, Mäestu E, Mengel E, Remmel L, Purge P, Tillmann V. Association between dietary calcium intake and adiposity in male adolescents. Nutrients 2019; 19:1454.^,44. Santos LC, Martini LA, Cintra IP, Fisberg M. Relationship between calcium intake and body mass index in adolescents. Arch Latinoam Nutr 2005; 88:345-9.^,55. Goldberg TB, Silva CC, Peres LN, Berbel MN, Heigasi MB, Ribeiro JM, et al. Calcium intake and its relationship with risk of overweight and obesity in adolescents. Arch Latinoam Nutr 2009; 59:14-21., other prospective studies show that this association is still inconclusive ⁶6. Berkey CS, Rockett HR, Willett WC, Colditz GA. Milk, dairy fat, dietary calcium, and weight gain. Arch Pediatr Adolesc Med 2005; 159:543-50.^,77. Huh SY, Rifas-Shiman SL, Rich-Edwards JW, Taveras EM, Gillman MW. Prospective association between milk intake and adiposity in preschool-aged children. J Am Diet Assoc 2010; 110:563-70.. Systematic reviews which suggested that dietary calcium or dairy intake is protective against excess weight in children, adolescents, and adults also had inconsistent results ⁸8. Louie JCY, Flood VM, Hector DJ, Rangan AM, Gill TP. Dairy consumption and overweight and obesity: a systematic review of prospective cohort studies. Obes Rev 2011; 12:582-92.^,99. Spence LA, Cifelli CJ, Millar GD. The role of dairy products in healthy weight and body composition in children and adolescents. Curr Nutr Food Sci 2011; 7:40-9..

Dietary calcium was first associated with obesity by McCarron ¹⁰10. McCarron DA. Dietary calcium as an antihypertensive agent. Nutr Rev 1984; 42:223-5. and later by Zemel. ¹¹11. Zemel MB. The role of calcium and dairy products in energy partitioning and weight management. Am J Clin Nutr 2004; 79 Suppl:907S-12S.. According to Zemel’s hypothesis, low calcium bioavailability increases the parathormone (PTH) and 1,25-Dihydroxyvitamin D3 synthesis as calcium-regulating mechanisms, leading to intracellular calcium concentration, which stimulates lipogenesis and suppresses lipolysis in the adipose tissue ¹¹11. Zemel MB. The role of calcium and dairy products in energy partitioning and weight management. Am J Clin Nutr 2004; 79 Suppl:907S-12S.^,1212. Major GC, Chaput JP, Ledoux M, St-Pierre S, Anderson GH, Zemel MB, et al. Recent developments in calcium-related obesity research. Obes Rev 2008; 9:428-45.. A meta-analysis of clinical control trials ¹³13. Gonzalez JT, Rumbold PL, Stevenson EJ. Effect of calcium intake on fat oxidation in adults: a meta-analysis of randomized, controlled trials. Obes Rev 2012; 13:848-57. showed that calcium intake around 1,000-1,400mg/day seems to enhance fat oxidation and is inversely correlated with PTH. Calcium intake of around 600-700mg/day is believed to boost the process of adipogenesis ¹⁴14. Jacqmain M, Doucet E, Després JP, Bouchard C, Tremblay A. Calcium intake, body composition, and lipoprotein-lipid concentrations in adults. Am J Clin Nutr 2003; 77:1448-52..

Brazilian nationwide surveys have shown that overweight and obesity increased in the country, reaching 20% of adolescents ¹⁵15. Instituto Brasileiro de Geografia e Estatística. Pesquisa de Orçamentos Familiares no Brasil, 2008/2009: antropometria e estado nutricional de crianças, adolescentes e adultos no Brasil. Rio de Janeiro: Instituto Brasileiro de Geografia e Estatística; 2010.. Meanwhile, the first Brazilian National Dietary Survey (INA) conducted in Brazil (2008-2009) ¹⁶16. Instituto Brasileiro de Geografia e Estatística. Pesquisa de Orçamentos Familiares, 2008-2009: análise do consumo alimentar pessoal no Brasil. Rio de Janeiro: Instituto Brasileiro de Geografia e Estatística; 2011. showed that, among Brazilian adolescents, the per capita dairy intake was lower than 100g/day and the average calcium intake was around 500mg/day ¹⁷17. Veiga GV, Costa RS, Araújo MC, Souza AM, Bezerra IN, Barbosa FS, et al. Inadequação do consumo de nutrientes entre adolescentes brasileiros. Rev Saúde Pública 2003; 47 Suppl 1:212S-21S.. Souza et al. ¹⁸18. Souza AM, Barufaldi LA, Abreu GA, Giannini DT, Oliveira CL, Santos MM, et al. ERICA: intake of macro and micronutrients of Brazilian adolescents. Rev Saúde Pública 2016; 50 Suppl 1:5s. found similar results in their analysis of a 2013 nationally representative school-based study with adolescents. Data from the 2017-2018 INA also showed that low calcium intake persists among Brazilian adolescents, estimating that 98.1% of them had inadequate intake of this micronutrient ¹⁹19. Verly-Jr E, Marchioni DM, Araujo MC, De Carli E, Oliveira DCRS, Yokoo EM, et al. Evolução da ingestão de energia e nutrientes no Brasil entre 2008-2009 e 2017-2018. Rev Saúde Pública 2021; 55 Supl 1:5s.. This low intake and the increased overweight and obesity provides an opportunity to examine the relationship between calcium intake and adiposity. We thus sought to assess the association between dietary calcium intake and adiposity indicators among a cohort of Brazilian adolescents.

Material and methods

Study population

The Adolescent Nutritional Assessment Longitudinal Study (ELANA) was carried out to assess changes in anthropometric and body composition measures in a cohort of adolescents in high school followed from 2010 to 2012 in the Metropolitan Region of Rio de Janeiro, Brazil ²⁰20. Moreira NF, Sichieri R, Reichenheim ME, Oliveira AS, Veiga GV. The associations of BMI trajectory and excessive weight gain with demographic and socio-economic factors: the Adolescent Nutritional Assessment Longitudinal Study cohort. Br J Nutr 2015; 114:2032-38..

The ELANA sample size was originally estimated at 1,200 adolescents, which would allow estimating the variation of one unit of body mass index (BMI) based on a 5% alpha error, 80% test power, mean BMI of 21.9kg/m², standard deviation (SD) of 3.0kg/m², and up to 20% lost to follow-up. The sample design was stratified by gender and type of school (public or private) with 300 individuals in each stratum. The repetition of the outcomes over time was considered in the sample size estimation. The sample size in each quintile of calcium intake could detect differences from 0.63 to 1.5 units between the outcomes, with statistical power ranging from 71 to 99% and a 5% alpha error.

All high school students (n = 1,131) of the two public and four private schools selected by convenience for the study were invited to participate in ELANA. Students with physical disabilities which prevented nutritional assessment, students on a diet, and pregnant female students were not eligible for the study. Moreover, 92 students (8.3%) could not be examined for lack of parental consent or refusal to participate in the study or for being absent from school. The ELANA cohort thus included 1,039 high school adolescents at baseline (2010, T₀). Our study included 962 students (aged 13-19 years) after excluding 77 (7%) individuals with incomplete information on food consumption or anthropometric/body composition measurements. In total, 779 adolescents were followed until 2011 (T₁) and 579 were followed until 2012 (T₂). Figure 1 shows the flowchart of the participants and the measurements taken at baseline and at follow-up. Further details on ELANA sample design and participant selection are available elsewhere ²⁰20. Moreira NF, Sichieri R, Reichenheim ME, Oliveira AS, Veiga GV. The associations of BMI trajectory and excessive weight gain with demographic and socio-economic factors: the Adolescent Nutritional Assessment Longitudinal Study cohort. Br J Nutr 2015; 114:2032-38..

The ELANA was approved by the Ethics Research Committee of the Institute of Social Medicine of the State University of Rio de Janeiro (IMS/UERJ, n. 0020.0.259.000-09). A written informed consent for the survey was provided by the adolescents and by their parents or guardians.

Anthropometric and body composition measurements

Anthropometric measures included body weight (kg), height (cm), and waist circumference (cm), assessed by well-trained staff and in accordance with standard protocols ²¹21. Lohman TG, Martorell R. Anthropometric standardization reference manual. Champaign: Human Kinetics Books; 1988. in 2010, 2011, and 2012.

Body weight was measured using a portable digital scale (Kratos, Cotia, Brazil) with 0.1kg of variation while the participant wore light clothes and no shoes. A stadiometer (Alturaexata, Belo Horizonte, Brazil) with a variation of 0.1cm was used to measure height twice, accepting a maximum variation of 0.5cm between the two measures. The average of both measurements was used in the analysis. BMI (kg/m²) was estimated. Excess weight was defined as BMI Z-score > +1 according to gender and age distributions ²²22. World Health Organization. Growth reference data for 5-19 years: body mass index-for-age, length/height-for-age and weight-for-height. Geneva: World Health Organization; 2007.. Waist circumference (WC) was assessed at the narrowest circumference of the trunk with an inelastic tape (Seca, Cotia, Brazil) with 0.1 cm accuracy. The WC was measured twice, allowing a maximum variation of 1.0cm between the two measurements, whose average was used in the analysis. Considering the absence of an international reference for WC distribution for adolescents, abdominal obesity was defined as values above the 90th percentile of WC distribution, specified by gender.

Body composition was measured at baseline and at the first year of follow-up (2011) using a bioelectrical impedance device (BIA RJL System, 101Q, Clinton Township, United States) with tetrapolar electrodes. Subjects were in the supine position and the electrodes were placed on the dorsal surface of their right foot and hand. A 50-kHz frequency electrical current was employed to measure resistance and reactance. Resistance values were used to estimate fat-free mass (FFM) using the equation for adolescents validated by Houtkooper et al. ²³23. Houtkooper LB, Going SB, Lohman TG, Roche AF, Van Loan M. Bioelectrical impedance estimation of fat-free body mass in children and youth: a cross-validation study. J Appl Physiol 1992; 72:366-73.. Fat mass (FM) was calculated as the difference between weight and fat-free mass. Body fat percentage (%BF) was obtained from: [(fat mass/total weight) * 100]. Excess body adiposity was defined as %BF ≥ 25% for boys and ≥ 30% for girls ²⁴24. Williams DP, Going SB, Lohman TG, Harsha DW, Srinivasan SR, Webber LS, et al. Body fatness and risk for elevated blood pressure, total cholesterol, and serum lipoprotein ratios in children and adolescents. Am J Public Health 1992; 82:358-63.. Since BMI is a global adiposity index, the body fat fraction was estimated using the fat mass index [FMI = fat mass/height² (kg/m²)] whereas the body fat-free fraction was estimated with the fat-free mass index [FFMI = fat-free mass/height² (kg/m²)].

Dietary intake assessment

Food intake was assessed using a qualitative self-administered Food Frequency Questionnaire (FFQ) referring to consumption in the last three months before the interview. The FFQ included 72 items and eight options for consumption frequency, ranging from “less than once a month or never” to “four or more times a day”. This is a reduced version of an FFQ validated for adolescents from Rio de Janeiro ²⁵25. Araujo MC, Yokoo EM, Pereira RA. Validation and calibration of a semiquantitative Food Frequency Questionnaire designed for adolescents. J Am Diet Assoc 2010; 110:1170-7.. To reduce the original FFQ, items with 90% of energy intake were selected and items consumed by less than 5% of the original group evaluated were excluded. Furthermore, several items were grouped into a broad category; different types of bread were simply listed as “bread”. The three food records which assessed the validation of the original FFQ were used to assess the relative validity of the short FFQ, showing similar results: deattenuated and adjusted Pearson’s correlation coefficients ranged from 0.16 to 0.42 for the reduced questionnaire (unpublished data) and from 0.17 to 0.47 for the complete version ²⁵25. Araujo MC, Yokoo EM, Pereira RA. Validation and calibration of a semiquantitative Food Frequency Questionnaire designed for adolescents. J Am Diet Assoc 2010; 110:1170-7..

For the analysis, calibrated FFQ energy and nutrient intake were estimated as described by Araujo et al. ²⁵25. Araujo MC, Yokoo EM, Pereira RA. Validation and calibration of a semiquantitative Food Frequency Questionnaire designed for adolescents. J Am Diet Assoc 2010; 110:1170-7.. In short, the database used for validating the original FFQ was used for the calibration process and, for each nutrient, linear regression models adjusted for age and gender were run, with the nutrient intake estimated from the mean of three food records as the dependent variable and from the FFQ as the independent variable. The intercept (α) and the slope (β) regression coefficients were computed to calibrate the energy and nutrient intakes estimated by the FFQ using the following equation:

The reported food frequency consumption was transformed into daily frequency and multiplied by the standard serving size associated with each item. The daily intake of energy, protein, carbohydrates, vitamins, and minerals - including calcium - was estimated using the Brazilian Food Composition Table (TBCA) ²⁶26. Instituto Brasileiro de Geografia e Estatística. Pesquisa de Orçamentos Familiares. Tabela de Composição Nutricional dos Alimentos Consumidos no Brasil. Rio de Janeiro: Ministério do Planejamento, Orçamento e Gestão; 2011..

Finally, dietary calcium intake was categorized into quintiles and the following cutoffs were estimated: 1^st quintile (Q1): < 492.0mg/day; 2^nd quintile (Q2): 492.0-567.9mg/day; 3^rd quintile (Q3): 568.0-670.9mg/day; 4^th quintile (Q4): 671.0-838.0mg/day; and 5^th quintile (Q5): > 838.0mg/day.

Other covariates

Age (as a continuous variable), gender, and type of school (public and private) were obtained by a structured self-administered questionnaire, applied in class under supervision of the nutritionists responsible for ELANA. Physical activity was measured using the short version of the International Physical Activity Questionnaire (IPAQ) validated for Brazilian adolescents ²⁷27. Guedes DP, Lopes CC, Guedes JERP. Reprodutibilidade e validade do questionário internacional de atividade física em adolescentes. Rev Bras Med Esporte 2005; 11:151-8.. Subjects were asked about their frequency and duration (in minutes) of walking, practice of specific moderate and vigorous activities, and sedentary habits (e.g., sitting) in the week before the interview. Physical activity was then categorized as: “very active” (vigorous: ≥ 5 days/week and ≥ 30 minutes per time and/or ≥ 3 days/week and ≥ 20 minutes per time + moderate and/or walking: ≥ 5 days/week and ≥ 30 minutes per time); “active” (vigorous: ≥ 3 days/week and ≥ 20 minutes per time and/or moderate or walking: ≥ 5 days/week and ≥ 30 minutes per time); “irregularly active A” (vigorous + moderate + walking: 5 days/week or 150 minutes/week); “irregularly active B” (those who did not match the frequency or duration of physical activity recommendations); or “sedentary” (those who do not practice physical activity for at least 10 minutes continuously per week) ²⁸28. Matsudo SM, Matsudo VR, Araújo T, Andrade D, Andrade E, Oliveira L, et al. Nível de atividade física da população do Estado de São Paulo: análise de acordo com o gênero, idade, nível socioeconômico, distribuição geográfica e de conhecimento. Rev Bras Ciênc Mov 2002; 10:41-50.. In this study, physical activity categories were divided into dichotomous categories: “active/highly active” (active + very active ) and “sedentary/less active” (“sedentary” + “irregularly active A” + “irregularly active B”).

Daily energy (Kcal/day) and vitamin D (mcg/day) intake and daily consumption frequency of fruits, vegetables, sugar-sweetened beverages, fast foods, and breakfast cereals were considered as potential confounders since they are usually associated with healthy or unhealthy dietary patterns and weight gain. All these variables were included in the longitudinal analysis to adjust the regression models.

Statistical analyses

Means and SD of continuous variables and proportions of categorical variables were estimated across quintiles of dietary calcium intake. For cross-sectional analysis, means of BMI, WC, FM, FFM, %BF, FMI, and FFMI were compared by analysis of variance (ANOVA) according to dietary calcium quintiles; Tukey’s Post Hoc test compared the 1st quintile with the 2nd to the 5th quintiles. The Tukey’s test compares the difference between each pair of means with appropriate adjustment for the multiple testing. The ANOVA F-test was applied to test for linear trend (increasing or decreasing) of means on each dependent variable across calcium intake quintiles.

For longitudinal analysis, linear mixed-effects models (LMM) were applied to assess the effect of dietary calcium intake on each adiposity measure (dependent variables) from baseline to the end of follow-up. Models were run for the whole sample and stratified by gender. The age (in years) of the adolescents at each point of follow-up was considered as the random effects in the LMM to guarantee that changes could be attributable to time instead of age. For each outcome variable, we constructed separated models including the random effects of both the intercepts and slopes. The correlation matrix was unstructured since this type generates full models, which are the best estimates of the regression parameters considering unbalanced and incomplete data (missing outcome variables), prevalent in longitudinal studies ²⁹29. Fitzmaurice GM, Laird NM, Ware JH. Longitudinal and clustered data. In: Fitzmaurice GM , Laird NM , Ware JH , editors. Applied longitudinal analysis. 2nd Ed. Hoboken: Wiley; 2011. p. 1-18.. Changes in time of the outcome variables across quintiles of dietary calcium intake were assessed by including the interaction term between age and dietary calcium (age * dietary calcium intake) on LMM. A significant interaction of age * dietary calcium intake indicates a differential growth rate for each adiposity measure over time.

The multicollinearity of confounder variables in regression models was assessed by the variance inflation factor (VIF). The VIF values for each confounder variable ranged from 1.13 to 3.14 and were lower (except for energy with VIF = 3.14) than the cutoff point of 2.5, which indicates considerable collinearity ³⁰30. Johnston R, Jones K, Manley D. Confounding and collinearity in regression analysis: a cautionary tale and an alternative procedure, illustrated by studies of British voting behaviour. Qual Quant 2018; 52:1957-76..

The Cochran-Mantel-Haenszel test was applied to assess differences between the proportions of cases at baseline and new cases at the end of follow-up regarding excess weight, abdominal obesity, and excess body fat in the calcium intake quintiles.

Statistical analyses were conducted using SPSS software (https://www.ibm.com/) and statistical significance was set at p-value > 0.05.

Results

Table 1 shows the characteristics of the study population at baseline according to quintiles of dietary calcium. At baseline, 53% of participants were women, 50% were from private schools, and 71% were active/highly active. The Q5 of calcium intake included more boys (32%) than girls (9%) and, regarding type of school, 24% of students in Q5 were from public schools whereas 16% were from private schools. Adolescents in Q5 were also classified according to categories of physical activity, in which most were active/highly active (23%) rather than sedentary/low active (15%). The mean age was 15.7 (SD = 0.9) years. Overall, the mean calcium intake was 700 (SD = 310) mg/day. Means of calcium, vitamin D, energy intake, consumption frequency of breakfast cereals, fast foods, fruits/vegetables, and sugar-sweetened beverages were higher in Q5 than in Q1 (Table 1).

Table 2 shows significant differences between most of the baseline anthropometric and body composition variables means at highest quintiles of calcium intake and at the first quintile. For the total sample, means of BMI (p for trend = 0.003), FM (p for trend < 0.01), %BF (p for trend < 0.01), and FMI (p for trend < 0.01) were lower at Q4 and Q5 than in Q1. Differences between means of FFM (p for trend < 0.01) and FFMI (p for trend = 0.001) from Q3 to Q5 were not significant compared to differences with Q1. Boys had lower means of BMI (p for trend = 0.003), WC (p for trend = 0.001), FM (p for trend = 0.01), and FMI (p for trend = 0.01) in Q5 than in Q1; they also had lower %BF means (p for trend < 0.01) in both Q4 and Q5. Among girls, means of adiposity indicators from Q2 to Q5 were equal compared to Q1, except for %BF means, which decreased across the quintiles of calcium intake (p for trend = 0.01) (Table 2).

Table 3 shows the adjusted anthropometric and body composition variable means (and SD) at the end of follow-up according to the quintiles of calcium intake at baseline, estimated from linear mixed models for longitudinal analysis. Changes in the means over time are estimated with the interaction term between age * dietary calcium intake. For the total sample and for boys, those at Q4 (p for interaction = 0.03) and Q5 (p for interaction = 0.02) quintiles had lower FMI means at the end of follow-up than at baseline. Girls at Q4 of calcium intake had the lowest WC (p for the interaction = 0.03) compared to those at Q1.

At baseline, the proportions of excess weight were lower at Q5 than at Q1 for the total sample (p = 0.02), for boys (p = 0.02), and for girls (p = 0.06). Accordingly, the proportion of abdominal obesity was lower at Q5 than Q1 for the total sample (p = 0.06) and for boys (p = 0.07). The proportion of excess body fat was lower at Q5 for the total sample (p < 0.01), boys (p = 0.002), and girls (p = 0.003). No association was observed between calcium intake and the incidence of excess weight, abdominal adiposity, and excess body fat at the end of the follow-up (Table 4).

Subjects lost to follow-up (n = 437) differ from those who completed the study regarding gender and type of school; however, since 76% (n = 779) of the cohort completed the first follow-up, contributing with at least two observations during the period, they were included in the longitudinal analyses with the LMM. Moreover, losses were more frequent in public schools (59.5%), likely because these students often have higher mobility between schools than students from private ones. No significant differences were observed between subjects that completed the study and those lost to follow-up regarding BMI (p = 0.31), WC (p = 0.07), and %BF (p = 0.19) at baseline (data not shown).

Discussion

In this study, at baseline, adolescents at Q4 and Q5 quintiles had lower means of body adiposity indicators (BMI, FM, %BF, and FMI) and higher means of lean mass indicators (FFM and FFMI) than those at the lowest quintile, with differences according to gender. Longitudinally, calcium intake was inversely associated with changes over time for FMI in boys at Q4 and Q5 and for WC in girls at Q4. Q5 had a lower proportion of excess weight and excess body fat than Q1 at baseline, but not at the end of the follow-up.

These results corroborate similar international ¹1. Abreu S, Santos R, Moreira CS, Vale S, Santos PC, Soares-Miranda L, et al. Association between dairy product intake and abdominal obesity in Azorean adolescents. Eur J Clin Nutr 2012; 66:830-5.^,22. Lee H-J, Cho J-i, Lee H-SH, Kim C-i, Cho E. Intakes of dairy products and calcium and obesity in Korean adults: Korean National Health and Nutrition Examination Surveys (KNHANES) 2007-2009. PLoS One 2014; 9:e99085.^,33 Jürimäe J, Mäestu E, Mengel E, Remmel L, Purge P, Tillmann V. Association between dietary calcium intake and adiposity in male adolescents. Nutrients 2019; 19:1454. and national studies ⁴4. Santos LC, Martini LA, Cintra IP, Fisberg M. Relationship between calcium intake and body mass index in adolescents. Arch Latinoam Nutr 2005; 88:345-9.^,55. Goldberg TB, Silva CC, Peres LN, Berbel MN, Heigasi MB, Ribeiro JM, et al. Calcium intake and its relationship with risk of overweight and obesity in adolescents. Arch Latinoam Nutr 2009; 59:14-21. that reported an inverse association between calcium intake and adiposity in adolescents. A systematic review of longitudinal studies (n = 10) ⁸8. Louie JCY, Flood VM, Hector DJ, Rangan AM, Gill TP. Dairy consumption and overweight and obesity: a systematic review of prospective cohort studies. Obes Rev 2011; 12:582-92. and a review conducted by Spence et al. ⁹9. Spence LA, Cifelli CJ, Millar GD. The role of dairy products in healthy weight and body composition in children and adolescents. Curr Nutr Food Sci 2011; 7:40-9. including clinical trials (n = 11) and cross-sectional (n = 23) and longitudinal studies (n = 13) also found that calcium and/or dairy helps prevent obesity. A meta-analysis including 33 studies, out of which nine were conducted with children and adolescents, found an inverse correlation between calcium intake and body weight gain among these age groups ³¹31. Li P, Fan C, Lu Y, Qi K. Effects of calcium supplementation on body weight: a meta-analysis. Am J Clin Nutr 2016; 104:1263-73., indicating that calcium intake could benefit obesity prevention more than obesity treatment.

The association between calcium intake and weight gain is likely inconsistent because of the variety of studies, which have large cultural differences in calcium/dairy consumption, different lengths of follow-up, limited outcome variables (usually based on BMI), and different methods for diet assessment and statistical analysis procedures.

Furthermore, a threshold of dietary calcium intake around 600-700mg/day ¹⁴14. Jacqmain M, Doucet E, Després JP, Bouchard C, Tremblay A. Calcium intake, body composition, and lipoprotein-lipid concentrations in adults. Am J Clin Nutr 2003; 77:1448-52. may start the biological mechanisms that regulate calcium levels, thus increasing intracellular calcium concentrations in the adipogenesis process into adiposity cells ¹¹11. Zemel MB. The role of calcium and dairy products in energy partitioning and weight management. Am J Clin Nutr 2004; 79 Suppl:907S-12S.^,1212. Major GC, Chaput JP, Ledoux M, St-Pierre S, Anderson GH, Zemel MB, et al. Recent developments in calcium-related obesity research. Obes Rev 2008; 9:428-45.. Individuals whose usual calcium intake exceeds this threshold could be on calcium homeostasis with no increases in adipogenesis rate in fat tissues. This condition could partly explain the controversial results in epidemiological studies since the definition of low or high calcium intake varies according to its original population and may be above or below this biological threshold.

The average calcium intake in Brazilian adolescents is low. The 2017-2018 INA estimated a 457mg intake ¹⁹19. Verly-Jr E, Marchioni DM, Araujo MC, De Carli E, Oliveira DCRS, Yokoo EM, et al. Evolução da ingestão de energia e nutrientes no Brasil entre 2008-2009 e 2017-2018. Rev Saúde Pública 2021; 55 Supl 1:5s., similarly to the threshold of Q1 estimated in our study (450mg/day) but below the threshold associated with the onset of adipogenesis into adiposity cells ¹⁴14. Jacqmain M, Doucet E, Després JP, Bouchard C, Tremblay A. Calcium intake, body composition, and lipoprotein-lipid concentrations in adults. Am J Clin Nutr 2003; 77:1448-52.. On the other hand, the mean calcium intake in Q5 was 1,174mg/day, close to the estimated average requirement (EAR) of calcium (1,100mg) for 9-18 year old children and adolescents ³²32. Institute of Medicine. Dietary reference intakes for calcium and vitamin D. Washington DC: The National Academies Press; 2011.. Even if half of the adolescents in the highest quintile do not reach the recommended intake, the Q5 cutoff limit (838.0mg/day) is higher than the threshold associated with the increased risk of overweight ¹⁴14. Jacqmain M, Doucet E, Després JP, Bouchard C, Tremblay A. Calcium intake, body composition, and lipoprotein-lipid concentrations in adults. Am J Clin Nutr 2003; 77:1448-52..

The findings of this cross-sectional analysis support Zemel’s hypothesis ¹¹11. Zemel MB. The role of calcium and dairy products in energy partitioning and weight management. Am J Clin Nutr 2004; 79 Suppl:907S-12S. considering that adolescents at Q5, especially boys, had lower means for most adiposity indicators and higher means of lean mass indicators. The longitudinal analysis found that high calcium intake was inversely associated with changes over time (calcium intake * age) for FMI in boys at Q4 and Q5 and for WC in girls at Q4. At baseline, boys were more likely to be classified at Q5 (32%) than girls (9%), which could explain the gender differences related to the association between calcium intake and adiposity indicators. Moreover, waist circumference alone is also a good indicator for visceral adipose depots ³³33. Katzmarzyk PT, Srinivasan SR, Chen W, Malina RM, Bouchard C, Berenson GS. Body mass index, waist circumference, and clustering of cardiovascular disease risk factors in a biracial sample of children and adolescents. Pediatrics 2004; 114:198-205.^,3434. Rodríguez G, Moreno LA, Blay MG, Blay VA, Garagorri JM, Sarría A, et al. Body composition in adolescents: measurements and metabolic aspects. Int J Obes Relat Metab Disord 2004; 28 Suppl 3:S54-8. in children and adolescents. Waist circumference values increased over time could be more sensitive to detect abdominal fat changes, especially in girls, as shown in this study ³⁵35. Garnett SP, Baur LA, Cowell CT. The prevalence of increased central adiposity in Australian schoolchildren 1985 to 2007. Obes Rev 2011; 12:887-96.

Our findings are similar to those of other longitudinal studies. The Avon Longitudinal Study of Parents and Children - ALSPAC (n = 2,455) found that children with high calcium intake (1,140mg/day) at 10 years old had a lower risk of excessive body weight and FM at 13 years old ³⁶36. Bigornia SJ, LaValley MP, Moore LL, Northstone K, Emmett P, Ness AR, et al. Dairy intakes at age 10 years do not adversely affect risk of excess adiposity at 13 years. J Nutr 2014; 144:1081-90. than children with low calcium intake (564mg/day). Our results are also similar to those of the IDEFICS study ³⁷37. Nappo A, Sparano S, Intemann T, Kourides YA, Lissner L, Molnar D, et al. Dietary calcium intake and adiposity in children and adolescents: cross-sectional and longitudinal results from IDEFICS/I.Family cohort. Nutr Metab Cardiovasc Dis 2019; 29:440-9., a large European cohort of children and adolescents (n = 6,696) which found an inverse association at baseline of total calcium intake with all the adiposity indices in boys and with the sum of skinfolds and FMI among girls. The authors of the cohort study also found that the prevalence of overweight/obesity decreased across the tertiles of calcium intake for both sexes. Prospectively, at baseline, boys at the highest tertile of calcium intake had significantly lower increase in BMI, WC, and FMI z-scores whereas girls had lower increase in WC z-score only. Moreover, boys had lower risk of overweight/obesity at the highest tertile of calcium intake than in the first ³⁷37. Nappo A, Sparano S, Intemann T, Kourides YA, Lissner L, Molnar D, et al. Dietary calcium intake and adiposity in children and adolescents: cross-sectional and longitudinal results from IDEFICS/I.Family cohort. Nutr Metab Cardiovasc Dis 2019; 29:440-9..

Similarly to our findings for boys, Barr ³⁸38. Barr SI. Calcium and body fat in peripubertal girls: cross-sectional and longitudinal observations. Obesity (Silver Spring) 2007; 15:1302-10. observed that girls with calcium intake below the median (773mg/day) had higher baseline body weight, total FM, percentage body fat, and percentage visceral fat than those above the median, with no changes after two years of follow-up. However, in a randomized controlled trial ³⁹39. Lappe JM, McMahon DJ, Laughlin A, Hanson C, Desmangles JC, Begley M, et al. The effect of increasing dairy calcium intake of adolescent girls on changes in body fat and weight. Am J Clin Nutr 2017; 105:1046-53., girls aged 13-14 years with usual calcium intake ≤ 600mg/day were supplemented with ≥ 1,200mg/day of calcium from low-fat milk and yogurt for 12 months, gaining amount of body fat similar to the control group (≤ 600mg/day). This shows no association between calcium or dairy and the decrease of body fat or weight gain in girls. These results lead us to ask whether the supplementation of calcium and dairy can counteract a physiological stimulus to weight and fat gain, as commonly observed during adolescence. To our understanding, they are not.

Dairy products are the major calcium source in a diet and rich in other nutritional components such as conjugated linoleic acid (CLA) and vitamin D, both involved in the lipolysis process ⁴⁰40. Dougkas A, Reynolds CK, Givens ID, Elwood PC, Minihane AM. Associations between dairy consumption and body weight: a review of the evidence and underlying mechanisms. Nutr Res Rev 2011; 24:72-95.. Moreover, other components such as lactose, cholecalciferol (vitamin D3), and casein phosphopeptides enhance their bioavailability ⁴¹41. Cashman KD. Calcium intake, calcium bioavailability and bone health. Br J Nutr 2002; 87 Suppl 2:S169-77.^,4242. Caroli A, Poli A, Ricotta D, Banfi G, Cocchi D. Invited review: dairy intake and bone health: a viewpoint from the state of the art. J Dairy Sci 2011; 94:5249-62.. This effect’s attribution to dietary calcium per se or to other components in dairy products, or even both, is still undefined. In fact, two meta-analyses ⁴³43. Dror DK. Dairy consumption and pre-school, school-age and adolescent obesity in developed countries: a systematic review and meta-analysis. Obes Rev 2014; 15:516-27.^,4444. Lu L, Xun P, Wan Y, Hi K, Cai W. Long-term association between dairy consumption and risk of childhood obesity: a systematic review and meta-analysis of prospective cohort studies. Eur J Clin Nutr 2016; 70: 414-23 suggested that dairy consumption is inversely and longitudinally associated with the risk of childhood and adolescence overweight/obesity.

This study has limitations. The period of follow-up may have been too short to capture small changes on anthropometric and body composition measures as those observed in the study. We did not consider the sexual maturation stage, but adolescents included in this study had a mean age of 15.7 years, indicating that most of them may have completed all five stages of pubertal development ⁴⁵45. Tanner JM. Growth and endocrinology of the adolescent. In: Gardner LI, editor. Endocrine and genetic diseases of childhood and adolescents. 2nd Ed. Philadelphia: WB Saunders; 1975. p. 14-64.; therefore, changes in body weight and composition may not result from the hormonal transformations of puberty. Moreover, linear mixed models included age as the random effect which seems to minimize the absence of sexual maturation stage information. The mixed linear models can generate better estimators of changes over time by minimizing the effect of common problems related to longitudinal designs, including the drop-outs during follow-up, unbalanced data, and unequal spacing between time intervals ²⁹29. Fitzmaurice GM, Laird NM, Ware JH. Longitudinal and clustered data. In: Fitzmaurice GM , Laird NM , Ware JH , editors. Applied longitudinal analysis. 2nd Ed. Hoboken: Wiley; 2011. p. 1-18.. Finally, the losses to follow-up did not bias our results because they did not differ regarding the main exposures and outcomes.

This study also has strengths. Firstly, we analyzed indicators of body fat besides BMI to assess changes in weight and fat gain during follow-up. Secondly, we applied a validated FFQ specially designed for adolescents including the representative main food sources of calcium in the Brazilian diet. Finally, we calibrated dietary data to improve the energy and nutrient estimates and minimize information bias usually associated with the FFQs.

Conclusions

This study’s results corroborate Zemel’s hypothesis about the effect of dietary calcium on adiposity in adolescents, especially among boys. To our knowledge, this is the first prospective study in Brazil to analyze the association between calcium intake and adiposity measures in adolescents. This study supports the hypothesis that dietary calcium intake is inversely associated with changes in adiposity among adolescents in a low/middle-income country with a high prevalence of calcium intake inadequacy. Further studies should consider that a threshold of low calcium intake could enhance adipogenesis.

Acknowledgments

We would like to thank the Brazilian National Research Council (CNPq; grant 47667/2011-9), the Rio de Janeiro State Research Foundation (FAPERJ; grants E26/110.847/2009, E26/110.626/2011 and E-26/110.774/2013), and the Brazilian Graduate Studies Coordinating Board (CAPES; grant 23038.007702/2011-5) for funding the project and the ELANA for graduate support fellowships.

References

Publication Dates

History