Next Article in Journal
Difference in Mortality Rates by Occupation in Japanese Male Workers Aged 25 to 64 Years from 1980 to 2015
Next Article in Special Issue
Discrete Event Simulation Model for Cost-Effectiveness Evaluation of Screening for Asymptomatic Patients with Lower Extremity Arterial Disease
Previous Article in Journal
Differences between Healthy-Weight and Overweight Serbian Preschool Children in Motor and Cognitive Abilities
Previous Article in Special Issue
Endothelial Function Assessment by Flow-Mediated Dilation Method: A Valuable Tool in the Evaluation of the Cardiovascular System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Association between Osteoporosis and Peripheral Artery Disease: A Population-Based Longitudinal Follow-Up Study in Taiwan

1
Department of Orthopedics, Fu Jen Catholic University Hospital, Fu Jen Catholic University, No. 69, Guizi Rd., Taishan Dist., New Taipei City 24352, Taiwan
2
Department of Family Medicine, Fu Jen Catholic University Hospital, Fu Jen Catholic University, No. 69, Guizi Rd., Taishan Dist., New Taipei City 24352, Taiwan
3
Graduate Institute of Business Administration, College of Management, Fu Jen Catholic University, No. 510, Zhongzheng Rd., Xinzhuang Dist., New Taipei City 242062, Taiwan
4
Artificial Intelligence Development Center, Fu Jen Catholic University, No. 510, Zhongzheng Rd., Xin-Zhuang Dist., New Taipei City 242062, Taiwan
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(18), 11327; https://doi.org/10.3390/ijerph191811327
Submission received: 18 July 2022 / Revised: 5 September 2022 / Accepted: 7 September 2022 / Published: 9 September 2022
(This article belongs to the Special Issue Vascular Disease and Health)

Abstract

:
Background: There are several possible links that have been used to claim that osteoporosis and peripheral artery disease (PAD) are associated; however, the solid evidence is not sufficient. This study aimed to use the Taiwan National Health Insurance Research Database (NHIRD) to determine if osteoporosis is associated with peripheral artery disease (PAD). Method: NHIRD records from 23 million patients were collected to recruit two matched cohort groups: 64,562 patients with and 64,562 patients without osteoporosis. To compare the crude hazard ratio (HR) and the incidence rate ratio between the two cohort groups for PAD, the Cox model was used. Result: With an adjusted HR of 1.18 (95% CI, 1.08–1.29), the osteoporosis cohort group had a significantly greater risk for PAD than the group without osteoporosis. The cumulative incidence of PAD in the cohort group was also statistically higher than it was in the group without osteoporosis (1.71% and 1.39%; p ≤ 0.0001, log-rank) over the 10-year follow-up period. In addition, the osteoporotic patients with ischemic stroke, chronic obstructive pulmonary disease (COPD), and congestive heart failure (CHF) had a significantly increased risk of PAD based on subgroup analysis. Conclusions: There was a positive association between osteoporosis and the development of PAD, as patients with osteoporosis had an increased incidence of PAD over time.

1. Introduction

Peripheral artery disease (PAD) has frequently been investigated to determine its prevalence, pathophysiology, frequency and occurrence, management, and prevention in previous literature. Despite a decline in mortality rates, PAD remains highly prevalent due to the aging population, with a prevalence of approximately 10 to 20% [1]. In addition, the economic burden of atherothrombotic disease is high and is expected to increase with the aging population, demonstrating an increased survival rate. Moreover, the direct costs associated with PAD are even higher than those for coronary artery disease because of the higher rates of polyvascular disease and the higher number of annual cardiovascular (CV) events and hospitalization rates [2]. On the other hand, osteoporosis is a chronic degenerative disease that has incidence rates that increase with age. With its association with age-related fragile fractures, osteoporosis may progress without symptoms, and it is now recognized as a major threat to the elderly population. Studies are becoming increasingly focused on the links between atherothrombotic disease and osteoporosis from many aspects, including epidemiology and pathophysiology [3,4,5]. Between the mineralization of bone tissue and the calcification of the blood vessel wall, there is abundant evidence that claims their association [6]. However, while osteoporosis represents a proven link between coronary heart disease and ischemic stroke [7,8], studies investigating the relationship between osteoporosis and PAD remain controversial. There is one small cohort study that discusses the association between PAD and osteoporosis, and the results found a weak and age-dependent association between PAD and osteoporosis that was only present in women [9]. A cross-sectional study in the USA reported that the prevalence of PAD was significantly higher in women with osteopenia (4.8%) and osteoporosis (11.8%) compared to in women with a normal body mineral density (BMD) (3.3%) (p < 0.001). Osteopenia (OR: 1.15) and osteoporosis (OR: 1.8) were found to be independent risk factors for the presence of PAD in women [10]. On the other hand, one study conducted in the Netherlands showed that women with a low femoral neck BMD had a significantly increased risk of PAD (OR: 1.49), but this was not the case in men [11]. Although both diseases have improved our understanding of the pathophysiological factors and associations shared by advanced science, the solid evidence-based studies investigating the etiologic role of osteoporosis associated with PAD are insufficient. Therefore, the purpose of this study was to observe and analyze NHIRD records to determine if osteoporosis is associated with PAD. The hypothesis for this study was that osteoporosis is correlated with an increased incidence of PAD.

