Skip to main content

WNT16 elevation induced cell senescence of osteoblasts in ankylosing spondylitis



WNT16 is critical for bone homeostasis, but the effect of WNT16 in ankylosing spondylitis (AS) is still unknown. Here, we investigated whether WNT16 influences bone formation and pathophysiological changes of AS in an in vitro model.


The bone tissue from the facet joints was obtained from seven disease control and seven AS patients. Primary osteoprogenitor cells of the facet joints were isolated using an outgrowth method. Isolated osteoprogenitor cells from both control and AS tissues were analyzed by microarray, RT-qPCR, immunoblotting, and immunohistochemistry. The bone-forming activity of osteoprogenitor cells was assessed by various in vitro assays. β-galactosidase staining and senescence-associated secretory phenotype (SASP) using RT-qPCR were used to assess cell senescence.


In microarray analysis, WNT16 expression was significantly elevated in AS osteoprogenitor cells compared to the control. We also validated that WNT16 expression was elevated in AS-osteoprogenitor cells and human AS-bone tissues. WNT16 treatment inhibited bone formation in AS-osteoprogenitor cells but not in the control. Intriguingly, AS-osteoprogenitor cells were stained markedly with β-galactosidase for cell senescence in WNT16 treatment. Furthermore, in an H2O2 stress-induced premature senescence condition, WNT16 treatment increased cell senescence in AS-osteoprogenitor cells and WNT16 treatment under the H2O2 stress condition showed an increase in p21 protein and SASP mRNA expression. The WNT16-induced SASP expression in AS-osteoprogenitor cells was reduced in WNT16 knockdown cultures.


WNT16 is highly expressed in AS and WNT16 treatment facilitated cell senescence in AS-osteoprogenitor cells during osteoblast differentiation accompanied by suppression of bone formation. The identified role of WNT16 in AS could influence bone loss in AS patients.


Ankylosing spondylitis (AS) is a chronic inflammatory arthritis characterized by spinal ankylosis due to excessive bone formation [1,2,3]. Clinically, AS has shown a paradoxical feature including bone loss as well as bone formation simultaneously [4,5,6,7,8,9,10]. Pathological bone changes in AS have been associated with the Wnt signaling pathway [11,12,13] raising intriguing questions of the molecular partners of Wnt signaling and bone changes in AS. Among the Wnt molecules, many studies have clearly shown that WNT16 signaling increases bone formation while WNT16 inhibition leads to decreased bone formation [14,15,16]. Although the role of WNT16 signaling and bone formation has been studied using various genetic murine models, the conclusions in human AS are still elusive.

The main cause of bone loss in patients with AS is unclear. It has been reported that the prevalence of AS patients with bone loss is about 25% and low vertebral bone mineral density in AS is closely associated with high modified Stoke Ankylosing Spondylitis Spinal Score (mSASSS) and C-reactive protein (CRP) [17, 18]. Chronic inflammation and immobilization have been speculated as risk factors in bone loss in AS patients. A recent study reported increased cellular senescence in mesenchymal stem cells (MSCs) derived from the bone marrow of patients with AS [19]. In this paper, oxidative stress (OS) facilitated senescence by an increase in the release of reactive oxygen species (ROS) and mitochondrial deregulation. Clinically, oxidative stress status and associated serum proteins were higher in AS and associated with low bone mineral density in patients with AS [20,21,22]. However, whether the senescence influences bone-forming activity in AS is still unknown.

Elevated expression of Wnt/β-catenin signaling has been reported to be associated with AS pathophysiology. Wnt molecules have been studied a lot with bone formation, but little is known about the association of Wnt with osteoporotic bone loss in AS. In the present study, we evaluated whether WNT16 could regulate bone formation and bone loss in AS. Moreover, the physiological changes caused by WNT16 were examined in AS-osteoprogenitor cells.

Materials and methods

Human facet joint bone tissue

The present study was conducted in accordance with the Helsinki Declaration and approved by the Institutional Review Board of Hanyang University, Seoul (IRB 2014-05-002), and Guri Hospital (IRB 2014-05-001). We obtained facet joints from seven AS patients who met modified New York criteria [23, 24] and underwent spinal surgery and from seven patients with non-inflammatory disease as the control. Clinical demographics are shown in Table 1.

Table 1 Clinical demographics in this study

Isolation of osteoprogenitor cells

The protocol for the isolation of osteoprogenitor cells was previously reported [25,26,27,28]. Briefly, osteoprogenitor cells were isolated until the fourth passage using the outgrowth method for one month in Dulbecco’s modified Eagle medium-high glucose (DMEM-HG; Hyclone, SH30243.01) containing 10% fetal bovine serum (FBS; Gibco, 16000-044) and 1% antibiotics (penicillin/streptomycin; Gibco, 1540-122). All isolated osteoprogenitor cells were checked for mycoplasma using a PCR-based method (Takara, 6601) before long-term storage in liquid nitrogen. Passages 3~7 were selected and used for further experiments.

