Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations.

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Genetic underpinning of the comorbidity between type 2 diabetes and osteoarthritis.

Multimorbidity is a rising public health challenge with important implications for health management and policy. The most common multimorbidity pattern is the combination of cardiometabolic and osteoarticular diseases. Here, we study the genetic underpinning of the comorbidity between type 2 diabetes and osteoarthritis. We find genome-wide genetic correlation between the two diseases and robust evidence for association-signal colocalization at 18 genomic regions. We integrate multi-omics and functional information to resolve the colocalizing signals and identify high-confidence effector genes, including FTO and IRX3, which provide proof-of-concept insights into the epidemiologic link between obesity and both diseases. We find enrichment for lipid metabolism and skeletal formation pathways for signals underpinning the knee and hip osteoarthritis comorbidities with type 2 diabetes, respectively. Causal inference analysis identifies complex effects of tissue-specific gene expression on comorbidity outcomes. Our findings provide insights into the biological basis for the type 2 diabetes-osteoarthritis disease co-occurrence.

Strategies to investigate and mitigate collider bias in genetic and Mendelian randomisation studies of disease progression.

Genetic studies of disease progression can be used to identify factors that may influence survival or prognosis, which may differ from factors that influence on disease susceptibility. Studies of disease progression feed directly into therapeutics for disease, whereas studies of incidence inform prevention strategies. However, studies of disease progression are known to be affected by collider (also known as "index event") bias since the disease progression phenotype can only be observed for individuals who have the disease. This applies equally to observational and genetic studies, including genome-wide association studies and Mendelian randomisation (MR) analyses. In this paper, our aim is to review several statistical methods that can be used to detect and adjust for index event bias in studies of disease progression, and how they apply to genetic and MR studies using both individual- and summary-level data. Methods to detect the presence of index event bias include the use of negative controls, a comparison of associations between risk factors for incidence in individuals with and without the disease, and an inspection of Miami plots. Methods to adjust for the bias include inverse probability weighting (with individual-level data), or Slope-Hunter and Dudbridge et al.'s index event bias adjustment (when only summary-level data are available). We also outline two approaches for sensitivity analysis. We then illustrate how three methods to minimise bias can be used in practice with two applied examples. Our first example investigates the effects of blood lipid traits on mortality from coronary heart disease, while our second example investigates genetic associations with breast cancer mortality.

Causal relationships between anthropometric traits, bone mineral density, osteoarthritis and spinal stenosis: A Mendelian randomisation investigation.

OBJECTIVE: Spinal stenosis is a common condition among older individuals, with significant morbidity attached. Little is known about its risk factors but degenerative conditions, such as osteoarthritis (OA) have been identified for their mechanistic role. This study aims to explore causal relationships between anthropometric risk factors, OA, and spinal stenosis using Mendelian randomisation (MR) techniques. DESIGN: We applied two-sample MR to investigate the causal relationships between genetic liability for select risk factors and spinal stenosis. Next, we examined the genetic relationship between OA and spinal stenosis with linkage disequilibrium score regression and Causal Analysis Using Summary Effect estimates MR method. Finally, we used multivariable MR (MVMR) to explore whether OA and body mass index (BMI) mediate the causal pathways identified. RESULTS: Our analysis revealed strong evidence for the effect of higher BMI (odds ratio [OR] = 1.54, 95%CI: 1.41-1.69, p-value = 2.7 × 10-21), waist (OR = 1.43, 95%CI: 1.15-1.79, p-value = 1.5 × 10-3) and hip (OR = 1.50, 95%CI: 1.27-1.78, p-value = 3.3 × 10-6) circumference on spinal stenosis. Strong evidence of causality was also observed for higher bone mineral density (BMD): total body (OR = 1.21, 95%CI: 1.12-1.29, p-value = 1.6 × 10-7), femoral neck (OR = 1.35, 95%CI: 1.09-1.37, p-value = 7.5×10-7), and lumbar spine (OR = 1.38, 95%CI: 1.25-1.52, p-value = 4.4 × 10-11). We detected high genetic correlations between spinal stenosis and OA (rg range: 0.47-0.66), with Causal Analysis Using Summary Effect estimates results supporting a causal effect of OA on spinal stenosis (ORallOA = 1.6, 95%CI: 1.41-1.79). Direct effects of BMI, BMD on spinal stenosis remained after adjusting for OA in the MVMR. CONCLUSIONS: Genetic susceptibility to anthropometric risk factors, particularly higher BMI and BMD can increase the risk of spinal stenosis, independent of OA status. These results may inform preventative strategies and treatments.

