We used the R package ClusterProfiler (version 4.0.5) [35] to perform enrichment analysis for the differentially expressed genes based on Gene Ontology (GO) biological processes and molecular functions. The enrichment significance threshold was FDR < 0.05. We further explored the functions of the skipped exons by Gene Set Enrichment Analysis (GSEA) [36] through ClusterProfiler using three pathway knowledgebases GO [37], KEGG [38] and Hallmarks [39]. An enrichment score was computed for each pathway to determine if it was enriched or depleted based on changes of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\hat{\mathrm \Psi }}}}$$\end{document}Ψ^ levels. The significance threshold was also FDR < 0.05.
