Microglia specific alternative splicing alterations in multiple sclerosis

Several aberrant alternative splicing (AS) events and their regulatory mechanisms are widely recognized in multiple sclerosis (MS). Yet the cell-type specific AS events have not been extensively examined. Here we assessed the diversity of AS events using web-based RNA-seq data of sorted CD15^-CD11b⁺ microglia in white matter (WM) region from 10 patients with MS and 11 control subjects. The GSE111972 dataset was downloaded from GEO and ENA databases, aligned to the GRCh38 reference genome from ENSEMBL via STAR. rMATS was used to assess five types of AS events, alternative 3'SS (A3SS), alternative 5'SS (A5SS), skipped exon (SE), retained intron (RI) and mutually exclusive exons (MXE), followed by visualizing with rmats2sashimiplot and maser. Differential genes or transcripts were analyzed using the limma R package. Gene ontology (GO) analysis was performed with the clusterProfiler R package. 42,663 raw counts of AS events were identified and 132 significant AS events were retained based on the filtered criteria: 1) average coverage >10 and 2) delta percent spliced in (ΔPSI) >0.1. SE was the most common AS event (36.36%), followed by MXE events (32.58%), and RI (18.94%). Genes related to telomere maintenance and organization primarily underwent SE splicing, while genes associated with protein folding and mitochondrion organization were predominantly spliced in the MXE pattern. Conversely, genes experiencing RI were enriched in immune response and immunoglobulin production. In conclusion, we identified microglia-specific AS changes in the white matter of MS patients, which may shed light on novel pathological mechanisms underlying MS.