- 1Department of Neurology, Liaocheng People’s Hospital and Liaocheng Hospital Affiliated to Shandong First Medical University, Liaocheng, China
- 2Clinical Laboratory, Liaocheng Veterans Hospital, Liaocheng, China
- 3Department of Neurosurgery, Liaocheng Traditional Chinese Medicine Hospital, Liaocheng, China
Background: The pathogenesis of Alzheimer’s disease (AD) is complex and multi-factorial. Increasing evidence has shown the important role of immune infiltration in AD. Thus the current study was designed to identify immune infiltration-related genes and to explore their diagnostic value in AD.
Methods: The expression data of AD patients were downloaded from the GEO database. The limma R package identified differentially expressed genes (DEGs) between AD and controls. The CIBERSORT algorithm identified differentially infiltrated immune cells (DIICs) between AD and controls. DIIC-correlated DEGs were obtained by Pearson correlation analysis. WGCNA was employed to identify DIIC-related modules. Next, LASSO, RFE, and RF machine learning methods were applied to screen robust DIIC-related gene signatures in AD, followed by the construction and validation of a diagnostic nomogram. Detection of the expression of related genes in the peripheral blood of Alzheimer’s disease and healthy volunteers by RT-PCR. In addition, the CTD database predicted chemicals targeting DIIC-related gene signatures in the treatment of AD.
Results: NK cells, M0 macrophages, activated myeloid dendritic cells, resting mast cells, CD8+ T cells, resting memory CD4+ T cells, gamma delta T cells, and M2 macrophages were differentially infiltrated between AD and controls. Pearson analysis identified a total of 277 DIIC-correlated DEGs between AD and controls. Thereafter, 177 DIIC-related genes were further obtained by WGCNA analysis. By LASSO, RFE and RF algorithms, CMTM2, DDIT4, LDHB, NDUFA1, NDUFB2, NDUFS5, RPL17, RPL21, RPL26 and NDUFAF2 were identified as robust gene signature in AD. The results of RT-PCR detection of peripheral blood samples from Alzheimer’s disease and healthy volunteers showed that the expression trend of ten genes screened was consistent with the detection results; among them, the expression levels of CMTM2, DDIT4, LDHB, NDUFS5, and RPL21 are significantly different among groups. Thus, a diagnostic nomogram based on a DIIC-related signature was constructed and validated. Moreover, candidate chemicals targeting those biomarkers in the treatment of AD, such as 4-hydroxy-2-nonenal, rosiglitazone, and resveratrol, were identified in the CTD database.
Conclusion: For the first time, we identified 10 immune infiltration-related biomarkers in AD, which may be helpful for the diagnosis of AD and provide guidance in the treatment of AD.
Introduction
Alzheimer’s disease (AD) is an irreversible neurodegenerative disease closely associated with aging (1). In the early stages, patients with the disease present with mild memory difficulties that progress to cognitive impairment and various cognitive dysfunctions that affect several areas of cognition, which seriously diminishes the quality of life (2). The prevalence of AD is increasing annually, and according to recent studies, the incidence of Alzheimer’s disease is expected to double and triple in Europe and worldwide by 2050, respectively (3). Despite many efforts to better understand the pathogenesis of AD, current clinical management and treatment are not very effective Early diagnosis and early intervention to slow down the progression of AD appear to be particularly crucial (4). Thus, further exploration of AD at a genetic level, identification of novel biomarkers, and refinement of the diagnostic model are of great significance to improving the intervention and treatment of AD patients.
In current clinical work, the diagnosis of AD is mainly based on the clinical manifestations of typical cognitive decline, with imaging manifestations such as atrophy of the hippocampus and other structures as an auxiliary diagnosis (5, 6). However, the typical clinical symptoms appear mostly in the middle and late stages of the disease process, which makes the early diagnosis of AD very challenging. Nowadays, with the rapid development of sequencing technology, disease sequencing data has increased dramatically (7). The bioinformatics parsing of sequencing data can identify biomarkers such as genes and targets that influence the development of transcendental development, which is of great value for the early diagnosis of AD (8, 9).
Recent studies have revealed that abnormal activation of neuroimmune cells and inflammatory response play a key role in the progression of AD (10). However, the diagnostic value of immune cell-related genes in AD remains unclear. A comprehensive understanding of the genetic and molecular mechanisms underlying immune cell involvement in AD could pave the way for developing targeted therapeutics and personalized treatment strategies. Furthermore, identifying reliable biomarkers and establishing diagnostic models may enable earlier and more accurate detection of AD, facilitating timely interventions and improving patient outcomes. The research team aims to discover novel genetic markers and molecular pathways associated with immune dysregulation in AD by analyzing comprehensive genomic data, transcriptomic profiles, and immune cell signatures. Additionally, the study intends to develop a robust diagnostic model capable of accurately differentiating AD patients from healthy individuals based on immune cell-related biomarkers.
Materials and methods
Data collection and processing
Three human AD microarray datasets, namely GSE85426, GSE63060, and GSE63061 were selected and downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). The details of the selected datasets are presented in Table 1. The selection criteria were as follows: (i) the datasets include both AD and control samples, (ii) The sample type used is blood, and (iii) The total number of samples included in the analysis is not less than 150. After eliminating the batch effects by the “sva” package in R (11), GSE85426, GSE63060 and GSE63061 were merged and used as AD training cohort. In addition, the E-MTAB-6094 dataset, including blood samples from 22 AD patients and 13 controls, was downloaded from the Array Express database (https://www.ebi.ac.uk/arrayexpress/) and used as an AD testing cohort. The “Limma” package (12) in R was used to identify DEGs between AD and control samples with FDR < 0.05 and |log2FC| > 0.5 as the cutoff threshold. The results were visualized in a volcanic plot and heatmap.
Evaluation of immune cell infiltration and their relationship with DEGs
The CIBERSORT algorithm was applied to calculate the proportion of 22 immune cells based on gene expression profiles of merging datasets. Student-t test was used to identify differentially infiltrated immune cells (DIICs) between AD and control samples in the training cohort with p < 0.05 as the cutoff threshold. Next, the Cor function of R was applied to perform correlation analysis between DIICs and DEGs to screen DIIC-correlated DEGs with p < 0.05 and |cor| > 0.3 as the cutoff threshold. The function of DIIC-correlated DEGs was analyzed by DAVID (13, 14), and GO terms and KEGG pathways with FDR < 0.05 were considered significantly enriched. In addition, the interactions among DIIC-correlated DEGs were evaluated in the STRING database (https://string-db.org/), followed by the construction of PPI network using DIIC-correlated DEGs with interaction score greater than 0.6 by Cytoscape software (15).
WGCNA analysis
WGCNA was carried out to build modules related to DIICs based on the gene expression profiles (16). After calculating the Pearson correlation coefficient, a similarity matrix was constructed. Then the similarity matrix was converted into an adjacency matrix, in which the optimal soft threshold was selected to build a scale-free network. Then adjacency matrix was converted to TOM for hierarchical clustering analysis of genes, followed by the identification of network modules by Dynamic TreeCut setting cut height as 0.995. Next, the correlations between modules and DIICs were calculated. To further narrow the scope of DIIC-correlated DEGs, Fisher’s exact test was performed by projecting DIIC-correlated DEGs into each module. DIIC-correlated DEGs in modules with p < 0.05 and FC >1 were selected for the following analysis.
Identification of robust DIIC-related biomarkers by machine learning
Three machine learning methods, including LASSO, RFE, and RF, were applied independently to screen DIIC-related biomarkers. Briefly, LASSO was performed by “lars” package in R (17) to screen the gene signatures under the optimal lambda with the slightest classification error. RF was performed by “randomForest” package in R (18), and gene features and their contributions were then analyzed by explaining the function of the “DALEX” R package (19, 20); SVM was performed by “caret” package in R at 100-fold cross-validation to determine the variables at the max accuracy (20). At last, the overlapped genes identified by LASSO, RF, and SVM were considered robust DIIC-related biomarkers in AD (21, 22). The expressions of DIIC-related biomarkers in AD and control samples from training and testing cohorts were compared and displayed in the bar charts.
Construction and validation of a diagnostic nomogram in AD
Then the overlapped genes were included in multivariate Cox regression analysis to construct the nomogram in training set by the “rms” package of R (23). The decision curves were plotted by the “rmda” package (24) of R to evaluate the clinical utility of the nomogram (25). In addition, the nomogram was tested in testing set. Furthermore, chemicals associated with AD and those biomarkers were searched in the CTD database (http://ctdbase.org/) and were employed to construct a chemical-gene network by Cytoscape software.
Quantitative real-time reverse transcriptions PCR
Total RNA was extracted from the peripheral blood of six patients with Alzheimer’s disease and six healthy adults from Liaocheng Hospital of Shandong First Medical University using Trizol reagent (AG21101, AG, CHINA). Inclusion criteria: people who met the diagnostic criteria for AD (NIA-AA criteria) published by the National Institute on Aging (NIA) and the Alzheimer’s Association (AA) in 2011 (Alzheimer’s disease group); people without cognitive impairment (control group).
Exclusion criteria: Patients with a clear history of stroke; Patients with psychiatric disorders; History of the immune system and severe hematologic disorders. Informed consent was obtained from the subjects or their legal guardians, and the Ethics Committee of Liaocheng Hospital, Shandong First Medical University, approved the study (NO.2023039).
RNA was reverse transcribed into cDNA with a commercial reverse transcription kit (AG11706, AG, CHINA), and quantitative real-time PCR was performed by the SYBR premixed Ex Taq kit (AG11718, AG, CHINA) and specific primers (Table 2). Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as the reference gene. Relative mRNA levels were expressed and calculated using the 2(-Delta CT) method. then subjected to one-way analysis of variance using Graphpad Prism8 software, and P value < 0.05 was considered a significant difference.
Results
Identification of DEGs in AD training cohort
Batch effects among GSE85426, GSE63060, and GSE63061 (Figure 1A) were first eliminated (Figure 1B) by “sva” package in R. A total of 404 DEGs, including 234 down-regulated and 170 up-regulated genes, were identified between AD (N = 374) and control (N = 328) samples in the merging AD training datasets (Table S1) and displayed in the volcanic plot (Figure 1C) and heatmap (Figure 1D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1147501/fimmu-14-1147501-HTML-r1/image_m/fimmu-14-1147501-g001.jpg)
Figure 1 Screening of significantly differentially expressed genes (A, B) The three data sets (GSE85426, GSE63060, GSE63061) remove the batch effect through the “sva” algorithm and merge them into one data set. The sample relationship before and after the batch effect removal. (C) The volcano map of DEGs. The green and red dots represent significantly lower and higher DEGs, respectively, the black horizontal line represents FDR<0.05, and the two vertical lines representlog2FC |>0.5. (D) Heatmap of DEGs. Black and white in the sample bar represent the AD and the control groups, respectively.
Identification of DEGs correlated with key infiltrated immune cells in AD
To identify key immune cells in AD, CIBERSORT was carried out to calculate the proportions of 22 immune cells in AD and control samples (Table S2). We found that compared to control samples, AD patients had significantly higher proportions of resting NK cells, M0 macrophages, activated myeloid dendritic cells, resting mast cells and lower proportions of CD8+ T cells, resting memory CD4+ T cells, gamma delta T cells, and M2 macrophages (Figure 2A), suggesting that those DIICs were key immune cells in AD. Next, 277 DEGs were detected to be significantly correlated with those DIICs (Table S3). Functional analysis showed that DIICs-correlated DEGs were significantly enriched into 26 GO terms (Figure 2B), such as SRP-dependent cotranslational protein targeting to membrane, nuclear-transcribed mRNA catabolic process, nonsense-mediated decay, translation, and 8 KEGG pathways (Figure 2C), including ribosome, oxidative phosphorylation, Huntington’s disease, Alzheimer’s disease, Parkinson’s disease, non-alcoholic fatty liver disease (NAFLD), proteasome and cardiac muscle contraction. Moreover, the PPI network, including 202 nodes and 1,096 interaction pairs, revealed close interactions among DIIC-correlated DEGs (Figure 2D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1147501/fimmu-14-1147501-HTML-r1/image_m/fimmu-14-1147501-g002.jpg)
Figure 2 Screening of DEGs significantly related to immune cells. (A) Distribution histogram of immune cells with significantly different distribution in different groups based on CIBERSORT evaluation. Immunologically significantly correlated DEGs, significantly correlated GO biological process, *P < 0.05, ***P < 0.001. (B) and KEGG signal pathway (C) bubble diagram, the horizontal axis represents the number of genes, the vertical axis represents the item name, the size of the dot represents the FoldEnrichment value, and the color represents significance. The larger the dot, the closer the color is and the higher the significance. (D) Interaction network diagram, color indicates the degree of significant difference.
Further screening of DIIC-correlated DEGs by WGCNA
A total of 277 DIIC-related DEGs were analyzed in WGCNA, the soft threshold power five was chosen to construct a scale-free network at R2 = 0.9 (Figure 3A). Then 7 modules were identified by Dynamic TreeCut (Figure 3B), and the multidimensional scaling analysis showed that genes in the same module were distributed in the same region (Figure 3C). Next, DIIC-correlated DEGs were projected to those modules. Fisher’s exact test revealed that DIIC-correlated DEGs were significantly enriched in blue, brown, green, and red modules (Table 3). The module-trait correlation analysis showed that the blue module was positively correlated with resting memory CD4 T cells (cor = 0.35) and delta gamma T cells (cor = 0.39), the brown module was negatively correlated with M0 macrophages (cor = -0.35), the green module was positively correlated with M0 macrophages (cor = 0.49) and negatively correlated with CD8 T cells (cor = -0.33), resting memory CD4 T cells (cor = -0.32), delta gamma T cells (cor = -0.34) and M2 macrophages (cor = -0.3), red module was negatively correlated with CD8 T cells (cor = -0.47) (Figure 3D). Thus, 118 DIIC-correlated DEGs in the blue module, 17 DIIC-correlated DEGs in brown module, 30 DIIC-correlated DEGs in green module, and 12 DIIC-correlated DEGs in the red module were used for the downstream analysis (Table S4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1147501/fimmu-14-1147501-HTML-r1/image_m/fimmu-14-1147501-g003.jpg)
Figure 3 WGCNA algorithm screening disease progression and immune-related modules. (A) Left figure: The weight parameter power selection diagram of the adjacency matrix. The horizontal axis represents the weight parameter power, and the vertical axis represents the square of the correlation coefficient between log (k) and log (p (k)) in the corresponding network. The red line represents the standard line where the square value of the correlation coefficient reaches 0.9. Right figure: Schematic diagram of average connectivity of genes under different power parameters. The red line indicates that under the weight parameter power of the adjacency matrix in the left figure, the average connectivity of network nodes is 1. (B) Gene dendrogram and module colors, each color represents a different module. (C) MDS diagram of DEGs contained in each module. (D) Correlation heat map between the proportion of sample immune cells and the partition of each module.
Construction and validation of a DIIC-related diagnostic model in AD
To screen robust DIIC-related biomarkers in AD, different machine-learning methods were employed to screen robust DIIC-related biomarkers in AD. Using RFE, 16 genes were identified as potential biomarkers for AD (Table S5, Figure 4A). By LASSO, 37 signatures were extracted as candidate biomarkers (Table S6, Figure 4B). 46 genes were identified by the RF algorithm (Table S5 and Figure 4C). At last, 10 overlapping genes, including CMTM2, DDIT4, LDHB, NDUFA1, NDUFB2, NDUFS5, RPL17, RPL21, RPL26 and NDUFAF2, were obtained and considered as robust biomarkers in AD. Except for DDIT4, the expressions of the other nine genes significantly decreased in AD compared with those in controls in the training set (Table S6, Figure 5A). Furthermore, the expression trend of these ten genes is similar in the test set, with significant differences in gene CMTM2, DDIT4, LDHB, NDUFS5, and RPL21 (Table S6, Figure 5B). Six Alzheimer patients and six healthy controls were included in the study (Table 4), and the expression trends of these genes verified by peripheral blood qPCR were consistent with the results of our previous bioinformatics analysis (Figure 6). Based on CMTM2, DDIT4, LDHB, NDUFA1, NDUFB2, NDUFS5, RPL17, RPL21, RPL26 and NDUFAF2, a nomogram was constructed (Figure 7A). The decision curves showed that the patients got the highest benefits from the combined nomogram model compared with the other single biomarker models (Figure 7B). The nomogram was also constructed in the testing cohort (Figure 7C), and the decision curves in the testing cohort further demonstrated the efficient clinical utility of the combined nomogram model (Figure 7D). Moreover, candidate chemicals targeting those biomarkers in AD were screened in the CTD database; a total of 28 chemicals, such as 4-hydroxy-2-nonenal, rosiglitazone, resveratrol, cannabidiol, colchicine were identified, followed by the construction of chemical-biomarker network composed of 94 chemical-biomarker pairs (Figure 8, Table S7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1147501/fimmu-14-1147501-HTML-r1/image_m/fimmu-14-1147501-g004.jpg)
Figure 4 Optimization algorithm for screening important immune-related markers. DEGs parameter diagram of RFE (A), LASSO algorithm (B) and RF (C) filtering features.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1147501/fimmu-14-1147501-HTML-r1/image_m/fimmu-14-1147501-g005.jpg)
Figure 5 The expression level distribution of characteristic genes in the combined training data set (A) and E-MTAB-6094 validation data set (B), *P < 0.05, ***P < 0.001.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1147501/fimmu-14-1147501-HTML-r1/image_m/fimmu-14-1147501-g006.jpg)
Figure 6 The relative expression levels of NDUFS5 (A), DDIT4 (B), LDHB(C), CMTM2 (D), RPL21 (E), NDUFB2 (F), NDUFA1 (G), RPL17 (H), RPL26 (I), and NDUFAF2 (J) in the model between control and AD samples identified by RT-PCR, GAPDH was used as a reference. * P < 0.05, ***P < 0.001, ns, non-significant.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1147501/fimmu-14-1147501-HTML-r1/image_m/fimmu-14-1147501-g007.jpg)
Figure 7 Construction and verification of Nomogram diagnostic model. (A) Nomogram model diagram based on the expression level of immune DEGs in the combined training data set based on ten features. (B) Model decision line diagram. (C) Nomogram model diagram of the expression level of immune DEGs in E-MTAB-6094 validation data set based on ten characteristics. (D) Model decision line diagram.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1147501/fimmu-14-1147501-HTML-r1/image_m/fimmu-14-1147501-g008.jpg)
Figure 8 Small molecule screening of chemical drugs related to characteristic immune DEGs. The linkage map of 10 characteristic genes and small chemical molecules. The square and yellow circles represent small molecules of characteristic genes and AD-related chemicals.
Discussion
Neuroimmune refers to the intricate interplay between the nervous and immune systems, encompassing the coordinated regulation of diverse cell types, including neurons, glial cells, and immune cells (26). Recent research has shed light on the pivotal role of aberrant activation of neuroimmune cells and heightened inflammatory responses in the progression of Alzheimer’s disease (AD) (10). This dysregulated activation of immune cells and the ensuing surge in inflammatory reactions lead to irreversible damage to neurons and the delicate functioning of the nervous system (10). However, the diagnostic value of immune cell-related genes in AD remains unclear. In the current study, a comprehensive analysis by CIBERSORT, WGCNA, and machine learning identified 10 immune cell-related genes in blood samples, followed by the construction of a reliable diagnostic nomogram and candidate chemicals in AD treatment. We hope to improve our understanding of the relationship between immune cells and AD at a genetic level and provide useful information for the treatment of AD patients.
The trends of expression levels of the ten hub genes obtained from this screening were identical in the training and validation sets, and to further validate the accuracy of our bioinformatics analysis, we verified their expression differences in peripheral blood in AD patients and healthy populations and obtained consistent results. Furthermore, evidence suggests that immune cell populations promise as diagnostic and prognostic biomarkers for AD (27). For example, a recent study demonstrated significant differences in the ratios of CD4+ T cells, NK cells, CD8+ T cells, and monocyte-macrophages between AD patients and healthy individuals at the single-cell level (28). These findings indicate that immune cell alterations and dysregulation may indicate disease status and progression.
CMTM2 is a member of a newly found gene family (chemokine-like factor-like MARVEL transmembrane domain-containing family members), the members which have been reported to be involved in various tumors, reproduction, and immunity (29). CMTM2 has been reported in spermiogenesis (30), hepatocellular carcinoma, HBV-related disorders (31), and gastric cancer (32). However, its role in immunity and the nervous system is unknown. In the current study, we found that the expression of CMTM2 significantly decreased in AD and was significantly correlated with M2 macrophages and gamma delta T cells (Table S3), which were DIICs between AD and controls. In AD, decreased M2 macrophage subset has been reported to indicate worse cognitive performance (33). Transplanting of M2 macrophages could reduce neuron loss, impair inflammation response and enhance cognitive ability (34). Gamma delta T cells produce IL-17, the accumulation of which was shown to be concomitated with the onset of cognitive decline in female AD mice model (35). Therefore, CMTM2 may regulate AD via M2 macrophages and gamma delta T cells, which needs further study.
The encoded product of DDIT4, RTP801 (also known as REDD1), is involved in Aβ-induced synaptic dysfunction by regulating Aβ’s cytotoxicity (36, 37), and knockdown of DDIT4 could improve cognitive ability and ameliorate neuroinflammation severity (38, 39). Those findings indicate that DDIT4 may be a promising therapeutic target in AD. In addition, our study found that the expression of DIIT4 was significantly correlated with resting NK cells. Consistent with our study, the proportion of resting NK cells was also significantly different between AD and control samples in other studies (40, 41), indicating an essential role of resting NK cells in AD. It has been reported that immune infiltration of NK cells in the brain may contribute to neuroinflammation in AD (42). However, the specific role of resting NK cells in the blood remains unclear. Given studies on DIIT4, it is inferred whether the role of DDIT4 in neuroinflammation is associated with resting NK cells.
In our study, the expression of LDHB was significantly decreased in AD samples and was correlated with resting mast cells and resting memory CD4+ T cells. LDHB is a vital player in lactate metabolism (43). Alteration of lactate metabolism was observed to have a close relationship with neuronal damage in AD mice model (44). In tumors, overexpression of LDHB could mitigate the effects of lactic acid on CD4 T cell-mediated cytokine production (45). However, to our knowledge, the role of LDHB on the CD4 T cells’ function in AD is unknown. Mast cells are involved in neuroprotection and neuroinflammation by secreting inflammatory mediators and cytokines (46), and activation of mast cells may accelerate the development of AD (47). Thus, LDHB may be involved in AD via lactate metabolism and neuroinflammation mediated by resting CD4+ T cells and resting memory cells.
RPL17, RPL21, and RPL26 are all ribosomal proteins. Among them, RPL26 was reported to be involved in methylation within AD neurons (48). However, the exact role of RPL26 in AD remains unclear. In the current study, we found that RPL17, RPL21, and RPL26 expressions correlated with M2 macrophages and gamma delta T cells. It will be, therefore, of interest in future studies to test the relationship between them and their role in modulating cognitive ability in AD.
NDUFA1, NDUFB2, NDUFS5, and NDUFAF2 are ubiquinone oxidoreductase subunits or complex assembly factors which play an important role in oxidative stress (49–51). Oxidative stress may be essential to AD development by promoting Aβ deposition, tau hyperphosphorylation, and the subsequent loss of synapses and neurons (52). However, the studies on those genes in AD are limited. Gly32Arg SNP is a mutational site in NDUFA1 Reported that may be associated with early-onset dementia (53). Interestingly, we found that those genes are also significantly correlated with M2 macrophages and gamma delta T cells, further indicating the critical role of M2 macrophages and gamma delta T cells in AD.
AD lacks effective treatment measures in clinical practice (54); Immunotherapy has garnered considerable attention as a promising treatment strategy aimed at modulating immune system function, reducing inflammation, and clearing pathological protein deposits to improve the disease condition and alleviate symptoms (55).
This study utilized gene expression profiling and bioinformatics methods to identify abnormal expression patterns of immune cell-related genes in individuals with Alzheimer’s disease (AD). Based on the these Hub genes, an effective diagnostic model of Nomogram was constructed and chemical drugs targeting these genes were predicted for the treatment of AD.These findings potentially reflect alterations within the immune system, including immune cell activation, cellular infiltration levels, and inflammatory responses. Investigating the mechanisms underlying immune cell activation, dysfunction, and interplay, as well as the regulatory networks of immune genes, holds promise for unraveling the pathogenesis of AD and establishing a theoretical foundation for developing novel therapeutic strategies. Nevertheless, certain limitations should be acknowledged in our current study. Firstly, although PCR validation confirmed the expression of immune cell-related biomarkers in the peripheral blood of AD patients and healthy volunteers, further comprehensive investigations are required to elucidate the specific mechanisms underlying the actions of these genes. Secondly, understanding the relationship between immune cells, biomarkers, and their interactions in regulating the development of AD necessitates further exploration.
Conclusion
In summary, our study identified CMTM2, DDIT4, LDHB, NDUFA1, NDUFB2, NDUFS5, RPL17, RPL21, RPL26, and NDUFAF2 as immune cell-related biomarkers and constructed a diagnostic model in AD. Our findings may provide novel ideas for managing and treating AD patients.
Data availability statement
The datasets in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, GSE85426, GSE63060 and GSE63061. https://www.ebi.ac.uk/arrayexpress/, E-MTAB-6094.
Ethics statement
The studies involving human participants were reviewed and approved by the Institutional Ethical Committee of the Liaocheng People’s Hospital and Liaocheng Hospital Affiliated to Shandong First Medical University (NO.2023039). The patients/participants provided their written informed consent to participate in this study.
Author contributions
XZ, MB, and GJ conceived and drafted the article, GZ and ZW provided bioinformatics analysis, GZ and HW provided samples and conducted RT-PCR test. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by Liaocheng Key Research and Development Plan(NO.2022YDSF26), Shandong Society of Geriatrics Scientific and Technological Research Project (NO. LKJGG2021W067).
Acknowledgments
The authors would like to thank GZ for his excellent technical support and XS for his critical review of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1147501/full#supplementary-material
References
1. Uddin MS, Tewari D, Mamun AA, Kabir MT, Niaz K, Wahed MII, et al. Circadian and sleep dysfunction in Alzheimer's disease. Ageing Res Rev (2020) 60:101046. doi: 10.1016/j.arr.2020.101046
2. Yan Y, Yang H, Xie Y, Ding Y, Kong D, Yu H. Research progress on Alzheimer's disease and resveratrol. Neurochem Res (2020) 45(5):989–1006. doi: 10.1007/s11064-020-03007-0
3. Scheltens P, De Strooper B, Kivipelto M, Holstege H, Chételat G, Teunissen CE, et al. Alzheimer's disease. Lancet (2021) 397(10284):1577–90. doi: 10.1016/S0140-6736(20)32205-4
4. Breijyeh Z, Karaman R. Comprehensive review on Alzheimer's disease: causes and treatment. Molecules (2020) 25(24):5789. doi: 10.3390/molecules25245789
5. Atri A. The Alzheimer's disease clinical spectrum: diagnosis and management. Med Clinics North America (2019) 103(2):263–93. doi: 10.1016/j.mcna.2018.10.009
6. Khan S, Barve KH, Kumar MS. Recent advancements in pathogenesis, diagnostics and treatment of Alzheimer's disease. Curr neuropharmacology (2020) 18(11):1106–25. doi: 10.2174/1570159x18666200528142429
7. Hrdlickova R, Toloue M, Tian B. Rna-seq methods for transcriptome analysis. Wiley Interdiscip Rev RNA (2017) 8(1):1364. doi: 10.1002/wrna.1364
8. Wang M, Song WM, Ming C, Wang Q, Zhou X, Xu P, et al. Guidelines for bioinformatics of single-cell sequencing data analysis in Alzheimer's disease: review, recommendation, implementation and application. Mol neurodegeneration (2022) 17(1):17. doi: 10.1186/s13024-022-00517-z
9. Graff-Radford J, Yong KXX, Apostolova LG, Bouwman FH, Carrillo M, Dickerson BC, et al. New insights into atypical Alzheimer's disease in the era of biomarkers. Lancet Neurol (2021) 20(3):222–34. doi: 10.1016/s1474-4422(20)30440-3
10. Jorfi M, Maaser-Hecker A, Tanzi RE. The neuroimmune axis of Alzheimer's disease. Genome Med (2023) 15(1):6. doi: 10.1186/s13073-023-01155-w
11. 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. Bioinf (Oxford England) (2012) 28(6):882–3. doi: 10.1093/bioinformatics/bts034
12. 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(7):e47. doi: 10.1093/nar/gkv007
13. Chen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling tumor infiltrating immune cells with cibersort. Methods Mol Biol (Clifton NJ) (2018) 1711:243–59. doi: 10.1007/978-1-4939-7493-1_12
14. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of Large gene lists using David bioinformatics resources. Nat Protoc (2009) 4(1):44–57. doi: 10.1038/nprot.2008.211
15. 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. doi: 10.1101/gr.1239303
16. Chen R, Mias GI, Li-Pook-Than J, Jiang L, Lam HY, Chen R, et al. Personal omics profiling reveals dynamic molecular and medical phenotypes. Cell (2012) 148(6):1293–307. doi: 10.1016/j.cell.2012.02.009
17. Goeman JJ. L1 penalized estimation in the cox proportional hazards model. Biom J (2010) 52(1):70–84. doi: 10.1002/bimj.200900028
18. Tolosi L, Lengauer T. Classification with correlated features: unreliability of feature ranking and solutions. Bioinformatics (2011) 27(14):1986–94. doi: 10.1093/bioinformatics/btr300
19. Sanz H, Valim C, Vegas E, Oller JM, Reverter F. Svm-rfe: selection and visualization of the most relevant features through non-linear kernels. BMC Bioinf (2018) 19(1):432. doi: 10.1186/s12859-018-2451-4
20. Deist TM, Dankers FJWM, Valdes G, Wijsman R, Hsu IC, Oberije C, et al. Machine learning algorithms for outcome prediction in (Chemo)Radiotherapy: an empirical comparison of classifiers. Med Phys (2018) 45(7):3449–59. doi: 10.1002/mp.12967
21. Zhao S, Zhang L, Ji W, Shi Y, Lai G, Chi H, et al. Machine learning-based characterization of cuprotosis-related biomarkers and immune infiltration in parkinson's disease. Front Genet (2022) 13:1010361. doi: 10.3389/fgene.2022.1010361
22. Lai G, Liu H, Deng J, Li K, Xie B. A novel 3-gene signature for identifying Covid-19 patients based on bioinformatics and machine learning. Genes (2022) 13(9):1602. doi: 10.3390/genes13091602
23. Wu J, Zhang H, Li L, Hu M, Chen L, Xu B, et al. A nomogram for predicting overall survival in patients with low-grade endometrial stromal sarcoma: a population-based analysis. Cancer Commun (London England) (2020) 40(7):301–12. doi: 10.1002/cac2.12067
24. Harrell FE Jr., Lee KL, Mark DB. Multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat Med (1996) 15(4):361–87. doi: 10.1002/(SICI)1097-0258(19960229)15:4<361::AID-SIM168>3.0.CO;2-4
25. Chi H, Peng G, Yang J, Zhang J, Song G, Xie X, et al. Machine learning to construct sphingolipid metabolism genes signature to characterize the immune landscape and prognosis of patients with uveal melanoma. Front Endocrinol (2022) 13:1056310. doi: 10.3389/fendo.2022.1056310
26. Hodo TW, de Aquino MTP, Shimamoto A, Shanker A. Critical neurotransmitters in the neuroimmune network. Front Immunol (2020) 11:1869. doi: 10.3389/fimmu.2020.01869
27. Heneka MT, Carson MJ, El Khoury J, Landreth GE, Brosseron F, Feinstein DL, et al. Neuroinflammation in Alzheimer's disease. Lancet Neurol (2015) 14(4):388–405. doi: 10.1016/s1474-4422(15)70016-5
28. Chen X, Holtzman DM. Emerging roles of innate and adaptive immunity in Alzheimer's disease. Immunity (2022) 55(12):2236–54. doi: 10.1016/j.immuni.2022.10.016
29. Guo X, Zhang S, Tan S, Bei C, Zhang H, Zhu X, et al. Downregulated Cmtm2 poses potential clinical significance in hepatocellular carcinoma. DNA Cell Biol (2020) 39(4):683–9. doi: 10.1089/dna.2019.5237
30. Sun CC, Li SJ, Hu W, Zhang J, Zhou Q, Liu C, et al. Comprehensive analysis of the expression and prognosis for E2fs in human breast cancer. Mol Ther (2019) 27(6):1153–65. doi: 10.1016/j.ymthe.2019.03.019
31. Chen S, Hu Q, Chen H, Zhang F-F, Duan L, Wang B, et al. Identification of serum Cmtm2 as a potential biomarker for hbv-related disorders. Dis Markers (2020) 2020:2032056. doi: 10.1155/2020/2032056
32. Qian C, Xu Z, Chen L, Wang Y, Yao J. Long noncoding rna Linc01391 restrained gastric cancer aerobic glycolysis and tumorigenesis Via targeting mir-12116/Cmtm2 axis. J Cancer (2020) 11(21):6264–76. doi: 10.7150/jca.48365
33. Hsieh S-W, Huang L-C, Chang Y-P, Hung C-H, Yang Y-H. M2b macrophage subset decrement as an indicator of cognitive function in Alzheimer's disease. Psychiatry Clin Neurosci (2020) 74(7):383–91. doi: 10.1111/pcn.13000
34. Zhu D, Yang N, Liu Y-Y, Zheng J, Ji C, Zuo P-P. M2 macrophage transplantation ameliorates cognitive dysfunction in amyloid-B-Treated rats through regulation of microglial polarization. J Alzheimers Dis (2016) 52(2):483–95. doi: 10.3233/JAD-151090
35. Brigas HC, Ribeiro M, Coelho JE, Gomes R, Gomez-Murcia V, Carvalho K, et al. Il-17 triggers the onset of cognitive and synaptic deficits in early stages of Alzheimer's disease. Cell Rep (2021) 36(9):109574. doi: 10.1016/j.celrep.2021.109574
36. Kim J-R, Lee S-R, Chung HJ, Kim S, Baek S-H, Kim JH, et al. Identification of amyloid beta-peptide responsive genes by cdna microarray technology: involvement of Rtp801 in amyloid beta-peptide toxicity. Exp Mol Med (2003) 35(5):403–11. doi: 10.1038/emm.2003.53
37. Morel M, Couturier J, Pontcharraud R, Gil R, Fauconneau B, Paccalin M, et al. Evidence of molecular links between pkr and mtor signalling pathways in abeta neurotoxicity: role of P53, Redd1 and Tsc2. Neurobiol Dis (2009) 36(1):151–61. doi: 10.1016/j.nbd.2009.07.004
38. Pérez-Sisqués L, Sancho-Balsells A, Solana-Balaguer J, Campoy-Campos G, Vives-Isern M, Soler-Palazón F, et al. Rtp801/Redd1 contributes to neuroinflammation severity and memory impairments in Alzheimer's disease. Cell Death Dis (2021) 12(6):616. doi: 10.1038/s41419-021-03899-y
39. Yi JH, Kwon H, Cho E, Jeon J, Lee J, Lee YC, et al. Redd1 is involved in amyloid B-induced synaptic dysfunction and memory impairment. Int J Mol Sci (2020) 21(24):9482. doi: 10.3390/ijms21249482
40. Tian Y, Lu Y, Cao Y, Dang C, Wang N, Tian K, et al. Identification of diagnostic signatures associated with immune infiltration in Alzheimer's disease by integrating bioinformatic analysis and machine-learning strategies. Front Aging Neurosci (2022) 14:919614. doi: 10.3389/fnagi.2022.919614
41. Liu Z, Li H, Pan S. Discovery and validation of key biomarkers based on immune infiltrates in Alzheimer's disease. Front Genet (2021) 12:658323. doi: 10.3389/fgene.2021.658323
42. Lu Y, Li K, Hu Y, Wang X. Expression of immune related genes and possible regulatory mechanisms in Alzheimer's disease. Front Immunol (2021) 12:768966. doi: 10.3389/fimmu.2021.768966
43. Urbańska K, Orzechowski A. Unappreciated role of ldha and ldhb to control apoptosis and autophagy in tumor cells. Int J Mol Sci (2019) 20(9):2085. doi: 10.3390/ijms20092085
44. Zhang M, Cheng X, Dang R, Zhang W, Zhang J, Yao Z. Lactate deficit in an Alzheimer disease mouse model: the relationship with neuronal damage. J Neuropathol Exp Neurol (2018) 77(12):1163–76. doi: 10.1093/jnen/nly102
45. Decking S-M, Bruss C, Babl N, Bittner S, Klobuch S, Thomas S, et al. Ldhb overexpression can partially overcome T cell inhibition by lactic acid. Int J Mol Sci (2022) 23(11):5970. doi: 10.3390/ijms23115970
46. Kempuraj D, Mentor S, Thangavel R, Ahmed ME, Selvakumar GP, Raikwar SP, et al. Mast cells in stress, pain, blood-brain barrier, neuroinflammation and Alzheimer's disease. Front Cell Neurosci (2019) 13:54. doi: 10.3389/fncel.2019.00054
47. Kempuraj D, Selvakumar GP, Thangavel R, Ahmed ME, Zaheer S, Raikwar SP, et al. Mast cell activation in brain injury, stress, and post-traumatic stress disorder and Alzheimer's disease pathogenesis. Front Neurosci (2017) 11:703. doi: 10.3389/fnins.2017.00703
48. Mastroeni D, Grover A, Delvaux E, Whiteside C, Coleman PD, Rogers J. Epigenetic changes in Alzheimer's disease: decrements in DNA methylation. Neurobiol Aging (2010) 31(12):2025–37. doi: 10.1016/j.neurobiolaging.2008.12.005
49. Mamelak AJ, Kowalski J, Murphy K, Yadava N, Zahurak M, Kouba DJ, et al. Downregulation of Ndufa1 and other oxidative phosphorylation-related genes is a consistent feature of basal cell carcinoma. Exp Dermatol (2005) 14(5):336–48. doi: 10.1111/j.0906-6705.2005.00278.x
50. Wang G, Tang C, Yan G, Feng B. Gene expression profiling of H9c2 cells subjected to H2o2-induced apoptosis with/without af-Hf001. Biol Pharm Bull (2016) 39(2):207–14. doi: 10.1248/bpb.b15-00601
51. Zhang G, Wang X, Li C, Li Q, An YA, Luo X, et al. Integrated stress response couples mitochondrial protein translation with oxidative stress control. Circulation (2021) 144(18):1500–15. doi: 10.1161/CIRCULATIONAHA.120.053125
52. Tönnies E, Trushina E. Oxidative stress, synaptic dysfunction, and Alzheimer's disease. J Alzheimers Dis (2017) 57(4):1105–21. doi: 10.3233/JAD-161088
53. Huttula S, Väyrynen H, Helisalmi S, Kytövuori L, Luukkainen L, Hiltunen M, et al. Ndufa1 P.Gly32arg variant in early-onset dementia. Neurobiol Aging (2022) 114:113–6. doi: 10.1016/j.neurobiolaging.2021.09.026
54. Cummings JL, Tong G, Ballard C. Treatment combinations for Alzheimer's disease: current and future pharmacotherapy options. J Alzheimer's Dis JAD (2019) 67(3):779–94. doi: 10.3233/jad-180766
Keywords: Alzheimer’s disease, immune infiltration, WGCNA, machine learning, diagnosis
Citation: Zhuang X, Zhang G, Bao M, Jiang G, Wang H, Li S, Wang Z and Sun X (2023) Development of a novel immune infiltration-related diagnostic model for Alzheimer’s disease using bioinformatic strategies. Front. Immunol. 14:1147501. doi: 10.3389/fimmu.2023.1147501
Received: 18 January 2023; Accepted: 03 July 2023;
Published: 20 July 2023.
Edited by:
Isaac Francos Quijorna, Spanish National Research Council (CSIC), SpainReviewed by:
Chao Cheng, Wuxi People’s Hospital of Nanjing Medical University, ChinaSusanta Pahari, Texas Biomedical Research Institute, United States
Copyright © 2023 Zhuang, Zhang, Bao, Jiang, Wang, Li, Wang and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zheng Wang, zhengyimingdao@gmail.com; Xiujuan Sun, salide111@163.com
†These authors have contributed equally to this work