Osteoarthritis: Insights Offered by the Study of Bone Mass Genetics.

PURPOSE OF REVIEW: This paper reviews how bone genetics has contributed to our understanding of the pathogenesis of osteoarthritis. As well as identifying specific genetic mechanisms involved in osteoporosis which also contribute to osteoarthritis, we review whether bone mineral density (BMD) plays a causal role in OA development. RECENT
FINDINGS: We examined whether those genetically predisposed to elevated BMD are at increased risk of developing OA, using our high bone mass (HBM) cohort. HBM individuals were found to have a greater prevalence of OA compared with family controls and greater development of radiographic features of OA over 8 years, with predominantly osteophytic OA. Initial Mendelian randomisation analysis provided additional support for a causal effect of increased BMD on increased OA risk. In contrast, more recent investigation estimates this relationship to be bi-directional. However, both these findings could be explained instead by shared biological pathways. Pathways which contribute to BMD appear to play an important role in OA development, likely reflecting shared common mechanisms as opposed to a causal effect of raised BMD on OA. Studies in HBM individuals suggest this reflects an important role of mechanisms involved in bone formation in OA development; however further work is required to establish whether the same applies to more common forms of OA within the general population.

Introduction

Osteoarthritis (OA) and osteoporosis are the two most prevalent musculoskeletal disorders globally, responsible for high levels of disability and healthcare utilisation, for example, through joint replacements and fractures, respectively [1]. OA is recognised as a disease of the whole joint and is considered a disorder of the bone-cartilage unit [2]. Cross-sectional studies consistently find that osteoporosis is inversely related to osteoarthritis (OA), as evidenced by positive associations between bone mineral density (BMD) and prevalent radiographic hip [3-7] and knee [8, 9] OA. Conversely, an inverse relationship has been observed between BMD and progression of knee OA [10, 11]. This may reflect a role of increased bone turnover, which is usually inversely related to BMD and is also associated with OA progression. For example, in the UK Chingford study, Bettica et al. observed that levels of the bone resorption marker, urinary NTX, were higher in progressive compared to that in non-progressive cases [12]. The suggestion that bisphosphonates, which suppress bone turnover and increase BMD, ameliorate some aspects of OA symptoms supports the suggestion that increased bone turnover contributes to OA pathogenesis [13]. This apparent paradox, whereby increased BMD contributes both positively and negatively to OA pathogenesis, may reflect distinct mechanisms according to the stage of OA development. For example, increased subchondral bone remodelling, involving the subchondral plate and subjacent cancellous bone, may contribute to OA initiation [14]. In contrast, subchondral sclerosis, a feature of relatively advanced OA, results from subsequent thickening and increased density of subchondral bone immediately adjacent to the joint and may contribute to further cartilage loss due to a loss of compliance leading to increased strains on the nearby cartilage [15].

In light of these well-established, albeit complex, relationships between BMD and OA, the study of bone genetics may offer important insights into the pathogenesis of OA. As well as identifying possible shared molecular mechanisms, genetic understanding may help in unravelling causal pathways. For example, studying whether those genetically predisposed to raised BMD are at increased risk of OA provides a stronger basis for making causal inferences compared with observational studies, inherently prone to residual confounding.

High Bone Mass and Osteoarthritis

OA Risk in Individuals with HBM

Over the last decade, we have studied our unique population of individuals with HBM to examine a possible causal role of raised BMD in the development of OA; these HBM individuals have a generalised increased BMD, more than 3 standard deviations above average for their age, which is clinically unexplained [16]. Our studies have been based on the premise that HBM individuals are genetically predisposed to their raised BMD, which manifests in young adulthood. Hence, raised BMD is expected to precede any subsequent development of OA, in keeping with a causal relationship. We initially found that individuals with HBM are at increased risk of developing OA, with a 4.5 times higher odds of partial or total hip replacement compared with that in family controls [17]. HBM individuals also had a greater prevalence of joint replacement at any site compared with general population data from the Health Survey for England [17]. We subsequently examined the risk of prevalent OA in the HBM population based on radiographs obtained from these individuals, using age- and sex-matched general population controls from Chingford and Hertfordshire cohort studies. An increased risk of hip and knee radiographic OA was observed [18, 19]. Body mass index (BMI) is an important confounder when examining BMD and OA relationships, since BMI is positively associated with both. This may well contribute to associations between HBM and OA, given our finding that BMI is significantly increased in HBM individuals [20]. Adjustment for BMI made little difference to the results for hip OA but attenuated the OR for the effect of HBM on knee OA from 2.38 to 1.62, suggesting the increased fat mass of HBM individuals partially contributes to the relationship between HBM and knee OA [20].

