ORIGINAL RESEARCH article

Front. Aging Neurosci., 20 January 2022

Sec. Alzheimer's Disease and Related Dementias

Volume 13 - 2021 | https://doi.org/10.3389/fnagi.2021.789698

A Systematic Analysis on the Genes and Their Interaction Underlying the Comorbidity of Alzheimer's Disease and Major Depressive Disorder

  • 1. School of Biomedical Engineering, Tianjin Medical University, Tianjin, China

  • 2. School of Public Health, Tianjin Medical University, Tianjin, China

Article metrics

View details

11

Citations

3,8k

Views

1,5k

Downloads

Abstract

Background:

During the past years, clinical and epidemiological studies have indicated a close relationship between Alzheimer's disease (AD) and other mental disorders like major depressive disorder (MDD). At the same time, a number of genes genetically associated with AD or MDD have been detected. However, our knowledge on the mechanisms that link the two disorders is still incomplete, and controversies exist. In such a situation, a systematic analysis on these genes could provide clues to understand the molecular features of two disorders and their comorbidity.

Methods:

In this study, we compiled the genes reported to be associated with AD or MDD by a comprehensive search of human genetic studies and genes curated in disease-related database. Then, we investigated the features of the shared genes between AD and MDD using the functional enrichment analysis. Furthermore, the major biochemical pathways enriched in the AD- or MDD-associated genes were identified, and the cross talks between the pathways were analyzed. In addition, novel candidate genes related to AD and MDD were predicted in the context of human protein-protein interactome.

Results:

We obtained 650 AD-associated genes, 447 MDD-associated genes, and 77 shared genes between AD and MDD. The functional analysis revealed that biological processes involved in cognition, neural development, synaptic transmission, and immune-related processes were enriched in the common genes, indicating a complex mechanism underlying the comorbidity of the two diseases. In addition, we conducted the pathway enrichment analysis and found 102 shared pathways between AD and MDD, which involved in neuronal development, endocrine, cell growth, and immune response. By using the pathway cross-talk analysis, we found that these pathways could be roughly clustered into four modules, i.e., the immune response-related module, the neurodevelopmental module, the cancer or cell growth module, and the endocrine module. Furthermore, we obtained 37 novel candidate genes potentially related to AD and MDD with node degrees > 5.0 by mapping the shared genes to human protein-protein interaction network (PPIN). Finally, we found that 37 novel candidate genes are significantly expressed in the brain.

Conclusion:

These results indicated shared biological processes and pathways between AD and MDD and provided hints for the comorbidity of AD and MDD.

Introduction

As the most common type of dementia, Alzheimer's disease (AD) affects ~10% of people of 65 years or older (Blennow et al., 2006). Although cognitive impairment is the core feature of this disease, neuropsychiatric symptoms appear commonly in patients with AD and are persistent (Fernández et al., 2010; Kales et al., 2015). Depression and depression symptoms are among the most common symptoms in patients with AD and may be the early symptoms of the dementing process (Lyketsos et al., 1997b). Depression or depressive symptoms may affect as much as 50% of patients with AD (Starkstein et al., 2005), and about 20–30% of patients with AD are influenced by major depressive disorder (MDD) (Ballard et al., 1996; Enache et al., 2011). The continuous increase of aging populations across the world has led to an increasing prevalence of cognitive impairment, which has become a major problem both clinically and socially (Prince et al., 2016; The Alzheimer's Association, 2020). Specifically, AD and MDD are both common diseases of the elderly and are often comorbidity (Lee and Lyketsos, 2003). Thus, exploring the association between AD and MDD will not only help us to elucidate the neurobiological mechanisms underlying the two diseases but also be useful for developing better therapeutic approaches for them.

However, the relationship between AD and MDD is complex and not well-understood. Depression and depressive symptoms not only are common and highly persistent in patients with AD but also may increase the risk of functional decline and behavioral disturbance (Migliorelli et al., 1995; Ballard et al., 1996; Lyketsos et al., 1997b; Olin et al., 2002). It is reported that depressed patients may be more likely to develop AD than people without the symptoms (Saczynski et al., 2010). In contrast, no connection between the two conditions has been found in some studies. For example, in a cohort including 4,615 elderly Canadians, the history of depression is not significantly associated with the risk of dementia (Lindsay et al., 2002); in another study based on people aged 85 years or over, it is reported that depression is associated with cognitive impairment, but depression alone may not increase the risk of cognitive decline (Vinkers et al., 2004).

At the molecular level, both AD and MDD are diseases influenced by multiple genetic factors. During the past years, a number of genes, such as apolipoprotein E (APOE), ATP binding cassette subfamily A member 7 (ABCA7), and others (Steinberg et al., 2015; Serrano-Pozo et al., 2021), have been found to be genetically associated with AD. Some studies found the shared genetic risk factors could partially explain the observed association between AD and MDD. A positive association between ApoE4 allele and late-onset depression has been reported in some studies (Krishnan et al., 1996; Kim et al., 2002). However, in other studies, no association between ApoE4 allele and depression in patients with AD was detected (Lyketsos et al., 1997a; Scarmeas et al., 2002). Since both AD and MDD are cognitive disorders, hippocampus plays a principal role in memory deficits in patients with AD and MDD. In patients with AD, the decrease in hippocampal volume is related to cognitive decline and other symptoms of AD neuropathology (Jack et al., 1999; van der Flier and Scheltens, 2009). Furthermore, a correlation between hippocampal volume loss in the presence of ApoE allele epsilon4 and decreased cerebrospinal fluid beta-amyloid, the biomarkers of AD, has been observed (Schuff et al., 2009). Structural studies also have found that hippocampal volume is smaller in patients with MDD compared with people without depression, and hippocampal atrophy is suggested as a biomarker for cognitive decline in patients with MDD (Taylor et al., 2014; Schmaal et al., 2017). Amyloid-β (Aβ) deposition is an important pathological change in the AD brain, and it has been found that lifelong major depression is associated with the deposition of Aβ in the brain (Wu et al., 2014). Abnormal neuroinflammation is related to the pathological features of AD and has a negative impact on cognitive function (Cai et al., 2019; Passamonti et al., 2019). Studies have revealed an association between neuroinflammation and the pathogenesis of MDD or AD (Barber, 2011); however, the impact of neuroinflammation is much larger and more complicated, and further exploration is necessary. The study found that late-life major depression with dementia was significantly associated with a neuropathologic diagnosis of AD (Sweet et al., 2004). However, some studies report no association between lower hippocampal volume and AD pathology in late-life depression (De Winter et al., 2017).

Therefore, the correlation between AD and depression is still not fully understood. Obviously, it is necessary to further explore the biological mechanisms underlying the comorbidity of AD and MDD, which may help us design better therapeutic approach for the diseases.

In current study, we comprehensively collected genes related to AD or MDD. We next employed the functional analysis to detect the major biological processes and biochemical pathways and further investigated the interactions among the significantly enriched pathways. Of significance, we predicted some novel candidate genes potentially related to AD and MDD by analyzing human protein-protein interaction network (PPIN). This study will provide valuable clues for understanding the molecular mechanisms underlying the comorbidity of AD and MDD.

Materials and Methods

Data Collection

The DisGeNET (https://www.disgenet.org/) is a database that comprehensively integrates human disease-associated genes and their variants (Piñero et al., 2020). It compiles data from expert curated repositories, text mining data extracted from scientific literature, experimentally validated data, and referred data. A gene with a gene-disease associations (GDA) score above 0.2 indicates a strong association with the disease. In the current study, only genes with a GDA score >0.2 were included in the analysis. A list of 146 AD-associated genes was retrieved from DisGeNET. Furthermore, in an earlier study (Hu et al., 2017), we collected 430 genes related to AD from reports deposited at PubMed (https://pubmed.ncbi.nlm.nih.gov/). These genes were identified from either candidate gene association studies (CGAS) or GWAS on human samples. This gene set was updated to include the genes reported in the recent publications, and a total of 587 genes associated with AD were obtained. Then, this gene list was merged with those retrieved from DisGeNET, a list of 650 unique AD-associated genes were left.

For MDD, we retrieved 266 genes with a GDA score >0.2 from the DisGeNET database. Additionally, Fan et al. (2020) compiled a list of 255 MDD-associated genes. By integrating the two gene sets, we obtained a total of 447 unique MDD-associated genes.

Functional Enrichment Analysis

The gene ontology (GO) and pathway enriched in the shared genes were annotated using the WebGestalt (http://www.webgestalt.org/option.php) (Zhang et al., 2005; Wang et al., 2013). The Benjamini and Hochberg (B&H) method for false discovery rate correction was applied to correct the original p-value. WebGestalt is an online-based system to integrate and visualize gene set in various biological contexts including GO. Only the GO terms with FDR value smaller than 0.05 were chosen as significantly enriched terms in our study. The biochemical pathways enriched in the 650 AD-associated genes and the 447 MDD-associated genes were identified by ToppGene (https://toppgene.cchmc.org/enrichment.jsp) (Chen et al., 2009), which is an online tool that integrates different pathway databases such as Kyoto Encyclopedia of Genes and Genomes (KEGG) (Aoki-Kinoshita and Kanehisa, 2007), Pathway Interaction Database (PID) (Schaefer et al., 2009), and Reactome (Croft et al., 2011). For the pathway analysis, only KEGG was chosen as the pathway database, and pathway with FDR < 0.05 was defined as significant pathway.

Cross-Talk Analysis Among Pathways

Here, we first obtained pathways enriched in the AD-associated genes or MDD-associated genes, respectively. We selected the overlapping pathways of AD and MDD to construct the pathway cross-talk network by integrating AD-associated genes and MDD-associated genes. As described in our previous report (Guo et al., 2021), we adopted two measurements (Jia et al., 2011; Liu et al., 2015), i.e., the Jaccard Coefficient (JC) and the Overlap Coefficient (OC) , with A and B being the lists of genes included in the two tested pathways, and |A| and |B| representing the number of genes contained in the two pathways. Finally, the cross talk between pathways was constructed and visualized using the Cytoscape 3.7.1 software (Kohl et al., 2011).

Prediction of New Candidate Genes Related to AD and MDD

In this study, the PPIN data were obtained from the Protein Interaction Network Analysis (PINA) database (Cowley et al., 2012) by pooling and curating the unique physical interaction information from six main public protein interaction databases, i.e., BioGRID, IntAct, DIP, MINT, MIPS/MPact, and HPRD. In the meantime, Menche et al. (2015) reported an interactome for human that contained 141,296 edges among 13,460 nodes. Later, we merged the two interactomes data by excluding the self-interaction and redundant pairs. We excluded ubiquitin-C, B, and D (i.e., UBC, UBB, and UBD) from this network because of the non-specific binding of ubiquitin to proteins for degradation (Ferrari et al., 2018). Finally, we obtained a relatively complete human physical interactome, which included 15,435 genes/protein and 218,161 interactions (Supplementary Material 1).

We mapped the 650 AD-associated genes to the PPIN and obtained an AD-specific network containing 9,089 genes/protein and 32,065 interactions. Similarly, by mapping the 447 MDD-associated genes to the PPIN, we obtained an MDD-specific network containing 6,541 genes/protein and 18,221 interactions. Then, the irrelevant interactions were excluded by merging the AD-specific network and MDD-specific network into a single disease network, which contained 10,203 genes/protein and 45,385 interactions. Finally, the AD-specific network was compared against the MDD-specific network, and the overlapping part was extracted, and the node degree of each node in this overlapping network was calculated with the disease network as background via the “Network Analyzer” (Assenov et al., 2008) plug-in in the Cytoscape software. After removing AD-associated genes and MDD-associated genes, the remaining genes could be expected to be novel candidate genes for both diseases.

Analysis of Temporal and Cell Type-Specific Expression

The brain is composed of different but interdependent regions, with which being made up of many different types of cells performing specific functions. The precise spatial and temporal regulation of gene expression in these cells plays key roles in the function of the brain. To explore the biological features of new candidate genes in a specific brain cell type, brain region, and developmental stages, the corresponding gene expression pattern was analyzed. Cell-type specific expression analysis (CSEA; http://genetics.wustl.edu/jdlab/csea-tool-2/) (Xu et al., 2014) was used to compare the expression pattern of new candidate genes in different brain cells or brain regions. The Fisher's exact test was performed across the cell types and developmental stages at a specificity index threshold (pSI) of 0.05. This threshold measures the possibility of the expression of a specific gene in the given cell type, brain region, and developmental stage.

Results

Susceptibility Genes of AD and MDD

By searching the available genetic studies on AD and collecting the genes curated in disease-related database DisGeNET, we compiled a list of 650 AD-associated genes (Supplementary Table S1; the gene list is referred to as ADgset, hereafter). Among them were those encoding apolipoprotein, ATP-binding cassette transporters, nicotinic acetylcholine receptors, adrenoceptors, serotonin receptors, dopamine degradation, dopamine receptors, complement receptors, etc. Similarly, we also compiled a list of 447 MDD-associated genes (Supplementary Table S2; the gene list is referred to as MDDgset, hereafter). The MDDgset included genes related to immune regulation, dopamine neurotransmission, serotonin neurotransmission, embryonic development, cellular stress response, etc.

There were 77 shared genes between the two diseases (Table 1), among which included five oxidative stress-related genes (i.e., HIF1A, HSPA4, MAOA, MAOB, SOD2), four interleukin genes (i.e., IL10, IL1B, IL33, IL6), four neural factors genes (i.e., BDNF, NGF, NPY, NTF3), two 5-hydroxytryptamine receptor genes (i.e., HTR2A, HTR6), two solute carrier genes (i.e., SLC6A3, SLC6A4), two estrogen receptor genes (i.e., ESR1, ESR2), two cancer-related genes (i.e., TNF, TP53), two LDL receptor-related protein coding genes (i.e., LRP1, LRP8), one transcription factors gene (i.e., TFCP2), one apolipoprotein gene (i.e., APOE), one dopamine receptor gene (i.e., DRD4), and others.

Table 1

Gene IDGene symbolGene nameGene source
5243ABCB1ATP binding cassette subfamily B member 1Literature/DisGeNet
1636ACEAngiotensin I converting enzymeLiterature/DisGeNet
155ADRB3Adrenoceptor beta 3Literature/DisGeNet
348APOEApolipoprotein ELiterature/DisGeNet
406ARNTLAryl hydrocarbon receptor nuclear translocator likeLiterature
540ATP7BATPase copper transporting betaLiterature/DisGeNet
613BCRBCR, RhoGEF and GTPase activating proteinLiterature/DisGeNet
627BDNFBrain derived neurotrophic factorLiterature/DisGeNet
6347CCL2C-C motif chemokine ligand 2Literature/DisGeNet
1029CDKN2ACyclin dependent kinase inhibitor 2ALiterature/DisGeNet
89832CHRFAM7ACHRNA7 (exons 5–10) and FAM7A (exons A-E) fusionLiterature
9575CLOCKClock circadian regulatorLiterature/DisGeNet
26047CNTNAP2Contactin associated protein 2Literature/DisGeNet
1312COMTcatechol-O-methyltransferaseLiterature/DisGeNet
1392CRHCorticotropin releasing hormoneDisGeNet
1565CYP2D6Cytochrome P450 family 2 subfamily D member 6Literature/DisGeNet
267012DAOAD-amino acid oxidase activatorLiterature/DisGeNet
1621DBHDopamine beta-hydroxylaseLiterature/DisGeNet
27185DISC1Disrupted in schizophrenia 1Literature/DisGeNet
1740DLG2Disks large MAGUK scaffold protein 2Literature/DisGeNet
1742DLG4Disks large MAGUK scaffold protein 4Literature/DisGeNet
1815DRD4Dopamine receptor D4Literature/DisGeNet
5610EIF2AK2Eukaryotic translation initiation factor 2 alpha kinase 2Literature
2099ESR1Estrogen receptor 1Literature/DisGeNet
2100ESR2Estrogen receptor 2Literature/DisGeNet
79068FTOAlpha-ketoglutarate dependent dioxygenaseLiterature
2692GHRHRGrowth hormone releasing hormone receptorLiterature
2784GNB3G protein subunit beta 3Literature/DisGeNet
2904GRIN2BGlutamate ionotropic receptor NMDA type subunit 2BLiterature
2932GSK3BGlycogen synthase kinase 3 betaLiterature/DisGeNet
3091HIF1AHypoxia inducible factor 1 subunit alphaLiterature/DisGeNet
3240HPHaptoglobinLiterature/DisGeNet
3308HSPA4Heat shock protein family A (Hsp70) member 4Literature/DisGeNet
3356HTR2A5-hydroxytryptamine receptor 2ALiterature/DisGeNet
3362HTR65-hydroxytryptamine receptor 6Literature
3586IL10Interleukin 10Literature/DisGeNet
3553IL1BInterleukin 1 betaLiterature/DisGeNet
90865IL33Interleukin 33Literature/DisGeNet
3569IL6Interleukin 6Literature/DisGeNet
3952LEPLeptinDisGeNet
4035LRP1LDL receptor related protein 1Literature/DisGeNet
7804LRP8LDL receptor related protein 8Literature/DisGeNet
4128MAOAMonoamine oxidase ALiterature/DisGeNet
4129MAOBMonoamine oxidase BDisGeNet
4133MAP2Microtubule associated protein 2DisGeNet
9175MAP3K13Mitogen-activated protein kinase 13Literature
4353MPOMyeloperoxidaseLiterature/DisGeNet
4524MTHFRMethylenetetrahydrofolate reductaseLiterature/DisGeNet
4599MX1MX dynamin like GTPase 1Literature
4803NGFNerve growth factorLiterature/DisGeNet
4852NPYNeuropeptide YLiterature/DisGeNet
4908NTF3Neurotrophin 3Literature/DisGeNet
4915NTRK2Neurotrophic receptor tyrosine kinase 2Literature/DisGeNet
84547PGBD1piggyBac transposable element derived 1Literature
65018PINK1PTEN induced kinase 1Literature/DisGeNet
5337PLD1Phospholipase D1Literature
5444PON1Paraoxonase 1Literature
5743PTGS2Prostaglandin-endoperoxide synthase 2Literature/DisGeNet
5649RELNReelinLiterature/DisGeNet
6285S100BS100 calcium binding protein BLiterature/DisGeNet
9037SEMA5ASemaphorin 5ALiterature
10280SIGMAR1Sigma non-opioid intracellular receptor 1Literature
23411SIRT1Sirtuin 1Literature
6531SLC6A3Solute carrier family 6 member 3Literature
6532SLC6A4solute carrier family 6 member 4Literature/DisGeNet
6648SOD2Superoxide dismutase 2Literature/DisGeNet
6671SP4Sp4 transcription factorLiterature
6688SPI1Spi-1 proto-oncogeneLiterature
6750SSTSomatostatinLiterature/DisGeNet
7024TFCP2Transcription factor CP2Literature/DisGeNet
7040TGFB1Transforming growth factor beta 1Literature
7124TNFTumor necrosis factorLiterature/DisGeNet
10452TOMM40Translocase of outer mitochondrial membrane 40Literature/DisGeNet
7157TP53Tumor protein p53Literature
7422VEGFAVascular endothelial growth factor ALiterature/DisGeNet
7434VIPR2Vasoactive intestinal peptide receptor 2Literature
91752ZNF804AZinc finger protein 804ALiterature

Susceptibility gene shared by Alzheimer's disease (AD) and major depressive disorder (MDD).

Biological Process Enrichment Analysis of Shared Genes

For these 77 shared genes, 74 were significantly enriched in the biological regulation, 47 in membrane, and 67 in the protein binding (Figure 1). Totally, 311 BP GO terms (Supplementary Table S3) were significantly enriched in the shared genes. The most significantly term is cognition, which is in line with the fact that patients with AD and MDD have cognitive impairments (Sierksma et al., 2010). Furthermore, biological processes related to neural development and synaptic transmission were enriched, such as neuron death, dopaminergic neuron differentiation, neuron migration, synaptic transmission, dopaminergic, chemical synaptic transmission, postsynaptic and synaptic transmission, and glutamatergic. Terms associated with drug reactions and metabolic processes were overrepresented, such as response to toxic substance, response to antibiotic, response to oxidative stress, reactive oxygen species metabolic process, hormone metabolic process, receptor metabolic process, steroid metabolic process, and amine metabolic process. Of significance, immune-related processes, such as neuroinflammatory response, inflammatory cell apoptotic process, and acute inflammatory response were significantly over-represented.

Figure 1

Figure 1

Functional enrichment analysis of the 77 shared genes. Gene ontology (GO) items belonging to biological process (BP), cellular component (CC), and molecular function (MF) are shown in gray, white, and black columns, respectively.

Pathway Enrichment Analysis of ADgset and MDDgset

For the ADgset, 153 significantly enriched pathways were identified (Supplementary Table S4), while for MDDgset, 146 significantly enriched pathways were identified (Supplementary Table S5). We obtained 102 shared pathways by integrating pathways enriched in ADgset and MDDgset (Table 2). Beside the shared pathways, there were 95 disease-specific pathways for the two disorders, which included 51 AD-specific pathways and 44 MDD-specific pathways. Among the 102 shared pathways, several pathways were involved in the neural function and neurotransmission, for example, cholinergic synapse, dopaminergic synapse, serotonergic synapse, neuroactive ligand-receptor interaction, and neurotrophin signaling pathway. Some pathways were related to neurological disorders, such as Huntington's disease and AD. At the same time, substance addiction-related pathways were identified, such as amphetamine addiction and cocaine addiction. Metabolism-related pathways such as drug metabolism (cytochrome P450) and tyrosine metabolism were also included, indicating these or related processes may play roles in the pathogenesis of AD and MDD. In addition, about one-third of the shared pathways were related to immune response, immune disease, and pathogens infection, such as inflammatory bowel disease, inflammatory mediator regulation of TRP channels, Th17 cell differentiation, T-cell receptor signaling pathway, Epstein-Barr virus infection, Herpes simplex infection, and HTLV-I infection, suggesting the immune response was involved in the pathogenesis of AD and MDD. Moreover, among the 102 shared pathways, some pathways were related to cell growth and/or survival, including PI3K-Akt signaling, Jak-STAT signaling pathway, MAPK signaling pathway, and mTOR signaling. Among the list, there were also pathways involved in cancer (e.g., bladder cancer, breast cancer, colorectal cancer, endometrial cancer, pancreatic cancer, prostate cancer, and small cell lung cancer), indicating the complicated pathological mechanisms of AD and MDD.

Table 2

ModulePathwayDegreeFDR in ADgsetFDR in MDDgsetNumber of genes
ADMDDShared
Immune responseHepatitis B798.37 × 10−051.08 × 10−0519174
HTLV-I infection772.43 × 10−037.03 × 10−0623247
Fluid shear stress and atherosclerosis761.74 × 10−111.19 × 10−0330135
Influenza A761.59 × 10−091.41 × 10−0930258
Tuberculosis721.43 × 10−138.03 × 10−0337135
TNF signaling pathway691.48 × 10−061.01 × 10−0419135
Chagas disease (American trypanosomiasis)684.21 × 10−113.32 × 10−0625157
Epstein-Barr virus infection664.77 × 10−053.79 × 10−0324155
Toll-like receptor signaling pathway665.89 × 10−056.90 × 10−0516133
T cell receptor signaling pathway631.20 × 10−023.10 × 10−0311103
Hepatitis C636.48 × 10−041.77 × 10−0416144
Measles629.96 × 10−051.46 × 10−0918226
Toxoplasmosis603.27 × 10−105.67 × 10−0425123
Amoebiasis565.45 × 10−085.39 × 10−0420115
NOD-like receptor signaling pathway558.97 × 10−062.40 × 10−0423164
Inflammatory mediator regulation of TRP channels522.88 × 10−032.06 × 10−0612153
Herpes simplex infection512.55 × 10−086.44 × 10−0529188
Autophagy - animal505.13 × 10−043.25 × 10−021691
IL-17 signaling pathway497.90 × 10−071.01 × 10−0418126
Th17 cell differentiation495.74 × 10−062.21 × 10−0518144
Cytokine-cytokine receptor interaction441.04 × 10−035.53 × 10−0625258
Leishmaniasis405.49 × 10−102.49 × 10−0420105
Inflammatory bowel disease (IBD)383.38 × 10−091.01 × 10−0418105
Th1 and Th2 cell differentiation362.03 × 10−031.48 × 10−021280
Rheumatoid arthritis363.38 × 10−093.14 × 10−0421116
RIG-I-like receptor signaling pathway352.54 × 10−031.18 × 10−021071
Legionellosis343.34 × 10−127.58 × 10−042083
African trypanosomiasis321.39 × 10−065.84 × 10−0711104
Malaria293.27 × 10−103.62 × 10−041787
Graft-versus-host disease267.18 × 10−061.61 × 10−021153
Endocytosis211.11 × 10−021.38 × 10−0221163
Antigen processing and presentation174.61 × 10−045.66 × 10−031282
Prion diseases64.35 × 10−048.85 × 10−03852
NeurodevelopmentNeurotrophin signaling pathway672.33 × 10−088.36 × 10−0423126
Focal adhesion641.36 × 10−027.74 × 10−0317143
Glioma584.66 × 10−033.27 × 10−069122
Cholinergic synapse581.18 × 10−035.62 × 10−1114221
Dopaminergic synapse552.20 × 10−041.87 × 10−18173210
Longevity regulating pathway533.82 × 10−056.98 × 10−0715153
Longevity regulating pathway—multiple species401.10 × 10−034.75 × 10−0910152
Calcium signaling pathway372.80 × 10−031.76 × 10−1318313
Alzheimer's disease352.37 × 10−177.61 × 10−0441156
Serotonergic synapse294.49 × 10−043.36 × 10−1915318
Amphetamine addiction246.77 × 10−032.06 × 10−149215
Long-term depression233.02 × 10−025.37 × 10−057101
Amyotrophic lateral sclerosis (ALS)239.80 × 10−062.41 × 10−0612114
Huntington's disease146.14 × 10−051.38 × 10−0223135
Neuroactive ligand-receptor interaction123.00 × 10−031.33 × 10−1924478
Cocaine addiction108.31 × 10−041.44 × 10−119166
Drug metabolism—cytochrome P45042.54 × 10−037.90 × 10−041093
Metabolism of xenobiotics by cytochrome P45031.11 × 10−024.78 × 10−02961
Catecholamine biosynthesis, tyrosine = > dopamine = > noradrenaline = > adrenaline13.55 × 10−027.82 × 10−04231
Cancer/cell growthPathways in cancer871.73 × 10−091.83 × 10−07483512
MAPK signaling pathway814.77 × 10−042.04 × 10−1025329
Proteoglycans in cancer791.74 × 10−072.37 × 10−0629226
PI3K-Akt signaling pathway758.50 × 10−071.24 × 10−0538287
HIF-1 signaling pathway742.55 × 10−081.64 × 10−0821183
Osteoclast differentiation706.00 × 10−041.77 × 10−0316124
Ras signaling pathway705.22 × 10−047.06 × 10−0423185
FoxO signaling pathway688.55 × 10−051.90 × 10−0418145
Apoptosis671.04 × 10−062.01 × 10−0222103
Rap1 signaling pathway662.53 × 10−033.48 × 10−0620223
Prostate cancer633.98 × 10−045.51 × 10−0513122
cAMP signaling pathway612.58 × 10−029.35 × 10−1516346
Pancreatic cancer614.66 × 10−033.27 × 10−069125
EGFR tyrosine kinase inhibitor resistance601.63 × 10−041.51 × 10−0713153
mTOR signaling pathway606.16 × 10−032.00 × 10−0315132
Jak-STAT signaling pathway591.75 × 10−022.41 × 10−0614193
MicroRNAs in cancer591.96 × 10−032.20 × 10−0226176
Colorectal cancer572.22 × 10−045.43 × 10−031173
Choline metabolism in cancer549.21 × 10−031.77 × 10−0411122
Small cell lung cancer508.61 × 10−032.77 × 10−021072
Breast cancer492.11 × 10−021.35 × 10−0413154
Non-small cell lung cancer452.59 × 10−022.15 × 10−04792
cGMP-PKG signaling pathway435.19 × 10−031.35 × 10−0516181
Endometrial cancer403.76 × 10−032.00 × 10−03872
AMPK signaling pathway365.00 × 10−078.96 × 10−0321102
p53 signaling pathway232.42 × 10−031.11 × 10−021072
Bladder cancer231.15 × 10−036.87 × 10−04873
Chemical carcinogenesis32.53 × 10−038.05 × 10−031181
EndocrineAGE-RAGE signaling pathway in diabetic complications764.26 × 10−071.06 × 10−0519146
Non-alcoholic fatty liver disease (NAFLD)712.82 × 10−121.32 × 10−0232116
Sphingolipid signaling pathway672.28 × 10−047.74 × 10−0316103
Endocrine resistance658.55 × 10−057.19 × 10−0915185
Insulin resistance653.23 × 10−074.06 × 10−0320103
Prolactin signaling pathway612.54 × 10−033.70 × 10−0510113
Phospholipase D signaling pathway592.78 × 10−043.48 × 10−0618181
Acute myeloid leukemia522.04 × 10−027.58 × 10−04781
Chronic myeloid leukemia498.82 × 10−038.62 × 10−04994
VEGF signaling pathway462.80 × 10−035.03 × 10−03972
GnRH signaling pathway453.64 × 10−021.07 × 10−069151
Apelin signaling pathway456.89 × 10−036.51 × 10−0614171
Adipocytokine signaling pathway439.10 × 10−063.01 × 10−031483
Type II diabetes mellitus342.42 × 10−036.12 × 10−03861
Type I diabetes mellitus272.60 × 10−081.51 × 10−041482
Regulation of lipolysis in adipocytes274.09 × 10−041.30 × 10−041093
Ovarian steroidogenesis221.94 × 10−072.00 × 10−031471
Aldosterone-regulated sodium reabsorption204.12 × 10−024.78 × 10−02540
Antifolate resistance201.96 × 10−045.43 × 10−03854
Renin secretion161.56 × 10−021.91 × 10−058112
Tyrosine metabolism49.19 × 10−031.65 × 10−03664
Not includedGABA biosynthesis, eukaryotes, putrescine = > GABANA1.33 × 10−025.02 × 10−03332

Pathways commonly enriched in ADgset and MDDgset.

Construction of the Interaction Network of Pathway-Pathway-Genes

Among the 102 shared pathways, 101 pathways, except the gamma-aminobutyric acid (GABA) biosynthesis, are included in the cross-talk analysis. Furthermore, we added the shared genes into the pathway-pathway interaction network to form the interaction network of pathway-pathway-shared genes to explore the biological function in which shared genes participated. Of the 77 shared genes, 56 were included in 99 shared pathways (except aldosterone-regulated sodium reabsorption, GABA biosynthesis and Th1 and Th2 cell differentiation). This network consisted of 2,767 edges (connections) and 157 nodes (pathways and shared genes) (Figure 2). Of the edges, 2,380 were connections between pathways, and the other 387 were connections between pathways and the shared genes. Among the 157 nodes, there were 101 pathways (Figure 2, triangular node) and 56 shared genes (Figure 2, circular nodes).

Figure 2

Figure 2

Pathway-pathway-shared genes network of 101 shared pathways and 56 shared genes of Alzheimer's disease (AD) and major depressive disorder (MDD). Gray triangular nodes represent pathways, red circular nodes represent genes located within two or more modules, and orange circular nodes represent genes located within one module. Edges represent connections between pathway-pathway and pathway-shared genes.

Based on the biological function, this network could be largely divided into four modules, i.e., immune, neural, cancer, and endocrine (Figure 2). The first module was mainly related to immunological regulation, immunological disorders, and pathogens infection (e.g., antigen processing and presentation, T-cell receptor signaling pathway and Th17 cell differentiation, tumor necrotic factor (TNF) signaling pathway, inflammatory bowel disease, influenza A, rheumatoid arthritis, Epstein-Barr virus infection, and herpes simplex infection). The second module primarily included pathways related to neuronal function, substance addiction, and neurological disorders (e.g., cholinergic synapse, dopaminergic synapse, serotonergic synapse, amphetamine addiction, cocaine addiction, Alzheimer's disease, and Huntington's disease). The major theme of the third module was cell growth and tumorigenesis (e.g., MAPK signaling pathway, PI3K-Akt signaling pathway, bladder cancer, breast cancer, pancreatic cancer, and prostate cancer). The last module was largely concentrated in the regulation of the endocrine system, like endocrine resistance and renin secretion. Obviously, these four modules were not independent of each other; instead, they were connected through one or more pathways, indicating the existence of an immune-neuronal-cell development-endocrine regulatory network underlying the mechanisms of AD and MDD. Meanwhile, some of the 56 shared genes mainly participated in the regulation within one module (Figure 2, orange circular nodes), but others participated in the regulation of two or more modules (Figure 2, red circular nodes). The shared genes between two or more modules may play essential roles in maintaining the balance of different systems.

To further explore the connections between disease-specific pathways and their linked genes, we constructed a network for AD-specific pathways, MDD-specific pathways, and their linked genes. Among the 95 disease-specific pathways (51 AD-specific and 44 MDD-specific pathways), 68 (32 AD-specific and 36 MDD-specific pathways) met the criterion of pathway cross talk. The constructed network contained 1,371 edges and 410 nodes (pathways and genes) (Figure 3). Of the edges, 386 were connections between pathways and 985 were connections between pathways and genes. The 410 nodes included 32 AD-specific pathways (Figure 3, blue triangular node), 36 MDD-specific pathways (Figure 3, gray triangular node), 186 AD-associated genes (Figure 3, red circular nodes), 117 MDD-associated genes (Figure 3, green circular nodes), and 39 shared genes (Figure 3, orange circular nodes). This network could be divided into two modules, namely, the upper left module that included the MDD pathways and their connected genes, and the lower right module that included the AD pathways and their connected genes. The two modules were linked by several shared genes.

Figure 3

Figure 3

Pathway-pathway-genes network of disease-specific pathways and their linked genes. Blue triangular nodes represent AD-specific pathways, gray triangular nodes represent MDD-specific pathways, red circular nodes represent AD-specific genes, green circular nodes represent MDD-specific genes, and orange circular node represents shared genes. Edges represent connections between pathway-pathway and pathway-genes.

Analysis of New Candidate Genes Related to AD and MDD

We extracted the overlapping part of the AD-specific network and the MDD-specific network, which contained 3,001 nodes and 4,900 interactions. Among the nodes, 192 were genes from ADgset, 106 were genes from MDDgset, and 74 were shared genes. There were 2,629 first neighbor nodes of the 77 shared genes, from which we identified the novel candidate genes via the “guilt-by-association” principle (Oliver, 2000). That is, a node tends to participate in the same or similar cellular functions if the majority of its neighbors in the interactome network associated with specific cellular functions (e.g., a certain disease or phenotype). Ranking in descending order of node degree, we found 37 novel candidate genes (Table 3) directly interacted with 6 or more of the 77 shared genes. Furthermore, we selected 7 novel candidate genes directly interacted with 9 or more of the 77 shared genes, including ELAV like RNA-binding protein 1 (ELAVL1), heat shock protein 90 alpha family class B member 1 (HSP90AB1), heat shock protein family A (HSPA8), MDM2 proto-oncogene (MDM2), SRC proto-oncogene, non-receptor tyrosine kinase (SRC), COP9 signalosome subunit 5 (COPS5), and HDAC1 histone deacetylase 1 (HDAC1). By mapping these seven genes with the most neighbors to the disease network, we could retrieve their specific interaction network (Figure 4), which consisted of 306 nodes (genes) and 522 edges (connections). Besides the seven novel candidate genes (Figure 4, triangular nodes), this network also included 168 AD-associated genes (Figure 4, red circular nodes) and 98 MDD-associated genes (Figure 4, green circular nodes) and 33 shared genes (Figure 4, orange circular nodes). These seven novel candidate genes were closely related to AD-associated genes, MDD-associated genes, and shared genes, indicating the credibility of our results. Above all, they were expected to provide guidance for clinical and basic medical researches.

Table 3

Gene IDGene symbolNode degreeGene name
1994ELAVL112ELAV like RNA binding protein 1
3326HSP90AB110Heat shock protein 90 alpha family class B member 1
3312HSPA810Heat shock protein family A
4193MDM210MDM2 proto-oncogene
6714SRC10SRC proto-oncogene, non-receptor tyrosine kinase
10987COPS59COP9 signalosome subunit 5
3065HDAC19Histone deacetylase 1
3725JUN8Jun proto-oncogene, AP-1 transcription factor subunit
60ACTB8Actin beta
2033EP3008E1A binding protein p300
2534FYN8FYN proto-oncogene, Src family tyrosine kinase
3066HDAC28Histone deacetylase 2
5578PRKCA8Protein kinase C alpha
5970RELA8RELA proto-oncogene, NF-kB subunit
4088SMAD38SMAD family member 3
4904YBX18Y-box binding protein 1
1051CEBPB7CCAAT enhancer binding protein beta
8841HDAC37Histone deacetylase 3
4609MYC7MYC proto-oncogene, bHLH transcription factor
6774STAT37Signal transducer and activator of transcription 3
203068TUBB7Tubulin beta class I
1869E2F16E2F transcription factor 1
2885GRB26Growth factor receptor bound protein 2
10524KAT56Lysine acetyltransferase 5
4067LYN6LYN proto-oncogene, Src family tyrosine kinase
4130MAP1A6Microtubule associated protein 1A
5594MAPK16Mitogen-activated protein kinase 1
4869NPM16Nucleophosmin 1
5499PPP1CA6Protein phosphatase 1 catalytic subunit alpha
5501PPP1CC6Protein phosphatase 1 catalytic subunit gamma
5566PRKACA6Protein kinase cAMP-activated catalytic subunit alpha
5925RB16RB transcriptional corepressor 1
6464SHC16SHC adaptor protein 1
4087SMAD26SMAD family member 2
10273STUB16STIP1 homology and U-box containing protein 1
6613SUMO26Small ubiquitin like modifier 2
7534YWHAZ6Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta

Novel candidate genes potentially related to AD and MDD.

Figure 4

Figure 4

Seven predicted novel candidate genes and known genes connected with them. Triangular nodes are predicted novel candidate genes, red nodes are AD-associated genes, green nodes are MDD-associated genes, and orange nodes are shared genes.

We further performed a spatio-temporal and cell type-specific expression analysis on the 37 new candidate genes. Of which, 28 genes could be detected in public available brain cell type expression dataset CSEA, and 6 of them (i.e., E2F1, RELA, MAPK1, FYN, CEBPB, STAT3) showed substantial overrepresentation (pSI = 0.05, Supplementary Table S6) in cone cells, oligodendrocytes, interneurons, and others (Figure 5A). In addition, all the 37 genes expressed in the brain regions of different development stages, and 20 genes (i.e., ACTB, CEBPB, COPS5, E2F1, ELAVL1, HDAC2, JUN, LYN, MAP1A, MAPK1, MYC, NPM1, PPP1CC, PRKCA, RB1, SMAD3, SRC, SUMO2, TUBB, YBX1) showed substantial overrepresentation under the threshold pSI = 0.05 (Supplementary Table S7). Striatum, hippocampus, cortex, and thalamus contained 10, 8, 8, and 8 novel candidate genes under pSI = 0.05, respectively (Figure 5B). Further, the number of genes was generally fewer after birth (Figure 5C). In addition, we searched the expression patterns of 37 novel candidate genes in Human Brain Transcriptome (HBT) database (www.humanbraintranscriptome.org). We found that the expression of most genes changed with different developmental stages, which are specifically manifested in a downward trend after birth ((Supplementary Figures 1A–D) and an upward trend after birth (Supplementary Figures 1E–H). Thus, most of these novel candidate genes were significantly expressed in the brain or brain cell lines and changed with different ages, indicating they could play important roles in the normal function of the brain.

Figure 5

Figure 5

Spatio-temporal expression of the 37 novel candidate genes. (A) Number of the novel candidate genes with substantial overrepresentation in brain cell types. (B) Number of the novel candidate genes with substantial overrepresentation in brain regions. (C) Number of the novel candidate genes expressed in specific developmental period.

Discussion

In the past decades, our capability in identifying genes related to complex disease like AD and MDD has been greatly improved. However, complex diseases are rarely a straight-forward dysfunction in a specific gene, but rather a consequence of the collective interplay of multiple genes (Menche et al., 2015). While previous studies have reported the neurobiology and susceptibility gene correlations between MDD and neurodegenerative diseases like AD, the underlying mechanisms still need to be explored. The analysis on the susceptibility gene of AD and MDD from the aspect of systems biology could provide insight into the comorbidity of these two diseases beyond studies based on single risk gene. In current research, we analyzed the functions of genes associated with AD and MDD and deciphered the interrelationship between these genes. In addition, we predicted the novel candidate genes potentially related to AD and MDD.

Based on 650 AD-associated genes and 447 MDD-associated genes, we retrieved 77 shared genes between AD and MDD. Then, GO biological processes enriched in the 77 shared genes and pathways enriched in AD-associated genes or MDD-associated genes were detected. Biological processes such as response to substances (e.g., toxic substances, hormones, oxidative stress, and alcohol), cognition, regulation of neurogenesis, aging, and regulation of developmental growth were enriched in the shared genes, indicating the significance of these processes in AD and MDD pathology. Among the 102 overlapping pathways, AD was the most significant one, indicating some AD susceptibility genes are MDD genes. Other pathways were related to neuroactive ligand-receptor interaction, synapse (e.g., serotonergic, dopaminergic, and cholinergic), and substance addictions (e.g., amphetamine and cocaine). Neuroactive ligand-receptor interaction, as a neurodevelopment-related pathway, has been implicated in AD and MDD. Synaptic dysfunction can lead to two important symptoms of AD and MDD, i.e., cognitive impairment and memory loss, indicating the importance of synapse-related pathways in pathophysiology of these two disorders (Boyle et al., 2010; Sierksma et al., 2010). Consistently, several shared pathways were related to immune responses (e.g., Th1 and Th2 cell differentiation and Th17 cell differentiation) and immune diseases (e.g., inflammatory bowel disease and type I diabetes). For example, Th17 participated in neuroinflammation and neurodegeneration of AD by regulating inflammatory cytokines and signaling (Zhang et al., 2013). In addition, a predominant Th1/Th17 inflammatory response in placental immune cells plays a crucial role in regulating depressive symptoms during pregnancy (Leff-Gelman et al., 2016). It has been showed the involvement of neuroinflammation in AD and MDD pathological mechanisms, in which the complement molecules exert indispensable efforts (Lian et al., 2016; Yao and Li, 2020). Studies have shown that the microbial imbalance of the intestinal flora may explain the association between inflammatory bowel disease and cognitive function of AD or MDD (Bonaz and Bernstein, 2013; Chen et al., 2016; Gareau, 2016; Spielman et al., 2018). Furthermore, it has been indicated that diabetes can cause deterioration of cognitive function and increase the risk of dementia and depression (Szatmári et al., 2017). Moreover, infection-related pathways were also enriched in our disease-related genes, such as Epstein-Barr virus (EBV) infection and herpes simplex virus (HSV) infection. Previous studies have reported HSV-1 infection may impair the hippocampal function and is a negative factor of AD (Li Puma et al., 2019). EBV has also been reported to be related to the cognitive decline in patients with AD (Shim et al., 2017). In addition, patients with MDD have altered levels and patterns of EBV antibodies indicating a therapeutic intervention targeting EBV available for MDD individuals (Jones-Brando et al., 2020). In our study, virus infection like HSV and EBV infection included in the shared pathway list, indicating that virus infection may increase the risk of AD and MDD. In the future, we should explore the impact of viral infections on the clinical treatment of AD and MDD. Moreover, some pathways involved in cancers (e.g., prostate cancer, glioma, and pancreatic cancer) and their developmental processes (e.g., apoptosis) were also existed in 102 shared pathways. Researchers have revealed that multitudinous cancer survivors had lower AD and MDD risk (Caruso et al., 2017; Frain et al., 2017). Our shared pathways include some cancer and related pathways, suggesting that cancer and/or cancer treatment may affect AD and MDD progression. In the pathway cross-talk analysis, the pathway-pathway-shared genes interaction network containing 2,767 edges (connections) and 157 nodes (pathways and shared genes). Pathways are tightly connected in this network and may collectively influence the occurrence of AD and MDD. In sharp contrast with the network between shared pathways, in the network of disease-specific pathways, the pathways between AD and MDD were only sparsely connected through a few shared genes.

We further identified 37 new genes highly interconnected with the 77 shared genes, which were also potentially related to AD and MDD. Of which, the node degrees of ELAVL1, HSP90AB1, HSPA8, MDM2, SRC, COPS5, and HDAC1 were the highest with degrees >8.0. ELAVL1 (also known as HuR/HuA), an RNA-binding protein, belongs to the ELAV/Hu family and is ubiquitously expressed in all human tissues (Fan and Steitz, 1998). ELAVL1 plays a critical role in stress-induced synaptic dysfunction (He et al., 2019) and contributes to neuroinflammation in ELAVL1-knockout mouse model (Chen et al., 2017). It has been showed that neuroinflammation and stress are involved in AD and MDD pathological mechanisms (Bolos et al., 2017; Ignácio et al., 2019; Sotiropoulos et al., 2019). ELAVL1 may promote the progress of AD and MDD through regulating neuroinflammatory and stress-mediated synaptic function. Whether it can be used as a potential therapeutic target for AD, MDD, and their comorbidity in the future remains to be further explored. HSP90AB1 and HSPA8 are heat shock proteins (HSPs), and disruption of HSPs expression has been described as a possible mechanism in the etiology and progression of AD (Gorenberg and Chandra, 2017). HSPA8 is a member of HSP70 family proteins that may affect the action of antidepressants and thus their therapeutic efficacy in treatment of MDD (Pae et al., 2007). Our results suggest that HSPs are expected to be potential intervention targets for the treatment of the comorbidity underlying AD and MDD. MDM2 gene is a proto-oncogene, and its expression mediates important cancer pathology (Oliner et al., 2016). Our pathway cross-talk analysis showed the close interlink between the neural and cancer models. The involvement of MDM2 in comorbidity of AD and MDD may be related to tumorigenesis. COPS5 is an evolutionarily conserved and multifunctional protein and its overexpression in mice brain significantly increased amyloid-β protein levels in the cortex and hippocampus (Wang et al., 2015). As we all know, amyloid-β takes part in the pathogenesis of both AD and MDD (Blennow et al., 2006; Pomara et al., 2012). The activation of HDAC1, a histone deacetylase 1, has shown therapeutic potential against functional decline in brain aging and neurodegeneration (Pao et al., 2020). It has been reported that HDAC1 is involved in fluoxetine-mediated antidepressant mechanism in the LPS-induced depression mouse model (Li et al., 2021). Although there are few reports that indicate its role of HDAC1 in the comorbidity processes of AD and MDD, histone modifications are widely involved in the gene expression regulation. Our PPI results show that HDAC1 is closely interlinked with shared genes of AD and MDD, implying that HDAC1 has potential regulatory functions on the shared genes. To verify the results, we further checked the expression pattern of these 37 novel candidate genes in brain tissues and cell lines. As expected, most of the novel candidate genes were significantly expressed in the brain.

There were several limitations in this study. For instance, our disease susceptibility genes were mainly collected from currently available studies and public databases, which should be deemed as incomplete as our exploration on the molecular features of AD and MDD is still ongoing. In future, more genes involved in these disorders may be identified, and some false-positive genes may be excluded. Similarly, the protein-protein interaction databases and the human interactome constructed in this work should be treated in the same way.

Conclusion

Based on the available susceptibility genes related to AD or MDD, we analyzed the biochemical pathways and the interactions between genes and pathways potentially involved in the comorbidity of the two disorders. We found that biological processes related to metabolism, immune response, cancer, and neural disease were enriched in 77 shared genes. And there existed four interconnected modules between 102 shared pathways, namely, immune, neuro, cancer, and endocrine. Of significance, we identified 37 interesting novel candidate genes related to AD and MDD based on the PPI analysis. Such analysis of AD- and MDD-related genes will not only improve our understanding of the mechanisms or their comorbidity but also help us to explore new biomarkers for the comorbidity of these diseases.

Funding

This study was supported in part by grants from National Key Research and Development Program of China (No. 2016YFC0906300) and National Natural Science Foundation of China (Nos. 91746205 and 31271411).

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.

Statements

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

PG and SC designed the experiments. PG, HW, and JW performed the experiments and data analysis. PG, YW, and JW wrote the manuscript. All authors have read and approved the final manuscript.

Acknowledgments

The authors would like to thank Yanshi Hu and Ting Fan for the help with the data collection in previous study.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2021.789698/full#supplementary-material

Supplementary Figure 1

The expression pattern of some novel candidate genes in specific developmental period. Some genes are specifically manifested in a downward trend after birth (A–D) and an upward trend after birth (E–H). The y-axis represents the intensity of gene expression, and x-axis represents age. Periods 1–7 represent postconceptional age and 8–15 represent postnatal years. NCX, neocortex; HIP, hippocampus; AMY, amygdala; ST, striatum; MD, mediodorsal nucleus of the thalamus; CBC, cerebellar cortex.

Supplementary Material 1

The background PPIN constructed in the current study.

Supplementary Table S1

List of 650 genes associated with Alzheimer's Disease.

Supplementary Table S2

List of 447 genes associated with major depressive disorder.

Supplementary Table S3

The GO biological processes significantly enriched in the 77 shared genes between AD and MDD.

Supplementary Table S4

153 pathways significantly enriched in the 650 AD-associated genes.

Supplementary Table S5

146 pathways significantly enriched in the 447 MDD-associated genes.

Supplementary Table S6

Expression of novel candidate genes in brain cells.

Supplementary Table S7

Expression of novel candidate genes in brain region and developmental stages.

References

  • 1

    Aoki-KinoshitaK. F.KanehisaM. (2007). Gene annotation and pathway mapping in KEGG. Methods Mol. Biol. 396, 7191. 10.1007/978-1-59745-515-2_6

  • 2

    AssenovY.RamírezF.SchelhornS. E.LengauerT.AlbrechtM. (2008). Computing topological parameters of biological networks. Bioinformatics24, 282284. 10.1093/bioinformatics/btm554

  • 3

    BallardC.BannisterC.SolisM.OyebodeF.WilcockG. (1996). The prevalence, associations and symptoms of depression amongst dementia sufferers. J. Affect. Disord. 36, 135144. 10.1016/0165-0327(95)00072-0

  • 4

    BarberR.. (2011). Inflammatory signaling in Alzheimer disease and depression. Cleve. Clin. J. Med. 78(Suppl. 1), S4749. 10.3949/ccjm.78.s1.08

  • 5

    BlennowK.de LeonM. J.ZetterbergH. (2006). Alzheimer's disease. Lancet368, 387403. 10.1016/S0140-6736(06)69113-7

  • 6

    BolosM.PereaJ. R.AvilaJ. (2017). Alzheimer's disease as an inflammatory disease. Biomol. Concepts8, 3743. 10.1515/bmc-2016-0029

  • 7

    BonazB. L.BernsteinC. N. (2013). Brain-gut interactions in inflammatory bowel disease. Gastroenterology144, 3649. 10.1053/j.gastro.2012.10.003

  • 8

    BoyleL. L.PorsteinssonA. P.CuiX.KingD. A.LynessJ. M. (2010). Depression predicts cognitive disorders in older primary care patients. J. Clin. Psychiatry71, 7479. 10.4088/JCP.08m04724gry

  • 9

    CaiM.LeeJ. H.YangE. J. (2019). Electroacupuncture attenuates cognition impairment via anti-neuroinflammation in an Alzheimer's disease animal model. J. Neuroinflam.16:264. 10.1186/s12974-019-1665-3

  • 10

    CarusoR.NanniM. G.RibaM.SabatoS.MitchellA. J.CroceE.et al. (2017). Depressive spectrum disorders in cancer: prevalence, risk factors and screening for depression: a critical review. Acta Oncol. 56, 146155. 10.1080/0284186X.2016.1266090

  • 11

    ChenC. H.LinC. L.KaoC. H. (2016). Irritable bowel syndrome is associated with an increased risk of dementia: a nationwide population-based study. PLoS ONE11:e0144589. 10.1371/journal.pone.0144589

  • 12

    ChenJ.BardesE. E.AronowB. J.JeggaA. G. (2009). ToppGene suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 37, W305311. 10.1093/nar/gkp427

  • 13

    ChenJ.MartindaleJ. L.CramerC.GorospeM.AtasoyU.DrewP. D.et al. (2017). The RNA-binding protein HuR contributes to neuroinflammation by promoting C-C chemokine receptor 6 (CCR6) expression on Th17 cells. J. Biol. Chem. 292, 1453214543. 10.1074/jbc.M117.782771

  • 14

    CowleyM. J.PineseM.KassahnK. S.WaddellN.PearsonJ. V.GrimmondS. M.et al. (2012). PINA v2.0: mining interactome modules. Nucleic Acids Res. 40, D862865. 10.1093/nar/gkr967

  • 15

    CroftD.O'KellyG.WuG.HawR.GillespieM.MatthewsL.et al. (2011). Reactome: a database of reactions, pathways and biological processes. Nucleic Acids Res39, D691697. 10.1093/nar/gkq1018

  • 16

    De WinterF. L.EmsellL.BouckaertF.ClaesL.JainS.FarrarG.et al. (2017). No association of lower hippocampal volume with Alzheimer's disease pathology in late-life depression. Am. J. Psychiatry174, 237245. 10.1176/appi.ajp.2016.16030319

  • 17

    EnacheD.WinbladB.AarslandD. (2011). Depression in dementia: epidemiology, mechanisms, and treatment. Curr. Opin. Psychiatry24, 461472. 10.1097/YCO.0b013e32834bb9d4

  • 18

    FanT.HuY.XinJ.ZhaoM.WangJ. (2020). Analyzing the genes and pathways related to major depressive disorder via a systems biology approach. Brain Behav.10:e01502. 10.1002/brb3.1502

  • 19

    FanX. C.SteitzJ. A. (1998). HNS, a nuclear-cytoplasmic shuttling sequence in HuR. Proc. Natl. Acad. Sci. U.S.A. 95, 1529315298. 10.1073/pnas.95.26.15293

  • 20

    FernándezM.GobarttA. L.Bala,ñáM. (2010). Behavioural symptoms in patients with Alzheimer's disease and their association with cognitive impairment. BMC Neurol. 10:87. 10.1186/1471-2377-10-87

  • 21

    FerrariR.KiaD. A.TomkinsJ. E.HardyJ.WoodN. W.LoveringR. C.et al. (2018). Stratification of candidate genes for Parkinson's disease using weighted protein-protein interaction network analysis. BMC Genomics19:452. 10.1186/s12864-018-4804-9

  • 22

    FrainL.SwansonD.ChoK.GagnonD.LuK. P.BetenskyR. A.et al. (2017). Association of cancer and Alzheimer's disease risk in a national cohort of veterans. Alzheimers Dement.13, 13641370. 10.1016/j.jalz.2017.04.012

  • 23

    GareauM. G.. (2016). Cognitive function and the microbiome. Int. Rev. Neurobiol. 131, 227246. 10.1016/bs.irn.2016.08.001

  • 24

    GorenbergE. L.ChandraS. S. (2017). The role of co-chaperones in synaptic proteostasis and neurodegenerative disease. Front. Neurosci. 11:248. 10.3389/fnins.2017.00248

  • 25

    GuoP.CaoC.MaY.WangJ. (2021). Network-based analysis on genetic variants reveals the immunological mechanism underlying Alzheimer's disease. J. Neural. Transm. 128, 803816. 10.1007/s00702-021-02337-9

  • 26

    HeZ. X.SongH. F.LiuT. Y.MaJ.XingZ. K.YinY. Y.et al. (2019). HuR in the medial prefrontal cortex is critical for stress-induced synaptic dysfunction and depressive-like symptoms in mice. Cereb. Cortex29, 27372747. 10.1093/cercor/bhz036

  • 27

    HuY. S.XinJ.HuY.ZhangL.WangJ. (2017). Analyzing the genes related to Alzheimer's disease via a network and pathway-based approach. Alzheimers Res. Ther. 9:29. 10.1186/s13195-017-0252-z

  • 28

    IgnácioZ. M.da SilvaR. S.PlissariM. E.QuevedoJ.RéusG. Z. (2019). Physical exercise and neuroinflammation in major depressive disorder. Mol. Neurobiol. 56, 83238335. 10.1007/s12035-019-01670-1

  • 29

    JackC. R.Jr.PetersenR. C.XuY. C.O'BrienP. C.SmithG. E.IvnikR. J.et al. (1999). Prediction of AD with MRI-based hippocampal volume in mild cognitive impairment. Neurology52, 13971403. 10.1212/WNL.52.7.1397

  • 30

    JiaP.KaoC.-F.KuoP.-H.ZhaoZ. (2011). A comprehensive network and pathway analysis of candidate genes in major depressive disorder. BMC Syst. Biol. 5:S12. 10.1186/1752-0509-5-S3-S12

  • 31

    Jones-BrandoL.DickersonF.FordG.StallingsC.OrigoniA.KatsafanasE.et al. (2020). Atypical immune response to Epstein-Barr virus in major depressive disorder. J. Affect. Disord. 264, 221226. 10.1016/j.jad.2019.11.150

  • 32

    KalesH. C.GitlinL. N.LyketsosC. G. (2015). Assessment and management of behavioral and psychological symptoms of dementia. BMJ350:h369. 10.1136/bmj.h369

  • 33

    KimD. H.PayneM. E.LevyR. M.MacFallJ. R.SteffensD. C. (2002). APOE genotype and hippocampal volume change in geriatric depression. Biol. Psychiatry51, 426429. 10.1016/S0006-3223(01)01272-0

  • 34

    KohlM.WieseS.WarscheidB. (2011). Cytoscape: software for visualization and analysis of biological networks. Methods Mol. Biol. 696, 291303. 10.1007/978-1-60761-987-1_18

  • 35

    KrishnanK. R.TuplerL. A.RitchieJ. C.Jr.McDonaldW. M.KnightD. L.NemeroffC. B.et al. (1996). Apolipoprotein E-epsilon 4 frequency in geriatric depression. Biol. Psychiatry40, 6971. 10.1016/0006-3223(95)00424-6

  • 36

    LeeH. B.LyketsosC. G. (2003). Depression in Alzheimer's disease: heterogeneity and related issues. Biol. Psychiatry54, 353362. 10.1016/S0006-3223(03)00543-2

  • 37

    Leff-GelmanP.Mancilla-HerreraI.Flores-RamosM.Cruz-FuentesC.Reyes-GrajedaJ. P.García-Cuétara MdelP.et al. (2016). The Immune System And The Role Of Inflammation In Perinatal Depression. Neurosci. Bull. 32, 398420. 10.1007/s12264-016-0048-3

  • 38

    Li PumaD. D.PiacentiniR.LeoneL.GironiK.MarcocciM. E.De ChiaraG.et al. (2019). Herpes simplex virus type-1 infection impairs adult hippocampal neurogenesis via Amyloid-beta protein accumulation. Stem Cells37, 14671480. 10.1002/stem.3072

  • 39

    LiW.AliT.ZhengC.LiuZ.HeK.ShahF. A.et al. (2021). Fluoxetine regulates eEF2 activity (phosphorylation) via HDAC1 inhibitory mechanism in an LPS-induced mouse model of depression. J. Neuroinflam.18:38. 10.1186/s12974-021-02091-5

  • 40

    LianH.LitvinchukA.ChiangA. C.AithmittiN.JankowskyJ. L.ZhengH. (2016). Astrocyte-microglia cross talk through complement activation modulates amyloid pathology in mouse models of Alzheimer's disease. J. Neurosci. 36, 577589. 10.1523/JNEUROSCI.2117-15.2016

  • 41

    LindsayJ.LaurinD.VerreaultR.HébertR.HelliwellB.HillG. B.et al. (2002). Risk factors for Alzheimer's disease: a prospective analysis from the Canadian Study of Health and Aging. Am. J. Epidemiol. 156, 445453. 10.1093/aje/kwf074

  • 42

    LiuM.FanR.LiuX.ChengF.WangJ. (2015). Pathways and networks-based analysis of candidate genes associated with nicotine addiction. PLoS ONE10:e0127438. 10.1145/2818302

  • 43

    LyketsosC. G.BakerL.WarrenA.SteeleC.BrandtJ.SteinbergM.et al. (1997a). Depression, delusions, and hallucinations in Alzheimer's disease: no relationship to apolipoprotein E genotype. J. Neuropsychiatry Clin. Neurosci. 9,6467. 10.1176/jnp.9.1.64

  • 44

    LyketsosC. G.SteeleC.BakerL.GalikE.KopunekS.SteinbergM.et al. (1997b). Major and minor depression in Alzheimer's disease: prevalence and impact. J. Neuropsychiatry Clin. Neurosci. 9,556561.

  • 45

    MencheJ.SharmaA.KitsakM.GhiassianS. D.VidalM.LoscalzoJ.et al. (2015). Disease networks. uncovering disease-disease relationships through the incomplete interactome. Science347:1257601. 10.1126/science.1257601

  • 46

    MigliorelliR.TesónA.SabeL.PetracchiM.LeiguardaR.StarksteinS. E. (1995). Prevalence and correlates of dysthymia and major depression among patients with Alzheimer's disease. Am. J. Psychiatry152, 3744. 10.1176/ajp.152.1.37

  • 47

    OlinJ. T.SchneiderL. S.KatzI. R.MeyersB. S.AlexopoulosG. S.BreitnerJ. C.et al. (2002). Provisional diagnostic criteria for depression of Alzheimer disease. Am. J. Geriatr. Psychiatry10, 125128. 10.1097/00019442-200203000-00003

  • 48

    OlinerJ. D.SaikiA. Y.CaenepeelS. (2016). The role of MDM2 amplification and overexpression in tumorigenesis. Cold Spring Harb. Perspect. Med. 6:a026336. 10.1101/cshperspect.a026336

  • 49

    OliverS.. (2000). Guilt-by-association goes global. Nature403, 601603. 10.1038/35001165

  • 50

    PaeC. U.MandelliL.SerrettiA.PatkarA. A.KimJ. J.LeeC. U.et al. (2007). Heat-shock protein-70 genes and response to antidepressants in major depression. Prog. Neuropsychopharmacol. Biol. Psychiatry31, 10061011. 10.1016/j.pnpbp.2007.02.011

  • 51

    PaoP. C.PatnaikD.WatsonL. A.GaoF.PanL.WangJ.et al. (2020). HDAC1 modulates OGG1-initiated oxidative DNA damage repair in the aging brain and Alzheimer's disease. Nat. Commun.11:2484. 10.1038/s41467-020-16361-y

  • 52

    PassamontiL.TsvetanovK. A.JonesP. S.Bevan-JonesW. R.ArnoldR.BorchertR. J.et al. (2019). Neuroinflammation and functional connectivity in Alzheimer's disease: interactive influences on cognitive performance. J. Neurosci. 39, 72187226. 10.1523/JNEUROSCI.2574-18.2019

  • 53

    PiñeroJ.Ramírez-AnguitaJ. M.Saüch-PitarchJ.RonzanoF.CentenoE.SanzF.et al. (2020). The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 48, D845d855. 10.1093/nar/gkz1021

  • 54

    PomaraN.BrunoD.SarrealA. S.HernandoR. T.NierenbergJ.PetkovaE.et al. (2012). Lower CSF amyloid beta peptides and higher F2-isoprostanes in cognitively intact elderly individuals with major depressive disorder. Am. J. Psychiatry169, 523530. 10.1176/appi.ajp.2011.11081153

  • 55

    PrinceM.AliG. C.GuerchetM.PrinaA. M.AlbaneseE.WuY. T. (2016). Recent global trends in the prevalence and incidence of dementia, and survival with dementia. Alzheimers Res. Ther. 8:23. 10.1186/s13195-016-0188-8

  • 56

    SaczynskiJ. S.BeiserA.SeshadriS.AuerbachS.WolfP. A.AuR. (2010). Depressive symptoms and risk of dementia: the Framingham Heart Study. Neurology75, 3541. 10.1212/WNL.0b013e3181e62138

  • 57

    ScarmeasN.BrandtJ.AlbertM.DevanandD. P.MarderK.BellK.et al. (2002). Association between the APOE genotype and psychopathologic symptoms in Alzheimer's disease. Neurology58, 11821188. 10.1212/WNL.58.8.1182

  • 58

    SchaeferC. F.AnthonyK.KrupaS.BuchoffJ.DayM.HannayT.et al. (2009). PID: the pathway interaction database. Nucleic Acids Res. 37, D674679. 10.1093/nar/gkn653

  • 59

    SchmaalL.HibarD. P.SämannP. G.HallG. B.BauneB. T.JahanshadN.et al. (2017). Cortical abnormalities in adults and adolescents with major depression based on brain scans from 20 cohorts worldwide in the ENIGMA Major Depressive Disorder Working Group. Mol. Psychiatry22, 900909. 10.1038/mp.2016.60

  • 60

    SchuffN.WoernerN.BoretaL.KornfieldT.ShawL. M.TrojanowskiJ. Q.et al. (2009). MRI of hippocampal volume loss in early Alzheimer's disease in relation to ApoE genotype and biomarkers. Brain132, 10671077. 10.1093/brain/awp007

  • 61

    Serrano-PozoA.DasS.HymanB. T. (2021). APOE and Alzheimer's disease: advances in genetics, pathophysiology, and therapeutic approaches. Lancet Neurol. 20, 6880. 10.1016/S1474-4422(20)30412-9

  • 62

    ShimS. M.CheonH. S.JoC.KohY. H.SongJ.JeonJ. P. (2017). Elevated epstein-barr virus antibody level is associated with cognitive decline in the Korean elderly. J. Alzheimers Dis. 55, 293301. 10.3233/JAD-160563

  • 63

    SierksmaA. S.van den HoveD. L.SteinbuschH. W.PrickaertsJ. (2010). Major depression, cognitive dysfunction and Alzheimer's disease: is there a link?Eur. J. Pharmacol. 626, 7282. 10.1016/j.ejphar.2009.10.021

  • 64

    SotiropoulosI.SilvaJ. M.GomesP.SousaN.AlmeidaO. F. X. (2019). stress and the etiopathogenesis of alzheimer's disease and depression. Adv. Exp. Med. Biol. 1184, 241257. 10.1007/978-981-32-9358-8_20

  • 65

    SpielmanL. J.GibsonD. L.KlegerisA. (2018). Unhealthy gut, unhealthy brain: the role of the intestinal microbiota in neurodegenerative diseases. Neurochem. Int. 120, 149163. 10.1016/j.neuint.2018.08.005

  • 66

    StarksteinS. E.JorgeR.MizrahiR.RobinsonR. G. (2005). The construct of minor and major depression in Alzheimer's disease. Am. J. Psychiatry162, 20862093. 10.1176/appi.ajp.162.11.2086

  • 67

    SteinbergS.StefanssonH.JonssonT.JohannsdottirH.IngasonA.HelgasonH.et al. (2015). Loss-of-function variants in ABCA7 confer risk of Alzheimer's disease. Nat. Genet. 47, 445447. 10.1038/ng.3246

  • 68

    SweetR. A.HamiltonR. L.ButtersM. A.MulsantB. H.PollockB. G.LewisD. A.et al. (2004). Neuropathologic correlates of late-onset major depression. Neuropsychopharmacology29, 22422250. 10.1038/sj.npp.1300554

  • 69

    SzatmáriS.Orbán-KisK.MihályI.LázárA. S. (2017). [Diabetes, dementia, depression, distress]. Ideggyogy Sz70, 293299. 10.18071/isz.70.0293

  • 70

    TaylorW. D.McQuoidD. R.PayneM. E.ZannasA. S.MacFallJ. R.SteffensD. C. (2014). Hippocampus atrophy and the longitudinal course of late-life depression. Am. J. Geriatr. Psychiatry22, 15041512. 10.1016/j.jagp.2013.11.004

  • 71

    The Alzheimer's Association (2020). 2020 Alzheimer's disease facts and figures. Alzheimer's Dement16, 391460. 10.1002/alz.12068

  • 72

    van der FlierW. M.ScheltensP. (2009). Alzheimer disease: hippocampal volume loss and Alzheimer disease progression. Nat. Rev. Neurol.5, 361362. 10.1038/nrneurol.2009.94

  • 73

    VinkersD. J.GusseklooJ.StekM. L.WestendorpR. G.van der MastR. C. (2004). Temporal relation between depression and cognitive impairment in old age: prospective population based study. BMJ329:881. 10.1136/bmj.38216.604664.DE

  • 74

    WangJ.DuncanD.ShiZ.ZhangB. (2013). WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013. Nucleic Acids Res. 41, W7783. 10.1093/nar/gkt439

  • 75

    WangR.WangH.CarreraI.XuS.LakshmanaM. K. (2015). COPS5 protein overexpression increases amyloid plaque burden, decreases spinophilin-immunoreactive puncta, and exacerbates learning and memory deficits in the mouse brain. J. Biol. Chem. 290, 92999309. 10.1074/jbc.M114.595926

  • 76

    WuK. Y.HsiaoI. T.ChenC. S.ChenC. H.HsiehC. J.WaiY. Y.et al. (2014). Increased brain amyloid deposition in patients with a lifetime history of major depression: evidenced on 18F-florbetapir (AV-45/Amyvid) positron emission tomography. Eur. J. Nucl. Med. Mol. Imaging41, 714722. 10.1007/s00259-013-2627-0

  • 77

    XuX.WellsA. B.O'BrienD. R.NehoraiA.DoughertyJ. D. (2014). Cell type-specific expression analysis to identify putative cellular mechanisms for neurogenetic disorders. J. Neurosci.34, 14201431. 10.1523/JNEUROSCI.4488-13.2014

  • 78

    YaoQ.LiY. (2020). Increased serum levels of complement C1q in major depressive disorder. J. Psychosom. Res. 133:110105. 10.1016/j.jpsychores.2020.110105

  • 79

    ZhangB.KirovS.SnoddyJ. (2005). WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. 33, W741748. 10.1093/nar/gki475

  • 80

    ZhangJ.KeK. F.LiuZ.QiuY. H.PengY. P. (2013). Th17 cell-mediated neuroinflammation is involved in neurodegeneration of aβ1-42-induced Alzheimer's disease model rats. PLoS ONE8:e75786. 10.1371/journal.pone.0075786

Summary

Keywords

Alzheimer's disease, major depressive disorder, comorbidity, pathway cross-talk, functional enrichment analysis

Citation

Guo P, Chen S, Wang H, Wang Y and Wang J (2022) A Systematic Analysis on the Genes and Their Interaction Underlying the Comorbidity of Alzheimer's Disease and Major Depressive Disorder. Front. Aging Neurosci. 13:789698. doi: 10.3389/fnagi.2021.789698

Received

05 October 2021

Accepted

20 December 2021

Published

20 January 2022

Volume

13 - 2021

Edited by

Francesco Panza, University of Bari Aldo Moro, Italy

Reviewed by

Takahisa Kanekiyo, Mayo Clinic Florida, United States; Yuyan Cheng, UCLA Health System, United States

Updates

Copyright

*Correspondence: Yaogang Wang Ju Wang

This article was submitted to Alzheimer's Disease and Related Dementias, a section of the journal Frontiers in Aging Neuroscience

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics