Introduction
Low bone mineral density (BMD), whether osteopenia or osteoporosis, has become a worldwide interest as the population is ageing globally(Reference Friedman and Mendelson1). A progressive decline in BMD begins as peak bone mass (PBM) is reached between 20 and 30 years of age, which occurs initially independent of sex hormone levels(Reference Hendrickx, Boudin and Van Hul2). Thus, both ageing and PBM attained during childhood and adolescence can be considered predictive factors of bone health at an older age(Reference Hendrickx, Boudin and Van Hul2,Reference Rizzoli, Bianchi and Garabedian3) .
Bone mass gain is influenced by various factors, such as sex, ethnicity, hormonal status and lifestyle-related factors, including diet, physical activity, cigarette smoking and alcohol intake(Reference Rizzoli, Bianchi and Garabedian3,Reference Curtis, Litwic and Cooper4) . Several studies of children and young adults have also proposed that consumption of carbonated soft drinks is associated with low BMD(Reference McGartland, Robson and Murray5,Reference Alghadir, Gabr and Al-Eisa6) and even increased risk of fractures(Reference Wyshak7,Reference Manias, McCabe and Bishop8) . In one experimental study, a significant delay in alveolar bone formation was observed in cola consuming rats(Reference Teofilo, Leonel and Lamano9). Furthermore, colas, but not other carbonated beverages, are related to low BMD in older women(Reference Tucker, Morita and Qiao10) or increased upper limb fractures in children(Reference Ma and Jones11). However, the results of the previous studies were different according to sex(Reference McGartland, Robson and Murray5), and the association between cola and BMD has been found mainly in women(Reference Wyshak7,Reference Tucker, Morita and Qiao10,Reference Wyshak and Frisch12) . Displacement of beneficial nutrient-rich foods such as milk(Reference McGartland, Robson and Murray5,Reference Kristensen, Jensen and Kudsk13) , decreased physical activity owing to watching television or computer while drinking colas(Reference Ma and Jones11), and the presence of phosphoric acid and caffeine in colas that interfere with calcium absorption and excretion(Reference Amato, Maravilla and Montoya14,Reference Rapuri, Gallagher and Kinyamu15) have been suggested as possible mechanisms. However, a clear mechanism has not yet been proposed, and a set of contributing factors may exist(Reference Fung, Arasaratnam and Grodstein16).
Most studies on the relationship between cola intake and BMD have been conducted with the Western populations, with very few studies involving Asian populations. Considering the clinical importance of PBM gain during childhood and adolescence and the increasing trend of consumption of carbonated soft drink in South Korea during this critical period for bone health(17,Reference Han, Kim and Powell18) , examining whether cola is related to low BMD in the growing young population is worthwhile. Thus, we aimed to examine the associations between cola consumption and BMD in Korean adolescents and young adults.
Methods
Data collection and study participants
The present study was based on data acquired from the Korea National Health and Nutrition Examination Survey (KNHANES) from 2008 to 2011. The KNHANES is a nationwide and population-based cross-sectional health survey conducted periodically by the Korean Ministry of Health and Welfare. This survey consists of three distinctive surveys: health interviews; health examinations; nutrition. Data were collected via household interviews and by direct standardised physical examinations in specially equipped mobile examination centres. The survey methods have been described elsewhere(Reference Lim, Kim and Rhee19).
BMD was initially measured in participants who were older than 19 years, but the target population was expanded to those who were older than 10 years since July 2009. A total of 20 808 participants were assessed for the body composition using dual-energy X-ray absorptiometry (DXA) from 2008 to 2011. Among these participants, we limited our study population to adolescents and young adults aged 12–25 years, the age when BMD has been reported to increase and reach a plateau(Reference Hendrickx, Boudin and Van Hul2,Reference Kang, Hong and Chung20) . We excluded subjects who were diagnosed with asthma, liver cirrhosis, renal failure, thyroid disease or any kind of cancer that can affect BMD. Subjects who were being treated for osteoporosis were also excluded. To adjust for the sex hormone effect, premenarchal female participants were also excluded. Finally, a total of 1222 males and 1277 females aged 12–25 years were enrolled as the study population for analysis.
BMD measurements
BMD (g/cm2) and body composition were measured using DXA (Hologic Discovery, Hologic Inc., Bedford, MA, USA). As the posterior–anterior spine and whole body except the head are recommended skeletal sites for BMD measurements in paediatric populations(Reference Crabtree, Arabi and Bachrach21), we used BMD at the whole body which did not include the head, whole femur, femoral neck and lumbar spine for analysis.
Dietary assessment
The nutritional survey was conducted via a face-to-face interview with dietitians. The usual dietary habits were assessed using a 24-h dietary recall questionnaire, and the frequency of sixty-three food-item consumption was also investigated. The frequency of food consumption was based on dietary habits during the past 1 year. We assessed the daily intake amount of energy and other nutrients based on the calculated data using food frequency questionnaires presented in the KNHANES. We also evaluated whether calcium, phosphorus and energy intake were sufficient according to the dietary reference intake for Koreans(22). Since the frequency of cola intake was not investigated in KNHANES, we estimated the level of cola consumption by summing the amount of cola intakes regardless of the brand, during a day using a 24-h dietary recall questionnaire. We also divided the study population into two groups based on their reports of cola consumption using a 24-h dietary recall data to examine the relationship between cola intake and BMD. Subjects who drank any kind of cola were categorised as cola drinkers. Subjects who did not reportedly drink any cola beverages were classified as non-cola drinkers. The amount of milk consumed daily was estimated using the same method. The frequency of food consumption was originally categorised into ten groups in the KNHANES: three times a day; twice a day; once a day; 4–6 times a week; 2–3 times a week; once a week; 2–3 times a month; once a month; 6–11 times a year; nearly non-existent. We regrouped these ten groups into five categories to examine the distribution of coffee and milk consumption. ‘Rarely’ (represented subjects who reported consumption nearly none or 6–11 times a year), ‘1–3 times/month’, ‘1–3 times/week’, ‘4–7 times/week’, ‘≥twice a day’.
Measurement of covariates
Height and weight were measured using a standard method to the nearest 0⋅1 cm for height and 0⋅1 kg for weight. Waist circumference was measured at the midpoint between the lowest level of the ribs and iliac crest at the mid-axillary line. We used the calculated values of body mass index (BMI) presented in the original data. Lifestyle-related components were assessed using a self-administered questionnaire. Physical activity was categorised as insufficient exercise, which included no exercise and low-intensity exercise, and sufficient exercise, which included moderate-to-vigorous-intensity exercise. Low-intensity exercise was defined as walking for more than 5 d/week for more than 30 min at a time in the last week. Moderate-to-vigorous-intensity exercise was defined as participating in moderate-intensity activities, such as swimming slowly or badminton, for more than 5 d/week for more than 30 min at a time or participating in vigorous-intensity activities, such as running, climbing or swimming fast for more than 3 d/week for more than 20 min at a time in the last week. Smoking status was divided into never smokers and ever smokers. Subjects who had smoked 100 cigarettes or more in their lifetime were categorised as ever smokers. Those who had never smoked or smoked less than 100 cigarettes in their lifetime were defined as never smokers. Alcohol consumption was assessed according to the frequency and amount of alcohol intake at a time. High-risk alcohol consumption was defined as drinking alcohol more than twice a week, with the average amount per drink of ≥7 glasses for males, or ≥5 glasses for females. We also assessed the serum 25(OH)D level.
Statistical analyses
All analyses were conducted using the SPSS software for Windows version 21.0 (SPSS Inc., Chicago, IL, USA). We used sampling weights based on the complex sample design of KNHANES in all analyses to obtain representatives of the Korean population. χ 2 statistics and one-way analysis of variance were used to assess the difference between the baseline characteristics based on cola intake and sex. As the association between cola and BMD has been reported differently based on sex in previous studies, we examined the association stratified by sex first. Then, we divided the study population into two groups: cola drinkers and non-cola drinkers. To examine the association between calcium intake and cola consumption, we calculated the tertiles of daily calcium intake in both sexes. Multiple linear regression was used to examine the association between cola intake and BMD. Since BMD values were not normally distributed in this study population, they were log-transformed to approximate a normal distribution for the analyses. To examine the independent association between cola intake and BMD, we adjusted the potential influence of age, BMI, physical activity, smoking, alcohol consumption, daily calcium intake, daily phosphorus intake, total daily energy intake and serum 25(OH)D level. We additionally adjusted oestrogen exposure in the female population. Since the sex hormone level was not investigated in KNHANES, we estimated the degree of oestrogen exposure by subtracting the age of menarche from the present age. A P-value of less than 0⋅05 was considered statistically significant in all analyses.
Ethical standards disclosure
The present study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving human subjects were approved by the Institutional Review Board at Sungkyunkwan University. All the participants in this survey signed an informed consent form, and legal guardian consent was obtained for participants aged under 14.
Results
Characteristics of study participants
Table 1 shows the baseline characteristics of the study population based on cola intake and sex. The numbers of non-cola drinkers and cola drinkers were 1086 and 136 males and 1113 and 164 females. The differences in most clinical variables such as age, body composition parameters, lifestyle-related factors and serum 25(OH)D levels between cola drinkers and non-cola drinkers were not significant in both sexes, but differences in nutritional factors based on cola intake were significant. Daily energy, fat, carbohydrate and protein intake were significantly higher in male cola drinkers, and daily fat and protein intake were significantly higher in female cola drinkers than in their non-cola drinking female and male counterparts, respectively. However, calcium intake was significantly lower in cola drinkers than in non-cola drinkers in males, which was not observed in females. In addition, no significant differences were observed in the daily phosphorus intake based on cola intake status in either sex.
Table 1. Baseline characteristics of the study population (n 2499) stratified by cola intake and sex
se, standard error.
Other beverage consumption and calcium intake according to cola consumption
Table 2 shows the frequency of coffee and milk consumption based on cola intake and sex. Among males, 26·9 % cola drinkers v. 18·3 % non-cola drinkers and among females, 35·8 % cola drinkers v. 26·0 % non-cola drinkers drank coffee more than four times a week. There were was no statistically significant difference in coffee intake among the groups. In addition, among males, 35·5 % cola drinkers v. 53·8 % non-cola drinkers and among females, 34·3 % cola drinkers v. 43·7 % non-cola drinkers drank milk more than four times a week; the differences were statistically significant. However, no significant difference in the daily amount of milk consumption was observed in either sex based on cola intake (data not shown).
Table 2. Beverage consumption in the study population (n 2499) stratified by cola intake and sex
se, standard error.
Table 3 shows the daily amount of cola consumption according to daily calcium intake. The cut-off values of calcium intake in cola drinkers and non-cola drinkers were, respectively, 363·7 and 608·5 mg in males and 280·2 and 499·1 mg in females. Males who consumed less calcium among tertiles showed higher cola consumption on comparing the lowest tertile group v. the middle tertile group and the middle tertile group v. the highest tertile group. However, statistical significance was not observed among those associations. Females in the lowest tertile of calcium consumption had higher cola consumption than those in the middle and highest tertiles. However, the association between calcium intake and cola consumption was not statistically significant in female cola drinkers and non-drinkers.
Table 3. Daily amount of cola consumption according to calcium intakea in each sex
se, standard error.
a T1, lowest tertile group; T2, middle tertile group; T3, highest tertile group. Cut-off values were 363·7 and 608·5 mg/d in males and 280·2 and 499·1 mg/d in females, respectively.
Relationship between cola intake and BMD
Table 4 shows the association between cola intake and BMD at various sites. When adjusted for age and BMI only, the association between cola intake and BMD was not significant at all the measured sites. However, when lifestyle-related factors were adjusted more, the association between cola intake and BMD was inversely significant for the whole body in males. After adjusting the possible confounding factors regarding nutrients, male cola drinkers tended to have significantly lower BMD at the whole body, whole femur and femoral neck with 4, 5 and 5 %, respectively, compared with that of non-cola drinkers by applying the exponential function to the presented β. Lumbar spine BMD was also lower in male cola drinkers than in non-cola drinkers. However, statistical significance was not observed in this inverse association. In females, whole body BMD was inversely associated with cola intake when all the possible confounding factors were adjusted for, although no statistical significance was observed. In addition, no significant association was observed at other skeletal sites in females. When all confounding factors were adjusted as for model 3 except phosphorus, the coefficients were not found to be significantly different in all measured skeletal sites (data not shown).
Table 4. Multivariate-adjusteda association between cola intake and bone mineral densityb at various sites
CI, confidence interval.
a Model 1: adjusted for age and body mass index.Model 2: adjusted for age, body mass index, physical activity, smoking and alcohol consumption.Model 3: adjusted for age, body mass index, physical activity, smoking, alcohol consumption, daily calcium intake, daily phosphorus intake, daily protein intake, serum 25(OH)D level and total daily energy intake. Oestrogen exposure was additionally adjusted for females.
b Log-BMD was used for multiple linear regression analysis.
Discussion
In the present study, we observed that male cola drinkers had lower BMD at the whole body, whole femur and femoral neck than non-cola drinkers after adjusting for possible confounding factors. We also found that cola drinkers had less frequent milk consumption than non-cola drinkers, but the average cola consumption was not statistically different according to calcium intake in both sexes.
In an observational study of Northern Irish adolescents aged 12 and 15 years, a significant inverse association between total carbonated soft drink intake and BMD at the dominant heel was observed in girls. However, no consistent associations were found in boys(Reference McGartland, Robson and Murray5). Similarly, a significant association between cola intake and increased risk of fracture was also reported in girls(Reference Wyshak7). In addition, a significant association between cola intake and lower hip BMD was reported in adults, which was confined to female participants(Reference Tucker, Morita and Qiao10). In agreement with these previous studies, we observed that BMD was inversely associated with cola intake. We also observed various inverse associations between cola intake and BMD depending on the measured skeletal sites, which was consistent with previous findings(Reference McGartland, Robson and Murray5,Reference Tucker, Morita and Qiao10) . Considering that each bone is composed of different materials and the trabecular bone is more metabolically active than cortical bone(Reference Rauch and Schoenau23), our findings might be influenced by various susceptibilities to nutritional influences based on the type of bone.
However, contrary to previous results, only male participants showed a significant inverse association between cola intake and BMD in the present study. When comparing the average amount of cola consumption according to sex in cola drinkers, female participants were found to consume 263·4 ± 20·7 g/d and male participants were found to consume 314·5 ± 24·9 g/d (data not shown). Since female cola drinkers consumed relatively small amounts of cola, the effects of cola on BMD may have been weakened. In addition, while the average daily calcium intake was 800–1000 mg/d in the previous studies on Western populations, in the present study, it was 544·5 ± 13·7 mg/d in males and 425·5 ± 9·6 mg/d in females, respectively (data not shown). Approximately 75 % of Korean adolescents report inadequate calcium intake, and its prevalence is increasing with age in children and adolescents(Reference Im, Kim and Lee24). The aforementioned previous studies examined subjects with sufficient calcium intake; therefore, the results of studies involving subjects with inadequate calcium intake may show different results. In addition, due to the lack of information in the original data, we were unable to assess hormonal interactions in the male participants. Although we adjusted for oestrogen exposure in females, the absence of full assessment of pubertal status in both sexes might weaken the associations observed in the present study. However, considering BMD increases in proportion to age(Reference Moretto de Oliveria, Cristiane da Silva and Kurokawa25) and reaches a plateau after completion of growth at about 16–18 years of age in females and 17–20 years of age in males(Reference Davies, Evans and Gregory26), the effect of pubertal status on BMD might have been adjusted to some degree by adjustment for chronological age. Additional research is needed to examine the different associations between cola intake and BMD based on sex and pubertal status.
We found that cola intake was associated with low BMD. However, we could not investigate how it is directly linked to low BMD. Some hypotheses have been proposed to explain the relationship between cola intake and BMD. Phosphoric acid and caffeine in cola beverages have been suggested to have deleterious effects on the bone(Reference Wyshak and Frisch12,Reference Amato, Maravilla and Montoya14,Reference Rapuri, Gallagher and Kinyamu15,Reference Fernando, Martha and Evangelina27) . However, no attenuated inverse association between cola intake and BMD was reported after adjusting for the calcium to phosphorus ratio or total caffeine intake(Reference Tucker, Morita and Qiao10), which suggests that other components of cola might play a role in this regard. In addition, since there were no significant associations between total intake of carbonated beverages or carbonated beverages except cola and BMD (data not shown), it can be said that certain components only contained in cola affect the bone.
Decreased milk consumption was found to be related to increased carbonated beverage intake in Irish adolescents(Reference McGartland, Robson and Murray5). Similar to these previous findings, we found a significant difference in milk consumption between cola drinkers and non-cola drinkers. However, there was no difference in the daily amount of milk consumption and phosphorus intake between cola drinkers and non-cola drinkers in either sex (data not shown). Furthermore, the inverse associations between cola and BMD did not differ with adjustment for phosphorus intake. Regarding the effect of phosphoric acid from cola and other foods, we could not determine the source of phosphorus since the KNHANES recorded only the total daily amount of phosphorus intake and not the sources. Based on our analysis after adjusting for daily total phosphorus intake, the inverse association observed between cola consumption and BMD did not seem to originate from phosphoric acid.
In the present study, calcium consumption was significantly lower in cola drinkers, which is consistent with previous findings(Reference McGartland, Robson and Murray5,Reference Tucker, Morita and Qiao10) . Although the statistical significance was not observed in females (P = 0·09), cola drinkers of both sexes consumed less calcium than non-cola drinkers. However, the average amount of cola consumption was not different among groups based on daily calcium intake and sex. Further studies are needed to elucidate whether specific components in cola play a role in lowering BMD and whether certain dietary habits are linked to cola consumption.
In the present study, the average amount of cola consumption was 35 g/d in female participants, and 40 g/d in male participants (data not shown), which was lower than the 230 g/d in girls and 320 g/d in boys, as reported in the Irish study on 12- and 15-year-olds(Reference McGartland, Robson and Murray5). The relatively low amount of cola consumption found in the present study might have attenuated the relationship between cola and BMD. However, the consumption of carbonated soft drinks has been steadily increasing in South Korea, reaching 29·6 l of annual consumption per capita in 2014(17). This indicates that the Korean population consumes an average of approximately 0·4 cup of carbonated soft drink per day. In addition, the average frequency of carbonated soft drink consumption in South Korean adolescents and young adults is 1–3 times a week with a tendency to decrease with age, suggesting that adolescents are the main consumers of carbonated beverages in South Korea(28). In a survey of males and females in their teens to 40s, cola was ranked the primarily purchased soda among all carbonated soft drinks (42 %), followed by cider, fruit-flavoured soda and other sodas(17). Moreover, Korean cola consumption per capita exceeded the worldwide consumption in 2012, according to a cola company report(29). Considering the trends of increase in the consumption of carbonated soft drinks including colas in South Korea, the inverse association between cola intake and BMD observed in the present study may be a matter of concern.
To the best of our knowledge, this is the first study to investigate the association between cola intake and BMD in South Korea. However, the causal relationship between cola intake and BMD could not be assessed because of the cross-sectional nature of the study. In addition, since we used self-reported questionnaires and 24-h dietary recall data, we were unable to investigate and quantify the total lifetime cola consumption. Moreover, the average amount of beverage consumption, including cola, in the present study, might not represent the usual consumption of the study participants. People who occasionally drink cola might have not had cola on the day of the examination, which would attenuate the association between cola and BMD found in the present study. In addition, we did not classify the colas according to the degree of caffeine content, which made it impossible to assess the effect of caffeine on BMD. Moreover, the hidden confounding effects of certain medications, such as corticosteroids, or certain diseases that were not excluded, such as hypogonadism or parathyroid disorders, might have affected the results. We also could not assess parathyroid hormone (PTH) levels in our study population since PTH levels were only assessed among those over 50 years of age in the KNHANES. With regard to statistical methods, we examined an extensive series of covariates affecting bone health to minimise the effects of possible confounders. However, this might produce a misleadingly high explanatory power of the regression model since several covariates were adjusted for in the analysis. Despite these limitations, since the present study was based on nationwide and population-based data in South Korea, our findings can be applied to the general population of South Korean adolescents and young adults.
In conclusion, the present study demonstrates that cola intake was inversely associated with BMD in Korean male adolescents and young adults. Considering the increase in the trend of consumption in carbonated beverages in South Korea, further studies are required to elucidate the causality between carbonated beverages and lower BMD and whether the mechanism underlying the effect of carbonated beverages on bone density is mediated by sex.
Acknowledgements
We thank professor Yun Mi Song and professor Jungkwon Lee in the Department of Family Medicine, Sungkyunkwan University School of Medicine for their comments that greatly improved the manuscript.
This research received no specific grant from any funding agency, commercial or not-for-profit sectors.
J.H.Y. developed the research questions and made substantial contributions to the study design. Y.A.K. performed the main analyses and contributed to manuscript writing and submission. All authors were involved in the review and revision of the manuscript and approved the final manuscript.
The authors declare that they have no conflicts of interest.
Open access