Microarray analysis

Passage 3 of control- and AS-osteoprogenitors was used and performed using the Affymetrix GeneChip Human Gene 2.0 ST array. For microarray analysis, background adjustment, normalization, and median polish summarization of CEL files were performed with the robust multi-average method using Partek Genomics Suite software version 7.0. Differential expression was calculated from the normalized data using ANOVA and visualized as a heatmap for Wnt genes.


Primary osteoprogenitor cells were lysed using TRIzol, and 1.0 μg of the total RNA was subjected to reverse transcription using a cDNA kit (Thermo Fisher Scientific, Waltham, MA, USA). Relative gene transcript levels were determined by RT-qPCR using SYBR Green, and qPCR was performed in triplicate for each sample. The mean expression level of each gene was normalized to GAPDH. The primers used in the study are listed in Table 2.

Table 2 RT-qPCR primers

Immunoblot assays

Stimulated osteoprogenitor cells were lysed with 1X RIPA buffer containing phosphatase (5870S, Cell Signaling Technology, USA) and protease (535140, Calbiochem, USA) inhibitors. Proteins were quantified using the Bradford assay. Total cell protein (10–30 μg) was subjected to immunoblotting using the following antibodies: β-catenin (9562S), active β-catenin (19807S), phos-AKT (Ser473) (4060), total-AKT (4691), phos-ERK (9101), total-ERK (9102), and GAPDH (2118) from Cell Signaling Technology (Danvers, MA, USA). CDKN1A (p21; A19094) and OCN (ab13420) antibodies were obtained from Abaclonal and abcam, respectively. WNT16 (sc-271897), CASP3 (sc-7272), and OPN (sc-21742) antibodies were obtained from Santa Cruz Biotechnology. Goat anti-rabbit IgG (111-035-003) and goat anti-mouse IgG (115-035-003) secondary antibodies were obtained from Jackson ImmunoResearch (West Grove, PA, USA).

Cell viability

The water-soluble tetrazolium salt (WST) assay was performed using the EZ-CYTOX kit (Dogen Bio, EZ-1000) according to the manufacturer’s instructions. Primary osteoprogenitor cells were seeded into 96-well plates at 1×103 cells/well and treated with WNT16 at the indicated dose and time point (days). At the indicated day, 10 μl WST solution was added to the cells and incubated for 1 h. Absorbance at 450 nm was measured using a microplate reader (Thermo Fisher Scientific, Inc.).


The detailed materials and methods for immunohistochemistry (IHC) analysis were previously described [25]. Briefly, facet joint tissues were fixed in 10% formalin for a week, decalcified with 10% formic acid for a week, and embedded in paraffin. Tissue slides (5 μm) were treated with Neo-Clear for deparaffination, followed by dehydration with ethanol solutions (100 to 50% of ethanol), permeabilization with 0.25% Triton X-100, and blocking with BLOXALL (Vector Lab, SP-6000) for 1 h. The slides were then incubated with primary WNT16 antibodies (sc271897, Santa Cruz Biotechnology; 1:100 ratio in antibody diluent) overnight, followed by washing with PBS-T, incubation with biotinylated secondary antibodies for 1 h, washing with PBS-T, incubation with ABC kit components as specified by the manufacturer (Vector Lab, PK-6102), and incubation using a DAB substrate kit (Vector Lab, sk4100) for 1 to 5 min. The slides were counterstained with haematoxylin (Merck, 1.05174.0500) for 10 s, dehydrated though ethanol and Neo-Clear, and coverslipped using a permanent mounting medium (Vector Lab, H-5000). Microscopy was used to visualize and image stain cells (Leica).

Human WNT16 serum level measurements

Sera were collected from 30 healthy male donors and 36 patients with AS who met the modified New York criteria [23]. Serum collection methods and storage conditions were previously described [25]. We measured WNT16 serum level using a human WNT16 ELISA kit (CUSABIO, CSB-EL026132HU) according to the manufacturer’s instructions.

In vitro bone formation

The protocol for the use of osteoprogenitor cells for assessing bone-forming activity was previously reported [25,26,27,28]. Briefly, osteoprogenitor cells were seeded in 96-well plates at 1×104 cells/well and incubated for 24 h. These cells were changed to osteogenic medium containing DMEM-HG (Hyclone, SH30243.01) containing 10% fetal bovine serum (FBS; Gibco, 16000-044), 1% antibiotics (penicillin/streptomycin; Gibco, 1540-122), 50 μM ascorbic acid (Sigma, A4544), 10 mM β-glycerol phosphate (Santa Cruz Biotechnology, sc-220452A;), and 100 nM dexamethasone (Sigma, D2915). The medium was changed every 3 or 4 days and bone-forming activity was assessed at the indicated days. The matrix maturation of bone formation was assessed by alkaline phosphatase (ALP) activity (BioVision, K412-500), ALP (Sigma, 85L2), and collagen using Picro-Sirius Red (Abcam, ab150681) staining. The matrix mineralization of bone formation was assessed using Alizarin Red S (ARS; ScienCell, 8678), hydroxyapatite (HA; Lonza, PA-1503), and Von Kossa (5% silver nitrate solution; SAMCHUN Chemical, S0228) staining. Stained wells were imaged via microscopy (Nikon eclipse Ti-U, Nikon). The overall image of a well for 96 well plates and its enlarged changes were seen in Suppl. Figure 1.

