The influence of genetic susceptibility and calcium plus vitamin D supplementation on fracture risk

Anaphase-Promoting Complex Subunit 1 Associates with Bone Mineral Density in Human Osteoporotic Bone

Genome-wide association studies (GWAS) are one of the most common approaches to identify genetic loci that are associated with bone mineral density (BMD). Such novel genetic loci represent new potential targets for the prevention and treatment of fragility fractures. GWAS have identified hundreds of associations with BMD; however, only a few have been functionally evaluated. A locus significantly associated with femoral neck BMD at the genome-wide level is intronic SNP rs17040773 located in the intronic region of the anaphase-promoting complex subunit 1 (ANAPC1) gene (p = 1.5 × 10-9). Here, we functionally evaluate the role of ANAPC1 in bone remodelling by examining the expression of ANAPC1 in human bone and muscle tissues and during the osteogenic differentiation of human primary mesenchymal stem cells (MSCs). The expression of ANAPC1 was significantly decreased 2.3-fold in bone tissues and 6.2-fold in muscle tissue from osteoporotic patients as compared to the osteoarthritic and control tissues. Next, we show that the expression of ANAPC1 changes during the osteogenic differentiation process of human MSCs. Moreover, the silencing of ANAPC1 in human osteosarcoma (HOS) cells reduced RUNX2 expression, suggesting that ANAPC1 affects osteogenic differentiation through RUNX2. Altogether, our results indicate that ANAPC1 plays a role in bone physiology and in the development of osteoporosis.

Genome-wide polygenic risk score for major osteoporotic fractures in postmenopausal women using associated single nucleotide polymorphisms

BACKGROUND

Osteoporosis is highly polygenic and heritable, with heritability ranging from 50 to 80%; most inherited susceptibility is associated with the cumulative effect of many common genetic variants. However, existing genetic risk scores (GRS) only provide a few percent predictive power for osteoporotic fracture.

METHODS

We derived and validated a novel genome-wide polygenic score (GPS) comprised of 103,155 common genetic variants to quantify this susceptibility and tested this GPS prediction ability in an independent dataset (n = 15,776).

RESULTS

Among postmenopausal women, we found a fivefold gradient in the risk of major osteoporotic fracture (MOF) (p < 0.001) and a 15.25-fold increased risk of severe osteoporosis (p < 0.001) across the GPS deciles. Compared with the remainder of the GPS distribution, the top GPS decile was associated with a 3.59-, 2.48-, 1.92-, and 1.58-fold increased risk of any fracture, MOF, hip fracture, and spine fracture, respectively. The top GPS decile also identified nearly twofold more high-risk osteoporotic patients than the top decile of conventional GRS based on 1103 conditionally independent genome-wide significant SNPs. Although the relative risk of severe osteoporosis for postmenopausal women at around 50 is relatively similar, the cumulative incident at 20-year follow-up is significantly different between the top GPS decile (13.7%) and the bottom decile (< 1%). In the subgroup analysis, the GPS transferability in non-Hispanic White is better than in other racial/ethnic groups.

CONCLUSIONS

This new method to quantify inherited susceptibility to osteoporosis and osteoporotic fracture affords new opportunities for clinical prevention and risk assessment.

Estrogen receptor alpha and NFATc1 bind to a bone mineral density-associated SNP to repress WNT5B in osteoblasts

