Associations between cola consumption and bone mineral density in Korean adolescents and young adults: a cross-sectional study using data from the Korea National Health and Nutrition Examination Survey, 2008-2011.

The present study aimed to elucidate the relationship between cola consumption and bone mineral density (BMD) in Korean adolescents and young adults. We used data from the Korea National Health and Nutrition Examination Survey 2008-2011. A total of 2499 adolescents and young adults aged 12-25 years were included. The study participants were classified as cola drinkers and non-cola drinkers according to 24-h dietary recall data. BMD was measured using dual X-ray absorptiometry. In the male population, whole body, whole femur and femoral neck BMD in cola drinkers were lower than that of non-cola drinkers by 4% (95% CI -0⋅071, -0⋅007), 5% (-0⋅092, -0⋅012) and 5% (-0⋅090, -0⋅001), respectively. In both sex groups, cola drinkers had less frequent milk consumption than non-cola drinkers. However, there were no significant differences in cola consumption according to calcium intake in both sexes. In conclusion, cola intake and BMD were inversely associated with Korean male adolescents and young adults. Considering the importance of peak bone mass attainment at adolescents and the increasing trend in carbonated beverage consumption in South Korea, further studies are needed to elucidate the causality between cola intake and lower BMD.

Introduction

Low bone mineral density (BMD), whether osteopenia or osteoporosis, has become a worldwide interest as the population is ageing globally(. 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(. Thus, both ageing and PBM attained during childhood and adolescence can be considered predictive factors of bone health at an older age(.

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(. Several studies of children and young adults have also proposed that consumption of carbonated soft drinks is associated with low BMD( and even increased risk of fractures(. In one experimental study, a significant delay in alveolar bone formation was observed in cola consuming rats(. Furthermore, colas, but not other carbonated beverages, are related to low BMD in older women( or increased upper limb fractures in children(. However, the results of the previous studies were different according to sex(, and the association between cola and BMD has been found mainly in women(. Displacement of beneficial nutrient-rich foods such as milk(, decreased physical activity owing to watching television or computer while drinking colas(, and the presence of phosphoric acid and caffeine in colas that interfere with calcium absorption and excretion( have been suggested as possible mechanisms. However, a clear mechanism has not yet been proposed, and a set of contributing factors may exist(.

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(, 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(.

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(. 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(, 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(. 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

	Males					Females
	Non-cola drinker (n 1086)		Cola drinker (n 136)		P-value	Non-cola drinker (n 1113)		Cola drinker (n 164)		P-value
	Mean or proportion	se	Mean or proportion	se		Mean or proportion	se	Mean or proportion	se
Age (yr)	18·7	0·2	19·4	0·4	0·09	18·7	0·2	18·6	0·4	0·74
Height (cm)	171·9	0·3	172·8	0·7	0·29	160·5	0·2	160·5	0·5	0·97
Weight (kg)	66·0	0·6	66·6	1·4	0·66	54·2	0·4	55·2	1·2	0·43
Waist circumference (cm)	76·1	0·4	76·6	1·0	0·66	69·4	0·3	70·4	1·0	0·38
Body mass index (kg/m2)	22·2	0·2	22·3	0·4	0·86	21·0	0·1	21·4	0·5	0·42
Lena body mass (kg)	50·9	0·4	51·7	0·9	0·43	36·2	0·2	36·3	0·5	0·92
Bone mineral density (g/cm2)
Whole body	0·990	0·005	0·980	0·010	0·46	0·890	0·004	0·890	0·010	0·64
Whole femur	0·979	0·005	0·975	0·015	0·82	0·879	0·004	0·889	0·001	0·39
Femoral neck	0·885	0·005	0·880	0·016	0·75	0·769	0·004	0·789	0·010	0·13
Lumbar spine	0·926	0·006	0·931	0·014	0·74	0·921	0·005	0·947	0·013	0·07
Physical activity (%)					0·77					0·39
Insufficient exercise	55·2	3·1	58·1	9·6		69·1	3·0	76·2	7·2
Sufficient exercise	44·8	3·1	41·9	9·6		30·9	3·0	23·8	7·2
Ever smoker (%)	39·6	1·9	36·6	6·0	0·64	12·5	1·2	16·6	3·2	0·24
High-risk alcohol consumption (%)	9·8	1·3	9·3	3·0	0·88	2·4	0·5	1·7	1·1	0·61
Serum 25(OH)D (ng/ml)	16·4	0·2	16·3	0·7	0·94	15·0	0·3	14·4	0·5	0·35
Daily energy intake (kcal)	2356·4	48·1	2774·8	122·8	<.001	1767·4	29·2	1942·8	38·2	0·05
Daily fat intake (g)	62·8	1·9	78·8	5·8	0·01	46·6	1·3	55·8	3·6	0·02
Daily carbohydrate intake (g)	344·2	6·1	379·6	14·8	0·03	271·9	4·5	276·3	10·1	0·69
Daily protein intake (g)	86·3	2·2	100·7	5·6	0·02	62·2	1·2	70·7	3·8	0·05
Daily calcium intake (mg)	557·0	16·3	478·6	25·3	0·01	431·4	10·6	393·7	19·0	0·09
Daily phosphorus intake (mg)	1339·9	28·1	1344·8	61·2	0·94	987·5	17·0	986·3	45·3	0·98

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

	Males					Females
	Non-cola drinker (n 1086)		Cola drinker (n 136)		P-value	Non-cola drinker (n 1113)		Cola drinker (n 164)		P-value
	%	se	%	se		%	se	%	se
Coffee consumption					0·25					0·28
Rarely	42·4	2·0	31·8	5·6		41·5	2·0	34·7	5·0
1–3 times/month	17·0	1·6	18·6	4·1		13·8	1·2	11·8	3·0
1–3 times/week	22·3	1·7	22·7	4·5		18·7	1·5	17·7	3·3
4–7 times/week	14·0	1·5	19·6	4·2		17·6	1·5	26·5	5·2
≥twice a day	4·3	0·9	7·3	2·7		8·4	1·1	9·3	2·8
Milk consumption					<0·001					0·02
Rarely	5·0	0·9	13·0	3·7		13·8	1·5	9·4	2·9
1–3 times/month	10·1	1·3	6·0	2·0		10·8	1·2	9·5	2·4
1–3 times/week	31·1	2·1	45·5	5·8		31·8	1·8	46·9	4·1
4–7 times/week	46·0	2·0	28·2	5·0		39·1	1·9	32·4	4·1
≥twice a day	7·8	1·0	7·3	2·6		4·6	0·7	1·9	1·2

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 intake in each sex

	Daily cola consumption (g)		P for trend
	Mean	se
Males (n 1222)			0·62
T1	54·8	9·6
T2	52·9	10·6
T3	42·5	9·3
Females (n 1277)			0·56
T1	46·6	7·0
T2	36·5	5·9
T3	40·4	10·5

se, standard error.

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-adjusted association between cola intake and bone mineral density at various sites

	Whole body		Whole femur		Femoral neck		Lumbar spine
	β	95 % CI	β	95 % CI	β	95 % CI	β	95 % CI
Males (n 1222)
Model 1	−0·011	−0·031, 0·008	−0·013	−0·043, 0·016	−0·015	−0·048, 0·018	−0·007	−0·035, 0·020
Model 2	−0·032	−0·061, −0·002	−0·029	−0·072, 0·014	−0·026	−0·074, 0·023	−0·015	−0·058, 0·028
Model 3	−0·039	−0·071, −0·007	−0·052	−0·092, −0·012	−0·046	−0·090, −0·001	−0·027	−0·071, 0·017
Females (n 1277)
Model 1	0·001	−0·013, 0·016	0·004	−0·017, 0·025	0·018	−0·008, 0·044	0·020	−0·002, 0·037
Model 2	0·006	−0·020, 0·032	0·000	−0·034, 0·035	0·009	−0·027, 0·046	0·025	−0·005, 0·055
Model 3	−0·008	−0·043, 0·028	0·019	−0·018, 0·056	0·021	−0·033, 0·076	0·023	−0·011, 0·056

CI, confidence interval.

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.

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(. Similarly, a significant association between cola intake and increased risk of fracture was also reported in girls(. In addition, a significant association between cola intake and lower hip BMD was reported in adults, which was confined to female participants(. 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(. Considering that each bone is composed of different materials and the trabecular bone is more metabolically active than cortical bone(, 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(. 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( 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(, 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(. However, no attenuated inverse association between cola intake and BMD was reported after adjusting for the calcium to phosphorus ratio or total caffeine intake(, 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(. 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(. 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(. 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(. 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(. 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(. Moreover, Korean cola consumption per capita exceeded the worldwide consumption in 2012, according to a cola company report(. 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.