Cross-talk of subchondral bone osteoblasts and articular cartilage chondrocytes : a new insight in understanding osteoarthritis pathogenesis

Osteoarthritis (OA) is the most common musculoskeletal disorder and represents a major health burden to society. In the course of the pathological development of OA, articular cartilage chondrocytes (ACCs) undergo atypical phenotype changes characterized by the expression of hypertrophic differentiation markers. Also, the adjacent subchondral bone shows signs of abnormal mineral density and enhanced production of bone turnover markers, indicative of osteoblast dysfunction. Collectively these findings indicate that the pathological changes typical of OA, involve alterations of the phenotypic properties of cells in both the subchondral bone and articular cartilage. However, the mechanism(s) by which these changes occur during OA development are not completely understood. The purpose of this project was to address the question of how subchondral bone osteoblasts (SBOs) and ACCs interact with each other with respect to regulation of respective cells’ phenotypic properties and in particular the involvement of mitogen activated protein kinase (MAPK) signalling pathways under normal and OA joint condition. We also endeavoured to test the influence of cross-talk between SBOs and ACCs isolated from normal and OA joint on matrix metalloproteinase (MMP) expression. For this purpose tissues from the knees of OA patients and normal controls were collected to isolate SBOs and ACCs. The cellular cross-talk of SBOs and ACCs were studied by means of both direct and indirect co-culture systems, which made it possible to identify the role of both membrane bound and soluble factors. Histology, immunohistochemistry, qRT-PCR, zymography, ELISA and western blotting were some of the techniques applied to distinguish the changes in the co-cultured vs. non co-cultured cells. The MAPK signalling pathways were probed by using targeted MAPK inhibitors, and their activity monitored by western blot analysis using phospho MAPK specific antibodies. Our co-culture studies demonstrated that OA ACCs enhanced the SBOs differentiation compared to normal ACCs. We demonstrated that OA ACCs induced these phenotypic changes in the SBOs via activating an ERK1/2 signalling pathway. The findings from this study thus provided clear evidence that OA ACCs play an integral role in altering the SBO phenotype. In the second study, we tested the influence of normal SBOs and OA SBOs on ACCs phenotype changes. The results showed that OA SBOs increased the hypertrophic gene expression in co-cultured ACCs compared to normal SBOs, a phenotype which is considered as pathological to the health and integrity of articular cartilage. It was demonstrated that these phenotype changes occurred via de-activation of p38 and activation of ERK1/2 signaling pathways. These findings suggest that the pathological interaction of OA SBOs with ACCs is mediated by cross-talking between ERK1/2 and p38 pathways, resulting in ACCs undergoing hypertrophic differentiation. Subsequent experiments to determine the effect on MMP regulation, of SBOs and ACCs cross-talk, revealed that co-culturing OA SBOs with ACCs significantly enhanced the proteolytic activity and expression of MMP-2 and MMP-9. In turn, co-culture of OA ACCs with SBOs led to abundant MMP-2 expression in SBOs. Furthermore, we showed that the addition of ERK1/2 and JNK inhibitors reversed the elevated MMP-2 and MMP-9 production which otherwise resulted from the interactions of OA SBOs-ACCs. Thus, this study has demonstrated that the altered interactions between OA SBOs-ACCs are capable of triggering the pathological pathways leading to degenerative changes seen in the osteoarthritic joint. In conclusion, the body of work presented in this dissertation has given clear in vitro evidence that the altered bi-directional communication of SBOs and ACCs may play a role in OA development and that this process was mediated by MAPK signalling pathways. Targeting these altered interactions by the use of MAPK inhibitors may provide the scientific rationale for the development of novel therapeutic strategies in the treatment and management of OA.

Size and habit of mineral particles in bone and mineralized callus during bone healing in sheep

Bone healing is known to occur through the successive formation and resorption of various tissues with different structural and mechanical properties. To get a better insight into this sequence of events, we used environmental scanning electron microscopy (ESEM) together with scanning small-angle X-ray scattering (sSAXS) to reveal the size and orientation of bone mineral particles within the regenerating callus tissues at different healing stages (2, 3, 6, and 9 weeks). Sections of 200 µm were cut from embedded blocks of midshaft tibial samples in a sheep osteotomy model with an external fixator. Regions of interest on the medial side of the proximal fragment were chosen to be the periosteal callus, middle callus, intercortical callus, and cortex. Mean thickness (T parameter), degree of alignment (ρ parameter), and predominant orientation (ψ parameter) of mineral particles were deduced from resulting sSAXS patterns with a spatial resolution of 200 µm. 2D maps of T and ρ overlapping with ESEM images revealed that the callus formation occurred in two waves of bone formation, whereby a highly disordered mineralized tissue was deposited first, followed by a bony tissue with more lamellar appearance in the ESEM and where the mineral particles were more aligned, as revealed by sSAXS. As a consequence, degree of alignment and mineral particle size within the callus increased with healing time, whereas at any given moment there were structural gradients, for example, from periosteal toward the middle callus.

