Wnt/β-catenin signaling modulates brain development and function and its deregulation underlies pathological changes occurring in neurodegenerative and neurodevelopmental disorders. Since one of the main effects of Wnt/β-catenin signaling is the modulation of target genes, in the present work we examined global transcriptional changes induced by short-term Wnt3a treatment (4 h) in primary cultures of rat hippocampal neurons. RNAseq experiments allowed the identification of 170 differentially expressed genes, including known Wnt/β-catenin target genes such as Notum, Axin2, and Lef1, as well as novel potential candidates Fam84a, Stk32a, and Itga9. Main biological processes enriched with differentially expressed genes included neural precursor (GO:0061364, p-adjusted = 2.5 × 10−7), forebrain development (GO:0030900, p-adjusted = 7.3 × 10−7), and stem cell differentiation (GO:0048863 p-adjusted = 7.3 × 10−7). Likewise, following activation of the signaling cascade, the expression of a significant number of genes with transcription factor activity (GO:0043565, p-adjusted = 4.1 × 10−6) was induced. We also studied molecular networks enriched upon Wnt3a activation and detected three highly significant expression modules involved in glycerolipid metabolic process (GO:0046486, p-adjusted = 4.5 × 10−19), learning or memory (GO:0007611, p-adjusted = 4.0 × 10−5), and neurotransmitter secretion (GO:0007269, p-adjusted = 5.3 × 10−12). Our results indicate that Wnt/β-catenin mediated transcription controls multiple biological processes related to neuronal structure and activity that are affected in synaptic dysfunction disorders.
Comisión Nacional de Investigación Científica y TecnológicaFondo Nacional de Desarrollo Científico y Tecnológico11403531. Introduction
The Wnt signaling cascade plays an essential role during embryogenesis and adult tissue homeostasis. Wnts are lipid modified secreted glycoproteins that signal through three major cellular pathways: the Planar Cell Polarity, the Wnt/Ca2+, and the Wnt/β-catenin signaling pathway, also referred to as the canonical Wnt signaling pathway [1–3]. The canonical cascade initiates with the binding of a Wnt ligand to Frizzled (FZD) receptors and LRP5/6 coreceptors located at the cellular membrane [4]. Wnt binding leads to the inhibition of the β-catenin destruction complex consisting of Axin, adenomatous polyposis coli (APC) [5], casein kinase 1 (CK1), and glycogen synthase kinase 3β (GSK3β) [6], which ultimately results in the stabilization of β-catenin protein in the cytosol and its subsequent nuclear translocation where it interacts with members of the T-cell factor/lymphoid enhancing factor (TCF/LEF) family of transcription factors to enhance transcription of Wnt/β-catenin target genes [3]. Conversely, in the absence of Wnt ligand activation, Axin and APC facilitate the sequential phosphorylation of β-catenin by CK1 and GSK3β [6] tagging this protein for ubiquitination and subsequent proteasome mediated degradation [7].
Throughout mammalian brain development the activity of the Wnt cascade is spatially confined to specialized regions such as the olfactory bulb, frontal cortex, hippocampal formation, and the cerebellum [8–11]. In these brain domains Wnt/β-catenin signaling participates in diverse biological processes including neurogenesis [12], axonal remodeling and patterning [13, 14], and development and maturation of functional synapses within the CNS [15–20]. Indeed, Wnt1, Wnt3a, Wnt7a, and Wnt8 are ligands known to activate Wnt/β-catenin signaling and are involved in brain development and synaptogenesis [21, 22]. Wnt3a is essential in early development of hippocampal structures and participates in the establishment of long term potentiation events [23, 24]. Wnt7a and Wnt8a have also been shown to regulate excitatory synaptic formation [17, 25]. Furthermore a recent study suggests that LRP6, Wnt/β-catenin signaling coreceptor, is critical for the development of functional synapses in vivo [25]. Therefore, given its multiple roles in synaptic function and brain homeostasis, Wnt/β-catenin signaling is a functional and positional candidate to understand complex prevalent neurological conditions in the human population.
At the presynaptic region, canonical Wnt ligands such as Wnt7a and Wnt3a enhance the clustering and recycling of synaptic vesicles (SVs) in primary cultures of rat hippocampal neurons [26]. Consistently, loss of Wnt7a function inhibits SVs clustering, an effect that is mimicked by loss-of-function of Dishevelled 1 (DVL1) signaling downstream of Wnt ligands [19]. Interestingly, the Wnt7a/Dvl1 double mutant shows defects in spine morphogenesis and excitatory synaptic neurotransmission [17], which parallels behavioral abnormalities with a disrupted presynaptic assembly and excitatory/inhibitory balance. Wnt/β-catenin signaling also seems to trigger neurotransmitter release and SV trafficking by modulating SV-associated phosphoproteins. While Wnt7a and Wnt3a enhance the clustering [27] and phosphorylation [28] of synapsin 1 at the synaptic button prior to neurotransmitter release, Dvl1 is involved in neurotransmitter release through direct binding to synaptotagmin I in differentiated neurons [29].
Experience driven plasticity is highly dependent on proper synaptic transmission and is mainly modulated by Ca2+ related pathways. In this regard, Wnt noncanonical and canonical pathways have been extensively related to Ca2+ homeostasis and signaling [19, 28, 30, 31]. For instance, ligands such as Wnt3a [28], Wnt5a [30], and Wnt7a [19] have all been shown to increase Ca2+ influxes to stimulate excitatory synaptic strength in hippocampal neurons or in peripheral nerves to alter pain sensitivity [32]. Other mechanisms modulating the activity of the synaptic terminal involve the function of cell adhesion proteins, most notably transsynaptic cadherin-β-catenin interactions that have an essential function during the recruitment and clustering of SVs to synapses [33–37]. Importantly, the effect on synaptic function, plus the biological control of neuronal polarity and axon outgrowth/navigation to their final synaptic targets, is mainly accomplished through rapid posttranslational changes affecting the cytoskeletal machinery [38, 39].
While Wnt/β-catenin ligands have a profound effect on the modulation of synaptic terminals; however the transcriptional program elicited by the signaling cascade in neurons has received little attention. Therefore, given that the integration of differentially expressed genes into biological networks provides a wider vision for the transcriptional landscape of a particular cell-type and considering that others and we have previously observed that Wnt/β-catenin target gene expression is rapidly induced (2–6 h) upon activation of the signaling cascade [40–44], here we use genome-wide transcriptional data to identify target genes and main biological processes directed by transient Wnt3a exposure (4 h) in primary cultures of rat hippocampal neurons. We expect that understanding of the genetic program maintained by β-catenin/TCF-LEF complexes in hippocampal cells will provide molecular mechanisms to couple activity, structure, and function of synaptic networks and circuits in this important brain region, which is severely affected in neurodegenerative and neurodevelopmental disorders.
2. Materials and Methods2.1. Primary Cultures of Rat Hippocampal Neurons
Hippocampal neurons were dissociated and maintained as described before [28]. Briefly, neurons were taken from 18-day pregnant Sprague-Dawley rats and maintained for 14 days in vitro (DIV) on 12-well culture plates (500,000 cells/well) coated with poly-L-lysine (Sigma) and supplemented with neurobasal/B27 media (Gibco). The culture was placed on a shelf in a 37°C humidified CO2 incubator, and the medium was changed every 2 days.
2.2. Wnt3a Purification
Wnt3a is a specific Wnt/β-catenin signaling agonist that can be efficiently recovered from conditioned medium from Wnt3a-secreting L-cells (ATCC, Rockville, MD; see also the Wnt homepage at http://web.stanford.edu/group/nusselab/cgi-bin/wnt/). Wnt3a purification was carried out as previously described [28, 57, 58]. The presence of the Wnt3a protein was detected with an anti-Wnt3a antibody (R & D Systems, Minneapolis, MN). Purity was analyzed by SDS-PAGE (8%), stained with Coomassie Blue G250, and analyzed through densitometry by using software ImageJ [59].
2.3. RNAseq Experiments
Rat hippocampal neurons were transiently exposed with Wnt3a for 4 h and then total RNA was extracted (3 control and 3 Wnt3a samples matched by experiment). RNA integrity was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). RNA integrity numbers (RIN) obtained were above 7.5 and suitable for sequencing. RNA was processed using the Illumina TruSeq Stranded mRNA kit to generate 280 bp size paired-end libraries. Libraries were sequenced using the Illumina HiSeq 2000 benchtop sequencer. Raw reads (average depth of 100.62 M 2 × 80) were aligned to the Rnor6 build for rat genome reference, using Bioconductor R Package (v3.1; [60]). Differential expression was statistically analyzed with DESeq v1.20 following author recommendations [61]. Heterogeneity across the samples was examined visualizing reads distribution plots before and after DESeq normalization. Pearson’s correlation analysis showed r2 values above 0.93 between all samples, within the standards (i.e., 0.92 and 0.98) suggested by ENCODE [62]. The complete pipeline of our study is presented in Supplementary Figure 1 in Supplementary Material available online at http://dx.doi.org/10.1155/2016/4672841. Summary statistics for the sequencing process and alignment are provided in Supplementary Tables 1 and 2, respectively. Likewise, read normalization results are presented in Supplementary Figure 2.
2.4. Gene Ontology (GO) Enrichment Analysis
We tested the gene ontology (GO; http://geneontology.org/) [63] structure and annotations using the package Ontologizer [64], considering categories with less than 500 members to avoid associations to major categories that are less informative (i.e., signaling) and excluding the ones “Inferred by Electronic Annotation” (IEA), from “Reviewed Computational Analysis” (RCA) and with “No biological Data available” (ND), which are characterized by a high rate of false positives [65, 66]. We used the parent-child-union algorithm to call for overrepresentation adjusting the p values with Benjamini-Hochberg multiple testing correction, to avoid false positives [67]. We considered GO terms significantly overrepresented with an adjusted p value below 1 × 10−6 for further analyses.
2.5. Functional Protein Association Network (FPAN)
To assay for known interactions between the differentially expressed genes we retrieved high confidence functional interactions from the STRING 10.0 database (search tool for the retrieval of interacting genes/proteins; http://string-db.org/) [68] which contains curated interactions from different evidence sources (i.e., genomic context, coexpression, and curated literature). We kept only “highest confidence” interactions (i.e., the interactions with a combined score >0.9 provided by STRING). The final FPAN was composed of 9,443 nodes (rat genes) and 309,728 nonredundant edges (interactions). All network and subnetworks were visualized with Cytoscape software [69].
2.6. Module Search
All gene-wise p values obtained with DESeq were introduced in to the FPAN as a floating-point attribute. A fold change greater than 1.25 (25%) was used as a threshold for inclusion (n = 2,892 genes). Module search was carried out with the Cytoscape JActiveModules plugin [70] with a gene overlap threshold of 20% following the method reported previously [66, 71]. Briefly, the program searches for significant highly connected subnetworks or modules enriched with expression information (p value). Starting from one random node (Monte Carlo procedure) modules grow in comparison with an expected background distribution created by the software from 10,000 internal randomizations, obtaining a specific standard deviation score (S score). Modules with S>3 (3 standard deviation above the mean of randomized scores) and with a gene number between 10 and 50 were considered significant [72]. To acquire a mean S score and standard deviation (SD) for each resulting module the search was performed 10 times. Finally, the same procedure was conducted with permuted p values over the entire genes present in the FPAN (permuted analysis). Statistical differences between permuted and real (Wnt3a) analyses were assessed through one-sided Student’s t-test.
3. Results3.1. Differentially Expressed Wnt/β-Catenin Target Genes and Ontological Categories
To observe a direct transcriptional effect of Wnt/β-catenin signaling in primary cultures of hippocampal neurons and thus avoid noise due to activation of secondary unrelated pathways [40–44], we compared the whole transcriptome of cells transiently incubated with purified Wnt3a (4 h, n=3) with cells from the same batch under control conditions (Wnt3a vehicle; n=3; see Material and Methods). We identified 170 differentially expressed genes with a nominal p value of <0.05 (Figure 1). Raw read counts per sample, in ascending order according to their respective fold change (FC), are shown in Figures 1(a) and 1(b). Smear plots between treatment and control conditions revealed that differentially expressed genes (red dots) displayed an inverse relationship between the observed FC and the mean of normalized counts (Figure 1(c)). Overall, we found that 1% of the total rat genes included in our analyses (170 out of 16,252; Figure 1(d)) showed a differential response upon Wnt3a stimulation and that 88% of these genes were upregulated (149 genes) (Supplementary Table 3). Functional characterization of differentially expressed genes revealed that transcription factors (18.1%), proteins with nucleic acid binding activity (14.8%), and signaling related molecules (8.7%) represent main protein categories activated upon Wnt3a stimulation (Figure 1(e)).
Transcriptome analysis, differentially expressed genes, and their functional classification. (a) Rat genes total amount of raw reads is displayed for each of the samples sequenced (C1: Control 1, W1: Wnt1, C2: Control 2, W2: Wnt2, C3: Control 3, and W3: Wnt3). Dark to light blue color scale denotes a higher to lower number of counts per gene, respectively. (b) Fold change observed for the cluster analysis of the 170 differentially expressed rat genes in red to green scale for up- to downregulation, respectively. (c) Smear plot. log 2 fold change (y axis) values versus normalized read counts per gene (x axis) between treatment and control conditions. A total of 170 confident transcript calls over the threshold of 0.05 p value (red dots) are shown. (d) Comparison of the 170 genes found differentially expressed and the total genes analyzed. (e) For all genes over threshold, we determined their protein class with the Panther software (v11.0).
We observed that the strongest upregulated gene in Wnt3a treated neurons was Fam84a (Gene ID: 313969; nominal p value = 2.2 × 10−6; FC = 2.4; Table 1). Likewise, among the top 20 upregulated genes we found 11 Wnt/β-catenin target genes (Table 1 and Supplementary Table 3), including the palmitoleoyl-protein carboxylesterase Notum (ID: 303743; nominal p value = 1.5 × 10−5; FC = 14.1), Lef1 (ID: 161452; nominal p value = 1.7 × 10−5; FC = 2.1), Axin2 (ID: 29134; nominal p value = 1.2 × 10−4; FC = 3.5), and the chemokine (C-X-C motif) ligand 3 Cxcl3 (ID: 171551; nominal p value = 5.2 × 10−4; FC = 4.6). Interestingly, novel Wnt/β-catenin targets genes ranked between known targets and included the serine/threonine kinase 32A Stk32a (ID: 364858; nominal p value = 1.3 × 10−3; FC = 1.8), the integrin subunit alpha 9 Itga9 (ID: 586004; nominal p value = 1.7 × 10−3; FC = 2.0), and the transmembrane protein 72 TMEM72 (ID: 362424; nominal p value = 2.0 × 10−3; FC = 1.7), among others (Table 1). These results indicate that after transient treatment with Wnt3a we can readily identify nascent RNA messages derived from activating the transcriptional program controlled by Wnt/β-catenin signaling in these hippocampal neurons. As to noncanonical Wnt signaling components being expressed in primary neurons as a consequence of Wnt3a treatment, we intersected the 170 nominally differentially expressed genes with the reactome pathway “beta catenin independent Wnt signaling” including 143 genes belonging to PCP and/or Ca2+ noncanonical pathways (DOI: 10.3180/REACT_172694.1). Only Lef1, Tcf7, and Prkg2 were found in both gene sets (Supplementary Figure 3a). While Lef1 and Tcf7 are widely acknowledged as major canonical components, Prkg2 is recognized only as a functional modulator of the canonical pathway (Table 2) through the inhibition of GSK3β activity [48]. Thus we conclude that Wnt3a did not enhance the expression of noncanonical components.
Top 20 nominal significant Wnt3a-upregulated genes in hippocampal neurons.
Symbol ID
FC
p value
Wnt pathway
References
Fam84a
2.4
2.2 × 10−6
—
—
Notum
14.1
1.5 × 10−5
Wnt target gene
[45]
Lef1
2.1
1.7 × 10−5
Wnt target gene
[46]
Axin2
3.5
1.2 × 10−4
Wnt target gene
[47]
Prkg2
2.7
1.4 × 10−4
Functional
[48]
Cxcl3
4.6
5.2 × 10−4
Wnt target gene
[49]
Ahr
2.1
1.0 × 10−3
Wnt target gene
[50]
Stk32a
1.8
1.3 × 10−3
—
—
Itga9
2.0
1.7 × 10−3
—
—
Tnfresf19
2.4
1.8 × 10−3
Wnt target gene
[51]
Tmem72
1.7
2.0 × 10−3
—
—
Fras1
1.9
2.2 × 10−3
—
—
Fam167a
1.8
2.2 × 10−3
—
—
Gata2
2.8
2.3 × 10−3
Wnt target gene
[52]
Id2
1.3
3.0 × 10−3
Wnt target gene
[53]
Msx2
2.1
3.9 × 10−3
Wnt target gene
[54]
Sp5
10.3
4.2 × 10−3
Wnt target gene
[55]
Cd83
2.1
4.4 × 10−3
—
—
Cgnl1
1.5
4.9 × 10−3
—
—
Hunk
1.5
5.3 × 10−3
Wnt target gene
[56]
Functional: known to modulate Wnt pathway transcriptional activity. FC: fold change.
GO categories overrepresented in the 170 differentially expressed genes.
Category
GO ID
GO Name
GIP
GIS
p-adjusted (BH)
BP
GO:0061351
Neural precursor cell proliferation
122
14
1.6 × 10−10
GO:0030900
Forebrain development
266
18
2.2 × 10−09
GO:0048863
Stem cell differentiation
267
18
2.3 × 10−09
GO:0007420
Brain development
418
22
1.0 × 10−06
GO:0048732
Gland development
289
18
1.6 × 10−06
GO:0048864
Stem cell development
225
16
1.6 × 10−06
GO:0021872
Forebrain generation of neurons
58
9
3.3 × 10−06
GO:0021953
Central nervous system neuron differentiation
138
12
7.0 × 10−06
GO:0045165
Cell fate commitment
173
13
9.6 × 10−06
MF
GO:0043565
Sequence-specific DNA binding
452
22
4.1 × 10−06
BP: biological processes; CC: cellular components; MF: molecular function; GO ID: gene ontology term ID; GIP: genes in population; GIS: genes in study; BH: Benjamini-Hochberg correction multiple testing correction.
At a global level, we found that differentially expressed genes after Wnt3a treatment define gene ontology (GO) categories that have essential roles during brain development and homeostasis (Table 2). We observed that major biological process (BP) categories enriched with differentially expressed genes included neural precursor cell proliferation (GO:0061351; adjusted p value = 1.6 × 10−10), forebrain development (GO:0030900; adjusted p value = 2.2 × 10−9), and stem cell differentiation (GO:0048863; adjusted p value = 2.3 × 10−9). Similarly, sequence-specific DNA binding activity (GO:0043565; adjusted p value = 4.1 × 10−6) was the main molecular function (MF) category overrepresented in our analysis. All 170 differential expressed genes and their functional relationships are shown in Figure 2. A high degree of overlap was observed between previous GO terms with 37 genes (7 of which were among top 20 upregulated; Figure 2, hexagonal nodes) responsible for the enrichment given by 89 ontological associations (Figure 2, color edges). Additionally, we detected 11 high confidence interactions connecting 21 out of the 170 genes using the entire rat functional protein association network (FPAN) extracted from STRING database (see Materials and Methods).
Functional relationships of Wnt differentially expressed genes. Nominal Wnt-induced differentially expressed genes are presented (n=170) alongside their mayor functional interactions. Node color represents the observed fold change (FC) upon Wnt3a treatment in a gradient from green (downregulation) to red (upregulation) values (i.e., green: FC < 1; red FC > 1; white: FC = 1). Similarly, node size is proportional to the -log10p value assigned by DESeq to the gene. Central boxes denote the top 5 overrepresented GO categories and their respective gene associations (color edges) which are presented for neural precursor (GO:0061364), forebrain development (GO:0030900), stem cell differentiation (GO:0048863), brain development (GO:0007420), and gland development (GO:0048732) categories. Direct functional interactions extracted from STRING are shown with black edges. Hexagon-shaped nodes show top 25 genes differentially expressed. Diamond-shaped node highlight top 1 differentially expressed gene (Fam68a). Black and grey node borders denote top 20 Wnt target or functional evidence, respectively.
3.2. Molecular Interaction Networks Activated by Wnt/β-Catenin Signaling
While differential gene expression analysis allows for the identification of disease altered genes in specific tissues by analyzing one gene at the time, this approach does not take into account the functional relationships between genes [73]. Therefore, since the convergence or aggregation of additional expression signals in a set of genes can become significant in a network context, we searched for highly connected subnetworks or modules displaying significant changes in expression upon Wnt signaling activation. We performed a module search using 2,892 upregulated genes (DESeq FC > 1.25) on the entire rat FPAN composed of a total of 9,443 genes with at least 1 high confidence interaction defined by STRING database (edge score above 900, see Materials and Methods) [68]. We observed a significantly higher number of modules (p value = 5.00 × 10−4) using actual Wnt3a transcriptomic data (average = 17.00; SD = 0.982) in comparison with modules arising from chance using permuted data (average = 11.00, SD = 0.987) (Figure 3(a)). While the average score of the first fourth modules derived from Wnt3a data was significantly higher compared to their permuted counterpart (M1–M4; Figure 3(b)), the average score for M4 was below the highest random score and therefore M4 was not considered for further analysis. Altogether, these results show that the number and structure of modules remained consistent across iterations and that they could not be reached by chance, supporting the idea that Wnt3a-derived subnetworks M1, M2, and M3 are indeed biologically meaningful. Components and interactions of M1, M2, and M3 subnetworks are provided in Figure 4. GO analysis further revealed that M1 was highly enriched with genes belonging to the glycerolipid metabolic process category (Table 3; GO:0046486; adjusted p value = 4.5 × 10−19) and with lipase activity (GO:0016298; adjusted p value = 5.2 × 10−9); M2 was mainly involved in learning or memory processes (Table 3; GO:0007611; adjusted p value = 4.0 × 10−5) with genes related to cell adhesion molecule binding activity (GO:0050839; adjusted p value = 1.1 × 10−6); and M3 was overrepresented with genes belonging to neurotransmitter secretion process (Table 3; GO:0007269; adjusted p value = 5.3 × 10−12) and with syntaxin-1 binding activity (GO:0017075; adjusted p value = 1.2 × 10−12). The complete list of the 93 genes belonging to M1, M2, and M3 subnetworks is presented in Supplementary Table 4 with their observed FC and p value.
Top 3 GO categories overrepresented in Wnt3a-induced M1, M2, and M3 subnetworks.
Module
Category
GO ID
GO Name
GIP
GIS
p-adjusted (BH)
M1
BP
GO:0046486
Glycerolipid metabolic process
193
18
4.5 × 10−19
GO:0006644
Phospholipid metabolic process
166
17
1.1 × 10−18
GO:0006650
Glycerophospholipid metabolic process
134
16
1.2 × 10−18
MF
GO:0016298
Lipase activity
59
8
5.2 × 10−09
GO:0016788
Hydrolase activity, acting on ester bonds
430
14
1.2 × 10−08
GO:0004620
Phospholipase activity
47
7
1.7 × 10−08
M2
BP
GO:0007611
Learning or memory
181
7
4.0 × 10−05
GO:0050890
Cognition
200
7
4.4 × 10−05
CC
GO:0043256
Laminin complex
8
4
2.0 × 10−07
GO:0005610
Laminin-5 complex
4
3
2.9 × 10−06
GO:0005605
Basal lamina
18
4
2.9 × 10−06
MF
GO:0050839
Cell adhesion molecule binding
131
7
1.1 × 10−06
GO:0001948
Glycoprotein binding
74
5
1.7 × 10−05
GO:0005178
Integrin binding
76
5
1.7 × 10−05
M3
BP
GO:0007269
Neurotransmitter secretion
82
10
5.3 × 10−12
GO:0006887
Exocytosis
176
12
5.3 × 10−12
GO:0016079
Synaptic vesicle exocytosis
37
8
1.2 × 10−11
CC
GO:0044456
Synapse part
373
14
2.6 × 10−11
GO:0043679
Axon terminus
131
10
7.7 × 10−11
GO:0044306
Neuron projection terminus
139
10
9.4 × 10−11
MF
GO:0017075
Syntaxin-1 binding
16
7
1.2 × 10−12
GO:0000149
SNARE binding
69
8
8.0 × 10−10
GO:0019905
Syntaxin binding
46
7
1.7 × 10−09
BP: biological processes; CC: cellular components; MF: molecular function; GO ID: gene ontology term ID; GIP: genes in population; GIS: genes in study; BH: Benjamini-Hochberg correction multiple testing correction.
Module search results over the rat FPAN. (a) The number of significant modules (size between 10 and 50 genes with a SD score >3) obtained with nominal data (Wnt3a, dark grey) and with the same values permuted across the FPAN (permuted data, light grey). (b) Individual SD score comparison of the top 10 modules obtained in the corresponding module searches presented in (a). Asterisk denotes the Student t-test p value significance (p∗<0.05; p∗∗<0.001; p∗∗∗<0.0001).
Molecular networks induced by Wnt3a treatment. Gene composition (nodes) and interactions (edges) are shown for modules M1 (upper right, green background), M2 (middle bottom, light blue background), and M3 (upper right, orange background). Node color represents the observed fold change (FC) upon Wnt3a treatment in a gradient from green (downregulation) to red (upregulation) values (i.e., green, FC < 1; red, FC > 1; white, FC = 1). Similarly, node size is proportional to the -log10p value assigned to the gene behavior. For each module the main GO category associated is shown: M1 = glycerolipid metabolism process (GO:0046486), M2 = learning and memory (GO:0007611), and M3 = neurotransmitter secretion (GO:0007269). Module structure: M1 is composed of 44 genes with 132 interactions; M1 is composed of 25 genes with 32 interactions; and M1 is composed of 30 genes with 62 interactions (see also Supplementary Table 4).
4. Discussion
The transcriptional program directed by Wnt/β-catenin signaling has been examined mainly in cancer models and has allowed for the identification of novel target genes and a better understanding of the oncogenic properties of this pathway [74–76]. Conversely, much remains to be elucidated about the transcriptional program of the cascade in the central nervous system. In this context, initial attempts to discover novel Wnt/β-catenin target genes functioning in the brain relied on in silico analyses of proximal gene promoters enriched with TCF/LEF binding sites [77, 78]. More recently, bioinformatics and transcriptomic approaches have been used in thalamic neurons [79], or human neural progenitor cells [41], to identify Wnt/β-catenin target genes involved in neuronal excitation or neurodegenerative diseases, including Alzheimer’s disease. In this regard, multiple lines of evidence support a functional role for Wnt/β-catenin signaling in prevalent neurological disorders related to synaptic dysfunction including autism, Alzheimer’s disease, or epilepsy [80–83]. Nevertheless, although these studies have resulted in hundreds of potential novel candidates, the expression of those target genes shows minimal overlap (less than 2%) in different cell backgrounds (Supplementary Figure 3b), likely due to cell-type specific transcriptional and/or translational mechanisms.
Ultimately the goal of Wnt/β-catenin signaling in neuronal cells is to modulate changes in gene expression that are manifested in diverse cellular processes, including neurogenesis, axonal pathfinding, dendritic development, synaptic formation, and plasticity [39, 84]. These changes are likely accomplished by transcription factors that either interact or are activated downstream of β-catenin transcriptional complexes. Interestingly, we found that a large proportion (18.1%, n=28) of Wnt/β-catenin targets in hippocampal neurons are genes coding for transcription factors (Figure 1). Indeed, within the top 20 Wnt3a upregulated genes we observed 6 transcription factors: Lef1, Ahr, Gata2, Id2, Msx2, and Sp5, that are target genes or functional partners with components of the signaling cascade and that play a role in synaptic development. For instance TCF/LEF transcription factors are instrumental in Wnt/β-catenin target gene expression [85], which include target genes participating in synaptic development and function [79]. Similarly, Ahr participates in the expression and membrane presentation of NMDA receptors [86]. Gata2 overexpression in hippocampal neurons leads to a decrease in synaptic spines density and depressive behavior in rats [87]. Id family genes including Id2 are upregulated in Rett syndrome [88], a neurodevelopmental disorder which features altered synaptic plasticity [89]; Id2 is also asymmetrically expressed in human embryonic brain hemispheres suggesting a role for this transcription factor in cortical specialization [90]. Additionally, several of these proteins also modulate Wnt pathway activation in other tissues [91–94]. While further experimental validation including proteomics and TCF/LEF transcription factor occupancy at a genome-wide level (ChIP-seq) is needed to confirm these potential candidates as Wnt/β-catenin targets in neuronal and nonneuronal cells, we note that none of the early Wnt3a-upregulated genes seems to be the result of downstream transcription factors or secondary transcriptional waves (Supplementary Figure 3c).
GO analysis results from differential expression data revealed that Wnt/β-catenin signaling enhances transcription of genes involved in brain development, particularly forebrain development, stem cell differentiation, and neural precursor cell proliferation (Figure 2). In this regard, a role for Wnt/β-catenin components in forebrain development has been clearly established [95] and recently this issue has received considerable attention mainly associated with the onset of autism [96]. Similarly, Wnt signaling is well known for enhancing cell renewal of neural stem cells and is a key modulator of neurogenesis [97]. Notably, Axin2, a classic Wnt target gene [47] and β-catenin destruction complex scaffold protein [98], has been recently shown to control the switch of intermediate progenitors from proliferative to differentiate status in the developing cerebral cortex [99]. Indeed, enhancing Axin expression in neuronal progenitors leads to an enlarged neocortex, an excess of excitatory synapses and autistic like behaviors [100].
Network based approaches have been widely used in the study of neurological diseases where usually the number of genetic markers that exceeds the significance threshold is very small; thus most of the identified markers are neglected [73]. The integration of these below-threshold markers into biologically significant networks has successfully identified novel modules of interaction in multiple sclerosis, Alzheimer’s disease, autism, and other disorders [66, 71, 96, 101]. Our network analysis of Wnt3a differentially expressed pathways extends the notion that Wnt/β-catenin signaling has an essential function in neurotransmitter secretion, learning, and memory (Figure 4, modules M2 and M3). We also found a highly significant module comprising components of glycerolipid metabolism that is differentially expressed in neurons exposed to Wnt3a (module M1). M1 is enriched with several lipid-modifying enzymes families, such as PI(3)Ks, PLAs, PLCs, and PLDs whose products are readily found in synaptic membranes enhancing vesicle docking and release [102]. For instance, PLAs produce arachidonic acid, a precursor which will be further processed by COX2, a known Wnt/β-catenin target gene [103], during the generation of several eicosanoids which mediate the inflammatory response [104]. Arachidonic acid released from the postsynaptic terminal is known to potentiate synaptic transmission by inhibiting presynaptic potassium channels [105]. Mutations in PLAs have been described in both Alzheimer’s disease and autism [104, 106]. Likewise, the production of phosphatidic acid by PLDs modulates several aspects of Wnt/β-catenin signaling in cancer and is recognized as a critical regulator of cell proliferation and tumorigenesis [107].
The M2 subnetwork, involved in learning or memory, is enriched in extracellular matrix components such as integrins and laminins. First, Itga9 (integrin subunit alpha 9) is among the 10 top upregulated genes reported in this study. Integrins are differentially expressed in specific regions in the adult brain [108] and they interact with Reln (Reelin) to activate cortical lamination [109]. Interestingly, the Reelin pathway, which crosstalks with Wnt/β-catenin signaling throughout brain development [110], has been associated with Alzheimer’s disease [111] and autism [112]. Second, the enrichment of laminin subunits is particularly interesting since laminins are necessary for synaptic ultrastructure [113] and have been proven to prevent beta amyloid aggregation [114]. Third, neuronal cell adhesion molecules Neuroligin (Nlgn) and Neurexin (Nrxn) are essential for synaptic structure and function in the synaptic cleft and Nlgn3 and Nrxn1 have been associated with autism [115]. Finally, Wnt3a-derived data also allowed the identification of M3 subnetwork related to neurotransmitter secretion enriched with transcription factors belonging to the ATF/cAMP response element binding proteins (CREB) family of proteins, whose function is involved in synaptic plasticity and memory [116, 117]. Similarly, M3 includes several genes whose products are involved in membrane fusion, synaptic exocytosis, and presynaptic dynamics [118], including SNARES (Vamp2, Vamp3, and Stx1A), STX binding proteins (Stxbp1 and Stxbp5), and Rab signaling molecules (Rab3a and Rims1), which may account for synaptic effects induced by the crosstalk between canonical Wnt pathway and calcium signaling components and that may be important in synaptic dysfunction disorders.
5. Conclusion
Transcriptomics and network analyses are highly useful for identifying novel targets that can be used to better understand cellular changes during the onset or progression of neurodegenerative and neurodevelopmental disorders. Our results provide novel insights into the early transcriptional program and molecular networking directed by Wnt/β-catenin signaling in hippocampal neurons and warrant further investigation.
Ethical Approval
All procedures involving experimentation on animals were approved by the Bioethical Committee of Universidad Andres Bello and were conducted in accordance with the guidelines of the National Fund for Scientific and Technological Research (FONDECYT; Chile).
Disclosure
Current address of Miguel E. Ávila is Instituto de Ciencias Naturales, Universidad de las Américas, Santiago, Chile.
Competing Interests
The authors declare no conflict of interests.
Authors’ Contributions
Eduardo Pérez-Palma, Víctor Andrade, Giorgia D. Ugarte, and Giancarlo V. De Ferrari conceived and designed the experiments. Eduardo Pérez-Palma, Víctor Andrade, Bernabé I. Bustos, Camilo Villaman, Matías A. Medina, Miguel E. Ávila, and Giorgia D. Ugarte performed experiments or analyzed the data. Eduardo Pérez-Palma, Víctor Andrade, Mario O. Caracci, and Giancarlo V. De Ferrari wrote the paper. All authors read and approved the final manuscript. Eduardo Pérez-Palma and Víctor Andrade contributed equally to this work.
Acknowledgments
This study was supported by CONICYT regular FONDECYT 1140353 grant to Giancarlo V. De Ferrari.
MoonR. T.KohnA. D.De FerrariG. V.KaykasA.WNT and β-catenin signalling: diseases and therapies20045969170110.1038/nrg14272-s2.0-4344584730StaalF. J. T.LuisT. C.TiemessenM. M.WNT signalling in the immune system: WNT is spreading its wings2008885815931861788510.1038/nri23602-s2.0-4874910193718617885CleversH.NusseR.Wnt/β-catenin signaling and disease20121496119212052268224310.1016/j.cell.2012.05.0122-s2.0-8486198605322682243HeX.SemenovM.TamaiK.ZengX.LDL receptor-related proteins 5 and 6 in Wnt/β-catenin signaling: arrows point the way200413181663167710.1242/dev.011172-s2.0-2342430345LipinaT. V.Kaidanovich-BeilinO.PatelS.WangM.ClapcoteS. J.LiuF.WoodgettJ. R.RoderJ. C.Genetic and pharmacological evidence for schizophrenia-related Disc1 interaction with GSK-320116532342482068711110.1002/syn.208392-s2.0-7865065764920687111HartM. J.De Los SantosR.AlbertI. N.RubinfeldB.PolakisP.Downregulation of β-catenin by human Axin and its association with the APC tumor suppressor, β-catenin and GSK3β1998810573581960164110.1016/S0960-9822(98)70226-X2-s2.0-00324929549601641AberleH.BauerA.StappertJ.KispertA.KemlerR.β‐Catenin is a target for the ubiquitin–proteasome pathway199716133797380410.1093/emboj/16.13.37972-s2.0-0030978351GroveE. A.ToleS.LimonJ.YipL.-W.RagsdaleC. W.The hem of the embryonic cerebral cortex is defined by the expression of multiple Wnt genes and is compromised in Gli3-deficient mice1998125122315232595841302-s2.0-00318699269584130HouartC.CaneparoL.HeisenbergC.-P.BarthK. A.Take-UchiM.WilsonS. W.Establishment of the telencephalon during gastrulation by local antagonism of Wnt signaling20023522552651216074410.1016/S0896-6273(02)00751-12-s2.0-003713045912160744MarettoS.CordenonsiM.DupontS.BraghettaP.BroccoliV.HassanA. B.VolpinD.BressanG. M.PiccoloS.Mapping Wnt/β-catenin signaling during mouse development and in colorectal tumors20031006329933041262675710.1073/pnas.04345901002-s2.0-003745293212626757ZhouC.-J.ZhaoC.PleasureS. J.Wnt signaling mutants have decreased dentate granule cell production and radial glial scaffolding abnormalities200424112112610.1523/jneurosci.4071-03.20042-s2.0-0842322748LieD.-C.ColamarinoS. A.SongH.-J.DésiréL.MiraH.ConsiglioA.LeinE. S.JessbergerS.LansfordH.DearieA. R.GageF. H.Wnt signalling regulates adult hippocampal neurogenesis20054377063137013751625196710.1038/nature041082-s2.0-2584443211916251967HollisE. R.ZouY.Expression of the Wnt signaling system in central nervous system axon guidance and regeneration20125510.3389/fnmol.2012.00005TakeichiM.AbeK.Synaptic contact dynamics controlled by cadherin and catenins20051542162211581737810.1016/j.tcb.2005.02.0022-s2.0-1684436246615817378De Ferrari andG. V.InestrosaN. C.Wnt signaling function in Alzheimer's disease200033111210.1016/s0165-0173(00)00021-72-s2.0-0033850796De FerrariG. V.MoonR. T.The ups and downs of Wnt signaling in prevalent neurological disorders200625577545755310.1038/sj.onc.1210064Ahmad-AnnuarA.CianiL.SimeonidisI.HerrerosJ.FredjN. B.RossoS. B.HallA.BrickleyS.SalinasP. C.Signaling across the synapse: a role for Wnt and Dishevelled in presynaptic assembly and neurotransmitter release2006174112713910.1083/jcb.2005110542-s2.0-33745606159InestrosaN. C.ArenasE.Emerging roles of Wnts in the adult nervous system201011277862001095010.1038/nrn27552-s2.0-7554908726120010950CianiL.BoyleK. A.DickinsE.SahoresM.AnaneD.LopesD. M.GibbA. J.SalinasP. C.Wnt7a signaling promotes dendritic spine growth and synaptic strength through Ca 2+/Calmodulin-dependent protein kinase II20111082610732107372167030210.1073/pnas.10181321082-s2.0-7996061469921670302De FerrariG. V.AvilaM. E.MedinaM. A.Pérez-PalmaE.BustosB. I.AlarcónM. A.Wnt/β-catenin signaling in Alzheimer's disease201413574575410.2174/18715273126661312231139002-s2.0-84905591087ItasakiN.JonesC. M.MercurioS.RoweA.DomingosP. M.SmithJ. C.KrumlaufR.Wise, a context-dependent activator and inhibitor of Wnt signalling200313018429543051290044710.1242/dev.006742-s2.0-014172548312900447CaricasoleA.FerraroT.IacovelliL.BarlettaE.CarusoA.MelchiorriD.TerstappenG. C.NicolettiF.Functional characterization of WNT7A signaling in PC12 cells: interaction with a FZD5∙LRP6 receptor complex and modulation by Dickkopf proteins200327839370243703110.1074/jbc.m3001912002-s2.0-0141484619LeeS. M. K.ToleS.GroveE.McMahonA. P.A local Wnt-3a signal is required for development of the mammalian hippocampus20001273457467106311672-s2.0-003396640510631167ChenJ.ParkC. S.TangS.-J.Activity-dependent synaptic Wnt release regulates hippocampal long term potentiation200628117119101191610.1074/jbc.m5119202002-s2.0-33744958556SharmaK.ChoiS.-Y.ZhangY.NielandT. J. F.LongS.LiM.HuganirR. L.High-throughput genetic screen for synaptogenic factors: identification of lrp6 as critical for excitatory synapse development2013551330134110.1016/j.celrep.2013.11.0082-s2.0-84890229121CerpaW.GodoyJ. A.AlfaroI.FaríasG. G.MetcalfeM. J.FuentealbaR.BonanscoC.InestrosaN. C.Wnt-7a modulates the synaptic vesicle cycle and synaptic transmission in hippocampal neurons20082839591859271809670510.1074/jbc.M7059432002-s2.0-4194914220218096705HallA. C.LucasF. R.SalinasP. C.Axonal remodeling and synaptic differentiation in the cerebellum is regulated by WNT-7a signaling200010055255351072199010.1016/S0092-8674(00)80689-32-s2.0-003459892110721990AvilaM. E.SepúlvedaF. J.BurgosC. F.Moraga-CidG.ParodiJ.MoonR. T.AguayoL. G.OpazoC.De FerrariG. V.Canonical Wnt3a modulates intracellular calcium and enhances excitatory neurotransmission in hippocampal neurons201028524189391894710.1074/jbc.m110.1030282-s2.0-77953295258KishidaS.HamaoK.InoueM.HasegawaM.MatsuuraY.MikoshibaK.FukudaM.KikuchiA.Dvl regulates endo- and exocytotic processes through binding to synaptotagmin200712149611721265410.1111/j.1365-2443.2006.01030.x2-s2.0-3384577941517212654Varela-NallarL.AlfaroI. E.SerranoF. G.ParodiJ.InestrosaN. C.Wingless-type family member 5A (Wnt-5a) stimulates synaptic differentiation and function of glutamatergic synapses20101074921164211692108463610.1073/pnas.10100111072-s2.0-7865050642521084636CerpaW.GambrillA.InestrosaN. C.BarriaA.Regulation of NMDA-receptor synaptic transmission by Wnt signaling20113126946694712171561110.1523/JNEUROSCI.6311-10.20112-s2.0-7995996079121715611SimonettiM.AgarwalN.StösserS.BaliK. K.KaraulanovE.KambleR.PospisilovaB.KurejovaM.BirchmeierW.NiehrsC.HeppenstallP.KunerR.Wnt-fzd signaling sensitizes peripheral sensory neurons via distinct noncanonical pathways20148311041212499195610.1016/j.neuron.2014.05.0372-s2.0-8490359675324991956BamjiS. X.ShimazuK.KimesN.HuelskenJ.BirchmeierW.LuB.ReichardtL. F.Role of β-catenin in synaptic vesicle localization and presynaptic assembly20034047197311462257710.1016/S0896-6273(03)00718-92-s2.0-024257331614622577IwaiY.HirotaY.OzakiK.OkanoH.TakeichiM.UemuraT.DN-cadherin is required for spatial arrangement of nerve terminals and ultrastructural organization of synapses20021933753881190621010.1006/mcne.2001.10812-s2.0-003621452811906210JünglingK.EulenburgV.MooreR.KemlerR.LessmannV.GottmannK.N-cadherin transsynaptically regulates short-term plasticity at glutamatergic synapses in embryonic stem cell-derived neurons200626266968697810.1523/jneurosci.1013-06.20062-s2.0-33745959647SunY.AigaM.YoshidaE.HumbertP. O.BamjiS. X.Scribble interacts with β-catenin to localize synaptic vesicles to synapses20092014339034001945819710.1091/mbc.E08-12-11722-s2.0-6774910353219458197TogashiH.AbeK.MizoguchiA.TakaokaK.ChisakaO.TakeichiM.Cadherin regulates dendritic spine morphogenesis200235177891212361010.1016/S0896-6273(02)00748-12-s2.0-003701444612123610OlivaC. A.VargasJ. Y.InestrosaN. C.Wnt signaling: role in LTP, neural networks and memory201312378680010.1016/j.arr.2013.03.0062-s2.0-84892522551SalinasP. C.Wnt signaling in the vertebrate central nervous system: from axon guidance to synaptic function2012422-s2.0-84862238182GujralT. S.MacBeathG.A system-wide investigation of the dynamics of Wnt signaling reveals novel phases of transcriptional regulation201054e1002410.1371/journal.pone.0010024WexlerE. M.RosenE.LuD.OsbornG. E.MartinE.RaybouldH.GeschwindD. H.Genome-wide analysis of a Wnt1-regulated transcriptional network implicates neurodegenerative pathways20114193, article ra6510.1126/scisignal.20022822-s2.0-80053523654NuñezF.BravoS.CruzatF.MontecinoM.De FerrariG. V.CampbellM.Wnt/β-catenin signaling enhances cyclooxygenase-2 (COX2) Transcriptional activity in gastric cancer cells201164e1856210.1371/journal.pone.0018562UgarteG. D.VargasM. F.MedinaM. A.LeónP.NecuñirD.ElorzaA. A.GutiérrezS. E.MoonR. T.LoyolaA.De FerrariG. V.Wnt signaling induces transcription, spatial proximity, and translocation of fusion gene partners in human hematopoietic cells201512615178517892633377610.1182/blood-2015-04-6384942-s2.0-8494365410626333776MedinaM. A.UgarteG. D.VargasM. F.AvilaM. E.NecuñirD.ElorzaA. A.GutiérrezS. E.De FerrariG. V.Alternative RUNX1 promoter regulation by Wnt/β-catenin signaling in leukemia cells and human hematopoietic progenitors201623171460146710.1002/jcp.25258ChangM. V.ChangJ. L.GangopadhyayA.ShearerA.CadiganK. M.Activation of wingless targets requires bipartite recognition of DNA by TCF200818231877188110.1016/j.cub.2008.10.047HovanesK.LiT. W. H.MunguiaJ. E.TruongT.MilovanovicT.Lawrence MarshJ.HolcombeR. F.WatermanM. L.β-Catenin-sensitive isoforms of lymphoid enhancer factor-1 are selectively expressed in colon cancer200128153571132627610.1038/882642-s2.0-003503352111326276JhoE.-H.ZhangT.DomonC.JooC.-K.FreundJ.-N.CostantiniF.Wnt/β-catenin/Tcf signaling induces the transcription of Axin2, a negative regulator of the signaling pathway2002224117211831180980810.1128/MCB.22.4.1172-1183.20022-s2.0-003614820811809808KawasakiY.KugimiyaF.ChikudaH.KamekuraS.IkedaT.KawamuraN.SaitoT.ShinodaY.HigashikawaA.YanoF.OgasawaraT.OgataN.HoshiK.HofmannF.WoodgettJ. R.NakamuraK.ChungU.-I.KawaguchiH.Phosphorylation of GSK-3β by cGMP-dependent protein kinase II promotes hypertrophic differentiation of murine chondrocytes20081187250625151855119510.1172/JC1352432-s2.0-4674913483518551195SegditsasS.SieberO.DeheragodaM.EastP.RowanA.JefferyR.NyeE.ClarkS.Spencer-DeneB.StampG.PoulsomR.SuraweeraN.SilverA.IlyasM.TomlinsonI.Putative direct and indirect Wnt targets identified through consistent gene expression changes in APC-mutant intestinal adenomas from humans and mice20081724386438751878285110.1093/hmg/ddn2862-s2.0-5704916240518782851Gerbal-ChaloinS.DuméA.-S.BriolottiP.KlieberS.RauletE.DuretC.FabreJ.-M.RamosJ.MaurelP.Daujat-ChavanieuM.The WNT/β-catenin pathway is a transcriptional regulator of CYP2E1, CYP1A2, and aryl hydrocarbon receptor gene expression in primary human hepatocytes201486662463410.1124/mol.114.0947972-s2.0-84912041758ButtittaL.TanakaT. S.ChenA. E.KoM. S. H.FanC.-M.Microarray analysis of somitogenesis reveals novel targets of different WNT signaling pathways in the somitic mesoderm20032581911041278168510.1016/S0012-1606(03)00116-72-s2.0-003768476912781685GustafsonB.SmithU.Activation of canonical wingless-type MMTV integration site family (Wnt) signaling in mature adipocytes increases β-catenin levels and leads to cell dedifferentiation and insulin resistance201028518140311404110.1074/jbc.m110.1028552-s2.0-77951519413RockmanS. P.CurrieS. A.CiavarellaM.VincanE.DowC.ThomasR. J. S.PhillipsW. A.Id2 is a target of the β-catenin/T cell factor pathway in colon carcinoma200127648451134511910.1074/jbc.m1077422002-s2.0-0035976964ZhaiY.IuraA.YeasminS.WieseA. B.WuR.FengY.FearonE. R.ChoK. R.MSX2 is an oncogenic downstream target of activated WNT signaling in ovarian endometrioid adenocarcinoma20113040415241622149930010.1038/onc.2011.1232-s2.0-8005361450421499300FujimuraN.VacikT.MachonO.VlcekC.ScalabrinS.SpethM.DiepD.KraussS.KozmikZ.Wnt-mediated down-regulation of Sp1 target genes by a transcriptional repressor Sp520072822122512371709053410.1074/jbc.M6058512002-s2.0-3384762006117090534ReedK. R.KorobkoI. V.NinkinaN.KorobkoE. V.HopkinsB. R.PlattJ. L.BuchmanV.ClarkeA. R.Hunk/Mak-v is a negative regulator of intestinal cell proliferation2015151, article no. 11010.1186/s12885-015-1087-22-s2.0-84928239089KishidaS.YamamotoH.KikuchiA.Wnt-3a and Dvl induce neurite retraction by activating rho-associated kinase20042410448745011512186610.1128/MCB.24.10.4487-4501.20042-s2.0-294259279815121866WillertK.BrownJ. D.DanenbergE.DuncanA. W.WeissmanI. L.ReyaT.YatesJ. R.IIINusseR.Wnt proteins are lipid-modified and can act as stem cell growth factors200342369384484521271745110.1038/nature016112-s2.0-003773772612717451SchneiderC. A.RasbandW. S.EliceiriK. W.NIH Image to ImageJ: 25 years of image analysis2012976716752293083410.1038/nmeth.20892-s2.0-8486320584922930834HuberW.CareyV. J.GentlemanR.AndersS.CarlsonM.CarvalhoB. S.BravoH. C.DavisS.GattoL.GirkeT.GottardoR.HahneF.HansenK. D.IrizarryR. A.LawrenceM.LoveM. I.MacDonaldJ.ObenchainV.OleśA. K.PagèsH.ReyesA.ShannonP.SmythG. K.TenenbaumD.WaldronL.MorganM.Orchestrating high-throughput genomic analysis with Bioconductor201512211512110.1038/nmeth.3252AndersS.HuberW.Differential expression analysis for sequence count data20101110, article R10610.1186/gb-2010-11-10-r106ConsortiumE. P.An integrated encyclopedia of DNA elements in the human genome20124897414577410.1038/nature11247Gene Ontology ConsortiumThe Gene Ontology (GO) database and informatics resource200432D258D26110.1093/nar/gkh036GrossmannS.BauerS.RobinsonP. N.VingronM.Improved detection of overrepresentation of Gene-Ontology annotations with parent-child analysis20072322302430311784839810.1093/bioinformatics/btm4402-s2.0-3644894614617848398Warde-FarleyD.DonaldsonS. L.ComesO.ZuberiK.BadrawiR.ChaoP.FranzM.GrouiosC.KaziF.LopesC. T.MaitlandA.MostafaviS.MontojoJ.ShaoQ.WrightG.BaderG. D.MorrisQ.The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function201038W214W22010.1093/nar/gkq537Pérez-PalmaE.BustosB. I.VillamánC. F.AlarcónM. A.AvilaM. E.UgarteG. D.ReyesA. E.OpazoC.De FerrariG. V.Overrepresentation of glutamate signaling in Alzheimer's disease: network-based pathway enrichment using meta-analysis of genome-wide association studies201494e9541310.1371/journal.pone.00954132-s2.0-84899676806WestfallP. H.YoungS. S.1993Wiley-InterscienceSzklarczykD.FranceschiniA.KuhnM.SimonovicM.RothA.MinguezP.DoerksT.StarkM.MullerJ.BorkP.JensenL. J.MeringC. v.The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored201039supplement 1D561D56810.1093/nar/gkq973ClineM. S.SmootM.CeramiE.Integration of biological networks and gene expression data using Cytoscape20072102366238210.1038/nprot.2007.324IdekerT.OzierO.SchwikowskiB.SiegelA. F.Discovering regulatory and signalling circuits in molecular interaction networks200218supplement 1S233S24010.1093/bioinformatics/18.suppl_1.s2332-s2.0-0000801240BaranziniS. E.GalweyN. W.WangJ.KhankhanianP.LindbergR.PelletierD.WuW.UitdehaagB. M. J.KapposL.PolmanC. H.MatthewsP. M.HauserS. L.GibsonR. A.OksenbergJ. R.BarnesM. R.Pathway and network-based analysis of genome-wide association studies in multiple sclerosis20091811207820901928667110.1093/hmg/ddp1202-s2.0-6554908428919286671BaranziniS. E.GalweyN. W.WangJ.KhankhanianP.LindbergR.PelletierD.WuW.UitdehaagB. M.KapposL.PolmanC. H.MatthewsP. M.HauserS. L.GibsonR. A.OksenbergJ. R.BarnesM. R.Pathway and network-based analysis of genome-wide association studies in multiple sclerosis200918112078209010.1093/hmg/ddp120ParikshakN. N.GandalM. J.GeschwindD. H.Systems biology and gene networks in neurodevelopmental and neurodegenerative disorders201516844145810.1038/nrg39342-s2.0-84937520584WillertJ.EppingM.PollackJ. R.BrownP. O.NusseR.The role of serine proteases and serine protease inhibitors in the migration of gonadotropin-releasing hormone neurons20022, article no. 11710.1186/1471-213x-2-12-s2.0-0142254155TaneyhillL.PennicaD.Identification of Wnt responsive genes using a murine mammary epithelial cell line model system20044, article no. 110.1186/1471-213x-4-12-s2.0-2642571109HatzisP.Van Der FlierL. G.Van DrielM. A.GuryevV.NielsenF.DenissovS.NijmanI. J.KosterJ.SantoE. E.WelborenW.VersteegR.CuppenE.Van De WeteringM.CleversH.StunnenbergH. G.Genome-wide pattern of TCF7L2/TCF4 chromatin occupancy in colorectal cancer cells2008288273227441826800610.1128/MCB.02175-072-s2.0-4214918357618268006ArrázolaM. S.Varela-NallarL.ColombresM.ToledoE. M.CruzatF.PavezL.AssarR.AravenaA.GonzálezM.MontecinoM.MaassA.MartínezS.InestrosaN. C.Calcium/calmodulin-dependent protein kinase type IV is a target gene of the Wnt/β-catenin signaling pathway2009221365866710.1002/jcp.219022-s2.0-70449699789HödarC.AssarR.ColombresM.AravenaA.PavezL.GonzálezM.MartínezS.InestrosaN. C.MaassA.Genome-wide identification of new Wnt/β-catenin target genes in the human genome using CART method2010111, article 34810.1186/1471-2164-11-3482-s2.0-77954279005WisniewskaM. B.NagalskiA.DabrowskiM.MisztalK.KuznickiJ.Novel β-catenin target genes identified in thalamic neurons encode modulators of neuronal excitability2012131, article 63510.1186/1471-2164-13-6352-s2.0-84869037450De FerrariG. V.MoonR. T.The ups and downs of Wnt signaling in prevalent neurological disorders20062557754575531714329910.1038/sj.onc.12100642-s2.0-3384535034617143299CaracciM. O.ÁvilaM. E.De FerrariG. V.Synaptic Wnt/GSK3β signaling hub in autism2016201610960375110.1155/2016/9603751De FerrariG. V.AvilaM. E.MedinaM. A.Pérez-PalmaE.BustosB. I.AlarcónM. A.Wnt/β-catenin signaling in Alzheimer's disease201413574575410.2174/18715273126661312231139002-s2.0-84905591087HuangC.FuX.-H.ZhouD.LiJ.-M.The role of Wnt/β-catenin signaling pathway in disrupted hippocampal neurogenesis of temporal lobe epilepsy: a potential therapeutic target?20154071319133210.1007/s11064-015-1614-12-s2.0-84938211518InestrosaN. C.Varela-NallarL.Wnt signaling in the nervous system and in Alzheimer's disease20146164742454915710.1093/jmcb/mjt0512-s2.0-8489682886724549157CadiganK. M.WatermanM. L.TCF/LEFs and Wnt signaling in the nucleus2012411a00790610.1101/cshperspect.a0079062-s2.0-84868540323LinC.-H.ChenC.-C.ChouC.-M.WangC.-Y.HungC.-C.ChenJ. Y.ChangH.-W.ChenY.-C.YehG. C.LeeY.-H.Knockdown of the aryl hydrocarbon receptor attenuates excitotoxicity and enhances NMDA-induced BDNF expression in cortical neurons200911137777891971205510.1111/j.1471-4159.2009.06364.x2-s2.0-7034990038919712055ChoiM.WangS. E.KoS. Y.KangH. J.ChaeS. Y.LeeS. H.KimY.-S.DumanR. S.SonH.Overexpression of human GATA-1 and GATA-2 interferes with spine formation and produces depressive behavior in rats2014910e10925310.1371/journal.pone.01092532-s2.0-84908587621PeddadaS.YasuiD. H.LaSalleJ. M.Inhibitors of differentiation (ID1, ID2, ID3 and ID4) genes are neuronal targets of MeCP2 that are elevated in Rett syndrome20061512200320141668243510.1093/hmg/ddl1242-s2.0-3374560927616682435WengS.-M.McLeodF.BaileyM. E. S.CobbS. R.Synaptic plasticity deficits in an experimental model of rett syndrome: long-term potentiation saturation and its pharmacological reversal201118031432110.1016/j.neuroscience.2011.01.0612-s2.0-79953167470SunT.PatoineC.Abu-KhalilA.VisvaderJ.SumE.CherryT. J.OrkinS. H.GeschwindD. H.WalshC. A.Early asymmetry of gene transcription in embryonic human left and right cerebral cortex200530857291794179810.1126/science.11103242-s2.0-20544435515WeidingerG.ThorpeC. J.Wuennenberg-StapletonK.NgaiJ.MoonR. T.The Sp1-related transcription factors sp5 and sp5-like act downstream of Wnt/β-catenin signaling in mesoderm and neuroectoderm patterning200515648950010.1016/j.cub.2005.01.0412-s2.0-15744365467ProcházkováJ.KabátkováM.BryjaV.UmannováL.BernatíkO.KozubíkA.MachalaM.VondrácěkJ.The interplay of the aryl hydrocarbon receptor and β-catenin alters both AhR-dependent transcription and wnt/β-catenin signaling in liver progenitors201112223493602160219110.1093/toxsci/kfr1292-s2.0-8005170206521602191MimotoM. S.KwonS.GreenY. S.GoldmanD.ChristianJ. L.GATA2 regulates Wnt signaling to promote primitive red blood cell fate2015407111110.1016/j.ydbio.2015.08.0122-s2.0-84947866159DongH.-J.JangG.-B.LeeH.-Y.ParkS.KimJ.NamJ.HongI.The Wnt/β-catenin signaling/Id2 cascade mediates the effects of hypoxia on the hierarchy of colorectal-cancer stem cells201662296610.1038/srep22966Harrison-UyS. J.PleasureS. J.Wnt signaling and forebrain development201247a00809410.1101/cshperspect.a0080942-s2.0-84868231671KrummN.O'RoakB. J.ShendureJ.EichlerE. E.A de novo convergence of autism genetics and molecular neuroscience2014372951052438778910.1016/j.tins.2013.11.0052-s2.0-8489370871324387789KalaniM. Y. S.CheshierS. H.CordB. J.BababeygyS. R.VogelH.WeissmanI. L.PalmerT. D.NusseR.Wnt-mediated self-renewal of neural stem/progenitor cells20081054416970169751895754510.1073/pnas.08086161052-s2.0-5594912167918957545KikuchiA.Roles of axin in the Wnt signalling pathway199911117777881061728010.1016/S0898-6568(99)00054-62-s2.0-003275658410617280FangW.-Q.ChenW.-W.FuA. Y.IpN.Axin directs the amplification and differentiation of intermediate progenitors in the developing cerebral cortex20137946656792397259610.1016/j.neuron.2013.06.0172-s2.0-8488267865223972596FangW.-Q.ChenW.-W.JiangL.LiuK.YungW.-H.FuA. K. Y.IpN. Y.Overproduction of upper-layer neurons in the neocortex leads to autism-like features in mice2014951635164310.1016/j.celrep.2014.11.0032-s2.0-84915783802JohnsonM. R.ShkuraK.LangleyS. R.Delahaye-DuriezA.SrivastavaP.HillW. D.RackhamO. J.DaviesG.HarrisS. E.Moreno-MoralA.RotivalM.SpeedD.PetrovskiS.KatzA.HaywardC.PorteousD. J.SmithB. H.PadmanabhanS.HockingL. J.StarrJ. M.LiewaldD. C.ViscontiA.FalchiM.BottoloL.RossettiT.DanisB.MazzuferiM.FoerchP.GroteA.HelmstaedterC.BeckerA. J.KaminskiR. M.DearyI. J.PetrettoE.Systems genetics identifies a convergent gene network for cognition and neurodevelopmental disease201619222323210.1038/nn.4205RohrboughJ.BroadieK.Lipid regulation of the synaptic vesicle cycle2005621391501568521910.1038/nrn16082-s2.0-1324427540315685219NuñezF.BravoS.CruzatF.MontecinoM.de FerrariG. V.Wnt/β-catenin signaling enhances cyclooxygenase-2 (COX2) transcriptional activity in gastric cancer cells201164e1856210.1371/journal.pone.00185622-s2.0-79953758074TamijiJ.CrawfordD. A.The neurobiology of lipid metabolism in autism spectrum disorders2011182981122134637710.1159/0003231892-s2.0-7995241959321346377CartaM.LanoreF.RebolaN.SzaboZ.Da SilvaS. V.LourençoJ.VerraesA.NadlerA.SchultzC.BlanchetC.MulleC.Membrane lipids tune synaptic transmission by direct modulation of presynaptic potassium channels20148147877992448608610.1016/j.neuron.2013.12.0282-s2.0-8489668992924486086Sanchez-MejiaR. O.MuckeL.Phospholipase A2 and arachidonic acid in Alzheimer's disease20101801878479010.1016/j.bbalip.2010.05.0132-s2.0-77953279191KangD. W.ChoiK.-Y.MinD. S.Phospholipase D meets Wnt signaling: a new target for cancer therapy201171229329710.1158/0008-5472.can-10-24632-s2.0-78751480479PinkstaffJ. K.DetterichJ.LynchG.GallC.Integrin subunit gene expression is regionally differentiated in adult brain1999195154115562-s2.0-0033103707DulabonL.OlsonE. C.TaglientiM. G.EisenhuthS.McGrathB.WalshC. A.KreidbergJ. A.AntonE. S.Reelin binds α3β1 integrin and inhibits neuronal migration200027133441093932910.1016/S0896-6273(00)00007-62-s2.0-003362427910939329ReinerO.SapirT.Similarities and differences between the Wnt and reelin pathways in the forming brain2005311–311713410.1385/mn:31:1-3:1172-s2.0-24644487688HerzJ.ChenY.Reelin, lipoprotein receptors and synaptic plasticity20067118508591705381010.1038/nrn20092-s2.0-3375033356917053810ZhangH.LiuX.ZhangC.MundoE.MacciardiF.GraysonD. R.GuidottiA. R.HoldenJ. J. A.Reelin gene alleles and susceptibility to autism spectrum disorders200279101210171239995610.1038/sj.mp.40011242-s2.0-003642885212399956EglesC.ClaudepierreT.ManglapusM. K.ChampliaudM.-F.BrunkenW. J.HunterD. D.Laminins containing the β2 chain modulate the precise organization of CNS synapses20073432882981718970110.1016/j.mcn.2006.11.0042-s2.0-3384727215417189701MorganC.InestrosaN. C.Interactions of laminin with the amyloid β peptide. Implications for Alzheimer's disease200134559760110.1590/s0100-879x20010005000062-s2.0-0035344128SüdhofT. C.Neuroligins and neurexins link synaptic function to cognitive disease200845572159039111892351210.1038/nature074562-s2.0-5404909194118923512LiJ.SutterC.ParkerD. S.BlauwkampT.FangM.CadiganK. M.CBP/p300 are bimodal regulators of Wnt signaling20072692284229410.1038/sj.emboj.76016672-s2.0-34247615104PasiniS.CoronaC.LiuJ.GreeneL. A.ShelanskiM. L.Specific downregulation of hippocampal ATF4 reveals a necessary role in synaptic plasticity and memory201511218319110.1016/j.celrep.2015.03.0252-s2.0-84929293304SudhofT. C.RizoJ.Synaptic vesicle exocytosis2011312a00563710.1101/cshperspect.a005637