When analysing the individual radiographic sub-phenotypes of OA (osteophytes, joint space narrowing [JSN], subchondral sclerosis), an increased odds of osteophytes in HBM individuals was observed at both hip and knee joints; plus an increased odds of subchondral sclerosis was observed at the hip. However, there was no difference in odds of JSN between individuals with and without HBM at either joint [18, 19]. The same pattern of increased odds of osteophytes but not JSN in individuals with HBM compared to unaffected relatives was also observed in the (non-weight-bearing) distal interphalangeal (DIP) and carpometacarpal (CMC) joints of the hand [21]. Recently, we were able to examine the relationship between HBM and subsequent OA incidence/progression. We found that HBM is associated with an increase in knee osteophyte and JSN score over an average of 8 years, as well as knee pain and functional limitation [22]. Similar results have also been observed at the hip (A Hartley et al., unpublished data). Therefore, contrary to findings from conventional observational studies, it would seem that raised BMD contributes to OA pathogenesis across all phases of disease.

There are several potential mechanisms for the increased risk of OA which we observed in HBM individuals. If raised BMD is on the causal pathway for OA (Fig. 1a), this would be consistent with a role of increased subchondral density and stiffness as an initiating event in OA development, first proposed by Radin [23]. However, this explanation is now thought less likely following detailed study of temporo-spatial changes during OA development [14]. On the other hand, genetic factors related to BMD might lead to OA via a separate, pleiotropic pathway. This could be due to a genetic variant operating either through multiple pathways or through a particular pathway having direct effects on both phenotypes. For example, as discussed below, the Wnt signalling pathway regulates both osteoblasts and chondrocytes, leaving the possibility that the same genetic influences on BMD lead to OA via separate effects on the cartilage (Fig. 1b). Alternatively, to the extent that HBM reflects an underlying tendency towards increased osteoblast activity and bone formation, the same tendency might increase the risk of developing OA as a consequence of greater susceptibility to osteophyte formation (Fig. 1c). The possible application of Mendelian randomisation (MR) to distinguish these scenarios is discussed under ‘Future Work’.

Fig. 1

Potential mechanisms for the increased risk of OA which we observed in HBM individuals. (A) Causal: raised BMD in HBM increases the risk of OA as a consequence of a causal effect of BMD on OA—in this scenario, genetic variants from any pathway that influences BMD would also be associated with OA. (B) Pleiotropy 1: pathways are affected which influence both osteoblast function and chondrocyte function, leading to an increase in BMD and OA risk, respectively—in this scenario, we might expect only genetic variants from specific BMD-influencing pathways to be associated with OA. (C) Pleiotropy 2: HBM is associated with a ‘bone-forming’ tendency, which results in increased BMD and a higher risk of OA due to greater susceptibility to osteophyte formation—in this scenario, though BMD is not the trait that causally effects OA (a proximal intermediate trait is), it is less likely that there will be genetic variants associated with BMD that are not also associated with OA, making this scenario difficult to distinguish from scenario (A)

HBM and the ‘Bone-Forming’ Phenotype

There are several lines of evidence showing that individuals from the HBM cohort have a tendency towards excess bone formation. For example, an increased prevalence of pelvic radiographic enthesophytes (formation of extra bone at tendon or ligament insertions) was observed in these individuals [24]. In addition, mutations in LRP5 are thought to lead to HBM through over-stimulation of bone formation [25, 26]. Similarly, the SMAD9 mutation, which we recently identified in association with HBM, is predicted to cause over-stimulation of osteoblast differentiation [27]. Our HBM population also possesses an over-representation of BMD-related alleles, several of which are implicated in increased osteoblast activity, including the novel BMD-associated loci NPR3 and SPON1 [28]. Consistent with the suggestion of a phenotype resulting from excess bone formation, HBM cases show evidence of increased periosteal bone formation as evidenced by greater periosteal circumference as assessed by peripheral quantitative computer tomography (pQCT) [29].

