- 1State Key Laboratory of Functions and Applications of Medicinal Plants, Guizhou Medical University, Guiyang, Guizhou, China
- 2Key Laboratory of Optimal Utilization of Natural Medicinal Resources, Guizhou Medical University, Guiyang, Guizhou, China
- 3Department of Anatomy, School of Basic Medical Sciences, Guizhou Medical University, Guiyang, China
- 4Guiyang Maternal and Child Health-Care Hospital, Guiyang, Guizhou, China
- 5Key Laboratory of Urological Diseases in Gansu Province, Gansu Nephro-Urological Clinical Center, Institute of Urology, Lanzhou University Second Hospital, Lanzhou, Gansu, China
Alzheimer's disease (AD) is a chronic neurodegenerative disease, and its underlying genes and treatments are unclear. Abnormalities in copper metabolism can prevent the clearance of β-amyloid peptides and promote the progression of AD pathogenesis. Therefore, the present study used a bioinformatics approach to perform an integrated analysis of the hub gene based on cuproptosis that can influence the diagnosis and treatment of AD. The gene expression profiles were obtained from the Gene Expression Omnibus database, including non-demented (ND) and AD samples. A total of 2,977 cuproptosis genes were retrieved from published articles. The seven hub genes associated with cuproptosis and AD were obtained from the differentially expressed genes and WGCNA in brain tissue from GSE33000. The GO analysis demonstrated that these genes were involved in phosphoribosyl pyrophosphate, lipid, and glucose metabolism. By stepwise regression and logistic regression analysis, we screened four of the seven cuproptosis genes to construct a diagnostic model for AD, which was validated by GES15222, GS48350, and GSE5281. In addition, immune cell infiltration of samples was investigated for correlation with these hub genes. We identified six drugs targeting these seven cuproptosis genes in DrugBank. Hence, these cuproptosis gene signatures may be an important prognostic indicator for AD and may offer new insights into treatment options.
Introduction
Alzheimer's disease is a neurodegenerative disease characterized by aging and irreversibility. It is predicted that dementia prevalence will double in Europe and triple globally by 2050 (1). The number of Americans with Alzheimer's dementia is estimated at 6.5 million over the age of 65 (2). The number of people with AD is forecast to grow to 13.8 million by 2060 without medical breakthroughs to prevent, slow, or treat it (2). According to the Alzheimer's Disease Foundation, Alzheimer's disease is the fifth leading cause of death among Americans aged 65 and older (3). It should be noted that about 5% of patients with AD inherit the disease in a chromosomal-dominant manner, and carriers can develop the disease as late as 40–60 years of age (3). Heritable familial AD (FAD) is most commonly caused by presenilin and the amyloid precursor protein (4). A major pathological feature of AD is the deposit of extracellular amyloid plaques and neurofibrillary tangles (NTFs) in the brain (5). According to the consensus, tau and amyloid β-protein (Aβ) are both essential to the neurodegenerative process in AD, and Aβ is located upstream of tau in the pathway (5). There is evidence indicating that amyloid plaque is not as detrimental as the soluble oligomeric form of Aβ, which causes synaptic dysfunction (6). However, the exact pathogenic form of Aβ and the molecular mechanism by which pathogenic Aβ causes subsequent synaptic and neurotoxic processes are yet unclear (7). Currently, the U.S. Food and Drug Administration (FDA) has licensed six medications (including donepezil, rivastigmine, galantamine, minocycline, and donepezil) to temporarily treat Alzheimer's symptoms without affecting the underlying brain changes or changing the course of the illness (3). Consequently, finding novel, potent targeted agents is urgently needed to further AD diagnosis and care.
In addition to Aβ and NTFs, neuroinflammation plays a significant role in the development of AD (1). It has been shown that Aβ causes microglia to release cytokines that pass through the blood-brain barrier with increasing age, leading to an accumulation of Aβ (8–10). Aβ-activated microglia secrete neurotoxic cytokines and chemokines (TNF, IL-6, IL-1, and CCL2), which cause neurological dysfunction and death (11). Clinical studies suggest anti-inflammatory cytokines (IL-2, IL-4, and IL-33) can modulate microglia activation and mitigate AD pathology (12, 13). The permeability of immune cells and molecules that pass through the blood-brain barrier increases with age, contributing to neurodegeneration in AD (14). Although a relationship between AD and neuroinflammation has been identified, it is unclear whether it is a cause or a consequence of the disease. Therefore, further research into the relationship between AD and inflammation could enhance our understanding of AD and contribute to developing more effective AD treatments.
A metal with redox activity, copper is involved in several metabolic processes in the brain (15). It serves as the active site for several copper enzymes, including cytochrome oxidase, ceruloplasmin, SOD1, and lysinase (16). The uptake and secretion of intracellular and extracellular copper ions are mediated by CTR1 and ATP7A/B (17). There is evidence that ATP7B (K832R) mutations may increase the risk of AD (18). A large number of copper ions accumulate in neurons in the brains of patients with Alzheimer's disease, primarily in amyloid plaques and tangles (19). It has been shown that copper ions can reduce copper ions in neurons and tissues of the brain when APP expression is increased (20, 21). The deficiency of intracellular copper ions can promote the production of Aβ, while the accumulation of extracellular copper ions can promote the aggregation of Aβ (22, 23). Copper ions can affect the division and proliferation of neutrophils in the peripheral circulation and have a positive correlation with the severity of AD (24). The neutrophil can cross the blood-brain barrier and aggregate near Aβ plaques in the AD mouse model (25). Depletion of neutrophils enhanced cognitive performance, reduced microglia, and reduced Aβ 1–42 levels in 3xTg-AD mouse brain homogenates (26). Thus, regulation of intracellular and extracellular copper iron transport is still under research for the treatment of AD. As one of the causes of neuronal ROS, copper ions combined with Tau protein can lead to the production of H2O2 in vitro. It was shown that intracellular delivery of copper ions reduced intracellular Tau phosphorylation in a mouse animal model of AD (APP/PS1) (27). Therefore, the relationship between copper ion-related metabolism, immune cells, and molecules and AD is not negligible.
Through the analysis of differentially expressed genes (DEGs) and Weighted gene co-expression network analysis (WGCNA) between ND and AD on the GSE33000 dataset, we identified seven hub genes related to cuproptosis and AD. The biological processes and pathways of seven hub genes have been analyzed using GO and gene set enrichment analysis (GSEA) in patients with AD. We developed a diagnostic model for AD based on stepwise regression and logistic regression analyses and validated it in the GSE15222, GSE48350, and GSE5281 datasets, respectively. It is evident from the AUC value of ROC that the model has good diagnostic performance and may be useful in the diagnosis of AD. Symptoms of neuroinflammation can be accompanied by tangles of tau, which can result from interactions with amyloid plaques. Thus, the significance of immune infiltrating cells for these hub genes was explored in AD samples. Finally, we retrieved drugs targeting seven hub genes from the DrugBank database, which have implications for the treatment of AD. The flowchart of the research is shown in Figure 1.
Materials and methods
Data acquisition
In total, 2,977 genes associated with cuproptosis were found in the published literature (28). In this study, all expression profiles were obtained from public databases. Gene expression data were obtained from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). GSE33000 expression data were obtained using the GPL10558 platform, containing 155 ND samples and 310 AD samples. GSE15222 expression data was used with the GPL2700 platform, containing 40 ND samples and 31 AD samples. GSE48350 expression data was used with the GPL570 platform, containing 27 ND samples and 21 AD samples. GSE5281 expression data was used with the GPL570 platform, containing 11 ND samples and 22 AD samples.
Differential expression analysis
Here, we used the “limma” package (version 3.40.6) for differential analysis to obtain differential genes between the AD and ND groups (29). Specifically, we obtained the expression profile dataset, performed multiple linear regression using the lmFit function, and further used the eBays function to compute moderated t-statistics, moderated F-statistics, and log-odds of differential expression by empirical Bayes moderation of the standard errors toward a common value, finally obtaining the differential significance of each gene. The differentially expressed genes (DEGs) between AD and AD were filtered with the threshold |logFC| > 1 and adj.p.val < 0.05.
Weighted gene co-expression network analysis
Weighted gene co-expression network analysis is a method for analyzing gene expression patterns of multiple samples (30). It can cluster genes by similar gene expression patterns, form modules, and analyze the relationships between modules and the clinical information of patients (30). First, we calculated the MAD (median absolute deviation) of each gene separately using the gene expression profile, and then we eliminated the top 50% of genes with the smallest MAD, removed the outlier genes and samples using the goodSamplesGenes method of the R package WGCNA, and further constructed a scale-free co-expression network using WGCNA. First, Pearson's correlation matrices and the average linkage method were both performed for all pair-wise genes. Then, a weighted adjacency matrix was constructed using a power function Amn = |Cmn|β (Cmn = Pearson's correlation between Genem and Genen; Amn = adjacency between Genem and Genen). β was a soft-thresholding parameter that could emphasize strong correlations between genes and penalize weak correlations. After choosing the power of 7, the adjacency was transformed into a topological overlap matrix (TOM), which could measure the network connectivity of a gene defined as the sum of its adjacency with all other genes for network generation, and the corresponding dissimilarity (1-TOM) was calculated. To classify genes with similar expression profiles into gene modules, average linkage hierarchical clustering was conducted according to the TOM-based dissimilarity measure with a minimum size (gene group) of 30 for the gene dendrogram. The sensitivity is 3. To further analyze the module, we calculated the dissimilarity of module eigengenes, chose a cut line for the module dendrogram, and merged some modules. In addition, we merged modules with distances < 0.25 and finally obtained 13 co-expression modules. Notably, the gray module was considered a set of genes that could not be assigned to any module. Finally, we calculated the correlation between module vectors and gene expression to obtain MM (MM >0.8), and 206 genes with high connectivity in the clinically significant module were identified as hub genes.
Identification of hub genes
To obtain genes related to cuproptosis genes and AD, DEGs, intersections of genes obtained by WGCNA and cuproptosis genes were taken using the “VeenDiagram” package in R software. The differential expression of the hub gene in ND and AD was represented using violin plots. The hypothesis tests used were the t-test and the Mann–Whitney U-test. The former was used if the data conformed to a normal distribution, and the latter if not. Significance was defined as p < 0.05.
Enrichment analysis
To investigate the biological mechanisms affecting the hub gene for AD, we performed a functional enrichment analysis. We first analyzed the biological process (BP) of gene ontology (GO) in which these genes are involved and presented the final results as a chord plot using the “GOplot” package in the R software. We obtained the GSEA software (version 4.3) from the GSEA (http://software.broadinstitute.org/gsea/index.jsp) website, divided the samples into high (≥50%) and low (< 50%) expression groups based on the expression levels of hub genes, and downloaded the c2.cp.kegg.v7.5.1.symbols.gmt subset from the Molecular Signatures Database (http://www.gsea-msigdb.org/gsea/downloads.jsp); the c2.cp.kegg.v7.5.1.symbols.gmt subset was downloaded to evaluate the relevant pathways and molecular mechanisms based on gene expression profiles and phenotypic groupings, setting a minimum gene set of 5 and a maximum gene set of 5,000, with 1,000 resamples. Screening conditions of p < 0.05 and FDR < 0.25 were considered to be statistically significant.
We used a query of multiple protein names (“IFI30,” “CLIC1,” “LYZ,” “PYGL,” “PLA1A,” “ALOX5AP,” and “A4GALT”) and organisms (“Homo sapiens”) to search the STRING website (https://string-db.org/). Following that, we set the following main parameters, namely, network type (“full STRING network”), network edge meaning (“evidence”), active interaction sources (“experiments”), the minimum required interaction score [“low confidence (0.150)”], and the maximum number of interactors to show (“no more than 50 interactors” in the first shell). Finally, the available experimentally determined binding proteins were obtained and visualized by Cytoscape software (version 3.9.1).
Signature for patients with AD
A logistic model is a statistical model that simulates the probability of an event by making the logarithm of the event a linear combination of one or more independent variables and is often used in disease diagnosis (31). In this study, logistic regression with two response variables was used, with 1 representing the AD sample and 0 representing the ND sample. Stepwise regression analysis was used to eliminate factors that were not significant for the response variable and only those that were significant were retained to simplify the model. Stepwise regression iteratively adds or removes variables from the model until the statistical value of the Akaike information criterion (AIC) is minimized. After that, logistic regression was used to fit the relationship between these significant factors and the response variables.
Finally, we performed ROC analysis using the R package pROC (version 1.18.0) to obtain the AUC. Specifically, we obtained the patient's risk score, performed ROC analysis using the roc function of pROC, and evaluated the AUC and confidence interval using the ci function of pROC to obtain the final AUC result of 0.91.
Immune infiltration
CIBERSORT, a deconvolution algorithm, was used to estimate 22 types of infiltration immune cells in patients with ND and AD by using the normalized gene expression (32). In this study, we analyzed the abundance of infiltrating immune cells in ND and AD tissues using CIBERSORT of the LM22 gene file, including plasma cells, naïve B cells, memory B cells, memory CD4+ T cells (activated and resting), naïve CD4+ T cells, CD8+ T cells, follicular helper T cells, γδ T cells, Tregs, macrophages (M0, M1, and M2), mast cells (activated and resting), dendritic cells (activated and resting), NK cells (activated and resting), neutrophils, monocytes, and eosinophils.
Drugs from the DrugBank
DrugBank (www.drugbank.ca) is a comprehensive, freely available web resource containing detailed drug, drug target, drug action, and drug interaction information on FDA-approved drugs and experimental drugs undergoing the FDA approval process (33). In this research, the “Targets” module of the DrugBank database was used to analyze the seven hub genes of drug targeting.
Statistical analysis
All statistical analyses in this study were performed using the R language. The t-test and the Mann–Whitney U-test were selected according to whether the data conformed to a normal distribution. Statistical significance was defined as p < 0.05.
Results
Identification of hub genes for AD and cuproptosis
To identify genes associated with AD, we analyzed a total of 299 differential genes between ND and AD in the GSE33000 dataset using the limma package with the screening criteria of adj.p.val < 0.05 and |logFC|> 1, including 139 upregulated genes and 160 downregulated genes (Supplementary Table S1). Expression differential genes were presented in a volcano plot (Figure 2A), and heat maps of the top 10 differentially expressed genes that were upregulated and downregulated, respectively, were plotted (Figure 2B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1064639/fneur-13-1064639-HTML/image_m/fneur-13-1064639-g002.jpg)
Figure 2. Differentially expression gene between ND and AD samples. (A) Volcano plot showing the significant genes found by limma analysis. Red genes represent significantly high expression in AD, blue genes indicate significantly high expression in ND, and gray genes mean no changes. (B) The top 20 genes significantly expressed in AD or ND samples, showing by the heatmap.
After removing outlier genes and samples using the goodSamplesGenes method in the WGCNA package, the expression profiles of 7,595 genes and 465 samples were taken from GSE33000 for constructing a weighted gene co-expression network (soft threshold force: 7, scale independence: 0.86, average linkage: 120.74, Figures 3A, B). The 13 different co-expression modules were obtained through dynamic tree cutting (module merging threshold: 0.25, minimum module: 30, Figure 3C). The correlation analysis was then performed for each module with clinical traits. It was shown that the midnight blue (r = 0.64, p = 8.6e-53) and dark red module (r = −0.62, p = 3.2e-55) had the highest positive and negative correlation with AD, respectively (Figure 3D). Thus, we selected the midnight blue module (containing 453 genes) with the highest correlation coefficients for further analysis. GS- and MM-related scatterplots showed that these genes were highly correlated with both modules and phenotypes (cor = 0.90, p < 0.001; Figure 3E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1064639/fneur-13-1064639-HTML/image_m/fneur-13-1064639-g003.jpg)
Figure 3. The WGCNA results. (A) Analysis of the scale-free fit index for various soft-thresholding powers (β). (B) The mean connectivity for the soft-thresholding powers. (C) Clustering dendrograms of genes, with dissimilarity based on the topological overlap, together with assigned module colors. (D) Correlations between different modules and clinical traits. (E) Correlation of module membership and gene significance in the midnight blue module.
By taking the intersection of 299 differential genes in GSE33000, 453 genes in the midnight blue module of WGCNA, and 2,977 genes related to cuproptosis genes, we finally obtained the seven hub genes of cuproptosis associated with AD (Figure 4A). The results of the violin plot showed that the seven hub genes were highly expressed in patients with AD (p < 0.05, Figure 5). In addition, the expression of these genes in GSE15222, GSE48350, and GSE5281 is shown in Supplementary Figure S1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1064639/fneur-13-1064639-HTML/image_m/fneur-13-1064639-g004.jpg)
Figure 4. Hub genes and GO analysis. (A) The intersections of the DE, cuproptosis genes, and midnight blue genes. (B) Biological processes in which the seven hub genes were involved.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1064639/fneur-13-1064639-HTML/image_m/fneur-13-1064639-g005.jpg)
Figure 5. The hub genes expression in the ND and AD samples of GSE33000. (A) IFI30, (B) CLIC1, (C) LYZ, (D) PYGL, (E) PLA1A, (F) ALOX5AP, and (G) A4GALT.
Enrichment for the hub gene
To further investigate the molecular mechanism of the seven hub genes in AD, we attempted to screen out the targeting of the hub gene binding protein. Based on the STRING tool, we found a total of binding proteins, which were validated by experimental data (Supplementary Figure S2). We performed an enrichment analysis to investigate the potential biological role of these hub genes. The GO analysis showed that four of the seven hub genes were involved in biological processes (BP), including phosphoribosyl pyrophosphate biosynthetic process, phosphoribosyl pyrophosphate metabolic process, lipoxin metabolic process, lipoxin biosynthetic process, cellular lipid metabolic process, and glycogen catabolic process (Figure 4B).
The results of the GSEA analysis showed that these hub genes were associated with neurodegenerative diseases (AD and Parkinson's disease), oxidative phosphorylation, leukocyte transendothelial migration, leishmania infection, hematopoietic cell lineage, complement and coagulation cascades, citrate cycle (TCA cycle), B cell receptor signaling pathway, systemic lupus erythematosus, cytokine receptor interaction, ECM receptor interaction, and Jak stat signaling pathway (Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1064639/fneur-13-1064639-HTML/image_m/fneur-13-1064639-g006.jpg)
Figure 6. GSEA revealed the enriched pathways of the hub genes. (A) A4GALT, (B) ALOX5AP, (C) CLIC1, (D) IFI30, (E) LYZ, (F) PLA1A, and (G) PYGL.
Construction and validation of diagnostic models
Logistic regression was used to perform a multi-gene prediction model based on GSE33000. Logistic regression was used to perform a multi-gene prediction model based on GSE33000. The prediction model was constructed for four of the seven cuproptosis genes, including IFI30, PLA1A, ALOX5AP, and A4GALT. The results showed that the constructed model had good diagnostic performance, and the area under the ROC curve was 0.91 (Figure 7A). We used the GSE48350, GSE15222, and GSE5281 datasets for validation, and the results showed that the area under the ROC curve was 0.75, 0.76, and 0.91, respectively (Figures 7B–D). The formula was as follows: [expression level of IFI30 × (−0.7449)] + [expression level of PLA1A × 0.6557] + [Expression level of ALOX5AP × 1.1483)] + [expression level of A4GALT × 1.4343] + 0.8891.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1064639/fneur-13-1064639-HTML/image_m/fneur-13-1064639-g007.jpg)
Figure 7. Receiver operating characteristic (ROC) curves and corresponding AUC values for the four expression cohorts. (A) GSE33000, (B) GSE48350, (C) GSE15222, and (D) GSE5281.
Immune infiltration in patients with AD
Immune cells, extracellular matrix, and other factors in the body are important for clinical therapeutic sensitivity and disease diagnosis. In this study, we used CIBERSORT to compare the proportion of 22 immune infiltrating cells in ND and AD samples (Figure 8A). The immune cell infiltration was compared between AD and ND samples in the boxplot (Figure 8B). The results showed that memory B cells, naïve CD4+ T cells, resting memory CD4+ T cells, NK cells resting, macrophages M2, activated mast cells, eosinophils, and neutrophils were significantly higher in patients with AD. In contrast, plasma cells, CD8+ T cells, activated memory CD4+ T cells, follicular helper T cells, regulatory T cells (Tregs), activated NK cells and resting mast cells were significantly lower in patients with AD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1064639/fneur-13-1064639-HTML/image_m/fneur-13-1064639-g008.jpg)
Figure 8. Immune cell infiltration between ND and AD samples. (A) The percentage of 22 immune cells in each sample. (B) Different levels of immune infiltrating cells between ND and AD samples.
Subsequently, the relationship between hub genes and immune infiltration was analyzed. The hub genes were significantly negatively associated with plasma cells, CD8+ T cells, activated memory CD4+ T cells, follicular helper T cells, Tregs, activated NK cells, and resting mast cells, while the opposite was true with resting memory CD4+ T cells, resting NK cells, monocytes, macrophages M1, macrophages M2, and neutrophils (Figure 9). The results suggest that the hub genes may play an important role in the immune microenvironment.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1064639/fneur-13-1064639-HTML/image_m/fneur-13-1064639-g009.jpg)
Figure 9. Correlation between the hub genes and immune infiltrating cells. **p < 0.01, ***p < 0.001.
Drug from DrugBank
In this study, we searched the DrugBank database for drugs targeting seven cuproptosis genes (Table 1). The FDA has approved six of these drugs; several others are being investigated intermittently; and one drug has been discontinued. There are no studies available for the A4GALT, PLA1A, and PYGL genes. It is believed that copper (DB09130) does have an effect on the CLIC1 gene and is mostly used in emergency contraception, non-traceable elements, and dietary supplements, although the exact mechanism is unknown due to the wide spectrum of enzymes that use copper ions as co-factors. Acting as a ligand for CLIC1 target genes, Artenimol (DB11638) also acts as a regulator of the cell cycle and inserts into membranes to form chloride channels at appropriate pH levels (34). In tissues and body fluids, sucrose (DB02772) acts as a nutritional supplement to Lysozyme C (encoded by the LYZ gene), which has an important role in enhancing immune responses (35, 36). As a feed additive, arsanilic acid (DB03006) is a toxic compound containing arsenic that induces blindness in animals. It has received approval for use in veterinary medicine for treating intestinal diseases in pigs and poultry (36). It is known that propyl alcohol (DB03175) is one of the best targets for lysozyme C and is primarily used for skin disinfection as well as a preservative that can be used in both clinical and domestic settings (35). Known as part of the non-essential amino acid family, aspartic acid (DB00128) can act on lysozyme C, but it needs to be studied further. Rose bengal (DB11182) is used as a diagnostic agent and is recommended for eye examinations, cornea, or conjunctiva (37) (Supplementary Table S2).
Discussion
In this study, seven cuproptosis genes were screened for association with AD using the public dataset GSE33000. GO analysis showed that they were involved in the phosphoribosyl pyrophosphate and lipoxin biological processes, as well as the cellular lipid metabolic process and glycogen catabolic process. There is no doubt that elevated levels of PRPP in humans are associated with an excess of uric acid production and accumulation, which is negatively correlated with Alzheimer's disease prevalence (38). Neuroinflammation may be beneficial for patients with early-onset Alzheimer's disease, and it may help reverse or at least slow tau protein accumulation in the brain, preventing dementia (25). By helping macrophage differentiation and activation, lipoxin, derived from arachidonic acid, could reduce the immune response in early patients with AD by releasing cytokines that decrease inflammation (39, 40). Based on these cuproptosis genes, prediction models were developed and validated using the GSE48350, GSE5281, and GSE15222 datasets, respectively. Furthermore, GSEA demonstrated that all of these pivotal genes were associated with neurodegenerative diseases, particularly AD, which validated these gene selections to some extent (Figure 6).
The diagnostic models developed from four of the seven cuproptosis genes (including IFI30, PLA1A, ALOX5AP, and A4GALT) can serve as a guide to the clinical diagnosis of AD. An antigen-presenting cell (APC)-expressing gamma-interferon-induced lysosomal thiol reductase (GILT, encoded by IFI30) reduces disulfide bonds through endocytic proteins, presenting immunogenic peptides bound to the major histocompatibility complex (MHC) class II (41). Consistent with our findings that the IFI30 gene is highly expressed in patients with AD, GILT is highly expressed in microglia surrounding Aβ and is involved in Aβ clearance (42, 43). Abnormalities of lipid metabolism in the brain are characteristic of AD (44, 45). It is possible to reduce the risk of AD by converting phosphatidylcholine to lysophosphatidylcholine-DHA in liver tissues and transporting it into the brain across the blood-brain barrier (46, 47). The 5-lipoxygenase-activating protein (FLAP, encoded by ALOX5AP) is widely distributed in the central nervous system and functions to regulate the activation of the 5-lipoxygenase enzyme (48). It has been shown that selective pharmacological inhibition of FLAP significantly reduces Aβ levels and deposition in amyloid precursor protein (APP) transgenic mice (49) and regulates endogenous tau metabolism in vivo (50). A4GALT is involved in regulating the synthesis of glycosphingolipids, which is associated with several neurodegenerative diseases, including AD and Parkinson's disease (36, 37). Neuroinflammation is more prevalent in patients with mild cognitive impairment and AD. There is research showing that a combination therapy that reduces amyloid plaque formation and limits neuroinflammation may be more effective than treating either alone (51).
Here, we further explored the level of infiltration of immune infiltrating cells in patients with AD. The results showed that the M0 macrophages were not significantly different in both AD and ND, while the M1/M2 macrophages were significantly lower in AD. According to the results, increased activation and differentiation of M0 macrophages into M2 macrophages reduces inflammation and contributes to delaying the onset of AD. Tregs, the subset of CD4+ T cells, are essential for maintaining immune homeostasis and downregulating patients with AD. The results have shown that early depletion of Tregs is associated with accelerated cognitive impairment and that restoring Tregs reduces Aβ deposition and improves cognition, but the results of this study remain controversial (52). Neutrophils are also upregulated in AD. A significant increase in neutrophils was found in AD brains and AD model mice, and hyperactivity of neutrophils is associated with AD (53). It was found that neutrophils may play a role in AD as they produce large quantities of reactive oxygen species, thereby causing AD and cognitive decline through the LFA-1 integrin (25).
The gamma interferon-inducible lysosomal thiol reductase (GILT, encoded by IFI30) is expressed in antigen-presenting cells (APCs), such as dendritic cells, monocytes/macrophages, and B cells, and it is highly expressed in microglia in the AD brain (42). The chloride intracellular channel 1 (CLIC1) protein, a potential marker of neurodegenerative processes, is significantly increased in peripheral blood mononuclear cells of patients with AD (54). The LYZ gene, a significant member of non-specific immunity, is upregulated in the cerebrospinal fluid (CSF) of patients with AD and inhibits the appearance of toxic Aβ oligomers (55). Abnormal brain glucose catabolic processes are associated with the formation of amyloid plaques in the brain and the onset of memory loss (56). In contrast, the relationship between PYGL (glucose metabolizing enzyme) and AD remains unclear.
Finally, the drugs targeting the above genes were retrieved from the DrugBank database. An excess of copper (DB09130) acts as a catalyst for a variety of biological processes (57), accumulating in neurofibrillary tangles and regulating APP gene expression (45, 46). Aβ peptide interaction with copper and other metals is thought to promote gain-of-function activity and lead to neurotoxicity (15). As an immunomodulator and neuroinflammatory drug, Artenimol (DB11638) inhibits neuronal apoptosis, modulates Tau autophagy, and protects AD mice from neuronal damage (58–60). There is no doubt that sucrose (DB02772) is a valuable nutritional supplement, but excessive consumption can lead to type 2 diabetes, and type 2 diabetes is associated with a high incidence of Alzheimer's disease (61). There is an increase in LYZ in the brains of transgenic mice and humans with AD, pointing to new therapeutic strategies to slow its progression (55). Propyl alcohol (DB03175) acts on the LYZ target, but its function still needs further study. In addition to its functionality in protein synthesis, aspartic acid (DB00128) is also involved in the urea cycle and gluconeogenesis and affects AD indirectly (62, 63). There is a xanthene dye, rose bengal (DB11182), which has been used to treat colon cancer (64). Furthermore, rose bengal inhibits the toxic effects caused by Aβ aggregation (65) and has a therapeutic effect on Tau aggregation (66).
Conclusion
In summary, we investigated seven cuproptosis genes and AD-related hub genes by using bioinformatics techniques. The biological role of these genes in AD development was explored. Further experiments are needed to confirm the function. Based on logistic regression analysis, we constructed a diagnostic model that can diagnose patients with AD by detecting the expression of several genes in the brain tissues. Additionally, immune cells expressed themselves more strongly in AD, indicating that they may be crucial to the immunological microenvironment. However, further research is needed to explore their specific effects. Currently, only a few drugs targeting these pivotal genes are predicted to alleviate AD, suggesting that additional drugs need to be developed.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
EZ and FD contributed to the conception, design of the study, and wrote the first draft of the manuscript. TC and SL organized the database. EZ, FD, and CX performed the statistical analysis. XS wrote sections of the manuscript. All authors contributed to the manuscript revision and read and approved the submitted version.
Funding
This study was supported by the Doctoral Scientific Research Foundation of Guizhou Medical University (No. YJ20066), the National Natural Science Foundation of China Cultivation Project of Guizhou Medical University (No. 20NSP051), the Scientific and Technological Project of Guizhou Province [No. ZK (2022) 398], and the Fund of Science and Technology Project of Guizhou provincial health commission (No. Gzwkj2023-564).
Acknowledgments
We thank all the patients who participated in the study. We also thank Dr. Jianming Zeng and his bioinformatics team for their reference codes.
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/fneur.2022.1064639/full#supplementary-material
Supplementary Figure 1. Expression of hub genes in the AD and ND samples GSE15222, GSE48350, and GSE5281.
Supplementary Figure 2. Protein-protein interaction of seven hub genes.
Supplementary Table 1. Differentially expressed genes between AD and ND samples in GSE33000.
Supplementary Table 2. Information on drugs targeting these seven hub genes.
References
1. Scheltens P, Strooper BD, Kivipelto M, Holstege H, Chételat G, Teunissen CE, et al. Alzheimer's disease. Lancet. (2021) 397:1577–90. doi: 10.1016/S0140-6736(20)32205-4
2. Rajan KB, Weuve J, Barnes LL, McAninch EA, Wilson RS, Evans DA. Population estimate of people with clinical Alzheimer's disease and mild cognitive impairment in the United States (2020-2060). Alzheimers Dement. (2021) 17:1966–75. doi: 10.1002/alz.12362
3. Alzheimer's disease facts and figures. Alzheimers Dement. (2022) 18:700–89. doi: 10.1002/alz.12638
4. Tanzi RE, Bertram L. Twenty years of the Alzheimer's disease amyloid hypothesis: a genetic perspective. Cell. (2005) 120:545–55. doi: 10.1016/j.cell.2005.02.008
5. Lewis J, Dickson DW, Lin WL, Chisholm L, Corral A, Jones G, et al. Enhanced neurofibrillary degeneration in transgenic mice expressing mutant tau and APP. Science. (2001) 293:1487–91. doi: 10.1126/science.1058189
6. Masters CL, Selkoe DJ. Biochemistry of amyloid β-protein and amyloid deposits in Alzheimer disease. Cold Spring Harb Perspect Med. (2012) 2:a006262. doi: 10.1101/cshperspect.a006262
7. Benilova I, Karran E, De Strooper B. The toxic Aβ oligomer and Alzheimer's disease: an emperor in need of clothes. Nat Neurosci. (2012) 15:349–57. doi: 10.1038/nn.3028
8. Zhu M, Wang X, Sun L, Schultzberg M, Hjorth E. Can inflammation be resolved in Alzheimer's disease? Ther Adv Neurol Disord. (2018) 11:1756286418791107. doi: 10.1177/1756286418791107
9. Blasko I, Veerhuis R, Stampfer-Kountchev M, Saurwein-Teissl M, Eikelenboom P, Grubeck-Loebenstein B. Costimulatory effects of interferon-gamma and interleukin-1beta or tumor necrosis factor alpha on the synthesis of Abeta1-40 and Abeta1-42 by human astrocytes. Neurobiol Dis. (2000) 7:682–9. doi: 10.1006/nbdi.2000.0321
10. Hu J, Akama KT, Krafft GA, Chromy BA, Van Eldik LJ. Amyloid-beta peptide activates cultured astrocytes: morphological alterations, cytokine induction and nitric oxide release. Brain Res. (1998) 785:195–206. doi: 10.1016/S0006-8993(97)01318-8
11. DiSabato DJ, Quan N, Godbout JP. Neuroinflammation: the devil is in the details. J Neurochem. (2016) 139 Suppl 2:136–53. doi: 10.1111/jnc.13607
12. Alves S, Churlaud G, Audrain M, Michaelsen-Preusse K, Fol R, Souchet B, et al. Interleukin-2 improves amyloid pathology, synaptic failure and memory in Alzheimer's disease mice. Brain. (2017) 140:826–42. doi: 10.1093/brain/aww330
13. Kiyota T, Okuyama S, Swan RJ, Jacobsen MT, Gendelman HE, Ikezu T, et al. CNS expression of anti-inflammatory cytokine interleukin-4 attenuates Alzheimer's disease-like pathogenesis in APP+PS1 bigenic mice. FASEB J. (2010) 24:3093–102. doi: 10.1096/fj.10-155317
14. 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:276. doi: 10.1186/s12974-018-1313-3
15. Roberts BR, Ryan TM, Bush AI, Masters CL, Duce JA. The role of metallobiology and amyloid-β peptides in Alzheimer's disease. J Neurochem. (2012) 120(Suppl 1):149–66. doi: 10.1111/j.1471-4159.2011.07500.x
16. Barnham KJ, Masters CL, Bush AI. Neurodegenerative diseases and oxidative stress. Nat Rev Drug Discov. (2004) 3:205–14. doi: 10.1038/nrd1330
17. Scheiber IF, Mercer JFB, Dringen R. Metabolism and functions of copper in brain. Prog Neurobiol. (2014) 116:33–57. doi: 10.1016/j.pneurobio.2014.01.002
18. Mercer SW, Wang J, Burke R. In vivo modeling of the pathogenic effect of copper transporter mutations that cause Menkes and Wilson Diseases, motor neuropathy, and susceptibility to Alzheimer's disease. J Biol Chem. (2017) 292:4113–22. doi: 10.1074/jbc.M116.756163
19. Ayton S, Lei P, Bush AI. Metallostasis in Alzheimer's disease. Free Radic Biol Med. (2013) 62:76–89. doi: 10.1016/j.freeradbiomed.2012.10.558
20. White AR, Reyes R, Mercer JF, Camakaris J, Zheng H, Bush AI, et al. Copper levels are increased in the cerebral cortex and liver of APP and APLP2 knockout mice. Brain Res. (1999) 842:439–44. doi: 10.1016/S0006-8993(99)01861-2
21. Maynard CJ, Cappai R, Volitakis I, Cherny RA, White AR, Beyreuther K, et al. Overexpression of Alzheimer's disease amyloid-beta opposes the age-dependent elevations of brain copper and iron. J Biol Chem. (2002) 277:44670–6. doi: 10.1074/jbc.M204379200
22. Ma QF, Li YM, Du JT, Kanazawa K, Nemoto T, Nakanishi H, et al. Binding of copper (II) ion to an Alzheimer's tau peptide as revealed by MALDI-TOF MS, CD, and NMR. Biopolymers. (2005) 79:74–85. doi: 10.1002/bip.20335
23. Zhou L-X, Du J-T, Zeng Z-Y, Wu W-H, Zhao Y-F, Kanazawa K, et al. Copper (II) modulates in vitro aggregation of a tau peptide. Peptides. (2007) 28:2229–34. doi: 10.1016/j.peptides.2007.08.022
24. Scali C, Prosperi C, Bracco L, Piccini C, Baronti R, Ginestroni A, et al. Neutrophils CD11b and fibroblasts PGE(2) are elevated in Alzheimer's disease. Neurobiol Aging. (2002) 23:523–30. doi: 10.1016/S0197-4580(01)00346-3
25. Zenaro E, Pietronigro E, Della Bianca V, 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:880–6. doi: 10.1038/nm.3913
26. Stock AJ, Kasus-Jacobi A, Pereira HA. The role of neutrophil granule proteins in neuroinflammation and Alzheimer's disease. J Neuroinflammation. (2018) 15:240. doi: 10.1186/s12974-018-1284-4
27. Crouch PJ, Hung LW, Adlard PA, Cortes M, Lal V, Filiz G, et al. Increasing Cu bioavailability inhibits Abeta oligomers and tau phosphorylation. Proc Natl Acad Sci U S A. (2009) 106:381–6. doi: 10.1073/pnas.0809057106
28. Tsvetkov P, Coy S, Petrova B, Dreishpoon M, Verma A, Abdusamad M, et al. Copper induces cell death by targeting lipoylated TCA cycle proteins. Science. (2022) 375:1254–61. doi: 10.1126/science.abf0529
29. 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:e47. doi: 10.1093/nar/gkv007
30. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. (2008) 9:559. doi: 10.1186/1471-2105-9-559
31. Tolles J, Meurer WJ. Logistic regression: relating patient characteristics to outcomes. JAMA. (2016) 316:533–4. doi: 10.1001/jama.2016.7653
32. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. doi: 10.1038/nmeth.3337
33. Wishart DS, Feunang YD, Guo AC, Lo EJ, Marcu A, Grant JR, et al. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res. (2018) 46:D1074–82. doi: 10.1093/nar/gkx1037
34. Finaurini S, Basilico N, Corbett Y, D'Alessandro S, Parapini S, Olliaro P, et al. Dihydroartemisinin inhibits the human erythroid cell differentiation by altering the cell cycle. Toxicology. (2012) 300:57–66. doi: 10.1016/j.tox.2012.05.024
35. Overington JP, Al-Lazikani B, Hopkins AL. How many drug targets are there? Nat Rev Drug Discov. (2006) 5:993–6. doi: 10.1038/nrd2199
36. Imming P, Sinning C, Meyer A. Drugs, their targets and the nature and number of drug targets. Nat Rev Drug Discov. (2006) 5:821–34. doi: 10.1038/nrd2132
37. Lee YC, Park CK, Kim MS, Kim JH. In vitro study for staining and toxicity of rose bengal on cultured bovine corneal endothelial cells. Cornea. (1996) 15:376–85. doi: 10.1097/00003226-199607000-00008
38. Tana C, Ticinesi A, Prati B, Nouvenne A, Meschi T. Uric acid and cognitive function in older individuals. Nutrients. (2018) 10:975. doi: 10.3390/nu10080975
39. Strosznajder AK, Wójtowicz S, Jezyna MJ, Sun GY, Strosznajder JB. Recent insights on the role of PPAR-β/δ in neuroinflammation and neurodegeneration, and its potential target for therapy. Neuromolecular Med. (2021) 23:86–98. doi: 10.1007/s12017-020-08629-9
40. Gerrits E, Brouwer N, Kooistra SM, Woodbury ME, Vermeiren Y, Lambourne M, et al. Distinct amyloid-β and tau-associated microglia profiles in Alzheimer's disease. Acta Neuropathol. (2021) 141:681–96. doi: 10.1007/s00401-021-02263-w
41. West LC, Cresswell P. Expanding roles for GILT in immunity. Curr Opin Immunol. (2013) 25:103–8. doi: 10.1016/j.coi.2012.11.006
42. Satoh J-I, Kino Y, Yanaizu M, Ishida T, Saito Y. Microglia express gamma-interferon-inducible lysosomal thiol reductase in the brains of Alzheimer's disease and Nasu-Hakola disease. Intractable Rare Dis Res. (2018) 7:251–7. doi: 10.5582/irdr.2018.01119
43. Chai K, Zhang X, Chen S, Gu H, Tang H, Cao P, et al. Application of weighted co-expression network analysis and machine learning to identify the pathological mechanism of Alzheimer's disease. Front Aging Neurosci. (2022) 14:837770. doi: 10.3389/fnagi.2022.837770
44. Di Paolo G, Kim T-W. Linking lipids to Alzheimer's disease: cholesterol and beyond. Nat Rev Neurosci. (2011) 12:284–96. doi: 10.1038/nrn3012
45. Penke B, Paragi G, Gera J, Berkecz R, Kovács Z, Crul T, et al. The role of lipids and membranes in the pathogenesis of Alzheimer's disease: a comprehensive view. Curr Alzheimer Res. (2018) 15:1191–212. doi: 10.2174/1567205015666180911151716
46. Ylilauri MPT, Voutilainen S, Lönnroos E, Virtanen HEK, Tuomainen T-P, Salonen JT, et al. Associations of dietary choline intake with risk of incident dementia and with cognitive performance: the Kuopio Ischaemic Heart Disease Risk Factor Study. Am J Clin Nutr. (2019) 110:1416–23. doi: 10.1093/ajcn/nqz148
47. Semba RD. Perspective: the potential role of circulating lysophosphatidylcholine in neuroprotection against Alzheimer disease. Adv Nutr. (2020) 11:760–72. doi: 10.1093/advances/nmaa024
48. Peters-Golden M, Brock TG. 5-lipoxygenase and FLAP. Prostaglandins Leukot Essent Fatty Acids. (2003) 69:99–109. doi: 10.1016/S0952-3278(03)00070-X
49. Chu J, Praticò D. Involvement of 5-lipoxygenase activating protein in the amyloidotic phenotype of an Alzheimer's disease mouse model. J Neuroinflammation. (2012) 9:127. doi: 10.1186/1742-2094-9-127
50. Chu J, Lauretti E, Di Meco A, Praticò D. FLAP pharmacological blockade modulates metabolism of endogenous tau in vivo. Transl Psychiatry. (2013) 3:e333. doi: 10.1038/tp.2013.106
51. Pascoal TA, Benedet AL, Ashton NJ, Kang MS, Therriault J, Chamoun M, et al. Microglial activation and tau propagate jointly across Braak stages. Nat Med. (2021) 27:1592–9. doi: 10.1038/s41591-021-01456-w
52. Dansokho C, Ait Ahmed D, Aid S, Toly-Ndour C, Chaigneau T, Calle V, et al. Regulatory T cells delay disease progression in Alzheimer-like pathology. Brain. (2016) 139:1237–51. doi: 10.1093/brain/awv408
53. Katayama H. Anti-interleukin-17A and anti-interleukin-23 antibodies may be effective against Alzheimer's disease: role of neutrophils in the pathogenesis. Brain Behav. (2020) 10:e01504. doi: 10.1002/brb3.1504
54. Carlini V, Verduci I, Cianci F, Cannavale G, Fenoglio C, Galimberti D, et al. CLIC1 protein accumulates in circulating monocyte membrane during neurodegeneration. Int J Mol Sci. (2020) 21:E1484. doi: 10.3390/ijms21041484
55. Sandin L, Bergkvist L, Nath S, Kielkopf C, Janefjord C, Helmfors L, et al. Beneficial effects of increased lysozyme levels in Alzheimer's disease modelled in Drosophila melanogaster. FEBS J. (2016) 283:3508–22. doi: 10.1111/febs.13830
56. Yang Y, Tapias V, Acosta D, Xu H, Chen H, Bhawal R, et al. Altered succinylation of mitochondrial proteins, APP and tau in Alzheimer's disease. Nat Commun. (2022) 13:159. doi: 10.1038/s41467-021-27572-2
57. Gaggelli E, Kozlowski H, Valensin D, Valensin G. Copper homeostasis and neurodegenerative disorders (Alzheimer's, prion, and Parkinson's diseases and amyotrophic lateral sclerosis). Chem Rev. (2006) 106:1995–2044. doi: 10.1021/cr040410w
58. Zhao Y, Long Z, Liu Y, Luo M, Qiu Y, Idris NFB, et al. Dihydroartemisinin ameliorates decreased neuroplasticity-associated proteins and excessive neuronal apoptosis in APP/PS1 mice. Curr Alzheimer Res. (2020) 17:916–925. doi: 10.2174/1567205017666201215124746
59. Xia L, Pang Y, Li J, Wu B, Du Y, Chen Y, et al. Dihydroartemisinin induces O-GlcNAcylation and improves cognitive function in a mouse model of tauopathy. J Alzheimers Dis. (2021) 84:239–48. doi: 10.3233/JAD-210643
60. Gao Y, Cui M, Zhong S, Feng C, Nwobodo AK, Chen B, et al. Dihydroartemisinin ameliorates LPS-induced neuroinflammation by inhibiting the PI3K/AKT pathway. Metab Brain Dis. (2020) 35:661–72. doi: 10.1007/s11011-020-00533-2
61. Cao D, Lu H, Lewis TL Li L. Intake of sucrose-sweetened water induces insulin resistance and exacerbates memory deficits and amyloidosis in a transgenic mouse model of Alzheimer disease *. J Biol Chem. (2007) 282:36275–82. doi: 10.1074/jbc.M703561200
62. Chen PE, Geballe MT, Stansfeld PJ, Johnston AR, Yuan H, Jacob AL, et al. Structural features of the glutamate binding site in recombinant NR1/NR2A N-methyl-D-aspartate receptors determined by site-directed mutagenesis and molecular modeling. Mol Pharmacol. (2005) 67:1470–84. doi: 10.1124/mol.104.008185
63. Kreamer BL, Siegel FL, Gourley GR. A novel inhibitor of beta-glucuronidase: L-aspartic acid. Pediatr Res. (2001) 50:460–6. doi: 10.1203/00006450-200110000-00007
64. Qin J, Kunda N, Qiao G, Calata JF, Pardiwala K, Prabhakar BS, et al. Colon cancer cell treatment with rose bengal generates a protective immune response via immunogenic cell death. Cell Death Dis. (2017) 8:e2584. doi: 10.1038/cddis.2016.473
65. Lee JS, Lee BI, Park CB. Photo-induced inhibition of Alzheimer's β-amyloid aggregation in vitro by rose bengal. Biomaterials. (2015) 38:43–9. doi: 10.1016/j.biomaterials.2014.10.058
Keywords: cuproptosis, Alzheimer's disease, diagnostic, drug, immune
Citation: Zhang E, Dai F, Chen T, Liu S, Xiao C and Shen X (2023) Diagnostic models and predictive drugs associated with cuproptosis hub genes in Alzheimer's disease. Front. Neurol. 13:1064639. doi: 10.3389/fneur.2022.1064639
Received: 08 October 2022; Accepted: 29 December 2022;
Published: 26 January 2023.
Edited by:
Alfredo Costa, University of Pavia, ItalyReviewed by:
Yanlin Bi, Qingdao Municipal Hospital, ChinaQisheng Su, Guangxi Medical University, China
Copyright © 2023 Zhang, Dai, Chen, Liu, Xiao and Shen. 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: Erdong Zhang, zhanged15@lzu.educn; Xiangchun Shen,
sxc@gmc.edu.cn
†These authors have contributed equally to this work