Dynamic FoxP2 levels in male zebra finches are linked to morphology of adult-born Area X medium spiny neurons

Page created by Guy Weber
 
CONTINUE READING
Dynamic FoxP2 levels in male zebra finches are linked to morphology of adult-born Area X medium spiny neurons
www.nature.com/scientificreports

            OPEN           Dynamic FoxP2 levels in male zebra
                           finches are linked to morphology
                           of adult-born Area X medium spiny
                           neurons
                           Jennifer Kosubek-Langer* & Constance Scharff

                           The transcription factor FOXP2 is crucial for the formation and function of cortico-striatal circuits.
                           FOXP2 mutations are associated with specific speech and language impairments. In songbirds,
                           experimentally altered FoxP2 expression levels in the striatal song nucleus Area X impair vocal learning
                           and song production. Overall FoxP2 protein levels in Area X are low in adult zebra finches and decrease
                           further with singing. However, some Area X medium spiny neurons (MSNs) express FoxP2 at high
                           levels (FoxP2high MSNs) and singing does not change this. Because Area X receives many new neurons
                           throughout adulthood, we hypothesized that the FoxP2high MSNs are newly recruited neurons, not yet
                           integrated into the local Area X circuitry and thus not active during singing. Contrary to our expectation,
                           FoxP2 protein levels did not predict whether new MSNs were active during singing, assayed via
                           immediate early gene expression. However, new FoxP2high MSNs had more complex dendrites, higher
                           spine density and more mushroom spines than new FoxP2low MSNs. In addition, FoxP2 expression levels
                           correlated positively with nucleus size of new MSNs. Together, our data suggest that dynamic FoxP2
                           levels in new MSNs shape their morphology during maturation and their incorporation into a neural
                           circuit that enables the maintenance and social modulation of adult birdsong.

                           The forkhead box P2 transcription factor (FOXP2) is linked to speech and language disorders. Heterozygous
                           FOXP2 mutations in humans affect both the coordination of fine orofacial movements and language percep-
                           tion1–3. Because songbirds – like humans – need to learn most of their communicative vocalizations, they offer
                           a unique model to study the role of FoxP2 (for nomenclature FOXP2/FoxP2 see Methods) for vocal learning
                           and for the maintenance of learned vocalizations as adults4. Studying the relationship between FoxP2 and vocal
                           learning in songbirds may inform the neurogenetic mechanism underlying the speech deficits in patients carrying
                           FOXP2 mutations for the following reasons. The FoxP2 protein coding sequence is highly conserved between
                           humans and songbirds as are the brain expression patterns, notably in the cerebellum and striatum5–7. Moreover,
                           genetic manipulations of FoxP2 expression levels in the striatal song nucleus Area X during the critical phase of
                           song learning lead to inaccurate and incomplete imitation of the tutor’s song and more variable vocal produc-
                           tion3,8–10. This phenotype bears similarities to the specific speech deficits called developmental verbal dyspraxia,
                           DVD (or childhood apraxia of speech), that patients carrying FOXP2 mutations suffer from. The core-phenotype
                           of DVD consists of altered precision, consistency and sequencing of movements underlying speech in the absence
                           of neuromuscular deficits11. In addition, altered FoxP2 levels in adult Area X affect the dopaminergic modulation
                           of corticostriatal signaling important to song variability and affect song maintenance12,13, stressing the fact that
                           tight regulation of FoxP2 expression is a prerequisite for correct neural transmission in differentiated neural cir-
                           cuits. Additional effects of Foxp2 and its disruption on the embryonic development and the function of neural cir-
                           cuits have been described in mice14–23. Further evidence for the biological relevance of tight regulation of FoxP2
                           expression levels comes from the following studies. FoxP2 expression levels in Area X transiently increase during
                           song learning but are lower in adults7,24. Singing decreases overall FoxP2 levels in Area X but not in the surround-
                           ing striatum and the degree of FoxP2 down regulation correlates with the amount of produced song25–27. This
                           relationship is missing in deafened birds, pointing to an important role of auditory feedback for singing-driven
                           FoxP2 down regulation27. How does singing affect FoxP2 expression at the cellular level? Medium spiny neurons

                           Department of Animal Behavior, Institute of Biology, Freie Universität Berlin, Berlin, Germany. *email: kosubekla@
                           zedat.fu-berlin.de

