- Open Access
Analysis of brain region-specific co-expression networks reveals clustering of established and novel genes associated with Alzheimer disease
Alzheimer's Research & Therapy volume 12, Article number: 103 (2020)
Identifying and understanding the functional role of genetic risk factors for Alzheimer disease (AD) has been complicated by the variability of genetic influences across brain regions and confounding with age-related neurodegeneration.
A gene co-expression network was constructed using data obtained from the Allen Brain Atlas for multiple brain regions (cerebral cortex, cerebellum, and brain stem) in six individuals. Gene network analyses were seeded with 52 reproducible (i.e., established) AD (RAD) genes. Genome-wide association study summary data were integrated with the gene co-expression results and phenotypic information (i.e., memory and aging-related outcomes) from gene knockout studies in Drosophila to generate rankings for other genes that may have a role in AD.
We found that co-expression of the RAD genes is strongest in the cortical regions where neurodegeneration due to AD is most severe. There was significant evidence for two novel AD-related genes including EPS8 (FDR p = 8.77 × 10−3) and HSPA2 (FDR p = 0.245).
Our findings indicate that AD-related risk factors are potentially associated with brain region-specific effects on gene expression that can be detected using a gene network approach.
Neurodegenerative diseases, such as Alzheimer disease (AD), Parkinson disease (PD), Huntington disease (HD), and amyotrophic lateral sclerosis, impair or damage neurons. Although many sub-cellular similarities between neurodegenerative diseases have been identified , the regional differences between them are quite profound [2,3,4,5]. For example, neuronal cell death from HD is primarily localized to the basal ganglia, whereas both AD and PD result in cell death throughout the brain . Furthermore, PD causes the most severe cell death in the substantia nigra  whereas AD most heavily affects the hippocampus, the frontal cortex, and the temporal lobe . These studies highlight the importance of studying gene expression signatures and relationships of AD-associated genes in different brain regions. For instance, an increased correlation in gene expression among two AD-associated genes in the brain structures such as the cortex as compared to other brain regions suggests either a functional relationship or cell/sub-region-specific expression biases towards cell types where the disease tend to originate or progress most rapidly.
Altered functional connectivity between brain regions has been demonstrated for several neuropsychiatric diseases including schizophrenia, depression, and AD using functional magnetic resonance imaging [6,7,8]. Brain imaging and neuropathological studies indicate that the hippocampus, which has a role in memory formation, is one of the first structures showing a marked neuronal loss in AD and, compared to other regions, suffers the largest relative reduction in volume by the latter stages of the disease . Regional specificity is also evident by longitudinal patterning of the AD-related tau and amyloid-β proteins that aggregate into neurofibrillary tangles and senile plaques, respectively . In the early stages of AD, a small number of tangles typically form in the brain stem and then spread aggressively to most of the cerebrum by the latest stage [11, 12]. Amyloid plaques form in the opposite pattern, beginning primarily in the outer cortex and spreading inward and then to the brain stem . Notably, very few protein aggregates form in the cerebellum even at the most severe stages of AD.
Differences in AD severity between the regions of the brain may be a consequence of a variety of factors. One such factor is the tissue-specific expression patterns of genes throughout the body, which is a relevant consideration for the human brain given its vast complexity and compartmentalization . An additional factor may also be the changing cell type fractions observed between major regions of the brain [14, 15]. Large-scale multi-omic approaches have been able to assist in understanding these complicated roots of neuropsychiatric disease [16, 17]. Furthermore, they have been able to identify novel disease-related gene candidates [18,19,20].
In this study, we integrated network-based correlation methods with existing genome-wide association study (GWAS) data and gene expression data derived from the brain in order to identify additional AD-related genes using a network methodology. In addition to identifying several novel biologically relevant genes for AD, we show that the strength of the correlations among previously established AD genes increases when the networks are restricted to the sub-regions of the brain that are most impacted by AD.
Acquisition of GWAS data and curation of AD genes
We obtained summarized results from a GWAS for AD risk conducted using 16,175 AD cases and 17,175 controls of European ancestry that were obtained as previously described [21, 22]. Association evidence with each gene was derived from p values for the association with individual single nucleotide polymorphisms (SNPs) corrected for multiple testing using an approximation that has been shown to be a conservative adjustment for recombination hotspots, linkage disequilibrium, and gene size . This correction can be expressed as:
where N is the number of SNPs existing within a gene. Analyses for this study were also predicated on a group of reproducible AD (RAD) genes which were previously curated from the literature .
Acquisition, labeling, and processing of brain expression data
Measurements of gene expression in the human brain were acquired from the Allen Brain Atlas (ABA). This database contains microarray data from 3702 single tissue samples extracted from six neuropathologically healthy brains (ages 24, 26, 31, 39, 49, and 57). Data derived from each sample consist of an expression vector containing expression measurements from 45,000 probes in the extracted tissue. Each sample was annotated at three different levels of granularity, which are defined for the purpose of this study as low-, mid-, or high-level structures in terms of region specificity. Principal components (PCs) of the expression vectors across all samples were computed using the prcomp method from the R programming package . Then, each sample was annotated according to the brain region based on the hierarchical labeling scheme described above. A scatterplot of the first and second PCs was produced to ascertain whether the expression vectors of samples displayed batch effects related to either the region or the brains from which samples were derived.
Determining gene expression within each brain region
Because some genes are queried by multiple probes, the mean expression of all probes mapping to each gene was computed, resulting in a gene × sample expression matrix. The expression of each gene in each brain region was adjusted using a mixed effect model approach to account for repeated sampling of both individual brains and regions in the ABA dataset . Mixed effect model specification is contained in Additional file 1.
Construction of region-specific brain co-expression networks
We constructed a co-expression network based upon correlations between all pairs of genes across the cerebrum, cerebellum, and the brain stem. In this instance, each node of the network is a gene, and each edge between genes is the absolute value of the Pearson correlation coefficient between a pair of genes. Due to the high impact of AD on the cerebrum, two additional correlation networks were created by subdividing the cerebrum into two non-overlapping subsets of regions based upon the relative time point in the course of the disease the region typically displays AD-related protein aggregation . In total, there were 79 regions of the cerebrum in the ABA dataset that showed some neuropathological evidence of AD at Braak stage 1, which we refer to as the early-stage regions, and 37 regions of the cerebrum that showed more pronounced AD pathology at Braak stage 3, which we refer to as late-stage regions. Correlations between all gene pairs were computed separately for early- and late-stage regions. In order to compare correlations between sets of genes of interest across networks, we normalized the correlations within each network. For this, we applied a novel metric, referred to as median ranking by correlation (MRC), that is derived using a “leave-one-out” strategy to normalize the distribution of correlations into a uniform distribution of ranks that is comparable across networks. Details of the MRC procedure are provided in Additional file 1.
Verifying consistency of gene rankings across correlation networks
Due to the small number of brain specimens in the ABA dataset, we tested the consistency of gene rankings by constructing a gene × gene correlation network for each brain. This procedure uses the same correlation approach described above, except that the expression levels of genes in each region were determined based on the measurements from a single brain at a time. Next, each non-RAD gene was ranked by its correlation to the RAD seed genes within each of the six individual correlation networks. Finally, a Kendall Tau rank correlation matrix was derived based upon all possible combinations of these six ranked lists.
Network-based ranking of novel AD genes
Based on the observation that the RAD genes tend to be highly correlated, we hypothesized that other genes showing a high correlation with established AD genes are likely to be AD-related genes. Therefore, a summed absolute Pearson correlation with all the RAD genes was computed for each non-RAD gene. Next, a percentile rank of each non-RAD gene-based upon these sums was computed and converted to Z-scores, which we refer to as network scores. If N genes are ranked, then the percentile rank for each gene is percentile = (rank)/(N + 1), ranging from 0 to 1. These percentiles form a uniform distribution, which are converted to Z-scores using qnorm(Percentile, lower.tail = F) in R. These network scores were then combined with genetic association Z-scores derived by a GWAS for AD risk including approximately 30,000 individuals using the Stouffer method implemented  in the meta-analysis tool METAL that was modified to equally weight both scores :
Further ranking was performed by integrating phenotypic information from gene orthologs in Flybase to focus on genes which when knocked out in a model organism result in an AD-related phenotype including premature aging, defective memory, defective aging, and oxidative stress . We also included genes which have external experimental evidence for influencing AD-related processes in human cell lines and brain. Linking fly and human gene orthologs was accomplished using the Drosophilia RNAi Screening Integrative Ortholog Prediction Tool (DIOPT) . This robust consensus mapping approach has been utilized as a functional validation strategy in studies of neurologically relevant phenotypes . The significance of network scores was determined based on the false discovery rate (FDR).
The principal component analysis revealed a potential batch effect in the six brain samples with respect to the gene expression in the three high-level brain structures (Fig. 1), noting that this analysis does not account for the non-independent gene co-expression. Further analysis revealed that most RAD genes tended to have fairly static expression in the cerebellum and brain stem regardless of the changes in the mid-level structure (Fig. 2). However, the expression for several of these genes appears more variable across the mid-level structures in the cerebrum.
Comparison of the co-expression of RAD genes across the high-level brain regions revealed higher correlation ranks (CRs) in the cerebrum (0.748) than in the brain stem (0.648) and cerebellum (0.574, Table 1). These differences appear to be due largely to a few genes including APOE and MAPT which showed much greater co-expression in the cerebrum (CR = 0.745 and 0.863, respectively) than in the cerebellum (CR = 0.280 and 0.542, respectively) and brain stem (CR = 0.216 and 0.337, respectively). Surprisingly, the CR for APP was much higher in the cerebellum (0.99) than in the brain stem (0.505) and cerebrum (0.376). Multiple RAD genes including PSEN2, EPHA1, LMX1B, TPBG, CLU, AKAP9, ZNF804B, PDGFRL, and ABCA7 were not meaningfully co-expressed with other RAD genes in any of the structures.
The MRC of the RAD genes was appreciably and nearly significantly higher for the late-stage network (0.733) than the early-stage region network (0.615, Wilcoxon signed rank p = 0.052), but the CRs for many individual genes including APOE, APP, and MAPT were similar across these two networks (Table 2). The comparison of correlation networks in the cerebrum constructed for each individual showed that RAD genes tend to have low variability in CR among individuals within this dataset (Fig. 3). The patterns of co-expression across the individual brains are moderately high and consistent with the CR values between 0.455 and 0.652 (Fig. 4).
In order to predict novel AD genes based upon the above observations, a network score was produced for each non-RAD gene using the cerebrum correlation network in which clustering of the RAD genes was strongest. These network scores were then combined with GWAS Z-scores resulting in re-ordered AD gene rankings. A normal approximation was used to evaluate the significance of the combined scores. These results were filtered using gene knockout information from Flybase to limit the focus to genes which have functional evidence for producing AD-related phenotypes. Of the remaining 654 genes after the final filtering step, there was significant evidence for two novel AD-related genes including EPS8 (FDR p = 8.77 × 10−3) and HSPA2 (FDR p = 0.245) (Table 3). Several previously reported AD genes also had high rankings but were not significant after FDR correction including ADAM10 (FDR p = 0.40) and HDAC1 (FDR p = 0.79). Only one of the top-ranked genes, RCAN1, when knocked out resulted in as many as three AD-related phenotypes in flies; however, the statistical support was modest (FDR p = 0.79).
Previous studies using correlation or other network strategies have increased discovery and understanding of the functional roles of novel disease-related genes across many biological contexts [18, 32,33,34]. In this study, we applied an integrative network strategy to capture complex relationships between RAD genes across the relevant regions of the brain and to aid the discovery of novel AD-related genes. This approach entailed integration of AD GWAS data, gene expression measures in multiple brain regions, and phenotypic information (i.e., memory and aging-related outcomes) from gene knockout studies in Drosophila . By separating the regions of the brain according to the established patterns of AD-related pathology including neurodegeneration and protein aggregation, we showed that the correlation of expression between previously established AD genes is highest in regions severely impacted by AD, noting gene expression data were derived from brains without AD pathology. In addition, we identified potential novel AD genes by numerically combining results from co-expression analysis of established AD genes and other genes in relevant brain regions with summary statistics from a large AD GWAS.
The most robust novel gene identified by our approach is EPS8. This gene encodes epidermal growth factor receptor substrate 8 which is involved in actin cytoskeleton regulation and is abundantly expressed in many brain regions . The accumulation of filamentous actin (F-actin) is associated with tau-induced neurodegeneration in Drosophila and mouse tauopathy models . The deletion of Eps8 in mice leads to a reduction in hippocampal synaptic plasticity and impaired cognitive performance . Three genes encoding heat shock proteins (HSPA2, HSPA6, and HSPB1) also emerged among our top findings. Notably, HSPA2 was also identified as related to AD in a recent network analysis in an independent dataset . Heat shock proteins have a major role in handling misfolded proteins including amyloid-β . Although the expression of heat shock protein genes has been well studied in AD , there is little evidence for the association of AD risk with polymorphisms in any members of this gene family .
Several other top-ranked genes in our study have directly or indirectly been linked to AD. ADAM10 encodes disintegrin and metalloproteinase 10 which is a synaptic enzyme that has been previously shown to limit amyloid-β1-42 peptide formation in AD. A variant in ADAM10 recently achieved genome-wide significance in one of the largest genetic studies of AD containing more than 95,000 individuals [42, 43]. The catalase protein encoded by CAT binds with amyloid and inhibition of this interaction has been reported to protect cells from toxic protein aggregation [44, 45]. Several genes in the HDAC family have been reported to impair memory in animal models, and inhibitors of several members of the HDAC gene family, including HDCA1 identified for the first time in our study as an AD candidate gene, have been gaining support as a therapeutic approach for treating AD [46,47,48]. In humans, loss of HDAC5 impairs memory function  and variants in HDAC9 have been associated with a dual outcome of neurofibrillary tangles and amyloid angiopathy . We also obtained mild evidence supporting a role for the gene encoding acetylcholinesterase (ACHE). This is a noteworthy finding in light of inconsistent and generally negative reports of association for AD with ACHE and related genes encoding choline acetyltransferase (CHAT) and butyrylcholinesterase (BCHE), despite the fact that AD is characterized by an extensive loss of cholinergic neurons from the basal forebrain area and the wide use of cholinesterase inhibitors to treat the early stages of cognitive decline . Expression of RCAN1, which encodes the regulator of calcineurin 1 and the only gene when knocked out resulted in three AD-related phenotypes in Drosophila, is increased in AD brain , and overexpression of the human RCAN1.1S isoform inserted in mice promotes early age-dependent memory and synaptic plasticity deficits and mitochondrial dysregulation leading to tau pathology .
A major motivation for our approach was to determine if the brain region-specific effects exhibited by AD can be detected using a correlation network approach. Recent work indicates that cell type compositions of the brain regions are highly variable in aging brains, so the cross-regional analysis is able to capture important properties such as changing cell fractions that may explain why the biological symptoms of AD are not uniformly present throughout the brain . The high MRC of the RAD genes in the cerebrum supports this notion, given that the cerebrum tends to be the most major structure in the brain affected by AD . Further evidence for this is also provided by the low MRC of the RAD genes in the other brain regions (brain stem, cerebellum) where the effect of AD is far less severe. Notably, these patterns appear to be consistent in our study of cognitively healthy individuals (Fig. 4).
Our findings also highlight several interesting patterns among several well-established RAD genes. We observed that expression of APOE and MAPT is highly correlated with other RAD genes in the cerebrum to the other RAD genes, but much less in the cerebellum and brain stem which is consistent with our observation of the RAD gene set as a whole. While most RAD genes are not highly correlated in the cerebellum, we observed a strong correlation among a few RAD genes, most notably, APP, which had a high CR in the cerebellum (0.99), but not in the cerebrum (0.38). APP is expressed across most regions of the brain, as evidenced in the Gene Tissue Expression (GTEx) portal . One possible explanation for a higher correlation of expression for a few RAD genes such as APP in the cerebellum is that they have an important role throughout the brain, whereas other RAD genes have a more localized role in cerebral function and health. A clearer understanding of this pattern will require a focused analysis of gene co-expression within specific regions in the cerebrum.
Interpretation of our results has several caveats. First, we analyzed a dataset that has few individuals but a high number of brain regions in which expression was measured. However, the expression patterns were consistent across individual brains in the dataset. If we had chosen instead a publicly available dataset containing a larger number of individuals but expression measurements in fewer regions, we would not have observed the high variation in the expression of the RAD genes across regions of the cerebrum. This underscores the need for larger samples of brains with expression data in more precisely defined regions. Second, the present study did not include any brains from AD individuals. Although we utilized known Braak staging to characterize regions, it is necessary to compare gene expression in the brains showing progressively severe AD pathology to determine whether the patterns observed in this study are related to AD. In addition, expression patterns for only a few genes remained significant after correction for multiple testing likely due to the relatively small sample size. Finally, because no brains from persons with AD were included in this study, we were unable to evaluate whether any of the gene co-expression patterns we identified differ between those with and without AD or are correlated with the degree of AD pathology in specific brain regions. However, the purpose of this study was to investigate whether established AD risk genes are co-expressed in the brain and whether their co-expression varies by brain region, specifically regions that are affected early in the disease compared to regions that are spared until later in the disease process. Because our findings were observed in a study of brains from relatively young individuals most of whom were autopsied several decades before typical onset of AD symptoms after age 65, our findings do provide insight about the coordinated expression of genes that are known to have a role in AD and specifically in regions of the brain that are temporally affected by the disease. Finally, the expression patterns for only a few genes remained significant after correction for multiple testing likely due to the relatively small sample size.
This work establishes a strong case for many potential follow-up investigations. Analysis of the expression in more fine-grained brain regions in larger samples including individuals with pathologically confirmed AD at various stages will allow more concise conclusions about the joint influences of multiple genes on the progression of AD from preclinical to late stages. Although highly granular regional expression data from AD brains is not readily available, efforts are in progress by the AMP-AD consortium to profile the expression of various regions of AD brains . Validated differences in cross-regional correlation patterns between healthy and AD brains would improve the understanding of the mechanisms underlying the progression of AD and inform strategies for developing more effective therapeutic targets.
Allen Brain Atlas
False discovery rate
Genome-wide association study
Median ranking by correlation
Reproducible Alzheimer disease
Single nucleotide polymorphism
Rubinsztein DC. The roles of intracellular protein-degradation pathways in neurodegeneration. Nature. 2006;443(7113):780–6.
Elbaz A, Carcaillon L, Kab S, Moisan F. Epidemiology of Parkinson’s disease. Rev Neurol (Paris). 2016;172(1):14–26.
Rowland LP. Amyotrophic lateral sclerosis. Curr Opin Neurol. 1994;7(4):310–5.
Wenk GL. Neuropathologic changes in Alzheimer’s disease: potential targets for treatment. J Clin Psychiatry. 2006;67(Suppl 3):3–7.
Browne SE, Bowling AC, MacGarvey U, Baik MJ, Berger SC, Muqit MM, et al. Oxidative damage and metabolic dysfunction in Huntington’s disease: selective vulnerability of the basal ganglia. Ann Neurol. 1997;41(5):646–53.
Hohenfeld C, Werner CJ, Reetz K. Resting-state connectivity in neurodegenerative disorders: is there potential for an imaging biomarker? Neuroimage Clin. 2018;18:849–70.
Gur RE, McGrath C, Chan RM, Schroeder L, Turner T, Turetsky BI, et al. An fMRI study of facial emotion processing in patients with schizophrenia. Am J Psychiatry. 2002;159(12):1992–9.
Yao Z, Wang L, Lu Q, Liu H, Teng G. Regional homogeneity in depression and its relationship with separate depressive symptom clusters: a resting-state fMRI study. J Affect Disord. 2009;115(3):430–8.
Mu Y, Gage FH. Adult hippocampal neurogenesis and its role in Alzheimer’s disease. Mol Neurodegener. 2011;6:85.
Ross CA, Poirier MA. Protein aggregation and neurodegenerative disease. Nat Med. 2004;10(Suppl):S10–7.
Braak H, Del Tredici K. The pathological process underlying Alzheimer’s disease in individuals under thirty. Acta Neuropathol. 2011;121:171–81.
Braak H, Thal DR, Ghebremedhin E, Del Tredici K. Stages of the pathologic process in Alzheimer disease: age categories from 1 to 100 years. J Neuropathol Exp Neurol. 2011;70:960–9.
Pierson E, GTEx Consortium, Koller D, Battle A, Mostafavi S, Ardlie KG, et al. Sharing and specificity of co-expression networks across 35 human tissues. PLoS Comput Biol. 2015;11(5):e1004220.
Mostafavi S, Gaiteri C, Sullivan SE, White CC, Tasaki S, Xu J, et al. A molecular network of the aging human brain provides insights into the pathology and cognitive decline of Alzheimer’s disease. Nat Neurosci. 2018;21(6):811–9.
Collado-Torres L, Burke EE, Peterson A, Shin J, Straub RE, Rajpurohit A, et al. Regional heterogeneity in gene expression, regulation, and coherence in the frontal cortex and hippocampus across development and schizophrenia. Neuron. 2019;103(2):203–16.
Moore JH, Asselbergs FW, Williams SM. Bioinformatics challenges for genome-wide association studies. Bioinformatics. 2010;26(4):445–55.
Shen L, Kim S, Risacher SL, Nho K, Swaminathan S, West JD, et al. Whole genome association study of brain-wide imaging phenotypes for identifying quantitative trait loci in MCI and AD: a study of the ADNI cohort. Neuroimage. 2010;53(3):1051–63.
Novarino G, Fenstermaker AG, Zaki MS, Hofree M, Silhavy JL, Heiberg AD, et al. Exome sequencing links corticospinal motor neuron disease to common neurodegenerative disorders. Science. 2014;343(6170):506–11.
Vanunu O, Magger O, Ruppin E, Shlomi T, Sharan R. Associating genes and protein complexes with disease via network propagation. PLoS Comput Biol. 2010;6(1):e1000641.
Cowen L, Ideker T, Raphael BJ, Sharan R. Network propagation: a universal amplifier of genetic associations. Nat Rev Genet. 2017;18(9):551–62.
Lancour D, Naj A, Mayeux R, Haines JL, Pericak-Vance MA, Schellenberg GD, et al. One for all and all for one: improving replication of genetic studies through network diffusion. PLoS Genet. 2018;14(4):e1007306.
Lambert J-C, Ibrahim-Verbaas CA, Harold D, Naj AC, Sims R, Bellenguez C, et al. Meta-analysis of 74,538 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat Genet. 2013;45(12):1452–8.
Segrè AV, DIAGRAM Consortium; MAGIC investigators, Groop L, Mootha VK, Daly MJ, Altshuler D. Common Inherited variation in mitochondrial genes Is not enriched for associations with Type 2 diabetes or related glycemic traits. Plos Genetics. 2010;6(8):e1001058.
R Development Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2019.
Sheiner LB, Grasela TH. An introduction to mixed effect modeling - concepts, definitions, and justification. J Pharmacokinet Biopharm. 1991;19(3):S11–24.
Braak H, Braak E. Neuropathological stageing of Alzheimer-related changes. Acta Neuropathol. 1991;82(4):239–59.
Stouffer SA, Suchman EA, DeVinney KC, Star SA, Williams RM Jr. The American soldier: volume I, adjustment during army life. Princeton: Princeton University Press; 1949.
Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26(17):2190–1.
Campanella C, Pace A, Caruso Bavisotto C, Marzullo P, Marino Gammazza A, Buscemi S, Palumbo Piccionello A. Heat shock proteins in Alzheimer's disease: role and targeting. Int J Mol Sci. 2018;19(9):2603.
Hu Y, Flockhart I, Vinayagam A, Bergwitz C, Berger B, Perrimon N, Mohr SE. An integrative approach to ortholog prediction for disease-focused and other functional studies. BMC Bioinformatics. 2011;12:357.
Satizabal CL, Adams HHH, Hibar DP, White CC, Knol MJ, Stein JL, et al. Genetic architecture of subcortical brain structures in 38,851 individuals. Nat Genet. 2019;51(11):1624–36.
Goh KI, Cusick ME, Valle D, Childs B, Vidal M, Barabasi AL. The human disease network. Proc Natl Acad Sci U S A. 2007;104(21):8685–90.
Liu M, Liberzon A, Kong SW, Lai WR, Park PJ, Kohane IS, Kasif S. Network-based analysis of affected biological processes in type 2 diabetes models. PLoS Genet. 2007;3(6):e96.
Zhang B, Gaiteri C, Bodea LG, Wang Z, McElwee J, Podtelezhnikov AA, et al. Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer’s disease. Cell. 2013;153(3):707–20.
Disanza A, Carlier MF, Stradal TE, Didry D, Frittoli E, Confalonieri S, et al. Eps8 controls actin-based motility by capping the barbed ends of actin filaments. Nat Cell Biol. 2004;6(12):1180–8.
Fulga TA, Elson-Schwab I, Khurana V, Steinhilb ML, Spires TL, Hyman BT, Feany MB. Abnormal bundling and accumulation of F-actin mediates tau-induced neuronal degeneration in vivo. Nat Cell Biol. 2007;9(2):139–48.
Wang YT, Huang CC, Lin YS, Huang WF, Yang CY, Lee CC, et al. Conditional deletion of Eps8 reduces hippocampal synaptic plasticity and impairs cognitive function. Neuropharmacology. 2017;112(Pt A):113–23.
Petyuk VA, Chang R, Ramirez-Restrepo M, Beckmann ND, Henrion MYR, Piehowski PD, et al. The human brainome: network analysis identifies HSPA2 as a novel Alzheimer’s disease target. Brain. 2018;141(9):2721–39.
Campanella C, Pace A, Caruso Bavisotto C, Marzullo P, Marino Gammazza A, Buscemi S, Palumbo PA. Heat shock proteins in Alzheimer’s disease: role and targeting. Int J Mol Sci. 2018;19(9).
Hamos JE, Oblas B, Pulaski-Salo D, Welch WJ, Bole DG, Drachman DA. Expression of heat shock proteins in Alzheimer’s disease. Neurology. 1991;41(3):345–50.
Kunkle BW, Grenier-Boley B, Sims R, Bis JC, Damotte V, Naj AC, et al. Genetic meta-analysis of diagnosed Alzheimer’s disease identifies new risk loci and implicates Aβ, tau, immunity and lipid processing. Nat Genet. 2019;51(3):414–30.
Colciaghi F, Marcello E, Borroni B, Zimmermann M, Caltagirone C, Cattabeni F, et al. Platelet APP, ADAM 10 and BACE alterations in the early stages of Alzheimer disease. Neurology. 2004;62(3):498–501.
Milton NG. Amyloid-beta binds catalase with high affinity and inhibits hydrogen peroxide breakdown. Biochem J. 1999;344(Pt 2):293–6.
Habib LK, Lee MT, Yang J. Inhibitors of catalase-amyloid interactions protect cells from beta-amyloid-induced oxidative stress and toxicity. J Biol Chem. 2010;285(50):38933–43.
Yang SS, Zhang R, Wang G, Zhang YF. The development prospection of HDAC inhibitors as a potential therapeutic direction in Alzheimer’s disease. Transl Neurodegener. 2017;6:19.
Janczura KJ, Volmar CH, Sartor GC, Rao SJ, Ricciardi NR, Lambert G, et al. Inhibition of HDAC3 reverses Alzheimer’s disease-related pathologies in vitro and in the 3xTg-AD mouse model. Proc Natl Acad Sci U S A. 2018;115(47):E11148–57.
Agis-Balboa RC, Pavelka Z, Kerimoglu C, Fischer A. Loss of HDAC5 impairs memory function: implications for Alzheimer’s disease. J Alzheimers Dis. 2013;33(1):35–44.
Kilgore M, Miller CA, Fass DM, Hennig KM, Haggarty SJ, Sweatt JD, Rumbaugh G. Inhibitors of class 1 histone deacetylases reverse contextual memory deficits in a mouse model of Alzheimer’s disease. Neuropsychopharmacology. 2010;35(4):870–80.
Chung J, Zhang X, Allen M, Wang X, Ma Y, Beecham G, et al. Genome-wide pleiotropy analysis of neuropathological traits related to Alzheimer’s disease. Alzheimers Res Ther. 2018;10(1):22.
Hampel H, Mesulam MM, Cuello AC, Farlow MR, Giacobini E, Grossberg GT, et al. The cholinergic system in the pathophysiology and treatment of Alzheimer’s disease. Brain. 2018;141(7):1917–33.
Wu Y, Ly PT, Song W. Aberrant expression of RCAN1 in Alzheimer’s pathogenesis: a new molecular mechanism and a novel drug target. Mol Neurobiol. 2014;50(3):1085–97.
Wong H, Levenga J, Cain P, Rothermel B, Klann E, Hoeffer C. RCAN1 overexpression promotes age-dependent mitochondrial dysregulation related to neurodegeneration in Alzheimer’s disease. Acta Neuropathol. 2015;130(6):829–43.
Carithers LJ, Moore HM. The Genotype-Tissue Expression (GTEx) Project. Biopreserv Biobank. 2015;13(5):307–8.
Wang M, Roussos P, McKenzie A, Zhou X, Kajiwara Y, Brennand KJ, et al. Integrative network analysis of nineteen brain regions identifies molecular signatures and networks underlying selective regional vulnerability to Alzheimer’s disease. Genome Med. 2016;8(1):104.
We thank Drs. Claire Hu and Norbert Perrimon for sharing detailed Flybase information related to the curation and mapping of phenotype-associated gene knockouts in Drosophila.
Availability of data and supporting materials
Summarized genome-wide association study results are available from the National Institute on Aging Genetics of Alzheimer’s Disease Data Storage Site (https://www.niagads.org/).
This work was supported by the National Institute of Aging grants U01-AG032984, RF1-AG057519, R01-AG048927, and P30-AG13846.
Ethics approval and consent to participate
This study was approved by the Boston University Institutional Review Board.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lancour, D., Dupuis, J., Mayeux, R. et al. Analysis of brain region-specific co-expression networks reveals clustering of established and novel genes associated with Alzheimer disease. Alz Res Therapy 12, 103 (2020). https://doi.org/10.1186/s13195-020-00674-7
- Alzheimer disease
- Gene network analysis
- Brain regions
- Genome-wide association study