2. Materials and Methods

2.1. Data Sources

This study used the National Health Insurance Research Database (NHIRD), which is a large, computerized database that represents more than 99% of the Taiwanese population (approximately 23 million individuals) [12]. The database diagnosis codes are based on the International Classification of Diseases 9th Revision Clinical Modification (ICD-9-CM), with the 10th Revision (ICD-10-CM) being implemented in 2016. The NHIRD has often been used in previous studies identifying and using data on patients, as it includes abundant information on basic demographics, disease diagnoses, medical expenses, prescriptions, surgical operations, and medication services [13]. This study was approved by the Fu Jen Catholic University Institutional Review Board (IRB) and Ethics Committee under approval number No. C108094. Due to the NHIRD’s ethical guidelines, the personal information of the patients was anonymized before accessing the database, and thus, the IRB waived the requirements for informed consent.

2.2. Study Population

This research was a retrospective, case-compared, database cohort study for the evaluation of the different risks of developing PAD in an osteoporosis group and in a nonosteoporosis group. If patients were (1) older than 50 years of age and (2) had been diagnosed with osteoporosis (ICD-9-CM Code: 733.0) by a physician from 1 January 2006 to 31 December 2008, they were included in this study. If a patient had at least three different medical claims issued in an out-patient setting or at least one claim in an in-patient setting, they were defined as having osteoporosis. BMD measurement is the gold standard for osteoporosis diagnosis in the spine and hip and measures how different the mineral and calcium levels are in a particular area of the bone. A value of 2.5 or more in the standard deviation (SD) indicates osteoporosis. As shown in Figure 1, a total of 125,762 osteoporosis cases were identified after excluding 44,229 patients that were younger than 50 years old; 721 patients that had missing data; and 103,805 patients that had PAD, osteoporosis, or osteoporotic fractures, including vertebral compression fractures and hip fractures, from 2002 to 2007. A total of 2,541,021 patients without osteoporosis were included after the same exclusion criteria were applied. The nonosteoporosis cohort was verified as maintaining osteoporosis-free status throughout the entire follow-up period to minimize any allocation bias. The index date was determined from the first osteoporosis date. In order to avoid reverse causation bias, patients who were diagnosed with PAD within one year after the index date were excluded from the study. A 1:1 ratio was also used for propensity score matching along with the index year, age, and sex to reduce sampling bias. The final osteoporosis cohort contained 64,562 patients, with the comparison cohort containing the same number of patients after propensity score matching. Both cohort groups were followed up until the occurrence of PAD, death, or December 2018, whichever occurred first.

2.3. Outcome Measure and Confounding Variables

The incidence of PAD was the primary outcome measurement. At least three different medical statements issued in an out-patient setting or at least one claim in an in-patient setting was the definition for diagnosing PAD (ICD-9-CM codes: 440.0, 440.2, 440.3, 440.8, 440.9, 443.9, 444.0, 444.2, 444.8, 447.8, and/or 447.9). During the first year of the follow up, patients with PAD were excluded from the overall risk calculation to avoid protopathic bias. Table 1 shows that the two cohorts were balanced for index date, age, and sex after matching. Comorbidity data were collected for hypertension, hyperlipidemia, chronic obstructive pulmonary disease (COPD), diabetes mellitus (DM), liver disease, peripheral vascular disease (PVD), congestive heart failure (CHF), stroke, rheumatoid arthritis (RA), chronic kidney disease (CKD), gout, and ischemic stroke. In addition, different kinds of anti-osteoporosis drugs (including selective estrogen-receptor modulator, denosumab, bisphosphonates, and teriparatide) were added as confounding variables.

2.4. Statistical Analysis

