Single-cell multiomics reveals disrupted glial gene regulatory programs in Alzheimer’s disease via interpretable machine learning

Recent development of single-cell technology across multiple omics platforms has provided new ways to obtain holistic views of cells to study disease pathobiology. Alzheimer’s disease (AD) is the most common form of dementia worldwide, yet the detailed understanding of its cellular and molecular mechanisms remains limited. In this study, we analyzed paired single-cell transcriptomic (scRNA-seq) and chromatin accessibility (scATAC-seq) data from the Seattle Alzheimer’s Disease Brain Cell Atlas (SEA-AD) Consortium to investigate the molecular mechanisms of AD at a cell-subpopulation-specific resolution focusing on glial cells. We benchmarked various multi-omics integration methods using diverse metrics and built an analytic workflow that enabled effective batch correction and cross-modality alignment, creating a unified cell state space. Through integrative analysis of 26 human brain samples, we uncovered AD-associated gene expression and pathway changes in glial subpopulations and highlighted important transcriptomic and epigenomic signatures via functional inference and interpretable machine learning paradigms, discovering the profound involvement of the Solute Carrier proteins (SLC) family genes in multiple glial cell types. We also identified glial cell-specific regulatory programs mediated by key transcription factors such as JUN and FOSL2 in astrocytes, the Zinc Finger (ZNF) family genes in microglia, and the SOX family of transcription factors in oligodendrocytes. Our study provides a comprehensive workflow and a high-resolution view of how glial regulatory programs are disrupted in AD. Our findings offer novel insights into disease-related changes in gene regulation and suggest potential targets for further research and therapy.