Integrating single-nucleus sequence profiling to reveal the transcriptional dynamics of Alzheimer’s disease, Parkinson’s disease, and multiple sclerosis

Background Alzheimer's disease (AD), Parkinson's disease (PD), and multiple sclerosis (MS) are three nervous system diseases that partially overlap clinically and genetically. However, bulk RNA-sequencing did not accurately detect the core pathogenic molecules in them. The availability of high-quality single cell RNA-sequencing data of post-mortem brain collections permits the generation of a large-scale gene expression in different cells in human brain, focusing on the molecular features and relationships between diseases and genes. We integrated single-nucleus RNA-sequencing (snRNA-seq) datasets of human brains with AD, PD, and MS to identify transcriptomic commonalities and distinctions among them. Methods The snRNA-seq datasets were downloaded from Gene Expression Omnibus (GEO) database. The Seurat package was used for snRNA-seq data processing. The uniform manifold approximation and projection (UMAP) were utilized for cluster identification. The FindMarker function in Seurat was used to identify the differently expressed genes. Functional enrichment analysis was carried out using the Gene Set Enrichment Analysis (GSEA) and Gene ontology (GO). The protein‐protein interaction (PPI) analysis of differentially expressed genes (DEGs) was analyzed using STRING database (http://string-db.org). SCENIC analysis was performed using utilizing pySCENIC (v0.10.0) based on the hg19-tss-centered-10 kb-10species databases. The analysis of potential therapeutic drugs was analyzed on Connectivity Map (https://clue.io). Results The gene regulatory network analysis identified several hub genes regulated in AD, PD, and MS, in which HSPB1 and HSPA1A were key molecules. These upregulated HSP family genes interact with ribosome genes in AD and MS, and with immunomodulatory genes in PD. We further identified several transcriptional regulators (SPI1, CEBPA, TFE3, GRHPR, and TP53) of the hub genes, which has important implications for uncovering the molecular crosstalk among AD, PD, and MS. Arctigenin was identified as a potential therapeutic drug for AD, PD, and MS. Conclusions Together, the integrated snRNA-seq data and findings have significant implications for unraveling the shared and unique molecular crosstalk among AD, PD, and MS. HSPB1 and HSPA1A as promising targets involved in the pathological mechanisms of neurodegenerative diseases. Additionally, the identification of arctigenin as a potential therapeutic drug for AD, PD, and MS further highlights its potential in treating these neurological disorders. These discoveries lay the groundwork for future research and interventions to enhance our understanding and treatment of AD, PD, and MS. Supplementary Information The online version contains supplementary material available at 10.1186/s12967-023-04516-6.


Introduction
In face of aging world population and the absence of effective cures, central nervous system (CNS) diseases pose a significant economic burden on society.Alzheimer's disease (AD), Parkinson's disease (PD), and multiple sclerosis (MS) are chronic and progressive diseases of the CNS, characterized by the loss of neurons in the brain or the spinal cord [1].Although a variety of efforts at the molecular level have attempted to elucidate the basic biological pathologies contributing to these diseases, the principal causes and cures for these diseases remain elusive, which hinders the discovery of disease-modifying therapies.
AD is the most common neurodegenerative disease, identified by the presence of extracellular aggregates of amyloid β (Aβ) peptides and intraneuronal tau neurofibrillary tangles in the brain.Patients with AD primarily exhibit the impairments in short-term memory and cognitive processing [2].PD is the second most common neurodegenerative disease, and is clinically defined by the symptoms of akinesia, rigidity, and tremor.Cytosolic Lewy bodies aggregated by α-synuclein and the loss of dopaminergic neurons in the substantia nigra pars compacta are the major neuropathological features [3].MS is a progressive neuroinflammatory disease with distinct lesion characteristics in the cortical grey versus subcortical white matter, and neurodegeneration at chronic stages.Clinically, the symptoms of MS are diverse, involving the impairments of movement, sensation, vision, and cognition [4].Although AD, PD, and MS differ in many clinical and pathological aspects, it is possible that they share cross-molecular features.Growing evidence indicates that brain endothelial dysfunction might play a significant role in the neurobiology of AD, PD, and MS [5,6], which causes derangement of the mitochondrial machinery [7], suppressing glutamate reuptake by astrocytes and resulting in glutamate-mediated toxicity [8].Immune cell infiltration from the periphery into the CNS as well as the inflammatory responses mediated by reactive astrocytes and by activated microglia in the CNS also have been implicated in AD, PD, and MS [9].Elucidating the shared and distinct molecular crosstalk could provide a biological basis for the treatments of AD, PD, and MS, which remains controversial in the field.
The establishment of single-nucleus RNA-sequencing (snRNA-seq) databases over the last decade permits global bioinformatics analyses of gene expression in different cells.Transcriptomic profiling, through snRNAseq of patient-derived tissues, can address confounding by cellular composition, providing previously unattainable insight into cell-type-specific transcriptomic pathology [10,11].Existing data resources have not yet been fully exploited to understand the causal disease pathways in AD, PD, and MS.With this in mind, we integrated several existing human snRNA-seq datasets to gain a comprehensive view of the cell-type-specific transcriptional changes of these diseases in the CNS.Our results showed the shared and distinct transcriptional changes in multiple cell types among AD, PD, and MS.We hypothesized that HSPB1 may be the core molecule of the shared pathological mechanism of the three diseases, which is the result of the induction by the blood-brain barrier (BBB) in response to cellular stress, providing insight into the nervous system diseases with unique pathogenic processes.Additionally, arctigenin has shown potential as a therapeutic drug for AD, PD, and MS.

Resources of single-nucleus RNA sequencing data
The selection of datasets was driven by the specific disease contexts and the relevant anatomical sites, which mainly based on the specific disease pathology.The complete snRNA-seq data sets used in this study were downloaded from the NCBI Gene Expression Omnibus database (GEO; https:// www.ncbi.nlm.nih.gov/ geo/) under accession numbers GSE138852 [12], GSE174367 [13], GSE157783 [14], and GSE118257 [15] (Table 1; Additional file 2: Table S1).Dataset GSE138852 contains entorhinal cortex samples from control and AD brains (n = 6 per group).GSE174367 contains prefrontal cortex samples from control (n = 8) and late-stage AD brains (n = 12).GSE157783 contains midbrain samples from idiopathic PD (n = 5) and control brains (n = 6).GSE118257 contains white matter samples from patients with MS (n = 12) and controls (n = 9).

Quality control, normalization, batch correction, dimensionality reduction, and clustering
All bioinformatic analyses were performed with Seurat (version 4.0; https:// satij alab.org/ seurat/) [16] in R software (version 4.0.2) for data processing and analysis.We obtained the gene expression matrix and metadata information of GSE138852, GSE174367, and GSE118257 datasets from GEO database.The metadata information includes cell-specific details such as cell origin, tissue source, disease status, and other relevant annotations.The expression matrix was directly transformed into a Seurat object using the CreateSeuratObject() function.Subsequently, we systematically incorporated the metadata information into this Seurat object.For the GSE157783 dataset, we downloaded three files from GEO, namely the cell barcode file, gene feature file, and expression matrix file.To integrate these files, we employed the Read10X() function, which resulted in an expression matrix with genes represented as rows and cells as columns.Then, the CreateSeuratObject() function was employed to generate a Seurat object.Finally, we utilized the merge function to integrate four Seurat objects together, facilitating further analysis.Cells with less than 200 genes, mitochondrial counts greater than 5%, and genes expressed in less than three cells were filtered out for quality control.The data was first normalized by functions Normalize-Data and ScaleData functions.Then FindVariable function was applied to select the top 2000 variable genes.We conducted principal component analysis (PCA) using the top 2000 variable genes.We selected the top 10 principal components (PCs) for sub-clustering and set the resolution parameter to 0.1.To mitigate batch effects and nonbiological technical biases, we employed the Harmony package in Seurat (https:// github.com/ immun ogeno mics/ harmo ny).Default parameter values were utilized for Harmony settings.The expression of known marker genes was used as a reference for annotation of different cell types.The results were visualized by Uniform Manifold Approximation and Projection (UMAP).

Differential expression analysis and functional enrichment analysis
To identify the genes that were upregulated in AD, PD, and MS brains across different cell types, the FindMarker function (Logfc.threshold= 1, p < 0.05, Only.pos = T, pct.1 − pct.2 > 0.2, Wilcoxon test) in Seurat was used.In addition, we defined the active (MS_Active) and chronic active (MS_CA) samples as the disease group to identify the truly upregulated genes in MS.Functional enrichment analysis was carried out using the Gene Set Enrichment Analysis (GSEA) [17] and hypergeometric tests with the clusterProfiler R package.Gene ontology (GO) enrichment analyses was performed using clusterProfiler R package with the input species set to Homo sapiens.The protein-protein interaction (PPI) analysis of differentially expressed genes (DEGs) was analyzed using STRING database (http:// string-db.org).

Gene regulatory network analysis
Single-cell regulatory network inference and clustering (SCENIC) analysis was performed using utilizing pySCE-NIC (v0.10.0) based on the hg19-tss-centered-10 kb-10species databases (https:// github.com/ aerts lab/ pySCE NIC) [18].Default parameters were used for the SCENIC workflow, and the raw count matrix from all the samples was used as the input.A three-step process was used for the analysis.Firstly, we calculated the co-expression modules and evaluated the weight between transcriptional factors (TFs) and their target genes using GRN-Boost, and then TFs with direct targets (regulons) were identified using RcisTarget.Finally, the activity of each regulon in each cell was evaluated using AUCell.For visualization, the average regulon activity (AUC) scores for each cell type were calculated, and a rank plot of regulons was drawn using ggplot2.

Connectivity map bioinformatics analysis
The connectivity map was analyzed on Connectivity Map (https:// clue.io) [19,20] with default settings.The differentially up-regulated genes in AD, PD, and MS were uploaded to the Connectivity Map database for analysis
The UMAP plot displayed the same expression patterns of the markers (Additional file 1: Fig. S1).We estimated cell-type proportions across disease groups and observed a decrease in oligodendrocytes and a marked increase in microglia occurred in the active (MS_Active) and chronic active lesion areas (MS_CA) of MS patients, and the midbrain of PD relative to the control groups; these changes were not observed in the AD group (Fig. 2D).OPCs were also significantly reduced in the active and CA lesion areas of MS, which is consistent with demyelination in the MS brain (Fig. 2D).

Defining potential shared gene regulatory network among AD, PD, and MS
To identify the shared gene regulatory network among AD, PD, and MS, we analyzed the association degree of DEGs (AD vs. control, PD vs. control, MS_Active and MS_CA vs. Ctrl_NWM) connectivity among the three diseases (Fig. 3A, Additional file 3: Table S2, Additional file 4: Table S3, Additional file 5: Table S4).Most genes upregulated in AD were expressed in excitatory neurons, astrocytes, and pericytes (Fig. 3A, Additional file 3: Table S2).PD-associated upregulated genes were found in almost all endothelial cells and astrocytes (Fig. 3A, Additional file 4: Table S3).Most of the upregulated genes in MS were found in endothelial cells, microglia, OPCs, and astrocytes (Fig. 3A, Additional file 5: Table S4).Upregulated DEGs that shared at least two cell types among three diseases were defined as hub-genes (Fig. 3A).Interestingly, the hub gene HSPB1 was identified to be significantly upregulated in the astrocytes, pericytes, OPCs, and excitatory cells of AD brains, and in endothelial cells of PD and MS brains (Fig. 3A, Additional file 3: Table S2, Additional file 4: Table S3, Additional file 5: Table S4).HSPA1A was found upregulated in oligodendrocytes, OPCs, pericytes, astrocytes, and excitatory cells of AD brains, and of oligodendrocytes, pericytes, and endothelial cells of PD brains (Fig. 3A, Additional file 3: Table S2, Additional file 4: Table S3).DNAJA1 was upregulated in pericytes of AD and in endothelial cells of PD (Fig. 3A, Additional file 3: Table S2, Additional file 4: Table S3).Additionally, SERPINH1, AEBP1, and DNAJB1 were similarly upregulated in both AD and MS (Fig. 3A, Additional file 3: Table S2, Additional file 5: Table S4).
To summarize the integrated ranking in the pathway, we performed GO analysis using GSEA (Fig. 3B).We prioritized the signaling pathways and identified five major pathways in each disease (p < 0.05, normalized enrichment score (NES) > 0, Fig. 3B).In the line part, the horizontal axis is the sorted gene, and the vertical axis is the corresponding Running Enrichment Score.There is a peak in the line chart, which is the Enrichment score of this gene set, and the corresponding gene is the core of the gene set.Collectively, two common pathways were identified in AD and PD: chaperone-mediated protein folding (GO: 0061077) and protein folding (GO: 0006457) (Fig. 3B), consistent with the well-established pathological mechanism of protein folding abnormalities in AD and PD [22].Interestingly, MS-related biological processes were mainly related to the transport of proteins to the endoplasmic reticulum (ER) (GO: 0006613, GO: 0072599, GO: 0045047, GO: 0006614) (Fig. 3B).
Both GO biological process and molecular functions analysis of the hub-genes revealed a significant enrichment in the protein folding and regulation of inclusion body assembly, which may module of the process of the pathological protein aggregation (Fig. 4A).HSPB1 and HSPA1A were both involved in multiple protein foldingrelated biological processes (Fig. 4B).STRING analysis of DEGs has unveiled the significant involvement of the heat shock protein (HSP) family in the pathogenesis of neurodegenerative diseases.In the brains of individuals with AD, HSPB1 and HSPA1A interacted with DNAJB1 in pericytes and excitatory cells, and with DNAJA1 in pericytes (Fig. 4C, Additional file 3: Table S2).In PD, HSPB1 and HSPA1A interacted with DNAJA1 in endothelial cells (Fig. 4D, Additional file 4: Table S3).Moreover, in the brains of individuals with MS, HSPB1 interacted with DNAJB1 in endothelial cells (Fig. 4E, Additional file 5: Table S4).We simultaneously observed that a series of ribosomal proteins play important roles in the pericytes of AD brains (Fig. 4C).For instance, RPL4, RPS7, RPL18, RPL18A, and RPS28 were all upregulated in the pericytes (Fig. 4C, Additional file 3: Table S2).These ribosomal proteins also interact directly or indirectly with HSPA1A (Fig. 4C).In PD, the HSP family interacted with IRF1 in endothelial cells, potentially mediating the regulation of inflammatory responses (Fig. 4D, Additional file 4: Table S3).In MS, HSPB1 not only demonstrated its role in upregulated proteins of astrocytes but also exhibited close interactions with upregulated ribosomal proteins of microglia (Fig. 4E, Additional file 5: Table S4).

Transcriptional dynamics and the regulators of AD, PD, and MS
The SCENIC was used to map the gene regulatory networks governing the different diseases and identify potential TFs modulating the DEGs in disease samples.We initially identified the modules of TFs which regulate the cell heterogeneity (Fig. 5A, Additional file 6: Table S5).Some modules of TFs were highly specific in AD, PD, and MS (M1, M4, and M7) (Fig. 5B, Additional file 6: Table S5).Using the regulon activity scores (RAS), we identified some regulons that were active in specific cell types among three diseases, regulating cell type-specific functions (Fig. 5B).Regulons activity in the M1 module was higher in the microglia of among three diseases (Fig. 5B).In the M4 module, the active regulons were mostly found in pericytes, especially in patients with PD and MS (Fig. 5B).Most of regulons in the M7 module were in the oligodendrocytes of PD and MS patients, indicating that PD and MS patients may share highly active TFs (Fig. 5B).Information regarding the modules and regulons is provided in Additional file 6: Table S5.
To obtain a specific correspondence between regulons and each cell type, we further analyzed the top 20 specific regulons in different cell types using the regulon specificity score (RSS) (Fig. 6A, Additional file 7: Table S6).We identified some shared specific regulons among three diseases by the Specificity score as candidate transcription factors underlying the gene expression differences in different cells (Fig. 6A).CEBPA and SPI1 in microglia, GRHPR in oligodendrocytes, and TFE3 and TP53 in astrocytes were shared by AD, PD, and MS (Fig. 6A).By integrating the transcription profiles previously identified upregulated in the three disease groups, we conducted an analysis to identify the candidate target genes regulated by these shared TFs (Fig. 6B).We identified several target genes that showed up-regulation in the disease groups, including PSAP and CNTN2 in AD, as well as CHI3L1 in PD (Fig. 6B, Additional file 8: Table S7).Interestingly, HSPBAP1, which interacts with hub-gene HSPB1 [23], was identified as a target gene regulated by both CEBPA and SPI1 (Fig. 6B).The corresponding motifs are showed in Fig. 6B.These analyses identified the upstream regulons that drive cell-type-specific state transitions toward disease.

Discovery of repurposable drugs
Connectivity Map (CMap) (https:// clue.io) [20] is a genome-wide transcriptional expression dataset of selected human cell lines processed by bioactive small molecules, including many drugs.CMap is used to discover functional links between drugs, genes, and diseases through transient signatures of common gene expression changes.With CMap, we identified that 30 shared drugs were inversely correlated with AD, PD, and MS (Fig. 7, Additional file 9: Table S8, Additional file 10: Table S9, Additional file 11: Table S10).Calcium channel blocker, topoisomerase inhibitor, MEK inhibitor, DNA methyltransferase inhibitor, and adenosine receptor agonist were inversely correlated with the up-regulated genes in AD, PD, and MS (Table 2).Arctigenin is a phenylpropanoid dibenzylbutyro lactone lignan compound which showed a potential therapeutic agent for AD, PD, and MS in our study (Fig. 7).Arctigenin has been shown having neuroprotective effects in vivo and in vitro of AD models [24].Previous studies showed that arctigenin effectively ameliorated memory impairment in AD mice model by targeting the production and clearance of β-amyloid [25].Arctigenin could improve the movement behaviors and upregulate dopamine and γ-aminobutyric acid levels in a PD mice model [26].Arctigenin may have anti-inflammatory and immunosuppressive properties via inhibiting Th17 cells in MS [27].These findings suggested that arctigenin holds promise as a potential therapeutic agent for AD, PD, and MS.

Discussion
AD, PD, and MS are CNS diseases that differ clinically.Previous studies have noted the importance of the common pathological changes in these diseases, such as small vessel disease [5] and inflammation [9].However, the molecular crosstalk between these three diseases remains largely unknown.This study aimed to identify the key molecules which regulate the pathogenic pathways in AD, PD, and MS.
We observed that the dataset obtained from the white matter of MS patients had a higher proportion of glial cells compared to neuronal populations.The nature of the tissue and the inherent heterogeneity of the cell type present pose challenges in achieving a completely unbiased representation.The result of the higher proportion of glial populations may lead to an overemphasis on glialspecific gene expression patterns while potentially masking or downplaying certain neuronal-specific signals.Furthermore, the altered proportions in the dataset could impact the functional interpretation of the results.It is important to consider the potential contributions of glial cells to the observed gene expression patterns and their implications for the underlying biological processes.
There was significant overlap in differentially expressed genes and pathways between AD and PD in our study.In the GSEA analysis, the genes of AD and PD were co-enriched into two identical pathways, chaperonemediated protein folding and protein folding, which has previously been shown to be associated with the abnormal protein folding in AD and PD.In our study, we found that HSPA1A, which encodes the major heat shock protein in the HSP70 family, was upregulated in certain cell types, including oligodendrocytes, OPCs, pericytes, and astrocytes in AD, as well as oligodendrocytes, pericytes, and endothelial cells in PD.Our findings align with prior research on AD and PD.Some studies showed that HSPA1A is upregulated in human entorhinal cortex samples and layer III pyramidal cells of AD, contributing to protein folding abnormalities, and altered synaptic transmission [28,29].A proteomics analysis of extracellular vesicles isolated from cerebrospinal fluid in AD also revealed elevated expression of HSPA1A compared to the mild cognitive impairment and control groups [30].Single-cell transcriptomics also uncovered the upregulation of HSPA1A (HSP70) in endothelial cells of patients with PD, which was further confirmed in PD patients' blood specimens and peripheral blood mononuclear cells [31,32].Furthermore, HSPA1A protein regulates the processing and generation of APP (amyloid precursor protein) and the production and aggregation of Aβ [33,34].HSPA1A exerted a significant restorative effect on neuronal morphology and functional status in the temporal cortex and hippocampal region in transgenic mouse [35].
The upregulation of HSPA1A can also play a preventive and decelerating role in PD-like neurodegeneration through its chaperone activity, which effectively inhibits α-synuclein aggregation and microglia activation [36,37].Taken together, our findings strongly support that HSPA1A could serve as a potential therapeutic target for AD and PD.
In our study, MS-related biological processes were mainly related to the transport of proteins to the endoplasmic reticulum (ER) (GO:0006613, GO:0072599, GO:0045047, GO:0006614), which is responsible for the proper folding and processing of polypeptide chains into functional proteins in cells.ER stress occurs when misfolded or unfolded proteins accumulate in the ER due to exogenous or endogenous factors.In MS, the destruction of the BBB leads to an escalation of pro-inflammatory damage, leading to cell damage.This is followed by oxidative damage and ER stress [38], leading to apoptosis or repair of CNS cells.These results are largely controlled by the unfolded protein response [39], which is consistent with our hypothesis.The preliminary findings indicate a strong association between ER function and MS, highlighting the need for additional investigation in this area.
The most prominent discovery arising from our analysis was the consistent upregulation of HSPB1 (HSP27) across among AD, PD, and MS.An increase of HSPB1 protein in the chronic active lesions of MS brains was reported [40], as we found in the present study.HSPB1 phosphorylation also was reported to be elevated in MS and AD [41,42].HSPB1 codes for a small heat shock protein that probably maintains folding competence in denatured proteins [43,44].In response to environmental stresses such as heat shock, HSPB1 protein is upregulated [45].The molecular chaperone activity of HSPB1 may regulate a broad spectrum of biological processes, including the phosphorylation of neurofilament proteins and their transport along axons [46].In AD, amyloid plaques and neurofibrillary tangles (composed of hyperphosphorylated tau) are present extracellularly and intracellularly, respectively [47].A transgenic mouse model of AD showed HSPB1 localization in plaques [48].By crossing the APPswe/PS1dE9 mouse model with the mouse model of HSPB1 overexpressed, spatial learning and electrophysiological parameters improved significant lye [49].HSPB1 can delay, but not prevent fibril formation by interacting with hyperphosphorylated tau [50]; this interaction is enhanced with increased tau phosphorylation [51].Moreover, overexpression of HSPB1 reduces the cellular deposition and cytotoxicity of α-synuclein in PD [52].HSPB1 binds to the surface of α-synuclein fibrils, thereby reducing their hydrophobicity [53].Based on the cumulative evidence from previous studies, we put forth the hypothesis that the observed upregulation of HSPB1 in our study strongly suggests its protective function in mitigating the formation of pathological protein aggregates, particularly in response to toxic stimuli or stressful conditions associated with neurological disorders.
We observed that the upregulation of HSPB1 occurred in different cells in the three diseases.HSPB1 exhibited upregulation in endothelial cells in PD and MS, as well as in astrocytes, pericytes, OPCs, and excitatory cells in AD.Initially, we attributed it to a non-specific cellular response.However, we have observed that the Threshold normalized connectivity score (normalized_cs) was set at < 0, and p < 0.05 upregulation of HSPB1 is consistently present in the constituent cells of the BBB, including endothelial cells, pericytes, and astrocytes, across AD, PD, and MS.This finding indicates that HSPB1 plays a role in the disruption of the BBB in three diseases.A previous study suggested that the endothelium-targeted overexpression of HSPB1 ameliorated BBB disruption after ischemic brain injury [54,55].The immune response of HSPB1 was also found to be present in the temporal cortex of patients with epilepsy and was mainly confined to vascular walls and glial cells [56].As a selective physical barrier, the BBB plays a protective role in maintaining the environmental balance in the brain.One hypothesis we consider is that conceptual BBB damage may play an important role in the pathogenesis of AD, PD, and MS.BBB dysfunction contributes to the onset and progression of AD [57,58], PD [59], and MS [60] as an upstream or downstream events.However, clinical studies in this area have shown that the BBB is less pronounced in patients with AD, PD, and MS.The BBB damage is more associated with cerebrovascular disease [61].The recent interest in AD immunotherapy has provided a reassessment of the BBB damage, as intact BBB is a potential barrier to effective treatment of anti-amyloid immunoglobulin [62].At the same time, dysfunctional BBB may be a risk factor for immune-mediated toxicity, including autoimmune encephalitis.Another reason for revisiting the role of BBB is that the current understandings of the pathogenesis of AD, PD, and MS are incomplete.Our findings suggest that HSPB1 is initially induced in regions where BBB is disrupted in response to cellular stress, particularly in endothelial cells and pericytes associated with blood vessels, as well as astrocytes.This endogenous response may serve as a protective mechanism against vascular and BBB injury.Further investigation is required to understand the mechanism by which HSPB1 maintains BBB structure and function.PPI analysis of DEGs revealed that the upregulated members of the HSP family may engage in interactions with specific proteins within various cells implicated in neurodegenerative diseases.In our study, the pericytes of AD brains appeared to be a crucial focal point.HSP family genes were upregulated in pericytes and interacted with the upregulated ribosomal protein family in AD brains.Different ribosomal proteins within the ribosome bear distinct roles and functions, collaborating to ensure accurate and efficient protein synthesis.Interestingly, previous study has also reported the elevated expression of ribosomal proteins in the capillaries of AD brains, which is remarkably like what we observed in pericytes [63].The researchers proposed that ribosomal function is augmented in the cerebral blood vessels of AD patients, leading to an aberrant protein translation network.The clinically used anti-AD drugs donepezil and tacrine were reported to inhibit ribosome biosynthesis [64].These viewpoints suggested the significant involvement of ribosome biosynthesis in the pathogenesis of AD and the importance of it as a therapeutic target.Increased ribosome biosynthesis also occurred in the microglia of MS, which indirectly interacted with HSPB1.This suggests the potential involvement of microglial protein synthesis in the pathogenesis of MS.In the endothelial cells of PD, we observed a concurrent upregulation and interaction between the HSP family and the inflammation-related gene IRF1.IRF1 is involved in immune responses, cell apoptosis, and tumorigenesis.Previous study has found that IRF1 protein is upregulated in α-Syn overexpressed SH-SY5Y cells and the substantia nigra of A53T α-Syn transgenic mice.Furthermore, α-Syn overexpression facilitates the translocation of IRF1 from the cytoplasm to the nucleus, mediating neuroimmune responses [65].Our research demonstrated the upregulation of IRF1 and HSP family genes in endothelial cells, revealing that besides immune cells, protein processing and immuneinflammatory responses in endothelial cells are also noteworthy.In this study, although we focused on the HSP family genes and their interacting genes, their prominent roles in distinct cell types appeared to vary across different diseases.This observation suggested that different cells may play diverse roles in the context of various diseases.
The transcriptional regulatory analysis identified several modules that are highly related to the pathogenesis of AD, PD, and MS, mainly in microglial, oligodendrocytes, and pericytes.We found that multiple top 20 specific TFs in microglia, oligodendrocytes, and astrocytes were shared across the AD, PD, and MS.Combined with the hub genes identified above, we found that these shared TFs can regulate multiple hub genes, including PSAP and CNTN2 in AD, and CHI3L1 in PD.The identification of these shared TFs has opened new avenues for studying their underlying pathological mechanisms.
We further identified that arctigenin could be a potential therapeutic drug for AD, PD, and MS.Arctigenin has also been shown having neuroprotective effect in vivo and in vitro of AD and PD models [24,26].Arctigenin could obviously attenuate the decrease of cell survival rates in SH-SY5Y cells with PD phenotypes by acting against cell apoptosis through the decrease of Bax/Bcl-2 and caspase-3, and by reducing the surplus reactive oxygen species production and downregulating the mitochondrial membrane potential [26].These results shed light on the potential therapeutic value of arctigenin as a treatment option for AD, PD, and MS.

Conclusions
In summary, we performed a large-scale snRNA-seq integration of diverse neurodegenerative diseases (AD, PD, and MS) from post-mortem patients.We showed that the shared and distinct molecular networks of the three diseases, which were significantly enriched for various well-known neurodegeneration-related pathobiological biological processes.HSPB1 was identified as a key molecule in the pathogenic mechanisms of these three diseases and is associated with the BBB.Arctigenin shows promise as a potential therapeutic agent for AD, PD, and MS.Further investigation into the pathogenic mechanisms associated with HSPB1 is required to gain a more comprehensive understanding of the involvement of the BBB in neurodegenerative diseases.

Fig. 1 Fig. 2
Fig.1The workflow of the study.After literature search, 64 cases from different single-nucleus sequence studies were selected.These studies were done on post-mortem human brain tissue, which were categorized into entorhinal cortex, prefrontal cortex, midbrain, and white matter, severely affected regions by AD, PD, and MS pathology, respectively

Fig. 3 Fig. 4
Fig. 3 Network of DEGs and enrichment analysis.A The network of upregulated genes across different cell types in AD, PD, and MS.B Enrichment plots from Gene set enrichment analysis (GSEA), the top five biological pathways sorted by normalized enrichment score across AD (left panel), PD (middle panel) and MS (right panel) are shown.The color of broken line represents different pathways

Fig. 5
Fig. 5 SCENIC analysis showing distinct and shared regulons across different cell types.A Heatmap of the identified TFs with hierarchical clustering.Different colors represent different modules.B Regulons activity scores of different cell types among modules (M1, M4, and M7)

Fig. 6
Fig. 6 Transcriptional regulatory network among AD, PD, and MS.A Specificity scores of regulons in microglia, oligodendrocytes, and astrocytes.The top 20 genes with higher activity are noted.B Network of TFs shared by the three diseases and their target genes (regulons), accompanied by the corresponding motifs of the TFs

Fig. 7
Fig.7 The shared 30 therapeutic drugs of AD, PD, and MS.Drug sensitivity analyses were performed using Connectivity Map (CMap).Threshold normalized connectivity score (normalized_cs) was set at < 0, and p < 0.05

Table 1
Summary of snRNA-seq datasets used in this study

Table 2
The names, Moa, and targets of the 30 drugs