Osteoarthritis-Like Changes in Bardet–Biedl Syndrome Mutant Ciliopathy Mice (Bbs1M390R/M390R): Evidence for a Role of Primary Cilia in Cartilage Homeostasis and Regulation of Inflammation

Osteoarthritis (OA) is a debilitating inflammation related disease characterized by joint pain and effusion, loss of mobility, and deformity that may result in functional joint failure and significant impact on quality of life. Once thought of as a simple “wear and tear” disease, it is now widely recognized that OA has a considerable metabolic component and is related to chronic inflammation. Defects associated with primary cilia have been shown to be cause OA-like changes in Bardet–Biedl mice. We examined the role of dysfunctional primary cilia in OA in mice through the regulation of the previously identified degradative and pro-inflammatory molecular pathways common to OA. We observed an increase in the presence of pro-inflammatory markers TGFβ-1 and HTRA1 as well as cartilage destructive protease MMP-13 but a decrease in DDR-2. We observed a morphological difference in cartilage thickness in Bbs1^M390R/M390R mice compared to wild type (WT). We did not observe any difference in OARSI or Mankin scores between WT and Bbs1^M390R/M390R mice. Primary cilia appear to be involved in the upregulation of biomarkers, including pro-inflammatory markers common to OA.

Introduction

Osteoarthritis (OA) is a debilitating disease characterized by joint pain and effusion, loss of mobility, and deformity that may result in functional joint failure and significant impact on quality of life (1–3). It is also one of the most common chronic diseases in the United States, with recent estimates suggesting that more than 27 million adults suffer from clinical OA in this country. Vulnerable populations include the elderly, where 80% of people above 65 years of age are symptomatic for OA (Felson et al., 1988; Guccione, 1994; Lawrence et al., 2008), and the obese. As upward obesity trends continue it is probable that an increase in clinical OA will also be observed (1–3). There is currently no cure, and treatment options are minimal, with pain management and surgical joint replacement procedures being the only reprieve.

Once thought of as a passive “wear and tear” disease, it is now widely recognized that OA has a considerable metabolic component. It has been shown that whether the onset of the disease is due to aging, obesity, joint misalignment, acute injury, or genetic predisposition, all result in the activation of a common molecular pathway (Murwantoko et al., 2004; Lawrence et al., 2008; Polur et al., 2010; Holt et al., 2012; Larkin et al., 2013). Previous work has successfully identified key mediators in this pathway, including transforming growth factor beta 1 (TGF-β1), high temperature receptor A-1 (HTRA1), discoidin domain-containing receptor 2 (DDR-2), matrix metalloprotease 13 (MMP-13), and other matrix-degrading enzymes that are secreted by chondrocytes and contribute to the degradation of articular cartilage. Stress to the chondrocyte leads to an increase in TGF-β1-mediated chondroplasia but also has the potential to stimulate osteogenesis. Perhaps to offset this phenomena, HTRA1 is secreted to cleave TGF-β1. Being a non-specific serine protease, HTRA-1 also cleaves proteins of the pericellular matrix and brings the collagenous matrix in contact with the surface of chondrocytes, binding to DDR2 and subsequently increasing production of MMP-13, which specifically cleaves the type II collagen, among other proteins, making up the majority of articular cartilage.

Activation of MMP-13 gene expression appears to occur through a number of molecular pathways that work through either inflammation or primary cilia. That is not to say there are not some common themes. Stress-inducible nuclear protein 1 (NUPR1) has been shown to regulate MMP-13 expression in vitro (Yammani and Loeser, 2014). Yammani and Loeser (2014) showed that NUPR1, expressed in cartilage, is required for expression of MMP-13 via IL-1β. This might be a pathway for the catabolic effects of OA to be mediated through inflammation. This is especially interesting in light of the study done by Xu et al. (2015) in which they analyzed differential expression of genes in cartilage involved in OA and rheumatoid arthritis (RA). While these researchers identified multiple genes associated with the regulation of MMPs, the predominant ones were associated with inflammation. This might give greater credence for the role of early inflammatory signals (i.e., AGEs, IL-1) in the initiation and progression of OA. While more obvious, a similar role for inflammation appears to be present in RA. Araki et al. (2016) reported that histone methylation and the binding of signal transducer activator of transcription 3 (STAT3) was associated with RA and OA. They report that histone H3 methylation is associated with elevated expression of MMP-1, -3, -9, and -13. However, STAT3 was shown to increase expression, either spontaneous or IL-6 activated, of MMPs 1, 3, and 13 but not 9. As previously indicated, primary cilia appear to also be involved in OA. Sugita et al. (2015) reported that transcription factor hairy and enhancer of split-1 (Hes1) is involved in the upregulation of expression of MMP-13. Normally Hes1 acts as a transcriptional repressor but under the influence of calcium/calmodulin-dependent protein kinase 2 (CaMK2) it becomes a transcriptional activator, thus upregulating MMP-13 expression (Ju et al., 2004). Thus Hes1 acts to increase expression of MMP-13. It is of particular interest to note that HES1 acts through Notch signaling pathway (Kageyama et al., 2007). Notch has previously been shown to modulate sonic hedgehog signaling and work through primary cilia (Ezratty et al., 2011; Kong et al., 2015). In an apparent unrelated mechanism, Niebler et al. (2015) showed that the transcription factor AP-2e is intimately involved in the upregulation of MMP-13 as OA progresses.

Bardet–Biedl syndrome (BBS, MIM 209900) is a pleiotropic genetically heterogeneous disorder characterized by obesity, retinopathy, polydactyly, renal anomalies including polycystic kidney disease, intellectual disabilities and hypo-genitalism (Chiang et al., 2004; Kaushik et al., 2009; Marion et al., 2011). To date, 21 BBS genes have been identified (Katsanis et al., 2000; Slavotinek and Biesecker, 2000; Mykytyn et al., 2001; Nishimura et al., 2001, 2005; Ansley et al., 2003; Badano et al., 2003; Chiang et al., 2004, 2006; Li et al., 2004; Stoetzel et al., 2006a,b, 2007). Although the cellular functions of BBS proteins are not yet fully understood, evidence from a variety of organisms, including BBS mouse models, demonstrate that BBS proteins play a role in cilia assembly and/or function, as well as intracellular vesicle trafficking. BBS mouse models lack spermatozoa flagella (Fan et al., 2004; Badano et al., 2005). Knock down of multiple BBS genes in zebrafish have been shown to interfere with function of nodal cilia, as well as to result in delay of intracellular vesicle transport. C. elegans homologs of Bbs1, Bbs2, Bbs7 and Bbs8 have been shown to be expressed in ciliated cells and Bbs8 was found in the ciliary basal body but not in the microtubule-based ciliary axoneme (Ansley et al., 2003). BBS1, BBS4, BBS5, BBS7, and BBS8 have also been reported to localize to the centrosome and/or basal body of widely used mammalian cell lines (Badano et al., 2005). The depletion of BBS4 through RNA interference has been shown to lead to the disruption of cytoplasmic microtubule arrays, the mis-localization of some pericentriolar proteins including pericentriolar material 1 protein (PCM1), and cell division arrest and apoptosis in cultured cells (Kim et al., 2004). Recent studies indicate that the seven most conserved BBS proteins (BBS1, BBS2, BBS4, BBS5, BBS7, BBS8, and BBS9), along with BBS18 form a stable complex known as the BBSome. The BBSome associates with centriolar satellites in the cytoplasm and transitions to a membrane-associated form at the base of the cilium. Depletion of some components of the BBSome affects ciliogenesis by interfering with membrane trafficking to the primary cilium (Nachury et al., 2007). Studies involving renal cells have shown a predisposition of primary cilia to act as surface receptors (Marion et al., 2011) and function as both mechano- and chemoreceptors for specific environments (Poole et al., 1985, 1997, 2001; Kouri et al., 1996, 1998; Jensen et al., 2004; McGlashan et al., 2006; Kaushik et al., 2009). For a recent reviews of BBS genes and associated diseases see the following papers (Heon et al., 2016; Kaur et al., 2017; Weihbrecht et al., 2017). We chose to work with mice with the Bbs1^M390R/M390R since it accounts for 25% of all human BBS cases and is directly involved with the BBSome.

