| 【Abstract】 Background Anxiety disorder represents a prevalent psychiatric condition with a persistently rising global prevalence, significantly compromising patients' quality of life and social functioning. Currently, there is a clinical scarcity of objective and reliable biomarkers for early diagnosis and treatment evaluation of anxiety disorders. Recent research suggests that plasma proteins may serve as potential biomarkers for psychiatric diseases; however, the causal relationship between specific plasma proteins and anxiety disorders remains unclear. Objective This study aims to identify plasma proteins causally associated with anxiety disorders and elucidate their underlying biological pathways, thereby providing insights for the discovery of biomarkers and potential therapeutic targets for this condition. Methods This study utilized protein quantitative trait locus (pQTL) data for 4,907 plasma proteins derived from the deCODE study comprising 35,559 Icelandic individuals, and genome-wide association study (GWAS) data encompassing 50,486 anxiety disorder cases and 330,460 healthy controls. Bidirectional Mendelian randomization analyses were conducted using inverse variance weighting (IVW) as the primary method, supplemented by MR-Egger, weighted median, simple mode, and weighted mode approaches.. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted for the identified proteins. Sensitivity analyses, including Cochran’s Q test, MR-Egger intercept test, MR-PRESSO global test, and leave-one-out analysis, were performed to assess the robustness of the findings. Results Ten plasma proteins were identified as significantly associated with anxiety disorders. Among these, SPATA9 (OR=0.856, 95% CI: 0.784~0.934, P<0.001) and PDE5A (OR=0.911, 95% CI: 0.864~0.961, P<0.001) were identified as protective factors, while CRYGD (OR=1.209, 95% CI: 1.095~1.334, P<0.001), BTN3A3 (OR=1.045, 95% CI: 1.018~1.073, P<0.001), SERPINB13 (OR=1.102, 95% CI: 1.040~1.168, P<0.001), ERBB4 (OR=1.283, 95% CI: 1.109~1.484, P<0.001), LSAMP (OR=1.096, 95% CI: 1.037~1.158, P<0.005), ICOSLG (OR=1.283, 95% CI: 1.104~1.490, P<0.005), DNAJB11 (OR=1.172, 95% CI: 1.076~1.277, P<0.001), and TREML1 (OR=1.115, 95% CI: 1.054~1.179, P<0.001) were identified as risk factors. Sensitivity analyses confirmed the robustness of the results, with no evidence of heterogeneity (Cochran’s Q test, P>0.05) or pleiotropy (MR–Egger intercept test, P>0.05). Enrichment analysis indicated that these proteins were significantly enriched in biological processes such as T-cell signaling, lymphocyte proliferation, cell membrane structure, and synaptic function, as well as in the intestinal immune network for IgA production and the ErbB signaling pathway. Conclusion This study identified ten plasma proteins associated with the risk of anxiety disorders through MR analysis. These proteins are involved in multiple biological processes, including neurodevelopment and immune regulation, providing evidence for potential biomarkers and therapeutic targets for anxiety disorders. |