Genetic factors and estrogen deficiency contribute to the development of osteoporosis. The single-nucleotide polymorphism (SNP) rs2887571 is predicted from genome-wide association studies (GWASs) to associate with osteoporosis but has had an unknown mechanism. Analysis of osteoblasts from 110 different individuals who underwent joint replacement revealed that the genotype of rs2887571 correlates with WNT5B expression. Analysis of our ChIP-sequencing data revealed that SNP rs2887571 overlaps with an estrogen receptor alpha (ERα) binding site. Here we show that 17β-estradiol (E2) suppresses WNT5B expression and further demonstrate the mechanism of ERα binding at the enhancer containing rs2887571 to suppress WNT5B expression differentially in each genotype. ERα interacts with NFATc1, which is predicted to bind directly at rs2887571. CRISPR-Cas9 and ChIP-qPCR experiments confirm differential regulation of WNT5B between each allele. Homozygous GG has a higher binding affinity for ERα than homozygous AA and results in greater suppression of WNT5B expression. Functionally, WNT5B represses alkaline phosphatase expression and activity, decreasing osteoblast differentiation and mineralization. Furthermore, WNT5B increases interleukin-6 expression and suppresses E2-induced expression of alkaline phosphatase during osteoblast differentiation. We show that WNT5B suppresses the differentiation of osteoblasts via receptor tyrosine kinase-like orphan receptor 1/2 (ROR1/2), which activates DVL2/3/RAC1/CDC42/JNK/SIN3A signaling and inhibits β-catenin activity. Together, our data provide mechanistic insight into how ERα and NFATc1 regulate the non-coding SNP rs2887571, as well as the function of WNT5B on osteoblasts, which could provide alternative therapeutic targets for osteoporosis.

Habitual use of fish oil supplements, genetic predisposition, and risk of fractures: a large population-based study

BACKGROUND

Epidemiologic studies have suggested an inverse association between circulating concentrations of long-chain ω-3 PUFAs and fracture risk. However, whether supplementation of long-chain ω-3 PUFA (i.e. fish oil) is associated with fracture risk, and whether the association is modified by genetic predisposition to fracture risk remain unclear.

OBJECTIVES

To evaluate the associations of habitual fish oil supplement use with fracture risk, and to explore the potential effect modification by genetic predisposition.

METHODS

This study included 492,713 participants from the UK Biobank who completed a questionnaire on habitual fish oil supplement use between 2006 and 2010. HRs and 95% CIs for fractures were estimated from multivariable Cox proportional hazards models. A weighted fracture-genetic risk score (GRS) was derived from 14 validated single nucleotide polymorphisms.

RESULTS

During a median follow-up of 8.1 y, 12,070 incident fractures occurred among participants free of fracture at baseline (n = 441,756). Compared with nonuse, habitual use of fish oil supplements was associated with a lower risk of total fractures (HR = 0.93; 95% CI: 0.89, 0.97), hip fractures (HR = 0.83; 95% CI: 0.75, 0.92), and vertebrae fractures (HR = 0.85; 95% CI: 0.72, 0.99). The inverse association for total fractures was more pronounced among participants having a higher fracture-GRS than among those with a lower fracture-GRS (P-interaction <0.001). Among participants with a history of fracture at baseline (n = 50,957), fish oil use was associated with a lower risk of total recurrent fractures (HR = 0.88; 95% CI: 0.82, 0.96) and vertebrae recurrent fractures (HR = 0.64; 95% CI: 0.46, 0.88) but not with hip fracture recurrence.

CONCLUSIONS

Our findings suggest that habitual fish oil supplement use is associated with lower risks of both incident and recurrent fractures. The inverse associations of fish oil use with total fractures appeared to be more pronounced among individuals at higher genetic risk of fractures than those with lower genetic risk.

Relationship between Serum Nutritional Factors and Bone Mineral Density: A Mendelian Randomization Study

PURPOSE

Multiple risk factors have been implicated in the development of osteoporosis. This study examined potential associations between serum nutritional factors and bone mineral density (BMD).

METHODS

Six nutritional factors were selected as exposures. Outcomes included total body BMD (n = 66 945); BMD at the forearm (FA), femoral neck (FN) and lumbar spine (LS) (n = 8143, n = 32 735, and n = 28 498, respectively); estimated heel BMD (HL eBMD) (n = 394 929); and HL eBMD stratified by sex (n = 206 496). A 2-sample Mendelian randomization approach was adopted to estimate the association between serum nutritional factors and BMD. The threshold for adjusted P value was 1.39 × 10-3.

RESULTS

Serum calcium levels were inversely associated with LS BMD (effect = -0.55; 95% CI, -0.86 to -0.24; P = 0.001), whereas serum selenium levels were positively correlated with HL eBMD (effect = 0.22; 95% CI, 0.10 to 0.33; P = 1.70 × 10-4). Regarding nominal significance, there was a positive association between serum selenium levels and FA BMD. Nominally significant results were also obtained for serum retinol as well as vitamin E levels and HL eBMD. Moreover, sex-specific effects of serum retinol and vitamin E levels on BMD were observed in men.

