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 formation in 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.

Muscle atrophy and fiber type alterations are well-characterized physiological adaptations to microgravity with both understood to be primarily regulated by differential gene expression (DGE). While microgravity-induced DGE has been extensively investigated, adaptations to microgravity due to alternative splicing (AS) have not been studied in a mammalian model. We sought to comprehensively elucidate the transcriptomic underpinnings of microgravity-induced muscle phenotypes in mice by evaluating both DGE and changes in AS due to extended spaceflight. Tissue sections and total RNA were isolated from the gastrocnemius and quadriceps, postural and phasic muscles of the hind limb, respectively, of 32-week-old female BALB/c mice exposed to microgravity or ground control conditions for nine weeks. Immunohistochemistry disclosed muscle type-specific physiological adaptations to microgravity that included i) a pronounced reduction in muscle fiber cross-sectional area in both muscles and ii) a prominent slow-to-fast fiber type transition in the gastrocnemius. RNA sequencing revealed that DGE and AS varied across postural and phasic muscle types with preferential employment of DGE in the gastrocnemius and AS in the quadriceps. Gene ontology analysis indicated that DGE and AS regulate distinct molecular processes. Various non-differentially expressed transcripts encoding musculoskeletal proteins (Tnnt3, Tnnt1, Neb, Ryr1, and Ttn) and muscle-specific RNA binding splicing regulators (Mbnl1 and Rbfox1) were found to have significant changes in AS that altered critical functional domains of their protein products. In striking contrast, microgravity-induced differentially expressed genes were associated with lipid metabolism and mitochondrial function. Our work serves as the first comprehensive investigation of coordinate changes in DGE and AS in large limb muscles across spaceflight. We propose that substantial remodeling of pre-mRNA by AS is a major component of transcriptomic adaptation of skeletal muscle to microgravity. The alternatively spliced genes identified here could be targeted by small molecule splicing regulator therapies to address microgravity-induced changes in muscle during spaceflight.