As well as the suggestion that increased BMD in HBM individuals likely arises from increased bone formation, the same may apply to their OA phenotype, given the stronger association of HBM with osteophytes as opposed to JSN. OA is recognised to have several sub-phenotypes, and it may be that increased BMD in our HBM population is responsible for a predominantly osteophytic form of disease. That said, though osteophytes are a recognised feature of OA, their precise role in OA pathogenesis remains unclear. Though osteophytes are often observed in the context of cartilage loss and subsequent biomechanical changes, they may have an adaptive role in joint stabilisation rather than actively contributing to further deterioration. In addition, the extent to which osteophytes contribute to clinical features such as pain remains to be established. However, our recent finding that adjustment for knee osteophyte severity attenuated the relationship between HBM and knee pain scores provides some evidence that osteophytes contribute to pain in OA [22]. Moreover, HBM appears to be associated with true symptomatic OA, in light of their increased risk of joint replacement as discussed above [17].

Shared Genetic Factors Between Bone and OA

The study of bone genetics has proved useful in identifying shared regulatory pathways which might contribute to the association between OA and BMD/osteoporosis. The bone and cartilage share several developmental and regulatory pathways, including the Wnt pathway and transforming growth factor beta (TGF-β) superfamily [2]. Genetic disorders affecting the Wnt pathway have been implicated in both OA and osteoporosis (see below); likewise loci in several genes linked to the TGF-β superfamily have been identified as being related to the risk of OA, including GDF5, BMP2, TGFβ1, SMAD3 and ASPN [30-34]. Hence, common genetic influences on BMD and OA might be expected, as a consequence of genetic variation in these shared pathways. This question was examined specifically by Hackinger et al., who studied the genetic overlap between OA and BMD, based on genome-wide data from GEFOS (BMD) and arcOGEN (OA) [30]. Significant genetic correlation was observed between hip and/or knee OA and lumbar spine BMD, but not femoral neck BMD. SMAD3 was subsequently identified as one of the loci contributing to both traits. One could hypothesise that genetic correlation limited to lumbar spine BMD could reflect artefactual elevation of BMD in the presence of OA by spinal osteophytes. However, there has also been evidence for a genetic correlation between total body BMD measured in childhood and combined OA [30]; since total body BMD measurement during childhood is unaffected by OA, these findings suggest that there may be a true genetic correlation between BMD and OA, reflecting shared genetic aetiology. Shared genetic risk factors can indicate either that there are distinct effects of the variants on two separate traits or that one trait is causal for the other. Both scenarios would result in genetic correlation between the two traits, and further work is necessary to determine which is more likely.

Wnt Signalling

Several genetic variants associated with BMD relate to the Wnt signalling pathway. LRP5, which as described above causes familial HBM [25, 26], is a Wnt co-receptor, whilst a common polymorphism at the wnt16 locus is related to both BMD and fracture risk [35]. Wnt signalling also contributes to cartilage homeostasis and OA risk. For example, DOT1L, a histone methyltransferase, is involved in the downregulation of Wnt signalling as well as cartilage preservation [36]. Suggestive of a role of Wnt inhibition in reducing OA risk, a SNP in the DOT1L gene was associated with hip joint space width, an indirect measure of cartilage thickness, with genome-wide significance [37], and this same SNP has been associated with risk of hip OA [37, 38]. Further Wnt pathway inhibitors have been suggested to contribute to OA pathogenesis. For example, variants in the gene for the Wnt inhibitor, secreted frizzled-related protein 3, have been found to be associated with hip OA [39], and a novel small molecule Wnt inhibitor is being evaluated to treat knee OA [40].

The Wnt signalling pathway is complex, with different factors playing a predominant role in different tissue types. However, there is some evidence that certain factors within the Wnt pathway affect both the bone and cartilage. For example, mice with an LRP5 deletion display increased cartilage degeneration compared to wild-type mice [41], although analyses by a different group reported the opposite [42]. In studies of the contribution of Wnt16 to OA aetiology, Nalesso et al. observed more severe OA in Wnt16 knockout mice compared to that in wild type, suggesting Wnt16 has a role in preserving the cartilage after joint injury [43]. In contrast, Van den Bosch et al. found that Wnt16 overexpression in the synovium of a mouse model led to increased expression of markers of cartilage degeneration [44], whereas Tornqvist et al. found no change in OA severity in transgenic mice overexpressing Wnt16 in osteoblasts despite the increase in subchondral bone mass [45].

