Novel β-catenin target genes identified in thalamic neurons encode modulators of neuronal excitability

Background LEF1/TCF transcription factors and their activator β-catenin are effectors of the canonical Wnt pathway. Although Wnt/β-catenin signaling has been implicated in neurodegenerative and psychiatric disorders, its possible role in the adult brain remains enigmatic. To address this issue, we sought to identify the genetic program activated by β-catenin in neurons. We recently showed that β-catenin accumulates specifically in thalamic neurons where it activates Cacna1g gene expression. In the present study, we combined bioinformatics and experimental approaches to find new β-catenin targets in the adult thalamus. Results We first selected the genes with at least two conserved LEF/TCF motifs within the regulatory elements. The resulting list of 428 putative LEF1/TCF targets was significantly enriched in known Wnt targets, validating our approach. Functional annotation of the presumed targets also revealed a group of 41 genes, heretofore not associated with Wnt pathway activity, that encode proteins involved in neuronal signal transmission. Using custom polymerase chain reaction arrays, we profiled the expression of these genes in the rat forebrain. We found that nine of the analyzed genes were highly expressed in the thalamus compared with the cortex and hippocampus. Removal of nuclear β-catenin from thalamic neurons in vitro by introducing its negative regulator Axin2 reduced the expression of six of the nine genes. Immunoprecipitation of chromatin from the brain tissues confirmed the interaction between β-catenin and some of the predicted LEF1/TCF motifs. The results of these experiments validated four genes as authentic and direct targets of β-catenin: Gabra3 for the receptor of GABA neurotransmitter, Calb2 for the Ca2+-binding protein calretinin, and the Cacna1g and Kcna6 genes for voltage-gated ion channels. Two other genes from the latter cluster, Cacna2d2 and Kcnh8, appeared to be regulated by β-catenin, although the binding of β-catenin to the regulatory sequences of these genes could not be confirmed. Conclusions In the thalamus, β-catenin regulates the expression of a novel group of genes that encode proteins involved in neuronal excitation. This implies that the transcriptional activity of β-catenin is necessary for the proper excitability of thalamic neurons, may influence activity in the thalamocortical circuit, and may contribute to thalamic pathologies.

Background β-catenin is an armadillo family protein that serves as a gene expression regulator in canonical Wnt signaling, in addition to its function in cell adhesion [1,2]. The canonical Wnt pathway involves the inhibition of GSK3α/β (i.e., a negative regulator of β-catenin), dissociation of the β-catenin destruction complex that contains APC and Axin proteins, and subsequent accumulation of the cytoplasmic pool of β-catenin that can translocate into the nucleus and activate LEF1/TCF transcription factors [3,4]. The crucial role of this pathway in cell differentiation and proliferation is well established, and mutations in its components lead to severe malformations in developing embryos and cancer in adults [5].
The nuclear mediators of canonical Wnt signaling include the transcription factors LEF1, TCF7 (also known as TCF1), TCF7L1 (also known as TCF3), and TCF7L2 (also known as TCF4), which are members of the high mobility group (HMG) family [6][7][8]. The HMG DNA binding domain of LEF1/TCFs recognizes the WWCAAAG consensus sequence. The N-terminus of LEF1/TCF recruits β-catenin, which does not bind to DNA itself but has a strong transactivation domain. The remaining regulatory domains and alternatively spliced C-termini of the proteins vary between these members and provide a molecular basis for the diversity and context dependence of LEF1/TCFs function. Lef1/Tcf genes are highly expressed during mouse embryogenesis, but their transcripts are virtually undetectable postnatally [9,10], except for intestinal and hair follicle stem cells [11,12]. Surprisingly, high levels of Lef1 and Tcf7l2 expression have been observed in the thalamus in the adult brain [13][14][15][16].
Most of the approximately 100 known β-catenin targets are involved in development and cancer. In differentiated neurons, only a few genes have been shown to be activated by β-catenin; several are already known Wnt target genes that are not neuron-specific [35][36][37]. As we recently demonstrated, a new target, Cacna1G, encodes voltage-gated Ca 2+ channels [16]. In this work, we identified in silico and experimentally validated a novel group of β-catenin-LEF1/TCF targets in thalamic neurons, consisting of genes that encode proteins that are important for neuronal function, including voltage-and ligand-gated ion channels and the Ca 2+ -binding protein calretinin. This indicates a role for β-catenin and LEF1/ TCF transcription factors in the maintenance of neuronal excitability.