Previous work has visualized the presence of primary cilia on chondrocytes, with each chondrocyte containing a single primary cilium (Scherft and Daems, 1967; Hart, 1968; Wilsman, 1978; Meier-Vismara et al., 1979; Wilsman et al., 1980; Kaushik et al., 2009). One study demonstrated cartilage abnormalities associated with defective primary cilia (Kaushik et al., 2009). Studies employing Oak Ridge polycystic disease mouse chondrocytes show that the primary cilia in chondrocytes are not involved in the initial mechanosensation, but are associated with signal processing in response to increased loads (Mykytyn and Sheffield, 2004). We hypothesized that dysfunctional primary cilia result in OA in Bbs1^M390R/M390R mice as they age through intra-cellular pathways resulting in up-regulation of the previously identified HTRA1-DDR2-MMP13 degradative pathway common to OA.

Materials and Methods

Tissue Samples

Knees from wild type (WT) n = 44 (n = 8 10–14 weeks old; n = 24 15–18 weeks old; n = 12 21+ weeks old) and Bbs1^M390R/M390R n = 61 (n = 18 10–14 weeks old; n = 39 15–18 weeks old; n = 4 21+ weeks old) mutant mice were obtained from the laboratory of Val Sheffield at the University of Iowa. Knees were harvested and fixed in 4% paraformaldehyde prior to decalcification and embedding according to a previously reported procedure (Spencer et al., 2015). Knees were serially sectioned and sections analyzed histologically after staining with Safranin O and fast green as previously reported (Spencer et al., 2015). Using a light microscope equipped with a digital camera, photographs of each knee joint were taken at 100× and 200× magnifications. The articular cartilage in two representative sections from each stained slide were analyzed using the OARSI scoring system (Glasson et al., 2010) to quantify the pathological state of each joint, with a score of zero representing unaltered articular cartilage and six representing severe OA.

Immunohistochemical Analysis

Immunohistochemistry (IHC) was performed on slides representing serial sections of mouse knee joints from all animals. Separate slides were stained with antibodies against HTRA1, DDR-2, MMP-13, and TGF-β1. Each slide was deparaffinized and then blocked with 5% bovine serum albumin for 1 h. Primary antibodies against HTRA-1 (ab38611) (Abcam, Cambridge, MA, United States), DDR-2 (SC-8989) (Santa Cruz Biotechnology, Santa Cruz, CA, United States), MMP-13 (AB8120) (Chemicon, Temecula, CA, United States) and TGF-β1 (Abcam- ab92486, Rabbit Polyclonal IgG) were used. All antibodies were diluted 1:200, applied to specimens, and incubated overnight at 4°C. On the second day, samples were rinsed with PBS and then incubated with an avidin/biotin ABC mix (Vectastain elite ABC Kit). Slides were rinsed again with PBS and incubated with a species appropriate biotinylated secondary antibody. After a third rinse, a color reaction was initiated using a peroxidase substrate (Vector Labs, NovaRED). Negative controls were prepared by staining without the addition of primary antibody. Differences in staining intensity were compared qualitatively with WT controls. Blind counting of stained cells was performed using ImageJ (NIH, Bethesda, MD, United States). IHC and histological analysis focused on different parts of the joint to demonstrate the universal nature of the OA in the Bbs1^M390R/M390R model. Animals were grouped into three general age categories (10–14 weeks, 15–18 weeks, and 21+ weeks) for comparison purposes. The age categories selected indicate young, middle age, and old mice. For the purposes of IHC and histological analysis, i.e., cell counting (see section “ImageJ Analysis” below), the number of mice used in each category was 10–14 weeks: WT n = 8, Bbs1^M390R/M390R mutant n = 18, 15–18 weeks: WT n = 24, Bbs1^M390R/M390R mutant n = 39, 21+ weeks: WT n = 12, Bbs1^M390R/M390R mutant n = 4.

Zebrafish Immunohistochemistry

Zebrafish IHC techniques were employed as previously described (Suli et al., 2014). Antibodies used for visualization: Primary cilia: anti-mouse acetylated tubulin (1:500, Sigma, T7451) with anti-mouse IgG Alexa488 (1:400, Invitrogen, A-28175); Primary DDR2: anti-mouse DDR2 (1:200, Santa Cruz SC-8989) with anti-rabbit IgG Alexa405 (1:400, Invitrogen, A31556).

Larvae were fixed with 4% paraformaldehyde either for 2 h at room temperature, washed three times, 20 min each, with PBST (0.1% Tween in PBS) and incubated 1 h in distilled water. They were placed in block solution [1% bovine serum albumin, 1% dimethyl sulfoxide (DMSO), and 0.02% sodium azide in PBST, 10% normal goat serum] for 1 h and then incubated with primary antibody overnight at 4°C. After four 20 min washes with PBST, they were incubated with secondary antibody for 3 h at room temperature, washed four times, 10 min each, in PBST and cleared in 50% glycerol/PBS. Embryos were imaged using an Olympus FV1000 confocal microscope.

Osteoarthritis Scoring

