Introduction
Osteoporosis (OP) is a common chronic skeletal disorder characterised by decreased bone mineral density (BMD) and deterioration of bone microarchitecture.(Reference Costa, da Silva and Brito1,Reference Rahman, Usman and Sheikh2) With the increase in life expectancy and the aging population worldwide, OP brought a significant financial burden and public health challenge.(Reference Shen, Chyu and Wang3,Reference Jaul and Barron4) The prevalence of OP is increasing in the United States.(Reference Ayers, Kansagara and Lazur5) Estrogen levels decline has been reported to contribute to the pathogenesis of BMD reduction.(Reference Tabatabaei-Malazy, Salari and Khashayar6) In women, with the reduction in estrogen production caused by postmenopausal, the secretion of osteoprotegerin (OPG) reduces, which leads to skeletal disorder and deterioration in BMD, and is associated with an increased risk of fracture, and further results in an increasing mortality.(Reference Tey, Chew and How7,Reference Alejandro and Constantinescu8) Therefore, the risk factors for BMD reduction in postmenopausal women are more worthy to be explored and identified.
Caffeine is a bioactive compound in coffee, and the relationship between caffeine consumption and human health has been widely reported.(Reference Cappelletti, Piacentino and Sani9) Studies indicated that caffeine intake negatively mediated bone homeostasis to cause bone loss and even OP, and the possible mechanism may be increased urinary calcium (Ca) excretion and reduced endogenous Ca absorption, leading to a negative calcium balance.(Reference Wang, Zhou and Guan10,Reference Brun, Brance and Lombarte11) Postmenopausal women as the major population of Ca loss at particularly increased risk for OP,(Reference Hallstrom, Wolk and Glynn12) and however, views on the relationship between caffeine consumption and BMD in postmenopausal women have still been inconsistent. Rapuri et al. (Reference Rapuri, Gallagher and Kinyamu13) found that caffeine intake can accelerate bone loss at the spine in elderly postmenopausal women, while other studies found that there was no significant overall association between high v. low consumption coffee and risk of fracture in women.(Reference Poole, Kennedy and Roderick14,Reference Hallstrom, Byberg and Glynn15)
Clinical trials are the gold standard in exploring relationship between caffeine and BMD, and however, they are usually impractical due to the caffeine metabolism status cannot be assessed in a short period of time. Metabolomics can measure the metabolites in biofluids and may capture the metabolic products of food and better reflect the dietary exposure after metabolism, which may be related to disease pathogenesis.(Reference Yuliana, Hunaefi and Goto16) Recently, Chau et al. (Reference Chau, Au and Li17) identified twelve serum metabolites were significantly associated with self-reported habitual coffee intake in community-dwelling Chinese adults, and four of these twelve coffee-associated metabolites were significantly associated with BMD. However, there are few studies focus on the association between caffeine metabolites and BMD in postmenopausal women.
Given the equivocal effect of caffeine and its metabolites on the bone, this study aims to explore the relationship between caffeine consumption and caffeine metabolites and BMD in postmenopausal women, in order to provide some references for early identification of high-risk populations and prevention of OP.
Methods
Study design and populations
Data in this cross-sectional study were extracted from the National Health and Nutrition Examination Survey (NHANES) database in 2009–14. NHANES is an ongoing study conducted by the Centers for Disease Control and Prevention (CDC) to assess the nutritional and health status of the non-institutionalized population in the United States. The regular data collection of the NHANES is carried out on approximately 5000 persons from 15 areas since 1999 and examined in 2-year periods. The NHANES is a multi-stage stratified sampling database, and data used in our statistical analyses were available on a public link address: https://wwwn.cdc.gov/nchs/nhanes/. Interviews in participants’ homes were conducted by trained professionals from the National Center for Health Statistics (NCHS), and extensive physical examinations were conducted at mobile exam centres (MECs).
A total of 4286 postmenopausal women who diagnosed with OP or received anti-osteoporosis therapy were initially included. The exclusion criteria were missing the information of caffeine consumption or caffeine metabolites, and having an extreme energy intake (<500 or >5000 kcal/d). Finally, 4042 of them were eligible for caffeine analyses and that 583 for caffeine metabolites analyses. The database protocol was approved by the institutional review board (IRB) of NCHS. Since the database was publicly available, the requirement of ethical approval for this study was waived by the IRB of the 958th Hospital of the PLA Army.
Measurement of BMD
In NHANES, total body scans of participants were performed for BMD (gm/cm2) measurement by the fast mode using dual-energy X-ray absorptiometry (DXA) with Hologic QDR 4500A fan-beam densitometers (Hologic Inc., Bedford, MA, USA).(Reference Baker, Weber and Neogi18) Details of DXA examination protocol were shown in the Body Composition Procedure Manual of the NHANES.(Reference Kaibori, Kawaguchi and Yokoigawa19) BMD levels in the database were measured based on the lumbar spine and femoral regions of total femur, femur neck, trochanter, and intertrochanter. In the current study, we used the femur neck BMD. The average of first through fourth lumbar vertebra was calculated as BMD of lumbar spine.(Reference Khalil, Chen and Lee20) The left hip was routinely scanned for regions of proximal femur. The right hip was scanned when the participant self-reported a left hip replacement, a fractured left hip, or a pin in the left hip. Participants who had fractures, replacements or pins in both hips were excluded from the DXA scan.(Reference Lim, Kim and Lim21) Women who had a BMD z-score ≥1 and ≤ −1⋅28 measured at either lumbar spine or total hip, respectively, were diagnosed with low BMD.
Dietary caffeine intake assessment
Dietary caffeine intake was collected via the 24-h dietary recalls, in which participants reported individual foods and drinks consumed during the midnight-to-midnight 24-h period prior to the in-person dietary interview. Investigators of the two 24-h dietary recalls interviews were well trained. The first dietary recall interview is conducted in-person in the MEC, and the second interview is by telephone 3–10 d after the first one. NHANES conducted the coding of interview data and conversion to total nutrient intakes by using the USDA Food and Nutrient Database for Dietary Studies, 5.0 (FNDDS 5.0) (http://www.ars.usda.gov/ba/bhnrc/fsrg). The FNDDS 5.0 nutrient values were based on the USDA National Nutrient Database for Standard Reference, release 24 (http://www.ars.usda.gov/nutrientdata).
Caffeine metabolites detection
Detailed urine samples collection and processing instructions are presented in the NHANES Laboratory/Medical Technologists Procedures Manual.(Reference Sato, Murai and Ueda22) Utilised an ultra-high-performance liquid chromatography coupled with tandem mass spectrometry for quantitative analysis of caffeine and fourteen of its metabolites, including MethyluricAcid1, MethyluricAcid3, MethyluricAcid7, MethyluricAcid13, MethyluricAcid17, MethyluricAcid37, MethyluricAcid137, Methylxanthine1, Methylxanthine3, Methylxanthine7, Methylxanthine13, Methylxanthine17, Methylxanthine37, and Methylxanthine137, were measured for the participants.
Variables collection
We collected variables from the NHANES database including age, poverty income ratio (PIR), race, educational level, smoking, drinking, vigorous activity, moderate activity, height, weight, body mass index (BMI), Ca intake and supplements, vitamin D (VD) intake and supplements, energy intake, protein intake, sugar intake, fat intake, carbohydrate intake, total cholesterol (TC), low-density lipoprotein (LDL), high-density lipoprotein (HDL), triglyceride (TG), diabetes mellitus (DM), hypertension, coronary heart disease (CHD), angina pectoris, stroke, previous fracture, parental history of OP, glucocorticoids, bone resorption inhibitors, estrogens, corticosteroids, OP treatment, and hypercholesterolaemia.
Physical examinations were conducted and the information was collected in the MECs. Laboratory examinations used the blood samples of participants, and the methodological details of the laboratory analyses have been described on the NHANES website. DM was defined as fasting blood glucose ≥7⋅0 mmol/l or glycosylated haemoglobin (HbAlc) ≥6⋅5 % or self-reported DM or receiving hypoglycaemic therapy. Hypertension was defined as a mean blood pressure exceeding 140/90 mmHg for systolic pressure and diastolic pressure, respectively. Other chronic diseases were diagnosed according to the self-report or the medications.
Statistical analyses
Non-normal distributed data were described by median and quartiles [M (Q1, Q3)] and the Mann–Whitney U rank test was used for comparison. Categorical data were expressed as frequency and constituent ratio [N (%)], and the chi-square (χ 2) test was used for comparison. The 2-year MEC exam weight (wtmec2yr) was used, which is needed for all NHANES analyses of 2009–10, 2011–12, and 2013–14. Detailed instructions for combining datasets from NHANES cycles are provided in the NHANES Analytic Guidelines.
Weighted univariate linear regression and stepwise regression analyses were used to screen the covariates. Weighted univariate and multivariate linear regression analyses were used to explore the association between caffeine intake and caffeine metabolites and BMD. The evaluation index was estimated value (β) with 95 % confidence intervals (CIs). Model 1 was the crude model. Model 2 adjusted for age, race, PIR, and BMI. Model 3 adjusted for age, race, PIR, BMI, vigorous activity, moderate activity, DM, previous fracture, VD intake, bone resorption inhibitors, estrogens, parental history of OP, and total energy intake. Subgroup analysis of age was also performed to explore these relationships. P < 0⋅05 was considered significant.
Statistical analyses were using SAS 9.4 (SAS Institute., Cary, NC, USA) and R 4.0 (Math Soft, Seattle, WA, USA). Missing variables were managed by multiple imputation using the ‘mice’ package by R software.
Result
Characteristics of postmenopausal females
Fig. 1 shows the flowchart of participants screening. A total of 4286 postmenopausal females with OP or received anti-osteoporosis therapy were initially included. We excluded those who missing the information of caffeine intake (n 123) or caffeine metabolites measurement (n 3689), and have extreme energy intake (n 121) for caffeine analyses, n 14 for caffeine metabolites analyses. Finally, 4042 of them were eligible for the caffeine analyses and that 583 for caffeine metabolites analyses.
Flowchart of data screening.
The characteristics of postmenopausal females are shown in Table 1. For the caffeine group, the median age was 61⋅1 years old, and the median caffeine intake was 130⋅3 mg. The median BMD was 0⋅7 gm/cm2. For the caffeine metabolites group, the median age was 59⋅9 years old, and the median caffeine intake was 135⋅7 mg. Similarly, the median BMD was 0⋅7 gm/cm2.
Characteristics of postmenopausal women
M, median; Q1, 1st quartile; Q3, 3rd quartile; PIR, poverty income ratio; BMI, body mass index; Ca, calcium; VD, vitamin D; BMD, bone mineral density; TC, total cholesterol; LDL, low-density lipoprotein; HDL, high-density lipoprotein; TG, triglyceride; T-score, T-score = (individual BMD – reference BMD)/reference sd; DM, diabetes mellitus; CHD, coronary heart disease; OP, osteoporosis.
Associations between caffeine intake and caffeine metabolites and BMD
Fig. 2 shows fit curves of the caffeine intake and caffeine metabolites and BMD. The difference sum between each point and the regression line is the minimum, indicating an approximate linear trend between caffeine intake and caffeine metabolites and BMD.
The fit curves of caffeine and its metabolites intake and BMD.
Tables 2 and 3 respectively shows the relationships between caffeine intake and caffeine metabolites and BMD. After adjusting for the covariates, we found that caffeine intake was not significantly associated with BMD [β = −0⋅012, 95 %CI (−0⋅049, 0⋅026)]. However, caffeine metabolites, including MethyluricAcid3 (β = −4⋅036 × 10−3), MethyluricAcid7 (β = −0⋅199 × 10−3), MethyluricAcid37 (β = −3⋅995 × 10−3), Methylxanthine3 (β = −0⋅150 × 10−3), and Methylxanthine37 (β = −0⋅379 × 10−3), were negatively associated with BMD.
Association between caffeine intake and BMD
PIR, poverty income ratio; BMI, body mass index; Ca, calcium; VD, vitamin D; BMD, bone mineral density; CI, confidence interval; TC, total cholesterol; LDL, low-density lipoprotein; HDL, high-density lipoprotein; TG, triglyceride; DM, diabetes mellitus; OP, osteoporosis.
Model 1: crude mode.
Model 2: adjusted for age, race, PIR, and BMI.
Model 3: adjusted for age, race, PIR, BMI, vigorous activity, moderate activity, DM, previous fracture, VD intake, bone resorption inhibitors, estrogens, parental history of OP, and total energy intake.
a The actually estimated value was that in table multiply by 10−3.
Association between caffeine metabolites and BMD
PIR, poverty income ratio; BMI, body mass index; Ca, calcium; VD, vitamin D; BMD, bone mineral density; CI, confidence interval; TC, total cholesterol; LDL, low-density lipoprotein; HDL, high-density lipoprotein; TG, triglyceride; DM, diabetes mellitus; CHD, coronary heart disease; OP, osteoporosis.
Model 1: crude mode.
Model 2: adjusted for age, race, PIR, and BMI.
Model 3: adjusted for age, race, PIR, BMI, vigorous activity, moderate activity, DM, previous fracture, VD intake, bone resorption inhibitors, estrogens, parental history of OP, and total energy intake.
a The actually estimated value was that in table multiply by 10−3.
Fig. 3 shows the residual analysis of caffeine and its metabolites, and the results show that all the points are symmetrically distributed around y = 0, and most of the points in the residual analysis diagram of caffeine/caffeine metabolites fall in the horizontal band (−0⋅5, +0⋅5) without any trend and are completely randomly distributed in the band, indicating that the regression equation adopted is a good fit for the sample data.
The residual analysis of caffeine and its metabolites.
Relationship between caffeine metabolites and BMD in age subgroups
The association between caffeine metabolites and BMD in age subgroup is shown in Table 4. After adjusting for covariates including race, PIR, BMI, vigorous activity, moderate activity, DM, previous fracture, VD intake, bone resorption inhibitors, estrogens, parental history of OP, and total energy intake, in women aged <65 years old, MethyluricAcid37 (β = −4⋅198 × 10−3) and Methylxanthine37 (β = −0⋅455 × 10−3) was negatively associated with the BMD, while in those who aged ≥65 years old, the negative relationships were found between MethyluricAcid3 (β = −4⋅890 × 10−3) and Methylxanthine3 (β = −0⋅153 × 10−3).
Association between caffeine metabolites and BMD in age subgroups
BMD, bone mineral density; CI, confidence interval.
Adjusted for race, PIR, BMI, vigorous activity, moderate activity, DM, previous fracture, VD intake, bone resorption inhibitors, estrogens, parental history of OP, and total energy intake.
a The actually estimated value was that in table multiply by 10−3.
Discussion
In this study, we explored the relationships between caffeine and its metabolites and BMD in postmenopausal women. The results showed that caffeine intake was not significantly associated with the BMD. However, caffeine metabolites, including MethyluricAcid3, MethyluricAcid7, MethyluricAcid37, Methylxanthine3, and Methylxanthine37, were negatively associated with BMD. In addition, the type of caffeine metabolites which negatively associated with BMD were different between women aged <65 years old and those who aged ≥65 years old.
High caffeine intake has been cited as a risk factor for OP among postmenopausal women,(Reference Rapuri, Gallagher and Kinyamu13) but views on the relationship between caffeine consumption and BMD in postmenopausal women are various and equally inconclusive. Hallström et al. (Reference Hallstrom, Byberg and Glynn15) demonstrated that a modestly increased risk of OP fractures was linked to a daily caffeine intake of 330 mg (equivalent to 600 ml of coffee) or more, especially in females with a low intake of Ca. Rapuri et al. (Reference Rapuri, Gallagher and Kinyamu13) found that caffeine intake amounted >300 mg/d (approximately 514 g, or 18 oz, brewed coffee) contributed to a bone loss at the spine in postmenopausal women. Oppositely, Bijelic et al. (Reference Bijelic, Milicevic and Balaban23) pointed that there was no evidence can prove that a daily coffee intake ≥3 cups is a risk factor for OP in postmenopausal women. A cross-sectional analysis among Korean postmenopausal women also showed that compared with those who had lower coffee consumption, women with higher level of coffee intake had lower odds for OP, and as coffee consumption increased, so did the BMD of femoral neck and lumbar spine.(Reference Choi, Choi and Park24) In our study, there was no significant association between caffeine intake and BMD in postmenopausal women. The median caffeine consumption of postmenopausal women was approximately 130 mg. However, the relationship between caffeine intake and BMD is not conclusive, and further exploration is still needed.
As a matter of fact, caffeine can result in an increased excretion of Ca in the urine, and that cannot be fully compensated even after 24 h.(Reference Bijelic, Milicevic and Balaban23) Caffeine consumption can also lead to a reduction in the interstitial Ca absorption.(Reference Li, Dai and Wu25) The possible mechanism between caffeine intake and BMD reduction may be that caffeine increases the urinary excretion of Ca to disturb bone metabolism, resulting in decreased stores available for bone deposition, and this derangement in Ca metabolism can alter the normal development of bone and reduce both BMD and bone volume.(Reference Lacerda, Matuoka and Macedo26) Additionally, another study demonstrated that caffeine induced inhibitory effects on the cell viability, proliferation, migration, and pluripotency of bone marrow mesenchymal stem cells (BMSCs) in vitro, and excessive caffeine could induce OP via the suppression of osteogenesis and the promotion of adipogenesis of BMSCs.(Reference Hua, Zou and Ma27) According to our results, we speculated that caffeine may affect the bone metabolism though regulating Ca levels so that contribute to a BMD reduction.
Caffeine metabolites, including MethyluricAcid3, MethyluricAcid7, MethyluricAcid37, Methylxanthine3, and Methylxanthine37, were negatively associated with BMD in this study. Similarly, Chau et al. (Reference Chau, Au and Li17) identified twelve serum metabolites significantly associated with self-reported habitual coffee intake in community-dwelling Chinese adults, and four of these twelve coffee-associated metabolites were significantly associated with BMD. Although the specific mechanisms by which these caffeine metabolites we identified affect BMD have not been reported, they have also been reported to be significantly associated with a number of other diseases. Glover et al. (Reference Glover, Caudle and Del Giudice28) found significant, inverse associations between six xanthine metabolic products of caffeine and testosterone. Ngueta et al. (Reference Ngueta29) showed that the odds of hypertension decreased across quartiles of MethyluricAcid3, MethyluricAcid7, Methylxanthine3, and Methylxanthine7. The metabolically active products of caffeine metabolism including theophylline and theobromine of the xanthine class have direct effects on gonadotropin-induced steroidogenesis.(Reference Williams, Horner and Catt30) Further studies on the underlying mechanisms of the effects of caffeine metabolites on BMD in postmenopausal woman needed to be conducted.
Subgroup analysis based on age showed that in females aged <65 years old, MethyluricAcid37 and Methylxanthine37 were negatively associated with the BMD, while in those who aged ≥65 years old, MethyluricAcid3 and Methylxanthine3 were more noticeable. Some epidemiological studies indicated that caffeine only had adverse effects on bone mass among postmenopausal women,(Reference Poole, Kennedy and Roderick14,Reference Rodan and Martin31,Reference Ilich, Brownbill and Tamborini32) but no negative association between caffeine consumption and bone mass was seen in young women.(Reference Wetmore, Ichikawa and LaCroix33,Reference Conlisk and Galuska34) Obviously, the different level of age-related estrogen is the major difference between old and young women. There has been plentiful evidence that showed estrogen had the ability to reduce bone loss,(Reference Davis, Ross and Johnson35) decrease risk of fractures,(Reference Prince, Dick and Beilby36) and increase Ca absorption(Reference Gennari and Agnusdei37) that in contrast to the effects of caffeine on bone. We presumed that the negative roles of caffeine on bone may exhibit in absence or low level of estrogen in vivo, since estrogen has been proved to have multiple positive impacts on the metabolism of bone-related cells.(Reference Wang, Yu and Zhai38–Reference Dai, Li and Quarles40) However, the role of age in the effect of caffeine metabolites including MethyluricAcid and Methylxanthine on BMD needed further explored.
Strengths and limitations
The study population was from the NHANES database so that the sample size was large and representative. There are few objective ways to measure caffeine intake, and in the present study, NHANES measured caffeine exposure through the food-frequency questionnaire (FFQ). However, individual's habitual diet could not be fully assessed since the data was completed by using a FFQ that is subject to recall bias and selective bias. This observational study was also unable to observe the dose of caffeine exposure over a long period of time. The causal relationship between caffeine consumption and caffeine metabolites and BMD in postmenopausal women cannot be determined given the study was a cross-sectional study. Further prospective cohort studies are needed to explore the roles of caffeine and its metabolites in OP among postmenopausal women.
Conclusion
Caffeine metabolites were negatively associated with BMD in postmenopausal women, and further study is still needed to find the underlying mechanisms.
Acknowledgements
This research received no specific grant from any funding agency, commercial, or not-for-profit sectors.
S. L., Y. Z., and Z. Z. designed the study. S. L. wrote the manuscript. J. Z., H. C., W. W., and F. Y. collected, analysed, and interpreted the data. Y. Z. and Z. Z. critically reviewed, edited, and approved the manuscript. All authors read and approved the final manuscript.
The authors declare that there is no conflict of interest.
Open access
