Elucidating dynamics and regulation of alternative splicing in osteogenic differentiation

Nearly all human multi-exonic genes undergo alternative splicing (AS) via regulation by RNA-binding proteins (RBPs), but few studies have examined the temporal dynamics of AS and its regulation during cell differentiation in the bone niche. We sought to evaluate how AS, under the control of RBPs, affects cell fate commitment during induced osteogenic differentiation of human bone marrow-derived multipotent stem/stromal progenitor cells (MSPCs). We generated a time-course RNA sequencing (RNA-seq) dataset representative of induced MSPC differentiation to osteoblasts. Our analysis revealed widespread AS changes, coordinated with differential RBP expression, at multiple time points, including many AS changes in non-differentially expressed genes. We also developed a computational approach to profile the dynamics and regulation of AS by RBPs using time-course RNA-seq data, by combining temporal patterns of exon skipping and RBP expression with RBP binding sites in the vicinity of regulated exons. In total we identified nine RBPs as potential key splicing regulators during MSPC osteogenic differentiation. Perturbation of one candidate,KHDRBS3, inhibited osteogenesis and bone formationin vitro, validating our computational prediction of “driver” RBPs. Overall, our work highlights a high degree of complexity in the splicing regulation of MSPC osteogenic differentiation. Our computational approach may be applied to other time-course data to explore dynamic AS changes and associated regulatory mechanisms in other biological processes or disease trajectories.