The articular cartilage in two representative sections from each stained slide was analyzed using the OARSI scoring system to quantify the pathological state of each joint, with a score of zero representing unaltered articular cartilage and six representing severe OA. Overall OARSI scoring was based on an osteoarthritic damage 0–6 subjective scoring system applied to all four quadrants of the knee (Glasson et al., 2010). The articular cartilage in two representative sections from each stained slide was also analyzed using a modified Mankin score to quantify the pathological state of each joint, with a score of zero representing unaltered articular cartilage and 14 representing severe OA. Overall Mankin scoring was based on a subset of scores including cartilage erosion score (0–6), chondrocyte periphery staining (0–2), spatial arrangement of chondrocytes (0–3), and background staining intensity (0–3) (Mankin et al., 1971; Xu et al., 2003; Bomsta et al., 2006). Statistical significance of the combined Mankin scores for the 28-day Bbs1^M390R/M390R mutant and wild-type surgery groups using a two-way ANOVA test conducted by the Department of Statistics at Brigham Young University. As in the case of IHC, animals were grouped into three general age categories (10–14 weeks, 15–18 weeks, and 21+ weeks) for comparison purposes. The age categories selected indicate young, middle age, and old mice. For the purposes of OARSI scoring analysis the number of mice used in each category was 10–14 weeks: WT n = 8, Bbs1^M390R/M390R mutant n = 18, 15–18 weeks: WT n = 24, Bbs1^M390R/M390R mutant n = 39, 21+ weeks: WT n = 12, Bbs1^M390R/M390R mutant n = 4.

ImageJ Analysis

The expression levels of HTRA1, DDR-2, MMP-13, and TGF-β1 were analyzed quantitatively by calculating the percentage of cells staining positive for the respective biomarkers and the total number of chondrocytes in a defined 200×900 pixel area of articular cartilage immediately distal to the tibial plateau (Larkin et al., 2013). All quantitative analysis was performed using ImageJ (National Institutes of Health, Bethesda, MD, United States). Cell counting was done by two independent investigators who were blinded to the strain of mouse (Bbs1^M390R/M390R mutant vs. WT). The quantitative results were subsequently analyzed statistically using an ANOVA test to detect differences in the mean percentages of positive staining for key OA biomarkers and mean chondrocyte counts between the Bbs1^M390R/M390R mutant and WT samples. The numbers of animals analyzed for each age group were previously indicated in the IHC methods section.

Adobe Photoshop Analysis

To compare the thickness of cartilage between Bbs1^M390R/M390R mutant vs. WT, cartilage measurements above the subchondral bone for each mouse specimen was captured using the Lasso Tool in Adobe Photoshop. This was done to measure the area between the apical surface of the cartilage and where the articular cartilage interfaces with subchondral bone. The technique was perfected in our laboratory and previously reported (Black et al., 2015). To help standardize the human measurement error, two different lab members independently measured the tissue samples.

Statistics

Statistical analysis was performed using a mixed models analysis of variance (ANOVA) on the OARSI or Mankin score as well as. The dependent variables were the OARSI, Mankin scores or biomarker staining for the knee. The independent variables were age and genotype along with their interaction. Blocking was done on each animal to account for their multiple measures. Post hoc t-tests were performed to determine differences in genotype biomarker expression level at each age, p-values of <0.05 were considered significant. The percent of chondrocyte biomarker expression level staining was compared using t-tests. A Bonferroni adjust was made for these four tests. Thus p-values of <0.0125 were considered significant.

Results

Co-expression of DDR-2 and Primary Cilia

We observed primary cilia (acetylated α-tubulin) in zebrafish chondrocytes residing in the mandible (Figure 1A) indicated in green staining. DDR-2, indicated by purple staining (Figure 1B), did not co-localize with primary cilia.

FIGURE 1

FIGURE 1. Normal primary cilia in zebrafish chondrocytes residing in the mandible, as indicated in green staining for acetylated α-tubulin (A). Ddr-2 is indicated by purple staining (B), demonstrating a lack of co-localization between primary cilia and Ddr-2.

Cartilage Morphology

In a histological comparison of 18-week-old WT and Bbs1^M390R/M390R mutant knee joints stained with Safranin-O, Fast Green we observed thicker cartilage (p < 0.01) in Bbs1^M390R/M390R mutant knees with noted surface erosion in some cases. See Figure 2 for representative pictures.

FIGURE 2

FIGURE 2. Histological comparison of WT and Bbs1^M390R/M390R mutant knee joints. (A) Knee joint stained with Safranin-O Fast Green from a 17-week-old male WT mouse. (B) Knee joint stained with Safranin-O, Fast Green from an 18-week-old male Bbs1^M390R/M390R mutant mouse. Note the surface erosion and exaggerated thickness of cartilage. Bbs1^M390R/M390R mutant knees exhibited significantly thicker cartilage compared to WT controls (p < 0.01). However, Bbs1^M390R/M390R mutant mice do not exhibit higher than normal OARSI or Mankin scores.

OARSI and Mankin Scores

Although Bbs1^M390R/M390R mutants tended to have higher OARSI scores than WT mice (Figure 3), the mixed models analysis showed no statistically significant difference was observed in OARSI or Mankin (data not shown) scoring between the two groups studied. The alteration from normal biomarker concentrations in Bbs1^M390R/M390R mutant mice, suggests that the metabolic changes occur prior to cartilage degradation.

FIGURE 3

FIGURE 3. OARSI score comparison of Bbs1^M390R/M390R mutant and WT mice at various ages. Articular cartilage degradation from both mutant and WT samples was objectively evaluated using the OARSI scoring system by two blind scorers. Samples were combined into generalized age groups of 10–14 (WT n = 8; Bbs1^M390R/M390R mutant n = 18), 15–18 (WT n = 24; Bbs1^M390R/M390R mutant n = 39) and 21+ (WT n = 12; Bbs1^M390R/M390R mutant n = 4) weeks. Mean OARSI scores of 2.25, 2.2, and 2.5 (10–14, 15–18, 21+ weeks, respectively) were not significantly different than those of WT mice.

MmMP-13, HTRA1, and TGF-β1 Expression in Bbs1^M390R/M390R Mutant Mice

T-tests of immunohistochemical staining revealed significant differences in the percent of chondrocytes staining positive for MMP-13 (p < 0.001) (Figure 4A vs. Figure 4E), HTRA1 (p < 0.001) (Figure 4B vs. Figure 4F), and TGF-β1 (p < 0.001) (Figure 4C vs. Figure 4G), but not DDR-2 (p > 0.05) (Figure 4D vs. Figure 4H) was not significant with the Bonferroni adjustment, in Bbs1^M390R/M390R mutant mice when compared to age matched WT mice. These findings remained consistent across all ages examined.

FIGURE 4

FIGURE 4. Representative images showing the results of the immunohistochemical and histological staining performed to analyze the presence of OA biomarkers. Although, Bbs1^M390R/M390R mutant and WT mouse knees did not exhibit a statistical difference in cartilage erosion scoring, we analyzed separately cell staining within the three major age (week) categories examined – 10–14 (WT n = 8; Bbs1^M390R/M390R mutant n = 18), 15–18 (WT n = 24; Bbs1^M390R/M390R mutant n = 39), 21+ (WT n = 12; Bbs1^M390R/M390R mutant n = 4). All tissue samples represented in this figure were age matched at 14 weeks. The patterns of cell staining for biomarkers examined did not differ significantly with age. (A) WT MMP-13; (B) WT HTRA-1; (C) WT TGF-β1; (D) WT DDR-2; (E) Bbs1^M390R/M390R mutant MMP13; (F) Bbs1^M390R/M390R mutant HTRA-1; (G) Bbs1^M390R/M390R mutant TGF-β1; (H) Bbs1^M390R/M390R mutant DDR-2.

When comparing all age groups, WT mice exhibited significantly less staining for HTRA1 and MMP-13 compared to Bbs1^M390R/M390R mutants. WT mice exhibited significantly more staining for TGF-β1 compared to Bbs1^M390R/M390R mutants, indicative of earlier stages of OA. This is consistent with previous observations of an inverse relationship between HRTA1 and TGF-β1 (Larkin et al., 2013). We did not observe a decrease in DDR2 expression in Bbs1^M390R/M390R mutant mice compared to WT controls (Figure 4D vs. Figure 4H), contrary to general observations of upregulation in OA cartilage compared to controls. Notwithstanding the verified catabolic role of DDR2, no protection was afforded Bbs1^M390R/M390R mutant cartilage by its decreased expression. Visualization of the primary cilia (acetylated α-tubulin) and DDR2 in embryonic mandibular zebrafish chondrocytes (Figure 1A) indicated that DDR2 did not co-localize with primary cilia. The data are summarized in Table 1.

TABLE 1

TABLE 1. Percent Bbs1^M390R/M390R mutant cells stained with biomarkers compared with WT.

Discussion

It is imperative to fully understand the metabolic pathway involved in the age-related onset and progression of OA if an effective treatment is to be developed. Common to the pathogenesis of human OA is the activation of the HTRA1-DDR2-MMP13 axis. HTRA1 is strongly expressed in the presence of stressors in murine OA models, and its activation augments the expression of DDR-2 (Polur et al., 2010). Upregulation of DDR2 results in excess binding of the receptor to its ligand, type II collagen (Leitinger et al., 2004), which in turn stimulates high levels of expression of the MMP-13 (Su et al., 2009). This ultimately leads to the destruction of articular cartilage. Aside from the wide array of evidence showing HTRA1-DDR2-MMP13 activity in OA of multiple modalities, the convergence of diverse noxious stimuli upon this axis is further supported in the protective effects exerted in DDR-2 hypomorphic strains of mice. When subjected to the destabilization of the medial meniscus (DMM) surgical procedure they demonstrate a significant decrease in the progression of OA compared to WT littermates (Xu et al., 2010). We aimed to demonstrate the significance of primary cilia as it pertains to the onset and progression of OA. It is noteworthy that ciliopathy related genes do not just function in cilia. For example, the BBSome plays a role in leptin receptor trafficking and insulin receptor trafficking to the plasma membrane in neurons in the hypothalamus. Thus, BBS genes may play a role in OA pathways even if not directly involving the primary cilia organelle. Curiously, WT mice consistently exhibited significantly more staining for DDR-2 than Bbs1^M390R/M390R mutants at all age observed age groups, inconsistent with previous observations correlating high levels of DDR-2 expression with more advanced OA. This may suggest a role for primary cilia in processing signals associated with DDR-2 activation and/or expression of the receptor.

In order to identify primary cilia as a component in the chondrocyte destruction pathway associated with OA, we used Bbs1^M390R/M390R mutant mice, which are known to have defective primary cilia (Figure 2B) (Mykytyn et al., 2004) and have been previously shown to exhibit OA-like phenotype (Kaushik et al., 2009). We compared Bbs1^M390R/M390R mutants to WT at different ages and in each instance found significant differences in the concentrations of key biomarkers previously shown to be associated with OA. Among these differences were increased concentrations of chondrocytes expressing HTRA-1 and a decreased concentration of chondrocytes expressing TGF-b1 and increased levels of MMP-13. All three observations are consistent with previous studies on biomarkers associated with OA (Murwantoko et al., 2004; Polur et al., 2010; Holt et al., 2012; Larkin et al., 2013). The increased presence of HTRA1 and MMP13 may help explain why the Bbs1^M390R/M390R mutants respond more vigorously to DMM surgery (Peterson et al., 2011).

TGFβ-1 functions largely in the SMAD signaling pathway; it is a cytokine involved in a number of cellular processes including cell growth, proliferation, differentiation, apoptosis, and the elaboration of extracellular matrix (van Beuningen, 2000). HTRA-1 has a degrading effect on the peri- and extracellular matrix, as well as an inhibitory effect on TGF-β1 (Oka et al., 2004; Tocharus et al., 2004). Thus HTRA-1 may contribute to OA not only through direct matrix destruction, but also by inhibiting possible protection to the matrix offered by TGF-β1. Expression of MMP-13, a protein with a degrading effect on extracellular matrix similar to HTRA-1, appears to be upregulated. This is consistent with alterations in MMP-13 expression as in OA cases resulting from other causes, such as aging, obesity, joint misalignment, and acute injury.

Our results indicate that age-related OA-associated expression of DDR-2 is altered in Bbs1^M390R/M390R mice. Curiously, DDR-2 expression appeared to be lower in Bbs1^M390R/M390R mutant mice than in WT mice at all ages. Whereas in cases of OA resulting from the aforementioned etiologies, DDR-2 expression is upregulated in mice presenting with OA. We believe a possible explanation for this could be that primary cilia might have a role in processing signals associated with DDR-2 activation and/or expression of the receptor. However, we have demonstrated that DDR2 does not co-localize to primary cilia (Figure 2A). Therefore, any interaction(s) between primary cilia and DDR2 would presumably be indirect.

We did not observe a statistical difference between Bbs1^M390R/M390R mutant and WT mice when using knee histological scoring methods (OARSI and Modified Mankin). We did observe a difference in the thickness of cartilage in Bbs1^M390R/M390R mutant mice compared to WT as others have previously noted (Kaushik et al., 2009).

Finally, much has been published regarding the potential role of inflammation in OA. More recently, the interaction of adipokines and inflammation in initiating OA has been demonstrated (Azamar-Llamas et al., 2017). Leptin, a peptide hormone involved in maintaining insulin sensitivity and contributing to the sensation of satiety, is expressed at very high levels in obese individuals. It appears to be correlated with OA as well, with intervention at the level of MMP-13 expression occurring. Downregulation of leptin mRNA translation via small interference RNA molecules inhibits MMP-13 expression in cultured osteoarthritic chondrocytes (Iliopoulos et al., 2007; Pallu et al., 2010). The situation appears to be most exacerbated in cases of extreme obesity where a strong positive correlation exists between the responsiveness of the MMP-13 gene to leptin and the BMI of osteoarthritic individuals. This is intriguing in light of the observation that primary cilia are intimately involved in leptin homeostasis (Seo et al., 2009; Guo et al., 2016).