In the quantification of bone mineralization, we followed the manufacturer’s instructions for ARS and HA staining. To analyze Von Kossa staining, each image was manually imported to the ImageJ, defined the area by CIRCLE or SQUARE; the image was segmented and then put on the “Analyze” and “Measure” to obtain the area percentage of each pre-defined class. Before calculation, the VON image value at day 0 was also prepared in the same way as above. Next, we collected the “Area” and “Percent” values, put them in Excel, and calculated the Relative Density values for the VON image at day 0 as a control.

Promoter assays

SaOS2 cells (Korean Cell Line Bank, 30085) were seeded in 60 mm plates at 2×105 cells/well and co-transfected with each alkaline phosphatase (ALP), osteoblast-specific element (OSE), bone sialoprotein (BSP), osteocalcin (OCN), TOP, or FOP promoter DNA (3 μg/well) and Renilla luciferase (0.3 μg/well) using Lipofectamine 3000 (Thermo Fisher, L3000015). A day after transfection, the cells were reseeded in 12-well plates at 5×104 cells/well and incubated for 24 h. The next day, cells were treated with WNT16 (50 ng/mL) for 24 h and lysed to measure luciferase activity (Promega, E1500), according to the manufacturer’s instructions. The luciferase activity data were normalized to Renilla luciferase activity.

Small interference RNA (siRNA) and knockdown efficiency testing

Five siRNA oligos against human WNT16 were designed and generated by GENOLUTION Inc. (Seoul, Korea). Oligo sequences for siRNA are listed in Table 3. To test WNT16 knockdown efficiency, FOB cells (ATCC, CRL-11372) were seeded in 6-well plates and transfected with one of the five siRNA oligos against WNT16 using Lipofectamine 3000 (Thermo Fisher, L3000015). After 48 h, the transfected FOB cells were lysed and analyzed by RT-PCR (Suppl. Figure 2). Selected siRNA oligos were used for further investigation.

Table 3 siRNA oligos against WNT16

Senescence staining with β-galactosidase

To test the effect of WNT16 on cell senescence, AS-osteoprogenitor cells were seeded in 96-well plates at 4×103 cells/well and incubated for 24 h. The next day, the cells were stimulated H2O2 (200 μM) for 2 h, washed with 1X PBS, and suspended in new growth medium with WNT16 (50 ng/mL) or vehicle. After 3 days, the senescence effect of WNT16 was assessed using a senescence detection kit (BioVision, K320), according to the manufacturer’s instructions. After staining, stained wells were imaged by microscopy (Nikon eclipse Ti-U). For senescence quantification, stained images were analyzed by counting senescence positive cells.

Statistical analysis

Data images were generated using GraphPad Prism 6 software (GraphPad Software, Inc.). Two-tailed Student’s t tests were used to compare data between two unpaired groups. All data are expressed as the median and interquartile range (n≥3). A p<0.05 was considered statistically significant.


WNT16 is highly expressed in the facet joint of AS

We performed microarray experiments with control and AS-osteoprogenitor cells and analyzed the data to attain molecular insights among the Wnt genes in those samples. Differential expression gene (DEG) analysis revealed that WNT1, WNT3, and WNT16 were highly expressed in AS-osteoprogenitor cells (Fig. 1A, red arrows). We validated these data using RT-qPCR and showed that WNT1 was not detected, but WNT3, WNT5, and WNT16 were significantly increased in AS (Fig. 1B). Since the functional role of WNT3 and WNT5 in AS had been reported elsewhere [3], we focused the role of WNT16 in AS-osteoprogenitor cells. WNT16 protein was found to be higher in AS-osteoprogenitor cells relative to the control, which was revealed by immunoblotting and its quantification (Fig. 1C). We confirmed that there was no significant difference in WNT16 serum levels between healthy and AS groups (Suppl. Figure 3), but we consistently found increased WNT16 expression in bone-lining cells (black arrows) and in the facet joint periosteum tissue (blue arrows) in AS-osteoprogenitor cells compared with the control (Fig. 1D). Collectively, WNT16 was significantly elevated in AS-osteoprogenitor cells and facet joint of AS.

Fig. 1
figure 1

