Alzheimer’s disease master regulators analysis: search for potential molecular targets and drug repositioning candidates
Alzheimer's Research & Therapy volume 10, Article number: 59 (2018)
Alzheimer’s disease (AD) is a multifactorial and complex neuropathology that involves impairment of many intricate molecular mechanisms. Despite recent advances, AD pathophysiological characterization remains incomplete, which hampers the development of effective treatments. In fact, currently, there are no effective pharmacological treatments for AD. Integrative strategies such as transcription regulatory network and master regulator analyses exemplify promising new approaches to study complex diseases and may help in the identification of potential pharmacological targets.
In this study, we used transcription regulatory network and master regulator analyses on transcriptomic data of human hippocampus to identify transcription factors (TFs) that can potentially act as master regulators in AD. All expression profiles were obtained from the Gene Expression Omnibus database using the GEOquery package. A normal hippocampus transcription factor-centered regulatory network was reconstructed using the ARACNe algorithm. Master regulator analysis and two-tail gene set enrichment analysis were employed to evaluate the inferred regulatory units in AD case-control studies. Finally, we used a connectivity map adaptation to prospect new potential therapeutic interventions by drug repurposing.
We identified TFs with already reported involvement in AD, such as ATF2 and PARK2, as well as possible new targets for future investigations, such as CNOT7, CSRNP2, SLC30A9, and TSC22D1. Furthermore, Connectivity Map Analysis adaptation suggested the repositioning of six FDA-approved drugs that can potentially modulate master regulator candidate regulatory units (Cefuroxime, Cyproterone, Dydrogesterone, Metrizamide, Trimethadione, and Vorinostat).
Using a transcription factor-centered regulatory network reconstruction we were able to identify several potential molecular targets and six drug candidates for repositioning in AD. Our study provides further support for the use of bioinformatics tools as exploratory strategies in neurodegenerative diseases research, and also provides new perspectives on molecular targets and drug therapies for future investigation and validation in AD.
Alzheimer’s disease (AD) is the most prevalent neurodegenerative disease and the major cause of dementia. In the United States about 10% of people over the age of 65 years have Alzheimer’s dementia, and the worldwide prevalence of the disease ranges from 4 to 8% [1, 2]. A total of roughly 46 million AD cases is estimated around the world and related cost are about U$800 billion per year [2, 3].
This neurodegenerative disease causes gradual loss of brain volume and synaptic dysfunction, leading to a progressive memory and reasoning impairment followed by global cognitive decline and, ultimately, dementia [4, 5]. AD is characterized by its histopathological hallmarks, which includes deposits of amyloid-β (Aβ) plaques and neurofibrillary tangles composed of hyperphosphorylated tau . Recent advances in our ability to detect AD pathophysiology using imaging biomarkers currently allow the identification of Aβ and tau pathology in living individuals [7, 8]. By contrast, few advancements have been made in terms of drug treatments, which currently are available only for ameliorating symptoms .
Sporadic AD, also called late-onset AD, represents the vast majority of cases (> 95%) and is recognized as a multifactorial, complex disease . Apolipoprotein E isoform ε4 (APOEε4) is the main susceptibility gene for AD, with a threefold increase in AD risk for one allele and 12-fold increase for two alleles . Genome-wide association studies have identified more than 20 AD risk genes and several disease-associated pathways . However, the AD risk genes identified so far are neither necessary nor sufficient for disease onset . Meanwhile, evidence suggest that nongenetic factors, such as cerebrovascular disease, diabetes, and obesity, also increase the risk of developing AD . Furthermore, gene expression profiling studies in AD brains have shown many genes working together in relevant altered biological pathways in the disease, leading to a growing acceptance that AD results from the impairment of several complex mechanisms at once that have not yet been fully elucidated .
In keeping with this, the high rate of failure in the development of AD-modifying therapies seems to be a consequence of the incomplete knowledge about the underlying mechanisms of the disease. Based on this, the use of new approaches to study the disease pathophysiology and search for alternative therapeutic targets are urgently required [9, 14, 15]. The use of integrative strategies, such as regulatory networks, for analyzing high-throughput expression data have produced significant knowledge towards the elucidation of biological mechanisms underlying complex diseases, such as cancer and obesity . Furthermore, it has been observed that regulatory networks inferred by reverse engineering algorithms can provide sufficient accuracy to estimate the impact of transcription factors (TFs) on phenotype transitions according to their transcriptional targets, and to identify the ones that are acting as master regulators (MRs) of diseases . Many approaches have shown that TFs can operate as key elements in the phenotypic determination by regulating large groups of transcriptional targets associated with complex cellular processes [17,18,19,20]. Therefore, the analysis of expression profiling data using a TF-centered regulatory networks approach seems an interesting strategy to study the mechanisms and common drivers associated with AD.
In this study, gene expression data available in the Gene Expression Omnibus repository (GEO; http://www.ncbi.nlm.nih.gov/geo/) was used to infer a transcriptional regulatory network, through reverse engineering, for the human hippocampus, a region that undergoes high rates of volume loss in AD. Afterwards, expression data from AD case-control studies of the same region were used to identify MRs potentially modulating phenotypic changes from a normal to a pathological scenario. Moreover, the prospection of new drug candidates to treat AD patients was carried out by a connectivity map approach using the inferred regulatory units of MR candidates.
Microarray data acquisition
A normal human brain expression dataset was obtained from the GEO database under the accession number GSE60862 . AD case-control microarray studies from hippocampal samples were acquired from GEO under accession numbers GSE5281 [22, 23], GSE29378 , GSE36980 , and GSE48350 . Table 1 summarizes the data information from the selected GEO datasets used in this study. Each expression dataset was treated and analyzed independently (Additional file 1: Figure S1).
Region-specific transcription network inference
The genome-wide region-specific transcriptional network (TN) centered on TFs and their predicted target genes were inferred using the normal brain hippocampus (HIP) expression data from GSE60862. The groups of inferred target genes associated with each TF are hereinafter referred as its regulatory unit. These computations were performed using the RTN package, which is designed to reconstruct and analyze TNs based on the mutual information (MI), a measure that evaluates dependencies between two random variables, using the ARACNe (Algorithm for the Reconstruction of Accurate Cellular Networks) method. Briefly, the regulatory structure of the network is derived by mapping significant associations between known TFs and all potential targets. Interactions below a minimum MI threshold are eliminated by a permutation step and unstable interactions are additionally removed by bootstrap to create a consensus bootstrap network. In a final step, the data processing inequality algorithm is applied with null tolerance to eliminate interactions that are likely to be mediated by another TF. Here, we used the package’s default number of permutations and number of bootstraps (1000 permutations and 100 bootstraps), but with a p value cutoff of 0.001. The resultant network will be hereinafter referred to as HIP-TN [18, 27, 28].
All computational analyses were performed in R statistical environment . Network figures were constructed with the RedeR graphical platform for exploration of biological networks , and other plots were constructed using ggplot2 .
Master regulators and gene set enrichment analysis
After the hippocampus transcriptional regulatory network (HIP-TN) inference, we applied the master regulator analysis (MRA) algorithm described by Carro et al.  to the regulatory units comprised of at least 100 targets. The algorithm computes the statistical significance of the overlap between the regulatory units in HIP-TN and the differentially expressed genes (false discovery rate (FDR)-adjusted p value < 0.05) obtained from each AD study, corrected for multiple comparisons. We then selected the regulatory units of the TFs showing significant enrichment of differentially expressed target genes in three or more studies, which we termed MR candidates.
Two-tail gene set enrichment analysis (GSEA) was also performed using the RTN package with 1000 permutations, as previously described . Briefly, the groups of target genes for each MR (regulatory units) were split into positive and negative mode of action targets using Pearson’s correlation. Next, the association of each subgroup was assessed by GSEA statistics in each ranked phenotype, resulting in independent enrichment scores (Es), with two enrichment distributions. Additionally, a differential enrichment was performed among subgroups (EsA-EsB) where maximum deviation from zero near opposite extremes is desirable for a clear association. Thus, a highly positive differential score implies that the regulatory unit is induced in the disease phenotype, while a highly negative differential score indicates that the regulatory unit is repressed in the disease phenotype. The two-tail GSEA computation p value cutoff was set to 0.05 and 1000 permutations were used.
The differentially expressed genes used in the MRA and the log fold change (logFC) metric used to obtain the ranked phenotypes required for the GSEA were computed using the Bioconductor package limma .
Connectivity map drug profiling approach
The previously identified MR candidate regulatory units were queried in the Connectivity Map online tool (The CMap build02; www.broadinstitute.org/cmap/) using the GSEA algorithm described by Lamb et al. . This tool compares queried signature with gene expression profile database of several cell lines after treatment with approximately 1000 compounds, most of which are FDA approved. Drugs whose signature opposes the disease signature are assumed to have a therapeutic potential.
For this, we first selected the MR candidates with two-tail GSEA p values less or equal to 0.01. Next, for each case-control study, the differentially expressed targets of these MR candidates (adjusted p value < 0.05) were filtered, grouped, tagged according to the logFC metric, converted to Affymetrix probe identifiers, and submitted as input for the cMap webtool. Then we obtained a connectivity map of drug-phenotype association for each case-control study.
Human hippocampus transcriptional regulatory network reconstruction
HIP-TN was computed from a normal brain gene expression dataset (GSE60862) using the reverse engineering ARACNe algorithm. Transcripts were classified as transcription factors when annotated in the Gene Ontology with the identifier GO:0003700 (transcription factor activity, sequence-specific DNA binding). Among a total of 20,311 transcripts in the dataset, 766 were annotated as TFs under GO:0003700. From these, 469 were classified as TFs with more than 25 inferred target genes. The resultant HIP-TN, comprising 132 regulatory units with more than 100 targets, was used for further analyses. Figure 1a shows the inferred HIP-TN inside the blue container, where each node symbolizes a TF regulatory unit and node sizes correspond to the number of predicted targets for each TF (Additional file 2: Table S1). Regulatory units with less than 100 targets are represented in black outside the blue container.
Hippocampus AD master regulator inference
Microarray gene expression from AD case-control studies available in GEO (GSE5281, GSE29378, GSE36980, and GSE48350) were used to obtain disease MR candidates considering the normal HIP-TN previously inferred. MRA was performed to evaluate HIP regulatory units enriched with genes differentially expressed between the two phenotypes (disease and control). Only regulatory units significantly enriched in at least three case-control studies were considered as MR candidates. These analyses resulted in the identification of 34 MR candidates (Fig. 1b) (Additional file 3: Table S2).
Two-tail GSEA was performed to infer the activation state of each MR candidate. The outcome of this analysis showed 14 MR candidates that were significantly repressed and 2 MR candidates that were significantly activated in AD (FDR adjusted p value ≤ 0.05) (Fig. 2a). This means that targets from the repressed MR candidates had predicted positive TF-target expression association under normal conditions but were decreased in the disease. To the contrary, targets with inferred negative TF-target association under normal conditions had increased expression in the pathology. For the activated MR candidates, the inferred positive or negative TF target expression associations do not reverse during AD. The remaining 18 MR candidates did not present statistically significant results regarding their activation states and thus were not considered for the next steps.
The AD subregulatory network graph in Fig. 2b and Additional file 4 (Table S3) shows the association pattern between those MR candidates with significant alteration of the activation state in the disease state compared with control. The nodes with the highest degrees of connectivity in this network correspond to the MRs ATF2 (activating transcription factor 2) and PARK2 (Parkin RBR E3 ubiquitin protein ligase). The number of common targets between any two MRs are represented by the connector line widths as assessed by the Jaccard coefficient and indicates that certain variations in target expression may be a result of the contiguous regulatory action of two or more TFs.
The connectivity map approach was used to search for drugs with therapeutic repurposing potential in AD. The 16 MR candidate regulatory units identified in the previous analysis were grouped, and their up- or downregulated differentially expressed targets were selected for each AD case-control studies and used as an input in the webtool (Fig. 3a). The consensus drugs are consistently present in at least two case-control studies (p ≤ 0.05). Six drugs were negatively associated with AD and assumed to have a therapeutic potential: Cefuroxime, Cyproterone, Dydrogesterone, Metrizamide, Trimethadione, and Vorinostat. Additionally, seven drugs were positively associated and thus considered AD mimetic: Calmidazolium, Ciclosporin, Disulfiram, Fluspirilene, Puromycin, Quipazine, and Spiperone (Fig. 3b) (Additional file 5: Table S4).
Based on reverse engineering coexpression regulatory network reconstruction for the human HIP, we identified a range of transcription factors that acts on large regulatory units, therefore being potentially important for the functionality of this region. Furthermore, from these regulatory units, we selected those which differentially expressed inferred target genes were overrepresented in AD versus control. Interestingly, among the ten largest hippocampal regulatory units, seven of them present differential expression when comparing AD versus control (Additional file 2: Table S1).
The 34 MR candidates selected in this study were annotated with the GO term “transcription factor activity, sequence-specific DNA binding”. Among them, only 5 MR candidates, namely KAT7, MTA3, RREB1, TSC22D1, and ZNF287, do not have this GO term assigned by a curator. Moreover, some of the MR candidates have also been associated with other transcription regulation functions, such as transcription corepressor activity (GO:0003714) and transcription coactivator activity (GO:0003713), bringing forth the possibility that the influence of each MR candidate on the expression levels of its inferred regulatory units may be related to expression regulation mechanisms other than the direct DNA binding activity. Nevertheless, it is important to note that, for the purpose of finding MRs for the disease, this distinction is not necessarily relevant once the sought MR function relates to the expression regulation as a broad and diverse regulatory phenomenon. The additional GO terms annotated for each MR candidate are described in Additional file 3 (Table S2).
Among the TFs inferred as MRs of the disease, several, such as ATF2 and PARK2, have already had their relationship with AD previously reported [35,36,37]. Indeed, both ATF2 and PARK2 showed a high degree of connectivity via their inferred targets in the AD regulatory subnetwork, thus acting as potential hubs, and predicted elements of great importance for network maintenance and robustness, with their regulatory units repressed in the disease.
The PARK2 gene encodes an E3 ubiquitin ligase and it is one of the genes involved in autosomal recessive juvenile parkinsonism . In addition to its function in the ubiquitin proteasome system, PARK2 is also involved in the regulation of gene expression, modulating genes associated with apoptosis or cellular stress reactions . Furthermore, PARK2 also acts as a direct transcriptional repressor of p53 promotor activity, thus modulating cell death pathways . Remarkably, PARK2 was also shown to directly mediate expression of two proteins related to the amyloidogenic pathway, Presenilin 1 and Presenilin 2, which are components of the γ-secretase complex . Mutations in the coding regions of these two proteins are related to AD familial cases .
ATF2 is a member of the ATF/CREB family that regulates gene expression through homodimerization or heterodimerization with several other protein partners. However, the role of each dimer in target regulation is very difficult to determine and the knowledge about them is still limited . ATF2 is activated by several cell-damaging stimuli, such as cisplatin-induced genotoxic stress and ultraviolet (UV) radiation exposure [43, 44]. This TF regulates the expression of genes involved in important cellular processes also altered in AD, such as inflammatory signaling, apoptotic pathway, DNA damage response, and cell cycling control, being regarded as an early stress response protein [42, 45,46,47,48]. In agreement with the results obtained in this study, reduced expression of ATF2 has been shown for the CA1 to CA4 hippocampal areas, granule cells of the dentate gyrus, and adjacent entorhinal cortex in AD patients [36, 37].
Furthermore, nuclear availability of ATF2 and PARK2 are strongly influenced by stressing factors. ATF2 translocation from the nucleus to the cytoplasm was found to be increased in situations of cellular stress and disease states, leading to cell death triggering by induced opening of mitochondrial membrane pores . PARK2 solubility is compromised by oxidative and nitrosative stress and aging, in some cases showing behavioral patterns equivalent to those PARK2 mutations correlated with Parkinson’s disease . Thus, the reduced nuclear availabilities of both TFs in response to severe stress may account for the target expression reductions identified in this study.
We also identified novel TFs that seem to be involved in AD: CNOT7 (CCR4-NOT transcription complex subunit 7), CSRNP2 (cysteine and serine rich nuclear protein 2), SLC30A9 (solute carrier family 30 member 9), and TSC22D1 (TSC22 domain family member 1). These MR candidates have also shown a high degree of connectivity in the AD subregulatory network, being also potentially important for this disease. In the following sections we discuss each of these MRs with a brief description of its known functions.
CNOT7 is a catalytic component of one of the major mRNA deadenylase complexes (CCR4-NOT). It has an antiproliferative function dependent both on its deadenylase activity and its association with BTG1 (BTG anti-proliferation factor 1) [49, 50]. SLC30A9, also called ZnT9, belongs to a family of zinc transporters. This protein contains a motif for interaction with nuclear receptors, apparently migrating to the nucleus in a cell cycle-dependent manner [51, 52]. It has been shown that SLC30A9 acts as a hormone-dependent nuclear receptor coactivator and also participates in the Wnt signaling pathway by interacting with β-catenin [53, 54]. TSC22D1 is the most studied among these transcription factors due to its tumor suppressor activity. It was isolated as a transforming growth factor (TGF)-β-induced transcript which encodes a leucine-zipper transcription factor and has transcriptional repressor activity [55, 56]. TSC22D1 has been shown to be a p53-positive regulator, inhibiting its degradation. Furthermore, it also inhibits cell proliferation, promoting apoptosis when overexpressed .
Although it was not possible to determine the activation state of several other MRs found, a handful of studies directly correlating MEF2A (myocyte enhancer factor 2A), STAT1 (signal transducer and activator of transcription 1), and YY1 (Yin and Yang 1 protein) TFs to AD are available in the literature. It was shown that YY1 is directly involved in the regulation of important AD-related genes, such as BACE1 (Beta-secretase 1) and APH1A (aph-1 homolog A, gamma-secretase subunit), which have binding sites for YY1 in their promoter regions [58,59,60]. STAT1 also has a role in controlling the gene expression of BACE1, binding to its promoter region, and can be upregulated by Aβ, characterizing a positive feedback loop that could lead to the progressive increase of production and further accumulation of Aβ [61,62,63]. Regarding MEF2A, Burton et al.  and Gonzalez et al.  have suggested that deregulation in the control of these TF activation pathways could be associated with increased risk of developing AD. Additionally, the genes MEF2C and CELF1, identified by genome-wide association studies as having a small effect on AD risk , were inferred as MEF2A and YY1 targets, respectively (data not shown), which reinforces the idea that these genes are part of a broad and complex context and that to discuss their roles in the whole scenario could be a much more constructive approach.
Neuronal loss and astrocytosis are well-known events related to AD, and both have been observed in postmortem brains of AD patients [66,67,68]. A reduction in the neuronal population is directly related to the progression of hippocampal atrophy, to the severity of the dementia [69, 70], and to the Braak stage of the disease . The presence of astrocytosis in AD has also been described, and it is thought to be related either to the proliferation of astrocytes to replace dying neurons, or to an increased activity of these cells in an effort to scavenge the toxic Aβ peptides [72, 73]. Although astrocytosis is known to be essential for tissue repair and early mitigation of lesions, it can also lead to further deleterious effects, either by amplifying the inflammatory response  or by diminishing the trophic support for neurons [74, 75].
To investigate whether our results could be related to these histopathological alterations, we conducted a preliminary GSEA to compare the expression levels of the 34 regulatory units (HIP-TN) in mouse neuron versus astrocyte data from the microarray dataset GSE9566 . We found that all 34 regulatory units were enriched with differentially expressed genes in astrocytes compared with neurons. Notably, 16 regulatory units followed the same pattern of activation found in the AD case-control analysis, whereas 18 showed nonsignificant states of activation (Additional file 6: Figure S2). These findings indicate that our results can be, at least in part, a reflex to an increase in the influence of astrocyte-related regulatory units in the overall signature of the disease, which may be compatible with the astrocyte hyperactivation and proliferation hypothesis in AD. Therefore, a reversion of the inferred transcriptional signature as a whole can be a promising strategy to alleviate deleterious effects potentially mediated by these responses.
The transition from a single-gene approach to a network-centric view is seen as a new path in the search for pharmacological strategies for complex diseases . In addition, drug repositioning has been shown to be a cheaper and faster alternative method for the development of new therapeutic regimens . The connectivity map proposal enables us to combine both of these paradigms by incorporating a data-driven method for exploring transcriptional profile alterations with drug effects on expression. We applied a connectivity map adaptation centered on transcription factor regulatory units and obtained six FDA-approved drug candidates that seem to revert AD phenotype (Cefuroxime, Cyproterone, Dydrogesterone, Metrizamide, Trimethadione, and Vorinostat). Interestingly, these drugs have several self-related or class-related neuroprotective effects previously reported in the literature. Notably, Cyproterone and Vorinostat have already been shown to be neuroprotective in AD models. Cyproterone is an antiandrogen that antagonizes androgen-mediated gene expression, although it exerts a testosterone-like neuroprotective effect against Aβ toxicity in primary neuronal cultures by an androgen receptor activation-dependent mechanism . Vorinostat is a histone deacetylase inhibitor (HDACi) used for cancer treatment, and it has been shown to restore memory deficits in an AD animal model and protects against Aβ toxicity in an AD cell model . This drug is currently at phase 1 clinical trial for assessment of its memory performance improvement capabilities in AD  (www.clinicaltrial.gov). Furthermore, there are several studies showing the role of HDACis in the reduction of inflammatory mediator expression, excitotoxicity, and oxidative stress, as well as enhancement of neurotrophic factor expression, which are relevant pathways for AD . Trimethadione is a T-type calcium channel inhibitor used as an anticonvulsant drug. It has been reported as a neuroprotective compound leading to both prevention of calcium homeostasis impairment, potentially associated with the onset of AD, and reduction of age-related degenerative effects in animal models [82, 83].
Cefuroxime is a second-generation cephalosporin antibiotic that can cross the blood-brain barrier, and Dydrogesterone is a progestogen usually administered in conditions associated with progesterone deficiency [84, 85]. Although neither of them has reported neuroprotective effects, there are several class-related central nervous system benefits associated with them in the literature [86,87,88]. Finally, Metrizamide, a radiocontrast shown to effectively inhibit the brain hexokinase, has a recent pharmacodynamic study exploring its effects on neuronal function .
Systems biology is an integrative, hypothesis-free approach based on biological component interactions and represents an interesting avenue to study complex diseases. Indeed, regulatory networks centered in TF have already been shown effective in identifying cancer drivers [17, 18, 32]. Furthermore, this approach is also gradually becoming the methodology of choice to study multifactorial complex neurodegenerative diseases [19, 90, 91]. Herein, employing a systems approach, we identified several TFs previously related to the disease as well as novel potential targets to be investigated. In addition, new therapeutic strategies using drug repositioning were prospected from the obtained transcriptional signatures. Nevertheless, further studies using both in vitro and in vivo models are required to fully evaluate the impact and benefits of these findings in AD.
APH-1 homolog A, gamma secretase subunit
Apolipoprotein E isoform ε4
Algorithm for the Reconstruction of Accurate Cellular Networks
Activating transcription factor 2
CCR4-NOT transcription complex subunit 7
Cysteine and serine rich nuclear protein 2
False discovery rate
Gene Expression Omnibus
Gene set enrichment analysis
Histone deacetylase inhibitor
Hippocampus transcriptional network
Log fold change
Myocyte enhancer factor 2A
Master regulator analysis
Parkin RBR E3 ubiquitin protein ligase
Solute carrier family 30 member 9
Signal transducer and activator of transcription 1
TSC22 domain family member 1
Yin and Yang 1 protein
Alzheimer’s Association. Alzheimer’s disease facts and figures. Alzheimers Dement. 2017;13:325–73.
Prince MJ: World Alzheimer report 2015: the global impact of dementia: an analysis of prevalence, incidence, cost and trends. Alzheimers Dis Int; 2015.
Group GBDNDC. Global, regional, and national burden of neurological disorders during 1990–2015: a systematic analysis for the global burden of disease study 2015. Lancet Neurol. 2017;16:877–97.
Dubois B, Hampel H, Feldman HH, Scheltens P, Aisen P, Andrieu S, Bakardjian H, Benali H, Bertram L, Blennow K, et al. Preclinical Alzheimer’s disease: definition, natural history, and diagnostic criteria. Alzheimers Dement. 2016;12:292–323.
McKhann GM, Knopman DS, Chertkow H, Hyman BT, Jack CR Jr, Kawas CH, Klunk WE, Koroshetz WJ, Manly JJ, Mayeux R, et al. The diagnosis of dementia due to Alzheimer’s disease: recommendations from the National Institute on Aging-Alzheimer’s Association workgroups on diagnostic guidelines for Alzheimer’s disease. Alzheimers Dement. 2011;7:263–9.
Reitz C, Mayeux R. Alzheimer disease: epidemiology, diagnostic criteria, risk factors and biomarkers. Biochem Pharmacol. 2014;88:640–51.
Leuzy A, Zimmer ER, Heurling K, Rosa-Neto P, Gauthier S. Use of amyloid PET across the spectrum of Alzheimer’s disease: clinical utility and associated ethical issues. Amyloid. 2014;21:143–8.
Zimmer ER, Leuzy A, Gauthier S, Rosa-Neto P. Developments in tau PET imaging. Can J Neurol Sci. 2014;41:547–53.
Appleby BS, Cummings JL. Discovering new treatments for Alzheimer’s disease by repurposing approved medications. Curr Top Med Chem. 2013;13:2306–27.
Karch CM, Goate AM. Alzheimer’s disease risk genes and mechanisms of disease pathogenesis. Biol Psychiatry. 2015;77:43–51.
Reitz C. Genetic diagnosis and prognosis of Alzheimer’s disease: challenges and opportunities. Expert Rev Mol Diagn. 2015;15:339–48.
Guerreiro RJ, Gustafson DR, Hardy J. The genetic architecture of Alzheimer’s disease: beyond APP, PSENs and APOE. Neurobiol Aging. 2012;33:437–56.
Huang Y, Mucke L. Alzheimer mechanisms and therapeutic strategies. Cell. 2012;148:1204–22.
Mei H, Feng G, Zhu J, Lin S, Qiu Y, Wang Y, Xia T. A practical guide for exploring opportunities of repurposing drugs for CNS diseases in systems biology. Methods Mol Biol. 2016;1303:531–47.
Brunden KR, Trojanowski JQ, Lee VM. Advances in tau-focused drug discovery for Alzheimer’s disease and related tauopathies. Nat Rev Drug Discov. 2009;8:783–93.
Santiago JA, Potashkin JA. A network approach to clinical intervention in neurodegenerative diseases. Trends Mol Med. 2014;20:694–703.
Carro MS, Lim WK, Alvarez MJ, Bollo RJ, Zhao X, Snyder EY, Sulman EP, Anne SL, Doetsch F, Colman H, et al. The transcriptional network for mesenchymal transformation of brain tumours. Nature. 2010;463:318–25.
Fletcher MN, Castro MA, Wang X, de Santiago I, O'Reilly M, Chin SF, Rueda OM, Caldas C, Ponder BA, Markowetz F, Meyer KB. Master regulators of FGFR2 signalling and breast cancer risk. Nat Commun. 2013;4:2464.
Aubry S, Shin W, Crary JF, Lefort R, Qureshi YH, Lefebvre C, Califano A, Shelanski ML. Assembly and interrogation of Alzheimer’s disease genetic networks reveal novel regulators of progression. PLoS One. 2015;10:e0120352.
Pfaffenseller B, da Silva Magalhães PV, De Bastiani MA, Castro MA, Gallitano AL, Kapczinski F, Klamt F. Differential expression of transcriptional regulatory units in the prefrontal cortex of patients with bipolar disorder: potential role of early growth response gene 3. Transl Psychiatry. 2016;6:e805.
Trabzuni D, Ryten M, Walker R, Smith C, Imran S, Ramasamy A, Weale ME, Hardy J. Quality control parameters on a large dataset of regionally dissected human control brains for whole genome expression studies. J Neurochem. 2011;119:275–82.
Liang WS, Dunckley T, Beach TG, Grover A, Mastroeni D, Walker DG, Caselli RJ, Kukull WA, McKeel D, Morris JC, et al. Gene expression profiles in anatomically and functionally distinct regions of the normal aged human brain. Physiol Genomics. 2007;28:311–22.
Liang WS, Reiman EM, Valla J, Dunckley T, Beach TG, Grover A, Niedzielko TL, Schneider LE, Mastroeni D, Caselli R, et al. Alzheimer’s disease is associated with reduced expression of energy metabolism genes in posterior cingulate neurons. Proc Natl Acad Sci U S A. 2008;105:4441–6.
Miller JA, Woltjer RL, Goodenbour JM, Horvath S, Geschwind DH. Genes and pathways underlying regional and cell type changes in Alzheimer’s disease. Genome Med. 2013;5:48.
Hokama M, Oka S, Leon J, Ninomiya T, Honda H, Sasaki K, Iwaki T, Ohara T, Sasaki T, LaFerla FM, et al. Altered expression of diabetes-related genes in Alzheimer’s disease brains: the Hisayama study. Cereb Cortex. 2014;24:2476–88.
Berchtold NC, Coleman PD, Cribbs DH, Rogers J, Gillen DL, Cotman CW. Synaptic genes are extensively downregulated across multiple brain regions in normal human aging and Alzheimer’s disease. Neurobiol Aging. 2013;34:1653–61.
Margolin AA, Wang K, Lim WK, Kustagi M, Nemenman I, Califano A. Reverse engineering cellular networks. Nat Protoc. 2006;1:662–71.
Steuer R, Kurths J, Daub CO, Weise J, Selbig J. The mutual information: detecting and evaluating dependencies between variables. Bioinformatics. 2002;18(Suppl 2):S231–40.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; 2018. [https://www.r-project.org/].
Castro MA, Wang X, Fletcher MN, Meyer KB, Markowetz F. RedeR: R/Bioconductor package for representing modular structures, nested networks and multiple levels of hierarchical associations. Genome Biol. 2012;13:R29.
Wickham H. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York; 2009.
Castro MA, de Santiago I, Campbell TM, Vaughn C, Hickey TE, Ross E, Tilley WD, Markowetz F, Ponder BA, Meyer KB. Regulators of genetic risk of breast cancer identified by integrative network analysis. Nat Genet. 2016;48:12–21.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.
Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, Lerner J, Brunet JP, Subramanian A, Ross KN, et al. The connectivity map: using gene-expression signatures to connect small molecules, genes, and disease. Science. 2006;313:1929–35.
Duplan E, Sevalle J, Viotti J, Goiran T, Bauer C, Renbaum P, Levy-Lahad E, Gautier CA, Corti O, Leroudier N, et al. Parkin differently regulates presenilin-1 and presenilin-2 functions by direct control of their promoter transcription. J Mol Cell Biol. 2013;5:132–42.
Pearson AG, Curtis MA, Waldvogel HJ, Faull RL, Dragunow M. Activating transcription factor 2 expression in the adult human brain: association with both neurodegeneration and neurogenesis. Neuroscience. 2005;133:437–51.
Yamada T, Yoshiyama Y, Kawaguchi N. Expression of activating transcription factor-2 (ATF-2), one of the cyclic AMP response element (CRE) binding proteins, in Alzheimer disease and non-neurological brain tissues. Brain Res. 1997;749:329–34.
Zhang CW, Hang L, Yao TP, Lim KL. Parkin regulation and neurodegenerative disorders. Front Aging Neurosci. 2015;7:248.
Unschuld PG, Dachsel J, Darios F, Kohlmann A, Casademunt E, Lehmann-Horn K, Dichgans M, Ruberg M, Brice A, Gasser T, Lucking CB. Parkin modulates gene expression in control and ceramide-treated PC12 cells. Mol Biol Rep. 2006;33:13–32.
da Costa CA, Sunyach C, Giaime E, West A, Corti O, Brice A, Safe S, Abou-Sleiman PM, Wood NW, Takahashi H, et al. Transcriptional repression of p53 by parkin and impairment by mutations associated with autosomal recessive juvenile Parkinson's disease. Nat Cell Biol. 2009;11:1370–5.
Checler F. Processing of the beta-amyloid precursor protein and its regulation in Alzheimer’s disease. J Neurochem. 1995;65:1431–44.
Watson G, Ronai ZA, Lau E. ATF2, a paradigm of the multifaceted regulation of transcription factors in biology and disease. Pharmacol Res. 2017;119:347–57.
Liu H, Deng X, Shyu YJ, Li JJ, Taparowsky EJ, Hu CD. Mutual regulation of c-Jun and ATF2 by transcriptional activation and subcellular localization. EMBO J. 2006;25:1058–69.
Hayakawa J, Mittal S, Wang Y, Korkmaz KS, Adamson E, English C, Ohmichi M, McClelland M, Mercola D. Identification of promoters bound by c-Jun/ATF2 during rapid large-scale gene activation following genotoxic stress. Mol Cell. 2004;16:521–35.
Lopes FM, Schröder R, da Frota ML, Zanotto-Filho A, Müller CB, Pires AS, Meurer RT, Colpo GD, Gelain DP, Kapczinski F, et al. Comparison between proliferative and neuron-like SH-SY5Y cells as an in vitro model for Parkinson disease studies. Brain Res. 2010;1337:85–94.
Lau E, Ronai ZA. ATF2—at the crossroad of nuclear and cytosolic functions. J Cell Sci. 2012;125:2815–24.
Canter RG, Penney J, Tsai LH. The road to restoring neural circuits for the treatment of Alzheimer’s disease. Nature. 2016;539:187–96.
Yankner BA, Lu T, Loerch P. The aging brain. Annu Rev Pathol. 2008;3:41–66.
Aslam A, Mittal S, Koch F, Andrau JC, Winkler GS. The Ccr4-NOT deadenylase subunits CNOT7 and CNOT8 have overlapping roles and modulate cell proliferation. Mol Biol Cell. 2009;20:3840–50.
Bogdan JA, Adams-Burton C, Pedicord DL, Sukovich DA, Benfield PA, Corjay MH, Stoltenborg JK, Dicker IB. Human carbon catabolite repressor protein (CCR4)-associative factor 1: cloning, expression and characterization of its interaction with the B-cell translocation protein BTG1. Biochem J. 1998;336(Pt 2):471–81.
Sim DL, Chow VT. The novel human HUEL (C4orf1) gene maps to chromosome 4p12-p13 and encodes a nuclear protein containing the nuclear receptor interaction motif. Genomics. 1999;59:224–33.
Sim DL, Yeo WM, Chow VT. The novel human HUEL (C4orf1) protein shares homology with the DNA-binding domain of the XPA DNA repair protein and displays nuclear translocation in a cell cycle-dependent manner. Int J Biochem Cell Biol. 2002;34:487–504.
Chen YH, Kim JH, Stallcup MR. GAC63, a GRIP1-dependent nuclear receptor coactivator. Mol Cell Biol. 2005;25:5965–72.
Chen YH, Yang CK, Xia M, Ou CY, Stallcup MR. Role of GAC63 in transcriptional activation mediated by beta-catenin. Nucleic Acids Res. 2007;35:2084–92.
Shibanuma M, Kuroki T, Nose K. Isolation of a gene encoding a putative leucine zipper structure that is induced by transforming growth factor beta 1 and other growth factors. J Biol Chem. 1992;267:10219–24.
Kester HA, Blanchetot C, den Hertog J, van der Saag PT, van der Burg B. Transforming growth factor-beta-stimulated clone-22 is a member of a family of leucine zipper proteins that can homo- and heterodimerize and has transcriptional repressor activity. J Biol Chem. 1999;274:27439–47.
Yoon CH, Rho SB, Kim ST, Kho S, Park J, Jang IS, Woo S, Kim SS, Lee JH, Lee SH. Crucial role of TSC-22 in preventing the proteasomal degradation of p53 in cervical cancer. PLoS One. 2012;7:e42006.
Lahiri DK, Ge YW, Rogers JT, Sambamurti K, Greig NH, Maloney B. Taking down the unindicted co-conspirators of amyloid beta-peptide-mediated neuronal death: shared gene regulation of BACE1 and APP genes interacting with CREB, Fe65 and YY1 transcription factors. Curr Alzheimer Res. 2006;3:475–83.
Nowak K, Lange-Dohna C, Zeitschel U, Gunther A, Luscher B, Robitzki A, Perez-Polo R, Rossner S. The transcription factor yin Yang 1 is an activator of BACE1 expression. J Neurochem. 2006;96:1696–707.
Qin W, Jia L, Zhou A, Zuo X, Cheng Z, Wang F, Shi F, Jia J. The -980C/G polymorphism in APH-1A promoter confers risk of Alzheimer’s disease. Aging Cell. 2011;10:711–9.
Cho HJ, Kim SK, Jin SM, Hwang EM, Kim YS, Huh K, Mook-Jung I. IFN-gamma-induced BACE1 expression is mediated by activation of JAK2 and ERK1/2 signaling pathways and direct binding of STAT1 to BACE1 promoter in astrocytes. Glia. 2007;55:253–62.
Sastre M, Walter J, Gentleman SM. Interactions between APP secretases and inflammatory mediators. J Neuroinflammation. 2008;5:25.
Hsu WL, Ma YL, Hsieh DY, Liu YC, Lee EH. STAT1 negatively regulates spatial memory formation and mediates the memory-impairing effect of Abeta. Neuropsychopharmacology. 2014;39:746–58.
Burton TR, Dibrov A, Kashour T, Amara FM. Anti-apoptotic wild-type Alzheimer amyloid precursor protein signaling involves the p38 mitogen-activated protein kinase/MEF2 pathway. Brain Res Mol Brain Res. 2002;108:102–20.
Gonzalez P, Alvarez V, Menendez M, Lahoz CH, Martinez C, Corao AI, Calatayud MT, Pena J, Garcia-Castro M, Coto E. Myocyte enhancing factor-2A in Alzheimer’s disease: genetic analysis and association with MEF2A-polymorphisms. Neurosci Lett. 2007;411:47–51.
Kril JJ, Hodges J, Halliday G. Relationship between hippocampal volume and CA1 neuron loss in brains of humans with and without Alzheimer’s disease. Neurosci Lett. 2004;361:9–12.
Verkhratsky A, Marutle A, Rodríguez-Arellano JJ, Nordberg A. Glial asthenia and functional paralysis: a new perspective on neurodegeneration and Alzheimer’s disease. Neuroscientist. 2015;21:552–68.
Verkhratsky A, Olabarria M, Noristani HN, Yeh CY, Rodriguez JJ. Astrocytes in Alzheimer’s disease. Neurotherapeutics. 2010;7:399–412.
Bobinski M, Wegiel J, Tarnawski M, Reisberg B, de Leon MJ, Miller DC, Wisniewski HM. Relationships between regional neuronal loss and neurofibrillary changes in the hippocampal formation and duration and severity of Alzheimer disease. J Neuropathol Exp Neurol. 1997;56:414–20.
Bobinski M, Wegiel J, Wisniewski HM, Tarnawski M, Reisberg B, Mlodzik B, de Leon MJ, Miller DC. Atrophy of hippocampal formation subdivisions correlates with stage and duration of Alzheimer disease. Dementia. 1995;6:205–10.
Rössler M, Zarski R, Bohl J, Ohm TG. Stage-dependent and sector-specific neuronal loss in hippocampus during Alzheimer’s disease. Acta Neuropathol. 2002;103:363–9.
Blasko I, Stampfer-Kountchev M, Robatscher P, Veerhuis R, Eikelenboom P, Grubeck-Loebenstein B. How chronic inflammation can affect the brain and support the development of Alzheimer’s disease in old age: the role of microglia and astrocytes. Aging Cell. 2004;3:169–76.
Liu C, Cui G, Zhu M, Kang X, Guo H. Neuroinflammation in Alzheimer’s disease: chemokines produced by astrocytes and chemokine receptors. Int J Clin Exp Pathol. 2014;7:8342–55.
Juraskova B, Andrys C, Holmerova I, Solichova D, Hrnciarikova D, Vankova H, Vasatko T, Krejsek J. Transforming growth factor beta and soluble endoglin in the healthy senior and in Alzheimer’s disease patients. J Nutr Health Aging. 2010;14:758–61.
Diniz LP, Tortelli V, Matias I, Morgado J, Bérgamo Araujo AP, Melo HM, Seixas da Silva GS, Alves-Leon SV, de Souza JM, Ferreira ST, et al. Astrocyte transforming growth factor Beta 1 protects synapses against Aβ oligomers in Alzheimer’s disease model. J Neurosci. 2017;37:6797–809.
Cahoy JD, Emery B, Kaushal A, Foo LC, Zamanian JL, Christopherson KS, Xing Y, Lubischer JL, Krieg PA, Krupenko SA, et al. A transcriptome database for astrocytes, neurons, and oligodendrocytes: a new resource for understanding brain development and function. J Neurosci. 2008;28:264–78.
Azmi AS, Wang Z, Philip PA, Mohammad RM, Sarkar FH. Proof of concept: network and systems biology approaches aid in the discovery of potent anticancer drug combinations. Mol Cancer Ther. 2010;9:3137–44.
Nguyen TV, Yao M, Pike CJ. Flutamide and cyproterone acetate exert agonist effects: induction of androgen receptor-dependent neuroprotection. Endocrinology. 2007;148:2936–43.
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:870–80.
NIH: US National Institutes of Health - ClinicalTrials. gov. 2012.
Kazantsev AG, Thompson LM. Therapeutic application of histone deacetylase inhibitors for central nervous system disorders. Nat Rev Drug Discov. 2008;7:854–68.
Kornfeld K, Evason K. Effects of anticonvulsant drugs on life span. Arch Neurol. 2006;63:491–6.
Wildburger NC, Lin-Ye A, Baird MA, Lei D, Bao J. Neuroprotective effects of blockers for T-type calcium channels. Mol Neurodegener. 2009;4:44.
Gold B, Rodriguez WJ. Cefuroxime: mechanisms of action, antimicrobial activity, pharmacokinetics, clinical applications, adverse reactions and therapeutic indications. Pharmacotherapy. 1983;3:82–100.
Schindler AE. Present and future aspects of dydrogesterone in prevention or treatment of pregnancy disorders: an outlook. Horm Mol Biol Clin Investig. 2016;27:49–53.
Singh M, Su C. Progesterone and neuroprotection. Horm Behav. 2013;63:284–90.
Tikka T, Usenius T, Tenhunen M, Keinanen R, Koistinaho J. Tetracycline derivatives and ceftriaxone, a cephalosporin antibiotic, protect neurons against apoptosis induced by ionizing radiation. J Neurochem. 2001;78:1409–14.
Zumkehr J, Rodriguez-Ortiz CJ, Cheng D, Kieu Z, Wai T, Hawkins C, Kilian J, Lim SL, Medeiros R, Kitazawa M. Ceftriaxone ameliorates tau pathology and cognitive decline via restoration of glial glutamate transporter in a mouse model of Alzheimer’s disease. Neurobiol Aging. 2015;36:2260–71.
Bertoni JM, Weintraub ST. Competitive inhibition of human brain hexokinase by metrizamide and related compounds. J Neurochem. 1984;42:513–8.
Satoh J, Illes Z, Peterfalvi A, Tabunoki H, Rozsa C, Yamamura T. Aberrant transcriptional regulatory network in T cells of multiple sclerosis. Neurosci Lett. 2007;422:30–3.
Zhang B, Xia C, Lin Q, Huang J. Identification of key pathways and transcription factors related to Parkinson disease in genome wide. Mol Biol Rep. 2012;39:10881–7.
FK is a recipient of a fellowship award from the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq).
The present study was supported by the Brazilian funds CNPq/MS/SCTIE/DECIT research about Neurodegenerative Disease (466989/2014–8) and INCT-TM/CNPq/FAPESP (#465458/2014–9).
Availability of data and materials
Datasets used in this study can be accessed via the NCBI GEO portal (www.ncbi.nlm.nih.gov/geo/). Intermediate data and codes generated are available from the corresponding author on request. All data results generated during this study are included in this article and its Additional files.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Methodology flowchart. (A) Publicly available expression profile from Gene Expression Omnibus (GEO) were retrieved for normal brain hippocampus (GSE60862) and AD case versus control hippocampus (GSE5281, GSE29378, GSE36980, GSE48350). (B) Normal brain dataset was submitted to reverse engineering TF-centered transcription network reconstruction using ARACNe algorithm. Inferred healthy hippocampus regulatory units were then employed to query the master regulators of AD using master regulator analysis. Finally, the master regulator candidates were investigated for their state of activation using two-tail GSEA and possible repurposing drugs using connectivity maps. (PDF 5913 kb)
Table S1. Hippocampus transcription factor-centered network nodes and edges information. (XLSX 40 kb)
Table S2. Master regulators analysis results. (XLSX 23 kb)
Table S3. Master regulator candidates subnetwork nodes and edges information. (XLSX 12 kb)
Table S4. Connectivity map of master regulator candidates. (XLSX 11 kb)
Figure S2. Activation state of MR candidates in AD case-control studies and mouse neuron versus astrocyte data. Tile plot representing the MR candidate state of activation (two-tail gene set enrichment analysis) for the AD case-control (GSE5281, GSE29378, GSE36980, and GSE48350) and mouse neuron versus astrocyte (GSE9566) expression datasets. (PDF 24 kb)
About this article
Cite this article
Vargas, D.M., De Bastiani, M.A., Zimmer, E.R. et al. Alzheimer’s disease master regulators analysis: search for potential molecular targets and drug repositioning candidates. Alz Res Therapy 10, 59 (2018). https://doi.org/10.1186/s13195-018-0394-7