Leveraging single cell multiomic analyses to identify gene regulatory networks that drive human articular cartilage cell fate.

Compromises to the integrity of articular cartilage, whether genetic, post-traumatic or age-related, lead to debilitating chronic degenerative diseases including osteoarthritis. In this study, our objective was to leverage multiomic and spatial transcriptomic datasets to identify gene regulatory networks that drive human articular cartilage cell fate, and build a foundation from which novel therapeutics to overcome degenerative disease could be developed. We jointly profiled the transcriptome and open chromatin regions in individual nuclei recovered from distal femora at 2 fetal timepoints and performed high-definition spatial transcriptomics at an additional timepoint. We established a human pluripotent stem cell platform to interrogate the function of computationally predicted transcription factors during human chondrocyte differentiation. We computationally predicted gene regulatory networks governing chondrocyte subsets comprising the human distal femur during development. Following functional analysis of two transcription factors predicted to function in the superficial zone, CREB5 and NFATC2, using our in vitro experimental platform, we found both have the potential to reprogram growth plate cartilage and induce features of articular cartilage. We further identified new biological roles for CREB5 related to ECM organization and taxis. We expect new regulatory networks we uncovered to be important for promoting cartilage health and treating disease, and our platform to be a useful tool for studying cartilage development and homeostasis in vitro. The ability to reprogram chondrocytes toward an articular-like fate has significant therapeutic potential to treat degenerative joint diseases.