The basic characteristics are presented as N (%) for the categorical variables according to the chi-square test, and the continuous variables are defined as the mean ± standard deviation (mean ± SD) according to a t-test. For the hazard ratios (aHRs) with 95% confidence intervals (CIs) and for the estimated incidences of the case and comparison groups, multivariate-adjusted Cox proportional hazard regression was calculated. A forest plot was applied to demonstrate the differences between the two cohorts in the subgroup analysis. Kaplan–Meier curve analyses were used to evaluate the cumulative incidences of patients with PAD, and the log-rank test was applied to determine the significance. All statistical analyses used the two-tailed p < 0.05 and were conducted with SAS version 9.4 (SAS Institute, Cary, NC, USA). The method used for feature selection was implemented using R software (version 3.4.3; R Foundation for Statistical Computing, Vienna, Austria).

3. Results

Table 1 displays the baseline characteristics and comorbidities of the study participants. There were 64,562 cases in the study and control groups. Their sex and age distributions were identical. Several comorbidities with significant differences between cohorts were found for hypertension, hyperlipidemia, chronic obstructive pulmonary diseases (COPDs), diabetes mellitus (DM), liver disease, congestive heart failure (CHF), rheumatoid arthritis (RA), chronic kidney disease (CKD), gout, and ischemic stroke. Regarding anti-osteoporosis drug usage, the osteoporosis group only had a high incidence of denosumab usage over their nonosteoporosis counterparts. On the other hand, the incidence of prescriptions of raloxifene, oral bisphosphate, zoledronic acid, and teriparatide were higher in the nonosteoporosis group than in the study group. There was no difference in ibandronate usage between both groups (p = 0.06).
As shown in Table 2, the osteoporosis cohort group showed a significant increase in the risk for PAD, with a crude hazard ratio (HR) and adjusted HR of 1.21 (95% confidence level (CI), 1.10–1.32) and 1.18 (95% CI, 1.08–1.29), respectively. Moreover, 1105 osteoporosis patients and 898 comparison-matched patients received diagnoses of PAD by a physician after the 10-year follow up. The incidence rate ratio was 1.19 (incidence rates of PAD per 100,000 patients per year were 172.9 in the osteoporosis cohort and 145.5 in the comparison cohort group). The crude HR was also statistically significant, showing a 21% increase in the risk of PAD in the osteoporosis group, with an HR of 1.21 (95% CI, 21.10–1.32). This significant increase in risk was sustained, even after the application of the multivariate Cox model adjusting for the aforementioned confounding factors, with an adjusted HR of 1.18 (95% CI, 1.08–1.29, p < 0.001).
Table 3 shows the adjusted HRs of PAD stratified by each patient characteristic in both groups using the multivariate Cox model analysis of the differential risk among the different groups. The results reveal that men and women have different levels of PAD risk, with a significant decrease in the risk for female patients (aHR = 0.658; 95% CI, 0.59–0.73; p < 0.001). For patients who were older in age (i.e., age 55), DM, RA, and CKD showed an increased risk of PAD compared to those who did not have these comorbidities. Regarding anti-osteoporosis drug usage, there was no significant difference in the risk of PAD. Figure 2 demonstrates the results of our multivariate Cox model analyzing the differential risk among different categorizations of PAD development stratified by the patient characteristics in the osteoporosis cohort. Osteoporosis patients with ischemic stroke had an increased risk of PAD in our study (aHR = 1.94; 95% CI, 1.40–2.67), as shown in Figure 2. In addition, the osteoporotic patients with COPD and CHF had a significantly increased risk of PAD. Our study also revealed an interesting finding that osteoporosis patients with oral bisphosphonate use had a decreased risk of PAD (aHR = 0.65; 95% CI, 0.43–0.97). On the contrary, an increased risk of PAD was noted in the study group that used ibandronate (aHR = 2.61; 95% CI, 1.08–6.31). A multivariate Cox model that excluded the follow-up events from the first year was constructed to derive the cumulative incidence function. In Figure 3, the cumulative incidence of PAD was 1.71% and 1.39% (p < 0.001, log-rank) for the osteoporosis and nonosteoporosis cohorts, respectively.

4. Discussion