WNT16 is highly expressed in facet joint of AS. A Heatmap diagram of differentially expressed probes mapped to Wnt genes in AS and control osteoprogenitor cells. In AS osteoprogenitor cells, WNT16 showed significant upregulated expression with a p value < 0.05 and fold change > 2 compared to the control (n=3 per each group). B Verification of mRNA level using qPCR of A for the Wnt family (n=6 per each group). C Immunoblotting for WNT16 protein level in control and AS-osteoprogenitor cells (n=5 per each group). D Immunohistochemistry (IHC) staining of WNT16 in spinal facet joint tissue from control or AS patients. Black arrows indicate bone-lining cells, and blue arrows indicate periosteum in facet joints tissues. Three representative images are shown (n=7 per each group). Data are presented as the median and interquartile range. *p<0.05, **p<0.01

WNT16 treatment decreased matrix mineralization during osteoblast differentiation and increased senescence in AS

We next investigated whether WNT16 treatment influences matrix maturation and mineralization of osteoprogenitor cells during osteogenic differentiation. Two doses (25 and 50 ng/ml) of WNT16 did not affect the proliferation rate of AS-osteoprogenitor cells (Suppl. Figure 4), but treatment with 50 ng/mL WNT16 inhibited ALP activity and staining status in AS-osteoprogenitor cells. To investigate the effect of WNT16, both control and AS-osteoprogenitor cells were treated with 50 ng/mL WNT16 during osteoblast differentiation. As shown in Fig. 2A–E, treatment with WNT16 significantly inhibited bone matrix maturation (ALP staining and activity) and matrix mineralization (ARS, VON, and HA staining) in AS-osteoprogenitor during osteogenic differentiation. There is no substantial observation in control-osteoprogenitor, but these inhibitory effects by WNT16 were pronounced in AS-osteoprogenitor cells. WNT16 treatment also reduced bone formation-related gene promoter activity in alkaline phosphatase (ALP), bone sialoprotein (BSP), and osteocalcin (OCN), as revealed by luciferase activity (Fig. 2F). Promoter activity of osteoblast-specific element (OSE) and mRNA expression of Osterix were no substantial changes in WNT16 treatment (data not shown). Intriguingly, AS-osteoprogenitor cells were markedly stained with SA-β-gal at 28 days of differentiation when compared with the control although the cell staining rates were not significantly different (Fig. 2G). As shown in the 2H, caspase 3 protein was upregulated in the vehicle group of AS-osteoprogenitor cells accompanied by gain of bone formation-related protein such as OCN and OPN during osteogenic differentiation, whereas the WNT16 treatment group inhibited the upregulation of those expressions and promoted CDKN1A (called as p21, cell cycle regulator as wells senescence indicator) expression in AS-osteoprogenitor. The RT-qPCR results of OCN, CDKN1A, and CASP3 were consistent with Fig. 2H (Fig. 2I). These findings show that treatment with WNT16 inhibits matrix maturation and mineralization during osteogenic differentiation and promotes cell senescence in AS-osteoprogenitor cells.

Fig. 2
figure 2

WNT16 treatment inhibited mineralization and promoted cell senescence of AS-osteoprogenitor cells during osteoblast differentiation. Both control and AS-osteoprogenitor cells were differentiated with osteogenic medium in the presence of vehicle or WNT16. Analysis of A ALP staining, B ARS staining, and VON staining, C HA staining, D intercellular ALP activity at indicated day, E quantitation of ARS and VON at day 21, and HA staining at day 28 (n=4 per each group). F SaOS2 cells were co-transfected with 3 μg ALP, BSP, OSE, or OCN promoter plasmid along with β-gal plasmid for 48 h followed by and WNT16 treatment for 24 h and then analyzed by luciferase assay (n=6 per each group). G SA-β-gal staining was performed at 28 days of osteoblast differentiation of control and AS-osteoprogenitor cells. Counting results of Fig. 3G (right panel) (n=4 per each group). AS-osteoprogenitor cells were differentiated into mature osteoblasts in the presence of vehicle and WNT16 (50ng/ml). Analysis of H immunoblotting for protein level and I RT-qPCR for mRNA level (n=3 per each group). Scale bar = 200 μm. Data are presented as the median and interquartile range. *p<0.05

WNT16 treatment facilitated cell senescence of AS-osteoprogenitor cells under H2O2 stimulation

We next sought to investigate AS without osteogenic differentiation and designed an experimental condition using H2O2-induced senescence (Fig. 3A). WNT16 treatment in AS-osteoprogenitor cells exhibited a marked increase in SA-β-gal staining, SA-β-gal positive cells rates, and hydrogen peroxidase activity (Fig. 3B–D). In these conditions, WNT16 treatment increased active β-catenin, β-catenin, and p21 protein expression in AS-osteoprogenitor cells and decreased phos-AKT without changes in phos-p38 or phos-ERK protein levels (Fig. 3E, F). Moreover, WNT16 treatment significantly increased p21 mRNA expression and senescence-associated secretory phenotype (SASP)-related genes such as IL-1α, IL-1β, IL-6, and IL-8 (Fig. 3G) in AS-osteoprogenitor cells. WNT16 treatment activated the β-catenin functional motif, as confirmed by promoter activity assays (Fig. 3H). In contrast, WNT16 knockdown using siRNA attenuated the increased phos-AKT and CDKN1A proteins by H2O2 stimulation (Fig. 3J). WNT16 knockdown groups significantly reduced the production of SASP mRNA expression response to H2O2 stimulation (Fig. 3K). Taken together, WNT16 treatment facilitated cell senescence of AS-osteoprogenitor cells under H2O2 stimulation.

