Bin He
Lifeng Yin1
Zhengxue Quan- 1Department of Orthopedics, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
- 2Department of General Practice, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
Background: Hypertension may have some association with osteoporosis. This Mendelian randomization (MR) study aimed to explore the causal effect of blood pressure (BP) on bone mineral density (BMD), fall, and fracture.
Methods: We used the genome-wide association study (GWAS) summary data among 330,956 European-descent individuals to identify 107 single-nucleotide polymorphisms (SNPs) as the instrumental variables of BP. MR analyses of these instruments were performed on 53,236 European individuals for the association with forearm BMD (FA-BMD), femoral neck BMD (FN-BMD), and lumbar spine BMD (LS-BMD); 451,179 European individuals for fall susceptibility; and up to 1.2 million individuals from European descent for fracture. Conventional inverse variance weighted (IVW) method was adopted to obtain the causal estimates of BP on different outcomes, while weighted median, MR-egger, and MR pleiotropy residual sum and outlier (MR-PRESSO) test were used for sensitivity analyses.
Results: Genetically high pulse pressure (PP) could significantly improve FA-BMD (beta-estimate: 0.038, 95% confidence interval [CI]: 0.013 to 0.063, SE:0.013, P-value=0.003<Bonferroni correction P) in the IVW analysis, indicating that 1-SD increase in PP was associated with the improvement in FA-BMD levels by 0.038 g/cm2 (95% CI: 0.013 to 0.063). This positive finding was also confirmed by weighted-median analysis (beta-estimate: 0.034, 95% CI: 0.000 to 0.067, SE:0.017, P-value=0.046) and MR-Egger analysis (beta-estimate: 0.117, 95% CI: 0.026 to 0.208, SE:0.046, P-value=0.011). However, there was no remarkable MR association between BP and other outcomes (i.e., FN-BMD, LS-BMD, fall, and fracture).
Conclusions: Our findings reveal a potentially causal relationship between high PP and improved FA-BMD, which may provide new sights for the treatment of osteoporosis.
Introduction
Osteoporosis is one common systemic skeletal disease characterized by decreased bone mineral density (BMD) and increased risk of fracture (1–3). Its prevention and treatment are still a big challenge and growing public health problem in the world (4–6). Genome-wide association study (GWAS) has demonstrated that BMD is a highly polygenic trait, and some genetic determinants of fracture act through low BMD (7–9).
Hypertension is one common and strongly heritable disease, with high mortality and morbidity (10–12). It may also increase the risk of stroke and coronary artery disease (12, 13). Several observational studies reported the association between hypertension and BMD, but their results were conflicting (14–16). Mendelian randomization (MR) has become an effective and powerful approach to establish the causal relationships between exposure phenotype and outcome phenotype (17, 18). The use of MR study can overcome the limitations of confounding factors and reverse causation bias that commonly occur in observational studies (19, 20).
In this two-sample MR analysis, we used single-nucleotide polymorphisms (SNPs) strongly associated with systolic blood pressure (SBP), diastolic blood pressure (DBP), and pulse pressure (PP) as instrumental variables. To our knowledge, this is the first two-sample MR study to explore the causal effect of blood pressure (BP) on forearm BMD (FA-BMD), femoral neck BMD (FN-BMD), lumbar spine BMD (LS-BMD), fall, and fracture.
Materials and Methods
This MR study was conducted based on the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline (Supplementary Table 1).
GWAS Summary Statistics of Blood Pressure
In the meta-analysis of genome-wide association study (GWAS) among 330,956 European-descent individuals, a threshold of P < 5*10−8 was used to denote genome-wide significance of blood pressure. The replication resources included a large blood pressure (BP) meta-analysis consortium (ICBP cohorts) and further cohorts with 1,000 Genomes data for GWAS (21). Blood pressure was defined as the mean arterial pressure with the lowest mean velocity index and estimated by the Omron device. Initially, 107 independent single-nucleotide polymorphisms (SNPs) were identified to have robust association with blood pressure at the GWAS threshold of statistical significance (P<5×10−8), including 24 SNPs for systolic blood pressure (SBP), 41 SNPs for diastolic blood pressure (DBP), and 42 SNPs for pulse pressure (PP) (Supplementary Table 2). These SNPs were all adjusted by the Principal Component Analysis in order to address the population stratification.
We excluded SNPs in strong linkage disequilibrium (LD) because they produce some bias. The clumping process (R2<0.001, window size =10,000 kb) was conducted with the European samples from the 1,000 genomes project, and we estimated LD between SNPs. Among the pairs of SNPs with r2≥0.001, the SNP with a larger association P value would be excluded. The SNPs that were absent from the LD reference panel were also removed. Thus, three SNPs about SBP, four SNPs about DBP, and five SNPs about PP were excluded due to LD. Finally, 21 SNPs for SBP, 37 SNPs for DBP, and 37 SNPs for PP were used as the instrumental variables (Supplementary Table 3). If SNPs were unavailable in the outcome dataset, the proxy SNPs in linkage disequilibrium (LD, r2>0.9) were used as the instrumental variables.
Outcome Data Sources
A large meta-analysis was conducted to explore the genetic variants associated with FA-BMD, FN-BMD, and LS-BMD among 53,236 individuals of European ancestry (22). BMD was measured at the trabecular structure of forearm (distal 1/3 of radius), femoral neck, and lumbar spine (L1-4). BMD was measured by dual X-ray absorptiometry, and low BMD was defined as the Z-score <−1.0. Each SNP was tested after adjusting for sex, age, age2, and weight (22). A genome-wide association analysis involved 89,076 fall cases and 362,103 controls from the UK Biobank Study of European ancestry, and revealed the genetic determinants of fall susceptibility. Fall cases were defined as participants who gave positive answer to the following question: “In the last year have you had any falls?” (23). Genetic determinants of fracture risk were revealed in one large GWAS meta-analysis in up to 1.2 million individuals combining the UK Biobank and 23andMe cohorts. Fractures were defined as a break in the continuity of the bone at any site except the fractures of skull, face, hands and feet, pathological fractures due to malignancy, atypical femoral fractures, periprosthetic and healed fracture codes. The diagnosis of fracture should be within the past five years. (1).
Statistical Analysis
We evaluated the causal effect of BP (SBP, DBP, and PP) on BMD (FA-BMD, FN-BMD, LS-BMD), fall, and fracture. MR estimates for instrumental variables were meta-analyzed by computing an inverse variance weighted (IVW) analysis for the primary analysis (24). We used Cochran’s Q analysis to assess the heterogeneity (25), where high heterogeneity indicated the presence of invalid genetic variants (26). For the sensitivity analysis, weighted-median analysis was conducted, which provided a valid estimate if at least 50% of weight originated from non-pleiotropic SNPs (27).
To assess the potential violation of these assumptions, MR-Egger analysis was used to assess the directional pleiotropy based on the intercept (28). The presence of pleiotropy was also assessed by the MR pleiotropy residual sum and outlier test (MR-PRESSO), during which outlying SNPs were excluded and the effect estimates were reassessed (29). The ethical approval and informed consent for each study included in the study can be found in the original publications. The differences with P<0.05 were considered statistically significant. In multiple testing, an adjusted P value after Bonferroni correction (P<0.05/4 = 0.0125) was considered statistically significant. All of the analyses were conducted in R V.4.0.4 by using the R packages of “MendelianRandomization” (30), “TwoSampleMR” (31), and “MR-PRESSO” (29).
Role of the Funding Source
The funders of this study had an important role in study design, data collection, data analysis, data interpretation, and writing of the report. All authors had full access to all data in the study and had final responsibility for the decision to submit for publication.
Results
Causal Effect of Blood Pressure on BMD
We evaluated the causal effect of blood pressure including SBP (Figure 1), DBP (Figure 2), and PP (Figure 3) on FA-BMD, FN-BMD, and LS-BMD in the MR analysis (Table 1). High PP was significantly associated with improved FA-BMD (beta-estimate: 0.038, 95% confidence interval [CI]: 0.013 to 0.063, standard error [SE]: 0.013, P-value=0.003<Bonferroni correction P, Figure 3) in the IVW analysis, and this positive finding was supported by weighted-median analysis (beta-estimate: 0.034, 95% CI: 0.000 to 0.067, SE:0.017, P-value=0.046) and MR-Egger analysis (beta-estimate: 0.117, 95% CI: 0.026 to 0.208, SE:0.046, P-value=0.011). However, PP showed no remarkable influence on FN-BMD or LS-BMD based on the results of IVW, weighted-median, and MR-Egger analyses.
Figure 1 Mendelian randomization estimates for the associations between SBP and outcomes. FA-BMD, forearm BMD; FN-BMD, femoral neck BMD; LS-BMD, lumbar spine BMD; CI, confidence interval.
Figure 2 Mendelian randomization estimates for the associations between DBP and outcomes. FA-BMD, forearm BMD; FN-BMD, femoral neck BMD; LS-BMD, lumbar spine BMD; CI, confidence interval.
Figure 3 Mendelian randomization estimates for the associations between PP and outcomes. FA-BMD, forearm BMD; FN-BMD, femoral neck BMD; LS-BMD, lumbar spine BMD; CI, confidence interval.
Table 1 Mendelian randomization estimates of blood pressure on outcomes.
According to primary IVW analyses, SBP showed no MR association with FA-BMD (beta-estimate: −0.015, 95% CI: −0.041 to 0.011, SE:0.013, P-value=0.258), FN-BMD (beta-estimate: 0.010, 95% CI: −0.001 to 0.021, SE:0.006, P-value=0.064), or LS-BMD (beta-estimate: 0.001, 95% CI: −0.014 to 0.016, SE:0.008, P-value=0.921, Figure 1), while DBP also demonstrated no obvious impact on FA-BMD (beta-estimate: 0.002, 95% CI: −0.026 to 0.030, SE:0.014, P-value=0.876), FN-BMD (beta-estimate: −0.012, 95% CI: −0.025 to 0.001, SE:0.007, P-value=0.077), or LS-BMD (beta-estimate: −0.004, 95% CI: −0.019 to 0.012, SE:0.008, P-value=0.666, Figure 2) during the IVW analyses. These results were also confirmed by the weighted-median analysis and MR-Egger analysis.
Causal Effect of Blood Pressure on Fall and Fracture
SBP, DBP, and PP showed null association with fall in the IVW (odds ratio [OR]: 1.002, 95% CI: 0.994 to 1.010; SE:0.004, P-value=0.620 for SBP, Figure 1; OR: 1.003, 95% CI: 0.991 to 1.015; SE:0.006, P-value=0.577 for DBP, Figure 2; OR: 0.996, 95% CI: 0.988 to 1.004; SE:0.004, P-value=0.368 for PP, Figure 3). Consistently, there was also no relationship between BP and fracture in the IVW analysis (OR: 1.003, 95% CI: 0.993 to 1.013; SE:0.005, P-value=0.258 for SBP, Figure 1; OR: 1.005, 95% CI: 0.993 to 1.017; SE:0.006, P-value=0.434 for DBP, Figure 2; OR: 1.007, 95% CI: 0.995 to 1.019; SE:0.006, P-value=0.246 for PP, Figure 3). These results were all confirmed by the weighted-median analysis and MR-Egger analysis.
Evaluation of Assumptions and Sensitivity Analyses
Little evidence of directional pleiotropy was found for all models except for the association between DBP and FN-BMD (MR-Egger intercept P-value=0.03) (Table 1). The estimates from the weighted-median approach and MR-Egger analysis were all consistent with those of IVW models (Table 1).
Among the instrument variables, MR-PRESSO method only identified one outlier (rs72799341) for the association between DBP and fall, and one outlier (rs12628032) for the association between PP and fall. After excluding these outliers, DBP and PP still revealed no causal effect on the incidence of fall (OR: 1.002, 95% CI: 0.990 to 1.014, SE:0.006, P-value=0.801 for DBP; OR: 0.993, 95% CI: 0.985 to 1.001, SE:0.004, P-value=0.072 for PP, Table 2).
Table 2 Mendelian randomization estimates between blood pressure and outcomes after excluding outliers detected by MR-PRESSO.
Discussion
In this two-sample MR analysis, we found the casual effect between high PP and improved FA-BMD (beta-estimate: 0.038, 95% CI: 0.013 to 0.063, SE:0.013, P-value=0.003<Bonferroni correction P), indicating that 1-SD increase in PP was associated with the improvement in FA-BMD levels by 0.038 g/cm2 (95% CI: 0.013 to 0.063). This positive finding was also confirmed by weighted median, MR-egger, and MR-PRESSO analyses. However, no causal association was seen between BP and other outcomes (i.e., FN-BMD, LS-BMD, fall, and fracture).
Hypertension has important association with alterations in calcium metabolism, including increased calcium loss, compensatory activation of parathyroid gland, and increased movement of calcium from the bones (32). Long-lasting impairment effect of hypertension on calcium homeostasis may result in age-related excessive reduction of BMD and fracture (32). Previous studies explored the association between BP and BMD, but reported conflicting results (14–16).
In one cross-sectional study involving 270 postmenopausal Turkish women, hypertension was found to be significant predictors of osteopenia in a multivariate analysis (OR: 2.541, 95% CI: 1.46–3.48, P=0.003) (14). A retrospective analysis of 586 postmenopausal women with a mean age of 60.8 ± 8.8 years old revealed that hypertension was associated with low spine BMD in postmenopausal women (16). In contrast, another cross-sectional study was conducted among 4,058 premenopausal and postmenopausal women aged 40 years or older (number=991 and 3,067, respectively), and the results revealed no link between BP (i.e., SBP and DBP) and BMD (15). These studies did not involve male patients and were indeed contradictory. These inconsistent results may result from small patient sample, confounding factors, and the limitation of study design.
Our two-sample MR analysis involved 53,236 European individuals for the association with BMD, 451,179 European individuals for fall susceptibility, and up to 1.2 million individuals from European descent for fracture. The results revealed the important causal effect of high PP on improved FA-BMD (P=0.003), which was confirmed by multiple sensitivity analyses. The association between PP and BMD was explored after adjusting for sex, age, age2, and weight (22). In addition, we found no obvious MR association between BP and other outcomes (i.e., FN-BMD, LS-BMD, fall, and fracture).
It is very interesting to confirm that genetically high PP shows strong MR association with improved BMD. PP is defined as SBP minus DBP. High PP largely results from large-artery stiffness, while decreased PP is caused by low stroke volume, such as congestive heart failure and aortic valve stenosis (33). Vascular smooth muscle cells (VSMC) have important roles in regulating arterial stiffness by overproducing various extracellular matrix components (e.g., collagen and elastin), which provide biomechanical, structural integrity, and signaling regulation of the extracellular matrix components to maintain vascular homeostasis (34, 35). However, patients with high PP have vascular calcification and increased arterial stiffness, which increase the expression of bone markers such as alkaline phosphatase (ALP) and type 1 collagen. These factors also improve the osteogenic differentiation and mineralization for bone formation and increased BMD (36).
Forearm BMD at the 1/3 radial site is commonly used to improve the prediction of hip fractures when considered together with FN-BMD. Forearm BMD represents the BMD of combined trabecular and cortical bone structures, indicating the better improvement of BMD than other sites of BMD after effective anti-osteoporosis treatments, which may account for the positive MR association only between PP and forearm BMD (37). In addition, 1,032 men and 1,701 women aged 50 years and older were included in the Dubbo Osteoporosis Epidemiology Study, and hypertension may be an independent risk factor for fragility fracture (hazard ratio, 1.49; 95% CI, 1.13–1.96) after adjusting for BMD and covariates (38). However, our two-sample MR analysis confirmed null association between BP and fracture.
To the best of our knowledge, this is the first two-sample MR study to find the positively casual association between PP and FA-BMD. The summary statistics of outcome phenotypes are retrieved from GWAS or genome-wide meta-analysis with huge sample size. This two-sample MR method allows for the estimation of the causal effects of BP on all outcomes while at the same time minimizing reverse causation bias and confounding factors. The intercepts for the MR-Egger analysis, except for the association between DBP and FN-BMD, indicate that there is no directional pleiotropy to influence other causal associations.
Several limitations should be taken into consideration. Firstly, there is some heterogeneity between DBP (or PP) and fall, which may be caused by the selection of some instrumental variables. Secondly, this MR study reveals the potential causal effect of PP on FA-BMD, but null association is observed between BP and FN-BMD (or LS-BMD). The factors of this inconsistency remain elusive. Thirdly, it is not feasible to perform the MR analysis based on different age stratums because of the limitation of GWAS summary statistics. Fourthly, the MR analysis is restricted in the European-ancestry population, which may limit the generalizability of our finding to other populations.
Conclusion
This two-sample MR reveals the potential causal effect of high PP on improved FA-BMD, suggesting the protective role of high PP for osteoporosis.
Data Availability Statement
The original data were available in GEnetic Factors for OSteoporosis Consortium (http://www.gefos.org/) and the UK Biobank (https://www.ukbiobank.ac.uk/). The instructions to download data from (http://www.gefos.org/) are presented in Datasheet 3.
Author Contributions
BH, ZQ, and QL conducted study design. BH, LFY and MZZ conducted data collection and statistical analysis. BH, ZQ, QL, and YO conducted data interpretation, manuscript preparation, and literature search. BH and QL conducted funds collection. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by Natural Science Foundation of Chongqing (cstc2019jcyj-msxmX0836), National Natural Science Foundation of China (81701382), Foundation of The First Affiliated Hospital of Chongqing Medical University (PYJJ2018-13), and Chongqing Yuzhong Nature Science Foundation of China (Grant No. 2018114).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors acknowledged the GEnetic Factors for OSteoporosis Consortium and the UK Biobank for contributing the data used in this work. We thanked all the genetics consortiums for making the GWAS summary data publicly available.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2021.716681/full#supplementary-material
References
1. Morris JA, Kemp JP, Youlten SE, Laurent L, Logan JG, Chai RC, et al. An Atlas of Genetic Influences on Osteoporosis in Humans and Mice. Nat Genet (2019) 51(2):258–66. doi: 10.1038/s41588-018-0302-x
2. Liu J, Curtis EM, Cooper C, Harvey NC. State of the Art in Osteoporosis Risk Assessment and Treatment. J Endocrinol Invest (2019) 42(10):1149–64. doi: 10.1007/s40618-019-01041-6
3. Leder BZ, Mitlak B, Hu MY, Hattersley G, Bockman RS. Effect of Abaloparatide vs Alendronate on Fracture Risk Reduction in Postmenopausal Women With Osteoporosis. J Clin Endocrinol Metab (2020) 105(3):938–43. doi: 10.1210/clinem/dgz162
4. Reid IR. A Broader Strategy for Osteoporosis Interventions, Nature Reviews. Endocrinology (2020) 16(6):333–9. doi: 10.1038/s41574-020-0339-7
5. Khosla S, Hofbauer LC. Osteoporosis Treatment: Recent Developments and Ongoing Challenges. Lancet Diabetes Endocrinol (2017) 5(11):898–907. doi: 10.1016/S2213-8587(17)30188-2
6. Compston JE, McClung MR, Leslie WD. Osteoporosis. Lancet (London England) (2019) 393(10169):364–76. doi: 10.1016/S0140-6736(18)32112-3
7. Richards JB, Zheng HF, Spector TD. Genetics of Osteoporosis From Genome-Wide Association Studies: Advances and Challenges, Nature Reviews. Genetics (2012) 13(8):576–88. doi: 10.1038/nrg3228
8. Trajanoska K, Rivadeneira F. The Genetic Architecture of Osteoporosis and Fracture Risk. Bone (2019) 126:2–10. doi: 10.1016/j.bone.2019.04.005
9. Yang TL, Shen H, Liu A, Dong SS, Zhang L, Deng FY, et al. A Road Map for Understanding Molecular and Genetic Determinants of Osteoporosis, Nature Reviews. Endocrinology (2020) 16(2):91–103. doi: 10.1038/s41574-019-0282-7
10. Muñoz M, Pong-Wong R, Canela-Xandri O, Rawlik K, Haley CS, Tenesa A. Evaluating the Contribution of Genetics and Familial Shared Environment to Common Disease Using the UK Biobank. Nat Genet (2016) 48(9):980–3. doi: 10.1038/ng.3618
11. Poulter NR, Prabhakaran D, Caulfield M. Hypertension. Lancet (London England) (2015) 386(9995):801–12. doi: 10.1016/S0140-6736(14)61468-9
12. Forouzanfar MH, Alexander L, Anderson HR, Bachman VF, Biryukov S, Brauer M, et al. Global, Regional, and National Comparative Risk Assessment of 79 Behavioural, Environmental and Occupational, and Metabolic Risks or Clusters of Risks in 188 Countries, 1990-2013: A Systematic Analysis for the Global Burden of Disease Study 2013. Lancet (London England) (2015) 386(10010):2287–323. doi: 10.1016/S0140-6736(15)00128-2
13. Karmali KN, Lloyd-Jones DM, van der Leeuw J, Goff DC Jr., Yusuf S, Zanchetti A, et al. Blood Pressure-Lowering Treatment Strategies Based on Cardiovascular Risk Versus Blood Pressure: A Meta-Analysis of Individual Participant Data. PloS Med (2018) 15(3):e1002538. doi: 10.1371/journal.pmed.1002538
14. Cakmak HA, Cakmak BD, Yumru AE, Aslan S, Enhos A, Kalkan AK, et al. The Relationships Between Blood Pressure, Blood Glucose, and Bone Mineral Density in Postmenopausal Turkish Women. Ther Clin Risk Manage (2015) 11:1641–8. doi: 10.2147/TCRM.S95017
15. Kaplan S, Smith SR, Zuckerman IH. Blood Pressure and Bone Mineral Density in Premenopausal and Postmenopausal Women. J Women’s Health (2002) (2010) 19(6):1209–15. doi: 10.1089/jwh.2009.1587
16. Yazici S, Yazici M, Korkmaz U, Engin Erkan M, Erdem Baki A, Erden I, et al. Relationship Between Blood Pressure Levels and Bone Mineral Density in Postmenopausal Turkish Women. Arch Med Sci: AMS (2011) 7(2):264–70. doi: 10.5114/aoms.2011.22077
17. Burgess S, Dudbridge F, Thompson SG. Combining Information on Multiple Instrumental Variables in Mendelian Randomization: Comparison of Allele Score and Summarized Data Methods. Stat Med (2016) 35(11):1880–906. doi: 10.1002/sim.6835
18. Dalbeth N, Topless R, Flynn T, Cadzow M, Bolland MJ, Merriman TR. Mendelian Randomization Analysis to Examine for a Causal Effect of Urate on Bone Mineral Density. J Bone Mineral Research: Off J Am Soc Bone Mineral Res (2015) 30(6):985–91. doi: 10.1002/jbmr.2434
19. Davey Smith G, Ebrahim S. What can Mendelian Randomisation Tell Us About Modifiable Behavioural and Environmental Exposures? BMJ (Clinical Res ed) (2005) 330(7499):1076–9. doi: 10.1136/bmj.330.7499.1076
20. Davies NM, Holmes MV, Davey Smith G. Reading Mendelian Randomisation Studies: A Guide, Glossary, and Checklist for Clinicians. BMJ (Clinical Res ed) (2018) 362:k601. doi: 10.1136/bmj.k601
21. Warren HR, Evangelou E, Cabrera CP, Gao H, Ren M, Mifsud B, et al. Genome-Wide Association Analysis Identifies Novel Blood Pressure Loci and Offers Biological Insights Into Cardiovascular Risk. Nat Genet (2017) 49(3):403–15. doi: 10.1038/ng.3768
22. Zheng HF, Forgetta V, Hsu YH, Estrada K, Rosello-Diez A, Leo PJ, et al. Whole-Genome Sequencing Identifies EN1 as a Determinant of Bone Density and Fracture. Nature (2015) 526(7571):112–7. doi: 10.1038/nature14878
23. Trajanoska K, Seppala LJ, Medina-Gomez C, Hsu YH, Zhou S, van Schoor NM, et al. Genetic Basis of Falling Risk Susceptibility in the UK Biobank Study. Commun Biol (2020) 3(1):543. doi: 10.1038/s42003-020-01256-x
24. Burgess S, Butterworth A, Thompson SG. Mendelian Randomization Analysis With Multiple Genetic Variants Using Summarized Data. Genet Epidemiol (2013) 37(7):658–65. doi: 10.1002/gepi.21758
25. Greco MF, Minelli C, Sheehan NA, Thompson JR. Detecting Pleiotropy in Mendelian Randomisation Studies With Summary Data and a Continuous Outcome. Stat Med (2015) 34(21):2926–40. doi: 10.1002/sim.6522
26. Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N, Thompson J. A Framework for the Investigation of Pleiotropy in Two-Sample Summary Data Mendelian Randomization. Stat Med (2017) 36(11):1783–802. doi: 10.1002/sim.7221
27. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent Estimation in Mendelian Randomization With Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol (2016) 40(4):304–14. doi: 10.1002/gepi.21965
28. Burgess S, Thompson SG. Interpreting Findings From Mendelian Randomization Using the MR-Egger Method. Eur J Epidemiol (2017) 32(5):377–89. doi: 10.1007/s10654-017-0255-x
29. Verbanck M, Chen CY, Neale B, Do R. Detection of Widespread Horizontal Pleiotropy in Causal Relationships Inferred From Mendelian Randomization Between Complex Traits and Diseases. Nat Genet (2018) 50(5):693–8. doi: 10.1038/s41588-018-0099-7
30. Yavorska OO, Burgess S. MendelianRandomization: An R Package for Performing Mendelian Randomization Analyses Using Summarized Data. Int J Epidemiol (2017) 46(6):1734–9. doi: 10.1093/ije/dyx034
31. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base Platform Supports Systematic Causal Inference Across the Human Phenome. eLife (2018) 7:e34408. doi: 10.7554/eLife.34408
32. Cappuccio FP, Meilahn E, Zmuda JM, Cauley JA. High Blood Pressure and Bone-Mineral Loss in Elderly White Women: A Prospective Study. Study of Osteoporotic Fractures Research Group. Lancet (London England) (1999) 354(9183):971–5. doi: 10.1016/S0140-6736(99)01437-3
33. McFarlane IM, Shin TH, Bhamra M, Alvarez MR, Leon SZ, Ozeri DJ, et al. The Relationship of Pulse Pressure and Bone Mineral Density in Adult USA Population: Analysis of the National Health and Nutritional Examination Survey. Rheumatol (Sunnyvale Calif) (2018) 8(2):240. doi: 10.4172/2161-1149.1000240
34. Cocciolone AJ, Hawes JZ, Staiculescu MC, Johnson EO, Murshed M, Wagenseil JE. Elastin, Arterial Mechanics, and Cardiovascular Disease, American Journal of Physiology. Heart Circulatory Physiol (2018) 315(2):H189–h205. doi: 10.1152/ajpheart.00087.2018
35. Bäck M, Gasser TC, Michel JB, Caligiuri G. Biomechanical Factors in the Biology of Aortic Wall and Aortic Valve Diseases. Cardiovasc Res (2013) 99(2):232–41. doi: 10.1093/cvr/cvt040
36. Chen Y, Zhao X, Wu H. Arterial Stiffness: A Focus on Vascular Calcification and Its Link to Bone Mineralization. Arteriosclerosis Thromb Vasc Biol (2020) 40(5):1078–93. doi: 10.1161/ATVBAHA.120.313131
37. Biver E, Durosier-Izart C, Chevalley T, van Rietbergen B, Rizzoli R, Ferrari S. Evaluation of Radius Microstructure and Areal Bone Mineral Density Improves Fracture Prediction in Postmenopausal Women. J Bone Mineral Research: Off J Am Soc Bone Mineral Res (2018) 33(2):328–37. doi: 10.1002/jbmr.3299
38. Yang S, Nguyen ND, Center JR, Eisman JA, Nguyen TV. Association Between Hypertension and Fragility Fracture: A Longitudinal Study. Osteoporosis International: A Journal Established as Result Cooperation Between the European Foundation for Osteoporosis and the National Osteoporosis Foundation USA (2014) 25(1):97–103. doi: 10.1007/s00198-013-2457-8
Keywords: blood pressure, bone mineral density, fall, fracture, Mendelian randomization study
Citation: He B, Yin L, Zhang M, Lyu Q, Quan Z and Ou Y (2021) Causal Effect of Blood Pressure on Bone Mineral Density and Fracture: A Mendelian Randomization Study. Front. Endocrinol. 12:716681. doi: 10.3389/fendo.2021.716681
Received: 29 May 2021; Accepted: 12 July 2021;
Published: 04 August 2021.
Edited by:
Melissa Orlandin Premaor, Federal University of Minas Gerais, BrazilReviewed by:
Enisa Shevroja, Centre Hospitalier Universitaire Vaudois (CHUV), SwitzerlandChao Xu, University of Oklahoma Health Sciences Center, United States
Copyright © 2021 He, Yin, Zhang, Lyu, Quan and Ou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiong Lyu, 425598350@qq.com; Zhengxue Quan, quanzx18@126.com