Spatial and temporal variations of mechanical properties and mineral content of the external callus during bone healing

After bone fracture, various cellular activities lead to the formation of different tissue types, which form the basis for the process of secondary bone healing. Although these tissues have been quantified by histology, their material properties are not well understood. Thus, the aim of this study is to correlate the spatial and temporal variations in the mineral content and the nanoindentation modulus of the callus formed via intramembranous ossification over the course of bone healing. Midshaft tibial samples from a sheep osteotomy model at time points of 2, 3, 6 and 9 weeks were employed. PMMA embedded blocks were used for quantitative back scattered electron imaging and nanoindentation of the newly formed periosteal callus near the cortex. The resulting indentation modulus maps show the heterogeneity in the modulus in the selected regions of the callus. The indentation modulus of the embedded callus is about 6 GPa at the early stage. At later stages of mineralization, the average indentation modulus reaches 14 GPa. There is a slight decrease in average indentation modulus in regions distant to the cortex, probably due to remodelling of the peripheral callus. The spatial and temporal distribution of mineral content in the callus tissue also illustrates the ongoing remodelling process observed from histological analysis. Most interestingly the average indentation modulus, even at 9 weeks, remains as low as 13 GPa, which is roughly 60% of that for cortical sheep bone. The decreased indentation modulus in the callus compared to cortex is due to the lower average mineral content and may be perhaps also due to the properties of the organic matrix which might be different from normal bone.

Altered cell interactions of subchondral bone osteoblasts and articular chondrocytes in osteoarthritis is through the mediation of ERK1/2 phosphorylation

Introduction: Osteoarthritis (OA) is the most common musculoskeletal disorder and represents a major health burden to society. In the course of the pathological development of OA, articular cartilage chondrocytes (ACCs) undergo a typical phenotype changes characterized by the expression of hypertrophic differentiation markers. Also, the adjacent subchondral bone shows signs of abnormal mineral density and enhanced production of bone turnover markers, indicative of osteoblast dysfunction. However, the mechanism(s) by which these changes occur during the OA development are not completely understood. Materials and Methods: ACCs and subchondral bone osteoblasts (SBOs) were harvested from OA and healthy patients for the cross-talk studies between normal and OA ACCs and SBOs. The involvement of mitogen activated protein kinase (MAPK) signalling pathway during the cell-cell interactions was analysed by zymography, ELISA and western blotting methods. Results: The direct and in-direct co-culture studies showed that OA (ACCs and SBOs) cells induced osteoarthritic changes of normal (ACC and SBOs) cells. This altered cell interaction induced by OA cells significantly aggravated the proteolytic activity, which resulted cartilage degeneration. The altered cell interaction appeared to significantly activate ERK 1/2 phosphorylation and inhibition of MAPK-ERK 1/2 pathway reversed the osteoarthrtitic phenotypic changes. Discussion and Conclusion: Our study has demonstrated that the altered bi-directional communication of SBOs and ACCs are critical for initiation and progression of OA related changes and that this process is mediated by MAPK signalling pathways. Targeting these altered interactions by the use of MAPK inhibitors may provide the scientific rationale for the development of novel therapeutic strategies in the treatment and management of OA related disorders.

The global burden attributable to low bone mineral density

Introduction The Global Burden of Disease Study 2010 estimated the worldwide health burden of 291 diseases and injuries and 67 risk factors by calculating disability-adjusted life years (DALYs). Osteoporosis was not considered as a disease, and bone mineral density (BMD) was analysed as a risk factor for fractures, which formed part of the health burden due to falls. Objectives To calculate (1) the global distribution of BMD, (2) its population attributable fraction (PAF) for fractures and subsequently for falls, and (3) the number of DALYs due to BMD. Methods A systematic review was performed seeking population-based studies in which BMD was measured by dual-energy X-ray absorptiometry at the femoral neck in people aged 50 years and over. Age- and sex-specific mean ± SD BMD values (g/cm2) were extracted from eligible studies. Comparative risk assessment methodology was used to calculate PAFs of BMD for fractures. The theoretical minimum risk exposure distribution was estimated as the age- and sex-specific 90th centile from the Third National Health and Nutrition Examination Survey (NHANES III). Relative risks of fractures were obtained from a previous meta-analysis. Hospital data were used to calculate the fraction of the health burden of falls that was due to fractures. Results Global deaths and DALYs attributable to low BMD increased from 103 000 and 3 125 000 in 1990 to 188 000 and 5 216 000 in 2010, respectively. The percentage of low BMD in the total global burden almost doubled from 1990 (0.12%) to 2010 (0.21%). Around one-third of falls-related deaths were attributable to low BMD. Conclusions Low BMD is responsible for a growing global health burden, only partially representative of the real burden of osteoporosis.

