Brain proteome‐wide association study linking‐genes in multiple sclerosis pathogenesis

OBJECTIVES: To identify genes that confer MS risk via the alteration of cis‐regulated protein abundance and verify their aberrant expression in human brain. METHODS: Utilizing a two‐stage proteome‐wide association study (PWAS) design, MS GWAS data (N = 41,505) was respectively integrated with two distinct human brain proteomes from the dorsolateral prefrontal cortex, including ROSMAP (N = 376) in the discovery stage and Banner (N = 152) in the confirmation stage. In the following, Bayesian colocalization analysis was conducted for GWAS and protein quantitative trait loci signals to prioritize candidate genes. Differential expression analysis was then used to verify the dysregulation of risk genes in white matter and gray matter for evidence at the transcription level. RESULTS: A total of 51 genes whose protein abundance had association with the MS risk were identified, of which 18 genes overlapped in the discovery and confirmation PWAS. Bayesian colocalization indicated six causal genes with genetic risk variants for the MS risk. The differential expression analysis of SHMT1 (P (FDR) = 4.82 × 10(−2)), FAM120B (P (FDR) = 8.13 × 10(−4)) in white matter and ICA1L (P (FDR) = 3.44 × 10(−2)) in gray matter confirmed the dysregulation at the transcription level. Further investigation of expression found SHMT1 significantly up‐regulated in white matter lesion, and FAM120B up‐regulated in both white matter lesion and normal appearing white matter. ICA1L was down‐regulated in both gray matter lesion and normal appearing gray matter. INTERPRETATION: Dysregulation of SHMT1, FAM120B and ICA1L may confer MS risk. Our findings shed new light on the pathogenesis of MS and prioritized promising targets for future therapy research.