Osteoporosis and atherosclerosis are two chronic degenerative diseases that share several biochemical pathways and risk factors. A previous study proved an association between osteoporosis and coronary artery disease, carotid atherosclerosis, and ischemic stroke [14,15,16]. However, studies investigating an association between PAD and osteoporosis are few and conflicting. There was a significant increase in the risk of developing PAD in patients with osteoporosis compared to their counterparts based on the multivariate Cox model (adjusted HR: 1.18). Additionally, based on the multivariate Cox model, the osteoporosis patients with ischemic stroke, COPD, and CHF had an increased risk of PAD in our study. In addition, the osteoporosis patients with oral bisphosphonate use had a decreased risk of PAD (aHR = 0.65; 95% CI, 0.43–0.97). On the contrary, an increased risk of PAD was noted in the study group that used ibandronate (aHR = 2.61; 95% CI, 1.08–6.31). Overall, a positive association between osteoporosis and subsequent PAD development was shown in this study.
The Rancho Bernardo Study [9] found a weak and age-dependent association between PAD and osteoporosis in women but not in men. Another study that had a larger sample size revealed that the prevalence of PAD was significantly higher in patients with osteopenia and osteoporosis compared to in patients with normal BMD in both genders [10]. In a large prospective cohort of 3998 Chinese men and women (65–92 years of age) in Hong Kong, the ABI correlated positively with hip BMD (correlation coefficient = 0.27; p < 0.001). However, after adjusting for confounding factors, the correlation was much weaker (correlation coefficient = 0.03; p < 0.05) [17]. This indicates that other factors may contribute to the association between peripheral vascular disease and osteoporosis. In our population-based longitudinal study, a significant increase in the risk of developing PAD was found in the osteoporosis cohort, even after adjusting for comorbidities and anti-osteoporosis drug usage. The results provide evidence of an association between the two diseases. Another current Taiwanese population-based study also used the database to compare osteoporosis and PAD [18], and they also concluded that there was a long-term risk of PAD in people with osteoporosis. There were consistent results between their study and ours, including those indicating that the PAD-free survival rate was significantly lower in the case group than in the control group. Three main differences set our study and the other study apart: First, their study used the Longitudinal Health Insurance Database 2005 (LHID2005). All of the registration and claims data for these 1,000,000 individuals were randomly sampled from the National Health Insurance Program and constitute the LHID2005. Our study used the whole population included in the NHIRD without random sampling. Second, in our study, patients who were diagnosed with PAD within one year after the index date were excluded to avoid reverse causation bias. Lastly, we included anti-osteoporosis drug usage to evaluate the drug effects.
As for the interconnections between osteoporosis and PAD, many possible pathogenesis pathways have been discussed [19]. First, when analyzing the osteoprotegerin (OPG), RANK, and the RANK-Ligand (RANKL) system, the OPG levels were significantly higher in patients with PAD than they were in patients without PAD [20]. In addition, one recent review article found that OPG levels were correlated with the presence, severity, and progression of PAD in eight articles, whereas in one article, OPG levels were not significantly elevated [21]. Second, many studies have proven that vitamin D deficiency is correlated with osteoporosis and PAD [22,23,24,25]. Third, previous studies have shown that inflammatory factors have aggravated atherogenesis and bone resorption [15]. High serum levels of C-reactive protein (CRP), interleukin (IL-6), and tumor necrosis factor-alpha (TNF-α) have all been associated with a higher risk developing PAD [26,27]. According to pathophysiologic studies, there is an association between osteoporosis and PAD. This is consistent with our study, which also provided an indication of an association. Nevertheless, a more well-designed method is needed to establish stronger evidenced-based data for a pathophysiologic association to clarify the detailed mechanisms and association.
This study has many advantages. A large sample size, a sufficient follow-up duration of up to ten years, and proper, well-balanced random matching between two cohort groups make this study valuable. However, this study has three limitations. First, due to the many physicians, specialists, and medical staff inputting data into the administrative database, possible coding errors and other coding problems may have occurred. The accuracy and coding were not verified, even though the definition of PAD and osteoporosis was based on the coding definition. However, to the best of our ability, we added two restrictions: the diagnosis (1) must have been made during two ambulatory care visits or during one or more in-patient visits and (2) needed to have been made by certain specialists. These case definitions for chronic health conditions have been widely used in the literature based on administrative data [28]. Second, NHIRD data were limited or unavailable for the proinflammatory cytokines IL-1, IL-6, and TNF-α; serial BMD value measurements; serum 25-hydroxyvitamin D; and serum parathyroid hormones; therefore, these factors could not be analyzed further. Additionally, other data, such as the patients’ BMI and lifestyle, individual dietary habits and physical activity, and data related to PAD and osteoporosis, were not recorded in the database; thus, these factors were also not analyzed. Third, although we added anti-osteoporosis drugs for analysis, different kinds of medications, such as aspirin, statin, anti-DM, and anticoagulants, were not included in the study. To eliminate the confounding factors, matching criteria may need to be added to these medications or to other interventions for propensity score calculation. A well thought-out in vivo or in vitro study should definitely be considered to analyze the pathophysiological link and causation between PAD and OP in the future since this study confirmed a positive association between the two diseases.

5. Conclusions

Patients with osteoporosis had a significant increased risk of developing PAD compared to patients without osteoporosis. A positive association between osteoporosis and the subsequent development of PAD could be observed. Patients with osteoporosis are recommended to receive a personalized risk screening for PAD to prevent the complications of PAD.

