Abstract
Early diagnosis and biomarker discovery to bolster the therapeutic pipeline for Parkinson's disease (PD) are urgently needed. In this study, we leverage the large-scale whole-blood total RNA-seq dataset from the Accelerating Medicine Partnership in Parkinson's Disease (AMP PD) program to identify PD-associated RNAs, including both known genes and novel circular RNAs (circRNA) and enhancer RNAs (eRNAs). There were 1,111 significant marker RNAs, including 491 genes, 599 eRNAs, and 21 circRNAs, that were first discovered in the PPMI cohort (FDR < 0.05) and confirmed in the PDBP/BioFIND cohorts (nominal p < 0.05). Functional enrichment analysis showed that the PD-associated genes are involved in neutrophil activation and degranulation, as well as the TNF-alpha signaling pathway. We further compare the PD-associated genes in blood with those in post-mortem brain dopamine neurons in our BRAINcode cohort. 44 genes show significant changes with the same direction in both PD brain neurons and PD blood, including neuroinflammation-associated genes IKBIP, CXCR2, and NFKBIB. Finally, we built a novel multi-omics machine learning model to predict PD diagnosis with high performance (AUC = 0.89), which was superior to previous studies and might aid the decision-making for PD diagnosis in clinical practice. In summary, this study delineates a wide spectrum of the known and novel RNAs linked to PD and are detectable in circulating blood cells in a harmonized, large-scale dataset. It provides a generally useful computational framework for further biomarker development and early disease prediction.