Results
Putative LEF1/TCF target genes identified in silico are enriched in known Wnt targets To find putative β-catenin target genes, we screened for LEF1/TCF family motifs in rat and human conserved regions within 10 kb upstream and downstream of transcription start sites (TSSs). We relied on the fact that enhancers are enriched in proximity to TSSs [39], and the clustering of TCF7L2-bound regions occurs within 10 kb of TSSs in colorectal cancer [40]. Using established bioinformatics tools with their default parameters, we identified 2,871 genes with at least one rat-to-human conserved LEF1/TCF motif and 851 genes that contain at least two LEF1/TCF motifs in the same conserved noncoding sequences (CNSs). These two lists were compared with experimentally established Wnt target genes cataloged at http://www.stanford.edu/group/nusselab/ cgi-bin/wnt/target_genes (accessed January 11, 2012; Figure 1). The degree of overlap between our lists and the list on the Wnt homepage shows that the Wnt targets are highly overrepresented among the genes that contain the predicted conserved LEF1/TCF binding site. This positively validated our bioinformatics approach to identifying LEF1/TCF targets. For further analysis, we selected only the genes with at least two LEF1/TCF motifs in one regulatory element, expecting fewer false-positive hits in this group. This presumption was based on the observation that more than 30% of the genes that contain pairs of TCF7L2 motifs in the same CNS display Wnt target-like patterns during mouse development [41]. We excluded genes with a "NULL" or "PREDICTED" description in the Ensembl database, creating a final list of 428 rat genes (Additional file 1), which we considered putative LEF1/TCF targets.
Putative LEF1/TCF targets are enriched in neuron-specific genes (GO analysis) We performed a gene ontology (GO) analysis of the 428 putative targets using DAVID Functional Annotation tools (http://david.abcc.ncifcrf.gov; accessed January 11, 2012) to identify functionally related groups of genes. As expected, our list was highly enriched in genes involved in development, patterning, and cell proliferation (Table 1). Other predictable categories of molecular function and biological processes were overrepresented: specific DNA binding (transcription factors), Wnt signaling, and cell migration, adhesion, and death. Unexpectedly, we found statistically significant enrichment in genes implicated in neuronal function, including genes for synaptic proteins and voltage-gated cation channels (VGCCs). Based on this analysis, we hypothesized that β-catenin-LEF1/TCF can regulate specific neuronal genetic programs. We collected the genes that were annotated with "synapse" and "gated-channel activity" terms and included a few other genes that failed to be annotated with the above GO terms, although they obviously belong to these categories (e.g., the γ-aminobutyric acid [GABA] receptor-associated protein gene Gabarap). Several genes were also excluded because they were not brain-specific (e.g., Chrnd, which encodes a neuromuscular cholinergic receptor). The final collection contained 41 putative neuronal LEF1/TCF targets grouped into five categories: nine VGCCs, 11 neurotransmitter receptors, eight proteins involved in synaptic vesicle organization, 10 synaptic structural proteins, and three other proteins implicated in synaptic transmission ( Table 2). All of these gene targets were chosen for experimental validation.
VGCC genes with conserved LEF1/TCF motifs compared with all other VGCC genes display preferential expression in the thalamus (RT-qPCR arrays) To initially validate our bioinformatics predictions, we concentrated on the VGCC group. The relative mRNA levels of the ensemble of VGCC genes were measured in three different parts of the forebrain: thalamus, cortex, and hippocampus. We expected that the VGCC genes that are postulated LEF1/TCF targets would be preferentially expressed in the thalamus because nuclear βcatenin and LEF1 and TCF7L2 transcription factors are present in thalamic neurons but not cortical or hippocampal neurons [13][14][15][16]. The comparative expression analysis was performed using custom-designed quantitative real-time polymerase chain reaction (RT-qPCR) arrays with six independent preparations from each brain region. The arrays contained all VGCC genes that have CNSs (Table 3), excluding a few genes that are known to be expressed only in non-neuronal tissues.
The expression of 53 of the 55 VGCC genes was detected in the forebrain; no signal was obtained for Kcnip1 and Kcnj1. Four of nine putative LEF1/TCF target VGCC genes were expressed at least two-fold higher in the thalamus than in the cortex and hippocampus, and none were expressed at a lower level. In the remaining group of 44 genes, four were expressed at a  Figure 1 Bioinformatics identification of putative LEF1/TCF target genes. The diagram presents the successive steps to select putative LEF1/TCF targets, beginning from the group of human-rat orthologs in the Ensembl database. Groups of genes are in rectangles, and computational procedures are in diamonds. In the blue frame, the crossing of our in silico-selected genes with known Wnt/β-catenin targets is shown. The p values (Fisher's Exact test) confirmed the enrichment of genes with at least one and at least two LEF1/TCF binding sites (predicted with Matinspector using Genomatix family V$LEFF) with the known targets.  higher level in the thalamus, and six were expressed at a lower level (Figure 2A). Fisher's Exact test, yielding a p value of 0.021, confirmed a correlation between thalamic expression and being a putative LEF1/TCF target. This validated the approach to identifying new LEF1/TCF targets based on combined in silico binding site prediction and high expression in the thalamus (i.e., the region with high levels of nuclear β-catenin and LEF1/TCF factors in neurons).
More than 20% of neuron-specific genes with conserved LEF1/TCF motifs are highly expressed in the thalamus (RT-qPCR arrays) We next expanded the analysis and profiled the expression of all 41 genes listed as putative neuronal targets of LEF1/TCF (Table 2) in the cortex, hippocampus, and thalamus. Another set of custom-designed RT-qPCR arrays and the previous six independent preparations from each brain region were used. Our aim was to determine which of the predicted targets are highly expressed in the thalamus, suggesting their actual regulation by β-catenin-LEF1/TCF. The expression of 40 of the 41 genes was detected in the forebrain; no signal was obtained for Htr3b. In the thalamus, nine genes (22.5%) were expressed at least two-fold higher than in the cortex and hippocampus ( Figure 2B), and the differences were statistically significant ( Figure 3). In the cortex vs. the other two regions, two genes (5%) were expressed at a higher level, and no single gene was highly expressed in the hippocampus ( Figure 2C).
The following genes had higher expression levels in the thalamus: four VGCC genes (Cacna1g, Cacna2d2, Kcna6, and Kcnh8, which were also observed in the previous experiment), four genes that represent neurotransmitter receptors (Drd3, Gabra3, Glra1, and Grid2), and the Calb2 gene that encodes the Ca 2+ buffer calretinin ( Figure 3). The genes from two other categories (i.e., the genes that encode synaptic vesicle proteins and structural synaptic proteins; Table 2) did not show specific thalamic expression.
β-catenin is associated with chromatin at the Cacna1g, Kcna6, Gabra3, Grid2, and Calb2 loci in the thalamus Gene profiling in the rat brain allowed us to observe a positive association between the relative expression of the neuronal genes with at least two conserved LEF1/ TCF motifs and the presence of β-catenin and LEF1/ TCF factors in the brain. To determine whether the β-catenin-LEF1/TCF complex might directly regulate Cacna1g, Cacna2d2, Kcna6, Kcnh8, Drd3, Gabra3, Glra1, Grid2, and Calb2, we analyzed the in vivo binding of β-catenin to LEF1/TCF motifs within the CNSs using a chromatin immunoprecipitation (ChIP) assay with designed primers ( Figure 4 and Table 4). We also examined the chromatin conformational state of the fragments that contained conserved LEF1/TCF motifs by precipitating them with an antibody specific for acetylhistone H3 (H3Ac; a hallmark of open chromatin [42,43]). In each ChIP assay, we compared four independent samples of chromatin isolated from the cortex, hippocampus, and thalamus.
We first assessed the acetylation status of histone H3 at the Gapdh promoter and Gapdh exon, an open and close chromatin region, respectively. We found high levels of H3Ac in the promoter while much lower levels in the first exon. This showed that our ChIP assays to monitor H3Ac was specific ( Figure 5A). We then analyzed the chromatin conformation of our genes of interest in fragments with conserved LEF1/TCF motifs. The chromatin fragments that were in close proximity to TSSs, Cacna1g-3 and Kcnh8-1, appeared to be in an open state. The same was observed for some fragments located distally from the TSSs (e.g., Grid2-1, -2, and −3), whereas other fragments precipitated with low efficiency (< 1%; e.g., Drd-1 and −2), indicating the closed conformation of chromatin ( Figure 5A). In most cases, no differences were found between the analyzed brain structures. However, some fragments (e.g., Gabra3-1, Cacna2d2-1, Cacna2d2-2, Cacna2d2-3, and Calb2-3) precipitated significantly more efficiently from the thalamic samples than from the cortex and hippocampus.
We then performed a ChIP assay with an anti-β-catenin antibody. To determine the background, normal immunoglobulin G (IgG) was used, which precipitated < 0.02% of the input. The signals for all of the examined fragments were at background levels in the case of the cortex and hippocampus ( Figure 5B). However, when the thalamic samples were used, fragments of Gabra3, Grid2, Cacna1g, Kcna6, and Calb2 precipitated with anti-β-catenin at levels of 0.04-0.1% ( Figure 5B), indicating the binding of β-catenin to these fragments. Moreover, for all of these The gene symbols are based on the Rat Genome Database (RGD). Table 3 List of genes that encode voltage-gated cation channels included in the custom "VGCC" arrays fragments, the differences in chromatin precipitation levels between the thalamic samples and other samples were statistically significant. This indicates that the aforementioned genes can be directly regulated by β-catenin and LEF1/TCF factors. Interestingly, no correlation was observed between the β-catenin-chromatin interaction ( Figure 5B) and acetylation status of histone H3 in the chromatin fragments ( Figure 5A).
Attenuation of β-catenin signaling leads to decreases in Cacna1g, Cacna2d2, Kcna6, Kcnh8, Gabra3, and Calb2 expression in cultured thalamic neurons Finally, we examined the effect of nuclear β-catenin removal in primary thalamic cultures on the expression of the nine genes identified by gene profiling. The cultures contained both neurons and glia (approximately 1:1), which is vital for the survival of thalamic neurons [38]. Thalamic neurons cultured in vitro maintain the nuclear localization of β-catenin [38]. To decrease its level, the cultures were treated with an adenovirus that carried Axin2, the product of which is a component of the βcatenin destruction complex and as such should reduce its cytoplasmic and nuclear pool. The control cultures were transduced with Gfp-expressing adenovirus. The percentage of β-catenin-positive neurons decreased from 40% to 10% in cultures with ectopic Axin2 expression ( Figure 6A). No nuclear β-catenin was detected in glia.
Afterward, the expression levels of the nine genes and control genes, Map2, Cacna1h, and Lef1, were measured by RT-qPCR. Similar expression levels of the neuronal marker Map2 in Gfp-and Axin2-expressing cultures verified the equal survival and condition of neurons ( Figure 6B). The level of Lef1, serving as a positive control as a classical target of β-catenin, decreased by~20%, confirming the impairment of β-catenin-dependent transcription in the treated cultures. The expression level of Cacna1h, a paralog of Cacna1g that lacks conserved LEF1/TCF motifs, was not modified. Of the nine examined genes, we observed a statistically significant decrease in the expression of six genes, Cacna1g, Cac-na2d2, Kcna6, Kcnh8, Gabra3, and Calb2, suggesting that these genes are authentic targets of β-catenin. The level of Glra1 expression increased, and Grid2 expression did not change. Drd3 transcription was not detected ( Figure 6B).
The actual decreases of the expression levels in given neurons were supposedly several times larger, considering that the changes in the level of nuclear β-catenin occurred in~30% of the cultured neurons. This percentage was estimated based on the following observations: (i) 40% of the neurons were β-catenin-positive in the control, and (ii) 12% of the neurons still maintained the nuclear localization of β-catenin after Axin2 adenovirus treatment ( Figure 6A). The most pronounced effects were observed for Cacna1g, Kcna6, and Calb2; together with the ChIP results ( Figure 5B), this indicates the high importance of β-catenin in the regulation of the expression of these genes in thalamic neurons. The ChIP and expression results were also consistent for Gabra3 and corroborated the regulation of this gene by β-catenin.

Discussion
The present study identified five new β-catenin target genes in thalamic neurons, in addition to previously described Cacna1g [16]. Three of them, Kcna6, Calb2, and Gabra3, were validated by ChIP in vivo and a lossof-function experiment in cultured neurons, confirming that they might be directly regulated by β-catenin. Two other genes, Cacna2d2 and Kcnh8, also displayed βcatenin-dependence in the latter experiment, although the binding of β-catenin to their regulatory elements was not found. Based on these data, we propose that βcatenin is a regulator of the electrophysiological properties of thalamic neurons in the adult brain.
Numerous genes that we selected in silico as potentially regulated by β-catenin belong to expected functional Table 3 List of genes that encode voltage-gated cation channels included in the custom "VGCC" arrays (Continued) 45 Kcnj5 categories: transcription regulation, cell proliferation, morphogenesis, motility, adhesion, differentiation, and programmed cell death. Similar clusters were observed by others in the genes bound by TCF7L2 in a human colorectal cancer cell line [40]. These results support the well-established role of Wnt/β-catenin in development.
Interestingly, the genes involved in neuronal differentiation exhibited the highest enrichment scores in our list. This was consistent with a study that identified β-catenin-LEF1/TCF targets based on a ChIP assay in NIH3T3 cells. Thirty percent of the target genes were implicated in developmental processes, and more than half of the targets from this group were involved in neuronal development [44]. Indeed, Wnt signaling has been particularly implicated in central nervous system development, from early brain patterning to embryonic and adult neurogenesis [30][31][32][33][45][46][47][48][49][50]. Although our in silico analysis corroborated the involvement of LEF1/TCFs in the regulation of well-known thalamus (x-axis) and hippocampus (y-axis). (Right plot) Hippocampus vs. thalamus (x-axis) and cortex (y-axis). Notice that this group of genes is not preferentially expressed in the cortex or hippocampus. n = 6 independent biological samples.
groups of genes, it also identified a group that has not been previously proposed to be a β-catenin-LEF1/TCF target. These were the genes of proteins involved in signal transmission in neurons, including voltage-gated ion channels, neurotransmitter receptors, synaptic vesicle proteins, and synaptic structural proteins. Moreover, we provided experimental evidence of the authentic regulation of some of these genes by β-catenin. The above gene clusters have not yet been identified, probably because the screenings for β-catenin target genes were performed on established cell lines or cancer cells [40,44,[51][52][53][54]. Additionally, studies of hippocampal neurons, in which β-catenin nuclear translocation was observed after NMDA stimulation, did not attempt to identify specific neuronal targets [35][36][37]. While examining the β-catenin-chromatin association and acetylation of histone H3, we did not observe any relationship between these two phenomena in the analyzed regions with the LEF1/TCF motif. This suggests that the interaction between the β-catenin-LEF1/TCF complex and DNA might not require the open conformation of chromatin. These results may also suggest that the β-catenin complex does not always increase histone acetylation, although it potentially has such an ability [55,56]. This is consistent with a recent study performed on embryonic stem cells, in which knockdown of Tcf 7 and Tcf 7l2 did not affect the active chromatin conformation of their targets [57]. We also noticed an interesting pattern of LEF1/TCF motif occurrence in the examined genes. The motifs were usually clustered downstream of the first exon and not in the promoter regions. This suggests that the transcription of these genes may be regulated by LEF1/TCF factors by gene looping, which has been demonstrated for the COX2 and MMP13 genes, in which LEF1/TCF binding sites were located in the 3 0 untranslated region [58].
The neuronal genes with conserved LEF1/TCF motifs that were highly expressed in the thalamus, the regulation of which by β-catenin was confirmed experimentally, encode proteins involved in neuronal excitability. Cav3.1 (encoded by Cacna1g), Cavα 2 δ2 (Cacna2d2), Kv1.6 (Kcna6), Kv12.1/ELK1 (Kcnh8), and GABA A receptor 3 (Gabra3) are all voltage-or ligand-gated ion channels [59][60][61]. As such, they underlie the cell membrane conductance of Ca 2+ , K + , and Cl -(in the case of the GABA receptor) ions and directly propagate, inhibit or modify electric signals [62][63][64]. Calretinin, in turn, is an intracellular Ca 2+ -binding protein [65,66] with diverse functions, including the modulation of intrinsic neuronal excitability [67]. We propose that β-catenin contributes to the proper excitability of thalamic neurons by regulating the expression of the above genes. However, more research is required to determine the real impact of β-catenin and LEF1/TCF factors on the expression of the identified genes and electrophysiology of the thalamus.
The other classes of putative neuronal targets of LEF1/ TCF (i.e., the genes that encode structural synaptic proteins, mainly with the PDZ domain, and synaptic vesicle proteins) did not show high expression in the thalamus. However, they still might be regulated by β-catenin and LEF1/TCF factors in some subtypes of neurons or under specific physiological conditions because the regulation of gene expression by β-catenin is very much contextdependent [8,44,68,69]. Particularly interesting would be the exploration of this possibility in future research because a membranous fraction of β-catenin interacts with  Figure 3 Expression analysis of putative LEF1/TCF targets in the forebrain. Volcano plots arrange genes along the dimensions of (x) mean expression fold difference between two brain structures and (y) p value (Student's t-test). A logarithmic scale is used. Red frames surround the plot area, in which the expression in the thalamus is at least two-fold higher than in the other structures, and the difference is statistically significant (p < 0.05). The genes inside the frames are considered likely LEF1/TCF targets in the thalamus. On every plot, the genes that met the criterion of a statistically significant two-fold expression difference between the two structures are labeled. Those that are higher in the thalamus vs. cortex and hippocampus are in green. n = 6 independent biological samples.
PDZ proteins in synapses and is implicated in synaptic vesicle localization [70][71][72][73]. The role of nuclear β-catenin in the regulation of PDZ and synaptic vesicle protein expression might complement the function of membranous β-catenin in neurons. We do not yet know whether variations in the nuclear level of β-catenin affect the expression of genes that encode VGCCs and neurotransmitter receptors and shape neuronal excitability in vivo. If so, then we could speculate that the inappropriate activity of β-catenin might affect proper signal transmission in thalamocortical circuits. Thalamocortical desynchronization underlies absence epilepsy [74], and many anticonvulsant drugs target voltage-gated channels (e.g., T-type Ca 2+ channels [63]). Specifically, the T-type voltage-gated channel Cav3.1 has been proposed to be implicated in absence seizures [75,76], in addition to the Cavα 2 δ2 regulatory subunit of voltage-gated channels [77] and GABA A receptor 3 [78]. Schizophrenia has also been associated with thalamic dysfunction [79][80][81][82][83][84]. Moreover, some variants of Tcf7l2 have been recently shown to be a risk factor in schizophrenia [23,85], and a group of synaptic genes involved in excitability has been found to be associated with the risk of schizophrenia [86]. Interestingly, Gabra3-deficient mice display impairments in sensorimotor gating, which is a feature of this disorder [87]. These results suggest a possible role for β-catenin-dependent gene expression in thalamic pathologies, but further in vivo studies are required to elucidate this issue.

Conclusions
We identified a novel group of genes regulated by β-catenin-LEF1/TCF that encode proteins that underlie the transmission of nerve impulses. These results point to a novel role for β-catenin in the thalamus, in which β-catenin is constantly present in the cell nuclei. The identified and confirmed genes are directly involved in neuronal excitability, suggesting that β-catenin and LEF1/TCF transcription factors maintain the proper activity of thalamocortical circuits. The presented results may implicate the involvement of β-catenin transcriptional activity in thalamic pathologies, such as absence epilepsy and schizophrenia.

In silico screening and validation of the method
For each human-rat orthologous gene pair in Ensembl (version 39), 20 kb flanks of the genomic sequence from −10 kb to +10 kb from TSSs in either species were aligned using the AVID global alignment algorithm. Conserved noncoding sequences, defined as sequence windows at least 100 bp long with at least 75% identity between human and rat, were selected as putative regulatory regions. The binding sites for LEF1/TCF were separately predicted using the position weight matrices in the V$LEFF family (Genomatix Matrix Family Library, version 6.2) with the MatInspector program for the human and rat sequences of each CNS pair. V$LEFF motifs with a nonzero number of instances in both the human and rat sequences of the same CNS pair (not necessarily the same number or in the same AVID-aligned position) were considered conserved. To validate the method, the group of genes with at least one or at least two conserved V$LEFF motifs were crossed with known Wnt target genes listed at http://www.stanford.edu/ group/nusselab/cgi-bin/wnt/target_genes (accessed January 11, 2011). The Wnt targets that had no CNSs, were indirect targets, or had confusing names were excluded from this analysis.

Gapdh -E TTGTGACAAAGTGGACATTG AACTTGCCGTGGGTAGAGTC
The positions of the primers are shown in Figure 4.
Gene enrichment in annotation terms was measured by determining the EASE score [91], which is a conservative correction of Fisher's Exact p value, for the proportions of genes that fell into each GO term.

Animal care
To perform this study, we used

Chromatin isolation, shearing, and immunoprecipitation
Three-month-old male Wistar rats were sacrificed by cervical dislocation, and the brains were removed and sectioned. Approximately 200 mg of cortical, hippocampal, and thalamic tissues obtained from two rats were chopped with a razor blade and fixed at room temperature for 30 min in 1 ml of fresh 3% formaldehyde in phosphate-buffered saline (PBS) supplemented with protease inhibitor cocktail (Roche) in an Eppendorf tube. Fixation was stopped by adding glycerol to a final concentration of 0.125 mM and incubation for a further 10 min to quench formaldehyde. The tissue was then spun down for 5 min at 1,000 × g and briefly homogenized in ice-cold PBS supplemented with protease inhibitor cocktail using an Eppendorf-fitting pestle homogenizer. After two cycles of washing in PBS (i.e., homogenization and centrifugation), the samples were lysed in 2.5 ml of lysis buffer (1% sodium dodecyl sulfate, 10 mM ethylenediamine tetraacetic acid [EDTA], 50 mM TrisHCl, pH 8.1, and protease inhibitor cocktail) for 30 min at 4°C and then disrupted in a Potter-Elvehjem homogenizer (~100 strokes). The homogenates were centrifuged for 10 min at 2,400 × g to pellet the nuclei, which were then frozen at −80°. After thawing, the samples were diluted in digestion buffer from the ChIP-IT Express Enzymatic kit (Active Motif ) and sonicated on ice for 10 × 20 s with 40 s breaks at a 50% duty cycle and 50% power using a Sonopuls Bandeline sonicator to disrupt the nuclei and preshear the chromatin. An enzymatic shearing protocol using the ChIP-IT Express Enzymatic kit was then followed. Shearing was conducted with 40   separated from debris by centrifugation at 18,000 × g. This treatment yielded 125-500 bp DNA fragments. A subsequent immunoprecipitation procedure was performed on an equivalent 30 mg of tissue with 10 μg of anti-rabbit IgG (Sigma-Aldrich), 5 μg of anti-acetylhistone H3 (rabbit polyclonal; Millipore), or 40 μg of anti-β-catenin (rabbit polyclonal; Santa Cruz Biotechnology) antibodies, according to the Millipore protocol using salmon sperm DNA protein A-agarose (Millipore). Immunoprecipitated DNA was used as a template for RT-qPCR with SYBR Green chemistry. The obtained data were analyzed using the relative quantification method with the 2 -ΔCT formula (DCT = CT input DNA -CT immunoprecipitated DNA ; CT is the cycle threshold) and are presented as a percentage of the input. The primers used in the ChIP assay are listed in Table 4.

Primary neuronal cultures
Dissociated primary thalamic cells were obtained from the brains of embryonic day 19 Wistar rat embryos as described previously [38]. The cells were plated on coverslips coated with poly-D-lysine (30 μg/ml; Sigma) at a density of 2.5 × 10 5 cells per well of a 24-well plate in Minimal Essential Medium (MEM; Gibco) supplemented with 10% fetal bovine serum and 0.2 mM glutamine (Sigma). The next day, the medium was replaced with Neurobasal (Gibco) supplemented with B27 (Gibco), 0.5 mM glutamine, 12.5 mM glutamate, and penicillin/ streptomycin (Sigma) mixed in a 1:1 ratio with the cortical neuron conditioned medium.

Adenoviral constructs and neuron transduction
The construction of recombinant adenoviral plasmids that expressed Gfp and Axin2 and the preparation and purification of the adenoviruses were described previously [38]. On day 4 in vitro, purified virus suspensions at a multiplicity of infection of 100 were added to each well of a 24-well plate for 6 h at 37°C. The virus solution was then removed and replaced with fresh medium. Seventy-two hours later, the transduced cultures were collected for RNA isolation or fixation. Drd3 mRNA was not detected (ND). n = 9 independent biological samples. Error bars indicate SD. ***p < 0.001, **p < 0.01, *p < 0.05 (Student's t-test).

Immunofluorescence and microscopic analysis
The immunofluorescence analysis of β-catenin localization was performed as described previously [38]. Briefly, the cells were incubated overnight at 4°C with anti-β-catenin rabbit antibody (1:250; Santa Cruz Biotechnology) in a humid chamber followed by incubation with anti-NeuN mouse antibody (1:150; Chemicon) for 1.5 h at room temperature. All of the antibodies were diluted in 2% bovine serum albumin. The slides were acquired under a confocal microscope (Zeiss LSM5 Exciter). The fluorescence intensity of intracellular βcatenin was analyzed using Cell Profiler Software. Based on a fluorescent histogram under the control conditions, two separate cell populations with low and high intensity were identified, and the intensity threshold between them was established. In an experimental variant, the cells with fluorescence greater than the threshold were counted as β-catenin-positive cells.

RNA isolation from brain tissue and neuronal cultures
To isolate RNA from fresh brains, the brains from male, 2-month-old Wistar rats were sectioned, and different structures were homogenized separately using a Potter-Elvehjem homogenizer. RNA was extracted with the RNeasy Lipid Tissue Mini kit with additional DNase treatment (Qiagen). For neuronal cultures, the cells from one to three wells of a 24-well plate were pooled together and homogenized with Qiashredders (Qiagen), and RNA was isolated with the RNeasy Plus Mini kit (Qiagen). cDNA was synthesized by reverse transcription (SuperScript III RnaseH-, Invitrogen).

Gene expression analysis by RT-qPCR
All RT-qPCR reactions were performed using the 7900HT system (Applied Biosystems). Gene profiling in the brain was performed with custom-designed TaqMan Low Density Arrays (Applied Biosystems), hereinafter referred to as RT-qPCR arrays. Approximately 1.5 μg of cDNA was loaded per array (384 assays). The obtained data were analyzed using the relative quantification method with the 2 -ΔCT formula (ΔCT = CT target -CT Gapdh ; CT is the cycle threshold). In thalamic cultures, gene expression levels were examined in individual qRT-PCR reactions. For Cac-na1g, Cacna2d2, Kcna6, Knch8, Drd3, Glra1, Gabra3, Grid2, and Calb1, commercial primers were used (Qiagen). For Map2, commercial primers and FAM dyelabeled TAqMan probes (Applied Biosystems) were used. Lef1 and Gapdh expression was measured with the following primers: Lef1 (forward, CCCACACGGACAGC-GACCTA; reverse, TAGGCTCCTFTTCCTTTCTCT), Gapdh (forward, TGACTCTACCCACGGCAAGTTCAA; reverse, ACGACATACTCAGCACCAGCATCA). SYBR Green chemistry was used, with the exception of Map2, for which TaqMan was used (Applied Biosystems). The obtained data were quantified using the relative standard method.

Statistical analysis
Three types of statistical tests were used. Two-tailed Fisher's Exact test was used to verify the enrichment of genes with LEF1/TCF binding sites with the known Wnt targets (Figure 1). This test was also performed to calculate p values for the relationships between high expression in the thalamus and being a putative LEF1/ TCF target (Figure 2; n = 6), and compare the percentage of β-catenin-positive neurons in two types of thalamic cultures ( Figure 6A). To calculate p values for differences in gene expression levels between two brain regions, we ran two-tailed paired Student's t-test, using ΔCT values for the results obtained with the RT-qPCR arrays (Figure 3; n = 6) or using relative values for the results obtained with individual RT-qPCRs ( Figure 6; n = 9).
To calculate p values in the ChIP experiment, in which we compared three brain regions at the same time, we used repeated measures ANOVA followed by Tukey's post hoc test ( Figure 5; n = 4).