Fig. 3
figure 3

WNT16 treatment facilitated cell senescence of AS-osteoprogenitor cells under H2O2-stimulation. A Experimental design for cell senescence treated with WNT16. AS-osteoprogenitor cells were exposed to H2O2 for 2 h and then incubated with vehicle or WNT16 treatment for 70 h (n=5 per each group). Analysis of B SA-β-gal staining, C counting results of SA-β-gal-positive cells, D assessment of osteoprogenitor cell lysates using hydrogen peroxidase assay (n=5 per each group), E immunoblotting for protein level, F protein quantification of Fig. 3E, and G RT-qPCR for mRNA level (n=4 per each group). H SaOS2 cells were transfected with TOP or FOP-flash in the presence β-galactosidase for 24 h, incubated with WNT16 treatment for 24 h, and then analyzed by luciferase assays (n=4 per each group). I Experimental design for cell senescence treated with siRNA against WNT16. AS-osteoprogenitor cells were transfected with 100 nM WNT16 siRNA #2 for 70 h, exposed to H2O2 for 2 h, and followed by analysis of H immunoblotting for WNT16 protein level and I qPCR for mRNA level (n=5 per each group). Data are presented as the median and interquartile range. *p<0.05, **p<0.01. Representative images are shown. Scale bar = 200 μm


In this study, WNT16 was significantly increased in AS-osteoprogenitor cells and was observed in AS facet joint tissues. During osteoblast differentiation, WNT16 treatment in AS-osteoprogenitor decreased matrix maturation and mineralization accompanied by increased cell senescence, but not in control-osteoprogenitor. Moreover, WNT16 treatment strikingly facilitated cell senescence, p21 protein, and SASP expression in AS-osteoprogenitor cells under H2O2-induced stimulation, while WNT16 knockdown attenuated this senescence effect. Collectively, high WNT16 expression in AS-osteoprogenitor cells revealed that inhibiting bone mineralization and inducing cellular senescence of AS-osteoprogenitor, indicating that it might influence the bone loss in patients with AS.

In general, WNT3 is well known to play a positive role in bone formation through β-catenin activation. Li et al. showed that WNT molecule serum levels (WNT3a, WNT4, WNT5a, WNT7b, and WNT10) were increased dependent on the progression status of spinal ankylosis. In particular, the result showed that WNT3 and WNT5 were associated with hyper-osteoblastic activity and spinal ectopic new bone formation, and their expressions were in response to TNF stimulation [3]. In this paper, we showed that WNT3 expression was highly expressed in AS-osteoprogenitor cells (Fig. 1A, B). However, reports of WNT3 serum level in AS groups have been inconsistent [3, 29, 30], suggesting that the functional and pathophysiological changes need further investigation and a constructive approach.

Many studies have indicated that WNT16 transgenic mice show an increase in bone mineral density and mass whereas WNT16 knockout mice show a decrease in cortical and cancellous bone mass [14,15,16]. Here, we suggest that human WNT16 plays a paradoxical role in bone formation in AS. Consistent with our finding, DKK1 knockdown in osteoblast of WNT16 transgenic mice showed the reduction of calcified nodule formation relative to the control [31], indicating that the functional role of WNT16 in mineralization may be different depending upon DKK1 expression. Also, human meta-analyses have reported low levels of DKK1 in AS [32], and high WNT16 expression in AS is considered to lead to restrained bone mineralization.

Many studies have reported that WNT16 expression is strongly associated with bone mineral density, suggesting that WNT16 is a therapeutic target for bone loss-related diseases such as osteopenia or osteoporosis. While these data are mainly from murine models, there are a few human study results. It is a questionable issue whether WNT16 plays a functional role in humans. Since cell senescence in old mice exhibits a decline in bone mass accompanied by Osterix protein decrease [33], we analyzed whether WNT16 treatment affects Osterix expression. However, we found no change in Osterix mRNA expression (data now shown).

WNT16 genome-wide association study (GWAS) meta-analyses have shown that polymorphisms of the WNT16 locus are strongly associated with cortical bone thickness, BMD, and risk of osteoporosis fracture in humans [15, 34,35,36,37,38]. Therefore, these studies indicate that WNT16 is an important determinant of cortical bone mass and risk of fracture. As discussed above, the functional relevance of WNT16 in the regulation of the bone and mineralization has been confirmed in both in vitro and in vivo studies. Importantly, there is a lack of information on the functional role of WNT16 in AS-associated bone loss; therefore, further study on the relationship between bone loss and WNT16 locus polymorphisms in patients with AS is needed.

