To further explore the function of the six identified SE in human brains, we designed a computational strategy to identify their downstream genes. For each SE, we first predicted 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 }}}}_{SE}$$\end{document}Ψ^SE using the genotypes in the transcribed region of the host gene from all available CMC samples, and then we stratified the individuals into two groups based on low and high \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 }}}}_{SE}$$\end{document}Ψ^SE.levels (Supplementary Fig. S2A–E). The splicing event inTBC1D5 was not analyzed because most samples clustered in the central region of all \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}Ψ^ values; thus, the number of remaining samples that could be stratified as high or low was insufficient for statistical analysis (Supplementary Fig. S2F). We identified the differentially expressed genes between the two groups for each SE independently. The number of the downstream differentially expressed genes identified for each SE ranged from 5 for LINC00665 to 471 for DRC1 (Supplementary Fig. S3). In total, 970
