Abstract
Alzheimer's Disease (AD) is marked by the accumulation of pathology, neuronal loss, and gliosis and frequently accompanied by cognitive decline. Understanding brain cell interactions is key to identifying new therapeutic targets to slow its progression. Here, we used systems biology methods to analyze single-nucleus RNA sequencing (snRNASeq) data generated from dorsolateral prefrontal cortex (DLPFC) tissues of 424 participants in the Religious Orders Study or the Rush Memory and Aging Project (ROSMAP). We identified modules of co-regulated genes in seven major cell types, assigned them to coherent cellular processes, and assessed which modules were associated with AD traits such as cognitive decline, tangle density, and amyloid-β deposition. Coexpression network structure was conserved in the majority of modules across cell types, but we also found distinct communities with altered connectivity, especially when compared to bulk RNASeq, suggesting cell-specific gene co-regulation. These coexpression modules can also capture signatures of cell subpopulations and be influenced by cell proportions. Using a Bayesian network framework, we modeled the direction of relationships between the modules and AD progression. We highlight two key modules, a microglia module (mic_M46), associated with tangles; and an astrocyte module (ast_M19), associated with cognitive decline. Our work provides cell-specific molecular networks modeling the molecular events leading to AD.