The impact of intense training on endogenous estrogen and progesterone concentrations and bone mineral acquisition in adolescent rowers

The effect of 18 months of training on the ovarian hormone concentrations and bone mineral density (BMD) accrual was assessed longitudinally in 14 adolescent rowers and 10 matched controls, aged 14–15 years. Ovarian hormone levels were assessed by urinary estrone glucuronide (E1G) and pregnanediol glucuronide (PdG) excretion rates, classifying the menstrual cycles as ovulatory or anovulatory. Total body (TB), total proximal femur (PF), femoral neck (FN) and lumbar spine (LS) (L2–4) bone mass were measured at baseline and 18 months using dual-energy X-ray densitometry. Results were expressed as bone mineral content (BMC), BMD and bone mineral apparent density (BMAD). Five rowers had anovulatory menstrual cycles compared with zero prevalence for the control subjects. Baseline TB BMD was significantly higher in the ovulatory rowers, with PF BMD, FN BMD and LS BMD similar for all groups. At completion, the LS bone accrual of the ovulatory rowers was significantly greater (BMC 8.1%, BMD 6.2%, BMAD 6.2%) than that of the anovulatory rowers (BMC 1.1%, BMD 3.9%, BMAD 1.6%) and ovulatory controls (BMC 0.5%, BMD 1.1%, BMAD 1.1%). No difference in TB, PF or FN bone accrual was observed among groups. This study demonstrated an osteogenic response to mechanical loading, with the rowers accruing greater bone mass than the controls at the lumbar spine. However, the exercise-induced osteogenic benefits were less when rowing training was associated with low estrogen and progesterone metabolite excretion.

Changes in cervical spine bone mineral density in response to flight training

BACKGROUND High magnitude loads and unusual loading regimes are two important determinants for increasing bone mass. Past research demonstrated that positive Gz-induced loading, providing high loads in an unaccustomed manner, had an osteogenic effect on bone. Another determinant of bone mass is that the bone response to loading is site specific. This study sought to further investigate the site specific bone response to loading, examining the cervical spine response, the site suspected of experiencing the greatest loading, to high performance flight. METHODS Bone mineral density (BMD) and bone mineral content (BMC) was monitored in 9 RAAF trainee fighter pilots completing an 8-mo flight training course on a PC-9 and 10 age-height-weight-matched controls. RESULTS At completion of the course, the pilots had a significant increase in cervical spine BMD and total body BMC. No significant changes were found for the control group. CONCLUSIONS This study demonstrated that the physical environment associated with flight training may have contributed to a significant increase in cervical spine bone mass in the trainee PC-9 pilots. The increase in bone mass was possibly a response to the strain generated by the daily wearing of helmet and mask assembly under the influence of positive sustained accelerative forces.

Genome-wide association study using extreme truncate selection identifies novel genes affecting bone mineral density and fracture risk

Osteoporotic fracture is a major cause of morbidity and mortality worldwide. Low bone mineral density (BMD) is a major predisposing factor to fracture and is known to be highly heritable. Site-, gender-, and age-specific genetic effects on BMD are thought to be significant, but have largely not been considered in the design of genome-wide association studies (GWAS) of BMD to date. We report here a GWAS using a novel study design focusing on women of a specific age (postmenopausal women, age 55-85 years), with either extreme high or low hip BMD (age- and gender-adjusted BMD z-scores of +1.5 to +4.0, n = 1055, or -4.0 to -1.5, n = 900), with replication in cohorts of women drawn from the general population (n = 20,898). The study replicates 21 of 26 known BMD-associated genes. Additionally, we report suggestive association of a further six new genetic associations in or around the genes CLCN7, GALNT3, IBSP, LTBP3, RSPO3, and SOX4, with replication in two independent datasets. A novel mouse model with a loss-of-function mutation in GALNT3 is also reported, which has high bone mass, supporting the involvement of this gene in BMD determination. In addition to identifying further genes associated with BMD, this study confirms the efficiency of extreme-truncate selection designs for quantitative trait association studies. © 2011 Duncan et al.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ~2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2×10-9). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3×10-12, and -0.16 SD per G allele, P = 1.2×10-15, respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3×10-9), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9×10-6 and rs2707466: OR = 1.22, P = 7.2×10-6). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16-/- mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5×10-13<P<5.9×10-4) at both femur and tibia, compared with their wild-type littermates. Natural variation in humans and targeted disruption in mice demonstrate that WNT16 is an important determinant of CBT, BMD, bone strength, and risk of fracture. © 2012 Zheng et al.

