Skip to main content
  • Research article
  • Open access
  • Published:

Correlation between bone turnover and metabolic markers with age and gender: a cross-sectional study of hospital information system data

Abstract

Background

Bone turnover and metabolic indicators are related to age and gender. Age and gender should be matched in subjects in disease control research of bone turnover and metabolism, but strict matching of gender and age increases the difficulty and cost of the research. Therefore, the aim of this study was to solve it is necessary to strictly match age and gender in clinical research in bone metabolism.

Methods

A cross-sectional study was conducted from the data were extracted from the HIS of ZhuJiang Hospital. Data relating to seven bone turnover and metabolic indicators from 1036 patients between January 2018 and October 2019 were analyzed.

Results

P1NP, β-CTx and 25(OH)D were significant different in individuals younger than 20 years of age. ALP was significantly higher in those under 20 years of age and lower at age 20–39 compared with other age groups. The concentrations of Ca and P were different among the groups aged 0–19, 20–39, and 40–59 years of age groups but exhibited no difference above 60 years of age. PTH expression was not dependent on age. P1NP, β-CTx and PTH concentrations were not significantly different between the genders within the same age group. ALP was significantly different between genders within the age range 20–59 years. Ca and 25(OH)D were significantly different between the genders for those older than 60. Serum P was significantly different in the two genders for those aged 40–79. Patients received both alfacalcidol and calcium treatment differently from the others in P1NP, β-CTx, Serum Ca, P and ALP.

Conclusion

P1NP and β-CTx were highly correlated with age. If these two indictors require analysis in a case control study, the patients and controls should be strictly matched by age under 20 years. The demarcation point for ALP was 40 years of age. Ca and P were strongly recommended strict matching according to age in disease research. The difference in P1NP, β-CTx, 25(OH)D and ALP between genders depends on age differences. Medication history should be considered in bone turnover and metabolic clinical research.

Peer Review reports

Background

The skeletal system operates at different levels of metabolic activity in different disease states which also varies with age [1]. Bone metabolism indicators have been studied well in detail menopausal women, children and adolescents. The majority of researchers believe that the age and gender of subjects should be matched in bone metabolism research. It remains unclear whether research of the control of bone metabolic diseases requires the strict matching of the age and sex of patients. Also unclear is which is more important to strictly match, age or gender. Strict matching involves great difficulty and high costs, especially in randomized controlled trials (RCTs). It is unclear whether strict matching is in fact required at every age. Since no epidemiological research is able currently to answer this question, especially that of infants and very young children aged 0–5 years old, we wished to investigate such relevant information from inpatient data as a reference. We suppose that there is association between markers of bone formation or turnover with age and gender. Therefore, the aim of this study was to answer the question of matching. We assessed the relationship between markers of bone formation or turnover with age and gender, fitted the most appropriate curve and inferred the apparent relationship between them. Thus, data for seven bone indicators related to bone formation, turnover and metabolism were extracted, and their association with age and gender analyzed. The markers measured were amino-terminal propeptide of type I collagen (P1NP), beta-isomerized C-terminal telopeptide of type I collagen (β-CTx), alkaline phosphatase (ALP), 25-hydroxyvitamin D (25(OH)D), parathyroid hormone (PTH), calcium (Ca) and phosphorus (P) in serum.

Methods

Study design

A cross-sectional study was conducted to assess the association between bone turnover and metabolic indicators with age and gender in inpatients. Because this study focused on age matching for disease studies, all patient data were extracted from the hospital information system (HIS) of ZhuJiang Hospital, Southern Medical University, from January 2018 to October 2019. Seven indicators from patients were retrieved. P1NP, β-CTx, 25(OH)D and PTH were analyzed at ZhuJiang Hospital, Guangzhou, China, using the method of ECLIA (Electro chemiluminescent immunoassay), and ALP was measured by the rapid colorimetric test. Serum calcium and P were measured with the colorimetric kit, according to the manufacturer’s guidelines. No additional medical procedures were performed on patients as all medical record data of patients who had already received treatment had already been recorded and stored within the computer system. ZhuJiang Hospital of Southern Medical University Ethics Committee waived the approval due to the retrospective nature of the study.

