Superficial zone protein (SZP), also known as lubricin and proteoglycan 4 (PRG4), plays an important role in the boundary lubrication of articular cartilage and is regulated by transforming growth factor (TGF)-β. Here, we evaluate the role of cell surface glycosaminoglycans (GAGs) during TGF-β1 stimulation of SZP/lubricin/PRG4 in superficial zone articular chondrocytes. We utilized primary monolayer superficial zone articular chondrocyte cultures and treated them with various concentrations of TGF-β1, in the presence or absence of heparan sulfate (HS), heparin, and chondroitin sulfate (CS). The cell surface GAGs were removed by pretreatment with either heparinase I or chondroitinase-ABC before TGF-β1 stimulation. Accumulation of SZP/lubricin/PRG4 in the culture medium in response to stimulation with TGF-β1 and various exogenous GAGs was demonstrated by immunoblotting and quantitated by enzyme-linked immunosorbent assay. We show that TGF-β1 and exogenous HS enhanced SZP accumulation of superficial zone chondrocytes in the presence of surface GAGs. At the dose of 1 ng/mL of TGF-β1, the presence of exogenous heparin inhibited SZP accumulation whereas the presence of exogenous CS stimulated SZP accumulation in the culture medium. Enzymatic depletion of GAGs on the surface of superficial zone chondrocytes enhanced the ability of TGF-β1 to stimulate SZP accumulation in the presence of both exogenous heparin and CS. Collectively, these results suggest that GAGs at the surface of superficial zone articular chondrocytes influence the response to TGF-β1 and exogenous GAGs to stimulate SZP accumulation. Cell surface GAGs modulate superficial zone chondrocytes' response to TGF-β1 and exogenous HS.

The expression of bone morphogenetic proteins (BMPs) and their cognate receptors (BMPRs) in osteochondromas has not been investigated. We determined the immunohistochemical localization and distribution of BMP-2/4, -6 and -7; BMP receptors BMPR-1A, BMPR-1B and BMPR-2; signal transducing proteins phosphorylated Smad1/5/8; and BMP antagonist noggin in the cartilaginous cap of solitary (SO) and multiple (MO) human osteochondromas and compared these with bovine growth plate and articular cartilage. The distribution and localization patterns for BMP-6, BMP-7, BMPR-1A and BMPR-2 were similar between the cartilaginous cap and the growth plate. BMP-2/4 and BMPR-1B were present throughout the growth plate. However, BMP-2/4 and phosphorylated Smad1/5/8 were mainly detected in proliferating chondrocytes of the cartilaginous cap. Also, BMPR-1B was found in hypertrophic chondrocytes of SO and proliferating chondrocytes of MO. Noggin was observed in resting chondrocytes and, to a lesser extent, in clustered proliferating chondrocytes in SO. On the other hand, noggin in MO was observed in proliferating chondrocytes. Since BMPs can stimulate proliferation and hypertrophic differentiation of chondrocytes, these findings suggest that there is an imbalance of BMP-2/4 and noggin interactions that may lead to abnormal regulation of chondrocyte proliferation and differentiation in the cartilaginous cap of human osteochondromas.

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

Our previous genome-wide association study (GWAS) for sagittal nonsyndromic craniosynostosis (sNCS) provided important insights into the genetics of midline CS. In this study, we performed a GWAS for a second midline NCS, metopic NCS (mNCS), using 215 non-Hispanic white case-parent triads. We identified six variants with genome-wide significance (P ≤ 5 × 10–8): rs781716 (P = 4.71 × 10–9; odds ratio [OR] = 2.44) intronic to SPRY3; rs6127972 (P = 4.41 × 10–8; OR = 2.17) intronic to BMP7; rs62590971 (P = 6.22 × 10–9; OR = 0.34), located ~ 155 kb upstream from TGIF2LX; and rs2522623, rs2573826, and rs2754857, all intronic to PCDH11X (P = 1.76 × 10–8, OR = 0.45; P = 3.31 × 10–8, OR = 0.45; P = 1.09 × 10–8, OR = 0.44, respectively). We performed a replication study of these variants using an independent non-Hispanic white sample of 194 unrelated mNCS cases and 333 unaffected controls; only the association for rs6127972 (P = 0.004, OR = 1.45; meta-analysis P = 1.27 × 10–8, OR = 1.74) was replicated. Our meta-analysis examining single nucleotide polymorphisms common to both our mNCS and sNCS studies showed the strongest association for rs6127972 (P = 1.16 × 10–6). Our imputation analysis identified a linkage disequilibrium block encompassing rs6127972, which contained an enhancer overlapping a CTCF transcription factor binding site (chr20:55,798,821–55,798,917) that was significantly hypomethylated in mesenchymal stem cells derived from fused metopic compared to open sutures from the same probands. This study provides additional insights into genetic factors in midline CS.