We measured WNT16 levels in collected serum and confirmed that there was no significant difference between the control and AS groups (data now shown). Our serum results are consistent with previous reports [3], but WNT16 expression was high in AS-osteoprogenitor cells and AS bone tissue compared to the control. Therefore, the difference in WNT16 level between serum and pathological tissue can be interpreted as the distinction between systemic and local disease.

A limitation of our paper is that we did not provide the direct regulatory mechanisms on WNT16-induced p21 expression in AS-osteoprogenitor cells. Second, we did not explore that how WNT16 increases oxidative stress and related protein expression in AS-osteoprogenitor cells. Third, further experiments are needed to explore the effect of WNT16 treatment on the reduction of bone formation-related genes.


Our present study shows that WNT16 is highly expressed in AS-osteoprogenitor cells and in destructive bone tissue. WNT16 treatment in AS-osteoprogenitor cells inhibits mineralization of osteoblast differentiation accompanied by promoting cell senescence. Therefore, elevated WNT16 might be a putative indicator for bone loss in AS.

Availability of data and materials

Microarray day used in this study are available from the corresponding author upon reasonable request.



Ankylosing spondylitis


Healthy controls


Alkaline phosphatase


Von Kossa




Osteoblast-specific element


Bone sialoprotein






Caspase 3


Senescence-associated secretory phenotype


Modified Stoke Ankylosing Spondylitis Spinal Score


C-reactive protein


Oxidative stress


Reactive oxygen species