Conclusion

Our findings support our initial hypothesis that BBS mutations are involved in the identified molecular pathway common to OA, with the exception of DDR-2. Since the BBSome plays a role in leptin receptor trafficking and insulin receptor trafficking to the plasma membrane in neurons in the hypothalamus, BBS genes may play a role in OA pathways even if not directly involving the primary cilia organelle. We have shown that MMP-13, HTRA-1, and TGF-β1 expression mimic their expression during all other investigated cases of OA. This suggests that primary cilia act as a surface receptor on chondrocytes and work to relay the condition of the cartilaginous matrix and its integrity to the associated chondrocyte. Primary cilia may also be an important link between inflammatory signals and OA. Future research should work to identify the specific role primary cilia play in the onset and progression of OA.

Strengths And Limitations

This paper adds comprehensive and important data to a paucity of information regarding the involvement of primary cilia in the pathogenesis of OA. Some significant and important findings in this work are that primary cilia are a link between pro-inflammatory bio-markers and OA and that the BBSome may play a role in this disease as well. While BBS is a pleiotropic disease with 21 known genes, the work associated with this paper employed the BBS mutant (Bbs1^M390R/M390R) which is associated with 25% of all human cases. Some weaknesses of this work include the lack of elucidating the exact molecular pathway(s) or link(s) between primary cilia and OA as well as no data from human BBS patients.

Institutional Safety

All U.S. National Institutes of Health guidelines for research involving Recombinant or Synthetic Nucleic Acid Molecules were followed. All chemical use was according to the U.S. Occupational Safety and Health Administration guidelines as overseen by the Brigham Young University Office of Risk Management.

Ethics Statement

This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved Institutional Animal Care and Use Committee (IACUC) protocols 5061426 and 6081818, of the University of Iowa and the Brigham Young University Institutional Animal Care and Use Committee protocol 13-0601. Animals were housed according to IACUC recommendations. Methods of euthanasia used were carbon dioxide inhalation followed by cervical dislocation, or anesthesia induced by ketamine/xylazine followed by transcardial perfusion. Humane endpoints were strictly observed, and every effort was made to minimize suffering.

Author Contributions

IS provided direct oversight for the project. MM, SG, DB, JC, BD, CW, BR, AS, MG, and JD performed all of the histological work including analyses associated therewith. DE was responsible for all statistical work assisted by MG. DE also assisted with experimental design. VS provided the BBS mice and assisted with experimental design. AS provided the zebrafish and oversaw all of that work, assisted by SG and JD. DK conceived the project and was the overall coordinator for it.

Funding

This work was supported, in part, by Brigham Young University Office of Research and Creative Activities Undergraduate and Mentoring Environment Grants. This work was supported in part by US National Institutes of Health grants R01 EY-011298 and RO1 EY-017168 (to VS) and the Roy J. Carver Charitable Trust.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors wish to thank Matthew Sterling for help in preparing this manuscript.

References

Ansley, S. J., Badano, J. L., Blacque, O. E., Hill, J., Hoskins, B. E., Leitch, C. C., et al. (2003). Basal body dysfunction is a likely cause of pleiotropic Bardet-Biedl syndrome. Nature 425, 628–633. doi: 10.1038/nature02030

PubMed Abstract | CrossRef Full Text | Google Scholar

Araki, Y., Tsuzuki, Wada T, Aizaki, Y., Sato, K., Yokota, K., Fujimoto, K., et al. (2016). Histone methylation and STAT-3 differentially regulate interleukin-6-induced matrix metalloproteinase gene activation in rheumatoid arthritis synovial fibroblasts. Arthritis Rheumatol. 68, 1111–1123. doi: 10.1002/art.39563

PubMed Abstract | CrossRef Full Text | Google Scholar

Azamar-Llamas, D., Hernández-Molina, G., Ramos-Ávalos, B., and Furuzawa-Carballeda, J. (2017). Adipokine contribution to the pathogenesis of osteoarthritis. Mediators Inflamm. 2017:5468023. doi: 10.1155/2017/5468023

PubMed Abstract | CrossRef Full Text | Google Scholar

Badano, J. L., Ansley, S. J., Leitch, C. C., Lewis, R. A., Lupski, J. R., and Katsanis, N. (2003). Identification of a novel Bardet-Biedl syndrome protein, BBS7, that shares structural features with BBS1 and BBS2. Am. J. Hum. Genet. 72, 650–658. doi: 10.1086/368204

PubMed Abstract | CrossRef Full Text | Google Scholar

Badano, J. L., Teslovich, T. M., and Katsanis, N. (2005). The centrosome in human genetic disease. Nat. Rev. Genet. 6, 194–205. doi: 10.1038/nrg1557

PubMed Abstract | CrossRef Full Text | Google Scholar

Black, C. S., Crepeau, P. K., Sheffield, I. D., Macdonald, J. R., Wootton, D. J., Baek, M., et al. (2015). Identification of the tidemark line of calcification in osteoarthritic cartilage using a stain for alkaline phosphatase. Osteoarthritis Cartilage 23:A139. doi: 10.1016/j.joca.2015.02.876

CrossRef Full Text | Google Scholar

Bomsta, B. D., Bridgewater, L. C., and Seegmiller, R. E. (2006). Premature osteoarthritis in the Disproportionate micromelia (Dmm) mouse. Osteoarthritis Cartilage 14, 477–485. doi: 10.1016/j.joca.2005.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiang, A. P., Beck, J. S., Yen, H. J., Tayeh, M. K., Scheetz, T. E., Swiderski, R. E., et al. (2006). Homozygosity mapping with SNP arrays identifies TRIM32, an E3 ubiquitin ligase, as a Bardet-Biedl syndrome gene (BBS11). Proc. Natl. Acad. Sci. U.S.A. 103, 6287–6292. doi: 10.1073/pnas.0600158103

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiang, A. P., Nishimura, D., Searby, C., Elbedour, K., Carmi, R., Ferguson, A. L., et al. (2004). Comparative genomic analysis identifies an ADP-ribosylation factor-like gene as the cause of Bardet-Biedl syndrome (BBS3). Am. J. Hum. Genet. 75, 475–484. doi: 10.1086/423903

PubMed Abstract | CrossRef Full Text | Google Scholar

Ezratty, E. J., Stokes, N., Chai, S., Shah, A. S., Williams, S. E., and Fuchs, E. (2011). A role for the primary cilium in Notch signaling and epidermal differentiation during skin development. Cell 145, 1129–1141. doi: 10.1016/j.cell.2011.05.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, Y., Esmail, M. A., Ansley, S. J., Blacque, O. E., Boroevich, K., Ross, A. J., et al. (2004). Mutations in a member of the Ras superfamily of small GTP-binding proteins causes Bardet-Biedl syndrome. Nat. Genet. 36, 989–993. doi: 10.1038/ng1414

PubMed Abstract | CrossRef Full Text | Google Scholar

Felson, D. T., Anderson, J. J., Naimark, A., Walker, A. M., and Meenan, R. F. (1988). Obesity and knee osteoarthritis. The Framingham Study. Ann. Intern. Med. 109, 18–24. doi: 10.7326/0003-4819-109-1-18

CrossRef Full Text | Google Scholar

Glasson, S. S., Chambers, M. G., Van, Den Berg WB, and Little, C. B. (2010). The OARSI histopathology initiative – recommendations for histological assessments of osteoarthritis in the mouse. Osteoarthritis Cartilage 18(Suppl. 3), S17–S23. doi: 10.1016/j.joca.2010.05.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Guccione, A. A. (1994). Arthritis and the process of disablement. Phys. Ther. 74, 408–414. doi: 10.1093/ptj/74.5.408

CrossRef Full Text | Google Scholar

Guo, D. F., Cui, H., Zhang, Q., Morgan, D. A., Thedens, D. R., Nishimura, D., et al. (2016). The BBSome controls energy homeostasis by mediating the transport of the leptin receptor to the plasma membrane. PLoS Genet. 12:e1005890. doi: 10.1371/journal.pgen.1005890

PubMed Abstract | CrossRef Full Text | Google Scholar

Hart, J. A. L. (1968). Cilia in articular cartilage. J. Anat. 103:222.

Google Scholar

Heon, E., Kim, G., Qin, S., Garrison, J. E., Tavares, E., Vincent, A., et al. (2016). Mutations in C8ORF37 cause Bardet Biedl syndrome (BBS21). Hum. Mol. Genet. 25, 2283–2294. doi: 10.1093/hmg/ddw096

PubMed Abstract | CrossRef Full Text | Google Scholar

Holt, D. W., Henderson, M. L., Stockdale, C. E., Farrell, J. T., Kooyman, D. L., Bridgewater, L. C., et al. (2012). Osteoarthritis-like changes in the heterozygous sedc mouse associated with the HtrA1-Ddr2-Mmp-13 degradative pathway: a new model of osteoarthritis. Osteoarthritis Cartilage 20, 430–439. doi: 10.1016/j.joca.2011.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Iliopoulos, D., Malizos, K. N., and Tsezou, A. (2007). Epigenetic regulation of leptin affects MMP-13 expression in osteoarthritic chondrocytes: possible molecular target for osteoarthritis therapeutic intervention. Ann. Rheum. Dis. 66, 1616–1621. doi: 10.1136/ard.2007.069377

PubMed Abstract | CrossRef Full Text | Google Scholar

Jensen, C. G., Poole, C. A., McGlashan, S. R., Marko, M., Issa, Z. I., Vujcich, K. V., et al. (2004). Ultrastructural, tomographic and confocal imaging of the chondrocyte primary cilium in situ. Cell Biol. Int. 28, 101–110. doi: 10.1016/j.cellbi.2003.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Ju, B. G., Solum, D., Song, E. J., Lee, K. J., Rose, D. W., Glass, C. K., et al. (2004). Activating the PARP-1 sensor component of the Groucho/TLE1 corepressor complex mediates a CaMKinase II delta-dependent neurogenic gene activation pathway. Cell 119, 815–829. doi: 10.1016/j.cell.2004.11.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Kageyama, R., Ohtsuka, T., and Kobayashi, T. (2007). The Hes gene family: repressors and oscillators that orchestrate embryogenesis. Development 134, 1243–1251. doi: 10.1242/dev.000786

PubMed Abstract | CrossRef Full Text | Google Scholar

Katsanis, N., Beales, P. L., Woods, M. O., Lewis, R. A., Green, J. S., Parfrey, P. S., et al. (2000). Mutations in MKKS cause obesity, retinal dystrophy and renal malformations associated with Bardet-Biedl syndrome. Nat. Genet. 26, 67–70. doi: 10.1038/79201

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaur, Y., de Souza, R. J., Gibson, W. T., and Meyre, D. (2017). A systematic review of genetic syndromes with obesity. Obes. Rev. 18, 603–634. doi: 10.1111/obr.12531

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaushik, A. P., Martin, J. A., Zhang, Q., Sheffield, V. C., and Morcuende, J. A. (2009). Cartilage abnormalities associated with defects of chondrocytic primary cilia in Bardet-Biedl syndrome mutant mice. J. Orthop. Res. 27, 1093–1099. doi: 10.1002/jor.20855

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J. C., Badano, J. L., Sibold, S., Esmail, M. A., Hill, J., Hoskins, B. E., et al. (2004). The Bardet-Biedl protein BBS4 targets cargo to the pericentriolar region and is required for microtubule anchoring and cell cycle progression. Nat. Genet. 36, 462–470. doi: 10.1038/ng1352

PubMed Abstract | CrossRef Full Text | Google Scholar

Kong, J. H., Yang, L., Dessaud, E., Chuang, K., Moore, D. M., Rohatgi, R., et al. (2015). Notch activity modulates the responsiveness of neural progenitors to sonic hedgehog signaling. Dev. Cell 33, 373–387. doi: 10.1016/j.devcel.2015.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Kouri, J. B., Argüello, C., Luna, J., and Mena, R. (1998). Use of microscopical techniques in the study of human chondrocytes from osteoarthritic cartilage: an overview. Microsc. Res. Tech. 40, 22–36. doi: 10.1002/(SICI)1097-0029(19980101)40:1<22::AID-JEMT4>3.0.CO;2-Z

PubMed Abstract | CrossRef Full Text | Google Scholar

Kouri, J. B., Jiménez, S. A., Quintero, M., and Chico, A. (1996). Ultrastructural study of chondrocytes from fibrillated and non-fibrillated human osteoarthritic cartilage. Osteoarthritis Cartilage 4, 111–125. doi: 10.1016/S1063-4584(05)80320-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Larkin, D. J., Kartchner, J. Z., Doxey, A. S., Hollis, W. R., Rees, J. L., Wilhelm, S. K., et al. (2013). Inflammatory markers associated with osteoarthritis after destabilization surgery in young mice with and without Receptor for Advanced Glycation End-products (RAGE). Front. Physiol. 4:121. doi: 10.3389/fphys.2013.00121

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawrence, R. C., Felson, D. T., Lawrence, R. C., Gabriel, S., Hirsch, R., Kwoh, C. K., et al. (2008). Estimates of the prevalence of arthritis and other rheumatic conditions in the United States. Part II. Arthritis Rheum. 58, 26–35. doi: 10.1002/art.23176