Patients

We selected as many patients as possible from the (HIS) of ZhuJiang Hospital, Southern Medical University. From January 2018 to October 2019, the records of 1036 subjects in the age range from infant to 90 years who had undergone a biochemical indicator test of bone turnover and metabolism in ZhuJiang Hospital of Southern Medical University were obtained. The levels of bone metabolism are different in age change, so all patients were grouped by age to hierarchically analyze the bone indicators of different age brackets, including 0–19,20-39,40-59,60–79 and over 80 years of age. A total of 156 patients aged under 20 were identified and analyzed separately. Levels of P1NP, β-CTx, ALP, 25(OH)D, PTH, Ca and P in serum had been recorded in each case during routine laboratory investigation. Kinds of disease of all the patient’s data in supplement table 1. To reduce information bias, the data were collected by four experienced researchers, who divided into two teams, one team collected data and the other team checked the data.

Statistical analysis

Statistical analysis was performed using IBM SPSS v24.0 (USA) software. Data are presented as median values (QR). Curve fitting and hierarchical analysis were conducted to quantify associations among the variables. Age-related changes in bone turnover and metabolic markers were modeled using third-degree polynomial functions, which gave the highest coefficients of determination (R2) in comparison with other models. Cubic and exponential curves were plotted from the calculated results. To determine potential factors for bone turnover and metabolic markers, multivariable logistic regression analysis was performed, and the associations between influence factors and outcomes are presented as odds ratios (ORs) and 95% CIs, after adjustment for confounders, including age, gender, types of anti-osteoporosis medicines, anti-osteoporosis medicines taking duration time, types of disease and bone turnover and metabolic detection season. Comparisons between groups were calculated using Chi-square, Mann-Whitney and Kruskal-Wallis tests. Use Max Likelihood analysis where missing data were taken into account. All statistical tests and confidence intervals were two-sided. P < 0.05 was considered statistically significant.

Results

For all patients, the seven bone turnover and metabolic markers namely P1NP, β-CTx, ALP, Ca, P, 25(OH)D and PTH were extracted then analyzed. All baseline characteristics are presented in Table 1.

Table 1 Bone turnover and metabolic indicators are displayed hierarchically by age

Correlation with age

Young patients had high levels of P1NP, β-CTx and 25(OH)D, the difference in their concentration in patients younger than 20 years compared with those who were older was significant, but not significant between groups for those over 20 years of age (Fig. 1a, b and g; Table 1). There was a significant difference in ALP concentrations for patients older than 20 years compared with those younger than 20 years age. ALP was significantly lower in patients aged 20–39 compared with other age groups. ALP concentration did not vary in patient groups over the age of 40 years (Fig. 1c; Table 1). The concentrations of Ca and P in different age groups were completely inconsistent, with differences among the 0–19, 20–39, and 40–59 years age groups, but not different over 60 years of age (Fig. 1e, f; Table 1). PTH expression was not dependent on age (Fig. 1d; Table 1).

Fig. 1
figure 1

a Different expression of serum P1NP concentration in different ages; b Different expression of serum β-CTx concentration in different ages; c Different expression of serum ALP concentration in different ages; d Different expression of serum PTH concentration in different ages; e Different expression of serum Ca concentration in different ages; f Different expression of serum P concentration in different ages; g Different expression of serum 25(OH)D concentration in different ages; p < 0.05 as assessed by Mann–Whitney U test

Because bone turnover and metabolic indicators were significantly different for those under 20 years of age, we analyzed and modeled those patients separately. In a total of 156 young patients, the relationship between P1NP and age was not a simple linear correlation. Regression analysis and curve fitting suggested that the relationship between P1NP and age was cubic polynomial. P1NP decreased with age below the age of 3–4, then increased over the range 4–14 years. Over14 year of age, P1NP then again declined with increasing age (Fig. 2a, b). β-CTx and ALP exhibited trends similar to those of P1NP, levels declining before the age of 4, then gradually increasing from 4 to 14 years, and gradually decreasing after 14 years until reaching the levels observed in adults (Fig. 3c, d, k, l). Serum Ca and P levels also fitted a cubic polynomial curve, although R2 was not greater than 0.3. (Fig. 2e-h). PTH and 25(OH)D levels did not exhibit a similar declining model but had a small rising trend with age (Fig. 2i-l).

Fig. 2
figure 2

B Box plots of serum P1NP concentration by age and gender in patients under 20 years; b Cubic polynomial curves of serum P1NP concentration by age and gender in patients under 20 years; c Box plots of serum β-CTx concentration by age and gender in patients under 20 years; d Cubic polynomial curves of serum β-CTx concentration by age and gender in patients under 20 years; e Box plots of serum Ca concentration by age and gender in patients under 20 years; f Cubic polynomial curves of serum Ca concentration by age and gender in patients under 20 years; g Box plots of serum P concentration by age and gender in patients under 20 years; h Cubic polynomial curves of serum P concentration by age and gender in patients under 20 years; i Box plots of serum 25(OH)D concentration by age and gender in patients under 20 years; j Cubic polynomial curves of serum 25(OH)D concentration by age and gender in patients under 20 years; k Box plots of serum ALP concentration by age and gender in patients under 20 years; l Cubic polynomial curves of serum ALP concentration by age and gender in patients under 20 years; m Box plots of serum PTH concentration by age and gender in patients under 20 years; n Exponential curve of serum PTH concentration by age in female patients under 20 years; p < 0.05 as assessed by Regression analysis

Fig. 3
figure 3

a Different expression of serum P1NP concentration in different ages and genders; b Different expression of serum β-CTx concentration in different ages and genders; c Different expression of serum ALP concentration in different ages and genders; d Different expression of serum PTH concentration in different ages and genders; e Different expression of serum Ca concentration in different ages and genders; f Different expression of serum P concentration in different ages and genders; g Different expression of serum 25(OH)D concentration in different ages and genders; p < 0.05 as assessed by Mann–Whitney U and Kruskal-Wallis tests

Correlation with gender

A stratified analysis was conducted according to the different age groups. The results indicated that serum P1NP, β-CTx and PTH were not significantly different between the two genders (Table 2, Fig. 3a, b and d). Serum ALP levels were significantly different for different genders for those aged 20–59 (Table 2, Fig. 3c). Serum Ca and 25(OH)D were significantly different in the two genders for patients older than 60 (Table 2, Fig. 3e, g). Serum P was significantly different in the two genders for those aged of 40–79 (Table 2, Fig. 3f). We additionally analyzed the bone markers of boys and girls aged between 0 to 19 years old. The age cut-off points for pre-pubertal and pubertal of boys are 12 years old and girls are 11 years old. We found only the levels of β-CTx (p = 0.006, Bonferroni Correction) has significant differences between pubertal boys and girls.

Table 2 Bone turnover and metabolic indicators are displayed hierarchically by gender and age

Correlation with anti-osteoporosis medication

Multiple logistic regression analyses found either type of anti-osteoporosis medicines or anti-osteoporosis medicines taking duration time are not the main factors to influence bone turnover and metabolic markers in these patients. However, nonparametric analysis (Bonferroni Correction) found only the group in which patients received both alfacalcidol and calcium treatment has difference than the other groups in P1NP, β-CTx, Serum Ca, P and ALP (Fig. 4a-e).

Fig. 4
figure 4

a Patients received both alfacalcidol and calcium treatment has a difference than other groups in P1NP; b Patients received both alfacalcidol and calcium treatment has a difference than other groups in β-CTx; c Patients received both alfacalcidol and calcium treatment has a difference than other groups in Serum Ca; d Patients received both alfacalcidol and calcium treatment has a difference than other groups in serum P; e Patients received both alfacalcidol and calcium treatment has a difference than other groups in ALP. p < 0.05 as assessed by Mann–Whitney U and Kruskal-Wallis tests

Another found

Multiple logistic regression analyses found the 25(OH)D test in summer will get a higher level than other seasons (OR:1.70, 95% CI:1.01–2.60, p = 0.044).

Discussion

The present study focused on bone metabolism indicators in patients of different age and gender using HIS data from one hospital, providing mathematical modeling to assist the future study of age and gender in disease research. Although multiple studies have been published that provide a reference for differences in children and adolescents [2,3,4,5,6,7,8], those studies did not provide a precise breakdown of age and gender, namely at what patient age does parameters require precise matching and which age does they not. Researchers believe that age matching should be performed on all subjects in disease control research, but strict age-matching substantially increases the cost and difficulty of conducting the study. The present research found that patients aged more older than 20 years demonstrate that P1NP, β-CTx and 25(OH)D in patient research do not require age matching. This will greatly help researchers reduce the cost and difficulty of their research.

The biological data for young individuals is very complicated, and so we established mathematical models to correlate data for people under the age of 20. We noted that P1NP and β-CTx levels were especially correlated with age. The concentrations of P1NP in infants and young children were higher although their values gradually decreased with age, but not in a simple linear relationship. Firstly, it significantly declined from 0 to 5 years, followed by a period of equilibrium and a small escalation phase between the ages of 6 and 15. Over 15, P1NP levels again rapidly declined by the age of 20 to values close to those of adults. Many studies have demonstrated that P1NP and β-CTx are associated with age, but none have analyzed the relationship in detail or curve fitted the relationship, as conducted here. According to statistical analysis, the relationship with P1NP and age is a complex third-degree polynomial function, and the curve fits well using a cubic relationship. There are two peaks in bone growth in childhood. The first peak appears in infancy and the second in early adolescence. The present study found that bone metabolic indicators grew at their highest rate during puberty, the fastest rate for bone minerals over the age 12–13. As an indicator related to rapid growth rate in healthy individuals, it is not surprising that serum ALP and PINP declines after puberty. β-CTx declined before the age of 4, then increased over the ages 4–14 years, gradually decreasing after the age of 14. In the present study, because R2 was less than 0.3, the trend was judged to not be significant, possibly due to the limited sample size. A positive result might have been possible had the sample size been larger. The data obtained from the HIS in the present study included a large number of records of infants and children which allowed us to construct a mathmatical model. Prior to this study, almost no bone metabolic data had been published regarding Chinese children or adolescents, especially babies. In one large study, the minimum age was 15–19 years [9]. Another study that researched a large cohort of healthy adults in China, focused on healthy individuals older than 20 years of age [10]. Thus, studies on children in China, both healthy and sick, are scarce [11, 12]. Studies of pediatric-specific diseases are valuable as reference data. However, due to the sample size and narrow-specificity of the diseases in question, the results are often generally only applicable in a limited fashion [13,14,15]. Nevertheless, the curve model of young patients obtained in this study provides a good basis for future research.

Bone metabolism is closely related to gender, an association which has been described in multiple studies [9, 16,17,18,19]. In the present research, we found significant differences in P1NP, β-CTx, 25(OH)D and ALP concentrations between genders at all ages. However, if study subjects are stratified by age, the differences in gender are not apparent. This suggests that differences in gender depend on the ages of individuals, indicating that gender does not need to be strictly matched if age has been matched. Conversely, if age is not strictly matched, it would be necessary to strictly match by gender. Clearly, in clinical research, it is easier to match gender than age, so matching gender only would be advantageous for future research studies.

The patients who both took alfacalcidol and calcium were analyzed separately. Patients taking both two medication have lower N1NP, β-CTx, 25(OH)D and ALP values than others. Taking two kinds of medications did have the correlations to bone turnover and metabolic marker, but we can’t actually judge the causal of the relationship. In bone turnover and metabolic clinical trials, patients who take both alfacalcidol and calcium should be excluded but only take alfacalcidol or calcium didn’t need to be excluded.

Conclusions