Adjusting for collider bias in genetic association studies using instrumental variable methods.

Genome-wide association studies have provided many genetic markers that can be used as instrumental variables to adjust for confounding in epidemiological studies. Recently, the principle has been applied to other forms of bias in observational studies, especially collider bias that arises when conditioning or stratifying on a variable that is associated with the outcome of interest. An important case is in studies of disease progression and survival. Here, we clarify the links between the genetic instrumental variable methods proposed for this problem and the established methods of Mendelian randomisation developed to account for confounding. We highlight the critical importance of weak instrument bias in this context and describe a corrected weighted least-squares procedure as a simple approach to reduce this bias. We illustrate the range of available methods on two data examples. The first, waist-hip ratio adjusted for body-mass index, entails statistical adjustment for a quantitative trait. The second, smoking cessation, is a stratified analysis conditional on having initiated smoking. In both cases, we find little effect of collider bias on the primary association results, but this may propagate into more substantial effects on further analyses such as polygenic risk scoring and Mendelian randomisation.

Age-specific effects of weight-based body size on fracture risk in later life: a lifecourse Mendelian randomisation study.

Musculoskeletal conditions, including fractures, can have severe and long-lasting consequences. Higher body mass index in adulthood is widely acknowledged to be protective for most fracture sites. However, sources of bias induced by confounding factors may have distorted previous findings. Employing a lifecourse Mendelian randomisation (MR) approach by using genetic instruments to separate effects at different life stages, this investigation aims to explore how prepubertal and adult body size independently influence fracture risk in later life.Using data from a large prospective cohort, univariable and multivariable MR were conducted to simultaneously estimate the effects of age-specific genetic proxies for body size (n = 453,169) on fracture risk (n = 416,795). A two-step MR framework was additionally applied to elucidate potential mediators. Univariable and multivariable MR indicated strong evidence that higher body size in childhood reduced fracture risk (OR, 95% CI: 0.89, 0.82 to 0.96, P = 0.005 and 0.76, 0.69 to 0.85, P = 1 × 10- 6, respectively). Conversely, higher body size in adulthood increased fracture risk (OR, 95% CI: 1.08, 1.01 to 1.16, P = 0.023 and 1.26, 1.14 to 1.38, P = 2 × 10- 6, respectively). Two-step MR analyses suggested that the effect of higher body size in childhood on reduced fracture risk was mediated by its influence on higher estimated bone mineral density (eBMD) in adulthood.This investigation provides novel evidence that higher body size in childhood reduces fracture risk in later life through its influence on increased eBMD. From a public health perspective, this relationship is complex since obesity in adulthood remains a major risk factor for co-morbidities. Results additionally indicate that higher body size in adulthood is a risk factor for fractures. Protective effect estimates previously observed are likely attributed to childhood effects.

Mendelian randomization provides evidence for a causal effect of higher serum IGF-1 concentration on risk of hip and knee osteoarthritis.

OBJECTIVES: How insulin-like growth factor-1 (IGF-1) is related to OA is not well understood. We determined relationships between IGF-1 and hospital-diagnosed hand, hip and knee OA in UK Biobank, using Mendelian randomization (MR) to determine causality. METHODS: Serum IGF-1 was assessed by chemiluminescent immunoassay. OA was determined using Hospital Episode Statistics. One-sample MR (1SMR) was performed using two-stage least-squares regression, with an unweighted IGF-1 genetic risk score as an instrument. Multivariable MR included BMI as an additional exposure (instrumented by BMI genetic risk score). MR analyses were adjusted for sex, genotyping chip and principal components. We then performed two-sample MR (2SMR) using summary statistics from Cohorts for Heart and Aging Research in Genetic Epidemiology (CHARGE) (IGF-1, N = 30 884) and the recent genome-wide association study meta-analysis (N = 455 221) of UK Biobank and Arthritis Research UK OA Genetics (arcOGEN). RESULTS: A total of 332 092 adults in UK Biobank had complete data. Their mean (s.d.) age was 56.5 (8.0) years and 54% were female. IGF-1 was observationally related to a reduced odds of hand OA [odds ratio per doubling = 0.87 (95% CI 0.82, 0.93)], and an increased odds of hip OA [1.04 (1.01, 1.07)], but was unrelated to knee OA [0.99 (0.96, 1.01)]. Using 1SMR, we found strong evidence for an increased risk of hip [odds ratio per s.d. increase = 1.57 (1.21, 2.01)] and knee [1.30 (1.07, 1.58)] OA with increasing IGF-1 concentration. By contrast, we found no evidence for a causal effect of IGF-1 concentration on hand OA [0.98 (0.57, 1.70)]. Results were consistent when estimated using 2SMR and in multivariable MR analyses accounting for BMI. CONCLUSION: We have found evidence that increased serum IGF-1 is causally related to higher risk of hip and knee OA.

