Biomarkers as therapy monitoring for postmenopausal osteoporosis: a systematic review

1Bone is highly dynamic with resorption and ossification to maintain tissue homeostasis. Bone alkaline phosphatase (bALP) and procollagen type I N propeptide (PINP) have been considered biomarkers of bone ossification, while serum cross-linked C-telopeptides of type I collagen (bCTx) and urinary cross-linked N-telopeptides of type I collagen (NTx) are indicators of bone resorption. Bone turnover markers (BTMs) highlight the dynamic balance of the bone tissue. Markers of ossification (bALP and PINP) derive from the procollagen metabolism or from osteoblasts. Markers of resorption (bCTx and NTx) are produced by osteoclasts or result from collagen degradation processes. BTMs are influenced by several endogenous factors, such as gender, age, ethnicity, fracture, and associated diseases. Exogenous factors, such as circadian rhythm, seasonal variation, diet, and exercise, also influence BTMs. Recently, many studies used BTMs to monitor the efficacy and safety of drugs influencing bone turnover and as therapy monitors in postmenopausal osteoporosis. Although the use of these biomarkers in clinical practice is common, their role as therapy monitors is still unclear. Indeed, no previous studies performed a systematical evaluation of their potential as therapy monitors in postmenopausal osteoporosis.

2The purpose of the study was to explore the potential of bALP, PINP, bCTx, and NTx in therapy monitoring for postmenopausal osteoporosis, investigating their association with bone mineral density (BMD) and the rate of adverse events.

3This study was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses: the PRISMA guidelines. The PICOTD algorithm was preliminarily set out: P (Problem): Postmenopausal osteoporosisI (Intervention): Therapy monitoringC (Comparison): bALP, PINP, bCTx, NTxO (Outcomes): BMD, rate of fractures, and adverse eventsT (Timing): Minimum 6 months of follow-upD (Design): RCTs

4P (Problem): Postmenopausal osteoporosis

5I (Intervention): Therapy monitoring

6C (Comparison): bALP, PINP, bCTx, NTx

7O (Outcomes): BMD, rate of fractures, and adverse events

8T (Timing): Minimum 6 months of follow-up

9D (Design): RCTs

10Two independent authors (**;**) performed the literature search in April 2021. The following databases were accessed: PubMed, Google Scholar, EMBASE, and Scopus. No time constrains were used for the search. The following keywords were used in combination: osteoporosis, treatment, management, drug, pharmacology, pharmacological, medicament, mineral, density, bone, BMD, bone alkaline phosphatase, ALP, procollagen type I N propeptide, PINP, serum cross-linked C-telopeptides of type I collagen, CTx, urinary cross-linked N-telopeptides of type I collagen, NTx, premenopausal, spine, pathological, fragility, fractures, hip, vertebral, disability, adverse events, Calcium, Vitamin D, PTH, osteoblast, and osteoclast. The same authors independently performed the initial screening. If the title and abstract matched the topic, the article full-text was accessed. A cross reference of the bibliographies was also performed to identify further studies.

11All randomized clinical trials (RCTs) comparing two or more pharmacological treatments for postmenopausal osteoporosis were accessed. According to the authors’ language capabilities, articles in English, French, German, Italian, Portuguese, and Spanish were eligible. Only level I studies, according to Oxford Centre of Evidence-Based Medicine, were considered for inclusion. Only articles reporting quantitative data under the outcomes of interest were eligible. Only clinical studies that reported the amount of bALP, PINP, bCTx, and NTx at last follow-up were included. Articles including patients with secondary osteoporosis were excluded. Studies concerning patients with tumors and/or bone metastases were also not included. Studies reporting data on patients with iatrogenic-induced menopausal and those on pediatric and/ or adolescent patients were not included. Combined therapies with multiple drugs were also not considered in the present study. Studies regarding selected patients undergoing immunosuppressive therapies or organ transplantation were also not considered. Studies with follow-up shorter than 6 months were not eligible, nor where those involving less than 10 patients. Studies reporting data of combined therapy with multiple anti-osteoporotic drugs were also not included. Missing data under these endpoints warranted the exclusion from the present work.

12Two authors (**;**) independently performed data extraction. Study generalities (author, year, journal, duration of the follow-up) and patient baseline demographic information were collected: number of samples and related mean age, mean body mass index (BMI), and mean bone mass index (BMD) of the spine, hip, and femur neck. Data concerning the following endpoints were collected at last follow-up: rate of vertebral, femoral, and hip osteoporotic fractures. Further, data concerning the following complications were collected: serious adverse events and those leading to study discontinuation, gastrointestinal events, musculoskeletal events, and mortality. Data concerning bALP, PINP, bCTx, and NTx were extracted at last follow-up. The ultimate aim was to assess association between biomarkers and clinical outcomes at last follow-up in terms of BMD, rate of pathological fractures, and adverse events.

13The methodological quality assessment was made through the risk of bias summary tool of the Review Manager Software (The Nordic Cochrane Collaboration, Copenhagen). The following risks of bias were evaluated: selection, detection, performance, reporting, attrition, and other sources of bias.

14The statistical analyses were performed by the main author (**). The IBM SPSS software version 25 was used to assess baseline data. The Shapiro-Wilk test was performed to investigate data distribution. For normal data, mean and standard deviation (SD) were calculated. For non-parametric data, median and interquartile range (IQR) were calculated. The Student T-test was used to assess significance for parametric data, while the Mann-Whitney U-test was used for non-parametric variables. Values of P < 0.05 considered statistically significant. A multivariate analysis was performed to assess associations between biomarkers and clinical outcomes at last follow-up. The STATA Software/MP (StataCorporation, College Station, TX, USA) was used for the statistical analyses. A multiple linear model regression analysis through the Pearson product-moment correlation coefficient (r) was used. The Cauchy–Schwarz formula was used for inequality: +1 was considered as positive linear correlation, and −1 a negative one. Values of 0.1< | r | < 0.3, 0.3< | r | < 0.5, and | r | > 0.5 were considered to have weak, moderate, and strong correlation, respectively. The overall significance was performed through the χ2 test, with values of P < 0.05 considered statistically significant.

15The literature search resulted in 1174 studies. Of them, 307 were duplicates. A further 749 articles were excluded because of nature of the study (N = 233), non-clinical studies (N = 301), secondary osteoporosis (N = 81), small population or short follow-up (N = 19), multiple therapies (N = 21), language limitations (N = 9), uncertain results (N = 13), and others (N = 72). Another 102 articles were excluded because data under the outcomes of interest were missing. Finally, 16 RCTs were eligible for inclusion in the present study (Fig. 1).

16The inclusion of only RCTs yields to low risk of selection bias. Many studies performed patients and assessor blinding, thus leading to moderate-low risk of detection and performance biases. The overall high quality of the included studies led also to an overall low-risk of attrition and reporting bias. Overall, the results of the evaluation of each risk of bias item for each individual study included in the present analysis was low to moderate, leading to a good assessment of the methodology. The risk of bias graph is shown in Fig. 2.

17A total of 14,446 patients were included. The median age was 67 (IQR 4.7), and the median BMI 25.4 (IQR 1.9). The median vertebral BMD was 0.82 (IQR 0.14), hip BMD 0.79 (IQR 0.1), and femur BMD 0.64 (IQR 0.02). The ANOVA test found optimal within-group variance concerning mean age, BMI, and BMD (P > 0.1). Generalities and patient baseline data of the included studies are shown in detail in Table 1.

18Greater bALP was associated with lower femoral BMD (r = − 0.87; P = 0.01). Greater NTx was associated with greater of occurrence of non-vertebral fractures (r = 0.98; P = 0.02). Greater NTx was associated with greater rate of therapy discontinuation (r = − 0.60; P = 0.04). There was evidence of positive association between PINP and CTx (r = −0.93; P = 0.0001). No other statistically significant associations were detected. These results are shown in Table 2. Added-variable plots of the statistically significant outcomes are displayed in Fig. 3.

19Our findings suggest that bALP and NTx may represent useful, valid, and reliable tools for therapy monitoring for postmenopausal osteoporotic patients. Higher bALP and NTx were associated to lower femoral BMD and higher rate of non-vertebral fractures, respectively. Furthermore, a positive association between NTx and the rate of adverse events leading to therapy discontinuation was evidenced. BTMs are implicated in bone turnover, and their level significantly varies during osteoporotic therapy. P1NP and bCTx did not show any statistically significant association with any of the considered variables in the present study. PINP is released following to the amino/carboxy-terminal extensions cleavage of the procollagen and can be found variably in the blood. bCTx is a form of the telopeptides of type I collagen released during collagen degradation. We were unable to find any significant association for these two BTMs; thus, their potential in therapy monitoring remains uncertain.

20Given their sensibility to reveal changes in bone turnover, BTMs gained popularity. BTMs’ variations related also to antiresorptive drugs, which produce a quick decrease of the bone resorption markers, followed by those of bone formation. Vice versa, anabolic drugs increase the level of bone formation markers, followed by those of bone resorption. BTMs’ changes are related to the risk of fragility fractures; thus, BTMs have been introduced to monitor therapy in osteoporotic patients. The effect of the therapy on BTMs strictly depends on the type of drug used. Antiresorptive drugs inhibit osteoclasts and cause a rapid reduction in resorption markers, followed by a reduction in bone formation markers. Indeed, the present study shows an association between BMD and bALP. BALP is a membrane-bound enzyme found in almost all tissues of the organism and can be easily measured in serum. bALP has been the first BTMs of bone turnover intensively investigated. It was initially used to monitor the efficacy and safety of some drugs acting on bone turnover, and subsequently acquired popularity to monitor therapy in osteoporotic patients. Bjarnason et al. found that bALP had stronger association than BMD to predict the risk of fragility fractures in patients undergoing raloxifene therapy. These results were unexpected, since BMD was considered a very reliable measure of the risk of fractures. Similar results were found by Iwamoto et al. and Gonnelli et al. evaluating the outcome of alendronate zoledronate and ibandronate on patients’ quality of life. Both studies evidenced a statistically significant association between the increase in BMD and a decrease in serum bALP levels. Comparable results were obtained by Delmas et al., evaluating the efficacy of raloxifene in preventing vertebral fracture in patients with postmenopausal osteoporosis. Muratore et al. found that the bALP and BMD variations were proportional to the dose of clodronate administered to patients. Overall, these findings encouraged the use of bALP to monitor therapy in patients undergoing pharmacological management of postmenopausal osteoporosis.

21Our analyses showed evidence of positive association between NTx and the rate of non-vertebral fractures. NTx and bCTX are two different forms of the telopeptide of type I collagen, which modulate the degradation process of collagen. These telopeptides are measurable in serum and in the urines and exhibit to circadian cycle variations. Iwamoto et al. demonstrated that alendronate reduced the urinary level of NTx. Garnero et al. found that the urinary excretion of NTx did not predict fractures, hypothesizing that it follows a different pattern of bone collagen degradation.

22This study shows limitations. The analyses were performed regardless to the drug type and administration. This enhanced the risk of bias of the present study, negatively affecting the reliability of our results. Furthermore, the heterogeneous daily administration of vitamin D and calcium represents another important limitation. We included only RCTs reporting quantitative data under the outcomes of interest, which were published in peer reviewed scientific. However, the role of BTMs has been poorly investigated, and none of the included articles did not aim to quantify directly the biomarkers’ variations. Results from this study should encourage future investigation to evaluate the potential of BTMs in a clinical setting, analyzing their variations as primary outcome. The biological variability of BTMs constitutes an important factor limiting their engagement in the management of osteoporosis.

23The analysis of BTMs in the investigation of their possible role in monitoring therapy demonstrates the need for studies that can validate their use in clinical practice. Our analysis supports the adoption of BTMs in therapy monitoring of postmenopausal osteoporosis patients. Further studies are needed to analyze variations of BTMs in relation to treatment as a primary outcome.