Mesenchymal stem cells


  1. 1.

    Sieper J, Poddubnyy D. Axial spondyloarthritis. Lancet. 2017;390(10089):73–84.

    Article  Google Scholar 

  2. 2.

    Ranganathan V, Gracey E, Brown MA, Inman RD, Haroon N. Pathogenesis of ankylosing spondylitis - recent advances and future directions. Nat Rev Rheumatol. 2017;13(6):359–67.

    CAS  Article  Google Scholar 

  3. 3.

    Li X, Wang J, Zhan Z, Li S, Zheng Z, Wang T, et al. Inflammation intensity-dependent expression of osteoinductive Wnt proteins is critical for ectopic new bone formation in ankylosing spondylitis. Arthritis Rheumatol. 2018;70(7):1056–70.

    CAS  Article  Google Scholar 

  4. 4.

    Klingberg E, Lorentzon M, Mellstrom D, Geijer M, Gothlin J, Hilme E, et al. Osteoporosis in ankylosing spondylitis - prevalence, risk factors and methods of assessment. Arthritis Res Ther. 2012;14(3):R108.

    Article  Google Scholar 

  5. 5.

    van der Weijden MA, Claushuis TA, Nazari T, Lems WF, Dijkmans BA, van der Horst-Bruinsma IE. High prevalence of low bone mineral density in patients within 10 years of onset of ankylosing spondylitis: a systematic review. Clin Rheumatol. 2012;31(11):1529–35.

    Article  Google Scholar 

  6. 6.

    Nigil Haroon N, Szabo E, Raboud JM, McDonald-Blumer H, Fung L, Josse RG, et al. Alterations of bone mineral density, bone microarchitecture and strength in patients with ankylosing spondylitis: a cross-sectional study using high-resolution peripheral quantitative computerized tomography and finite element analysis. Arthritis Res Ther. 2015;17:377.

    Article  Google Scholar 

  7. 7.

    Inman RD. Axial spondyloarthritis: current advances, future challenges. Journal of Rheumatic Diseases. 2021;28(2):55–9.

    Article  Google Scholar 

  8. 8.

    Akgol G, Kamanli A, Ozgocmen S. Evidence for inflammation-induced bone loss in non-radiographic axial spondyloarthritis. Rheumatology (Oxford). 2014;53(3):497–501.

    Article  Google Scholar 

  9. 9.

    Davey-Ranasinghe N, Deodhar A. Osteoporosis and vertebral fractures in ankylosing spondylitis. Curr Opin Rheumatol. 2013;25(4):509–16.

    Article  Google Scholar 

  10. 10.

    Stal R, van Gaalen F, Sepriano A, Braun J, Reijnierse M, van den Berg R, et al. Facet joint ankylosis in r-axSpA: detection and 2-year progression on whole spine low-dose CT and comparison with syndesmophyte progression. Rheumatology (Oxford). 2020;59(12):3776–83.

    Article  Google Scholar 

  11. 11.

    Lories RJ, Luyten FP, de Vlam K. Progress in spondylarthritis. Mechanisms of new bone formation in spondyloarthritis. Arthritis Res Ther. 2009;11(2):221.

    Article  Google Scholar 

  12. 12.

    Xie W, Zhou L, Li S, Hui T, Chen D. Wnt/beta-catenin signaling plays a key role in the development of spondyloarthritis. Ann N Y Acad Sci. 2016;1364:25–31.

    CAS  Article  Google Scholar 

  13. 13.

    Lories RJ, Corr M, Lane NE. To Wnt or not to Wnt: the bone and joint health dilemma. Nat Rev Rheumatol. 2013;9(6):328–39.

    CAS  Article  Google Scholar 

  14. 14.

    Tornqvist AE, Grahnemo L, Nilsson KH, Funck-Brentano T, Ohlsson C, Moverare-Skrtic S. Wnt16 overexpression in osteoblasts increases the subchondral bone mass but has no impact on osteoarthritis in young adult female mice. Calcif Tissue Int. 2020;107(1):31–40.

    Article  Google Scholar 

  15. 15.

    Zheng HF, Tobias JH, Duncan E, Evans DM, Eriksson J, Paternoster L, et al. WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk. PLoS Genet. 2012;8(7):e1002745.

    CAS  Article  Google Scholar 

  16. 16.

    Moverare-Skrtic S, Henning P, Liu X, Nagano K, Saito H, Borjesson AE, et al. Osteoblast-derived WNT16 represses osteoclastogenesis and prevents cortical bone fragility fractures. Nat Med. 2014;20(11):1279–88.

    CAS  Article  Google Scholar 

  17. 17.

    Ghozlani I, Ghazi M, Nouijai A, Mounach A, Rezqi A, Achemlal L, et al. Prevalence and risk factors of osteoporosis and vertebral fractures in patients with ankylosing spondylitis. Bone. 2009;44(5):772–6.

    CAS  Article  Google Scholar 

  18. 18.

    Kim JW, Chung MK, Lee J, Kwok SK, Kim WU, Park SH, et al. Low bone mineral density of vertebral lateral projections can predict spinal radiographic damage in patients with ankylosing spondylitis. Clin Rheumatol. 2019;38(12):3567–74.

    Article  Google Scholar 

  19. 19.

    Ye G, Xie Z, Zeng H, Wang P, Li J, Zheng G, et al. Oxidative stress-mediated mitochondrial dysfunction facilitates mesenchymal stem cell senescence in ankylosing spondylitis. Cell Death Dis. 2020;11(9):775.

    CAS  Article  Google Scholar 

  20. 20.

    Wang L, Gao L, Jin D, Wang P, Yang B, Deng W, et al. The relationship of bone mineral density to oxidant/antioxidant status and inflammatory and bone turnover markers in a multicenter cross-sectional study of young men with ankylosing spondylitis. Calcif Tissue Int. 2015;97(1):12–22.

    CAS  Article  Google Scholar 

  21. 21.

    Yazici C, Kose K, Calis M, Kuzuguden S, Kirnap M. Protein oxidation status in patients with ankylosing spondylitis. Rheumatology (Oxford). 2004;43(10):1235–9.

    CAS  Article  Google Scholar 

  22. 22.

    Karakoc M, Altindag O, Keles H, Soran N, Selek S. Serum oxidative-antioxidative status in patients with ankylosing spondilitis. Rheumatol Int. 2007;27(12):1131–4.

    CAS  Article  Google Scholar 

  23. 23.

    van der Linden S, Valkenburg HA, Cats A. Evaluation of diagnostic criteria for ankylosing spondylitis. A proposal for modification of the New York criteria. Arthritis Rheum. 1984;27(4):361–8.

    Article  Google Scholar 

  24. 24.

    Kiltz U, Braun J. Assessments of functioning in patients with axial spondyloarthritis. Journal of Rheumatic Diseases. 2020;27(1).

  25. 25.

    Jo S, Wang SE, Lee YL, Kang S, Lee B, Han J, et al. IL-17A induces osteoblast differentiation by activating JAK2/STAT3 in ankylosing spondylitis. Arthritis Res Ther. 2018;20(1):115.

    Article  Google Scholar 

  26. 26.

    Jo S, Lee EJ, Nam B, Kang J, Lee S, Youn J, et al. Effects of dihydrotestosterone on osteoblast activity in curdlan-administered SKG mice and osteoprogenitor cells in patients with ankylosing spondylitis. Arthritis Res Ther. 2020;22(1):121.

    CAS  Article  Google Scholar 

  27. 27.

    Park P-R, Jo S, Jin S-H, Kim T-J. MicroRNA-10b plays a role in bone formation by suppressing interleukin-22 in ankylosing spondylitis. Journal of Rheumatic Diseases. 2020;27(1).

  28. 28.

    Jo S, Won EJ, Kim MJ, Lee YJ, Jin SH, Park PR, et al. STAT3 phosphorylation inhibition for treating inflammation and new bone formation in ankylosing spondylitis. Rheumatology (Oxford). 2020.

  29. 29.

    Iaremenko O, Shynkaruk I, Fedkov D, Iaremenko K, Petelytska L. Bone turnover biomarkers, disease activity, and MRI changes of sacroiliac joints in patients with spondyloarthritis. Rheumatol Int. 2020;40(12):2057–63.

    CAS  Article  Google Scholar 

  30. 30.

    Klingberg E, Nurkkala M, Carlsten H, Forsblad-d'Elia H. Biomarkers of bone metabolism in ankylosing spondylitis in relation to osteoproliferation and osteoporosis. J Rheumatol. 2014;41(7):1349–56.

    Article  Google Scholar 

  31. 31.

    Jiang Z, Von den Hoff JW, Torensma R, Meng L, Bian Z. Wnt16 is involved in intramembranous ossification and suppresses osteoblast differentiation through the Wnt/beta-catenin pathway. J Cell Physiol. 2014;229(3):384–92.

    CAS  Article  Google Scholar 

  32. 32.

    Zhang L, Ouyang H, Xie Z, Liang ZH, Wu XW. Serum DKK-1 level in the development of ankylosing spondylitis and rheumatic arthritis: a meta-analysis. Exp Mol Med. 2016;48:e228.

    CAS  Article  Google Scholar 

  33. 33.

    Kim HN, Chang J, Shao L, Han L, Iyer S, Manolagas SC, et al. DNA damage and senescence in osteoprogenitors expressing Osx1 may cause their decrease with age. Aging Cell. 2017;16(4):693–703.

    CAS  Article  Google Scholar 

  34. 34.

    Garcia-Ibarbia C, Perez-Nunez MI, Olmos JM, Valero C, Perez-Aguilar MD, Hernandez JL, et al. Missense polymorphisms of the WNT16 gene are associated with bone mass, hip geometry and fractures. Osteoporos Int. 2013;24(9):2449–54.

    CAS  Article  Google Scholar 

  35. 35.

    Hendrickx G, Boudin E, Fijalkowski I, Nielsen TL, Andersen M, Brixen K, et al. Variation in the Kozak sequence of WNT16 results in an increased translation and is associated with osteoporosis related parameters. Bone. 2014;59:57–65.

    CAS  Article  Google Scholar 

  36. 36.

    Kemp JP, Medina-Gomez C, Estrada K, St Pourcain B, Heppe DH, Warrington NM, et al. Phenotypic dissection of bone mineral density reveals skeletal site specificity and facilitates the identification of novel loci in the genetic regulation of bone mass attainment. PLoS Genet. 2014;10(6):e1004423.

    Article  Google Scholar 

  37. 37.

    Koller DL, Zheng HF, Karasik D, Yerges-Armstrong L, Liu CT, McGuigan F, et al. Meta-analysis of genome-wide studies identifies WNT16 and ESR1 SNPs associated with bone mineral density in premenopausal women. J Bone Miner Res. 2013;28(3):547–58.

    CAS  Article  Google Scholar 

  38. 38.

    Medina-Gomez C, Kemp JP, Estrada K, Eriksson J, Liu J, Reppe S, et al. Meta-analysis of genome-wide scans for total body BMD in children and adults reveals allelic heterogeneity and age-specific effects at the WNT16 locus. PLoS Genet. 2012;8(7):e1002718.

    CAS  Article  Google Scholar 

