Abstract
Background
Although many studies reported a close relationship between depression and Alzheimer’s disease (AD), the underlying pathophysiological mechanism remains unclear. The present study aimed to investigate the mechanism of AD and major depressive disorder (MDD).
Keywords
Significance Statement
● The differentially expressed genes of AD-related datasets were significantly enriched in the long-term depression pathway. ● DYNC1H1, MAPRE3, TTBK2, ITGB1, and WASL might play essential roles in the pathophysiological mechanism of Alzheimer’s disease and major depressive disorder.
Introduction
Alzheimer’s disease (AD) is a central nervous system disease characterized by anterograde amnesia, accounting for 60 to 80 percent of all types of dementia. 1 Alzheimer’s Disease International (ADI) estimates that the number of people with dementia worldwide will increase to 152 million by 2050. The annual cost of dementia is estimated to be $2 trillion in 2030. 2 Alzheimer’s disease severely damages the ability of daily life and quality of life of patients and brings a substantial care burden and economic losses to caregivers and society. The US Food and Drug Administration (FDA) approved six drugs to alleviate the symptom development of AD patients,3,4 but the drug treatment effects were not satisfactory.
Major depressive disorder (MDD) is a complex disease characterized by mood depression and various physiological disorders, with high morbidity and suicide mortality. More than 264 million people worldwide suffer from depressive disorders, 5 and nearly 1 million people commit suicide each year due to mental disorder. 6 Globally, 163 million have MDD 5 ; the lifetime prevalence of MDD is approximately 15–18% and the 12-month MDD prevalence is 6%. 7 Genetic factors account for approximately 40% of the causes of MDD, 8 and chronic or nervous system diseases (e.g., AD, 9 Parkinson’s disease, 9 Huntington’s disease, 9 migraine, 10 glioblastoma multiforme, 11 and so forth) are also important factors affecting MDD. A recent network meta-analysis study has reported that MDD does exist in approximately 40% of AD cases. 12 Nevertheless, when AD coexists with MDD, the number of hospitalizations is considerably higher when compared to either MDD or AD alone.13,14 Most importantly, however, the mortality rate is increased which is experienced by individuals with comorbid MDD and AD. 15 Therefore, there is a high risk of comorbidities between MDD and AD.
A considerable amount of literature has revealed that depression may be a risk factor for dementia or part of the prodromal symptoms of dementia.16,17 Older patients with MDD exhibit significant cognitive deficits, which have been demonstrated in various cognitive domains in more than half of patients with late-life depression. 18 It has been reported that the prevalence of depression in patients with dementia is 29%.16,19 In addition, venlafaxine, a commonly used drug for the treatment of depression, has been reported to effectively improve the behavioral performance of rats in the Morris water maze test. 20 Although depression plays a significant role in the occurrence of AD, the pathological mechanism of the two comorbidities remains unclear.
In this study, we downloaded the gene expression profiles of AD, MDD, and MDD patients treated with venlafaxine for eight weeks. After determining the differentially expressed genes (DEGs) between diseased samples and healthy subjects, cluster analysis and functional enrichment analysis were carried out to explore the biological pathways in the processes of AD and MDD. Finally, the protein–protein interaction (PPI) network was constructed to find the potentially related genes between AD and MDD, and the hub genes were identified in MDD patients treated with venlafaxine. The identified hub genes may play an essential role in the genetics underpinning AD and MDD. These findings provide a reference for explaining the relationship between AD and MDD, but the specific pathological mechanism remains to be further studied.
Materials and Methods
Data Collection
Description of data used in this study.
①58 healthy human brain tissues included 25 hippocampal tissues and 33 frontal lobe tissues;
②40 brain tissues of AD patients included 19 hippocampal tissues and 21 frontal lobe tissues; ③24 healthy human brain tissues included 13 hippocampal tissues and 11 frontal lobe tissues; ④33 brain tissues of AD patients included 10 hippocampal tissues and 23 frontal lobe tissues; ⑤Patients with MDD were only treated with venlafaxine for 8 weeks.
Identification of Differentially Expressed Genes
We carried out a comprehensive bioinformatics analysis of the GEO dataset. The identification of DEGs was performed to compare the differences in gene expression between two conditions. The limma package 22 was used in GEO2R 23 to screen DEGs in the GEO dataset. To screen out appropriate number of DEGs and high-quality molecular markers of AD and MDD, we used different |log base 2 Fold Change| to screen DEGs and compared their results. When P < .05 and |log base 2 Fold Change | > 1, which identified 6 DEGs of AD-related datasets, 0 DEGs shared by AD and MDD were identified. When P < .05 and |log base 2 Fold Change | > .5, which identified 66 DEGs of AD-related datasets, 9 DEGs shared by AD and MDD were identified. When P < .05 and |log base 2 Fold Change | > .3, which identified 171 DEGs of AD-related datasets, 79 DEGs shared by AD and MDD were identified. Finally, we selected P < .05 and |log base 2 Fold Change | > .3 as the criteria for DEGs. 24 For visualization, the Venn diagram and the intersection plot of DEGs were generated using two web tools: (http://bioinformatics.psb.ugent.be/webtools/Venn/) and (http://www.bioinformatics.com.cn/plot_basic_upsetR_plot_009).
Functional Enrichment Analysis
To explore more comprehensive biological information for the above DEGs, GO enrichment analysis, including biological processes (BP), cellular component (CC), and molecular function (MF) terms, was performed on the R package cluster profiler (V3.16.0). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment was performed in the R packages cluster profiler (V3.16.0) and DOSE (V3.14.0). 25 The bubble map was drawn using the ggplot2 R package for the top 10 significant KEGG pathways based on p. adjust <.05. The GO enrichment results were divided into three functional categories, BP, CC, and MF, and the top 10 results of BP, CC, and MF (p. adjust <.05.) were visualized using the GO plot R package (V1.0.2).
Construction of the PPI Network and Identification of Hub Genes
The STRING 26 (http://string-db.org) online database was used to determine the interaction genes, and the interaction parameters were set to a maximum confidence of >.4. 27 To fully understand the functional interactions between proteins and select essential hub genes, the PPI network was visualized and analyzed with Cytoscape (version 3.7.2) software (http://www.cytoscape.org/), and hub genes were screened with cytoHubba, 28 a plug-in of Cytoscape software.
CytoHubba was used to rank genes in the PPI network from highest to lowest according to the maximum neighborhood component (MNC), density of maximum neighborhood component (DMNC), maximal clique centrality (MCC), and degree method scores (degree). The top 10 genes were selected in each algorithm, and the common genes were identified as hub genes. Finally, nine hub genes were obtained. 29
Validation of Hub Gene Signatures
To evaluate the clinical value of these 9 hub genes, we used one independent GEO cohort as a validation dataset. The GSE32280 dataset, including peripheral blood samples from 8 MDD patients before treatment and 8 MDD patients treated with venlafaxine for eight weeks, was used to demonstrate the expression of 9 hub genes. To show the expression patterns of 9 hub genes between MDD patients before treatment and MDD patients treated with venlafaxine, differences between each group were analyzed by T-test, and violin plot using the ggplot2 (V3.3.1) package in R (V4.0.0) was used for visualization.
Construction of Target Genes—Transcription Factor Regulatory Network
The miRNet database (https://www.mirnet.ca/) was used to predict TFs of significant target genes. Finally, Cytoscape was used to visualize the regulatory network of TF-target genes, which described the interaction between TFs and their potential targets in AD and MDD. 30
Construction of Target Genes-miRNA Regulatory Network
The miRNet database (https://www.mirnet.ca/) was used to predict miRNAs of significant target genes. Finally, Cytoscape was used to visualize the regulatory network of miRNA-target genes, which described the interaction between miRNAs and their potential targets in AD and MDD. 30
Receiver Operating Characteristic Curve Analysis
ROC curves were used to evaluate classifiers in bioinformatics applications. In order to further evaluate the predictive accuracy of hub genes, ROC analysis was performed to distinguish AD from normal control and MDD from normal control. Based on the expression profile by high-throughput sequencing data of hub genes, the pROC in R was used to generate the ROC curves of hub genes. The diagnostic value of hub gene was compared by measuring the area under ROC curve and AUC. 30
Results
Identification of Differentially Expressed Genes
In the hippocampal tissue of AD, there were 11188 DEGs (GSE48350 HIP, including 8128 upregulated genes and 3060 downregulated genes) in GSE48350 and 9539 DEGs (GSE5281 HIP, including 5605 upregulated genes and 3934 downregulated genes) in GSE5281. In the frontal lobe tissue of GSE48350, there were 6099 DEGs (GSE48350 FL, including 2279 upregulated genes and 3820 downregulated genes) and 10127 DEGs (GSE5281 FL, including 748 upregulated genes and 9379 downregulated genes) in GSE5281. There were 1866 DEGs (765 upregulated genes and 1101 downregulated genes) in GSE18309. After screening the DEGs of the GSE98793 dataset, a dataset related to MDD patients, including 64 MDD patients and 64 standard samples, we obtained 6846 DEGs, including 1904 upregulated genes and 4942 downregulated genes. Figure 1A shows the volcano plots of DEGs in each dataset drawn using the ggplot2 (V3.3.1) package in R (V4.0.0). Identification of DEGs among each GEO dataset. (A) The volcano plots present the expression pattern of DEGs in each dataset. Red: upregulated; Blue: downregulated; Black: normal. (B) The intersection plot of DEGs among these datasets. The vertical yellow line represents the intersection of DEGs (79 DEGs shared in AD and MDD datasets) from six datasets. (C) The Venn plot of DEGs (171 DEGs of AD-related datasets) among these datasets. HIP: hippocampal tissue; FL: frontal lobe tissue.
As shown in Figure 1B, the 79 DEGs (orange bar chart), the intersection of all the datasets (GSE48350 HIP, GSE5281 HIP, GSE48350 FL, GSE5281 FL, GSE18309, and GSE98793) in this study, were common to AD and MDD. Herein, 26 upregulated genes and 53 downregulated genes were found in AD, while 40 upregulated genes and 39 downregulated genes were found in MDD. In Figure 1C, the 171 DEGs are in common among all the AD datasets.
Functional Enrichment Analysis of AD-Related Differentially Expressed Genes
According to KEGG pathway enrichment analysis, AD-related DEGs were significantly enriched in the Rap1 signaling pathway, adrenergic signaling in cardiomyocytes, retrograde endocannabinoid signaling, glutamatergic synapse, cholinergic synapse, circadian entrainment, dilated cardiomyopathy (DCM), salivary secretion, and nicotine addiction. The long-term depression pathway was also enriched (Figures S1A and Table S1).
The DEG GO enrichment analysis results were divided into three functional categories (BP, CC, and MF). For BP, the DEGs were mainly enriched in adult behavior, synaptic vesicle cycle, regulation of membrane potential, regulation of postsynaptic membrane potential, synapse organization, hindbrain development, adult locomotory behavior, cognition, positive regulation of morphogenesis of an epithelium, and locomotory behavior (Figures S1B and Table S2). In terms of CC, the genes were mainly enriched in the glutamatergic synapse, presynapse, synaptic membrane, GABAergic synapse, ion channel complex, axon terminus, postsynaptic density, transmembrane transporter complex, asymmetric synapse, and transporter complex (Figures S1C and Table S2). In the MF group, the DEGs were significantly enriched in calmodulin-binding, ligand-gated anion channel activity, ion channel activity, substrate-specific channel activity, channel activity, passive transmembrane transporter activity, gated channel activity, cell adhesion molecule binding, extracellular ligand-gated ion channel activity, and alpha-actinin binding (Figures S1D and Table S2).
Functional Enrichment Analysis for the Differentially Expressed Genes Related to AD and Major Depressive Disorder
The top 10 KEGG pathways of the DEGs with overlapping AD and MDD were mainly enriched in vascular smooth muscle contraction, platelet activation, circadian entrainment, dilated cardiomyopathy (DCM), salivary secretion, gap junction, renin secretion, long-term depression, endocrine and other factor-regulated calcium reabsorption, and Hippo signaling pathway-multiple species (Figures 2A and Table S3). Functional characteristics analysis for the DEGs related to AD and MDD. (A) Top 10 KEGG pathway enrichment results. (B) The BP category of GO enrichment analyses. (C) The CC category of GO enrichment analyses. (D) The MF category of GO enrichment analyses.
The BP category of the GO analysis results showed that DEGs were significantly enriched in cytoplasmic mRNA processing body assembly, locomotory behavior, adult locomotory behavior, adult behavior, lamellipodial organization, Hippo signaling, regulation of microtubule-binding, regulation of mRNA catabolic process, regulation of Hippo signaling, and negative regulation of neuron death (Figures 2B and Table S4). For CC, these genes were enriched in the lamellipodial, cell leading edge, nuclear speck, voltage-gated potassium channel complex, endocytic vesicle, potassium channel complex, cell cortex, integrin complex, outer membrane, and spindle midzone (Figures 2C and Table S4). Moreover, they were significantly enriched in integrin binding, microtubule plus-end binding, dynein light intermediate chain binding, actin binding, transcription coactivator activity, alpha-tubulin binding, tubulin binding, kinesin binding, transcription cofactor binding, and tau protein binding in the MF categories (Figures 2D and Table S4).
PPI Network Analysis
To further explore the biological interaction between AD and MDD, the PPI network of DEGs was constructed based on the STRING database. The separated and partially connected nodes were removed. The PPI network was composed of 76 nodes and 44 edges, with an average local clustering coefficient of .38, and the PPI enrichment P-value was <.011 (Figure S2). The intersection of the four algorithms cytoHubba (MNC, DMNC, MCC, and Degree; Table S5) was used to determine hub genes (Figure 3E). Finally, nine hub genes were identified, namely, DYNC1H1, MAPRE3, TTBK2, ITGB1, WASL, PRKCA, ITSN1, CTTN, and REPS2 (Figures 3A–3E). Establishment of the protein-protein interaction network, TF-target gene network and miRNA-target gene network. (A-D) The hub genes were identified using four models (Degree, DMNC, MCC, and MNC) with the Cytoscape plug-in cytoHubba. The darker the red color of a point is, the higher the score of the gene in the algorithm is, indicating the gene is more important in the network. (E) A Venn diagram was used to identify the nine hub genes in AD and MDD (DYNC1H1, MAPRE3, TTBK2, ITGB1, WASL, PRKCA, ITSN1, CTTN, and REPS2). (F). Target gene-TF regulatory network between target genes and TFs. target genes are marked in red. (G). Target gene-miRNA regulatory network between target genes and miRNAs. target genes are marked in red.
Construction of TF—Target Regulatory Network
After combining the consequences of TF-target genes with the interactive network of TFs, the TFs that may control the target genes are shown in Figure 3F. The top five target genes were PPARA, ITGB1, IGF2BP2, MBP, and REPS2 in the TF-target regulatory network. PPARA interacts with 9 TFs such as RELA and others; ITGB1 interacts with 6 TFs such as FOXF2 and others; IGF2BP2 interacts with 3 TFs such as HMGA2 and others; MBP interacts with 3 TFs such as SOX10 and others; and REPS2 interacts with 2 TFs such as NFKB1 and others (Table S6).
Construction of miRNA—Target Regulatory Network
After combining the consequences of miRNA-target genes with the interaction network of miRNAs, the genes and miRNAs are shown in Figure 3G. The top five target genes were QKI, HIPK2, KMT2A, WASL, and IL6ST in the miRNA-target regulatory network. QKI interacts with 45 miRNAs such as hsa-mir-16-5p and others; HIPK2 interacts with 39 miRNAs such as hsa-mir-124-3p and others; KMT2A interacts with 33 miRNAs such as hsa-mir-155-5p and others; WASL interacts with 32 miRNAs such as hsa-mir-1-3p and others; and IL6ST interacts with 31 miRNAs such as hsa-mir-23b-3p and others (Table S7).
Changes in the Gene Expression Values of the Nine Hub Genes After Treatment With Venlafaxine in Major Depressive Disorder
To further validate the 9 potential hub genes in the relationship between AD and MDD, the expression of these genes was demonstrated in GSE32280 in a group of MDD patients treated with venlafaxine for eight weeks. The results showed that five of the nine hub genes, including DYNC1H1, MAPRE3, TTBK2, ITGB1, and WASL (Figures 4A–E), changed significantly in MDD patients after venlafaxine treatment (P<.05, T-test), while the gene expression values of PRKCA, ITSN1, CTTN, and REPS2 (Figures 4F–I) were not significantly different in MDD patients after treatment with venlafaxine (P>.05, T-test). As illustrated in Figure 4, compared with untreated control group, in the five hub genes whose expression changed significantly after venlafaxine treatment, the expression values of four hub genes were increased significantly and tended to return to normal levels, namely, DYNC1H1, MAPRE3, TTBK2, and WASL; and the expression value of ITGB1 was decreased significantly and tended to return to normal levels. More details of the five hub genes are listed in Table 2. This further supported the possibility that these five genes play an essential role in the process of AD and MDD. Changes in gene expression values of nine hub genes after treatment of MDD patients with venlafaxine. (A-I) represents the expression of DYNC1H1, MAPRE3, TTBK2, ITGB1, WASL, PRKCA, ITSN1, CTTN, and REPS2 genes in eight venlafaxine treated MDD patients (red) and eight control MDD patients (green). The hub genes in AD and MDD shared network. afrom the GeneCards database (www.genecards.org).
Receiver Operating Characteristic Curve Analysis
In addition, ROC curve analysis was implemented to calculate the capacity of the expression levels of hub genes to distinguish AD from normal control and MDD from normal control. The AUC value of DYNC1H1, MAPRE3, TTBK2, ITGB1, and WASL in AD indicated a moderate diagnostic efficiency (.7 > AUC >.5) (Figures 5A–E). The AUC value of DYNC1H1, and TTBK2 in MDD indicated a moderate diagnostic efficiency (.7 > AUC >.5) (Figures 5F and H). The results of ROC curve analysis suggested five hub genes exhibited a moderate diagnostic efficiency in AD and MDD. The AUC value of MAPRE3, ITGB1, and WASL in MDD suggested good diagnostic efficiency (.9 > AUC >.7) (Figures 5G, I, and J). The AUC value of combining MAPRE3 and ITGB1 was the highest in MDD and suggested good diagnostic efficiency (.9 > AUC >.7) (Figure 5K). ROC curve validated the sensitivity, specificity of five hub genes as a predictive biomarker for AD and MDD prognosis. (A). DYNC1H1 in AD; (B). MAPRE3 in AD; (C). TTBK2 in AD; (D). ITGB1 in AD; (E). WASL in AD; (F). DYNC1H1 in MDD; (G). MAPRE3 in MDD; (H). TTBK2 in MDD; (I). ITGB1 in MDD; (J). WASL in MDD.
Discussion
AD and MDD are common and frequently occurring diseases.2,31 At present, there is mounting evidence of a close relationship between AD and depression. It has been reported that the severity of depressive symptoms is positively correlated with whole brain cognitive impairment in patients with mild cognitive impairment (MCI), 32 while the elderly with severe depression and increasing depression symptoms have a higher risk of dementia. 33 The co-disease relationship between MDD and AD is exceptionally complex. At present, most of the relationship between AD and MDD is about clinical symptoms and co-disease, and more research on the common pathological mechanism of AD and MDD is needed. Therefore, the purpose of this study was to reveal the potential pathophysiological relationship between AD and MDD.
In this study, in order to mitigate the impact of sequencing batch effects and screen out molecular markers that can be shared by brain and blood tissues, 34 we restricted our choice of datasets to those produced from the same platform (GPL570 platform) and four gene expression profile datasets (GSE48350, GSE5281, GSE18309, and GSE98793) were used. Since there really is not suitable dataset related to brain tissue in MDD patients, we can only use the dataset related to blood tissue (https://www.ncbi.nlm.nih.gov/geo/). The samples of GSE18309 are blood samples with AD patients. GSE18309 acted as a bridge between AD and MDD’s samples. Although GSE18309 provides gene expression data from fewer samples compared to the other datasets, the data also provide a strong evidence for screening hub genes. For researchers, blood samples are the easiest to obtain, and we will examine the expression of hub genes in the blood tissues of patients with AD and MDD in the next experiment.
We used bioinformatics methods to integrate four gene expression profile datasets, which identified 171 DEGs (Table S8) of AD-related datasets; 79 DEGs (Table S9) shared by AD and MDD were identified. The main enrichment result of KEGG in AD-related DEGs was enriched in the long-term depression pathway provides a powerful reference in the cellular process for the relationship between AD and depression, which is indirectly consistent with previous studies showing that depression is closely related to AD.35-37 Seventy-nine DEGs shared by AD and MDD were mainly enriched in signaling pathways such as circadian entrainment, dilated cardiomyopathy (DCM), gap junction, and long-term depression, and GO enrichment results were enriched in adult behavior, endocytic vesicle, and protein binding. The role of circadian entrainment36,37 and gap junction38,39 in the KEGG pathway in AD and MDD has been reported, while the roles of dilated cardiomyopathy (DCM) 40 and the Hippo signaling pathway 41 in AD have been reported, but the roles of these pathways in MDD remain to be elucidated.
Subsequently, through the establishment of the PPI network by the STRING database and module analysis, nine hub genes were identified. Furthermore, after validating the expression of nine hub genes in MDD patients treated with venlafaxine, we found five genes (DYNC1H1, MAPRE3, TTBK2, ITGB1, and WASL) with significant changes in expression values (P < .05, T-test), which makes the expression level of these five hub genes in MDD patients closer to that in healthy people. Finally, we conducted the ROC curve analysis for five hub genes, and we found the ROC curves of MAPRE3 and ITGB1 looked good and the ROC curve of the combination of MAPRE3 and ITGB1 only improved slightly compared to either gene alone. We surmise that MAPRE3 and ITGB1 are in the same pathway and thus highly correlated in expression. The results indicate that these five hub genes had critical biological functions in the development of AD and MDD, and also suggest that antidepressant therapy may play a decisive role in the treatment of AD.
Dynein heavy chain 1 (DYNC1H1) dimer is the core of the cytoplasmic dynein 1 complex. The dynamic protein complex with housekeeping function is involved in the localization of the mitotic spindle, nuclear localization, golgi apparatus maintenance, and endosomal dynamics,42-44 while in the nervous system, the complex affects gene expression, development and regeneration, and the return of misfolded proteins and organelles from synapses to the cell body. 45 It has been reported that compared with the control group, the expression of DYNC1H1 in the CA1 region of the hippocampus in AD patients was significantly decreased, which may be due to the facts that the hyperphosphorylation of tau protein in the hippocampus of AD patients destroyed the stability of microtubules and affected the retrograde transport of DYNC1H1 along microtubules. 46 In addition, it has been reported that DYNC1H1 is associated with lower limb dominant spinal muscular atrophy (SMA-LED), 47 cortical dysplasia (MCD), 48 epileptic encephalopathy, 49 Parkinson’s disease, and hereditary spastic paraplegia. 50 Consequently, DYNC1H1 is likely to play an important role in AD and MDD.
MAPRE3 (microtubule-associated protein RP/EB family member 3) binds microtubule plus-ends to promote mature mushroom spines over immature thin spines in cultured neurons. 51 MAPRE3 has been reported to be associated with neuroblastoma. 52 The ADNP-SIRT1 complex consists of ADNP and SIRT1 and one of two major interaction points between ADNP and SIRT1 is at the level of EB1 (MAPRE1)/EB3 (MAPRE3)/Tau (MAPT) binding, shown to be significant for microtubule dynamics, axonal transport, synapse formation, and protection against tauopathy. The ADNP-SIRT1 complex expression is significantly reduced in brain of AD patients. So, MAPRE3 may participate in the occurrence of AD through the ADNP-SIRT1 complex. 53 Although MAPRE3 has not been reported in MDD, the expression value of MAPRE3 changed most obviously after the treatment of MDD patients with venlafaxine, which shows that MAPRE3 may be a crucial factor in the development of AD and MDD.
TTBK2, a kinase known to phosphorylate tau and tubulin, acts as a critical regulator of ciliogenesis. 54 TTBK2 is involved in tau/tubulin phosphorylation 55 and TDP-43 accumulation. 56 A study revealed that TTBK2 was directly or potentially associated with AD and neural development. 57 Although TTBK2 can phosphorylate two AD-related sites of tau, the primary phenotype of the mutant TTBK2 has nothing to do with AD but is associated with severe neurodegenerative disease, autosomal dominant spinocerebellar ataxia type 11 (SCA11). 58 TTBK2 is also associated with cancer progression.59,60 Combined with these previous reports, the results of our research further confirmed the association between the dysregulation of TTBK2 and the development of AD at the gene level. Nevertheless, to explain how this gene functions in MDD or the treatment of MDD, more experimental research is needed.
ITGB1 is a member of the integrin family. ITGB1 is also known as CD29. In mammals, most integrin subunits are expressed in different regions of the brain and are highly expressed in developing neurons. 61 Previous studies have shown that integrins, especially ITGB1 and ITGB3, play a role in nervous system development by regulating processes related to neuronal connections, such as axon growth and guidance, dendritic spine and synaptic formation and maintenance, and synaptic plasticity. 62 ITGB1 mRNA expression was changed in hippocampus and peripheral blood mononuclear cells (PBMCs) of AD patients, and researchers suggested using ITGB1 as an AD biomarker.34,63,64 In hippocampal neurons, ITGB1 and its downstream signaling pathway signals through the nonreceptor tyrosine kinase Arg (also known as Abelson tyrosine-protein kinase 2; ABL2) to regulate dendritic branches, synaptic plasticity, and behavior in the postnatal mouse hippocampus. 65 Furthermore, changes in the genetic code of integrin and other cell adhesion molecules (CAM) or proteins that mediate CAM signal transduction affect neural connections and are closely related to neurodevelopmental disorders, such as schizophrenia 66 and autism.67,68 Stress-associated epigenetic change of ITGB1 was correlated with depression. 69 In addition, ITGB1 70 is also associated with cancer progression.
WASL is also known as neural Wiskott–Aldrich syndrome protein (N-WASP). The expression of N-WASP is highest in nerve tissue. In developing hippocampal neurons, N-WASP is phosphorylated by FAK (focal adhesion kinase) to promote axonal growth 71 and is located in the growth cone by Nck1 and Cdc42 to promote expansion. 72 In addition to regulating axonal growth, N-WASP seems to be involved in the differentiation of neural stem cells by promoting an increase in filamentous foot processes. 73 In addition, N-WASP has been found to affect the myelin growth of Schwann cells, which wrap axons by using a giant lamellipodial sheet-like structure. 74 N-WASP, which regulates the reorganization and rearrangement of actin assembly, is highly expressed in brain tissue from AD patients, and may participate in the neurodegenerative aberrant sprouting in AD neurons. 75 Further research is required to determine the importance of WASL in AD and MDD.
Furthermore, the five hub genes in AD have been previously reported, while their roles in psychiatric disorders have been less reported. The five genes were novel DEGs identified in brain and blood samples from AD and MDD patients. Taken together, we speculate that AD and MDD are related not only in common genes and pathways but also in the function of proteins such as DYNC1H1, MAPRE3, TTBK2, ITGB1, and WASL.
Conclusions
In this study, bioinformatics analysis methods were used to screen and verify the shared and potential pathophysiologically related network hub genes in AD and MDD. DYNC1H1, MAPRE3, TTBK2, ITGB1, and WASL were identified as hub genes and they all associated with the progression of AD and MDD. Our findings may provide a reference for further exploration of the treatment of AD and MDD.
Supplemental Material
sj-pdf-1-aja-10.1177_15333175211046123 – Supplemental Material for Identification of Hub Genes Related to Alzheimer’s Disease and Major Depressive Disorder
Supplemental Material, sj-pdf-1-aja-10.1177_15333175211046123 for Identification of Hub Genes Related to Alzheimer’s Disease and Major Depressive Disorder by Yajing Cheng, Meiyue Sun, Feng Wang, Xin Geng and Fei Wang in American Journal of Alzheimer's Disease & Other Dementias
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by a Chinese National Natural Science Foundation Grant (No. 81671054, No. 81771135, No. 31771520, No. 91649102), the key project of Tianjin Natural Science Grant (19JCZDJC35600), Tianjin Natural Science Grant (19JCJQJC63500), and the Training program for young and middle-aged key Innovative Personnel (303078100407).
Supplemental Material
Supplemental material for this article is available online.
Appendix
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
