Assessment of the genetic and clinical determinants of fracture risk: genome wide association and mendelian randomisation study

Mendelian randomization in the bone field

Identification of causative risk factors amenable for modification is essential for the prevention and treatment of osteoporosis. Observational studies have identified associations between several potentially modifiable risk factors and osteoporosis. However, observational studies are susceptible to confounding, reverse causation bias, and measurement error, all of which limit their ability to provide causal estimates of the effect of exposures on outcomes, thereby reducing their ability to inform prevention and treatment strategies against bone loss and fractures. In addition, not all risk factors are suitable for an analysis in a randomized clinical trial. Mendelian randomization is a genetic epidemiological method that exploits genetic variants as unbiased proxies for modifiable exposures (e.g., biomarkers, adiposity measures, dietary factors, and behaviors) to determine the causal relationships between exposures and health outcomes. This technique has been used to provide evidence of causal associations of serum estradiol concentrations, smoking, body mass index, and type 2 diabetes with bone mineral density and the lack of associations of serum thyroid stimulating hormone, urate, C-reactive protein, and 25‑hydroxyvitamin D concentrations with bone mineral density in generally healthy populations. This review will briefly explain the concept of Mendelian randomization, the advantages and potential limitations of this study design, and give examples of how Mendelian randomization has been used to investigate questions relevant to osteoporosis.

The genetic architecture of osteoporosis and fracture risk

Osteoporosis and fracture risk are common complex diseases, caused by an interaction of numerous disease susceptibility genes and environmental factors. With the advances in genomic technologies, large-scale genome-wide association studies (GWAS) have been performed which have broadened our understanding of the genetic architecture and biological mechanisms of complex disease. Currently, more than ~90 loci have been found associated with DXA derived bone mineral density (BMD), over ~500 loci with heel estimated BMD and several others with other less widely available bone parameters such as bone geometry, shape, and microarchitecture. Notably, several of the pathways identified by the GWAS efforts correspond to pathways that are currently targeted for the treatment of osteoporosis. Overall, tremendous progress in the field of the genetics of osteoporosis has been achieved with the discovery of WNT16, EN1, DAAM2, and GPC6 among others. Assessment of the function and biological mechanisms of the remaining genes may further untangle the complex genetic landscape of osteoporosis and fracture risk. With this review we aimed to provide a general overview of the existing GWAS studies on osteoporosis traits and fracture risk.

Can Classification and Regression Tree Analysis Help Identify Clinically Meaningful Risk Groups for Hip Fracture Prediction in Older American Men (The MrOS Cohort Study)?

Although the WHO fracture risk algorithm (FRAX) is used to predict fracture, the utility of some simple machine‐learning methods, such as classification and regression trees (CARTs) should be evaluated to determine their efficacy in fracture prediction. Follow‐up time for the hip fracture analyses of 5977 community‐dwelling American men aged ≥65 years old was truncated to 10 years. There were 172 (2.9%) men who had an incident nontraumatic hip fracture. The CARTs were developed using hip BMD and common clinical risk factors as follows: model 1 = using classification with continuous variables of age, total hip BMD, and femoral neck BMD, or together with common clinical risk factors; and model 2 = using classification with continuous variables of age, total hip BMD, femoral neck BMD, FRAX score, osteoporosis by T‐score at the hip, and common clinical risk factors. The predictive performance of risk models derived from CARTs was compared with the basic classification of FRAX at 3% (basic model). From model 1, discriminators selected by CART were total hip BMD, age, and femoral neck BMD; no other clinical risk factors were selected. From model 2, discriminators selected by CART were FRAX score, femoral neck BMD, and age. Compared with the basic model using only a high‐risk group by FRAX ≥3%, no significantly improved predictive performance was demonstrated by model 1 or model 2 as identified by CART with the area under the receiver‐operating characteristic curve for each model of 0.714 (95% CI, 0.676 to 0.751) or 0.726 (95% CI, 0.690 to 0.762) versus 0.703 (95% CI, 0.667 to 0.740), respectively. The improved overall net reclassification improvement index was 0.02 (95% CI, –0.04 to 0.08) and 0.05 (95% CI, –0.01 to 0.10), respectively. Although a FRAX category is a good clinical indicator for hip fracture risk, a simple classification by age and BMD may provide an alternative way to estimate a clinical risk level of 3.0%. © 2019 The Authors. JBMR Plus is published by Wiley Periodicals, Inc. on behalf of the American Society for Bone and Mineral Research.

Genetic predisposition to increased serum calcium, bone mineral density, and fracture risk in individuals with normal calcium levels: mendelian randomisation study

OBJECTIVE

To determine if genetically increased serum calcium levels are associated with improved bone mineral density and a reduction in osteoporotic fractures.

DESIGN

Mendelian randomisation study.

SETTING

Cohorts used included: the UK Biobank cohort, providing genotypic and estimated bone mineral density data; 25 cohorts from UK, USA, Europe, and China, providing genotypic and fracture data; and 17 cohorts from Europe, providing genotypic and serum calcium data (summary level statistics).

PARTICIPANTS

A genome-wide association meta-analysis of serum calcium levels in up to 61 079 individuals was used to identify genetic determinants of serum calcium levels. The UK Biobank study was used to assess the association of genetic predisposition to increased serum calcium with estimated bone mineral density derived from heel ultrasound in 426 824 individuals who had, on average, calcium levels in the normal range. A fracture genome-wide association meta-analysis comprising 24 cohorts and the UK Biobank including a total of 76 549 cases and 470 164 controls, who, on average, also had calcium levels in the normal range was then performed.

RESULTS

A standard deviation increase in genetically derived serum calcium (0.13 mmol/L or 0.51 mg/dL) was not associated with increased estimated bone mineral density (0.003 g/cm2, 95% confidence interval -0.059 to 0.066; P=0.92) or a reduced risk of fractures (odds ratio 1.01, 95% confidence interval 0.89 to 1.15; P=0.85) in inverse-variance weighted mendelian randomisation analyses. Sensitivity analyses did not provide evidence of pleiotropic effects.

CONCLUSIONS

Genetic predisposition to increased serum calcium levels in individuals with normal calcium levels is not associated with an increase in estimated bone mineral density and does not provide clinically relevant protection against fracture. Whether such predisposition mimics the effect of short term calcium supplementation is not known. Given that the same genetically derived increase in serum calcium is associated with an increased risk of coronary artery disease, widespread calcium supplementation in the general population could provide more risk than benefit.

Meta-Analysis of Genomewide Association Studies Reveals Genetic Variants for Hip Bone Geometry

Hip geometry is an important predictor of fracture. We performed a meta-analysis of GWAS studies in adults to identify genetic variants that are associated with proximal femur geometry phenotypes. We analyzed four phenotypes: (i) femoral neck length; (ii) neck-shaft angle; (iii) femoral neck width, and (iv) femoral neck section modulus, estimated from DXA scans using algorithms of hip structure analysis. In the Discovery stage, 10 cohort studies were included in the fixed-effect meta-analysis, with up to 18,719 men and women ages 16 to 93 years. Association analyses were performed with ∼2.5 million polymorphisms under an additive model adjusted for age, body mass index, and height. Replication analyses of meta-GWAS significant loci (at adjusted genomewide significance [GWS], threshold p ≤ 2.6 × 10-8 ) were performed in seven additional cohorts in silico. We looked up SNPs associated in our analysis, for association with height, bone mineral density (BMD), and fracture. In meta-analysis (combined Discovery and Replication stages), GWS associations were found at 5p15 (IRX1 and ADAMTS16); 5q35 near FGFR4; at 12p11 (in CCDC91); 11q13 (near LRP5 and PPP6R3 (rs7102273)). Several hip geometry signals overlapped with BMD, including LRP5 (chr. 11). Chr. 11 SNP rs7102273 was associated with any-type fracture (p = 7.5 × 10-5 ). We used bone transcriptome data and discovered several significant eQTLs, including rs7102273 and PPP6R3 expression (p = 0.0007), and rs6556301 (intergenic, chr.5 near FGFR4) and PDLIM7 expression (p = 0.005). In conclusion, we found associations between several genes and hip geometry measures that explained 12% to 22% of heritability at different sites. The results provide a defined set of genes related to biological pathways relevant to BMD and etiology of bone fragility. © 2019 American Society for Bone and Mineral Research.

GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures

Bone area is one measure of bone size that is easily derived from dual-energy X-ray absorptiometry (DXA) scans. In a GWA study of DXA bone area of the hip and lumbar spine (N ≥ 28,954), we find thirteen independent association signals at twelve loci that replicate in samples of European and East Asian descent (N = 13,608 - 21,277). Eight DXA area loci associate with osteoarthritis, including rs143384 in GDF5 and a missense variant in COL11A1 (rs3753841). The strongest DXA area association is with rs11614913[T] in the microRNA MIR196A2 gene that associates with lumbar spine area (P = 2.3 × 10-42, β = -0.090) and confers risk of hip fracture (P = 1.0 × 10-8, OR = 1.11). We demonstrate that the risk allele is less efficient in repressing miR-196a-5p target genes. We also show that the DXA area measure contributes to the risk of hip fracture independent of bone density.

The effects of resistance training on bone mineral density and bone quality in type 2 diabetic rats

Resistance training (RT) has been known to be effective in maintaining and improving bone strength, which is based on bone mineral density (BMD) and bone quality. However, it is not clear whether RT is effective in improving bone strength in patients with type-2 diabetes mellitus (T2DM), who have a high risk of fracture. Therefore, we tested the effects of a 6-week RT regimen using percutaneous electrical stimulation in T2DM model rats, male Otsuka Long-Evans Tokushima Fatty (OLETF), and its control, Long-Evans Tokushima Otsuka (LETO). After 6 weeks of RT, tibial BMD in RT legs was significantly higher than that in control (CON) legs in both groups. In diaphyseal cortical bone, bone area/tissue area, and cortical thickness was significantly increased in RT legs compared with CON legs in both groups. Cortical porosity was highly observed in OLETF compared with LETO, but RT improved cortical porosity in both groups. Interestingly, trabecular number, trabecular thickness and trabecular space as well as BMD and bone volume/tissue volume in proximal tibial metaphyseal trabecular bone were significantly improved in RT legs compared with CON legs in both groups. In contrast, connectivity density and structural model index were not affected by RT. These results indicate that the 6-week RT regimen effectively increased BMD and improved bone quality in T2DM model rats as well as control rats. Therefore, RT may have the potential to improve bone strength and reduce fracture risk, even in patients with T2DM.

Human health in relation to exposure to solar ultraviolet radiation under changing stratospheric ozone and climate

The Montreal Protocol has limited increases in the UV-B (280-315 nm) radiation reaching the Earth's surface as a result of depletion of stratospheric ozone. Nevertheless, the incidence of skin cancers continues to increase in most light-skinned populations, probably due mainly to risky sun exposure behaviour. In locations with strong sun protection programs of long duration, incidence is now reducing in younger age groups. Changes in the epidemiology of UV-induced eye diseases are less clear, due to a lack of data. Exposure to UV radiation plays a role in the development of cataracts, pterygium and possibly age-related macular degeneration; these are major causes of visual impairment world-wide. Photodermatoses and phototoxic reactions to drugs are not uncommon; management of the latter includes recognition of the risks by the prescribing physician. Exposure to UV radiation has benefits for health through the production of vitamin D in the skin and modulation of immune function. The latter has benefits for skin diseases such as psoriasis and possibly for systemic autoimmune diseases such as multiple sclerosis. The health risks of sun exposure can be mitigated through appropriate sun protection, such as clothing with both good UV-blocking characteristics and adequate skin coverage, sunglasses, shade, and sunscreen. New sunscreen preparations provide protection against a broader spectrum of solar radiation, but it is not clear that this has benefits for health. Gaps in knowledge make it difficult to derive evidence-based sun protection advice that balances the risks and benefits of sun exposure.

Role of Epigenomics in Bone and Cartilage Disease

Phenotypic variation in skeletal traits and diseases is the product of genetic and environmental factors. Epigenetic mechanisms include information-containing factors, other than DNA sequence, that cause stable changes in gene expression and are maintained during cell divisions. They represent a link between environmental influences, genome features, and the resulting phenotype. The main epigenetic factors are DNA methylation, posttranslational changes of histones, and higher-order chromatin structure. Sometimes non-coding RNAs, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), are also included in the broad term of epigenetic factors. There is rapidly expanding experimental evidence for a role of epigenetic factors in the differentiation of bone cells and the pathogenesis of skeletal disorders, such as osteoporosis and osteoarthritis. However, different from genetic factors, epigenetic signatures are cell- and tissue-specific and can change with time. Thus, elucidating their role has particular difficulties, especially in human studies. Nevertheless, epigenomewide association studies are beginning to disclose some disease-specific patterns that help to understand skeletal cell biology and may lead to development of new epigenetic-based biomarkers, as well as new drug targets useful for treating diffuse and localized disorders. Here we provide an overview and update of recent advances on the role of epigenomics in bone and cartilage diseases. © 2019 American Society for Bone and Mineral Research.

Vitamin D testing and treatment: a narrative review of current evidence

Vitamin D testing and treatment is a subject of controversial scientific discussions, and it is challenging to navigate through the expanding vitamin D literature with heterogeneous and partially opposed opinions and recommendations. In this narrative review, we aim to provide an update on vitamin D guidelines and the current evidence on the role of vitamin D for human health with its subsequent implications for patient care and public health issues. Vitamin D is critical for bone and mineral metabolism, and it is established that vitamin D deficiency can cause rickets and osteomalacia. While many guidelines recommend target serum 25-hydroxyvitamin D (25[OH]D) concentrations of ≥50 nmol/L (20 ng/mL), the minimum consensus in the scientific community is that serum 25(OH)D concentrations below 25-30 nmol/L (10-12 ng/mL) must be prevented and treated. Using this latter threshold of serum 25(OH)D concentrations, it has been documented that there is a high worldwide prevalence of vitamin D deficiency that may require public health actions such as vitamin D food fortification. On the other hand, there is also reason for concern that an exploding rate of vitamin D testing and supplementation increases costs and might potentially be harmful. In the scientific debate on vitamin D, we should consider that nutrient trials differ from drug trials and that apart from the opposed positions regarding indications for vitamin D treatment we still have to better characterize the precise role of vitamin D for human health.

LRP5, Bone Density, and Mechanical Stress: A Case Report and Literature Review

The Wnt-β-catenin pathway receptor, low-density lipoprotein receptor-related protein 5 (LRP5), is a known regulator of bone mineral density. It has been hypothesized that specific human polymorphisms in LRP5 impact bone density, in part, by altering the anabolic response of bone to mechanical loading. Although experiments in animal models support this hypothesis, there is limited evidence that LRP5 polymorphisms can alter the anabolic response of bone to mechanical loading in humans. Herein, we report a young male who harbors a rare LRP5 missense mutation (A745V) and who provides potential proof of principle for this mechanotransduction hypothesis for low bone density. The subject had no history of fractures until age 18, a year into a career in competitive distance running. As he continued to run over the following 2 years, his mileage threshold to fracture steadily and rapidly decreased until he was diagnosed with severe osteoporosis (lumbar spine BMD Z-score of -3.2). By contextualizing this case within the existing LRP5 and mechanical stress literature, we speculate that this represents the first documented case of an individual in whom a genetic mutation altered the anabolic response of bone to mechanical stress in a manner sufficient to contribute to osteoporosis.

Recent Advances in the Genetics of Fractures in Osteoporosis

Genetic susceptibility, together with old age, female sex, and low bone mineral density (BMD) are amongst the strongest determinants of fracture risk. Tmost recent large-scale genome-wide association study (GWAS) meta-analysis has yielded fifteen loci. This review focuses on the advances in the research of genetic determinants of fracture risk. We first discuss the genetic architecture of fracture risk, touching upon different methods and overall findings. We then discuss in a second paragraph the most recent advances in the field and focus on the genetics of fracture risk and also of other endophenotypes closely related to fracture risk such as bone mineral density (BMD). Application of state-of-the-art methodology such as Mendelian randzation in fracture GWAS are reviewed. The final part of this review touches upon potential future directions in genetic research of osteoporotic fractures.

Mouse Models and Online Resources for Functional Analysis of Osteoporosis Genome-Wide Association Studies

Osteoporosis is a complex genetic disease in which the number of loci associated with the bone mineral density, a clinical risk factor for fracture, has increased at an exponential rate in the last decade. The identification of the causative variants and candidate genes underlying these loci has not been able to keep pace with the rate of locus discovery. A large number of tools and data resources have been built around the use of the mouse as model of human genetic disease. Herein, we describe resources available for functional validation of human Genome Wide Association Study (GWAS) loci using mouse models. We specifically focus on large-scale phenotyping efforts focused on bone relevant phenotypes and repositories of genotype-phenotype data that exist for transgenic and mutant mice, which can be readily mined as a first step toward more targeted efforts designed to deeply characterize the role of a gene in bone biology.

Using Mendelian Randomization to Decipher Mechanisms of Bone Disease

PURPOSE OF REVIEW

This review summarizes the basic principles of Mendelian randomization (MR) and provides evidence for the causal effect of multiple modifiable factors on bone outcomes.

RECENT FINDINGS

Several studies using MR approach have provided support for the causal effect of obesity on bone mineral density (BMD). Strikingly, studies have failed to prove a causal association between elevated 25(OH) D concentrations and higher BMD in community-dwelling individuals. The MR approach has been successfully used to evaluate multiple factors related to bone mineral density variation and/or fracture risk. The MR approach avoids some of the classical observational study limitations and provides more robust causal evidence, ensuring bigger success of the clinical trials. The selection of interventions based on genetic evidence could have a substantial impact on clinical practice.