Author Contributions

Conceptualization, D.-K.S.; methodology, D.-K.S., S.-H.H. and P.-C.Y.; software, Y.-F.K. and Y.-C.H.; validation, D.-K.S., S.-H.H., P.-C.Y., T.-L.L., Y.-F.K., Y.-C.H. and C.-C.J.; formal analysis, D.-K.S., S.-H.H., P.-C.Y., T.-L.L., Y.-F.K., Y.-C.H. and C.-C.J.; investigation, D.-K.S., S.-H.H. and P.-C.Y.; resources, T.-L.L., Y.-F.K., Y.-C.H., C.-C.J. and M.C.; data curation, Y.-F.K. and Y.-C.H.; writing—original draft preparation, D.-K.S., S.-H.H. and P.-C.Y.; writing—review and editing, D.-K.S., S.-H.H., P.-C.Y., T.-L.L., Y.-F.K., Y.-C.H., C.-C.J. and M.C.; visualization, Y.-C.H.; supervision, C.-C.J. and M.C.; project administration, C.-C.J. and M.C.; funding acquisition, none. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fu Jen Catholic University (A0110152) and the Ministry of Science and Technology (MOST109-2221-E030-011-MY3).

Institutional Review Board Statement

This study was approved by the Institutional Review Board and Ethics Committee of Fu Jen Catholic University (IRB approval number: No C108094). No consent was required for this study. All data collected in this study were encrypted using the same encryption algorithm to cross-link the data while protecting the privacy of the patients.

Informed Consent Statement

Not applicable.

Data Availability Statement

These data were available to us as staff of the Department of Orthopedics at Fu Jen Catholic University Hospital and of Fu Jen Catholic University using the National Hospital Research Database (NHIRD). These data are protected by the Ministry of Health and Welfare and patient privacy laws in Taiwan; no public links to these protected health information datasets are available. These data will be made available to others after receiving appropriate data privacy and human subject approval from the institution. Requests should be sent to [email protected].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bauersachs, R.; Zeymer, U.; Briere, J.-B.; Marre, C.; Bowrin, K.; Huelsebeck, M. Burden of Coronary Artery Disease and Peripheral Artery Disease: A Literature Review. Cardiovasc. Ther. 2019, 2019, 8295054. [Google Scholar] [CrossRef] [PubMed]
  2. Mahoney, E.M.; Wang, K.; Keo, H.H.; Duval, S.; Smolderen, K.G.; Cohen, D.J.; Steg, G.; Bhatt, D.L.; Hirsch, A.T. Vascular hospitalization rates and costs in patients with peripheral artery disease in the United States. Cardiovasc. Qual. Outcomes 2010, 3, 642–651. [Google Scholar] [CrossRef] [PubMed]
  3. Eastell, R.; Newman, C.; Crossman, D.C. Cardiovascular disease and bone. Arch. Biochem. Biophys 2010, 503, 78–83. [Google Scholar] [CrossRef] [PubMed]
  4. Cassar, A.; Poldermans, D.; Rihal, C.S.; Gersh, B.J. The management of combined coronary artery disease and peripheral vascular disease. Eur. Heart J. 2010, 31, 1565–1572. [Google Scholar] [CrossRef] [PubMed]
  5. Narula, N.; Olin, J.W.; Narula, N. Pathologic Disparities Between Peripheral Artery Disease and Coronary Artery Disease. Arterioscler. Thromb. Vasc. Biol. 2020, 40, 1982–1989. [Google Scholar] [CrossRef]
  6. Anagnostis, P.; Karagiannis, A.; Kakafika, A.I.; Tziomalos, K.; Athyros, V.G.; Mikhailidis, D.P. Atherosclerosis and osteoporosis: Age-dependent degenerative processes or related entities? Osteoporos Int. 2009, 20, 197–207. [Google Scholar] [CrossRef]
  7. Hu, X.; Ma, S.; Yang, C.; Wang, W.; Chen, L. Relationship between senile osteoporosis and cardiovascular and cerebrovascular diseases. Exp. Ther. Med. 2019, 17, 4417–4420. [Google Scholar] [CrossRef] [PubMed]
  8. Liu, N.; Chen, J.; Zhang, K.; Tang, Z. A community-based study of the relationship between coronary artery disease and osteoporosis in Chinese postmenopausal women. Coron. Artery Dis. 2016, 27, 59–64. [Google Scholar] [CrossRef]
  9. von Mühlen, D.; Allison, M.; Jassal, S.K.; Barrett-Connor, E. Peripheral arterial disease and osteoporosis in older adults: The Rancho Bernardo Study. Osteoporos. Int. 2009, 20, 2071–2078. [Google Scholar] [CrossRef]
  10. Baldwin, M.J.; Policha, A.; Maldonado, T.; Hiramoto, J.S.; Honig, S.; Conte, M.S.; Berger, J.; Rockman, C.B. Novel association between bone mineral density scores and the prevalence of peripheral artery disease in both sexes. Vasc. Med. 2017, 22, 13–20. [Google Scholar] [CrossRef]
  11. van der Klift, M.; Pols, H.A.P.; Hak, A.E.; Witteman, J.C.M.; Hofman, A.; de Laet, C.E.D.H. Bone mineral density and the risk of peripheral arterial disease: The Rotterdam Study. Calcif. Tissue Int. 2002, 70, 443–449. [Google Scholar] [CrossRef]
  12. Lin, A.C.-C.; Chao, E.; Yang, C.-M.; Wen, H.-C.; Ma, H.-L.; Lu, T.-C. Costs of staged versus simultaneous bilateral total knee arthroplasty: A population-based study of the Taiwanese National Health Insurance Database. J. Orthop. Surg. Res. 2014, 9, 59. [Google Scholar] [CrossRef]
  13. Hsiao, F.-Y.; Huang, Y.-T.; Yang, C.-L.; Huang, W.-F. Using Taiwan’s National Health Insurance Research Databases for Pharmacoepidemiology Research. J. Food Drug Anal. 2007, 15, 99–108. [Google Scholar] [CrossRef]
  14. Laroche, M.; Pécourneau, V.; Blain, H.; Breuil, V.; Chapurlat, R.; Cortet, B.; Sutter, B.; Degboé, Y. Osteoporosis and ischemic cardiovascular disease. Jt. Bone Spine 2017, 84, 427–432. [Google Scholar] [CrossRef]
  15. Zhang, Y.; He, B.; Wang, H.; Shi, J.; Liang, H. Associations between bone mineral density and coronary artery disease: A meta-analysis of cross-sectional studies. Arch. Osteoporos. 2020, 15, 24. [Google Scholar] [CrossRef]
  16. Lin, C.-H.; Chang, W.-C.; Kuo, C.-N.; Yu, H.-C.; Yang, C.-C.; Lin, Y.-W.; Hung, K.-S. A population-based five-year study on the risk of stroke in patients with osteoporosis in Taiwan. Bone 2015, 72, 9–13. [Google Scholar] [CrossRef]
  17. Wong, S.Y.-S.; Kwok, T.; Woo, J.; Lynn, H.; Griffith, J.F.; Leung, J.; Tang, Y.Y.N.; Leung, P.C. Bone mineral density and the risk of peripheral arterial disease in men and women: Results from Mr. and Ms Os, Hong Kong. Osteoporos. Int. 2005, 16, 1933–1938. [Google Scholar] [CrossRef]
  18. Liang, Y.S.; Yeh, K.C.; Pan, S.L. Osteoporosis and the long-term risk of peripheral artery disease: A population-based longitudinal follow-up study in Taiwan. Osteoporos. Int. 2022, 33, 1117–1123. [Google Scholar] [CrossRef]
  19. Gaudio, A.; Xourafa, A.; Rapisarda, R.; Castellino, P.; Signorelli, S.S. Peripheral artery disease and osteoporosis: Not only age-related (Review). Mol. Med. Rep. 2018, 18, 4787–4792. [Google Scholar] [CrossRef]
  20. Demková, K.; Kozárová, M.; Malachovská, Z.; Javorský, M.; Tkáč, I. Osteoprotegerin concentration is associated with the presence and severity of peripheral arterial disease in type 2 diabetes mellitus. Vasa 2018, 47, 131–135. [Google Scholar] [CrossRef]
  21. Kapetanios, D.; Karkos, C.; Giagtzidis, I.; Papazoglou, K.; Kiroplastis, K.; Spyridis, C. Vascular calcification biomarkers and peripheral arterial disease. Int. Angiol. 2016, 35, 455–459. [Google Scholar] [PubMed]
  22. Holick, M.F. Vitamin D deficiency. N. Engl. J. Med. 2007, 357, 266–281. [Google Scholar] [CrossRef]
  23. Narula, R.; Tauseef, M.; Ahmad, I.A.; Agarwal, K.; Ashok, A.; Anjana, A. Vitamin d deficiency among postmenopausal women with osteoporosis. J. Clin. Diagn. Res. 2003, 7, 336–338. [Google Scholar] [CrossRef] [PubMed]
  24. Oh, J.; Weng, S.; Felton, S.K.; Bhandare, S.; Riek, A.; Butler, B.; Proctor, B.M.; Petty, M.; Chen, Z.; Schechtman, K.B.; et al. 1,25(OH) 2 Vitamin D Inhibits Foam Cell Formation and Suppresses Macrophage Cholesterol Uptake in Patients With Type 2 Diabetes Mellitus. Circulation 2009, 120, 687–698. [Google Scholar] [CrossRef]
  25. Sharma, J.; Raggi, P.; Melamed, K.L.; Muntner, P.; Michos, E.D.; Uribarri, J.; Weber, C. Serum 25-hydroxyvitamin D levels and the prevalence of peripheral arterial disease: Results from NHANES 2001 to 2004. Arterioscler Thromb. Vasc. Biol. 2008, 28, 1179–1185. [Google Scholar]
  26. Tzoulaki, I.; Murray, G.D.; Tzoulaki, I.; Murray, G.D.; Lee, A.J.; Rumley, A.; Lowe, G.D.; Fowkes, F.G.R. C-reactive protein, interleukin-6, and soluble adhesion molecules as predictors of progressive peripheral atherosclerosis in the general population: Edinburgh Artery Study. Circulation 2005, 112, 976–983. [Google Scholar] [CrossRef]
  27. Pande, R.L.; Brown, J.; Buck, S.; Redline, W.; Doyle, J.; Plutzky, J.; Creager, M.A. Association of monocyte tumor necrosis factorαexpression and serum inflammatory biomarkers with walking impairment in peripheral artery disease. J. Vasc. Surg. 2015, 61, 155–161. [Google Scholar] [CrossRef] [Green Version]
  28. Tonelli, M.; Wiebe, N.; Fortin, M.; Guthrie, B.; Hemmelgarn, B.R.; James, M.T.; Klarenbach, S.W.; Lewanczuk, R.; Manns, B.J.; Ronksley, P.; et al. Methods for identifying 30 chronic conditions: Application to administrative data. BMC Med. Inform. Decis. Mak. 2015, 15, 31. [Google Scholar]
