Int J Biol Sci 2023; 19(2):610-624. doi:10.7150/ijbs.78785 This issue Cite
Research Paper
1. Department of Plastic Surgery, Affiliated Friendship Plastic Surgery Hospital of Nanjing Medical University, Nanjing Medical University, Nanjing, China.
2. The Research Center for Bone and Stem Cells, Department of Anatomy, Histology and Embryology, Nanjing Medical University, Nanjing, China.
3. Center for Experimental Medicine, The Third Xiangya Hospital, Central South University, Changsha, China.
4. Calcium Research Laboratory, McGill University Health Centre and Department of Medicine, McGill University, Montreal, Quebec H4A 3J1, Canada.
*These authors contribute equally to this work.
Emerging observational data suggest that vitamin D deficiency is associated with the onset and progression of knee osteoarthritis (OA). However, the relationship between vitamin D level and OA and the role of vitamin D supplementation in the prevention of knee OA are controversial. To address these issues, we analyzed the articular cartilage phenotype of 6- and 12-month-old wild-type and 1α(OH)ase-/- mice and found that 1,25(OH)2D deficiency accelerated the development of age-related spontaneous knee OA, including cartilage surface destruction, cartilage erosion, proteoglycan loss and cytopenia, increased OARSI score, collagen X and Mmp13 positive chondrocytes, and increased chondrocyte senescence with senescence-associated secretory phenotype (SASP). 1,25(OH)2D3 supplementation rescued all knee OA phenotypes of 1α(OH)ase-/- mice in vivo, and 1,25(OH)2D3 rescued IL-1β-induced chondrocyte OA phenotypes in vitro, including decreased chondrocyte proliferation and cartilage matrix protein synthesis, and increased oxidative stress and cell senescence. We also demonstrated that VDR was expressed in mouse articular chondrocytes, and that VDR knockout mice exhibited knee OA phenotypes. Furthermore, we demonstrated that the down-regulation of Sirt1 in articular chondrocytes of 1α(OH)ase-/- mice was corrected by supplementing 1,25(OH)2D3 or overexpression of Sirt1 in mesenchymal stem cells (MSCs) and 1,25(OH)2D3 up-regulated Sirt1 through VDR mediated transcription. Finally, we demonstrated that overexpression of Sirt1 in MSCs rescued knee OA phenotypes in 1α(OH)ase-/- mice. Thus, we conclude that 1,25(OH)2D3, via VDR-mediated gene transcription, plays a key role in preventing the onset of aging-related knee OA in mouse models by up-regulating Sirt1, an aging-related gene that promotes articular chondrocyte proliferation and extracellular matrix protein synthesis, and inhibits senescence and SASP.
Keywords: Vitamin D deficiency, vitamin D supplementation, vitamin D receptor, osteoarthritis, Sirt1
Osteoarthritis (OA) characterized by pain, progressive loss of articular cartilage and structural changes is extremely common in aged individuals [1]. Vitamin D deficiency has also been a common health problem. Over a billion people worldwide are vitamin D deficient or insufficient [2]. Vitamin D deficiency is closely associated with common chronic diseases such as metabolic bone disorders, diabetes and cardiovascular diseases [3]. Emerging observational data suggest that vitamin D insufficiency is associated with the onset and progression of knee OA [4-6]. However, studies examining the prevention of OA after vitamin D supplementation have been inconclusive [7]. A pilot randomized controlled trial in India revealed that 12-months of vitamin D treatment had small but statistically significant clinical benefits on pain and function in patients with knee OA compared with placebo [8]. In contrast, a subsequent randomized controlled trial in the United States showed that vitamin D supplementation over 24 months did not reduce knee pain or cartilage volume loss [9]. Subsequently, data from randomized controlled trials in Australia [10] and the UK [11] did not find a significant clinical benefit of vitamin D supplementation for knee OA. Well-designed animal studies are required to investigate the role of vitamin D in the initiation and progression of knee OA; however, adequate cell and animal models are lacking.
Vitamin D, mainly synthesized in skin tissue, is firstly metabolized by liver 25-hydroxylase, to 25-hydroxyvitamin D [25(OH)D], which is then converted to 1,25(OH)2D by the renal or extrarenal 1α-hydroxylase [1α(OH)ase]. 1,25(OH)2D then exerts its principal actions by binding to the vitamin D receptor (VDR). Although 1,25(OH)2D deficient [1α(OH)ase-/-] [12] and VDR deficient (VDR-/-) [13] mouse models have been generated over 2 decades and have provided considerable insight into the regulation of mineral and skeletal physiology by 1,25(OH)2D [14, 15], whether OA phenotypes occur in long bone joints of 1α(OH)ase-/- and VDR-/- mice has not been reported. We have previously reported that 1,25(OH)2D deficiency may play a role in the pathogenesis of temporomandibular joint OA [16]. However, it is unknown whether 1,25(OH)2D deficiency can induce the occurrence of knee OA via VDR and what the key downstream targets of 1,25(OH)2D might be in preventing the occurrence of knee OA.
Sirt1 is a mammalian homologue of Sir2, which is a histone deacetylase that regulates gene expression and protein function by deacetylating lysine residues in histone and non-histone proteins, including p65 and forkhead box protein O. Sirt1 exerts anti-aging and other protective effects through regulation of mitochondrial biogenesis, oxidative stress, and inflammation-pathways that are known to drive OA [17]. Sirt1 has a chondroprotective role by suppressing interleukin-1β (IL-1β)- and tumor necrosis factor α-induced expression of cartilage-degrading enzymes. In mice, Sirt1 haploinsufficiency [18] or Sirt1-inactivating mutations [19] results in delayed growth and increased spontaneous OA. Chondrocyte-specific deletion of Sirt1 resulted in increased severity of OA with aging and following joint injury [20]. These findings in mice are also consistent with the reduction of Sirt1 expression in human OA cartilage. However, it is unclear whether Sirt1 is a key downstream target for 1,25(OH)2D in preventing the occurrence of knee OA.
Previous studies have shown that cells isolated from the superficial zone of mouse and human articular cartilage express mesenchymal stem cell (MSC) markers and exhibit MSC characteristics, and MSCs have been a therapeutic option for the prevention of osteoarthritis [21-23]. Indeed, accumulation of senescent MSCs in joints contributes to the development of osteoarthritis, and recent studies have reported that the transplantation of MSCs or the rejuvenation of endogenous senescent MSCs could be a therapeutic option for osteoarthritis [21-23]. Prx1-expressing periosteal cells can differentiate into chondrocytes and osteoblasts in vitro and Prx1 has been used to study chondrogenesis and/or skeletal limb formation. Therefore, periosteal Prx1-expressing cells are most likely chondro-osteoprogenitor cells [24]. Recently, we generated a transgenic mouse model by overexpressing Sirt1 in MSCs using Prx1 as a promoter (Sirt1Tg) and found that overexpression of Sirt1 in MSCs protected against bone loss induced by Bmi1 deficiency in long bone and mandibles [25, 26] and improved the mandibular bone loss induced by 1,25(OH)2D deficiency [27] without altering mineral levels.
In our current study, to investigate the mechanism of 1,25(OH)2D deficiency in the pathogenesis of knee OA, we firstly compared the articular cartilage phenotype of 6- and 12-month-old wild-type and 1α(OH)ase-/- mice fed a high calcium/phosphate (rescue) diet to assess whether 1,25(OH)2D deficiency induced age-related spontaneous knee OA. We then treated 1α(OH)ase-/- mice or IL-1β-treated human articular chondrocytes with exogenous 1,25(OH)2D3 to assess whether treatment with 1,25(OH)2D3 rescued the knee OA phenotype caused by 1,25(OH)2D deficiency, and also the IL-1β-induced reduction of chondrocyte proliferation and the increases of oxidative stress and cellular senescence. We also compared the articular chondrocyte and cartilage phenotypes of 6-month-old wild-type and VDR-/- mice to assess whether VDR deficiency induced chondrocyte senescence and a knee OA phenotype. The mechanisms of 1,25(OH)2D3 in regulating Sirt1 via VDR were examined in chondrocytes using molecular techniques. Finally, we generated 1α(OH)ase-/- mice with Sirt1 overexpression in MSCs, and compared their knee phenotypes with 1α(OH)ase-/- mice on the rescue diet and wild-type littermates to assess whether overexpression of Sirt1 in MSCs could prevent 1,25(OH)2D deficiency-induced development of knee OA.
The generation and characterization of 1α(OH)ase-/- mice and vitamin D receptor knockout mice (VDR-/-) were previously described [28]. After weaning, male 1α(OH)ase-/- mice and wild-type (WT) littermates were given a rescue diet (RD) containing 2.0% calcium, 1.25% phosphorus and 20% lactose and mice aged 6 and 12 months were sacrificed for joint phenotype analysis. WT and 1α(OH)ase-/- mice on a normal diet were subcutaneously injected, respectively, with vehicle and 1,25(OH)2D3 (1 μg/kg B.W.) three times per week, from the time of weaning until analysis at the age of 12 months. WT and VDR-/- littermates were given the rescue diet after weaning and sacrificed at the age of 6 months for joint phenotype analysis. The transgenic mouse model of Sirt1 overexpression driven by Prx1 gene, which represents the mesenchymal lineage (Sirt1Tg) was generated as previously described [29]. Meanwhile, 1α(OH)ase+/- mice and Sirt1Tg mice were mated to produce offspring heterozygous at both loci (1α(OH)ase+/-Sirt1Tg), which were then mated to generate 1α(OH)ase-/-Sirt1Tg mice. After weaning, WT, Sirt1Tg, 1α(OH)ase-/- and 1α(OH)ase-/-Sirt1Tg mice were given the rescue diet and used for phenotype analysis at the age of 12 months. All animal experiments were performed in compliance with the guidelines approved by the Institutional Animal Care and Use Committee of Nanjing Medical University.
Serum calcium and phosphorus levels were analyzed by an autoanalyzer (Beckman Synchron 67; Beckman Instruments). Serum intact PTH was measured using an ELISA kit (Immutopics, Inc., San Clemente, CA). The serum levels of 25(OH)D3 (Nanjing Jiancheng Bioengineering Institute) and 1,25(OH)2D (Nanjing Jiancheng Bioengineering Institute) were measured as previously described [30].
At the time of euthanasia, mouse knee joints were dissected free of soft tissue and fixed in PLP fixative buffer for 48 hours at 4°C, decalcified with 10% EDTA and embedded in paraffin, after which 5-μm sections were stained for safranin O/fast green. The histological OA severity was evaluated using the Osteoarthritis Research Society International (OARSI) scoring system (grade 0-6). For each joint, OA grade was examined by a single observer under blinded conditions.
Knee joints of indicated groups were collected and examined using μCT and subchondral bone volume were examined as previously described [31].
Human articular chondrocytes were purchased from Procell (Wuhan, China). Primary mouse chondrocytes isolated from articular cartilage of newborn mice were cultured as described previously. Briefly, cartilage tissue was cut into 1 mm3 pieces, washed three times with PBS and digested with 0.2% type II collagenase (Worthington Biochemical, 4176) in high-glucose DMEM-F12 (Gibco) overnight at 37 °C. After digestion, the filtrate was passed through a 70-µm strainer, and cells were washed three times with PBS, and seeded in DMEM-F12 complete medium containing 100 U/ml penicillin, 100 µg/ml streptomycin and 10% FBS (Gibco) at 37 °C with 5% CO2. For 3D cultures of human and mouse chondrocytes, 400,000 cells were seeded in 96-well MicroWell round-bottom plates (Thermo Fisher Scientific). Pellets were obtained after centrifugation at 150g for 10 min. The chondrocyte pellets were maintained at 37 °C with 5% CO2 in 200µl of chondrocyte growth medium consisting of high-glucose DMEM-F12, 1% Pen/Strep and 10% FBS. Recombinant Human Interleukin-1β (IL-1β) (200-01B, PeproTech), at 10ng/ml, was added as an inflammatory stimulus, and 100nM 1,25(OH)2D3 and control vehicle were added into indicated wells. Pellets were harvested after 21 days for evaluation. Medium was changed every other day until the end of the experiment. SA-β-gal staining, crystal violet staining, EdU assay and dihydroethidium (DHE) staining were performed as previously described [32]. Ex527, a selective Sirt1 inhibitor was purchased from MedChemExpress (MCE).
For immunofluorescence staining, cultured cells and sections were incubated with rabbit anti- VDR (Abcam, ab3508) at 4 °C overnight followed by using Alexa Fluor Plus 594 goat anti-rabbit IgG (H+L) secondary antibody (Invitrogen, A32740, 1:500 dilution) to detect immunoreactivity. For Immunohistochemistry staining, dewaxed and rehydrated paraffin-embedded sections were incubated with 6% hydrogen peroxide to block endogenous peroxidase activity and then washed in PBS (pH 7.6). The slides were then incubated with the primary antibodies, including rabbit anti-VDR (Abcam, ab3508), rabbit anti-β-gal (Abcam, ab616), rabbit anti-Collagen II (Abcam, ab34712), mouse anti-Collagen X (Abcam, ab49945), rabbit anti-Mmp13 (Proteintech, 18165-1-AP) and mouse anti-p16 (Santa Cruz, sc-1661) overnight at 4 °C. After rinsing with PBS for 15 min, sections were incubated with secondary antibody (biotinylated goat anti-rabbit IgG and goat anti-mouse IgG, Sigma). Sections were then washed and incubated with the Vectastain Elite ABC reagent (Vector Laboratories) for 30 minutes. Staining was done using 3,3-diaminobenzidine (2.5 mg/ml) followed by counterstaining with Mayer's hematoxylin.
Chromatin immunoprecipitation (ChIP) was performed using the ChIP kit (Millipore, USA) as previously described [33]. Briefly, human chondrocytes were treated with vehicle or 10-7 M 1,25(OH)2D3 for 3 hours, and then cell samples were subjected to immunoprecipitation using either a control IgG or rabbit anti-VDR antibody (Abcam, ab3508). The co-precipitated chromatin was determined by qPCR for the presence of human Sirt1 promoter sequence using Sirt1 sense 5′-TTAGAGTGGCTTACAGGC-3′ and antisense 5′-ACATCTTCTGGCTTCCTT-3′ primers.
To generate Sirt1 promoter-activated luciferase reporter, -1442 to +97 bp of Sirt1 promoter and the deletion mutant were directly synthesized and cloned into pGL3-basic. Luciferase reporter assay was performed as previously described [33]. Briefly, human chondrocytes (passage 4) in 12-well plates were transfected with Sirt1-promoter or Sirt1-promoter mutant luciferase reporter plasmid. Twelve hours later, the medium was changed and the indicated concentration of vehicle or 1,25(OH)2D3 were added. After another 36 hours, luciferase activity was measured using the Dual-Luciferase assay kit (Promega). pRL-TK was co-transfected to normalize transfection efficiency.
Cell or tissue lysates were extracted for loading into 10% SDS-PAGE gels and immunoblotting was performed as previously described [34]. Primary antibodies, including rabbit anti-Cyp27b1/1α(OH)ase (Abcam, ab206655), rabbit anti-VDR (Abcam, ab3508), rabbit anti-Sirt1 (Millipore, 07-131), rabbit anti-p16 (Proteintech, 10883-1-AP) and rabbit anti-β-actin (Cell Signaling Technology, 8457S) were used for immunoblotting. Immunoreactive bands were visualized with ECL chemiluminescence (Bio-Rad) and analyzed by Image J.
RNA isolation from cells and joint tissue, cDNA synthesis and quantitative real-time PCR (qRT-PCR) were performed as previously described [32]. The PCR primer sequences used in this study are shown in the Supplementary Table 1.
Measured data were described as mean ± SD. The statistical analyses were performed using GraphPad Prism (Version 8.0). The number of animals used in the experiments in vivo was shown in indicated figure legends. Two-tailed Student's t-test was used to compare differences between groups. For multiple comparisons, one-way ANOVA analysis of variance followed by Tukey's post-hoc test was used. P values <0.05, <0.01, and <0.001 were considered statistically significant (*, **, ***).
To determine whether 1,25(OH)2D deficiency accelerated age-related knee OA, we compared the articular cartilage phenotype of 6- and 12-month-old wild-type and 1,25(OH)2D deficient [1α(OH)ase-/-] mice fed a high calcium/phosphate (rescue) diet using safranin O staining. The rescue diet normalized serum calcium, phosphorus and intact PTH levels, but did not alter increased 25(OH)D levels, while 1,25(OH)2D was undetectable in 1α(OH)ase-/- mice (Figure 1A). Three-month-old 1α(OH)ase-/- mice on the normal diet were 28% lighter than wild-type littermates (the average lifespan of 1α(OH)ase-/- mice on a normal diet was 3 months) (Figure S1A), while 6- and 12-month-old 1α(OH)ase-/- mice on the rescue diet were 10% and 12% lighter than control littermates, respectively. (Figure S1B). The cartilage superficial destruction, cartilage erosion, proteoglycan loss and cytopenia increased progressively with increasing age in both wild-type mice and 1α(OH)ase-/- mice on the rescue diet; however, these alterations were more dramatic in 1α(OH)ase-/- mice on the rescue diet compared with age-matched wild-type mice (Figure 1B). We then assessed the OARSI grade based on the sections stained with Safranin O/Fast green using the OARSI histological scoring system. The OARSI scores were progressively elevated with increasing age in both wild-type mice and 1α(OH)ase-/- mice on the rescue diet, but were significantly higher in 1α(OH)ase-/- mice than in age-matched wild-type mice (Figure 1C). Our previous study reported that 1,25(OH)2D deficiency induced aging-related bone loss [31]. To test this indirect effect induced by 1,25(OH)2D deficiency, we examined the subchondral bone volume using micro-CT, and found that subchondral bone volume was not significantly altered in 6- and 12-month-old 1,25(OH)2D-deficient mice as compared with wild-type controls (Fig. 1D, E). To determine whether the age-related OA-like phenotype occurring in wild-type mice was associated with down-regulation of 1α(OH)ase in articular cartilage, we examined the protein expression level of 1α(OH)ase in articular cartilage isolated from 6- and 12-month-old wild-type mice and found that the protein levels of 1α(OH)ase in articular cartilage were significantly down-regulated in 12-month-old wild-type mice compared with 6-month-old wild-type mice (Fig. 1F, G). These results demonstrate that the enzyme which generates active 1,25(OH)2D declines in articular cartilage with age and that 1,25(OH)2D deficiency accelerates the age-related development of knee OA.
1,25(OH)2D deficiency accelerates the development of age-related spontaneous knee OA. (A) Serum calcium, phosphorus, intact PTH, 25(OH)D and 1,25(OH)2D3 levels in 12-month-old wild-type (n=8) and 1α(OH)ase-/- (n=8) mice on a normal diet (ND) or a rescue diet (RD). N.S.= not significant; N.D.= not detected. **: p< 0.01, ***: p<0.001. Representative sectioned images of articular cartilage stained for (B) safranin O and (C) the quantification of OARSI grade in 6- and 12-month-old WT and 1α(OH)ase-/- mice on the rescue diet. n=8 mice per group. (D) Representative μCT images of knee joints and (E) the quantification of subchondral bone volume (BV/TV, %) in indicated groups of mice. (F) Western blot detection of 1α(OH)ase in articular cartilage isolated from 6- and 12-month-old WT mice (n=5 mice in each group). Kidney tissue from 6-month-old WT mice (n=2) were used as a positive control. (G) Quantification of 1α(OH)ase relative protein levels. *: p< 0.05, **: p< 0.01, ***: p<0.001.
Supplementation of 1,25(OH)2D3 rescues the knee OA phenotype caused by 1,25(OH)2D deficiency. Representative images of articular cartilage stained for (A) safranin O and immunostained for (C) collagen X, (E) Mmp13 and (G) p16 in 12-month-old wild-type mice on the rescue diet (RD) (n=4), 1α(OH)ase-/- mice on the RD (n=4), and 1α(OH)ase-/- mice treated with 1,25(OH)2D3 (n=5). Quantification of (B) OARSI grade, (D) collagen X+ cells, (F) Mmp13+ cells and (H) p16+ cells. (I) Quantification of mRNA levels for SASP including p16, p21, IL-1β, IL-6 and Mmp13. **: p< 0.01, ***: p<0.001.
To assess whether supplementation of exogenous 1,25(OH)2D3 can prevent knee OA caused by 1,25(OH)2D deficiency, 1α(OH)ase-/- mice were either fed the rescue diet after weaning or were injected with 1,25(OH)2D3 subcutaneously (1 μg/kg, three times a week) after weaning until 12 months of age; wild-type mice were fed the rescue diet as the control. Serum calcium, phosphorus, and PTH levels were normalized in 1α(OH)ase-/- mice fed on the rescue diet [33]. Knee phenotypes were analyzed in 12-month-old mice. The cartilage superficial destruction, cartilage erosion, proteoglycan loss and cytopenia in 1α(OH)ase-/- mice were largely rescued by supplementation of exogenous 1,25(OH)2D3 (Figure 2A). OARSI scores, the percentages of collagen Ⅹ+, Mmp13+, and p16+ articular chondrocytes, and the mRNA levels of p16, p21, IL-1β, IL-6 and Mmp13 were all significantly increased in 1α(OH)ase-/- mice, and were reduced by supplementation of exogenous 1,25(OH)2D3 (Figure 2B-I). These results suggest that supplementation with exogenous 1,25(OH)2D3 rescued knee OA phenotype induced by 1,25(OH)2D deficiency by inhibiting chondrocyte senescence and SASP.
1,25(OH)2D3 rescues the reduction of human chondrocyte proliferation and the increases of oxidative stress and cellular senescence induced by IL-1β. Human articular chondrocytes were cultured in vitro in the presence or absence of IL-1β with or without 1,25(OH)2D3 and stained (A) with crystal violet, (D) immunocytochemically for EdU, (F) with DHE and (H) cytochemically for SA-β-gal. n = 3 wells per condition. (B) Clonal expansion analysis. (C) CCK8 assay showing cell viability. (E) EdU-positive cells. (G) ROS relative levels (relative intensity of DHE fluorescence analyzed by Image J). (I) SA-β-gal positive cells. (J) Western blot detection of the cellular senescence marker p16 and (K) the quantification of protein levels of p16. n = 3 wells per condition. *: p< 0.05, **: p< 0.01, ***: p<0.001.
To determine the potential mechanism underlying the protective effects of 1,25(OH)2D3 on knee OA phenotypes, we studied the impact of in vitro application of IL-1β, a critical factor in the onset of OA [35-37]. Human articular chondrocytes were cultured in vitro in the presence or absence of IL-1β with or without 1,25(OH)2D3 and then cell proliferation, senescence and reactive oxygen species (ROS) levels were analyzed. ROS regulate cellular homeostasis and act as prime modulators of cellular dysfunction contributing to disease pathophysiology. We firstly found that the protein expression of VDR was significantly reduced in human articular chondrocytes treated with IL-1β, and this reduction was markedly rescued upon 1,25(OH)2D3 treatment (Figure S2A-B). In addition, results revealed that the density and viability of chondrocytes and percentage of EdU-positive chondrocytes were significantly reduced in IL-1β + vehicle-treated cultures relative to vehicle-treated cultures, whereas they were insignificantly altered in 1,25(OH)2D3-treated cultures relative to vehicle-treated cultures. However, EdU-positive chondrocytes were significantly increased in IL-1β + 1,25(OH)2D3-treated cultures relative to IL-1β + vehicle-treated cultures (Figure 3A-E). Furthermore, IL-1β was reported to induce an OA cell phenotype in vitro, that includes increased oxidative stress and cellular senescence [35]. Here we found that ROS levels, the percentage of SA-β-gal positive chondrocytes and protein levels of p16 were significantly increased in IL-1β + vehicle-treated cultures relative to vehicle-treated cultures, whereas they were insignificantly altered in 1,25(OH)2D3-treated cultures relative to vehicle-treated cultures. However, these parameters of oxidative stress and cellular senescence were significantly decreased in IL-1β + 1,25(OH)2D3-treated cultures relative to IL-1β + vehicle-treated cultures (Figure 3F-K). These results demonstrated that 1,25(OH)2D3 could rescue the reduction of chondrocyte proliferation and the increase of oxidative stress and cellular senescence induced by IL-1β.
We next assessed whether the development of knee OA induced by 1,25(OH)2D deficiency was mediated through VDR. We firstly confirmed the protein expression of VDR in articular chondrocytes by immunohistochemistry and Western blots in WT mouse articular tissues (Figure 4A) and in BMSCs (Figure 4B). To determine whether VDR in cartilage chondrocytes can respond to exogenous 1,25(OH)2D3, articular chondrocytes were treated with 10-8 M 1,25(OH)2D3, and the mRNA expression levels of VDR and CYP24A1, a known downstream target gene of VDR, were examined using RT-PCR; nuclear localization of VDR was also examined by immunofluorescence. The results showed that the mRNA levels of VDR and CYP24A1 were clearly increased, and VDR nuclear localization was also clearly increased in wild-type chondrocytes, but not in VDR-/- cells when they were treated with 1,25(OH)2D3 (Figure 4C, D). To determine whether VDR deficiency induces chondrocyte senescence and a knee OA phenotype, the phenotypes of micromass pellets from micromass cultures of articular chondrocytes derived from wild-type and VDR-/- mice, and the knee phenotypes from 6-month-old wild-type and VDR-/- mice on the rescue diet were analyzed. We found that serum calcium, phosphorus, and PTH levels were normalized in VDR-/- mice on the rescue diet relative to wild-type mice (Figs. S3A-C); however, safranin O and type II collagen, and the mRNA expression levels of col2a1 and aggrecan were decreased, whereas Mmp13+ and p16+ chondrocytes and their mRNA expression levels were increased in VDR-/- articular chondrocytes compared with wild-type ones (Figure 4E-I). Cartilage erosion, proteoglycan loss and cytopenia were observed in VDR-/- mice relative to their wild-type littermates (Figure 4J). The percentages of Mmp13+, Collagen Ⅹ+, and p16+ chondrocytes, and the mRNA expression levels of p16, IL-1β, IL-6, Mmp13, were all significantly increased, whereas the mRNA expression levels of aggrecan and col2a1 were markedly decreased in VDR-/- mice relative to their wild-type littermates (Figure 4K-O). These results demonstrated that VDR deficiency also induced knee OA by stimulating chondrocyte senescence and SASP, indicating that the development of knee OA induced by 1,25(OH)2D was mediated though VDR.
To uncover the molecular mechanism underlying how 1,25(OH)2D3 inhibits IL-1β-induced oxidative stress and cellular senescence, we treated human articular chondrocytes with vehicle or 1,25(OH)2D3 in the presence of IL-1β for 72 hours; subsequent KEGG pathways enrichment results showed that 1,25(OH)2D3 significantly regulated pathways such as “Cellular senescence” and “FoxO signaling pathway”, which are target pathways of Sirt1 (Figure S4A). To determine whether 1,25(OH)2D3 suppressed IL-1β-induced chondrocyte senescence by transcriptionally up-regulating Sirt1 in a VDR dependent manner, we examined the expression of Sirt1 in articular cartilage of WT and 1α(OH)ase-/- mice with and without applications of supplementary 1,25(OH)2D3. We found that exogenous 1,25(OH)2D3 rescued 1,25(OH)2D deficiency-induced reduction of Sirt1-positive chondrocytes (Figure S4B, C). We also found that the protein expression level of Sirt1 was significantly down-regulated in articular chondrocytes isolated from VDR-/- mice relative to those isolated from wild-type mice (Figure 5A). Human articular chondrocytes were treated with vehicle- and 1,25(OH)2D3 at indicated time points (0, 12, 24, and 48 hours for protein analyses, or 0, 3, 6 and 12 hours for mRNA analyses) in the presence of IL-1β. The results showed that 1,25(OH)2D3 up-regulated the protein and mRNA expression of Sirt1 in a time dependent manner (Figure 5B & C). To further determine whether 1,25(OH)2D3 regulates Sirt1 via the VDR at a transcriptional level, a VDRE-like sequence located upstream of the Sirt1 gene (http://jaspar.genereg.net/) was identified (Figure 5D, upper). The ChIP-PCR assay showed that VDR could directly bind to the Sirt1 promoter at the predicted binding site (Figure 5D, lower). Luciferase reporter assays showed that 1,25(OH)2D3 increased luciferase activity significantly in human articular chondrocytes transfected with a Sirt1 pro-Luc plasmid, but failed to activate the pro-mutant Luc reporter (Figure 5E & F). Furthermore, we found that 1,25(OH)2D3 inhibited ROS accumulation, as determined by dihydroethidium (DHE) staining (Figure 5G & H), increased chondrocyte proliferation (Figure 5I & J), decreased chondrocyte senescence (Figure 5K & L) and up-regulated p16 protein expression (Figure 5M & N) in human articular chondrocyte treated with IL-1β, and these parameters were largely blocked by treatment with the Sirt1 inhibitor, Ex527 (Figure 5G-N). These results demonstrated that 1,25(OH)2D3 suppresses IL-1β-induced chondrocyte senescence by VDR-mediated transcriptional up-regulation of Sirt1.
VDR deficiency induces chondrocyte senescence and a knee OA phenotype. (A) Representative micrographs of articular cartilage sections of 6-month-old wild-type mice (n=5) immunostained for VDR using two sources of VDR antibody (Abcam 3508 and Santa Cluz D6). (B) Western blot detection of VDR in monolayer-cultured wild-type mouse chondrocytes. Proteins isolated from wild-type BM-BMSCs and VDR-/- chondrocytes were loaded as positive and negative control, respectively. (C) Marked increase of Cyp24a1 and VDR mRNA levels in 1,25(OH)2D3-treated wild-type mouse chondrocytes, but not in VDR-/- cells demonstrated by RT-PCR. GAPDH was used as a loading control. (D) Immunofluorescence staining for VDR in mouse chondrocytes after treatment with vehicle or 1,25(OH)2D3 for 24 h. 1,25(OH)2D3-treated cells showed an increased localization of VDR in the nucleus compared with vehicle-treated cells. (E-H) Representative images of sections from wild-type (n=3) and VDR-/- (n=3) chondrocyte micromass cultures stained with (E) safranin O and immunostained for (F) collagen II, (G) Mmp13 and (H) p16 (scale bars, 200 µm). (I) Relative mRNA levels of col2a1, aggrecan, Mmp13 and p16 in macromass cultures above. n=3 wells per group. Representative micrographs of sections from 6-month-old wild-type (n=4) and VDR-/- (n=4) mice on the rescue diet stained (J) with safranin O and immunostained for (K) Mmp13, (L) collagen X and (M) p16. (N) Statistical analysis of Mmp13+, collagen X+ and p16+ cells. (O) Relative mRNA levels of p16, IL-1β, IL-6, Mmp13, aggrecan and col2a1 in wild-type and VDR-/- cartilage. *: p< 0.05, **: p< 0.01, ***: p<0.001.
1,25(OH)2D3 suppresses IL-1β-induced chondrocyte senescence by VDR-mediated transcriptional up-regulation of Sirt1. (A) Western blot detection of Sirt1 in articular chondrocytes isolated from WT (n=3) and VDR-/- (n=3) mice. (B) Western blot detection of Sirt1 in human articular chondrocyte cultures treated with IL-1β alone or combined with 1,25(OH)2D3 for indicated times. (C) Relative Sirt1 mRNA levels in human articular chondrocytes cultured with 1,25(OH)2D3 for 3, 6 and 12 hours. (D) Upper: VDRE-like elements in human Sirt1 promoter region highlighted in red. Lower: Chromatin immunoprecipitation (ChIP) with IgG antibody, H3 antibody (positive control) or VDR antibody were performed in human articular chondrocytes and relative enrichment of human Sirt1 promoter was determined using RT-PCR assay. (E) Human Sirt1 promoter or Sirt1 promoter mutant Luc-plasmid were transfected into human articular chondrocytes followed by vehicle or 1,25(OH)2D3 treatment for 12 hours and (F) relative luciferase activity was analyzed after 48 hours. n=3 replicates per condition. Human articular chondrocytes were cultured with vehicle, 1,25(OH)2D3, or 1,25(OH)2D3 plus Ex527 (a Sirt1 inhibitor) in the presence of IL-1β and stained (G) with dihydroethidium (DHE) to determine relative ROS levels, (I) cytochemically for EdU and (K) cytochemically for SA-β-gal. n=3 replicates per condition. Statistical analysis of (H) DHE fluorescence intensity, (J) EdU-positive cells and (L) SA-β-gal-positive cells. (M) Protein expression levels of p16 in human articular chondrocytes treated with vehicle or the Sirt1 inhibitor (Ex527) in the presence or absence of IL-1β and 1,25(OH)2D3. n=3 wells per condition. (N) Quantification of (M). *: p< 0.05, ***: p<0.001.
Due to the complexity and heterogeneity of knee OA pathogenesis, the combination of the removal of prevalent senescent cells and the rejuvenation of MSCs has been established as a potential therapeutic option [21-23, 38]. To assess whether overexpression of Sirt1 in MSCs could prevent 1,25(OH)2D deficiency-induced development of knee OA, we crossed Prx1-Sirt1 transgenic mice with 1α(OH)ase-/- mice to generate 1α(OH)ase-/- mice with Sirt1 overexpression in MSCs [1α(OH)ase-/- Sirt1Tg], and compared their knee phenotypes with 1α(OH)ase-/- mice on the rescue diet and wild-type littermates. Serum calcium, phosphorus, and PTH levels were normalized, whereas serum 1,25(OH)2D3 was undetectable in both 1α(OH)ase-/- and 1α(OH)ase-/- Sirt1Tg mice (Figure 6A). In addition, the percentage of Sirt1-positive chondrocytes was significantly increased in 1α(OH)ase-/- Sirt1Tg mice as compared with WT and 1α(OH)ase-/- mice (Figure 6B & C). The overexpression of Sirt1 in MSCs rescued all knee OA phenotypes observed in 12-month-old 1α(OH)ase-/- mice, including the significantly increased OARSI scores (Figure 6D & E), the increased percentages of collagen Ⅹ+ (Figure 6F & G), Mmp13+ (Figure 6H & I), and p16+ (Figure 6J & K) chondrocytes, and the increased mRNA levels of IL-1α, IL-1β, IL-6, Mmp3, and Mmp13 (Figure 6L). These results demonstrated that overexpression of Sirt1 in MSCs could partially prevent the development of 1,25(OH)2D deficiency-induced knee OA by inhibiting chondrocyte senescence and SASP.
The cause of aging-related OA is likely complex and the direct roles of VDR and 1,25(OH)2D in OA are largely unknown or controversial. In this study we demonstrated, for the first time, a direct link between the 1,25(OH)2D-VDR signaling pathway and the onset of age-related OA at tissue-, cellular-, and molecular-levels using multiple genetically-modified mouse lines and comprehensive analytical techniques, including histopathology, 3D cultures of human and mouse chondrocytes, immunofluorescence and immunohistochemistry staining, chromatin immunoprecipitation and dual-luciferase assay. Our key findings are 1) 1,25(OH)2D deficiency accelerated the development of age-related spontaneous knee OA; 2) 1,25(OH)2D3 supplementation rescued all knee OA phenotypes of 1α(OH)ase-/- mice in vivo, and 1,25(OH)2D3 rescued IL-1β-induced chondrocyte OA phenotypes in vitro; 3) VDR knockout mice exhibited knee OA phenotypes; 4) 1,25(OH)2D3 up-regulated Sirt1 through VDR mediated transcription in human articular chondrocytes; 5) Overexpression of Sirt1 in MSCs rescued knee OA phenotypes in 1α(OH)ase-/- mice. Our results indicate that age-related knee OA induced by 1,25(OH)2D deficiency occurs in a calcium-, phosphorus- and 25(OH)D-independent manner. Our results may help to explain the inconsistent results surrounding the association between vitamin D levels, as measured by serum 25(OH)D, and OA. Serum 25(OH)D levels are not decreased in 1α(OH)ase-/- mice relative to wild-type mice [12]; however, our previous studies demonstrated that 1α(OH)ase-/- mice display female infertility, male infertility, hypertension, premature aging, a high incidence of multiple tumors and osteoporosis [31, 33, 39-44]. These results suggest that sufficient 1,25(OH)2D rather than 25(OH)D levels are necessary for the prevention of vitamin D deficiency-induced diseases including knee OA.
Vitamin D supplementation may improve pain and function in patients with knee OA; however, the evidence from observational studies and randomized controlled trials are controversial [45, 46]. To determine whether supplementation of exogenous 1,25(OH)2D3 can prevent knee OA caused by 1,25(OH)2D deficiency, we performed in vivo supplementation of exogenous 1,25(OH)2D3 in 1α(OH)ase-/- mice and in vitro treatment of 1,25(OH)2D3 in 2D or 3D cultures of human articular chondrocytes with IL-1β induction. Our in vivo experiments revealed that all knee OA phenotypes observed in 1α(OH)ase-/- mice were rescued by supplementation with exogenous 1,25(OH)2D3. In vitro 2D models have used cytokines such as IL-1β to induce an OA phenotype which has allowed for the screening of chondroprotective compounds to attenuate the catabolic factors involved in articular chondrocyte degradation [47]. Both healthy and OA articular chondrocytes express the IL-1 receptor type I. The effect of IL-1β on chondrocytes can be generally described as catabolic and involves the upregulation of aggrecanases and MMPs, the induction of additional inflammatory mediators, as well as the downregulation of chondrogenic extracellular matrix synthesis [48]. A recent study implies that omentin-1 could have clinical relevance in the treatment of osteoarthritis by inhibiting IL-1β-induced chondrocyte senescence [49]. In our in vitro experiments, we found that IL-1β inhibited the proliferation of human articular chondrocytes, and stimulated oxidative stress and articular chondrocyte senescence. Treatment with 1,25(OH)2D3, but not with vehicle, rescued IL-1β-induced chondrocyte phenotypes. Consequently, supplementation with exogenous 1,25(OH)2D3 prevented the onset and progression of knee OA by stimulating chondrocyte proliferation and extracellular matrix synthesis, and by inhibiting chondrocyte differentiation and degradation of extracellular matrix proteins and chondrocyte senescence in articular cartilage (Figure 7). In a recent study, we reported that progressive down-regulation of the protein expression levels of 1α(OH)ase in kidney, intestine and bone of 3-, 9- and 18-month-old wild-type mice were associated with progressive bone loss [30]. In the current study, we found that the protein expression level of 1α(OH)ase in articular cartilage was also reduced in 12-month-old wild-type mice relative to 6-month-old wild-type mice. Consequently, whether vitamin D supplementation can prevent the onset and progression of knee OA may depend on the local 1α(OH)ase expression levels and activity that is required to convert supplementary vitamin D into active vitamin D.
Overexpression of Sirt1 in MSCs prevents 1,25(OH)2D deficiency-induced development of knee OA. (A) Serum calcium, phosphorus, intact PTH and 1,25(OH)2D3 levels in 12-month-old wild-type (n=5), 1ɑ(OH)ase-/- (n=5) and 1ɑ(OH)ase-/-Sirt1Tg mice (n=5) on the rescue diet (RD). (B) Representative micrographs of articular cartilage sections immunostained for Sirt1 in indicated groups of mice and (C) a quantitative analysis of the percentage of Sirt1-positive chondrocytes. Representative micrographs of knee sections stained with (D) safranin O, and immunostained for (F) collagen X, (H) Mmp13 and (J) p16. Evaluation for (E) OARSI score and quantification of (G) collagen X+, (I) Mmp13+ and (K) p16+ cells. (L) Relative mRNA levels of SASP factors including IL-1ɑ, IL-1β, IL-6, Mmp3, Mmp13 were analyzed using qPCR. n=5 mice per group. *: p< 0.05, **: p< 0.01, ***: p<0.001.
Model of mechanisms leading from 1,25(OH)2D3 to prevent age-related spontaneous knee osteoarthritis via VDR-mediated upregulation of Sirt1.
It has been reported that VDR is expressed in the articular cartilage of OA patients and human articular chondrocytes [50]. We therefore asked whether the development of knee OA induced by 1,25(OH)2D deficiency was mediated though VDR. To address this issue, we first confirmed that VDR is localized in mouse articular chondrocytes as demonstrated by both immunohistochemistry and Western blots. We also found that treatment with 1,25(OH)2D3 could up-regulate VDR expression and increase its nuclear translocation in human articular chondrocytes. Then we demonstrated that VDR deficiency resulted in the reduction of extracellular matrix protein synthesis and the increase of Mmp13 and p16 expression in 3D cultured mouse articular chondrocytes. VDR deficiency also induced knee OA phenotypes in vivo by enhancing extracellular matrix degradation and articular chondrocyte hypertrophic differentiation and stimulating chondrocyte senescence and SASP. Consequently, VDR deficiency could induce the development of knee OA, supporting the role of VDR mediated 1,25(OH)2D deficiency in the development of knee OA.
The expression of Sirt1 in articular cartilage is negatively associated with severity of knee OA [51], and cartilage-specific Sirt1-deficient mice develop accelerated OA progression under mechanical stress and aging [20]. Our current study demonstrated that the expression level of Sirt1 was dramatically reduced in articular chondrocytes from 1,25(OH)2D deficient mice, and this was rescued by supplementation with exogenous 1,25(OH)2D3. Treatment with 1,25(OH)2D3 significantly upregulated the expression of Sirt1 at both protein and mRNA levels in human articular chondrocytes in the presence of IL-1β, while the inhibitory effect of 1,25(OH)2D3 on IL-1β-induced chondrocyte senescence and p16 up-regulation were largely blocked by treatment with the Sirt1 inhibitor, Ex527. Furthermore, we demonstrated that 1,25(OH)2D3 up-regulated Sirt1 expression via VDR-mediated transcription in human articular chondrocytes. Consequently, Sirt1 appears to be a direct downstream target of 1,25(OH)2D.
In view of the fact that previous studies suggested that Sirt1 plays a preventive role against the development of OA [18, 20], we investigated whether overexpressing Sirt1 in MSCs could prevent 1,25(OH)2D deficiency-induced development of knee OA. We found that similar to supplementation with 1,25(OH)2D3, overexpressing Sirt1 in MSCs rescued 1,25(OH)2D-deficiency-induced knee OA by up-regulating Sirt1 expression in articular chondrocytes, inhibiting the differentiation of chondrocytes, degradation of extracellular matrix proteins, chondrocyte senescence and SASP in articular cartilage (Figure 7). Therefore, our results indicate that 1,25(OH)2D can prevent the development of knee OA by inhibiting extracellular matrix degradation, cellular senescence and SASP via Sirt1 mediation (Figure 7).
We recently showed that 1,25(OH)2D deficiency accelerates age-related bone loss [31], raising the question of whether 1,25(OH)2D deficiency-induced OA is actually driven by the loss of subchondral bone. In this study, we found that the subchondral bone volume of 1,25(OH)2D deficient mice was not significantly altered compared with WT controls. In contrast, VDR-knockout chondrocytes displayed cellular senescence and an OA phenotype in vitro, and 1,25(OH)2D also protected chondrocytes against an IL-1β-induced OA phenotype in vitro, including increased ROS levels and cellular senescence. These findings demonstrate that 1,25(OH)2D can prevent an OA phenotype via VDR by directly targeting chondrocytes. However, the clarification of the role of subchondral bone in contributing to age-related OA is a potential limitation of our study, and aged mice with specific deletion of VDR in subchondral bone and in chondrocytes should be constructed to further examine the direct contribution of vitamin D deficiency in subchondral bone and in chondrocytes to an OA phenotype with aging.
Overall, the results of this study provide a model that suggests that 1,25(OH)2D via the VDR prevents aging-related knee OA by up-regulating Sirt1 through enhanced VDR mediated gene transcription, activating p16/p21 signaling and decreasing oxidative stress, promoting articular chondrocyte proliferation and extracellular matrix protein synthesis, and reducing articular chondrocyte senescence and SASP (Figure 7). This study therefore not only identifies novel mechanisms of 1,25(OH)2D deficiency in accelerating age-related knee OA development, but may also provide experimental and theoretical evidence for potential utilization of 1,25(OH)2D3 or downstream targets, including Sirt1 activators to prevent age-related knee OA.
Supplementary figures and table.
This work was supported by grants from the National Key R&D Program of China (2018YFA0800800 to DM), the National Natural Science Foundation of China (81730066 and 81230009 to DM) and the Canadian Institutes of Health Research (PJT-152963 to DG). We thank GENEWIZ for their assistance in RNA-sequencing.
D.M., R.Y. and D.G. conceived the project. R.Y., J.C., J.Z., J.L. and R.Q. performed most of the experiments, analyzed, and compiled the data. N.L. helped with experiments. R.Y., D.M., and D.G. participated in writing or editing the paper.
The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.
The authors have declared that no competing interest exists.
1. Kraus VB, Blanco FJ, Englund M, Karsdal MA, Lohmander LS. Call for standardized definitions of osteoarthritis and risk stratification for clinical trials and clinical use. Osteoarthritis and cartilage. 2015;23:1233-41
2. Holick MF, Chen TC. Vitamin D deficiency: a worldwide problem with health consequences. Am J Clin Nutr. 2008;87:1080S-6S
3. Wang H, Chen W, Li D, Yin X, Zhang X, Olsen N. et al. Vitamin D and Chronic Diseases. Aging Dis. 2017;8:346-53
4. Cao Y, Winzenberg T, Nguo K, Lin J, Jones G, Ding C. Association between serum levels of 25-hydroxyvitamin D and osteoarthritis: a systematic review. Rheumatology (Oxford, England). 2013;52:1323-34
5. Tripathy SK, Gantaguru A, Nanda SN, Velagada S, Srinivasan A, Mangaraj M. Association of vitamin D and knee osteoarthritis in younger individuals. World J Orthop. 2020;11:418-25
6. Vaishya R, Vijay V, Lama P, Agarwal A. Does vitamin D deficiency influence the incidence and progression of knee osteoarthritis? - A literature review. J Clin Orthop Trauma. 2019;10:9-15
7. Bergink AP, Trajanoska K, Uitterlinden AG, van Meurs JBJ. Mendelian randomization study on vitamin D levels and osteoarthritis risk: a concise report. Rheumatology (Oxford, England). 2021;60:3409-12
8. Sanghi D, Mishra A, Sharma AC, Singh A, Natu SM, Agarwal S. et al. Does vitamin D improve osteoarthritis of the knee: a randomized controlled pilot trial. Clin Orthop Relat Res. 2013;471:3556-62
9. McAlindon T, LaValley M, Schneider E, Nuite M, Lee JY, Price LL. et al. Effect of vitamin D supplementation on progression of knee pain and cartilage volume loss in patients with symptomatic osteoarthritis: a randomized controlled trial. JAMA. 2013;309:155-62
10. Jin X, Jones G, Cicuttini F, Wluka A, Zhu Z, Han W. et al. Effect of Vitamin D Supplementation on Tibial Cartilage Volume and Knee Pain Among Patients with Symptomatic Knee Osteoarthritis: A Randomized Clinical Trial. JAMA. 2016;315:1005-13
11. Arden NK, Cro S, Sheard S, Dore CJ, Bara A, Tebbs SA. et al. The effect of vitamin D supplementation on knee osteoarthritis, the VIDEO study: a randomised controlled trial. Osteoarthritis and cartilage. 2016;24:1858-66
12. Panda DK, Miao D, Tremblay ML, Sirois J, Farookhi R, Hendy GN. et al. Targeted ablation of the 25-hydroxyvitamin D 1alpha -hydroxylase enzyme: evidence for skeletal, reproductive, and immune dysfunction. Proc Natl Acad Sci U S A. 2001;98:7498-503
13. Li YC, Pirro AE, Amling M, Delling G, Baron R, Bronson R. et al. Targeted ablation of the vitamin D receptor: an animal model of vitamin D-dependent rickets type II with alopecia. Proc Natl Acad Sci U S A. 1997;94:9831-5
14. Goltzman D. Vitamin D action: Lessons learned from genetic mouse models. Ann N Y Acad Sci. 2010;1192:145-52
15. Miao D, Goltzman D. Probing the Scope and Mechanisms of Calcitriol Actions Using Genetically Modified Mouse Models. JBMR Plus. 2021;5:e10434
16. Shen M, Luo Y, Niu Y, Chen L, Yuan X, Goltzman D. et al. 1,25(OH)2D deficiency induces temporomandibular joint osteoarthritis via secretion of senescence-associated inflammatory cytokines. Bone. 2013;55:400-9
17. Almeida M, Porter RM. Sirtuins and FoxOs in osteoporosis and osteoarthritis. Bone. 2019;121:284-92
18. Gabay O, Oppenhiemer H, Meir H, Zaal K, Sanchez C, Dvir-Ginzberg M. Increased apoptotic chondrocytes in articular cartilage from adult heterozygous SirT1 mice. Annals of the rheumatic diseases. 2012;71:613-6
19. Gabay O, Sanchez C, Dvir-Ginzberg M, Gagarina V, Zaal KJ, Song Y. et al. Sirtuin 1 enzymatic activity is required for cartilage homeostasis in vivo in a mouse model. Arthritis and rheumatism. 2013;65:159-66
20. Matsuzaki T, Matsushita T, Takayama K, Matsumoto T, Nishida K, Kuroda R. et al. Disruption of Sirt1 in chondrocytes causes accelerated progression of osteoarthritis under mechanical stress and during ageing in mice. Annals of the rheumatic diseases. 2014;73:1397-404
21. Deng L, Ren R, Liu Z, Song M, Li J, Wu Z. et al. Stabilizing heterochromatin by DGCR8 alleviates senescence and osteoarthritis. Nature communications. 2019;10:3329
22. Fu L, Hu Y, Song M, Liu Z, Zhang W, Yu FX. et al. Up-regulation of FOXD1 by YAP alleviates senescence and osteoarthritis. PLoS biology. 2019;17:e3000201
23. Ren X, Hu B, Song M, Ding Z, Dang Y, Liu Z. et al. Maintenance of Nucleolar Homeostasis by CBX4 Alleviates Senescence and Osteoarthritis. Cell reports. 2019;26:3643-56 e7
24. Kawanami A, Matsushita T, Chan YY, Murakami S. Mice expressing GFP and CreER in osteochondro progenitor cells in the periosteum. Biochem Biophys Res Commun. 2009;386:477-82
25. Wang H, Hu Z, Wu J, Mei Y, Zhang Q, Zhang H. et al. Sirt1 Promotes Osteogenic Differentiation and Increases Alveolar Bone Mass via Bmi1 Activation in Mice. Journal of bone and mineral research: the official journal of the American Society for Bone and Mineral Research. 2019;34:1169-81
26. Sun W, Qiao W, Zhou B, Hu Z, Yan Q, Wu J. et al. Overexpression of Sirt1 in mesenchymal stem cells protects against bone loss in mice by FOXO3a deacetylation and oxidative stress inhibition. Metabolism. 2018;88:61-71
27. Chen H, Hu X, Yang R, Wu G, Tan Q, Goltzman D. et al. SIRT1/FOXO3a axis plays an important role in the prevention of mandibular bone loss induced by 1,25(OH)2D deficiency. Int J Biol Sci. 2020;16:2712-26
28. Panda DK, Miao D, Bolivar I, Li J, Huo R, Hendy GN. et al. Inactivation of the 25-hydroxyvitamin D 1alpha-hydroxylase and vitamin D receptor demonstrates independent and interdependent effects of calcium and vitamin D on skeletal and mineral homeostasis. J Biol Chem. 2004;279:16754-66
29. Wang H, Hu Z, Wu J, Mei Y, Zhang Q, Zhang H, et al. Sirt1 Promotes Osteogenic Differentiation and Increases Alveolar Bone Mass via Bmi1 Activation in Mice
30. Qiao W, Yu S, Sun H, Chen L, Wang R, Wu X. et al. 1,25-Dihydroxyvitamin D insufficiency accelerates age-related bone loss by increasing oxidative stress and cell senescence. Am J Transl Res. 2020;12:507-18
31. Yang R, Chen J, Zhang J, Qin R, Wang R, Qiu Y. et al. 1,25-Dihydroxyvitamin D protects against age-related osteoporosis by a novel VDR-Ezh2-p16 signal axis. Aging cell. 2020;19:e13095
32. Yang R, Zhang J, Li J, Qin R, Chen J, Wang R. et al. Inhibition of Nrf2 degradation alleviates age-related osteoporosis induced by 1,25-Dihydroxyvitamin D deficiency. Free Radical Biology and Medicine. 2022;178:246-61
33. Chen L, Yang R, Qiao W, Zhang W, Chen J, Mao L. et al. 1,25-Dihydroxyvitamin D exerts an antiaging role by activation of Nrf2-antioxidant signaling and inactivation of p16/p53-senescence signaling. Aging cell. 2019;18:e12951
34. Miao D, Su H, He B, Gao J, Xia Q, Zhu M. et al. Severe growth retardation and early lethality in mice lacking the nuclear localization sequence and C-terminus of PTH-related protein. Proc Natl Acad Sci U S A. 2008;105:20309-14
35. Zhu X, Chen F, Lu K, Wei A, Jiang Q, Cao W. PPARgamma preservation via promoter demethylation alleviates osteoarthritis in mice. Annals of the rheumatic diseases. 2019;78:1420-9
36. Matsushita T, Sasaki H, Takayama K, Ishida K, Matsumoto T, Kubo S. et al. The overexpression of SIRT1 inhibited osteoarthritic gene expression changes induced by interleukin-1β in human chondrocytes. Journal of orthopaedic research: official publication of the Orthopaedic Research Society. 2013;31:531-7
37. Tong W, Geng Y, Huang Y, Shi Y, Xiang S, Zhang N. et al. In vivo Identification and Induction of Articular Cartilage Stem Cells by Inhibiting NF-κB Signaling in Osteoarthritis. Stem Cells. 2015;33:3125-37
38. Jeon OH, Kim C, Laberge RM, Demaria M, Rathod S, Vasserot AP. et al. Local clearance of senescent cells attenuates the development of post-traumatic osteoarthritis and creates a pro-regenerative environment. Nature medicine. 2017;23:775-81
39. Chen L, Yang R, Qiao W, Yuan X, Wang S, Goltzman D. et al. 1,25-Dihydroxy vitamin D prevents tumorigenesis by inhibiting oxidative stress and inducing tumor cellular senescence in mice. Int J Cancer. 2018;143:368-82
40. Sun H, Qiao W, Cui M, Yang C, Wang R, Goltzman D. et al. The Polycomb Protein Bmi1 Plays a Crucial Role in the Prevention of 1,25(OH)2 D Deficiency-Induced Bone Loss. Journal of bone and mineral research: the official journal of the American Society for Bone and Mineral Research. 2020;35:583-95
41. Sun W, Chen L, Zhang W, Wang R, Goltzman D, Miao D. Active vitamin D deficiency mediated by extracellular calcium and phosphorus results in male infertility in young mice. American journal of physiology Endocrinology and metabolism. 2015;308:E51-62
42. Sun W, Xie H, Ji J, Zhou X, Goltzman D, Miao D. Defective female reproductive function in 1,25(OH)2D-deficient mice results from indirect effect mediated by extracellular calcium and/or phosphorus. American journal of physiology Endocrinology and metabolism. 2010;299:E928-35
43. Xue Y, Karaplis AC, Hendy GN, Goltzman D, Miao D. Exogenous 1,25-dihydroxyvitamin D3 exerts a skeletal anabolic effect and improves mineral ion homeostasis in mice that are homozygous for both the 1alpha-hydroxylase and parathyroid hormone null alleles. Endocrinology. 2006;147:4801-10
44. Zhou C, Lu F, Cao K, Xu D, Goltzman D, Miao D. Calcium-independent and 1,25(OH)2D3-dependent regulation of the renin-angiotensin system in 1alpha-hydroxylase knockout mice. Kidney Int. 2008;74:170-9
45. Gao XR, Chen YS, Deng W. The effect of vitamin D supplementation on knee osteoarthritis: A meta-analysis of randomized controlled trials. Int J Surg. 2017;46:14-20
46. Hussain S, Singh A, Akhtar M, Najmi AK. Vitamin D supplementation for the management of knee osteoarthritis: a systematic review of randomized controlled trials. Rheumatol Int. 2017;37:1489-98
47. Samvelyan HJ, Hughes D, Stevens C, Staines KA. Models of Osteoarthritis: Relevance and New Insights. Calcif Tissue Int. 2021;109:243-56
48. Jenei-Lanzl Z, Meurer A, Zaucke F. Interleukin-1beta signaling in osteoarthritis - chondrocytes in focus. Cell Signal. 2019;53:212-23
49. Chai B, Zheng ZH, Liao X, Li KY, Liang JS, Huang YX. et al. The protective role of omentin-1 in IL-1beta-induced chondrocyte senescence. Artif Cells Nanomed Biotechnol. 2020;48:8-14
50. Tetlow LC, Woolley DE. Expression of vitamin D receptors and matrix metalloproteinases in osteoarthritic cartilage and human articular chondrocytes in vitro. Osteoarthritis and cartilage. 2001;9:423-31
51. Li Y, Xiao W, Wu P, Deng Z, Zeng C, Li H. et al. The expression of SIRT1 in articular cartilage of patients with knee osteoarthritis and its correlation with disease severity. J Orthop Surg Res. 2016;11:144
Corresponding authors: Renlei Yang, Ph.D., Department of Plastic Surgery, Affiliated Friendship Plastic Surgery Hospital of Nanjing Medical University, Nanjing, Jiangsu, 210029, The People's Republic of China. Tel & FAX: 011-86-25-8686-9377; E-mail: yang_renleicom; Dengshun Miao, M.D., Ph.D., Department of Plastic Surgery, Affiliated Friendship Plastic Surgery Hospital of Nanjing Medical University, Nanjing, Jiangsu, 210029, The People's Republic of China. Tel & FAX: 011-86-25-8686-9377; E-mail: dsmiaoedu.cn.
Received 2022-9-9
Accepted 2022-12-9
Published 2023-1-1