Download references




This work was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2019R1A2C2004214, 2020R1A2C1102386, and 2021R1A6A1A03038899).

Author information




SJ designed and performed the experiments. SW assisted with the laboratory work. M-AJ and HK analyzed microarray data analysis. BL and T-JK evaluated the human sera. Y-SP provided and evaluated the human bone samples. SJ and T-HK wrote the manuscript. T-HK conceived the project and was responsible for the overall design and oversight of the project. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Tae-Hwan Kim.

Ethics declarations

Ethics approval and consent to participate

The present study was conducted in accordance with the Helsinki Declaration and approved by the Institutional Review Board of Hanyang University, Seoul (IRB 2014-05-002), and Guri Hospital (IRB 2014-05-001). The written informed consent was obtained from all subjects.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Figure 1

. The overall image of a well for 96 well plates and its enlarged changes by Nikon microscopy.

Additional file 2: Supplementary Figure 2

. The result of siRNA efficiency test. (A) Total four human osteoblasts cell lines (MG63, FOB, SaOS2, U2OS) were extracted total RNA and analyzed by RT-PCR to test WNT16 level. (B) FOB cell were transfected with five siRNA oligo against human WNT16 and analyzed by RT-PCR for siRNA efficiency.

Additional file 3: Supplementary Figure 3

. Analysis of WNT16 level in collected sera. The human WNT16 levels were measures in the collected sera (HC, n=3; AS, n=36) using ELISA. N.S., not significant.

Additional file 4: Supplementary Figure 4

. WNT16 treatment inhibited ALP activity in AS-osteoprogenitors. AS-osteoprogenitors were plated in 96 wells plate and treated with 0, 25, or 50 ng/ml WNT16, and then analyzed as indicated by (A) WST assay for cell proliferation. (B) ALP and COL staining, and (C) ALP activity. Data are presented as the mean ± SD (n=3). *p<0.05.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Jo, S., Weon, S., Nam, B. et al. WNT16 elevation induced cell senescence of osteoblasts in ankylosing spondylitis. Arthritis Res Ther 23, 301 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: