Identification of immune cells infiltrating in hippocampus and key genes associated with Alzheimer’s disease
BMC Medical Genomics volume 16, Article number: 53 (2023)
Alzheimer’s disease (AD) is the most prevalent cause of dementia and is primarily associated with memory impairment and cognitive decline, but the etiology of AD has not been elucidated. In recent years, evidence has shown that immune cells play critical roles in AD pathology. In the current study, we collected the transcriptomic data of the hippocampus from gene expression omnibus database, and investigated the effect of immune cell infiltration in the hippocampus on AD, and analyzed the key genes that influence the pathogenesis of AD patients. The results revealed that the relative abundance of immune cells in the hippocampus of AD patients was altered. Of all given 28 kinds of immune cells, monocytes were the important immune cell associated with AD. We identified 4 key genes associated with both AD and monocytes, including KDELR1, SPTAN1, CDC16 and RBBP6, and they differentially expressed in 5XFAD mice and WT mice. The logistic regression and random forest models based on the 4 key genes could effectively distinguish AD from healthy samples. Our research provided a new perspective on immunotherapy for AD patients.
Alzheimer’s disease (AD) is an age-related neurodegenerative disorder that primarily involves memory decline and executive dysfunction. The main features of AD are abnormal aggregation of extracellular amyloid plaques and hyperphosphorylation of neuronal tau, which lead to synaptic loss and neuronal atrophy . Experts now believe that, like other common chronic diseases, AD is caused by a combinational factor , including age, environment, genetics, or specific susceptibility genes [3,4,5]. Cardiovascular disease, diabetes, obesity, and diet are generally considered to be factors that increase the risk of AD [6,7,8]. Activated microglia and astrocytes in AD patient brains usually have higher levels of inflammatory markers, which are generally distributed around amyloid plaques and neurofibrillary tangles [9, 10]. The role of immune response in the brain of AD patients may be bidirectional. On the one hand, pathogenic substances such as cell debris and protein aggregates can be eliminated by phagocytosis of microglia and astrocytes; On the other hand, persistent neuroinflammation is a chronic response of the innate immune system to neurological changes, and the sustained activation of glial cells causes harm to the nervous system . In cell culture studies, activated microglia could produce harmful substances, which may damage neurons [12,13,14]. Another feature of AD is the impairment of the blood brain barrier (BBB), and a compromised BBB might increase the permeability of immune cells and peripheral tissue molecules, which could lead to neurodegeneration . Both peripheral macrophages and neutrophils can infiltrate the brain of AD patients through the BBB and induce the activated innate immune response in AD patients [16,17,18]. In addition, activated T cells are also found in the brain of AD patients, where they could release inflammatory factors[19,20,21]. Amyloid β (Aβ), which aggregates alone, has been found to be a powerful complement activator . Activation of the complement system in AD patients results in the production of allergenic toxins that further promote inflammation , cytokine-induced APP production, and higher Aβ production due to increased APP amounts [24, 25]. Although numerous studies have shown that inflammation plays a vital role in the pathogenesis of AD, the identification of immune cells closely related to AD and the molecular mechanisms of AD pathogenesis requires further elucidation.
In this study, we assess the level of immune infiltration from the hippocampus based on the expression of given immune cell genes by single sample gene set enrichment analysis (ssGSEA), and revealed the differences in the immune infiltration of hippocampal tissue in AD and healthy samples. We identified key genes from highly correlated co-expression modules, which were closely associated with disease and immune cells. This study laid the foundation for further finding effective targets for curing AD and developing immunomodulatory regimens for effective treatment of AD.
Materials and methods
Data preprocessing and immune infiltration assessment
At first, we used “hippocampus” and “Alzheimer’s disease” as keywords to search the datasets in the GEO database, and we found the GSE5281 and GSE48350 datasets, which were both from the GPL570 platform. The ssGSEA could assess the infiltration of 28 immune cells for each AD and control sample through GSEA package . We retained the immune cells with significant differences as traits for subsequent analysis (p value < 0.05).
Weighted gene co-expression network analysis (WGCNA)
At first, we used the limma package to normalize the raw data of all samples, and then we removed genes containing NA. We used WGCNA package to construct a gene co-expression network to find key modules and module genes . Genes were clustered based on the phase dissimilarity machine. The division of modules was based on the high topological overlap of genes within the modules . We selected the modules associated with disease for subsequence analysis. For genes within modules, we further screened based on gene significance (GS) and module importance (MM). The genes with high MM and high GS were described as the central module genes, which were strongly associated with disease and candidate immune cell. In our study, the central module genes were the genes in the candidate module with |MM| > 0.8 and |GS| > 0.2.
Analysis of protein–protein interaction (PPI) network and identification of key genes
The Search Tool for the Retrieval of Interacting Genes (STRING) online tool  was used to analyze the PPI of central module genes with the default parameters. Then we used the cytoHubba plugin  of the Cytoscape (version 3.8.2) to identify the key genes . The cytoHubba provides 12 analysis algorithms to calculate hub genes in protein interaction network graphs, we used five of which to identify key genes in the PPI network, including Degree, Edge Penetration Component (EPC), Maximum Neighborhood Component (MNC), Density of Maximum Neighborhood Component (DMNC) and Maximum Group Centrality (MCC). We regarded the intersection of top10 genes, which were obtained by cytoHubba’s five algorithms as the key genes, the VennDiagram package was used to visualize these results .
Validation of the key genes
We constructed logistic regression model and random forest model by the intersection genes of cytoHubba’s five algorithms to explore the correlation between disease and key genes. We randomly divided all samples into test and training cohort according to the proportion of 3/7, we generated logistic regression model and random forest model in the training cohort and validated the performance of the models in the test cohort. The receiver operating characteristic (ROC) curves and confusion matrix were used to assess the validity of the models .
5.5-month-old heterozygous 5XFAD mice are housed in Tongji University Animal Center under standard conditions. Aβ42 began to accumulate in the brains of 5XFAD mice at 1.5 months of age . There are many Aβ plaques in the hippocampus at 5.5 months of age. In this study, 5.5-month-old 5XFAD mice were euthanized and hippocampal tissue was isolated for subsequent experiments. Both AD and control groups contained two female mice and two male mice.
RNA extraction and quantitative real-time PCR (qRT-PCR)
The total RNA of the hippocampus in all mice were extract by RNAiso Plus (9109, TaKaRa, China). According to the manufacturer’s instructions, qRT-PCR was performed by the AceQ Universal SYBR qPCR Master Mix (Q511-02, Vazyme, China). All genes’ expression levels were normalized to β-Actin by the comparative CT method (2−ΔΔCt). Table 1 showed the sequences of all RNA primers.
All data were visualized and analyzed by GraphPad Prism 8. T test was used to compare expression level between the AD and WT groups and p value < 0.05 were considered statistically significant.
Gene set enrichment analysis (GSEA)
Based on the median expression levels of key genes, we divided all samples into high and low expression groups, and GSEA was performed to explore hallmark pathways between the two groups . We used p value < 0.05 and p-adjust < 0.25 as the screening criterion for statistically significant.
Construction of mRNA-miRNA-lncRNA network
To further explore the miRNA and lncRNA regulatory networks associated with key genes, we constructed an mRNA-miRNA-lncRNA network based on the key genes screened in the previous result. At first, based on the “multiMiR” package , we used experimentally validated data to explore miRNAs associated with key genes. After obtaining miRNAs that interact with key genes, we used the starBase database to explore lncRNAs that interact with miRNAs . The lncRNAs that interact weakly with miRNAs are removed. Our screening criteria were that miRNA-lncRNA expression was negatively correlated in more than four cancers and validated by more than three clip-seq experiments. Finally, we used Cytoscape for visualization of the mRNA-miRNA-lncRNA network.
This study procedure was conducted methodically based on the steps outlined in the flow diagram (Fig. 1). Based on the search for keywords in Materials and Methods, we downloaded two datasets, GSE5281 and GSE48350, from the GEO database. As we mainly focused on the changes in transcriptome data of hippocampal tissue, we selected the data of hippocampal tissue, GSE5281 containing 23 samples and GSE48350 containing 62 samples. The R software was used to process the raw expression profiles of these two datasets, and the limma package was used to normalize the raw data . The batch effects of these two datasets were processed by the “sva” package .
Immune infiltration in the hippocampus of AD patients might be altered
As described in methods and materials, ssGSEA was performed on 29 AD samples and 56 control samples to assess the scores of 28 immune cells (Fig. 2A, B). Our results indicated that the scores of activated B cell, activated CD8 T cell, CD56 bright natural killer cell, effector memory CD8 T cell, eosinophil, immature B cell, macrophage, memory B cell, monocyte, myeloid derived suppressor cell, natural killer cell, natural killer T cell and type 17 T helper cell were significantly different between AD and healthy groups (p < 0.05), indicating that the level of immune cell infiltration might be altered in the hippocampus of AD patients.
Monocytes were the important immune cell associated with AD in the hippocampus
To identify disease-associated immune cell types associated with disease, we constructed gene co-expression modules using WGCNA. We first normalized the data from the datasets and subsequently removed genes containing NA. 2971 genes were eligible for further analysis. We built a scale-free (scale-free R2 > 0.85) co-expression network using soft threshold power β = 12 (Additional file 1: Fig. S1). These 2971 genes were clustered into 10 different color modules (Fig. 3A, B). Then, we analyzed the correlation between each module and immune cell types or sample types (AD and control) (Figure B). As a result, the green module was positively correlated with AD but negatively correlated with monocytes, and in contrast, the pink module was negatively correlated with AD but positively correlated with monocytes (p < 0.05). Additionally, monocytes exhibited highly correlation with both pink and green module. These results suggested that monocytes infiltrating the hippocampus might be the important immune cell associated with AD.
Pink and green central module genes were mainly involved in lipid metabolism, notch signaling pathway and material transport by Golgi-associated vesicles transport.
13 and 55 genes with high connectivity (|MM| > 0.8 and |GS| > 0.2) were screened from pink and green modules, respectively (Fig. 3C, D), which were considered as central module genes. To clarify the biological processes of the pink and green central module genes, we conducted GO, KEGG and Reactome enrichment analysis. According to our selection criteria, as for GO, the pink central module genes were mainly involved in misfolded protein reactions, proteasome-mediated proteolytic metabolic processes and low-density lipoprotein particle metabolism (Fig. 4A). The KEGG analysis suggested the pink central module genes were mainly involved in type 1 diabetes mellitus, legionellosis and endocrine and other factor-regulated calcium reabsorption (Fig. 4B). Reactome analysis demonstrated that pink central module genes were mainly involved in wnt signaling pathway and lipid metabolism (Fig. 4C). The same analysis was also performed on the green module genes. The GO analysis suggested that the 55 green central module genes were mainly enriched in histone modification and Golgi-associated vesicles transport (Fig. 4D). KEGG analysis revealed that the 55 green central module genes were mainly involved in thyroid hormone signaling pathway, notch signaling pathway, lysine degradation and C-type lectin receptor signaling pathway (Fig. 4E). As for Reactome analysis, green central module genes were mainly focused on notch signaling and the transport of substances between the Golgi and the endoplasmic reticulum (Fig. 4F). As a result, based on the frequency of terms, the pink central module genes were mainly affected lipid metabolism, and the green central module genes were mainly affected notch signaling pathway and material transport by Golgi-associated vesicles transport.
KEDLR1, SPTAN1, CDC16 and RBBP6 were identified as key genes associated with AD and monocytes
As for all the 68 central genes in the pink and green modules, we explored the PPI of these genes by STRING database, and the result was shown in Fig. 5A. The five algorithms of the cytoHubba, including EPC, MCC, MNC, DMNC and Degree, were used to process the PPI network to identify the top 10 genes (Table 2). KEDLR1, SPTAN1, CDC16 and RBBP6 were regarded as the key genes associated with monocytes and AD, which were the common genes identified by the five algorithms, respectively (Fig. 5B). Correlation analysis showed that KDELR1, SPTAN1 and RBBP6 were positively associated with AD and negatively associated with monocytes, while CDC16 was negatively associated with AD and positively associated with monocytes (Fig. 5C).
To validate correlation between KEDLR, SPTAN1, CDC16 and RBBP6 and AD occurrence, we constructed logistic regression model and random forest model. The area under curve (AUC) of logistic regression model was 0.789 (95% CI = 0.641–0.938), and the AUC of RF model was 0.828 (95% CI = 0.688–0.878) (Fig. 5D, F).
The results of confusion matrix were shown in Fig. 5E, G, and the accuracy and recall of the models were shown in Table 3. These results suggested the logistic regression model and random forest model based on KDELR1, SPTAN1, CDC16 and RBBP6 can distinguish AD patients from healthy samples. Then, we verified the expression values of these 4 genes between the two groups and found that they were significantly different in AD and healthy group (Fig. 6). In summary, multiple algorithms verified KDELR1, SPTAN1, CDC16 and RBBP6 were the key genes corelated with AD.
We also validated the relative mRNA levels of Kdelr1, Sptan1, Cdc16 and Rbbp6 in 5XFAD mice and WT mice. Compared with 5XFAD mice, the relatively mRNA levels of Kdelr1, Sptan1, Cdc16 and Rbbp6 were significantly increased in WT mice (Fig. 7).
GSEA revealed that lipid metabolism and immune response play important roles in AD
On the basis of the expression value of these 4 key genes, we performed GSEA to explore the potential pathways. We found that samples with high expression of KDELR1, SPTAN1, CDC16 and RBBP6 were enriched for adipogenesis, fatty acid metabolism, glycolysis, mTORc1 signaling, MYC targets V1 and protein secretion proteolysis (Fig. 8). In addition, samples with high expression of KDELR1 were enriched in four other gene sets including apical surface, hedgehog signaling, oxidative phosphorylation and UV response up (Fig. 8A), while samples with high expression of SPTAN1 were enriched in UV response down (Fig. 8B). Apical surface, Cholesterol homeostasis and UV response down were also significantly enriched in samples with high expression of CDC16 (Fig. 8C). Coagulation, interferon alpha response and interferon gamma response gene sets was significantly enriched in samples with high expression of RBBP6 (Fig. 8D). It has been shown that dysregulation of lipid metabolism is associated with aging, alterations in lipid rafts and brain lipid peroxidation levels. Our results showed the core role of KDELR1, SPTAN1, CDC16 and RBBP6 in the lipid metabolism and immune response.
Construction of key gene-related mRNA-miRNA-lncRNA network
For a further understanding of the role of key genes in AD occurrence, we built mRNA-miRNA-lncRNA network based the four key genes. Firstly, 12 miRNAs interacting with key genes were found in the multiMiR database based on four key genes KDELR1, SPTAN1, CDC16 and RBBP6, and then 38 lncRNAs interacting with 12 miRNAs were identified in the starBase database based on our screening criteria (Fig. 9). Thus, we obtained the mRNA-miRNA-lncRNA regulatory network of 4 key genes (containing 54 nodes and 60 edges). These interacting RNAs may be key mechanisms affecting the pathogenesis of AD.
Microglia are the brain-resident immune cells, and many studies regards Aβ-associated mononuclear phagocytes as microglia . There are now evidences that blood-derived monocytes can infiltrate the brain of AD patients through the BBB [45, 46]. In cell cultures incorporating Aβ42, the percentage of monocytes/macrophages (M/M) is significantly higher and M/M express chemokines to promote their migration through the BBB . Monocytes recruited in the brain can phagocytose Aβ in the brain parenchyma . In addition, not only Aβ in central nervous system can be removed, but also Aβ that spreads from the brain to the periphery can be captured and phagocytosed by peripheral monocytes. In this study, we analyzed transcriptomic data from hippocampus of AD patients in GEO database. We revealed the difference immune cell types in hippocampus between AD patients and healthy controls. In addition, we identified the pink and green modules are the key modules closely related to AD.
Based on the PPI network and cytoHubba, we identified 4 key genes associated with monocytes and AD, including KDELR1, SPTAN1, CDC16 and RBBP1, and found that these 4 genes differentially expressed in 5XFAD transgenic mice and WT mice. The GSEA and mRNA-miRNA-lncRNA network based on these 4 key genes further confirmed the possibility of these key genes affecting AD.
KDELR1, KDEL endoplasmic reticulum protein retention receptor 1. It could regulate integrated stress responses (ISR), and promote the naive T-cell survival in vivo , and regulates T-cell homeostasis through PP1 (protein phosphatase) . KDELR1 is also one of the candidate molecules associated with neurodevelopmental disorders , suggesting it may be one of the key molecules associated with the occurrence of AD. SPTAN1, spectrin alpha, non-erythrocytic 1, is essential for myelin formation . Patients with SPTAN1 mutations have also been found to present with peripheral neuropathy, severe dyslexia, and executive function difficulties . SPTAN1 is downregulated in the hippocampus of patients with medial temporal lobe epilepsy(MTLE), which is usually involved in drug-resistant seizures and cognitive deficits. Therefore, we believe that SPTAN1 is also a key potential molecule associated with Alzheimer’s disease. CDC16, cell division cycle 16, functions as a protein ubiquitin ligase. Together with other proteins, CDC16 forms a protein complex containing the Tre2-Bub2-Cdc16 (TBC) structural domain, the protein that belongs to the Rab-specific GTPase-activating protein (GAP) and is highly conserved in eukaryotes . The TBC and LysM Domain containing (TLDc) proteins containing the structural domain of TBC1 domain family member 24 (TBC1D24) are associated with neurodevelopmental disorders and are mainly involved in the oxidative stress response [56, 57]. Therefore, we speculate that CDC16 may also be one of the key molecules affecting neurodevelopment in AD. RBBP6, retinoblastoma binding protein 6. In various human cancers, RBBP6 is involved in the regulation of cell cycle and apoptosis . However, the role of RBBP6 has not been studied in AD, and it may be a new target related to AD pathology. What’s more, we performed the GSEA and mRNA-miRNA-lncRNA regulatory network to have a more comprehensive knowledge of the roles of key genes in AD.
To sum up, the current study initially assessed the abundance of immune cells in the hippocampus and identified monocytes were associated with AD. We identified and verified 4 key genes associated with the occurrence of AD by multiple methods, and also revealed the signaling pathways associated with immune response in AD, which might provide new insights for immunological studies in AD pathology.
Availability and data materials
All data generated or analyzed during this study are included in this published article. The datasets generated and/or analyzed during the current study are available in the GEO repository, https://www.ncbi.nlm.nih.gov/geo/, (GSE5281, GSE48350). We have deposited all the code used in this study and the input data for each figure on the Figshare website. All R code for the manuscript: https://doi.org/10.6084/m9.figshare.22014869.v2. All the input data for each figure: https://doi.org/10.6084/m9.figshare.22014806.v1.
Gene expression omnibus
Blood brain barrier
Single sample gene set enrichment analysis
Weighted gene co-expression network analysis
Kyoto encyclopedia of genes and genomes
Integrated stress responses
Search tool for the retrieval of interacting genes
Receiver operating characteristic
Gene set enrichment analysis
Area under curve
Integrated stress responses
Protein phosphatase 1
Medial temporal lobe epilepsy
TBC and LysM Domain containing
Lane CA, Hardy J, Schott JM. Alzheimer’s disease. Eur J Neurol. 2018;25(1):59–70.
Alzheimer’s A. Alzheimer’s disease facts and figures. Alzheimer’s Dement. 2016;12(4):459–509.
Silva MVF, Loures CdMG, Alves LCV, de Souza LC, Borges KBG, Carvalho MdG. Alzheimer’s disease: risk factors and potentially protective measures. J Biomed Sci. 2019;26(1):33.
Fluegge K. Environmental factors influencing the link between APOE ε4 and Alzheimer’s disease (AD) risk. Int Psychogeriatr. 2019;31(2):305–6.
Spinney L. Alzheimer’s disease: the forgetting gene. Nature. 2014;510(7503):26–8.
Stampfer MJ. Cardiovascular disease and Alzheimer’s disease: common links. J Intern Med. 2006;260(3):211–23.
Butterfield DA, Di Domenico F, Barone E. Elevated risk of type 2 diabetes for development of Alzheimer disease: a key role for oxidative stress in brain. Biochim Biophys Acta. 2014;1842(9):1693–706.
Akbaraly TN, Singh-Manoux A, Dugravot A, Brunner EJ, Kivimäki M, Sabia S. Association of midlife Diet with subsequent risk for Dementia. JAMA. 2019;321(10):957–68.
Navarro V, Sanchez-Mejias E, Jimenez S, Muñoz-Castro C, Sanchez-Varo R, Davila JC, et al. Microglia in Alzheimer’s disease: activated, dysfunctional or degenerative. Front Aging Neurosci. 2018;10:140.
Beach TG, Walker R, McGeer EG. Patterns of gliosis in Alzheimer’s disease and aging cerebrum. Glia. 1989;2(6):420–36.
Forloni G, Balducci C. Alzheimer’s Disease, oligomers, and inflammation. J Alzheimer’s Disease JAD. 2018;62(3):1261–76.
Lin L, Zheng LJ, Zhang LJ. Neuroinflammation, gut microbiome, and Alzheimer’s disease. Mol Neurobiol. 2018;55(11):8243–50.
Cao J, Hou J, Ping J, Cai D. Advances in developing novel therapeutic strategies for Alzheimer’s disease. Mol Neurodegener. 2018;13(1):64.
Ozben T, Ozben S. Neuro-inflammation and anti-inflammatory treatment options for Alzheimer’s disease. Clin Biochem. 2019;72:87–9.
Newcombe EA, Camats-Perna J, Silva ML, Valmas N, Huat TJ, Medeiros R. Inflammation: the link between comorbidities, genetics, and Alzheimer’s disease. J Neuroinflammation. 2018;15(1):276.
Hultman K, Strickland S, Norris EH. The APOE ɛ4/ɛ4 genotype potentiates vascular fibrin(ogen) deposition in amyloid-laden vessels in the brains of Alzheimer’s disease patients. J Cereb Blood Flow Metab. 2013;33(8):1251–8.
Fiala M, Liu QN, Sayre J, Pop V, Brahmandam V, Graves MC, et al. Cyclooxygenase-2-positive macrophages infiltrate the Alzheimer’s disease brain and damage the blood–brain barrier. Eur J Clin Invest. 2002;32(5):360–71.
Zenaro E, Pietronigro E, Bianca VD, Piacentino G, Marongiu L, Budui S, et al. Neutrophils promote Alzheimer’s disease—like pathology and cognitive decline via LFA-1 integrin. Nat Med. 2015;21(8):880–6.
Klegeris A, Walker DG, McGeer PL. Toxicity of human THP-1 monocytic cells towards neuron-like cells is reduced by non-steroidal anti-inflammatory drugs (NSAIDs). Neuropharmacology. 1999;38(7):1017–25.
Itagaki S, McGeer PL, Akiyama H. Presence of T-cytotoxic suppressor and leucocyte common antigen positive cells in Alzheimer’s disease brain tissue. Neurosci Lett. 1988;91(3):259–64.
Togo T, Akiyama H, Iseki E, Kondo H, Ikeda K, Kato M, et al. Occurrence of T cells in the brain of Alzheimer’s disease and other neurological diseases. J Neuroimmunol. 2002;124(1):83–92.
Lian H, Litvinchuk AA-O, Chiang AC, Aithmitti N, Jankowsky JL, Zheng H. Astrocyte-Microglia Cross talk through complement activation modulates amyloid Pathology in Mouse Models of Alzheimer’s Disease. J Neurosci. 2016;36(Electronic):1529–2401.
Poole S, Singhrao SK, Chukkapalli S, Rivera M, Velsko I, Kesavalu L, et al. Active invasion of porphyromonas gingivalis and infection-induced complement activation in ApoE−/− mice brains. J Alzheimers Dis. 2015;43:67–80.
Blasko I, Veerhuis R, Stampfer-Kountchev M, Saurwein-Teissl M, Eikelenboom P, Grubeck-Loebenstein B. Costimulatory effects of interferon-γ and Interleukin-1β or tumor necrosis factor α on the synthesis of Aβ1–40 and Aβ1–42 by human astrocytes. Neurobiol Dis. 2000;7(6):682–9.
Hu J, Akama KT, Krafft GA, Chromy BA, Van Eldik LJ. Amyloid-β peptide activates cultured astrocytes: morphological alterations, cytokine induction and nitric oxide release. Brain Res. 1998;785(2):195–206.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:559.
Li H, Sun Y, Chen R. Constructing and validating a diagnostic nomogram for multiple sclerosis via bioinformatic analysis. 3 Biotech. 2021;11(3):127.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Yu G, He Q-Y. ReactomePA: an R/Bioconductor package for reactome pathway analysis and visualization. Mol Biosyst. 2016;12(2):477–9.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43(Database issue):D447–D52.
Chin C-H, Chen S-H, Wu H-H, Ho C-W, Ko M-T, Lin C-Y. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8(Suppl 4):S11–S.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Chen H, Boutros PC. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics. 2011;12:35.
Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez J-C, et al. pROC: an open-source package for R and S + to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.
Oakley H, Cole SL, Logan S, Maus E, Shao P, Craft J, et al. Intraneuronal beta-amyloid aggregates, neurodegeneration, and neuron loss in transgenic mice with five familial Alzheimer’s disease mutations: potential factors in amyloid plaque formation. J Neurosci. 2006;26(40):10129–40.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005;102(43):15545-50.
Ru Y, Kechris KJ, Tabakoff B, Hoffman P, Radcliffe RA, Bowler R, et al. The multiMiR R package and database: integration of microRNA–target interactions along with their disease and drug associations. Nucleic Acids Res. 2014;42(17):e133–e.
Li J-H, Liu S, Zhou H, Qu L-H, Yang J-H. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein–RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(D1):D92–D7.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:1362–4962. (Electronic)).
Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics (Oxford England). 2012;28(6):882–3.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Kao Y-C, Ho P-C, Tu Y-K, Jou IM, Tsai K-J. Lipids and Alzheimer’s disease. Int J Mol Sci. 2020;21(4):1505.
Gold M, El Khoury J. β-amyloid, microglia, and the inflammasome in Alzheimer’s disease. Semin Immunopathol. 2015;37(6):607–11.
Malm TM, Koistinaho M, Pärepalo M, Vatanen T, Ooka A, Karlsson S, et al. Bone-marrow-derived cells contribute to the recruitment of microglial cells in response to β-amyloid deposition in APP/PS1 double transgenic Alzheimer mice. Neurobiol Dis. 2005;18(1):134–42.
Fiala M, Lin J, Ringman J, Kermani-Arab V, Tsao G, Patel A, et al. Ineffective phagocytosis of amyloid-β by macrophages of Alzheimer’s disease patients. J Alzheimers Dis. 2005;7:221–32.
Saresella M, Marventano I, Calabrese E, Piancone F, Rainone V, Gatti A, et al. A complex proinflammatory role for peripheral monocytes in Alzheimer’s disease. J Alzheimers Dis. 2014;38:403–13.
Fiala M, Cribbs DH, Rosenthal M, Bernard G. Phagocytosis of amyloid-β and inflammation: two faces of Innate Immunity in Alzheimer’s disease. J Alzheimers Dis. 2007;11:457–63.
Kamimura D, Arima Y, Tsuruoka M, Jiang J-j, Bando H, Meng J, et al. Strong TCR-mediated signals suppress integrated stress responses induced by KDELR1 deficiency in naive T cells. Int Immunol. 2016;28(3):117–26.
Kamimura D, Katsunuma K, Arima Y, Atsumi T, Jiang J-j, Bando H, et al. KDEL receptor 1 regulates T-cell homeostasis via PP1 that is a key phosphatase for ISR. Nat Commun. 2015;6(1):7474.
Pérez-Palma E, Saarentaus E, Ravoet M, De Ferrari GV, Nürnberg P, Isidor B, et al. Duplications at 19q13.33 in patients with neurodevelopmental disorders. Neurol Genet. 2018;4(1):e210–e.
Voas MG, Lyons DA, Naylor Stephen G, Arana N, Rasband Matthew N, Talbot WS. αII-Spectrin is essential for assembly of the nodes of ranvier in myelinated axons. Curr Biol. 2007;17(6):562–8.
Ylikallio E, Ritari N, Sainio M, Toppila J, Kivirikko S, Tyynismaa H, et al. De novo SPTAN1 mutation in axonal sensorimotor neuropathy and developmental disorder. Brain. 2020;143(12):e104–e.
Persike DS, Marques-Carneiro JE, Stein MLdL, Yacubian EMT, Centeno R, Canzian M, et al. Altered proteins in the hippocampus of patients with mesial temporal lobe epilepsy. Pharmaceuticals (Basel). 2018;11(4):95.
Wei Z, Zhang M, Li C, Huang W, Fan Y, Guo J, et al. Specific TBC domain-containing proteins control the ER-Golgi-plasma membrane trafficking of GPCRs. Cell Rep. 2019;28(2):554–66e4.
Balestrini S, Milh M, Castiglioni C, Lüthy K, Finelli MJ, Verstreken P, et al. TBC1D24 genotype-phenotype correlation: epilepsies and other neurologic features. Neurology. 2016;87(1):77–85.
Finelli MJ, Oliver PL. TLDc proteins: new players in the oxidative stress response and neurological disease. Mamm Genome. 2017;28(9–10):395–406.
Mbita Z, Hull R, Mbele M, Makhafola T, Dlamini Z. Expression analysis of RbBP6 in human cancers: a prospective biomarker. Anti-Cancer Drugs. 2019;30(8):767–73.
We appreciate the support of all lab members, and we would like to thank ‘GEO database’ for providing high quality biological data.
This work was supported by grants from the National Natural Science Foundation of China (31830111, 82171387), Key Research and Innovation Program of Shanghai Municipal Education Commission (2019-01-07-00-07-E00040).
Ethics approval and consent to participate
All experimental protocols were approved by the animal experimentation ethics committee of Tongji University (Approval-Number: TJAA08820101), all methods were carried out in accordance with relevant guidelines and regulations, and following the ARRIVE Guidelines Checklist as well.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Liu, C., Xu, S., Liu, Q. et al. Identification of immune cells infiltrating in hippocampus and key genes associated with Alzheimer’s disease. BMC Med Genomics 16, 53 (2023). https://doi.org/10.1186/s12920-023-01458-2
- Alzheimer’s disease
- Immune infiltration
- Key genes