PubMed Abstract | CrossRef Full Text | Google Scholar

Leitinger, B., Steplewski, A., and Fertala, A. (2004). The D2 period of collagen II contains a specific binding site for the human discoidin domain receptor, DDR2. J. Mol. Biol. 344, 993–1003. doi: 10.1016/j.jmb.2004.09.089

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J. B., Gerdes, J. M., Haycraft, C. J., Fan, Y., Teslovich, T. M., May-Simera, H., et al. (2004). Comparative and basal genomics identifies a flagellar and basal body proteome that includes the BBS5 human disease gene. Cell 117, 541–552. doi: 10.1016/S0092-8674(04)00450-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Mankin, H. J., Dorfman, H., Lippiello, L., and Zarins, A. (1971). Biochemical and metabolic abnormalities in articular cartilage from osteo-arthritic human hips. II. Correlation of morphology with biochemical and metabolic data. J. Bone Joint Surg. Am. 53, 523–537. doi: 10.2106/00004623-197153030-00009

PubMed Abstract | CrossRef Full Text | Google Scholar

Marion, V., Schlicht, D., Mockel, A., Caillard, S., Imhoff, O., Stoetzel, C., et al. (2011). Bardet-Biedl syndrome highlights the major role of the primary cilium in efficient water reabsorption. Kidney Int. 79, 1013–1025. doi: 10.1038/ki.2010.538

PubMed Abstract | CrossRef Full Text | Google Scholar

McGlashan, S. R., Jensen, C. G., and Poole, C. A. (2006). Localization of extracellular matrix receptors on the chondrocyte primary cilium. J. Histochem. Cytochem. 54, 1005–1014. doi: 10.1369/jhc.5A6866.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Meier-Vismara, E., Walker, N., and Vogel, A. (1979). Single cilia in the articular cartilage of the cat. Exp. Cell Biol. 47, 161–171. doi: 10.1159/000162933

PubMed Abstract | CrossRef Full Text | Google Scholar

Murwantoko, Yano, M., Ueta, Y., Murasaki, A., Kanda, H., Oka, C., et al. (2004). Binding of proteins to the PDZ domain regulates proteolytic activity of HtrA1 serine protease. Biochem. J. 381(Pt 3), 895–904. doi: 10.1042/BJ20040435

PubMed Abstract | CrossRef Full Text | Google Scholar

Mykytyn, K., Braun, T., Carmi, R., Haider, N. B., Searby, C. C., Shastri, M., et al. (2001). Identification of the gene that, when mutated, causes the human obesity syndrome BBS4. Nat. Genet. 28, 188–191. doi: 10.1038/88925

PubMed Abstract | CrossRef Full Text | Google Scholar

Mykytyn, K., Mullins, R. F., Andrews, M., Chiang, A. P., Swiderski, R. E., Yang, B., et al. (2004). Bardet-Biedl syndrome type 4 (BBS4)-null mice implicate Bbs4 in flagella formation but not global cilia assembly. Proc. Natl. Acad. Sci. U.S.A. 101, 8664–8669. doi: 10.1073/pnas.0402354101

PubMed Abstract | CrossRef Full Text | Google Scholar

Mykytyn, K., and Sheffield, V. C. (2004). Establishing a connection between cilia and Bardet-Biedl Syndrome. Trends Mol. Med. 10, 106–109. doi: 10.1016/j.molmed.2004.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Nachury, M. V., Loktev, A. V., Zhang, Q., Westlake, C. J., Peränen, J., Merdes, A., et al. (2007). A core complex of BBS proteins cooperates with the GTPase Rab8 to promote ciliary membrane biogenesis. Cell 129, 1201–1213. doi: 10.1016/j.cell.2007.03.053

PubMed Abstract | CrossRef Full Text | Google Scholar

Niebler, S., Schubert, T., Hunziker, E. B., and Bosserhoff, A. K. (2015). Activating enhancer binding protein 2 epsilon (AP-2 epsilon)-deficient mice exhibit increased matrix metalloproteinase 13 expression and progressive osteoarthritis development. Arthritis Res. Ther. 17:119. doi: 10.1186/s13075-015-0648-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishimura, D. Y., Searby, C. C., Carmi, R., Elbedour, K., Van, Maldergem L, Fulton, A. B., et al. (2001). Positional cloning of a novel gene on chromosome 16q causing Bardet-Biedl syndrome (BBS2). Hum. Mol. Genet. 10, 865–874. doi: 10.1093/hmg/10.8.865

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishimura, D. Y., Swiderski, R. E., Searby, C. C., Berg, E. M., Ferguson, A. L., Hennekam, R., et al. (2005). Comparative genomics and gene expression analysis identifies BBS9, a new Bardet-Biedl syndrome gene. Am. J. Hum. Genet. 77, 1021–1033. doi: 10.1086/498323

PubMed Abstract | CrossRef Full Text | Google Scholar

Oka, C., Tsujimoto, R., Kajikawa, M., Koshiba-Takeuchi, K., Ina, J., Yano, M., et al. (2004). HtrA1 serine protease inhibits signaling mediated by Tgfbeta family proteins. Development 131, 1041–1053. doi: 10.1242/dev.00999

PubMed Abstract | CrossRef Full Text | Google Scholar

Pallu, S., Francin, P. J., Guillaume, C., Gegout-Pottie, P., Netter, P., Mainard, D., et al. (2010). Obesity affects the chondrocyte responsiveness to leptin in patients with osteoarthritis. Arthritis Res. Ther. 12:R112. doi: 10.1186/ar3048

PubMed Abstract | CrossRef Full Text | Google Scholar

Peterson, D. M., Weipert, D., Airmet, K. W., Killingsworth, S., Klein, M. R., Beardsley, T. D., et al. (2011). Premature osteoarthritis in the temporomandibular and knee joints of bardet-biedle syndrome mice. Osteoarthritis Cartilage 19:S62. doi: 10.1016/S1063-4584(11)60147-7

CrossRef Full Text | Google Scholar

Polur, I., Lee, P. L., Servais, J. M., Xu, L., and Li, Y. (2010). Role of HTRA1, a serine protease, in the progression of articular cartilage degeneration. Histol. Histopathol. 25, 599–608. doi: 10.14670/HH-25.599

PubMed Abstract | CrossRef Full Text | Google Scholar

Poole, C. A., Flint, M. H., and Beaumont, B. W. (1985). Analysis of the morphology and function of primary cilia in connective tissues: a cellular cybernetic probe? Cell Motil. 5, 175–193.

PubMed Abstract | Google Scholar

Poole, C. A., Jensen, C. G., Snyder, J. A., Gray, C. G., Hermanutz, V. L., and Wheatley, D. N. (1997). Confocal analysis of primary cilia structure and colocalization with the Golgi apparatus in chondrocytes and aortic smooth muscle cells. Cell Biol. Int. 21, 483–494. doi: 10.1006/cbir.1997.0177

