Background As the interest in manned spaceflight increases, so does the requirement to understand the transcriptomic mechanisms that underlay the detrimental physiological adaptations of skeletal muscle to microgravity. While microgravity-induced differential gene expression (DGE) has been extensively investigated, the contribution of differential alternative splicing (DAS) to the plasticity and functional status of the skeletal muscle transcriptome has not been studied in an animal model. Therefore, by evaluating both DGE and DAS across spaceflight, we set out to provide the first comprehensive characterization of the transcriptomic landscape of skeletal muscle during exposure to microgravity. Methods RNA-sequencing, immunohistochemistry, and morphological analyses were conducted utilizing total RNA and tissue sections isolated from the gastrocnemius and quadriceps muscles of 30-week-old female BALB/c mice exposed to microgravity or ground control conditions for 9 weeks. Results In response to microgravity, the skeletal muscle transcriptome was remodeled via both DGE and DAS. Importantly, while DGE showed variable gene network enrichment, DAS was enriched in structural and functional gene networks of skeletal muscle, resulting in the expression of alternatively spliced transcript isoforms that have been associated with the physiological changes to skeletal muscle in microgravity, including muscle atrophy and altered fiber type function. Finally, RNA-binding proteins, which are required for regulation of pre-mRNA splicing, were themselves differentially spliced but not differentially expressed, an upstream event that is speculated to account for the downstream splicing changes identified in target skeletal muscle genes. Conclusions Our work serves as the first investigation of coordinate changes in DGE and DAS in large limb muscles across spaceflight. It opens up a new opportunity to understand (i) the molecular mechanisms by which splice variants of skeletal muscle genes regulate the physiological adaptations of skeletal muscle to microgravity and (ii) how small molecule splicing regulator therapies might thwart muscle atrophy and alterations to fiber type function during prolonged spaceflight.

Vitamin D is known to modulate human immune responses, and vitamin D deficiency is associated with increased susceptibility to infection. However, what constitutes sufficient levels or whether vitamin D is useful as an adjuvant therapeutic is debated, much in part because of inadequate elucidation of mechanisms underlying vitamin D’s immune modulatory function. Cathelicidin antimicrobial peptide (CAMP) has potent broad-spectrum activity, and the CAMP gene is regulated in human innate immune cells by active 1,25(OH)2D3, a product of hydroxylation of inactive 25(OH)D3 by CYP27B1-hydroxylase. We developed a CRISPR/Cas9-edited human monocyte-macrophage cell line containing the mCherry fluorescent reporter gene at the 3′ end of the endogenous CAMP gene. The High Throughput CAMP Assay (HiTCA) developed here is a novel tool for evaluating CAMP expression in a stable cell line that is scalable for a high-throughput workflow. Application of HiTCA to serum samples from a small number of human donors (n = 10) showed individual differences in CAMP induction that were not fully accounted for by the serum vitamin D metabolite status of the host. As such, HiTCA may be a useful tool that can advance our understanding of the human vitamin D-dependent antimicrobial response, which is being increasingly appreciated for its complexity.

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.