Scientific Reports |   (2020) 10:4787 | https://doi.org/10.1038/s41598-020-61740-6                                                           1
Dynamic FoxP2 levels in male zebra finches are linked to morphology of adult-born Area X medium spiny neurons
www.nature.com/scientificreports/                                                               www.nature.com/scientificreports

                           (MSNs), the most abundant cell type in the avian striatum, predominantly express FoxP2 at low levels (FoxP2low)
                           while a subset expresses FoxP2 at very high levels (FoxP2high). Both subtypes are not equally affected by singing;
                           the density of FoxP2high MSNs is not measurably different after singing, contrary to the decreasing density of
                           FoxP2low MSNs24. The authors hypothesized that the difference might be due to the neuronal age. Adult Area X
                           constantly receives new MSNs that originate at the ventricular zone28–33. FoxP2high MSNs colocalize more fre-
                           quently with a marker for new neurons than FoxP2low MSNs24. Recently we showed that new MSNs mature and
                           participate in singing activity – as measured by immediate early gene activation – within a timeframe of six
                           weeks34. Whether FoxP2 influences not only the function but also the integration of new neurons into existing
                           circuits is still an open question. Based on the results of Thompson et al. (ref. 21) we hypothesized that FoxP2high
                           MSNs are newly recruited into Area X and need to become FoxP2low MSNs before they can participate in singing.
                           To test this, we labelled neuronal progenitors in adult zebra finches. At different time points after these cells had
                           migrated into Area X, we quantified their expression levels of FoxP2 and whether they also expressed the imme-
                           diate early gene expression EGR-1 after singing.
                               In rodents, Foxp2 expression is associated with neurite outgrowth, neuronal morphology and synapse forma-
                           tion in cortico-striatal circuits18,19,35–37. Foxp2 expression levels vary in striatal MNSs and these differences may
                           be relevant for the morphology of striatal MSNs. Dopamine receptor 1 (D1) expressing MSNs38 express Foxp2
                           at higher levels than dopamine receptor 2 (D2) expressing MSNs23,35. These differences in Foxp2 levels may be
                           linked to anatomical differences between D1 and D2 MSNs35,39. Furthermore, in mice carrying humanized Foxp2
                           alleles (Foxp2hum/hum mice), Foxp2high MSNs are more numerous in the dorsal striatum and their MSNs have
                           longer dendrites than wildtype mice36,37. Based on the latter results we hypothesized that FoxP2 levels of new
                           MSNs in adult songbirds correlate with their neural morphology. To test this, we virally labelled neural progen-
                           itors and analyzed their FoxP2 expression, dendrite complexity and spine density after migration into Area X.

                           Results
                           Dynamic FoxP2 levels in new MSNs in Area X. To assess FoxP2 protein levels in individual newborn
                           neurons we labelled progenitor cells with Bromodeoxyuridine (BrdU) and detected BrdU+/FoxP2+ cells after
                           21, 31 and 42 days post BrdU injection (dpi) in Area X of adult male zebra finches (Fig. 1a,b). We found that
                           FoxP2 expression in Area X was very variable, with some neurons expressing FoxP2 at particularly high levels
                           and some at low levels (Fig. 1c). At 21 dpi and at 31 dpi, the mean pixel intensities of all BrdU+/FoxP2+ cells
                           formed a bimodal distribution, whereas at 42 dpi the distribution was unimodal and shifted to low FoxP2 expres-
                           sion (Fig. 1d). We classified all neurons that had expression intensities within the top 30% of the measured pixel
                           intensity distribution as FoxP2high neurons. Neurons within the bottom 30% of the measured pixel intensity dis-
                           tribution were considered as FoxP2low. Because we were interested in the two extremes of the expression levels in
                           this study, we did not analyze the new neurons with intermediate FoxP2 expression levels further (29.3% ± 9.4,
                           SD, see Methods). At 21 dpi and at 31 dpi 36.17% ± 6.02 (SEM) and 34.91% ± 2.53 (SEM) percent of all BrdU+
                           cells were FoxP2high neurons, respectively. At 42 dpi the percentage of FoxP2high cells had significantly decreased to
                           12.95% ± 2.87 (SEM) on average (p = 0.0077, Kruskal-Wallis test, Fig. 1f,g). The percentage of FoxP2low neurons
                           increased significantly from 21 dpi (34.59 ± 5.86, SEM) and 31 dpi (34.06 ± 3.35, SEM) to 42 dpi (59.19 ± 4.02,
                           SEM, p = 0.013, ANOVA, Fig. 1e,g). We also noticed that BrdU+ cells varied in their nucleus size. Quantification
                           revealed that the distribution of the nucleus size was shifted towards bigger nuclei at 31 dpi (data not shown).
                           BrdU+/FoxP2+ cells at 31 dpi had significantly bigger nuclei (7.37 µm ± 0.94 (SD) than BrdU+/FoxP2+ cell at
                           21 dpi (6.88 µm ± 0.9 (SD), p = 0.00025, chi-square = 75.358, df = 2) or 42 dpi (6.74 µm ± 0.59 (SD), p = 2 × 10−6,
                           chi-square = 75.358, df = 2, data not shown). Interestingly there was a significant positive relationship with a low
                           effect size between nucleus size and FoxP2 expression levels in all three experimental groups (21 dpi: r2 = 0.086,
                           p = 0.017, 31 dpi: r2 = 0.083, p = 5.2 × 10-7, 42 dpi = r2 = 0.039, p = 0.0012, Fig. 1h).

                           Singing induced activity of new MSNs is independent of FoxP2 levels.                         In adult zebra finches
                           FoxP2 expression levels in Area X are behaviourally regulated. Undirected singing leads to downregulation of
                           FoxP2 mRNA and protein24–26. Undirected singing is also associated with expression of the immediate early gene
                           EGR-1 in Area X, which is therefore often used as a molecular readout of the neuronal activity associated with
                           undirected song40–44. We hypothesized that FoxP2 levels in new FoxP2high neurons needed to be downregulated
                           before activation by singing could occur, resulting in EGR-1 expression. Consequently, we did not expect to find
                           BrdU+/FoxP2high/EGR-1+ MSNs in Area X. To test this, we analyzed BrdU+/EGR-1 ± /FoxP2+ cells in birds
                           that had sung before sacrifice after 21 dpi, 31 dpi and 42 dpi.
                               Contrary to our hypothesis we found BrdU+/FoxP2high/EGR-1+ cells in Area X (Fig. 2a) in all groups, with
                           large differences between birds. At 31 dpi, on average 41.45% ± 15.43 (SEM) of new neurons expressed FoxP2high
                           and were also activated by singing (BrdU+/FoxP2high/EGR-1+) and 40.28% ± 13.34 (SEM) of the new neu-
                           rons that expressed FoxP2high were not activated by singing (BrdU+/FoxP2high/EGR-1-) (Fig. 2b). At 42 dpi, the
                           more birds had sung the fewer FoxP2high new neurons were found, resulting in a significant negative relation-
                           ship between the number of BrdU+/FoxP2+/EGR-1+ neurons and the number of motifs sung before sacrifice
                           (r2 = 0.753, p = 0.025, Fig. 2c) which was not the case at 21 dpi (r2 = 0.168, p = 0.493, data not shown) nor at 31
                           dpi (r2 = 0.164, p = 0.425, data not shown).

                           FoxP2 levels in new MSNs influence dendritic arborization and spine density.                             Because of the
                           relationship of FoxP2 levels and nucleus size we wanted to further characterize the morphology of new neurons
                           expressing different FoxP2 levels. We used a lentiviral approach to express Green Fluorescent Protein (GFP) in
                           progenitor cells at their place of birth in the lateral wall of the lateral ventricle (Fig. 3a). After 31 and 42 dpi we
                           found labeled neurons (GFP+) in Area X and the surrounding striatum (Fig. 3b). GFP+ neurons exhibited the
                           typical morphology of medium spiny neurons with small somata and spiny dendrites (Fig. 3c) and expressed

Scientific Reports |   (2020) 10:4787 | https://doi.org/10.1038/s41598-020-61740-6                                                              2
Dynamic FoxP2 levels in male zebra finches are linked to morphology of adult-born Area X medium spiny neurons
www.nature.com/scientificreports/                                                             www.nature.com/scientificreports

                           Figure 1. Dynamic FoxP2 expression levels in new MSNs. (a) The song motor pathway (main nuclei outlined
                           in black) controls the vocal organ (syrinx). The anterior forebrain pathway (shown in green) forms a cortico
                           (HVCm, lMAN) -basal ganglia (Area X) –thalamo (DLM) –cortical (RA) loop. (b) Experimental schedule: adult
                           male zebra finches received BrdU on five consecutive days and were sacrificed 21, 31 or 42 days after injections
                           (dpi). (c) Examples of FoxP2high (top row, immunoreactivity shown in green) and FoxP2low (bottom row) new
                           (BrdU+ immunoreactivity shown in white) MSNs in Area X at 31 dpi. Nuclear expression of FoxP2 and BrdU
                           coincides with DAPI label in blue. (d) Density plots of FoxP2 pixel intensities of individual new MSNs in Area
                           X at different time points after BrdU injections. The color scheme indicates increasing pixel intensity from low
                           (blue) to high (yellow) intensity. Ticks at the bottom of each plot represent individual MSNs. (e) Percentage
                           of new FoxP2low MSNs significantly increases from 21/31 dpi to 42 dpi. (f) Percentage of new FoxP2high MSNs
                           significantly decreases from 21/31 dpi to 42 dpi. (g) Empirical cumulative distribution function (ECDF) of
                           FoxP2 pixel intensities of individual new MSNs in Area X. FoxP2 pixel intensities are similar at 21 and 31 dpi
                           and are lower at 42 dpi. (h) FoxP2 pixel intensities of new MSNs positively correlate with nucleus diameter. Each
                           dot represents one new MSNs. Sample size (d-h): 733 MSNs of 17 zebra finches. *p ≤ 0.05, **p ≤ 0.01. Scale
                           bar: 20 µm (g). RA, Robust nucleus of the arcopallium; LMAN, Lateral magnocellular nucleus of the anterior
                           nidopallium; NXIIts, tracheosyringeal part of the hypoglossal nucleus; DLM, Dorsal lateral nucleus of the
                           medial thalamus.

                           FoxP2 (Fig. 3d,e). First, we traced FoxP2high and FoxP2low new neurons at 31 dpi (Fig. 4a) and found that FoxP2high
                           neurons had more primary dendrites (p = 0.021, Mann Whitney test, Fig. 4b), a higher total branch length
                           (p = 0.003, Mann Whitney test, Fig. 4e) and thicker dendrites than FoxP2low neurons (p = 0.003, t-test, Fig. 4f).
                           Second, we analyzed the extent of dendritic arborization of GFP+/FoxP2+ neurons using a Sholl analysis (see

Scientific Reports |   (2020) 10:4787 | https://doi.org/10.1038/s41598-020-61740-6                                                          3
Dynamic FoxP2 levels in male zebra finches are linked to morphology of adult-born Area X medium spiny neurons
www.nature.com/scientificreports/                                                             www.nature.com/scientificreports

                           Figure 2. Singing-induced EGR-1 activation of new MSNs is independent of FoxP2 levels. (a) The white
                           arrow in all 4 panels points to a new (BrdU+ immunoreactivity, white) MSN that expresses FoxP2high
                           (immunoreactivity, green) and also EGR-1 (immunoreactivity, purple) after undirected singing. The blue DAPI
                           staining shows other cells that are not new, but express FoxP2 at low levels, some of which also express EGR-
                           1. Scale bar: 10 µm. (b) At 31 dpi new neurons can either be activated by undirected singing (EGR-1+, right
                           column) or not (EGR-1-, left column). In both cases, the new MSN can either express FoxP2low or FoxP2high. (c)
                           At 42 dpi the number of FoxP2/BrdU+/EGR-1+ neurons negatively correlate with the number of motifs sung
                           during the 90 min before sacrifice. Sample size (b): 108 MSNs of 6 zebra finches. Sample size (c): 156 MSNs of 6
                           zebra finches.

                           Figure 3. New MSNs that were GFP-labeled as progenitors at the ventricular zone and migrated to Area X
                           express FoxP2. (a) Injections with a lentivirus into the ventricular zone resulted in GFP expression in the neural
                           progenitors of MSN. (b) 31 dpi after viral injections, many GFP-labelled neurons can be observed in Area X and
                           the surrounding striatum. (c) New GFP-expressing MSNs in Area X. Somata in the dashed boxes are magnified
                           in (d) and (e). (d,e) Virally labelled MSNs (GFP immunoreactivity, white) express FoxP2 (immunoreactivity,
                           green). Scale bars: 10 µm (d, e), 50 µm (c), 1 mm (b).

Scientific Reports |   (2020) 10:4787 | https://doi.org/10.1038/s41598-020-61740-6                                                         4
Dynamic FoxP2 levels in male zebra finches are linked to morphology of adult-born Area X medium spiny neurons
www.nature.com/scientificreports/                                                             www.nature.com/scientificreports

                           Figure 4. FoxP2 expression levels are linked to distinct morphologies of new MSNs at 31 dpi. (a) Examples of
                           dendrite tracings of new FoxP2low and FoxP2high MSNs. The black dot marks the center of the soma. (b) Number
                           of primary dendrites is significantly increased in FoxP2high compared to FoxP2low new MSN. (c) Maximal
                           number of intersections between dendrites and Sholl circles is significantly higher in FoxP2high than in FoxP2low
                           new MSN. (d) Sholl analysis revealed that dendrites of new FoxP2high MSNs had more complex arborizations
                           as indicated by more intersections at 20 µm from the soma than new FoxP2low MSNs. Shown are mean ± SEM.
                           Data points of FoxP2high neurons were slightly shifted to the right for better visibility. (e) FoxP2high new MSNs
                           have a significantly higher total branch length than FoxP2low new MSNs. (f) At 31 dpi new FoxP2high MSNs have
                           thicker dendrites than new FoxP2low MSNs. Sample size (a-f): 52 MSNs of 4 zebra finches. Scale bars: 25 µm.
                           *p ≤ 0.05; **p ≤ 0.01.

                           Methods). At 31 dpi FoxP2high neurons had more intersections at 20 µm distance from the soma (p = 0.024, paired
                           t-test, Fig. 4a,d) and a higher number of maximal intersections (p = 0.019, t-test, Fig. 4c) than FoxP2low neurons,
                           reflecting more extensive dendritic arborizations in FoxP2high than FoxP2low neurons. Second, Third, we used a
                           semi-automated quantification approach to assess the number of dendritic spines in GFP+/FoxP2+ neurons at
                           31 dpi. FoxP2high neurons in Area X had more dendritic spines than FoxP2low neurons (p = 0.034, paired t-test,
                           Fig. 5a–c). Overall, FoxP2high neurons had more mushroom spines than FoxP2low neurons (p = 0.0186, Mann
                           Whitney test, Fig. 5d). There was no difference in the number of stubby spines (p = 0.819, Mann Whitney test) or
                           thin spines (p = 0.409, Mann Whitney test) between FoxP2high and FoxP2low neurons (Fig. 5d). Because of the dif-
                           ference in mushroom spine number, we assumed that FoxP2 expression levels might influence mushroom spine
                           head size, too. However, quantification revealed no difference in mushroom spine head size between FoxP2high
                           and FoxP2low neurons (p = 0.317, Mann Whitney test, Fig. 5e). In a last step we compared spine densities between
                           new MSNs at 31 dpi and 42 dpi. Since at 42 dpi only few new neurons were FoxP2high we included only FoxP2low
                           new neurons in this analysis (Fig. 5f). At 42 dpi, the spine density of FoxP2low new neurons was higher than in
                           FoxP2low neurons at 31 dpi (p = 2.517 × 10−4, t-test, Fig. 5f). This elevated spine density was largely due to an
                           increase in thin spines (p = 8.422 × 10−4, Mann Whitney test) and not mushroom spines (p = 0.39, t-test) or
                           stubby spines (p = 0.119, Mann Whitney test).

                           Discussion
                           In the present study we investigated the dynamics of FoxP2 expression in adult-born MSNs in the striatal song
                           nucleus Area X of adult male zebra finches. We show that the new MSN strongly expressed FoxP2 at their arrival
                           in Area X from the ventricular zone (VZ) where they were born 21 days prior. During this stage and at interme-
                           diate maturation stages (31 days) one third of new MSNs expressed FoxP2 at high levels. At the late maturation
                           stage (42 days) most new MSNs expressed FoxP2 at low levels (Fig. 6). Together with our previous data we con-
                           clude that reaching low FoxP2 levels is a sign that adult-born born MSN in Area X have reached maturity by 6
                           weeks after their generation in the VZ34.

Scientific Reports |   (2020) 10:4787 | https://doi.org/10.1038/s41598-020-61740-6                                                          5
www.nature.com/scientificreports/                                                             www.nature.com/scientificreports

                           Figure 5. FoxP2 levels are associated with dendritic spine density of new MSNs. (a) New FoxP2high MSNs
                           had significantly more dendritic spines than new FoxP2low MSNs at 31 dpi (shown are mean ± SEM). Lines
                           connect data from the same animal. (b) Confocal 12 µm projection showing an example of FoxP2low new MSN
                           dendrite with mushroom (filled arrow), stubby (dashed arrow) and thin spines (unfilled arrow). (c) Example of
                           FoxP2high new MSNs dendrite. (d) New FoxP2high MSNs have more mushroom spines than FoxP2low new MSN
                           at 31 dpi. (e) Mushroom spine head size is not different between MSNs with different FoxP2 levels. (f) New
                           FoxP2low MSNs at 42 dpi show overall more dendritic spines and more thin spines than new FoxP2low MSNs at
                           31 dpi. Sample size (a,d,e,f): 52 MSNs of 4 zebra finches. Sample size (f): 23 MSNs of 3 zebra finches. *p ≤ 0.05;
                           ***p ≤ 0.001. Scale bars: 5 µm (b,c).

                               Since our previous work demonstrated that new MSNs participate in singing-associated neural activity in
                           Area X34 here we asked if this was linked to FoxP2 levels. After singing FoxP2 mRNA and protein levels are lower
                           in Area X tissue25–27 than when birds were silent whereas the expression of the immediate early gene EGR-1
                           increases linearly the more birds sing undirected song40,41. Analyzing expression levels in individual neurons
                           revealed that only FoxP2low MSNs, but not Foxp2high MSNs seemed to be subject to singing induced FoxP2 down-
                           regulation24. We therefore hypothesized that FoxP2high MSNs were not yet connected into the circuit and therefore
                           not regulated by singing. Our data contradict this hypothesis. We show that the induction of the immediate early
                           gene EGR-1 in MSNs after singing was equally likely in FoxP2low and FoxP2high MSNs, suggesting that both were
                           functionally incorporated into the song circuit.
                               Previous work showed that the degree of FoxP2 downregulation correlates with the quantity of produced song
                           and depends on auditory feedback, which is relayed to Area X via the cortical song nucleus HVC27,45. In our study,
                           the relationship between FoxP2 downregulation and song quantity was present in new MSNs at 42 dpi but not
                           earlier and we therefore suggest that new MSNs start to receive auditory feedback signals between 31 and 42 dpi.
                               What might cause the age-dependent FoxP2 downregulation in new MSNs? One possibility is that intrinsic
                           mechanisms, depending more on cell age than on extracellular inputs, downregulate FoxP2 during maturation.
                           Another possibility is that EGR-1 gradually decreases FoxP2 levels during every singing event. This latter scenario
                           is consistent with the findings that the FoxP2 promoter contains EGR-1 binding sites4,46 and that EGR-1 expres-
                           sion is crucial for functional integration of new neurons in the adult rodent hippocampus47.
                               We also tested if varying FoxP2 levels affect neuronal morphology in adult male zebra finches and analyzed
                           dendrite complexity and spine density of virally labelled new FoxP2low and FoxP2high MSNs. We now show dif-
                           ferences between the morphology of adult generated MSNs that express FoxP2 at high or low levels; high FoxP2
                           levels were associated with greater dendrite complexity and higher dendritic spine density in comparison to
                           neurons with low FoxP2 levels (Fig. 6). Concerning spine density, our results are consistent with previous studies

Scientific Reports |   (2020) 10:4787 | https://doi.org/10.1038/s41598-020-61740-6                                                          6
www.nature.com/scientificreports/                                                               www.nature.com/scientificreports

                           Figure 6. Graphical summary of the main results. The fraction of adult generated FoxP2high MSNs from zebra
                           finch Area X decreases with maturation. FoxP2high new MSNs show higher spine densities, more mushroom
                           spines and a more complex dendritic arborization than FoxP2low new MSNs.

                           in juvenile zebra finches, since experimental FoxP2 knockdown decreased overall spine densities of new Area X
                           MSNs48. Moreover, our results are in line with findings in mice that link Foxp2 to neuronal outgrowth and spine
                           density in striatal neurons and their progenitors18,19,35–37. Additionally, we find similarities between mice and birds
                           on the level of spine types. We show that FoxP2low new MSNs 31 dpi have fewer mushroom spines than FoxP2high
                           new MSNs. In mice, striatal spiny neurons of heterozygous Foxp2 knockdown mice show specifically a decrease
                           of mushroom and branched spines whereas stubby and thin spines are not affected19. In birds and mammals,
                           dendritic spines of striatal MSNs receive both glutamatergic dopaminergic input from cortical/pallial regions and
                           the midbrain, respectively. In the case of new MSNs in Area X, we hypothesize that high FoxP2 expression levels
                           during their maturation might increase their capacity for receiving both inputs.
                               What mechanisms might account for morphological differences between FoxP2low and FoxP2high MSNs? They
                           may originate from differential target gene activation in FoxP2low and FoxP2high MSNs. FoxP2 has hundreds of
                           downstream target genes of which many are part of networks associated with neurite development35,49–52. One
                           specific candidate is the myocyte enhancer factor 2C (Mef2c), a negative regulator of synaptogenesis19. Foxp2
                           specifically promotes corticostriatal synaptogenesis via the repression of Mef2c. Whether a similar mecha-
                           nism shapes the integration of new MSNs into the avian corticostriatal network remains to be elucidated by
                           future studies. In zebra finch Area X, two direct FoxP2 targets are associated with neuronal outgrowth; the
                           very-low-density-lipoprotein receptor (VLDLR) and the contactin-associated protein-like 2 (CNTNAP2). Their
                           expression correlates positively with FoxP2 in juveniles and in singing adults53,54. Thus, in new Area X neurons,
                           VLDLR and CNTNAP2 would be expected to be highly activated during singing in FoxP2high MSN but not in
                           FoxP2low MSNs and may thus generate the diverging MSNs morphology we found.
                               We would like to propose some speculations regarding possible functions of two MSNs subpopulations that
                           differ in FoxP2 expression levels. We found that these populations differed in nucleus size, dendritic complexity
                           and spine density during an early time period of their integration into Area X. We do not know if the observed
                           morphological differences persist long-term because of a lack of markers that could distinguish former FoxP2high
                           MSNs from former FoxP2low MSNs in later maturation phases. If these two subpopulations persist long-term, they
                           might resemble striato-nigral and striato-pallidal MSNs of the direct and indirect pathway in mammals. These
                           MSNs subtypes are morphologically and neurochemically different. Direct pathway MSNs express the dopa-
                           mine receptor D1 and their dendrites are more complex than indirect pathway MSNs that express the dopamine
                           receptor D239,55. High Foxp2 levels in D1 MSNs and low Foxp2 in D2 MSNs have been proposed to be linked to
                           this anatomical dichotomy35. The avian direct and indirect pathway through the basal ganglia however is not
                           characterized by different MSN projections but rather by direct and indirect pallidal-like output neurons that
                           project from Area X to the thalamus56. To date, it is not known if MSN subtypes exclusively synapse on either
                           direct or indirect pallidal-like neurons57. Contrary to mammalian MSNs, more than half of the Area X MSNs
                           express multiple dopamine receptors58 so that they cannot be used as markers for indirect versus direct pathway
                           neurons. Investigating potential avian MSNs subtypes and the developmental role FoxP2 plays in those will be of
                           interest for future studies.
                               What might be the function of new MSNs in Area X and how is it affected by FoxP2 expression lev-
                           els? Our previous study showed that once matured, new MSNs have similar characteristics as older, res-
                           ident MSNs and are active during singing34. General MSNs function is feed forward inhibition within a
                           cortico-striatal-thalamo-cortical loop during singing59. We hypothesize that new MSNs in adult zebra finches are
                           entrained to produce a correct firing pattern in a plastic phase during their maturation and thus may counteract

Scientific Reports |   (2020) 10:4787 | https://doi.org/10.1038/s41598-020-61740-6                                                              7
www.nature.com/scientificreports/                                                               www.nature.com/scientificreports

                           song drift, as has been suggested before33,60. This process might be in influenced by varying FoxP2 expression
                           levels and the resulting morphological differences between putative subpopulations of new MSNs. FoxP2high new
                           neurons with a higher dendrite complexity and more dendritic spines might be more receptive to tuning than
                           FoxP2low new neurons. For further interpretations of our findings it will be crucial to gain additional knowledge
                           on the microcircuitry of Area X and on the role neurons play for its function.
                              In summary, FoxP2 expression levels vary in adult-born MSN at different maturation times after they have
                           been recruited to Area X. We show that the different FoxP2 expression levels correlate with neuronal morphology
                           and spine density. Varying FoxP2 expression levels during a specific time window might permit different target
                           gene activation important for correct incorporation and function of new MSNs in Area X.

                           Methods
                           FoxP2 nomenclature. We follow the nomenclature proposed by61, FOXP refers to the human gene, Foxp
                           refers to the mouse gene and FoxP refers to all other species. FOXP, Foxp2 and FoxP2 correspond to the protein
                           product.

                           Animals.       42 adult male zebra finches (Taeniopygia guttata, age >120 days) were used in the present study,
                           bred and housed at the Department of Animal Behaviour at Freie Universität Berlin. The colony was kept under a
                           12:12 h light:dark-cycle with food and water ad libitum. All experiments were reviewed and approved by the vet-
                           erinary department of the Freie Universität Berlin and by the ethics committee of the Regional Office for Health
                           and Social Affairs Berlin and were performed in accordance with relevant guidelines and regulations. The permit
                           numbers are G0116/13 and G0296/15.

                           Experiments. We conduced three experiments. In the first, we analyzed FoxP2 expression levels of new neu-
                           rons (BrdU+, see below) in Area X at three time points, e.g. at 21 days (5 birds, 166 neurons), 31 days (6 birds,
                           295 neurons) or 42 days (6 birds, 272 neurons) after BrdU injections (dpi). In the second experiment we analyzed
                           FoxP2 expression levels and the expression of the early growth response protein 1 (EGR-1) at 21 dpi (6 birds, 127
                           neurons), 31 dpi (6 birds, 108 neurons) and at 42 dpi (6 birds 156 neurons). In the third experiment we analyzed
                           FoxP2 expression levels, dendrite morphology and spine density of new neurons that were labelled via lentiviral
                           infection. In total we analyzed 52 neurons of 4 zebra finches (13 ± 3 neurons/bird, mean ± SD) at 31 dpi and 23
                           neurons of 3 zebra finches at 42 dpi (7.6 ± 0.5 neurons/bird, mean ± SD).

                           BrdU injections. For the analysis of FoxP2 levels and EGR-1 expression in newborn neurons 35 birds
                           received BrdU (50 µg/g) via intramuscular injections in the mornings for 5 consecutive days. Birds were assigned
                           to three groups with different survival times (21, 31, and 42 days after BrdU injection, dpi).

                           Song monitoring. For FoxP2 expression level analysis after BrdU treatment or lentiviral injections, 17 birds
                           were isolated in sound attenuated chambers for one night before sacrifice. In the following morning, birds were
                           kept from singing by the experimenter sitting nearby for 1.5 h after lights went on. For EGR-1 expression analysis
                           in new neurons after singing, 18 birds were kept in sound attenuated chambers for three nights and were perfused
                           in the morning of the 4th day 1.5 h after the lights went on. During those 1.5 h birds had to sing at least 150 motifs
                           to be included in the subsequent analysis of EGR-1 expression. Vocalizations were continuously monitored via
                           Sound Analysis Pro62.

                           Lentiviral Vector injection. To label progenitor cells at the lateral wall of the ventricle, the lentiviral expres-
                           sion vector pFUGW63 containing a GFP reporter gene was generated as described in Lois et al.63 and stereotac-
                           tically injected into the ventricular zone of 7 birds under isoflurane anesthesia. Titers ranged from 2 × 106 and
                           3 × 107 viral particles/µl. Birds were fixed in a stereotaxic head holder, with the beak in a 45° angle from the
                           vertical axis. In each hemisphere, approximately 200 µl of virus containing solution were injected into four sites.
                           Following coordinates relative to the bifurcation of the midsagittal sinus were used: anterior-posterior 3.8–4.1,
                           medial-latera −1.3/+1.3, dorsal-ventral −5.0, injection angle: 10° lateral.

                           Immunohistochemistry.             For immunohistochemical staining birds were overdosed with isoflurane and
                           immediately perfused transcardially with phosphate-buffered saline (PBS) followed by 4% paraformaldehyde
                           (PFA) in PBS. Brains were dissected, post-fixed in 4% PFA for one night and washed for another night in PBS.
                           Brains were cut sagitally or frontally into 50 µm sections using a vibrating microtome (VT1000S, Leica). BrdU
                           antigen retrieval required incubation in 2N HCl for 30 min at 37 °C and neutralization with borate buffer. GFP
                           signal was enhanced via antibody staining. All immunostainings were performed according to standard proto-
                           cols. The following primary antibodies were used: anti FoxP2 (goat, Abcam ab1307, dilution:1:1000), anti EGR-1
                           (rabbit, Santa Cruz sc-189, dilution:1:600), anti BrdU (rat, Bio-Rad MCA2060, dilution: 1:200), anti GFP (rabbit,
                           Abcam ab290, dilution:1:1000). Fluorescent secondary antibodies were the following: anti-rat-Alexa-Fluor-488
                           (Life Technologies, A21208, dilution: 1:200), anti-rabbit-Alexa-Fluor-568 (Life Technologies, A10042, dilution:
                           1:200), anti-goat-Alaxa-Fluor-647 (Life Technologies, A21447, dilution: 1:200). To visualize nuclei, all sections
                           were counterstained with 4′, 6-Diamidin-2-phenylindol (DAPI, Serva).

                           Confocal imaging and image processing and quantification. Z-Stacks of BrdU+ or GFP+ cells
                           in Area X were obtained with a SP8 confocal microscope (Leica). For FoxP2 scanning, all microscope settings
                           were kept constant. Scans of BrdU+ nuclei were performed using a 63x lens (digital zoom 2.0), an image size of
                           1024 × 1024 pixels and a z-stack size of 1 µm. Whole neurons (GFP+) were imaged using a 63x lens, an image size
                           of 2042 × 2042 pixel and a z-stack size of 1 µm. Acquired images were processed using the Fiji software package64.
                           Only neurons with spiny long dendrites were included in the analysis. The Rolling Ball Background Subtraction

Scientific Reports |   (2020) 10:4787 | https://doi.org/10.1038/s41598-020-61740-6                                                             8
www.nature.com/scientificreports/                                                                          www.nature.com/scientificreports

                           plugin was used to subtract background. We measured the mean pixel intensity of nuclear FoxP2 expression, by
                           positioning a circle of 4 µm in diameter (12.56 µm2) in the center of the BrdU+ nucleus. In total we analyzed the
                           intensity of the FoxP2 expression dependent fluorescence of 166 BrdU+ cells at 21 dpi (n = 5), 295 BrdU+ cells
                           at 31 dpi (n = 6) and 272 BrdU+ cells at 42 dpi (n = 6). FoxP2high were defined as cells that fell into the top 30% of
                           measured mean pixel intensities in one animal (i.e. if the highest mean pixel intensity in one animal was 240 we
                           counted all BrdU+ cells that had a FoxP2 mean pixel intensity between 168 and 240 as FoxP2high neurons). We
                           decided on the 30% value because it covered the FoxP2high expressing cells in the bimodal distribution of all FoxP2
                           intensities. We defined the neurons that fell into the bottom 30% of measured mean pixel intensities as FoxP2low.
                           As for the FoxP2high cutoff, the bottom 30% contained the low-intensity peak of the bimodal density distribution.
                           Because we particularly wanted to address the effect of high and low FoxP2 expression levels on neuronal prop-
                           erties, neurons with intermediate FoxP2 expression levels were not considered for further analysis. The Simple
                           Neurite Tracer plugin (Fiji) was used to trace individual neurons and we analyzed their total branch length and
                           number of primary dendrites. The traces were then used by the Sholl analysis plugin in Fiji65. We measured inter-
                           sections of dendrites with concentric circles that were placed every 10 µm starting from the center of the soma.
                           The maximal number of intersections per neuron was extracted from the Sholl analysis dataset. For dendritic
                           spine analysis images were deconvolved using the Tikhonov-Miller algorithm in the DeconvolutionLab plugin
                           in Fiji66. Prior to deconvolution an individual point spread function was generated for each image using the Born
                           and Wolf 3D optical model in the PSF Generator plugin in Fiji67. Semiautomated dendritic spine counts were per-
                           formed using the software NeuronStudio68 that uses a spine classification algorithm. For spine classification, the
                           default settings were used to classify spines as mushroom, stubby or thin spines: a head-to-neck ratio threshold
                           of 1.1 µm, a height-to-width ratio threshold of 2.5 µm and a minimum mushroom head size of 0.35 µm. A spine is
                           considered mushroom if the head-to-neck ratio is above the threshold and its head is larger than 0.35 µm. A spine
                           is considered stubby if its head-to-neck ratio and its heights-to-width ratio are below threshold. In all other cases
                           spines were classified as thin. On average, we analyzed spines densities on secondary dendrites along the length
                           of 118 µm ± 1.92 (mean, SEM) per neuron. In total, we analyzed spines of 52 individual neurons of 4 animals in
                           experimental group 31 dpi, and 23 neurons of 3 animals at 42 dpi. Additionally, we measured the dendrite diam-
                           eter of 44 new neurons in 4 animals (8–12 neurons per animal) using the line measuring tool in Fiji64. We took 5
                           measures on each of 3 secondary dendrites per neuron (in total 15 measures per neuron). The experimenter was
                           blind to FoxP2 levels of individual neurons during the whole quantification process, because cells were selected
                           for quantification based on their BrdU+ fluorescence or their EGR-1 fluorescence and FoxP2 fluorescence in a
                           different channel was quantified last. The datasets generated and analysed during the current study are available
                           from the corresponding author on request.

                           Statistics. The software R was used to analyze data69. Significance level was p < 0.05 for all tests. Plots were
                           generated using the ggplot package in R70. For the analysis of FoxP2high neurons we used a Kruskal-Wallis test fol-
                           lowed by a Dunn’s test for pairwise comparison. For the analysis of FoxP2low neurons we used ANOVA followed
                           by a Bonferroni’s multiple comparison test. The relationship of (a) FoxP2 levels and nuclear diameter as well as (b)
                           FoxP2+ neurons and singing were determined using a linear regression analysis. Dendritic spine data (all spines)
                           and Sholl data were analyzed using the paired Student’s t-test. Data from the spine type analysis, mushroom head
                           size, total dendrite length, number of primary dendrites and dendrite thickness were analyzed using the Mann
                           Whitney test. Number of maximal intersections was analyzed using the Student’s t-test. Choice of test was based
                           on previous analysis for normality using the Shapiro-Wilk-test and variance analysis using the F-test or Levene’s
                           test.

                           Received: 29 October 2019; Accepted: 29 February 2020;
                           Published: xx xx xxxx

                           References
                            1. Lai, C. S., Fisher, S. E., Hurst, J. A., Vargha-Khadem, F. & Monaco, A. P. A forkhead-domain gene is mutated in a severe speech and
                               language disorder. Nature 413, 519–523, https://doi.org/10.1038/35097076 (2001).
                            2. Deriziotis, P. & Fisher, S. E. Speech and Language: Translating the Genome. Trends Genet 33, 642–656, https://doi.org/10.1016/j.
                               tig.2017.07.002 (2017).
                            3. Morgan, A., Fisher, S. E., Scheffer, I. & Hildebrand, M. In GeneReviews (eds Adam, M. P. et al.) (2016).
                            4. Wohlgemuth, S., Adam, I. & Scharff, C. FoxP2 in songbirds. Curr Opin Neurobiol 28, 86–93, https://doi.org/10.1016/j.
                               conb.2014.06.009 (2014).
                            5. Teramitsu, I., Kudo, L. C., London, S. E., Geschwind, D. H. & White, S. A. Parallel FoxP1 and FoxP2 expression in songbird and
                               human brain predicts functional interaction. J Neurosci 24, 3152–3163, https://doi.org/10.1523/JNEUROSCI.5589-03.2004 (2004).
                            6. Lai, C. S., Gerrelli, D., Monaco, A. P., Fisher, S. E. & Copp, A. J. FOXP2 expression during brain development coincides with adult
                               sites of pathology in a severe speech and language disorder. Brain 126, 2455–2462, https://doi.org/10.1093/brain/awg247 (2003).
                            7. Haesler, S. et al. FoxP2 expression in avian vocal learners and non-learners. J Neurosci 24, 3164–3175, https://doi.org/10.1523/
                               JNEUROSCI.4369-03.2004 (2004).
                            8. Heston, J. B. & White, S. A. Behavior-Linked FoxP2 Regulation Enables Zebra Finch Vocal Learning. J Neurosci 35, 2885–2894,
                               https://doi.org/10.1523/JNEUROSCI.3715-14.2015 (2015).
                            9. Haesler, S. et al. Incomplete and inaccurate vocal imitation after knockdown of FoxP2 in songbird basal ganglia nucleus Area X.
                               PLoS Biol 5, e321, https://doi.org/10.1371/journal.pbio.0050321 (2007).
                           10. Norton, P., Barschke, P., Scharff, C. & Mendoza, E. Differential Song Deficits after Lentivirus-Mediated Knockdown of FoxP1,
                               FoxP2, or FoxP4 in Area X of Juvenile Zebra Finches. J Neurosci 39, 9782–9796, https://doi.org/10.1523/JNEUROSCI.1250-19.2019
                               (2019).
                           11. Vargha-Khadem, F. et al. Neural basis of an inherited speech and language disorder. Proc Natl Acad Sci USA 95, 12695–12700,
                               https://doi.org/10.1073/pnas.95.21.12695 (1998).
                           12. Day, N. F., Hobbs, T. G., Heston, J. B. & White, S. A. Beyond Critical Period Learning: Striatal FoxP2 Affects the Active Maintenance
                               of Learned Vocalizations in Adulthood. eNeuro 6, https://doi.org/10.1523/ENEURO.0071-19.2019 (2019).

Scientific Reports |   (2020) 10:4787 | https://doi.org/10.1038/s41598-020-61740-6                                                                                9
www.nature.com/scientificreports/                                                                            www.nature.com/scientificreports

                           13. Murugan, M., Harward, S., Scharff, C. & Mooney, R. Diminished FoxP2 levels affect dopaminergic modulation of corticostriatal
                               signaling important to song variability. Neuron 80, 1464–1476, https://doi.org/10.1016/j.neuron.2013.09.021 (2013).
                           14. Tsui, D., Vessey, J. P., Tomita, H., Kaplan, D. R. & Miller, F. D. FoxP2 regulates neurogenesis during embryonic cortical development.
                               J Neurosci 33, 244–258, https://doi.org/10.1523/JNEUROSCI.1665-12.2013 (2013).
                           15. Garcia-Calero, E., Botella-Lopez, A., Bahamonde, O., Perez-Balaguer, A. & Martinez, S. FoxP2 protein levels regulate cell
                               morphology changes and migration patterns in the vertebrate developing telencephalon. Brain Struct Funct 221, 2905–2917, https://
                               doi.org/10.1007/s00429-015-1079-7 (2016).
                           16. Clovis, Y. M., Enard, W., Marinaro, F., Huttner, W. B. & De Pietri Tonelli, D. Convergent repression of Foxp2 3′UTR by miR-9 and
                               miR-132 in embryonic mouse neocortex: implications for radial migration of neurons. Development 139, 3332–3342, https://doi.
                               org/10.1242/dev.078063 (2012).
                           17. Kast, R. J., Lanjewar, A. L., Smith, C. D. & Levitt, P. FOXP2 exhibits projection neuron class specific expression, but is not required
                               for multiple aspects of cortical histogenesis. Elife 8, https://doi.org/10.7554/eLife.42012 (2019).
                           18. Chiu, Y. C. et al. Foxp2 regulates neuronal differentiation and neuronal subtype specification. Dev Neurobiol 74, 723–738, https://
                               doi.org/10.1002/dneu.22166 (2014).
                           19. Chen, Y. C. et al. Foxp2 controls synaptic wiring of corticostriatal circuits and vocal communication by opposing Mef2c. Nat
                               Neurosci, https://doi.org/10.1038/nn.4380 (2016).
                           20. Groszer, M. et al. Impaired synaptic plasticity and motor learning in mice with a point mutation implicated in human speech
                               deficits. Curr Biol 18, 354–362, https://doi.org/10.1016/j.cub.2008.01.060 (2008).
                           21. French, C. A. et al. An aetiological Foxp2 mutation causes aberrant striatal activity and alters plasticity during skill learning. Mol
                               Psychiatry 17, 1077–1085, https://doi.org/10.1038/mp.2011.105 (2012).
                           22. French, C. A. et al. Differential effects of Foxp2 disruption in distinct motor circuits. Mol Psychiatry 24, 447–462, https://doi.
                               org/10.1038/s41380-018-0199-x (2019).
                           23. van Rhijn, J. R., Fisher, S. E., Vernes, S. C. & Nadif Kasri, N. Foxp2 loss of function increases striatal direct pathway inhibition via
                               increased GABA release. Brain Struct Funct 223, 4211–4226, https://doi.org/10.1007/s00429-018-1746-6 (2018).
                           24. Thompson, C. K. et al. Young and intense: FoxP2 immunoreactivity in Area X varies with age, song stereotypy, and singing in male
                               zebra finches. Front Neural Circuits 7, 24, https://doi.org/10.3389/fncir.2013.00024 (2013).
                           25. Teramitsu, I. & White, S. A. FoxP2 regulation during undirected singing in adult songbirds. J Neurosci 26, 7390–7394, https://doi.
                               org/10.1523/JNEUROSCI.1662-06.2006 (2006).
                           26. Miller, J. E. et al. Birdsong decreases protein levels of FoxP2, a molecule required for human speech. J Neurophysiol 100, 2015–2025,
                               https://doi.org/10.1152/jn.90415.2008 (2008).
                           27. Teramitsu, I., Poopatanapong, A., Torrisi, S. & White, S. A. Striatal FoxP2 is actively regulated during songbird sensorimotor
                               learning. PLoS One 5, e8548, https://doi.org/10.1371/journal.pone.0008548 (2010).
                           28. Alvarez-Buylla, A., Kirn, J. R. & Nottebohm, F. Birth of projection neurons in adult avian brain may be related to perceptual or
                               motor learning. Science 249, 1444–1446 (1990).
                           29. Alvarez-Buylla, A., Ling, C. Y. & Yu, W. S. Contribution of neurons born during embryonic, juvenile, and adult life to the brain of
                               adult canaries: regional specificity and delayed birth of neurons in the song-control nuclei. J Comp Neurol 347, 233–248, https://doi.
                               org/10.1002/cne.903470207 (1994).
                           30. Lipkind, D., Nottebohm, F., Rado, R. & Barnea, A. Social change affects the survival of new neurons in the forebrain of adult
                               songbirds. Behav Brain Res 133, 31–43, https://doi.org/10.1016/s0166-4328(01)00416-8 (2002).
                           31. Rochefort, C., He, X., Scotto-Lomassese, S. & Scharff, C. Recruitment of FoxP2-expressing neurons to area X varies during song
                               development. Dev Neurobiol 67, 809–817, https://doi.org/10.1002/dneu.20393 (2007).
                           32. Barnea, A. & Pravosudov, V. Birds as a model to study adult neurogenesis: bridging evolutionary, comparative and neuroethological
                               approaches. Eur J Neurosci 34, 884–907, https://doi.org/10.1111/j.1460-9568.2011.07851.x (2011).
                           33. Pytte, C. L. Adult Neurogenesis in the Songbird: Region-Specific Contributions of New Neurons to Behavioral Plasticity and
                               Stability. Brain Behav Evol 87, 191–204, https://doi.org/10.1159/000447048 (2016).
                           34. Kosubek-Langer, J., Schulze, L. & Scharff, C. Maturation, Behavioral Activation, and Connectivity of Adult-Born Medium Spiny
                               Neurons in a Striatal Song Nucleus. Front Neurosci 11, 323, https://doi.org/10.3389/fnins.2017.00323 (2017).
                           35. Vernes, S. C. et al. Foxp2 regulates gene networks implicated in neurite outgrowth in the developing brain. PLoS Genet 7, e1002145,
                               https://doi.org/10.1371/journal.pgen.1002145 (2011).
                           36. Enard, W. et al. A humanized version of Foxp2 affects cortico-basal ganglia circuits in mice. Cell 137, 961–971, https://doi.
                               org/10.1016/j.cell.2009.03.041 (2009).
                           37. Reimers-Kipping, S., Hevers, W., Paabo, S. & Enard, W. Humanized Foxp2 specifically affects cortico-basal ganglia circuits.
                               Neuroscience 175, 75–84, https://doi.org/10.1016/j.neuroscience.2010.11.042 (2011).
                           38. Stanley, G., Gokce, O., Malenka, R. C., Sudhof, T. C. & Quake, S. R. Continuous and Discrete Neuron Types of the Adult Murine
                               Striatum. Neuron, https://doi.org/10.1016/j.neuron.2019.11.004 (2019).
                           39. Gertler, T. S., Chan, C. S. & Surmeier, D. J. Dichotomous anatomical properties of adult striatal medium spiny neurons. J Neurosci
                               28, 10814–10824, https://doi.org/10.1523/JNEUROSCI.2660-08.2008 (2008).
                           40. Jarvis, E. D., Scharff, C., Grossman, M. R., Ramos, J. A. & Nottebohm, F. For whom the bird sings: context-dependent gene
                               expression. Neuron 21, 775–788 (1998).
                           41. Mello, C. V. & Ribeiro, S. ZENK protein regulation by song in the brain of songbirds. J Comp Neurol 393, 426–438 (1998).
                           42. Knapska, E. & Kaczmarek, L. A gene for neuronal plasticity in the mammalian brain: Zif268/Egr-1/NGFI-A/Krox-24/TIS8/ZENK?
                               Prog Neurobiol 74, 183–211, https://doi.org/10.1016/j.pneurobio.2004.05.007 (2004).
                           43. Hessler, N. A. & Doupe, A. J. Social context modulates singing-related neural activity in the songbird forebrain. Nat Neurosci 2,
                               209–211, https://doi.org/10.1038/6306 (1999).
                           44. Zengin-Toktas, Y. & Woolley, S. C. Singing modulates parvalbumin interneurons throughout songbird forebrain vocal control
                               circuitry. PLoS One 12, e0172944, https://doi.org/10.1371/journal.pone.0172944 (2017).
                           45. Schmidt, M. F. & Konishi, M. Gating of auditory responses in the vocal control system of awake songbirds. Nat Neurosci 1, 513–518,
                               https://doi.org/10.1038/2232 (1998).
                           46. Warren, W. C. et al. The genome of a songbird. Nature 464, 757–762, https://doi.org/10.1038/nature08819 (2010).
                           47. Veyrac, A. et al. Zif268/egr1 gene controls the selection, maturation and functional integration of adult hippocampal newborn
                               neurons by learning. Proc Natl Acad Sci USA 110, 7062–7067, https://doi.org/10.1073/pnas.1220558110 (2013).
                           48. Schulz, S. B., Haesler, S., Scharff, C. & Rochefort, C. Knockdown of FoxP2 alters spine density in Area X of the zebra finch. Genes
                               Brain Behav 9, 732–740, https://doi.org/10.1111/j.1601-183X.2010.00607.x (2010).
                           49. Spiteri, E. et al. Identification of the transcriptional targets of FOXP2, a gene linked to speech and language, in developing human
                               brain. Am J Hum Genet 81, 1144–1157, https://doi.org/10.1086/522237 (2007).
                           50. Vernes, S. C. et al. High-throughput analysis of promoter occupancy reveals direct neural targets of FOXP2, a gene mutated in
                               speech and language disorders. Am J Hum Genet 81, 1232–1250, https://doi.org/10.1086/522238 (2007).
                           51. Hickey, S. L., Berto, S. & Konopka, G. Chromatin Decondensation by FOXP2 Promotes Human Neuron Maturation and Expression
                               of Neurodevelopmental Disease Genes. Cell Rep 27, 1699–1711 e1699, https://doi.org/10.1016/j.celrep.2019.04.044 (2019).
                           52. Konopka, G. et al. Human-specific transcriptional regulation of CNS development genes by FOXP2. Nature 462, 213–217, https://
                               doi.org/10.1038/nature08549 (2009).

Scientific Reports |   (2020) 10:4787 | https://doi.org/10.1038/s41598-020-61740-6                                                                                 10
www.nature.com/scientificreports/                                                                              www.nature.com/scientificreports

                           53. Adam, I., Mendoza, E., Kobalz, U., Wohlgemuth, S. & Scharff, C. FoxP2 directly regulates the reelin receptor VLDLR developmentally
                               and by singing. Mol Cell Neurosci 74, 96–105, https://doi.org/10.1016/j.mcn.2016.04.002 (2016).
                           54. Adam, I., Mendoza, E., Kobalz, U., Wohlgemuth, S. & Scharff, C. CNTNAP2 is a direct FoxP2 target in vitro and in vivo in zebra
                               finches: complex regulation by age and activity. Genes Brain Behav 16, 635–642, https://doi.org/10.1111/gbb.12390 (2017).
                           55. Calabresi, P., Picconi, B., Tozzi, A., Ghiglieri, V. & Di Filippo, M. Direct and indirect pathways of basal ganglia: a critical reappraisal.
                               Nat Neurosci 17, 1022–1030, https://doi.org/10.1038/nn.3743 (2014).
                           56. Farries, M. A., Ding, L. & Perkel, D. J. Evidence for “direct” and “indirect” pathways through the song system basal ganglia. J Comp
                               Neurol 484, 93–104, https://doi.org/10.1002/cne.20464 (2005).
                           57. Pidoux, M., Bollu, T., Riccelli, T. & Goldberg, J. H. Origins of basal ganglia output signals in singing juvenile birds. J Neurophysiol
                               113, 843–855, https://doi.org/10.1152/jn.00635.2014 (2015).
                           58. Kubikova, L., Wada, K. & Jarvis, E. D. Dopamine receptors in a songbird brain. J Comp Neurol 518, 741–769, https://doi.org/10.1002/
                               cne.22255 (2010).
                           59. Perkel, D. J., Farries, M. A., Luo, M. & Ding, L. Electrophysiological analysis of a songbird basal ganglia circuit essential for vocal
                               plasticity. Brain Res Bull 57, 529–532, https://doi.org/10.1016/s0361-9230(01)00690-6 (2002).
                           60. Wilbrecht, L. & Kirn, J. R. Neuron addition and loss in the song system: regulation and function. Ann N Y Acad Sci 1016, 659–683,
                               https://doi.org/10.1196/annals.1298.024 (2004).
                           61. Kaestner, K. H., Knochel, W. & Martinez, D. E. Unified nomenclature for the winged helix/forkhead transcription factors. Genes Dev
                               14, 142–146 (2000).
                           62. Tchernichovski, O., Nottebohm, F., Ho, C. E., Pesaran, B. & Mitra, P. P. A procedure for an automated measurement of song
                               similarity. Anim Behav 59, 1167–1176, https://doi.org/10.1006/anbe.1999.1416 (2000).
                           63. Lois, C., Hong, E. J., Pease, S., Brown, E. J. & Baltimore, D. Germline transmission and tissue-specific expression of transgenes
                               delivered by lentiviral vectors. Science 295, 868–872, https://doi.org/10.1126/science.1067081 (2002).
                           64. Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. Nat Methods 9, 676–682, https://doi.org/10.1038/
                               nmeth.2019 (2012).
                           65. Ferreira, T. A. et al. Neuronal morphometry directly from bitmap images. Nat Methods 11, 982–984, https://doi.org/10.1038/
                               nmeth.3125 (2014).
                           66. Sage, D. et al. DeconvolutionLab2: An open-source software for deconvolution microscopy. Methods 115, 28–41, https://doi.
                               org/10.1016/j.ymeth.2016.12.015 (2017).
                           67. Kirshner, H., Aguet, F., Sage, D. & Unser, M. 3-D PSF fitting for fluorescence microscopy: implementation and localization
                               application. J Microsc 249, 13–25, https://doi.org/10.1111/j.1365-2818.2012.03675.x (2013).
                           68. Rodriguez, A., Ehlenberger, D. B., Dickstein, D. L., Hof, P. R. & Wearne, S. L. Automated three-dimensional detection and shape
                               classification of dendritic spines from fluorescence microscopy images. PLoS One 3, e1997, https://doi.org/10.1371/journal.
                               pone.0001997 (2008).
                           69. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2013).
                           70. Wickham, H. gglot2: Elegant Graphics for Data Analysis. (Springer, 2016).

                           Acknowledgements
                           J.K.L. received funds from the Elsa Neumann scholarship by the state of Berlin, the Christiane Nüsslein-Volhard
                           Foundation and the Max Planck School of Cognition.

                           Author contributions
                           J.K.L. and C.S. designed experiments, J.K.L. performed experiments and data analysis, J.K.L. and C.S. wrote the
                           manuscript.

                           Competing interests
                           The authors declare no competing interests.

                           Additional information
                           Correspondence and requests for materials should be addressed to J.K.-L.
                           Reprints and permissions information is available at www.nature.com/reprints.
                           Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and
                           institutional affiliations.
                                         Open Access This article is licensed under a Creative Commons Attribution 4.0 International
                                         License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
                           format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
                           ative Commons license, and indicate if changes were made. The images or other third party material in this
                           article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the
                           material. If material is not included in the article’s Creative Commons license and your intended use is not per-
                           mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
                           copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

                           © The Author(s) 2020

Scientific Reports |   (2020) 10:4787 | https://doi.org/10.1038/s41598-020-61740-6                                                                                     11
You can also read