Metabolomics analysis in adults with high bone mass identifies a relationship between bone resorption and circulating citrate which replicates in the general population.

OBJECTIVE: Bone turnover, which regulates bone mass, may exert metabolic consequences, particularly on markers of glucose metabolism and adiposity. To better understand these relationships, we examined cross-sectional associations between bone turnover markers (BTMs) and metabolic traits in a population with high bone mass (HBM, BMD Z-score ≥+3.2). DESIGN: ß-C-terminal telopeptide of type-I collagen (ß-CTX), procollagen type-1 amino-terminal propeptide (P1NP) and osteocalcin were assessed by electrochemiluminescence immunoassays. Metabolic traits, including lipids and glycolysis-related metabolites, were measured using nuclear magnetic resonance spectroscopy. Associations of BTMs with metabolic traits were assessed using generalized estimating equation linear regression, accounting for within-family correlation, adjusting for potential confounders (age, sex, height, weight, menopause, bisphosphonate and oral glucocorticoid use). RESULTS: A total of 198 adults with HBM had complete data, mean [SD] age 61.6 [13.7] years; 77% were female. Of 23 summary metabolic traits, citrate was positively related to all BTMs: adjusted ßß-CTX  = 0.050 (95% CI 0.024, 0.076), P = 1.71 × 10-4 , ßosteocalcin  = 6.54 × 10-4 (1.87 × 10-4 , 0.001), P = .006 and ßP1NP  = 2.40 × 10-4 (6.49 × 10-5 , 4.14 × 10-4 ), P = .007 (ß = increase in citrate (mmol/L) per 1 µg/L BTM increase). Inverse relationships of ß-CTX (ß = -0.276 [-0.434, -0.118], P = 6.03 × 10-4 ) and osteocalcin (-0.004 [-0.007, -0.001], P = .020) with triglycerides were also identified. We explored the generalizability of these associations in 3664 perimenopausal women (age 47.9 [4.4] years) from a UK family cohort. We confirmed a positive, albeit lower magnitude, association between ß-CTX and citrate (adjusted ßwomen  = 0.020 [0.013, 0.026], P = 1.95 × 10-9 ) and an inverse association of similar magnitude between ß-CTX and triglycerides (ß = -0.354 [-0.471, -0.237], P = 3.03 × 10-9 ). CONCLUSIONS: Bone resorption is positively related to circulating citrate and inversely related to triglycerides. Further studies are justified to determine whether plasma citrate or triglyceride concentrations are altered by factors known to modulate bone resorption, such as bisphosphonates.

Correlates of high-impact physical activity measured objectively in older British adults.

Background: Exposure to higher magnitude vertical impacts is thought to benefit bone health. The correlates of this high-impact physical activity (PA) in later life are unknown. Methods: Participants were from the Cohort for Skeletal Health in Bristol and Avon, Hertfordshire Cohort Study and MRC National Survey of Health and Development. Associations of demographic, behavioural, physiological and psychological factors with vertical acceleration peaks ≥1.5 g (i.e. high-impact PA) from 7-day hip-worn accelerometer recordings were examined using linear regression. Results: A total of 1187 participants (mean age = 72.7 years, 66.6% females) were included. Age, sex, education, active transport, self-reported higher impact PA, walking speed and self-rated health were independently associated with high-impact PA whereas BMI and sleep quality showed borderline independent associations. For example, differences in log-high-impact counts were 0.50 (P < 0.001) for men versus women and -0.56 (P < 0.001) for worst versus best self-rated health. Our final model explained 23% of between-participant variance in high impacts. Other correlates were not associated with high-impact activity after adjustment. Conclusions: Besides age and sex, several factors were associated with higher impact PA in later life. Our findings help identify characteristics of older people that might benefit from interventions designed to promote osteogenic PA.

Characterization of Vertical Accelerations Experienced by Older People Attending an Aerobics Class Designed to Produce High Impacts.

The purpose of this study was to establish the feasibility of using an aerobics class to produce potentially bone protective vertical impacts of ≥ 4g in older adults and to determine whether impacts can be predicted by physical function. Participants recruited from older adult exercise classes completed an SF-12 questionnaire, short physical performance battery, and an aerobics class with seven different components, performed at low and high intensity. Maximum g and jerk values were identified for each activity. Forty-one participants (mean 69 years) were included. Mean maximal values approached or exceeded the 4g threshold for four of the seven exercises. In multivariate analyses, age (-0.53; -0.77, -0.28) (standardized beta coefficient; 95% CI) and 4-m walk time (-0.39; -0.63, -0.16) were inversely related to maximum g. Aerobics classes can be used to produce relatively high vertical accelerations in older individuals, although the outcome is strongly dependent on age and physical function.

Protein identification for stroke progression via Mendelian Randomization in Million Veteran Program and UK Biobank.

Background: Individuals who have experienced a stroke, or transient ischemic attack, face a heightened risk of future cardiovascular events. Identification of genetic and molecular risk factors for subsequent cardiovascular outcomes may identify effective therapeutic targets to improve prognosis after an incident stroke. Methods: We performed genome-wide association studies (GWAS) for subsequent major adverse cardiovascular events (MACE) (Ncases=51,929, Ncntrl=39,980) and subsequent arterial ischemic stroke (AIS) Ncases=45,120, Ncntrl=46,789) after first incident stroke within the Million Veteran Program and UK Biobank. We then used genetic variants associated with proteins (pQTLs) to determine the effect of 1,463 plasma protein abundances on subsequent MACE using Mendelian randomization (MR). Results: Two variants were significantly associated with subsequent cardiovascular events: rs76472767 (OR=0.75, 95% CI = 0.64-0.85, p= 3.69×10-08) with subsequent AIS and rs13294166 (OR=1.52, 95% CI = 1.37-1.67, p=3.77×10-08) with subsequent MACE. Using MR, we identified 2 proteins with an effect on subsequent MACE after a stroke: CCL27 (effect OR= 0.77, 95% CI = 0.66-0.88, adj. p=0.05), and TNFRSF14 (effect OR=1.42, 95% CI = 1.24-1.60, adj. p=0.006). These proteins are not associated with incident AIS and are implicated to have a role in inflammation. Conclusions: We found evidence that two proteins with little effect on incident stroke appear to influence subsequent MACE after incident AIS. These associations suggest that inflammation is a contributing factor to subsequent MACE outcomes after incident AIS and highlights potential novel targets.

The identification of distinct protective and susceptibility mechanisms for hip osteoarthritis: findings from a genome-wide association study meta-analysis of minimum joint space width and Mendelian randomisation cluster analyses.

BACKGROUND: Hip minimum joint space width (mJSW) provides a proxy for cartilage thickness. This study aimed to conduct a genome-wide association study (GWAS) of mJSW to (i) identify new genetic determinants of mJSW and (ii) identify which mJSW loci convey hip osteoarthritis (HOA) risk and would therefore be of therapeutic interest. METHODS: GWAS meta-analysis of hip mJSW derived from plain X-rays and DXA was performed, stratified by sex and adjusted for age and ancestry principal components. Mendelian randomisation (MR) and cluster analyses were used to examine causal effect of mJSW on HOA. FINDINGS: 50,745 individuals were included in the meta-analysis. 42 SNPs, which mapped to 39 loci, were identified. Mendelian randomisation (MR) revealed little evidence of a causal effect of mJSW on HOA (ORIVW 0.98 [95% CI 0.82-1.18]). However, MR-Clust analysis suggested the null MR estimates reflected the net effect of two distinct causal mechanisms cancelling each other out, one of which was protective, whereas the other increased HOA susceptibility. For the latter mechanism, all loci were positively associated with height, suggesting mechanisms leading to greater height and mJSW increase the risk of HOA in later life. INTERPRETATIONS: One group of mJSW loci reduce HOA risk via increased mJSW, suggesting possible utility as targets for chondroprotective therapies. The second group of mJSW loci increased HOA risk, despite increasing mJSW, but were also positively related to height, suggesting they contribute to HOA risk via a growth-related mechanism. FUNDING: Primarily funded by the Medical Research Council and Wellcome Trust.

A Genome-Wide Association Study Meta-Analysis of Alpha Angle Suggests Cam-Type Morphology May Be a Specific Feature of Hip Osteoarthritis in Older Adults.

OBJECTIVE: To examine the genetic architecture of cam morphology using alpha angle (AA) as a proxy measure and conduct an AA genome-wide association study (GWAS) followed by Mendelian randomization (MR) to evaluate its causal relationship with hip osteoarthritis (OA). METHODS: Observational analyses examined associations between AA measurements derived from hip dual x-ray absorptiometry (DXA) scans from the UK Biobank study and radiographic hip OA outcomes and subsequent total hip replacement. Following these analyses, an AA GWAS meta-analysis was performed (N = 44,214) using AA measurements previously derived in the Rotterdam Study. Linkage disequilibrium score regression assessed the genetic correlation between AA and hip OA. Genetic associations considered significant (P < 5 × 10-8 ) were used as AA genetic instrument for 2-sample MR analysis. RESULTS: DXA-derived AA showed expected associations between AA and radiographic hip OA (adjusted odds ratio [OR] 1.63 [95% confidence interval (95% CI) 1.58, 1.67]) and between AA and total hip replacement (adjusted hazard ratio 1.45 [95% CI 1.33, 1.59]) in the UK Biobank study cohort. The heritability of AA was 10%, and AA had a moderate genetic correlation with hip OA (rg  = 0.26 [95% CI 0.10, 0.43]). Eight independent genetic signals were associated with AA. Two-sample MR provided weak evidence of causal effects of AA on hip OA risk (inverse variance weighted OR 1.84 [95% CI 1.14, 2.96], P = 0.01). In contrast, genetic predisposition for hip OA had stronger evidence of a causal effect on increased AA (inverse variance weighted ß = 0.09 [95% CI 0.04, 0.13], P = 4.58 × 10-5 ). CONCLUSION: Expected observational associations between AA and related clinical outcomes provided face validity for the DXA-derived AA measurements. Evidence of bidirectional associations between AA and hip OA, particularly for risk of hip OA on AA, suggests that hip shape modeling secondary to a genetic predisposition to hip OA contributes to the well-established relationship between hip OA and cam morphology in older adults.

Lowering of Circulating Sclerostin May Increase Risk of Atherosclerosis and Its Risk Factors: Evidence From a Genome-Wide Association Meta-Analysis Followed by Mendelian Randomization.

OBJECTIVE: In this study, we aimed to establish the causal effects of lowering sclerostin, target of the antiosteoporosis drug romosozumab, on atherosclerosis and its risk factors. METHODS: A genome-wide association study meta-analysis was performed of circulating sclerostin levels in 33,961 European individuals. Mendelian randomization (MR) was used to predict the causal effects of sclerostin lowering on 15 atherosclerosis-related diseases and risk factors. RESULTS: We found that 18 conditionally independent variants were associated with circulating sclerostin. Of these, 1 cis signal in SOST and 3 trans signals in B4GALNT3, RIN3, and SERPINA1 regions showed directionally opposite signals for sclerostin levels and estimated bone mineral density. Variants with these 4 regions were selected as genetic instruments. MR using 5 correlated cis-SNPs suggested that lower sclerostin increased the risk of type 2 diabetes mellitus (DM) (odds ratio [OR] 1.32 [95% confidence interval (95% CI) 1.03-1.69]) and myocardial infarction (MI) (OR 1.35 [95% CI 1.01-1.79]); sclerostin lowering was also suggested to increase the extent of coronary artery calcification (CAC) (ß = 0.24 [95% CI 0.02-0.45]). MR using both cis and trans instruments suggested that lower sclerostin increased hypertension risk (OR 1.09 [95% CI 1.04-1.15]), but otherwise had attenuated effects. CONCLUSION: This study provides genetic evidence to suggest that lower levels of sclerostin may increase the risk of hypertension, type 2 DM, MI, and the extent of CAC. Taken together, these findings underscore the requirement for strategies to mitigate potential adverse effects of romosozumab treatment on atherosclerosis and its related risk factors.

A Guide for Understanding and Designing Mendelian Randomization Studies in the Musculoskeletal Field.

Mendelian randomization (MR) is an increasingly popular component of an epidemiologist's toolkit, used to provide evidence of a causal effect of one trait (an exposure, eg, body mass index [BMI]) on an outcome trait or disease (eg, osteoarthritis). Identifying these effects is important for understanding disease etiology and potentially identifying targets for therapeutic intervention. MR uses genetic variants as instrumental variables for the exposure, which should not be influenced by the outcome or confounding variables, overcoming key limitations of traditional epidemiological analyses. For MR to generate a valid estimate of effect, key assumptions must be met. In recent years, there has been a rapid rise in MR methods that aim to test, or are robust to violations of, these assumptions. In this review, we provide an overview of MR for a non-expert audience, including an explanation of these key assumptions and how they are often tested, to aid a better reading and understanding of the MR literature. We highlight some of these new methods and how they can be useful for specific methodological challenges in the musculoskeletal field, including for conditions or traits that share underlying biological pathways, such as bone and joint disease. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Using multivariable Mendelian randomization to estimate the causal effect of bone mineral density on osteoarthritis risk, independently of body mass index.

OBJECTIVES: Observational analyses suggest that high bone mineral density (BMD) is a risk factor for osteoarthritis (OA); it is unclear whether this represents a causal effect or shared aetiology and whether these relationships are body mass index (BMI)-independent. We performed bidirectional Mendelian randomization (MR) to uncover the causal pathways between BMD, BMI and OA. METHODS: One-sample (1S)MR estimates were generated by two-stage least-squares regression. Unweighted allele scores instrumented each exposure. Two-sample (2S)MR estimates were generated using inverse-variance weighted random-effects meta-analysis. Multivariable MR (MVMR), including BMD and BMI instruments in the same model, determined the BMI-independent causal pathway from BMD to OA. Latent causal variable (LCV) analysis, using weight-adjusted femoral neck (FN)-BMD and hip/knee OA summary statistics, determined whether genetic correlation explained the causal effect of BMD on OA. RESULTS: 1SMR provided strong evidence for a causal effect of BMD estimated from heel ultrasound (eBMD) on hip and knee OA {odds ratio [OR]hip = 1.28 [95% confidence interval (CI) = 1.05, 1.57], p = 0.02, ORknee = 1.40 [95% CI = 1.20, 1.63], p = 3 × 10-5, OR per standard deviation [SD] increase}. 2SMR effect sizes were consistent in direction. Results suggested that the causal pathways between eBMD and OA were bidirectional (ßhip = 1.10 [95% CI = 0.36, 1.84], p = 0.003, ßknee = 4.16 [95% CI = 2.74, 5.57], p = 8 × 10-9, ß = SD increase per doubling in risk). MVMR identified a BMI-independent causal pathway between eBMD and hip/knee OA. LCV suggested that genetic correlation (i.e. shared genetic aetiology) did not fully explain the causal effects of BMD on hip/knee OA. CONCLUSIONS: These results provide evidence for a BMI-independent causal effect of eBMD on OA. Despite evidence of bidirectional effects, the effect of BMD on OA did not appear to be fully explained by shared genetic aetiology, suggesting a direct action of bone on joint deterioration.

Increased development of radiographic hip osteoarthritis in individuals with high bone mass: a prospective cohort study.

BACKGROUND: Individuals with high bone mass (HBM) have a greater odds of prevalent radiographic hip osteoarthritis (OA), reflecting an association with bone-forming OA sub-phenotypes (e.g. osteophytosis, subchondral sclerosis). As the role of bone mineral density (BMD) in hip OA progression is unclear, we aimed to determine if individuals with HBM have increased incidence and/or progression of bone-forming OA sub-phenotypes. METHODS: We analysed an adult cohort with and without HBM (L1 and/or total hip BMD Z-score > + 3.2) with pelvic radiographs collected at baseline and 8-year follow-up. Sub-phenotypes were graded using the OARSI atlas. Superior/inferior acetabular/femoral osteophyte and medial/superior joint space narrowing (JSN) grades were summed and Δosteophyte and ΔJSN derived. Pain and functional limitations were quantified using the WOMAC questionnaire. Associations between HBM status and change in OA sub-phenotypes were determined using multivariable linear/logistic regression, adjusting for age, sex, height, total body fat mass, follow-up time and baseline sub-phenotype grade. Generalised estimating equations accounted for individual-level clustering. RESULTS: Of 136 individuals, 62% had HBM at baseline, 72% were female and mean (SD) age was 59 (10) years. HBM was positively associated with both Δosteophytes and ΔJSN (adjusted mean grade differences between individuals with and without HBM ßosteophyte = 0.30 [0.01, 0.58], p = 0.019 and ßJSN = 0.10 [0.01, 0.18], p = 0.019). Incident subchondral sclerosis was rare. HBM individuals had higher WOMAC hip functional limitation scores (ß = 8.3 [0.7, 15.98], p = 0.032). CONCLUSIONS: HBM is associated with the worsening of hip osteophytes and JSN over an average of 8 years, as well as increased hip pain and functional limitation.

Genome-wide association study implicates novel loci and reveals candidate effector genes for longitudinal pediatric bone accrual.

BACKGROUND: Bone accrual impacts lifelong skeletal health, but genetic discovery has been primarily limited to cross-sectional study designs and hampered by uncertainty about target effector genes. Here, we capture this dynamic phenotype by modeling longitudinal bone accrual across 11,000 bone scans in a cohort of healthy children and adolescents, followed by genome-wide association studies (GWAS) and variant-to-gene mapping with functional follow-up. RESULTS: We identify 40 loci, 35 not previously reported, with various degrees of supportive evidence, half residing in topological associated domains harboring known bone genes. Of several loci potentially associated with later-life fracture risk, a candidate SNP lookup provides the most compelling evidence for rs11195210 (SMC3). Variant-to-gene mapping combining ATAC-seq to assay open chromatin with high-resolution promoter-focused Capture C identifies contacts between GWAS loci and nearby gene promoters. siRNA knockdown of gene expression supports the putative effector gene at three specific loci in two osteoblast cell models. Finally, using CRISPR-Cas9 genome editing, we confirm that the immediate genomic region harboring the putative causal SNP influences PRPF38A expression, a location which is predicted to coincide with a set of binding sites for relevant transcription factors. CONCLUSIONS: Using a new longitudinal approach, we expand the number of genetic loci putatively associated with pediatric bone gain. Functional follow-up in appropriate cell models finds novel candidate genes impacting bone accrual. Our data also raise the possibility that the cell fate decision between osteogenic and adipogenic lineages is important in normal bone accrual.

A Rare Mutation in SMAD9 Associated With High Bone Mass Identifies the SMAD-Dependent BMP Signaling Pathway as a Potential Anabolic Target for Osteoporosis.

Novel anabolic drug targets are needed to treat osteoporosis. Having established a large national cohort with unexplained high bone mass (HBM), we aimed to identify a novel monogenic cause of HBM and provide insight into a regulatory pathway potentially amenable to therapeutic intervention. We investigated a pedigree with unexplained HBM in whom previous sequencing had excluded known causes of monogenic HBM. Whole exome sequencing identified a rare (minor allele frequency 0.0023), highly evolutionarily conserved missense mutation in SMAD9 (c.65T>C, p.Leu22Pro) segregating with HBM in this autosomal dominant family. The same mutation was identified in another two unrelated individuals both with HBM. In silico protein modeling predicts the mutation severely disrupts the MH1 DNA-binding domain of SMAD9. Affected individuals have bone mineral density (BMD) Z-scores +3 to +5, mandible enlargement, a broad frame, torus palatinus/mandibularis, pes planus, increased shoe size, and a tendency to sink when swimming. Peripheral quantitative computed tomography (pQCT) measurement demonstrates increased trabecular volumetric BMD and increased cortical thickness conferring greater predicted bone strength; bone turnover markers are low/normal. Notably, fractures and nerve compression are not found. Both genome-wide and gene-based association testing involving estimated BMD measured at the heel in 362,924 white British subjects from the UK Biobank Study showed strong associations with SMAD9 (PGWAS = 6 × 10-16 ; PGENE = 8 × 10-17 ). Furthermore, we found Smad9 to be highly expressed in both murine cortical bone-derived osteocytes and skeletal elements of zebrafish larvae. Our findings support SMAD9 as a novel HBM gene and a potential novel osteoanabolic target for osteoporosis therapeutics. SMAD9 is thought to inhibit bone morphogenetic protein (BMP)-dependent target gene transcription to reduce osteoblast activity. Thus, we hypothesize SMAD9 c.65T>C is a loss-of-function mutation reducing BMP inhibition. Lowering SMAD9 as a potential novel anabolic mechanism for osteoporosis therapeutics warrants further investigation. © 2019 The Authors. Journal of Bone and Mineral Research published by American Society for Bone and Mineral Research.