Suggestive linkage of 2p22-25 and 11q12-13 with low bone mineral density at the lumbar spine in the Irish population

Osteoporosis is a disease characterized by low bone mineral density (BMD) and poor bone quality. Peak bone density is achieved by the third decade of life, after which bone is maintained by a balanced cycle of bone resorption and synthesis. Age-related bone loss occurs as the bone resorption phase outweighs the bone synthesis phase of bone metabolism. Heritability accounts for up to 90% of the variability in BMD. Chromosomal loci including 1p36, 2p22-25, 11q12-13, parathyroid hormone receptor type 1 (PTHR1), interleukin-6 (IL-6), interleukin 1 alpha (IL-1α) and type II collagen A1/vitamin D receptor (COL11A1/VDR) have been linked or shown suggestive linkage with BMD in other populations. To determine whether these loci predispose to low BMD in the Irish population, we investigated 24 microsatellite markers at 7 chromosomal loci by linkage studies in 175 Irish families of probands with primary low BMD (T-score ≤ -1.5). Nonparametric analysis was performed using the maximum likelihood variance estimation and traditional Haseman-Elston tests on the Mapmaker/Sibs program. Suggestive evidence of linkage was observed with lumbar spine BMD at 2p22-25 (maximum LOD score 2.76) and 11q12-13 (MLS 2.55). One region, 1p36, approached suggestive linkage with femoral neck BMD (MLS 2.17). In addition, seven markers achieved LOD scores > 1.0, D2S149, D11S1313, D11S987, D11S1314 including those encompassing the PTHR1 (D3S3559, D3S1289) for lumbar spine BMD and D2S149 for femoral neck BMD. Our data suggest that genes within a these chromosomal regions are contributing to a predisposition to low BMD in the Irish population.

COL1A1 C-propeptide cleavage site mutation causes high bone mass, bone fragility and jaw lesions: A new cause of gnathodiaphyseal dysplasia?

Gnathodiaphyseal dysplasia (GDD) is a rare autosomal dominant condition characterized by bone fragility, irregular bone mineral density (BMD) and fibro-osseous lesions in the skull and jaw. Mutations in Anoctamin-5 (ANO5) have been identified in some cases. We aimed to identify the causative mutation in a family with features of GDD but no mutation in ANO5, using whole exome capture and massive parallel sequencing (WES). WES of two affected individuals (a mother and son) and the mother's unaffected parents identified a mutation in the C-propeptide cleavage site of COL1A1. Similar mutations have been reported in individuals with osteogenesis imperfecta (OI) and paradoxically increased BMD. C-propeptide cleavage site mutations in COL1A1 may not only cause 'high bone mass OI', but also the clinical features of GDD, specifically irregular sclerotic BMD and fibro-osseous lesions in the skull and jaw. GDD patients negative for ANO5 mutations should be assessed for mutations in type I collagen C-propeptide cleavage sites.

Meta-analysis of genome-wide studies identifies MEF2C SNPs associated with bone mineral density at forearm

Background Forearm fractures affect 1.7 million individuals worldwide each year and most occur earlier in life than hip fractures. While the heritability of forearm bone mineral density (BMD) and fracture is high, their genetic determinants are largely unknown. Aim To identify genetic variants associated with forearm BMD and forearm fractures. Methods BMD at distal radius, measured by dualenergy x-ray absorptiometry, was tested for association with common genetic variants. We conducted a metaanalysis of genome-wide association studies for BMD in 5866 subjects of European descent and then selected the variants for replication in 715 Mexican American samples. Gene-based association was carried out to supplement the single-nucleotide polymorphism (SNP) association test. We then tested the BMD-associated SNPs for association with forearm fracture in 2023 cases and 3740 controls. Results We found that five SNPs in the introns of MEF2C were associated with forearm BMD at a genome-wide significance level (p<5×10-8) in meta-analysis (lead SNP, rs11951031[T] -0.20 SDs per allele, p=9.01×10-9). The gene-based association test suggested an association between MEF2C and forearm BMD ( p=0.003). The association between MEF2C variants and risk of fracture did not achieve statistical significance (SNP rs12521522[A]: OR=1.14 (95% CI 0.92 to 1.35), p=0.14). Meta-analysis also revealed two genome-wide suggestive loci at CTNNA2 and 6q23.2. Conclusions These findings demonstrate that variants at MEF2C were associated with forearm BMD, implicating this gene in the determination of BMD at forearm.

'Sink or swim': An evaluation of the clinical characteristics of individuals with high bone mass

Summary High bone mineral density on routine dual energy X-ray absorptiometry (DXA) may indicate an underlying skeletal dysplasia. Two hundred fifty-eight individuals with unexplained high bone mass (HBM), 236 relatives (41% with HBM) and 58 spouses were studied. Cases could not float, had mandible enlargement, extra bone, broad frames, larger shoe sizes and increased body mass index (BMI). HBM cases may harbour an underlying genetic disorder. Introduction High bone mineral density is a sporadic incidental finding on routine DXA scanning of apparently asymptomatic individuals. Such individuals may have an underlying skeletal dysplasia, as seen in LRP5 mutations. We aimed to characterize unexplained HBM and determine the potential for an underlying skeletal dysplasia. Methods Two hundred fifty-eight individuals with unexplained HBM (defined as L1 Z-score ≥ +3.2 plus total hip Z-score ≥ +1.2, or total hip Z-score ≥ +3.2) were recruited from 15 UK centres, by screening 335,115 DXA scans. Unexplained HBM affected 0.181% of DXA scans. Next 236 relatives were recruited of whom 94 (41%) had HBM (defined as L1 Z-score + total hip Z-score ≥ +3.2). Fifty-eight spouses were also recruited together with the unaffected relatives as controls. Phenotypes of cases and controls, obtained from clinical assessment, were compared using random-effects linear and logistic regression models, clustered by family, adjusted for confounders, including age and sex. Results Individuals with unexplained HBM had an excess of sinking when swimming (7.11 [3.65, 13.84], p < 0.001; adjusted odds ratio with 95% confidence interval shown), mandible enlargement (4.16 [2.34, 7.39], p < 0.001), extra bone at tendon/ligament insertions (2.07 [1.13, 3.78], p = 0.018) and broad frame (3.55 [2.12, 5.95], p < 0.001). HBM cases also had a larger shoe size (mean difference 0.4 [0.1, 0.7] UK sizes, p = 0.009) and increased BMI (mean difference 2.2 [1.3, 3.1] kg/m 2, p < 0.001). Conclusion Individuals with unexplained HBM have an excess of clinical characteristics associated with skeletal dysplasia and their relatives are commonly affected, suggesting many may harbour an underlying genetic disorder affecting bone mass.

Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture

Bone mineral density (BMD) is the most widely used predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and east Asian ancestry. We tested the top BMD-associated markers for replication in 50,933 independent subjects and for association with risk of low-trauma fracture in 31,016 individuals with a history of fracture (cases) and 102,444 controls. We identified 56 loci (32 new) associated with BMD at genome-wide significance (P < 5 × 10−8). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal stem cell differentiation, endochondral ossification and Wnt signaling pathways. However, we also discovered loci that were localized to genes not known to have a role in bone biology. Fourteen BMD-associated loci were also associated with fracture risk (P < 5 × 10−4, Bonferroni corrected), of which six reached P < 5 × 10−8, including at 18p11.21 (FAM210A), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.

Common variants in the region around Osterix are associated with bone mineral density and growth in childhood

Peak bone mass achieved in adolescence is a determinant of bone mass in later life. In order to identify genetic variants affecting bone mineral density (BMD), we performed a genome-wide association study of BMD and related traits in 1518 children from the Avon Longitudinal Study of Parents and Children (ALSPAC). We compared results with a scan of 134 adults with high or low hip BMD. We identified associations with BMD in an area of chromosome 12 containing the Osterix (SP7) locus, a transcription factor responsible for regulating osteoblast differentiation (ALSPAC: P = 5.8 × 10-4; Australia: P = 3.7 × 10-4). This region has previously shown evidence of association with adult hip and lumbar spine BMD in an Icelandic population, as well as nominal association in a UK population. A meta-analysis of these existing studies revealed strong association between SNPs in the Osterix region and adult lumbar spine BMD (P = 9.9 × 10-11). In light of these findings, we genotyped a further 3692 individuals from ALSPAC who had whole body BMD and confirmed the association in children as well (P = 5.4 × 10-5). Moreover, all SNPs were related to height in ALSPAC children, but not weight or body mass index, and when height was included as a covariate in the regression equation, the association with total body BMD was attenuated. We conclude that genetic variants in the region of Osterix are associated with BMD in children and adults probably through primary effects on growth.