Fig. 4

Regulation of m6A in alternative splicing during muscle fiber conversion. A The fraction of mRNA with alternative splicing events in EDL and SOL. B Cumulative frequency of log2FC (MFPKM) for m6A modified genes with alternative splicing. High (AS Counts ≥ upper quartile; medium (lower quartile < AS Counts < upper quartile); low (AS Counts ≤ lower quartile). C The distribution of SE-DAS sites. D Four-quadrant graph showing the distribution of skipping exons with a marked change in both PSI and m6A methylation levels in SOL and EDL. E Genome browser tracks showing MeRIP-seq (red) and RNA-seq (light-blue) data at gene loci in SOL and EDL. m6A sites was labeled with red box. F RBP binding motif enriched in the skipping exon or neighboring regions. G Expression heatmap of RBPs. Differentially expressed gene is labeled in red. H Chord plot illustrating the GO biological process terms enriched for SE-DMAS genes. Genes contributing to their respective enrichment are shown on the left, and enriched GO terms are shown on the right