Figure 1. Consort diagram showing the detailed steps for assembling the two study cohorts.
Figure 1. Consort diagram showing the detailed steps for assembling the two study cohorts.
Ijerph 19 11327 g001
Figure 2. Cox model derived crude and multivariate adjusted hazard ratios (HRs) of PAD development stratified by patient characteristics in the osteoporosis cohort.
Figure 2. Cox model derived crude and multivariate adjusted hazard ratios (HRs) of PAD development stratified by patient characteristics in the osteoporosis cohort.
Ijerph 19 11327 g002
Figure 3. Cumulative incidence of peripheral artery diseases, which was 1.71% and 1.39%, respectively, in the osteoporosis and nonosteoporosis cohorts (p < 0.0001, compared to the log-rank).
Figure 3. Cumulative incidence of peripheral artery diseases, which was 1.71% and 1.39%, respectively, in the osteoporosis and nonosteoporosis cohorts (p < 0.0001, compared to the log-rank).
Ijerph 19 11327 g003
Table 1. Demographics and comorbidities at baseline between the osteoporosis cohort and the age-, sex-, and index date-matched comparison cohort without osteoporosis.
Table 1. Demographics and comorbidities at baseline between the osteoporosis cohort and the age-, sex-, and index date-matched comparison cohort without osteoporosis.
Case
(N = 64,562)
Comparison
(N = 64,562)
p
n(%)n(%)
Year
200623,95537.123,89837.020.6959
200721,32533.0321,24332.9
200819,28229.8719,42130.08
Sex
Female50,45478.1550,55578.30.4958
Male14,10821.8514,00721.7
Age
49–5511,05917.1311,05917.131
≥5553,50382.8753,50382.87
Mean (SD)66.25 (10.34)66.25 (10.34)1
Baseline Comorbidity
Hypertension62029.6149717.7<0.001
Hyperlipidemia22,04434.1417,54527.18<0.001
COPD18,19028.1711,77218.23<0.001
Diabetes mellitus16,02524.8212,70019.67<0.001
Liver disease15802.457951.23<0.001
CHF692610.7351307.95<0.001
RA1260.2870.130.0075
CKD39696.1521853.38<0.001
Gout11,22717.39798012.36<0.001
Ischemic stroke 28094.3512611.95<0.001
Raloxifene8061.2510241.59<0.001
Denosumab23933.7115472.4<0.001
Oral bisphosphate12031.8622603.5<0.001
Zoledronic acid5080.796230.960.0006
Ibandronate3420.533910.610.0695
Teriparatide2440.383380.520.0209
SD: standard deviation; COPD: chronic obstructive pulmonary disease; CHF: congestive heart failure; RA: rheumatic arthritis; CKD: chronic kidney disease.
Table 2. Incidence of peripheral artery diseases and the crude and adjusted hazard ratios (HRs) derived from the Cox model for the osteoporosis cohort compared to the comparison nonosteoporosis cohort stratified by patient characteristics.
Table 2. Incidence of peripheral artery diseases and the crude and adjusted hazard ratios (HRs) derived from the Cox model for the osteoporosis cohort compared to the comparison nonosteoporosis cohort stratified by patient characteristics.
Case
(N = 64,562)
Comparison
(N = 64,562)
N%Incidence RateN%Incidence Rate
PAD11051.71172.98981.39145.5
Case vs. Comparison
Crude Hazard Ratio (95% CI)p-ValueAdjusted HR (95% CI)p-Value
PAD1.21 (1.10–1.32)<0.0011.18(1.08–1.29)<0.001
PAD: peripheral artery disease; rate incidence per 100,000 PYs.
Table 3. Overall Cox model-derived crude- and multivariate-adjusted hazard ratios (HRs) of PAD development stratified by patient characteristics in both cohorts.
Table 3. Overall Cox model-derived crude- and multivariate-adjusted hazard ratios (HRs) of PAD development stratified by patient characteristics in both cohorts.
PAD
aHR (95% CI)p-Value
Gender <0.001
Male1
Female0.658 (0.595–0.727)
Age <0.001
49–551
≥551.729 (1.497–1.995)
Hypertension 0.5842
No1
Yes1.045 (0.893–1.223)
Hyperlipidemia 0.3151
No1
Yes0.949 (0.857–1.051)
COPD 0.2621
No1
Yes0.942 (0.848–1.046)
Diabetes <0.001
No1
Yes1.291 (1.162–1.434)
Liver cirrhosis 0.7197
No1
Yes0.94 (0.669–1.32)
CHF 0.1491
No1
Yes1.12 (0.96–1.307)
RA 0.0034
No1
Yes2.664 (1.382–5.135)
CKD 0.0029
No1
Yes1.558 (1.163–2.086)
Gout 0.1728
No1
Yes1.086 (0.964–1.224)
Ischemic stroke or SE 0.3575
No1
Yes0.887 (0.687–1.145)
Raloxifene 0.3124
No1
Yes0.815 (0.548–1.212)
Denosumab 0.2987
No1
Yes1.132 (0.896–1.432)
Oral bisphosphate 0.2789
No1
Yes1.151 (0.892–1.485)
Zoledronic acid 0.8917
No1
Yes0.969 (0.612–1.533)
Ibandronate 0.8755
No1
Yes0.955 (0.538–1.697)
Teriparatide 0.1937
No1
Yes1.433 (0.833–2.465)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Syu, D.-K.; Hsu, S.-H.; Yeh, P.-C.; Lee, T.-L.; Kuo, Y.-F.; Huang, Y.-C.; Jiang, C.-C.; Chen, M. The Association between Osteoporosis and Peripheral Artery Disease: A Population-Based Longitudinal Follow-Up Study in Taiwan. Int. J. Environ. Res. Public Health 2022, 19, 11327. https://doi.org/10.3390/ijerph191811327