Sclerostin, a product of the SOST gene, was initially discovered as the underlying mutation in the HBM disorder sclerosteosis [46] and has since been targeted to treat postmenopausal osteoporosis in the form of romosozumab [47]. Sclerostin is produced by osteocytes and has been found to exert an important suppressive effect on osteoblastic bone formation [48]. Sclerostin binds to LRP5/LRP6 receptors to inhibit Wnt signalling and osteoblast activity. Several studies have examined a possible role of sclerostin in cartilage and OA development; however, studies of the contribution of sclerostin to OA have yielded somewhat conflicting results. Sost knockout mice have been found to show more severe cartilage lesions compared to wild-type mice [49]. This was supported by an earlier work of Chan et al. who provided further evidence that SOST is chondroprotective [50]. However, others have found no difference in cartilage lesions in SOST knockout mice [51], whilst most recently Zhou et al. identified more rapid OA development following surgical induction of OA in mice overexpressing Sost [52]. Sclerostin expression has been reported to be decreased in osteocytes of the subchondral bone in OA [50, 51, 53], and this decreased expression has been associated with subchondral bone thickening [50]. However, other studies have shown increased sclerostin staining in bone tissue and cartilage from tibial plateaus of patients who had undergone total knee replacement (TKR), compared to healthy control samples [54]. When investigating the relationship between circulating sclerostin and OA, one study found little evidence of an association between OA and serum sclerostin [55], whereas another reported reduced sclerostin levels in those with a higher radiographic severity of OA [56]. Consistent with a chondroprotective effect of sclerostin, in the HBM population described above, higher sclerostin levels were found to be associated with a lower risk of hip and thumb-base JSN [57].

Although these experimental studies provide some evidence that genes related to the Wnt pathway, such as LRP5, WNT16 and SOST, are involved in the pathogenesis of OA, much of the evidence is conflicting. Moreover, it is often difficult to distinguish separate effects of a specific pathway on the cartilage, from those mediated by changes in the bone. A further question remains as to whether biological pathways shared between the bone and cartilage contribute to OA by altering bone formation or resorption. As discussed above, increased subchondral bone turnover and bone loss have been implicated in OA progression. Moreover, as well as regulating osteoblasts, the Wnt pathway and sclerostin are known to influence osteoclastic bone resorption [58, 59].

Cathepsin K

Genetic variation in the osteoclast enzyme cathepsin K, which is required for bone resorption, was found to relate to the risk of developing OA in a recent GWAS, based on an association between a CTSK locus and prevalent OA [60]. This finding may reflect a protective effect of cathepsin K deficiency, and hence reduced bone resorption, on OA development, based on findings that a Ctsk knockout mouse showed delayed development of OA after surgery compared to wild-type mice [61]. Moreover, following a randomised controlled trial of a cathepsin K inhibitor, MIV-711, individuals with knee OA showed reduced MRI-assessed features of structural progression in the intervention groups, compared to placebo, as well as reduced levels of the cartilage breakdown marker uCTX-II [62]. However, rather than being mediated by altered bone resorption, these observations could reflect a direct influence of cathepsin K on cartilage: cathepsin K is also able to cleave type II collagen [63] and is expressed by human chondrocytes [63], and mRNA levels are increased in OA cartilage [64] and correlate with OA severity [65].

Use of Genetic Studies to Examine the Causal Effect of BMD on OA Risk

The causal relationship between BMD and OA has also been examined using genetic variants related to BMD in the general population. For example, in a study of knee OA cases and controls, Yerges-Armstrong et al. identified an association between four BMD-associated SNPs (identified by Estrada et al. in the GEFOS cohort [66]) and knee OA [67]. These four SNPs associated with knee OA were among those most strongly associated with BMD, supporting a causal relationship (one would expect, in the case of a true causal effect, for the SNP-outcome association to increase as the SNP-exposure association increases). Additionally, one of the four SNPs was annotated to the LRP5 gene implicated in familial HBM as described earlier [67]; the BMD increasing allele was associated with an increased odds of knee OA.

Mendelian Randomisation (MR) Studies

MR is a more formal analysis framework for examining causal relationships [68]. Though MR can be limited by lack of power due to weak instruments, several genome-wide association study (GWAS) signals for BMD have been identified, making it feasible to explore causal relationships between BMD and OA. For example, the GEFOS consortium identified 56 loci associated with hip and/or lumbar spine BMD [66], and 518 loci were identified using BMD estimated from heel ultrasound (eBMD) in the UK Biobank [35]. Though the latter represents an indirect estimate of BMD, relatively high genetic correlation exists between the two traits, suggesting eBMD represents a reasonable proxy measure for BMD [69]. MR relies on several assumptions, including the absence of pleiotropy (i.e. a given genetic instrument relates to the outcome solely via a causal effect mediated by the exposure), which several methods have been developed to address [70]. As discussed above, genetic influences on BMD may be associated with OA risk via a separate chondrocyte pathway, or a shared ‘bone-forming’ tendency, both of which represent pleiotropy. Similarly, pleiotropy might exist between genetic factors related to BMD and other characteristics known to influence OA, such as BMI.