CONCLUSION

Serum calcium and selenium levels influence BMD at specific skeletal sites. This implies that these nutritional factors play crucial roles in bone metabolism.

Vitamin D supplementation and fracture risk: Evidence for a U-shaped effect

During the last decade, a cascade of evidence has questioned the anti-fracture efficacy of vitamin D supplementation. In general, vitamin D status, reflected by serum 25-hydroxy-vitamin D [25(OH)D] concentrations, seems to predict fracture risk and bone mineral density (BMD). Despite the well-documented detrimental effect of vitamin D deficiency on bones, vitamin D monotherapy does not seem to reduce the risk of fractures. On the other hand, high vitamin D doses, either at monthly (60,000-100,000 IU) or daily intervals (>4000 IU), appear to be harmful with regard to falls, fracture risk and BMD, especially for people without vitamin D deficiency and at low fracture risk. Therefore, a U-shaped effect of vitamin D on the musculoskeletal system may be supported by the current evidence. Vitamin D supplementation could be of value, at daily doses of at least 800 IU, co-supplemented with calcium (1000-1200  mg/day), in elderly populations, especially those with severe vitamin D deficiency [25(OH)D <25-30  nmol/L (<10-12  ng/mL)], although its anti-fracture and anti-fall efficacy is modest. Good compliance and at least 3-5 years of therapy are required.

Effects of vitamin D2-fortified shiitake mushroom on bioavailability and bone structure

Bioavailability and bone loss inhibitory effects of vitamin D₂ derived from UV-irradiated shiitake mushroom were determined in vivo. The effect of the absence of ovaries on the bioavailability of vitamin D₂ and bone structure was also investigated. Sham operated (sham) and ovariectomized (OVX) rats were divided in 3 groups according to their diets, i.e. control: only vitamin D-deficient diets; UV(X): vitamin D-deficient diets with non-irradiated mushroom powder; UV(O): vitamin D-deficient diets with irradiated mushroom powder. The obtained results showed that vitamin D₂ from shiitake mushroom was able to increase bone mineral density and trabecular bone structure of femur bone as well as its bioavailability. The absence of estrogen induced adverse effects not only on bioavailability of vitamin D₂ but also on trabecular bone. In conclusion, vitamin D₂-fortified shiitake mushroom might help postmenopausal women increase vitamin D₂ bioavailability and retard trabecular bone loss. Abbreviations: OVX: ovariectomized; 25(OH)D: 25-hydroxyvitamin D; 1,25(OH)2D: 1,25-dihydroxyvitamin D; BMD: bone mineral density; micro-CT: micro computed tomography; RSM: response surface methodology; RP-HPLC: Reverse phase-high performance liquid chromatography; MS/MS: tandem mass spectrometry; E2: estradiol; NTx: N-terminal telopeptide of type I collagen; BV/TV: bone volume/total volume; BS/BV: bone surface/bone volume; Tb.Th: trabecular thickness; Tb.Sp: trabecular separation.

Estimating the Relative Excess Risk Due to Interaction in Clustered-Data Settings

The risk difference scale is often of primary interest when evaluating public health impacts of interventions on binary outcomes. However, few investigators report findings in terms of additive interaction, probably because the models typically used for binary outcomes implicitly measure interaction on the multiplicative scale. One measure with which to assess additive interaction from multiplicative models is the relative excess risk due to interaction (RERI). The RERI measure has been applied in many contexts, but one limitation of previous approaches is that clustering in data has rarely been considered. We evaluated the RERI metric for the setting of clustered data using both population-averaged and cluster-conditional models. In simulation studies, we found that estimation and inference for the RERI using population-averaged models was straightforward. However, frequentist implementations of cluster-conditional models including random intercepts often failed to converge or produced degenerate variance estimates. We developed a Bayesian implementation of log binomial random-intercept models, which represents an attractive alternative for estimating the RERI in cluster-conditional models. We applied the methods to an observational study of adverse birth outcomes in mothers with human immunodeficiency virus, in which mothers were clustered within clinical research sites.