P1NP and β-CTx concentrations are highly correlated with age, especially in individuals under 20 years of age. If these two indicators require analysis in a case-control study, patients and controls must be strictly matched by age or hierarchical research should be conducted according to age if subjects are under 20 years of age. In contrast, researchers do not need to match patients strictly if all subjects are older than 20. The threshold for ALP is 40 years of age. The situation for Ca and P is extremely complicated and we strongly recommend strictly matching according to age in disease research. The difference in P1NP, β-CTx, 25(OH)D and ALP concentrations between genders depends on the difference in age. These don’t require strict gender matching if age has been matched. Medication history should be considered in bone turnover and metabolic clinical research.

Availability of data and materials

The datasets used during the current study are available from the corresponding author on reasonable request.

Abbreviations

P1NP:

Amino-terminal propeptide of type I collagen

β-CTx:

Beta-isomerized C-terminal telopeptide of type I collagen

ALP:

Alkaline phosphatase

25(OH)D:

25-hydroxyvitamin D

PTH:

Parathyroid hormone

Ca:

Calcium

P:

Phosphorus

References

  1. Riggs BL, Melton LJ III, Robb RA, Camp JJ, Atkinson EJ, McDaniel L, Amin S, Rouleau PA, Khosla S. A population-based assessment of rates of bone loss at multiple skeletal sites: evidence for substantial trabecular bone loss in young adult women and men. J Bone Miner Res. 2008;23(2):205–14.

    Article  Google Scholar 

  2. Lieuw AFM, Sierra RI, Specker BL. Carboxy-terminal propeptide of human type I collagen and pyridinium cross-links as markers of bone growth in infants 1 to 18 months of age. J Bone Miner Res. 1995;10(6):849–53.

    Article  Google Scholar 

  3. Tobiume H, Kanzaki S, Hida S, Ono T, Moriwake T, Yamauchi S, Tanaka H, Seino Y. Serum bone alkaline phosphatase isoenzyme levels in normal children and children with growth hormone (GH) deficiency: a potential marker for bone formation and response to GH therapy. J Clin Endocrinol Metab. 1997;82(7):2056–61.

    CAS  PubMed  Google Scholar 

  4. Marwaha RK, Garg MK, Mithal A, Gupta S, Shukla M, Chadha A. Effect of vitamin D supplementation on bone turnover markers in children and adolescents from North India. Indian J Endocrinol Metab. 2019;23(1):27–34.

    Article  CAS  Google Scholar 

  5. Jorgensen NR, Mollehave LT, Hansen YBL, Quardon N, Lylloff L, Linneberg A. Comparison of two automated assays of BTM (CTX and P1NP) and reference intervals in a Danish population. Osteoporos Int. 2017;28(7):2103–13.

    Article  CAS  Google Scholar 

  6. Pirila S, Taskinen M, Turanlahti M, Kajosaari M, Makitie O, Saarinen-Pihkala UM, Viljakainen H. Bone health and risk factors of cardiovascular disease--a cross-sectional study in healthy young adults. PLoS One. 2014;9(10):e108040.

    Article  Google Scholar 

  7. Brett NR, Parks CA, Lavery P, Agellon S, Vanstone CA, Kaufmann M, Jones G, Maguire JL, Rauch F, Weiler HA. Vitamin D status and functional health outcomes in children aged 2-8 y: a 6-mo vitamin D randomized controlled trial. Am J Clin Nutr. 2018;107(3):355–64.

    Article  Google Scholar 

  8. Puntus T, Schneider B, Meran J, Peterlik M, Kudlacek S. Influence of age and gender on associations of body mass index with bone mineral density, bone turnover markers and circulating calcium-regulating and bone-active sex hormones. Bone. 2011;49(4):824–9.

    Article  CAS  Google Scholar 

  9. Li M, Li Y, Deng W, Zhang Z, Deng Z, Hu Y, Xia W, Xu L. Chinese bone turnover marker study: reference ranges for C-terminal telopeptide of type I collagen and procollagen I N-terminal peptide by age and gender. PLoS One. 2014;9(8):e103841.

    Article  Google Scholar 

  10. Lu HK, Zhang Z, Ke YH, He JW, Fu WZ, Zhang CQ, Zhang ZL. High prevalence of vitamin D insufficiency in China: relationship with the levels of parathyroid hormone and markers of bone turnover. PLoS One. 2012;7(11):e47264.

    Article  CAS  Google Scholar 

  11. Hu WW, Zhang Z, He JW, Fu WZ, Wang C, Zhang H, Yue H, Gu JM, Zhang ZL. Establishing reference intervals for bone turnover markers in the healthy shanghai population and the relationship with bone mineral density in postmenopausal women. Int J Endocrinol. 2013;2013:513925.

    PubMed  PubMed Central  Google Scholar 

  12. Xie Z, Xia W, Zhang Z, Wu W, Lu C, Tao S, Wu L, Gu J, Chandler J, Peter S, et al. Prevalence of vitamin D inadequacy among Chinese postmenopausal women: a Nationwide, multicenter, cross-sectional study. Front Endocrinol (Lausanne). 2018;9:782.

    Article  Google Scholar 

  13. Gascoin-Lachambre G, Trivin C, Brauner R, Souberbielle JC. Serum procollagen type 1 amino-terminal propeptide (P1NP) as an early predictor of the growth response to growth hormone treatment: comparison of intrauterine growth retardation and idiopathic short stature. Growth Hormon IGF Res. 2007;17(3):194–200.

    Article  CAS  Google Scholar 

  14. van Lierop AH, Hamdy NA, van Egmond ME, Bakker E, Dikkers FG, Papapoulos SE. Van Buchem disease: clinical, biochemical, and densitometric features of patients and disease carriers. J Bone Miner Res. 2013;28(4):848–54.

    Article  Google Scholar 

  15. Neumeyer AM, Gates A, Ferrone C, Lee H, Misra M. Bone density in peripubertal boys with autism spectrum disorders. J Autism Dev Disord. 2013;43(7):1623–9.

    Article  Google Scholar 

  16. Szulc P, Kaufman JM, Delmas PD. Biochemical assessment of bone turnover and bone fragility in men. Osteoporos Int. 2007;18(11):1451–61.

    Article  CAS  Google Scholar 

  17. Sodi R, Hazell MJ, Durham BH, Rees C, Ranganath LR, Fraser WD. The circulating concentration and ratio of total and high molecular weight adiponectin in post-menopausal women with and without osteoporosis and its association with body mass index and biochemical markers of bone metabolism. Clin Biochem. 2009;42(13–14):1375–80.

    Article  CAS  Google Scholar 

  18. Botella S, Restituto P, Monreal I, Colina I, Calleja A, Varo N. Traditional and novel bone remodeling markers in premenopausal and postmenopausal women. The Journal of Clinical Endocrinology & Metabolism. 2013;98(11):E1740–8.

    Article  CAS  Google Scholar 

  19. Martinez J, Olmos JM, Hernandez JL, Pinedo G, Llorca J, Obregon E, Valero C, Gonzalez-Macias J. Bone turnover markers in Spanish postmenopausal women: the Camargo cohort study. Clin Chim Acta. 2009;409(1–2):70–4.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

The authors gratefully acknowledge all participants of the study.

Funding

This study was supported by the National Natural Science Foundation of China (grant no. 81601397) and National Natural Science Foundation of China (grant no. 81771727). The funding body had no role in the design of the study and collection, analysis, interpretation of data or in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

J S, SS Z, Y Q and BB L collected the data. J W and J S conceived the study and analyzed the results. J W and QH Y supervised the study and prepared the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jing Wu.

Ethics declarations

Ethics approval and consent to participate

ZhuJiang Hospital of Southern Medical University Ethics Committee waived the approval due to the retrospective nature of the study.

Consent for publication

Not applicable.

Competing interests

The authors report no conflicts of interest in this work.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

General characteristics of the 36 kinds of diseases.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shao, J., Zhou, SS., Qu, Y. et al. Correlation between bone turnover and metabolic markers with age and gender: a cross-sectional study of hospital information system data. BMC Musculoskelet Disord 21, 603 (2020). https://doi.org/10.1186/s12891-020-03610-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-020-03610-w

Keywords