AMA Style

Syu D-K, Hsu S-H, Yeh P-C, Lee T-L, Kuo Y-F, Huang Y-C, Jiang C-C, Chen M. The Association between Osteoporosis and Peripheral Artery Disease: A Population-Based Longitudinal Follow-Up Study in Taiwan. International Journal of Environmental Research and Public Health. 2022; 19(18):11327. https://doi.org/10.3390/ijerph191811327

Chicago/Turabian Style

Syu, De-Kai, Shu-Hua Hsu, Ping-Chun Yeh, Tsung-Lin Lee, Yu-Feng Kuo, Yen-Chun Huang, Ching-Chuan Jiang, and Mingchih Chen. 2022. "The Association between Osteoporosis and Peripheral Artery Disease: A Population-Based Longitudinal Follow-Up Study in Taiwan" International Journal of Environmental Research and Public Health 19, no. 18: 11327. https://doi.org/10.3390/ijerph191811327

APA Style

Syu, D. -K., Hsu, S. -H., Yeh, P. -C., Lee, T. -L., Kuo, Y. -F., Huang, Y. -C., Jiang, C. -C., & Chen, M. (2022). The Association between Osteoporosis and Peripheral Artery Disease: A Population-Based Longitudinal Follow-Up Study in Taiwan. International Journal of Environmental Research and Public Health, 19(18), 11327. https://doi.org/10.3390/ijerph191811327

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop