# Early Transcriptional Changes Induced by Wnt/β-Catenin Signaling in Hippocampal Neurons

Authors: Eduardo Pérez-Palma; Víctor Andrade; Mario O. Caracci; Bernabé I. Bustos; Camilo Villaman; Matías A. Medina; Miguel E. Ávila; Giorgia D. Ugarte; Giancarlo V. De Ferrari
Journal: Neural Plasticity
 (2016)

---

## Abstract

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.


---

## Body

## 1. 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 Methods

### 2.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-wisep 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.


## 2.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-wisep 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. Results

### 3.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)).Figure 1
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).
(a) (b) (c) (d) (e)We observed that the strongest upregulated gene in Wnt3a treated neurons was Fam84a (Gene ID: 313969; nominalp 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.Table 1
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.Table 2
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 (Table2). 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).Figure 2
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 -log⁡10  p 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.Table 3
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.Figure 3
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 Studentt-test p value significance (p∗<0.05; p∗∗<0.001; p∗∗∗<0.0001).
(a) (b)Figure 4
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-log⁡10  p 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).


## 3.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)).Figure 1
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).
(a) (b) (c) (d) (e)We observed that the strongest upregulated gene in Wnt3a treated neurons was Fam84a (Gene ID: 313969; nominalp 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.Table 1
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.Table 2
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 (Table2). 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).Figure 2
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 -log⁡10  p 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.Table 3
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.Figure 3
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 Studentt-test p value significance (p∗<0.05; p∗∗<0.001; p∗∗∗<0.0001).
(a) (b)Figure 4
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-log⁡10  p 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.



---

*Source: 4672841-2016-12-27.xml*