Bipolar disorder (BD) presents a significant challenge due to its chronic and relapsing nature, with its underlying pathogenesis remaining elusive. This study employs Mendelian randomization (MR), a widely recognized genetic approach, to unveil intricate causal associations between proteins and BD, leveraging protein quantitative trait loci (pQTL) as key exposures. We integrate pQTL data from brain, cerebrospinal fluid (CSF), and plasma with genome-wide association study (GWAS) findings of BD within a comprehensive systems analysis framework. Our analyses, including two-sample MR, Steiger filtering, and Bayesian colocalization, reveal noteworthy associations. Elevated levels of AGRP, FRZB, and IL36A in CSF exhibit significant associations with increased BD_ALL risk, while heightened levels of CTSF and LRP8 in CSF, and FLRT3 in plasma, correlate with decreased BD_ALL risk. Specifically for Bipolar I disorder (BD_I), increased CSF AGRP levels are significantly linked to heightened BD_I risk, whereas elevated CSF levels of CTSF and LRP8, and plasma FLRT3, are associated with reduced BD_I risk. Notably, genes linked to BD-related proteins demonstrate substantial enrichment in functional pathways such as “antigen processing and presentation,” “metabolic regulation,” and “regulation of myeloid cell differentiation.” In conclusion, our findings provide beneficial evidence to support the potential causal relationship between IL36A, AGRP, FRZB, LRP8 in cerebrospinal fluid, and FLRT3 in plasma, and BD and BD_I, providing insights for future mechanistic studies and therapeutic development.