Two previous MR studies have examined causal relationships between BMD and OA, using common genetic variants associated with BMD. Evidence of a causal effect of BMD on OA was found using an instrument derived from a lumbar spine BMD GWAS [71]. In a recent MR study, Funck-Brentano et al. identified a causal effect of BMD on both hip and knee OA, using an instrument based on genetic association with femoral neck BMD [72]. To exclude pleiotropy involving BMI pathways, Funck-Brentano et al. performed sensitivity analyses following the removal of BMI-associated SNPs, after which equivalent results were obtained. A number of other sensitivity analyses were also performed to explore possible pleiotropy, including MR-Egger regression and weighted median analysis. However, when two traits such as BMD and OA are highly co-regulated, it is difficult to detect pleiotropy using these methods. For example, the weighted median method assumes that the majority of SNPs are non-pleiotropic, which might not be the case for highly correlated variables such as BMD and OA.

As well as shared pathways between BMD and OA, causal pathways may also exist between these conditions and BMI. For example, given the suggested role for osteocalcin linking bone to energy metabolism [73], bone turnover may be related to adiposity, providing a further pathway between BMD and OA (Fig. 2). Hence, in a subsequent study, we examined each of these pathways in turn, accounting for possible relationships mediated by BMI using multivariable MR. Although we found evidence of an independent causal effect between BMD and OA, we also found evidence of a causal effect of OA on BMD, suggesting either a bi-directional relationship or pleiotropy that remained unaccounted for [74]. Multivariable effect estimates were similar to univariable estimates, providing little evidence of confounding by BMI. Therefore, though initial MR studies suggest that BMD has a causal effect on OA, the limited power of conventional pleiotropy tests, and our subsequent finding of bi-directional relationships, makes it difficult to exclude pleiotropy as an alternative explanation for these findings.

Fig. 2

Hypothesised causal diagram for relationships between BMD, BMI and OA. In addition to the a priori causal pathway between BMD and OA suggested by Funck-Brentano [72], we propose a possible bi-directional causal pathway between OA and BMD. BMI is thought to act as a confounder via the biomechanical effects of excess weight on BMD and OA. In addition, bi-directional causal pathways might exist between BMD and BMI (through increased bone turnover) and between OA and BMI (through reduced physical activity)

Future Work

Looking to the future, several newer MR methods are emerging to detect and account for pleiotropy [75-78], the application of which should enable the causal relationship between BMD and OA to be examined more thoroughly. We have also started to develop methods to identify which pathways contribute to the shared genetics of these two traits. For example, if only a small number of pathway loci are shared between BMD and OA, this would imply that BMD is not causal for OA and that the branch point in the causal pathway that produces pleiotropic effects occurs relatively early on (equivalent to the scenario on Fig. 1b). Alternatively, if all pathways studied have proportionate effects on both BMD and OA, this suggests either scenario A or C in Fig. 1 that either BMD really is causal for OA or an intermediate phenotype proximal to BMD is the true causal risk factor for OA, consistent with the role of a ‘bone-forming’ phenotype underlying both raised BMD and osteophytic forms of OA in the HBM cohort (Fig. 1c). This area of study will be facilitated by the availability of large GWAS cohorts linked to information about OA sub-phenotypes including osteophytes, making it possible to establish whether similar relationships exist in the general population.

Conclusions

As well as providing insights into shared genetic mechanisms between BMD and OA, bone genetics has been applied to examine the nature of relationships between BMD and OA. For example, the finding that individuals from our HBM cohort are at increased risk of developing OA raises the possibility that BMD plays a causal role in the development of OA. Although initial MR studies provided support for a causal effect of BMD on OA, further investigation has shown this to be more complicated, with either bi-directional or shared biological pathways also likely. Given the nature of the HBM phenotype and its genetic determinants, this may reflect the role of a ‘bone-forming’ phenotype in the development of OA. Further studies are required to dissect out the genetic pathways which contribute to the bone-OA relationship and to establish whether a similar ‘bone-forming’ phenotype underlies this relationship in more common forms of OA within the general population.