PubMed Abstract | CrossRef Full Text | Google Scholar

Poole, C. A., Zhang, Z. J., and Ross, J. M. (2001). The differential distribution of acetylated and detyrosinated alpha-tubulin in the microtubular cytoskeleton and primary cilia of hyaline cartilage chondrocytes. J. Anat. 199(Pt 4), 393–405. doi: 10.1046/j.1469-7580.2001.19940393.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Scherft, J. P., and Daems, W. T. (1967). Single cilia in chondrocytes. J. Ultrastruct. Res. 19:546. doi: 10.1016/S0022-5320(67)80080-7

CrossRef Full Text | Google Scholar

Seo, S., Guo, D. F., Bugge, K., Morgan, D. A., Rahmouni, K., and Sheffield, V. C. (2009). Requirement of Bardet-Biedl syndrome proteins for leptin receptor signaling. Hum. Mol. Genet. 18, 1323–1331. doi: 10.1093/hmg/ddp031

PubMed Abstract | CrossRef Full Text | Google Scholar

Slavotinek, A. M., and Biesecker, L. G. (2000). Phenotypic overlap of McKusick-Kaufman syndrome with bardet-biedl syndrome: a literature review. Am. J. Med. Genet. 95, 208–215. doi: 10.1002/1096-8628(20001127)95:3<208::AID-AJMG5>3.0.CO;2-J

PubMed Abstract | CrossRef Full Text | Google Scholar

Spencer, K., Wilhelm, S. K., Kartchner, J. Z., Mecham, D., Reynolds, P. R., and Kooyman, D. L. (2015). Effect of pharmacological blocking of TLR-4 on osteoarthritis in mice. J. Arthritis 4:164. doi: 10.4172/2167-7921.1000164

CrossRef Full Text

Stoetzel, C., Laurier, V., Davis, E. E., Muller, J., Rix, S., Badano, J. L., et al. (2006a). BBS10 encodes a vertebrate-specific chaperonin-like protein and is a major BBS locus. Nat. Genet. 38, 521–524. doi: 10.1038/ng1771

PubMed Abstract | CrossRef Full Text | Google Scholar

Stoetzel, C., Laurier, V., Faivre, L., Mégarbané, A., Perrin-Schmitt, F., Verloes, A., et al. (2006b). BBS8 is rarely mutated in a cohort of 128 Bardet-Biedl syndrome families. J. Hum. Genet. 51, 81–84. doi: 10.1007/s10038-005-0320-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Stoetzel, C., Muller, J., Laurier, V., Davis, E. E., Zaghloul, N. A., Vicaire, S., et al. (2007). Identification of a novel BBS gene (BBS12) highlights the major role of a vertebrate-specific branch of chaperonin-related proteins in Bardet-Biedl syndrome. Am. J. Hum. Genet. 80, 1–11. doi: 10.1086/510256

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, J., Yu, J., Ren, T., Zhang, W., Zhang, Y., Liu, X., et al. (2009). Discoidin domain receptor 2 is associated with the increased expression of matrix metalloproteinase-13 in synovial fibroblasts of rheumatoid arthritis. Mol. Cell. Biochem. 330, 141–152. doi: 10.1007/s11010-009-0127-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Sugita, S., Hosaka, Y., Okada, K., Mori, D., Yano, F., Kobayashi, H., et al. (2015). Transcription factor Hes1 modulates osteoarthritis development in cooperation with calcium/calmodulin-dependent protein kinase 2. Proc. Natl. Acad. Sci. U.S.A. 112, 3080–3085. doi: 10.1073/pnas.1419699112

PubMed Abstract | CrossRef Full Text | Google Scholar

Suli, A., Guler, A. D., Raible, D. W., and Kimelman, D. (2014). A targeted gene expression system using the tryptophan repressor in zebrafish shows no silencing in subsequent generations. Development 141, 1167–1174. doi: 10.1242/dev.100057

PubMed Abstract | CrossRef Full Text | Google Scholar

Tocharus, J., Tsuchiya, A., Kajikawa, M., Ueta, Y., Oka, C., and Kawaichi, M. (2004). Developmentally regulated expression of mouse HtrA3 and its role as an inhibitor of TGF-beta signaling. Dev. Growth Differ. 46, 257–274. doi: 10.1111/j.1440-169X.2004.00743.x

PubMed Abstract | CrossRef Full Text | Google Scholar

van Beuningen H. M., Glansbeek, H. L., van der Kraan P. M., and van den Berg W. B. (2000). Osteoarthritis-like changes in the murine knee joint resulting from intra-articular transforming growth factor-beta injections. Osteoarthritis Cartilage 8, 25–33.

PubMed Abstract | Google Scholar

Weihbrecht, K., Goar, W. A., Pak, T., Garrison, J. E., DeLuca, A. P., Stone, E. M., et al. (2017). Keeping an eye on Bardet-Biedl syndrome: a comprehensive review of the role of Bardet-Biedl syndrome genes in the eye. Med. Res. Arch. 5, 1–21. doi: 10.18103/mra.v5i9.1526

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilsman, N. J. (1978). Cilia of adult canine articular chondrocytes. J. Ultrastruct. Res. 64, 270–281. doi: 10.1016/S0022-5320(78)90036-9

CrossRef Full Text | Google Scholar

Wilsman, N. J., Farnum, C. E., and Reedaksamit, D. K. (1980). Incidence and morphology of equine and murine chondrocytic cilia. Anat. Rec. 197, 355–361. doi: 10.1002/ar.1091970309

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, L., Flahiff, C. M., Waldman, B. A., Wu, D., Olsen, B. R., Setton, L. A., et al. (2003). Osteoarthritis-like changes and decreased mechanical function of articular cartilage in the joints of mice with the chondrodysplasia gene (cho). Arthritis Rheum. 48, 2509–2518. doi: 10.1002/art.11233

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, L., Servais, J., Polur, I., Kim, D., Lee, P. L., Chung, K., et al. (2010). Attenuation of osteoarthritis progression by reduction of the discoidin domain receptor 2 in mice. Arthritis Rheum. 62, 2736–2744. doi: 10.1002/art.27582

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y., Huang, Y., Cai, D., Liu, J., and Cao, X. (2015). Analysis of differences in the molecular mechanism of rheumatoid arthritis and osteoarthritis based on integration of gene expression profiles. Immunol. Lett. 168, 246–253. doi: 10.1016/j.imlet.2015.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Yammani, R. R., and Loeser, R. F. (2014). Stress-Inducible Nuclear Protein 1 Regulates Matrix Metalloproteinase 13 Expression in Human Articular Chondrocytes. Arthritis Rheumatol. 66, 1266–1271. doi: 10.1002/art.38391

PubMed Abstract | CrossRef Full Text | Google Scholar