Introduction

Evolutionary changes that occurred specifically in the human genome have been proposed to be significant drivers of human speciation. Recently, significant progress has been made in cataloging human-specific genetic changes, especially with respect to evolution of the human brain1. However, knowledge of how these human-specific genetic modifications impact corresponding molecular pathways is still lacking. For example, multiple human-specific gene duplications (HSGDs) have recently been identified to play a role in controlling human brain development2,3,4,5,6,7,8,9,10,11, but how these genetic modifiers affect specific signaling pathways in the developing or adult human brain remains poorly understood. This limits our understanding of how these changes led to the emergence of human-specific biological traits.

Previous work has demonstrated how during human evolution, gene duplication of SLIT-ROBO Rho-GTPase-activating protein 2 (SRGAP2A), which is highly expressed in the developing brain of all mammals12,13, resulted in emergence of the human-specific paralogs SRGAP2B and SRGAP2C8,14. Both are truncated open reading frames (including only the first 9 exons out of 22 in Srgap2a) leading to the expression of a truncated protein corresponding to the extended F-BAR domain (F-BARx), that lacks the last C-terminal 49 amino acids (F-BARΔ49) of SRGAP2A’s FBARx domain and replaces them with 7 unique C-terminal amino acids (Fig. 1B). In addition, SRGAP2C acquired a series of unique nonsynonymous base pair mutations selectively targeting five arginine residues compared to SRGAP2B6,8 (Fig. 1C). This truncation and these specific arginine mutations have been shown to reduce solubility of SRGAP2C and increase its ability to heterodimerize with SRGAP2A to form an insoluble complex15.

Figure 1
figure 1

SRGAP2A and the emergence of the human-specific paralogs SRGAP2B and SRGAP2C. (A) The ancestral copy SRGAP2A, present in rodents, primates, and humans, was duplicated in the human lineage to form the truncated copies SRGAP2B and SRGAP2C. (B) SRGAP2A contains three distinct protein domains. Through binding of Homer1 to an EVH1 site in the F-BAR domain, SRGAP2A promotes maturation of excitatory synapses, while binding of Gephyrin to the SH3 domain promotes inhibitory synapse maturation. The Rho-GAP domain is involved in limiting synapse density through Rac1. (C) Duplication of SRGAP2A in humans generated a truncated protein, F-BARΔ49, containing the F-BARx domain of SRGAP2A that lacks the last C-terminal 49 amino acids. Subsequent nonsynonymous base pair mutations converted arginines to non-polar residues which led to the formation of SRGAP2B and SRGAP2C. Binding of SRGAP2C to SRGAP2A inhibits all functions of SRGAP2A. (D) In rodents, synaptic development of cortical pyramidal neurons is complete at P21. At this stage, both synapse size and synaptic density reach levels of those observed in adulthood. In contrast, expression of SRGAP2C in mouse cortical pyramidal neurons binds and inhibits SRGAP2A resulting in protracted synaptic maturation and increased synaptic density.

Expression of SRGAP2C in mouse cortical pyramidal neurons leads to the emergence of human-specific traits of synaptic development, including an increase in synaptic density and protracted synaptic maturation of both excitatory and inhibitory synapses6,7. This phenotype is strikingly similar to that of a constitutive knockout mouse in which Srgap2a expression is genomically reduced, strongly suggesting that human-specific SRGAP2C functions largely by inhibiting SRGAP2A function6. Interestingly, the ancestral copy SRGAP2A limits excitatory (E) and inhibitory (I) synapse density through its Rac1-specific GAP domain, while promoting maturation of both E and I synapses through its ability to bind to the postsynaptic scaffolding protein Homer1 at E synapses through its class II EVH1 binding domain embedded in its F-BAR domain, and the postsynaptic protein Gephyrin through its SH3 domain at I synapses6,7,16,17 (Fig. 1A). Because SRGAP2C directly binds SRGAP2A through its truncated F-BARx domain, we previously hypothesized that this binding directly inhibits the function of SRGAP2A6. However, the mechanisms underlying the ability of SRGAP2C to inhibit all functions of SRGAP2A remained unknown.

Recent evidence suggests that both SRGAP2B and SRGAP2C are expressed in the human brain2. In the human population, SRGAP2C copy numbers are remarkably fixed, while SRGAP2B exhibits significant copy number variation (CNV) in the human population8. This suggests that in contrast to SRGAP2B, SRGAP2C has been rapidly fixed in the human population since its emergence during human brain evolution, and suggests that SRGAP2C might have played a unique role compared to SRGAP2B. However, the human-specific paralog SRGAP2B has not been characterized functionally. SRGAP2B is almost identical to SRGAP2C in all other respects: like SRGAP2C, it lacks the last 49 amino acids of its F-BARx domain and shares the same 7 unique C-terminal amino acids, but has its own unique point mutations that differ from those of SRGAP2C6,8.

Here we provide evidence showing that in cortical neurons, the truncation of the F-BARx domain present in both SRGAP2B and SRGAP2C leads to proteasome-mediated degradation of these proteins. Moreover, upon binding to SRGAP2A, SRGAP2C targets this hetero-dimer to the proteasome degradation pathway, thereby effectively reducing SRGAP2A protein levels in dendrites of cortical PNs. We show that SRGAP2C is uniquely more potent than SRGAP2B at a long-lasting increase of synaptic density into adult cortical PNs. Together, these results show how the emergence of the human-specific paralog SRGAP2C directly impacted the ancestral copy SRGAP2A, a critical regulator of synaptic development during human brain evolution.

Results

Human-specific partial duplication of SRGAP2A resulted in the emergence of two truncated paralogs (SRGAP2B and SRGAP2C) which both encode for a truncated protein consisting of the F-BARx domain of SRGAP2A lacking the last C-terminal 49 amino acids6,8 (F-BARΔ49; Fig. 1). In addition, SRGAP2C acquired a series of point mutations that led to five non-synonymous amino acid changes that in SRGAP2A are all arginine residues (Fig. 1). First, we examined the effects of the Δ49 truncation. Previous work using recombinant proteins and cell lines showed that this Δ49 truncation of the F-BARx domain generates intrinsically unstable and insoluble protein aggregates15. In order to investigate whether the Δ49 truncation caused the formation of similar aggregates in neurons, we generated a construct in which we reintroduced these 49 amino acids (full-length F-BARx domain) and compared it to F-BARx with the Δ49 truncation (F-BARΔ49), and SRGAP2C. We co-expressed mRFP-tagged versions of these constructs together with Venus (cell fill) by ex utero electroporation in mouse cortical pyramidal neurons (PNs) cultured for 18–21 days in vitro (DIV) (Figs. 2 and S1). Expression of F-BARx protein was clearly observed in soma and dendrites of cortical PNs. In contrast, F-BARΔ49 and SRGAP2C expression levels were generally very low (Figs. 2A and S2). We wondered whether the low expression of F-BARΔ49 and SRGAP2C was the result of an active degradation mechanism, such as the proteasome degradation pathway, in response to the insolubility of these proteins. We therefore treated with MG-132, a potent and widely used inhibitor of proteasome function, and used a live-cell confocal imaging approach to measure protein levels in the same neurons over time. Upon treatment with MG-132, we observed a rapid increase of F-BARΔ49 and SRGAP2C protein levels, while F-BARx levels increased only slightly (Fig. 2A,B). These results show that the Δ49 truncation generates an insoluble protein that in neurons leads to proteasome-mediated degradation.

Figure 2
figure 2

Proteasome degradation of SRGAP2 proteins. (A) Expression of RFP-tagged F-BARx, F-BARΔ49, or SRGAP2C in mouse cortical pyramidal neurons cultured for 18–21 days in vitro (DIV). Low levels of expression were observed for F-BARΔ49 and SRGAP2C, while addition of the proteasome inhibitor MG-132 resulted in a strong increase after 8 h of treatment. F-BARx is highly expressed throughout the experiment. Scale bar top panels: 25 µm, bottom panels: 5 µm. (B) Quantification of fluorescence intensity after MG-132 treatment. nF-BARx = 31 neurons, nF-BARΔ49 = 50 neurons, nSRGAP2C = 47 neurons; multiple t tests with Holm-Sidak correction for multiple comparisons (α = 0.05); ****p < 0.0001; mean ± SEM. (C–E) Co-expression of SRGAP2A and SRGAP2C in cultured mouse cortical neurons. SRGAP2A levels are reduced when co-expressed with SRGAP2C in both soma and dendritic spines. Scale bars: 5 µm. (E) Quantification of SRGAP2A fluorescence intensity. Soma: nControl = 110, nSRGAP2C = 68; Spines: nControl = 969, nSRGAP2C = 1259; Mann-Whitney test; ***p < 0.001; mean ± SEM. (F) Treatment of cultured mouse cortical neurons expressing both SRGAP2A-GFP and SRGAP2C-RFP with MG-132. Levels for both SRGAP2A and SRGAP2C increased upon treatment, and co-localization of SRGAP2A and SRGAP2C protein clusters were observed (red arrow). Scale bar: 10 µm. (G) Quantification of SRGAP2A fluorescent intensity when co-expressed with either F-BARx or SRGAP2C in the presence of MG-132. nF-BAR = 60 neurons, nSRGAP2C = 59 neurons; multiple t tests with Holm-Sidak correction for multiple comparisons (α = 0.05); ***p < 0.001; mean ± SEM.

Considering that SRGAP2C directly binds SRGAP2A, and that SRGAP2A plays a critical role in synaptic development6,7, we next asked how SRGAP2C expression affects SRGAP2A protein. We co-expressed EGFP-tagged or mTagBFP-HA-tagged SRGAP2A with or without mRFP-tagged SRGAP2C in mouse cortical neurons. Additionally, we co-expressed a cell fill, either Venus or mTagBFP-HA, for visualization of neuronal morphology. Interestingly, expression levels for SRGAP2A were significantly reduced in both soma (Fig. 2C) and dendritic spines (Fig. 2D) upon co-expression with SRGAP2C (Fig. 2E).

Since SRGAP2C binds to SRGAP2A through its F-BARx domain6,15, we next asked whether the reduction of SRGAP2A levels in response to ‘humanization’ of SRGAP2C expression is proteasome-dependent, i.e., if the mechanism leading to reduction of SRGAP2A function in presence of SRGAP2C is due to targeting of this heterodimer complex to the proteasome for degradation. We co-expressed SRGAP2A-GFP with either F-BARx-RFP or SRGAP2C-RFP and treated neurons with MG-132. While SRGAP2A levels were only slightly elevated upon treatment with MG-132 when co-expressed with the full-length F-BARx domain, we observed a significant increase when co-expressed with SRGAP2C (Fig. 2F,G) and frequently observed the formation and co-localization of SRGAP2A and SRGAP2C protein clusters (Fig. 2F).

We performed a series of similar experiments using the protein synthesis inhibitor Cycloheximide (Chx) to measure the rate of SRGAP2A protein turnover when co-expressed with SRGAP2C. Indeed, Chx treatment led to a reduction of protein fluorescence over time. Compared to F-BARx, co-expression with SRGAP2C significantly decreased the half-life of SRGAP2A (Fig. S3). Together, these observations show that the Δ49 truncation present in SRGAP2C is sufficient to induce proteasome-mediated degradation of both SRGAP2C and SRGAP2A when co-expressed together.

Besides the Δ49 truncation, SRGAP2C also acquired a number of point mutations resulting in non-synonymous changes of five arginines. These mutations are highly specific to SRGAP2C (not present in SRGAP2B) and have most likely been under selective pressure over the course of human evolution8. Furthermore, both SRGAP2B and C are expressed during human brain development2,8. Finally, because SRGAP2B displays significantly more copy number variation (CNVs) in the human population than SRGAP2C, it has been proposed that these two HSGDs might have different functions and, most importantly, might exert different strengths in their ability to act as modifiers of SRGAP2A function. To test whether the arginine mutations specific to SRGAP2C contribute to its ability to regulate synaptic development differentially than SRGAP2B in vivo, we used in utero cortical electroporation (IUCE) to express Venus (cell filler) and mRFP-tagged fusion protein of either F-BARx, F-BARΔ49, SRGAP2B, and SRGAP2C in layer 2/3 pyramidal cortical neurons in vivo (Figs. 3A, S4, and Table S1) and compared their abilities to phenocopy a partial SRGAP2A loss-of-function6. We analyzed the ability of these proteins to regulate excitatory synaptic density and maturation by quantifying spine density and spine head size of apical oblique dendrites at juvenile (P21 – Fig. 3) and adult (>P65 – Fig. 4) stages. Comparing the effect of SRGAP2C on synaptic development to F-BARx and F-BARΔ49 allowed us to determine how the distinct evolutionary events, from the truncation of the F-BARx domain (present in both SRGAP2B and C) to the progressive acquisition of five arginine conversions (present only in SRGAP2C), contributed to human-specific evolution of synaptic development.

Figure 3
figure 3

In vivo analysis of spine density and spine size at juvenile stage. (A) Schematic illustrating in utero electroporation approach targeting layer 2/3 cortical pyramidal neurons and subsequent quantification of dendritic spine density and head size. Scale bar: 50 µm. (B) Representative images showing apical oblique dendrites of neurons at P21 co-expressing either Venus and tdTomato (Control), or mRFP-tagged F-BARx, F-BARΔ49, SRGAP2B, or SRGAP2C. Scale bar: 2 µm. (C,D) Quantification of spine density (C) and spine head size (D). Quantifications are shown as frequency distributions. Spine density (segments): nControl = 26, nF-BARx = 25, nF-BARΔ49 = 34, nSRGAP2B = 25, nSRGAP2C = 25; Spine size (number of spines): nControl = 1693, nF-BARx = 2054, nF-BARΔ49 = 3267, nSRGAP2B = 2900, nSRGAP2C = 2448; one-way ANOVA with Holm-Sidak’s post hoc multiple comparison test; *p < 0.05, ****p < 0.0001 compared to control.

Figure 4
figure 4

In vivo analysis of spine density and spine size at the adult stage. (A) Representative images showing apical oblique dendrites of adult neurons co-expressing either Venus (shown) and tdTomato (Control), or mRFP-tagged F-BARx, F-BARΔ49, SRGAP2B, or SRGAP2C (not shown). Scale bar: 2 µm. (B,C) Quantification of spine density (B) and spine head size (C). Quantifications are shown as frequency distributions. Spine density (segments): nControl = 25, nF-BARx = 24, nF-BARΔ49 = 22, nSRGAP2B = 27, nSRGAP2C = 25; Spine size (number of spines): nControl = 1872, nF-BARx  = 2049, nF-BARΔ49 = 1843, nSRGAP2B = 2295, nSRGAP2C = 2630; one-way ANOVA with Holm-Sidak’s post hoc multiple comparison test; NS: p > 0.05, *p < 0.05, ***p < 0.001, ****p < 0.0001 compared to control.

We first analyzed how these proteins affected spine density. Compared to tdTomato controls, F-BARx slightly increased spine density at P21, while expression of SRGAP2C induced an increase in spine density at P21 that was significantly larger compared to either control or F-BARx expression (Fig. 3C). Interestingly, both F-BARΔ49 and SRGAP2B had intermediary phenotypes, with increased spine density compared to control and F-BARx, but significantly smaller compared to SRGAP2C.

When analyzing spine head size, which is linearly proportional to postsynaptic AMPA receptor content and therefore classically used as an index of synaptic maturation18,19,20, we found that for the control condition, spine head size at P21 was almost identical to that found at adult stages, showing that spines have reached maturity at P21 (Figs. 3D and 4C). In contrast, F-BARx expression led to reduced spine size compared to control at P21 (Fig. 3D). A similar reduction in spine size was observed for F-BARΔ49 and SRGAP2B, while a significantly larger reduction in spine size was observed for SRGAP2C.

We next analyzed whether changes in spine density and head size observed at P21 persisted into adulthood for all the constructs tested (Fig. 4), as previously reported for humanization of SRGAP2C expression or loss of SRGAP2A6,7. Interestingly, only SRGAP2C significantly increased spine density in adults (Fig. 4A,B). The increased spine density observed for F-BARΔ49 and SRGAP2B at P21 was no longer observed at adult stages. This persistent increase of synaptic density induced by SRGAP2C in adult layer 2/3 PNs cannot be due to differences in protein expression levels between constructs (quantified in Fig. S5). While spine head size was slightly smaller in adult PNs for all of the proteins tested when compared to control, the difference in size was exceedingly small (for SRGAP2C a reduction of 8.2%), suggesting that for all conditions dendritic spines did mature to sizes very similar to controls at P65.

Discussion

The first evolutionary step that led to the emergence of SRGAP2B and then SRGAP2C was an incomplete genomic duplication of the ancestral copy SRGAP2A generating a truncated protein consisting of a truncated extended F-BAR domain with a unique C-terminal region that lacks the last 49 amino acids and gained 7 unique amino acids (F-BARx)8,21. This led to the formation of a highly insoluble protein that is sensitive to aggregate formation15. We show that in cortical pyramidal neurons this truncated protein is targeted to the proteasome for degradation. Moreover, our data also show that SRGAP2C reduces SRGAP2A levels in a proteasome-dependent manner. Since SRGAP2C directly binds SRGAP2A6, our results demonstrate that SRGAP2C downregulates SRGAP2A function by binding to it and conferring to the heterodimers proteosome-sensitivity. SRGAP2A plays a critical role in regulating synaptic density and maturation6,7. Restricting the available amount of SRGAP2A by targeting it to the proteasome represents an effective mechanism through which SRGAP2C increases synaptic density and delays maturation. Interestingly, HSGDs have often resulted in the formation of truncated paralogs, and several of these truncated HSGDs have recently been implicated in regulating brain development2,3,4,8,22,23. This suggests that evolutionary changes affecting levels and patterns of expression not only occurred through modifications in regulators of gene expression at the genomic level24, but also at the level of protein interactions through the emergence of insoluble proteins.

Following the initial duplication event of SRGAP2A, the newly generated SRGAP2B paralog acquired multiple base pair mutations specifically targeting five arginine residues, leading to the emergence of SRGAP2C. Previous work has shown how conversion of these arginine residues leads to reduced membrane binding for SRGAP2C and increased stability of the SRGAP2A:SRGAP2C heterodimer15. Besides SRGAP2C, the human-specific paralog SRGAP2B is also expressed in the human brain2. One of the question raised by this finding is whether SRGAP2B and SRGAP2C are redundant in their ability to downregulate SRGAP2A function or whether SRGAP2C has a distinct function in regulating synaptic development. The latter seems more plausible since SRGAP2B displays more CNVs in the human population than SRGAP2C8, strongly suggesting that these two paralogs have been under different selective pressures since their evolutionary emergence in the human population. Our data show that SRGAP2C has a unique ability compared to SRGAP2B not only to delay synaptic maturation but also to increase synaptic density into adulthood. Considering that these proteins differ only in a small number of nonsynonymous base pair mutations, these results further underline the importance of the conversion of these five arginine residues during emergence of SRGAP2C.

Our results not only provide a molecular mechanism whereby SRGAP2B/C inhibits SRGAP2A function by targeting it for proteasome degradation but also, and for the first time, provide evidence that SRGAP2C (and its 5 specific arginine mutations) is more potent in its ability to inhibit synaptic maturation and maintain increased spine density throughout adulthood compared to SRGAP2B. This might explain why these two paralogs are under different selection pressures during recent human evolution and, more importantly, why emergence of SRGAP2C has played an important role during human brain evolution.

Methods

Animals

All animals were handled according to NIH guidelines and protocols approved by the institutional animal care and use committee (IACUC) at Columbia University, New York. Timed pregnant CD1 mice (Crl:CD1(ICR)) were obtained from Charles River and maintained in a 12 h light/dark cycle. Juveniles correspond to mice at P21. Adults were P65 or older.

Constructs

Plasmids for expression of SRGAP2A-RFP, SRGAP2C-RFP, SRGAP2B-RFP, F-BARΔ49 -RFP, 5R-RFP, F-BARx-RFP, and myristoylated-Venus (mVenus), were as previously described6. SRGAP2A-BFP was made by exchanging mVenus for Blue Fluorescent Protein (BFP).

Primary neuronal culture

Cortical neurons were ex utero electroporated (EUE) by injecting endotoxin-free DNA (500 ng/µl for cell filler and 750 ng/µl for each fusion protein construct) into the lateral ventricle of E15.5 mouse embryos and applying a current of 5 pulses of 20 V for 100 ms with 200 ms intervals using a square wave electroporator (ECM 830, BTX). Cortices were subsequently dissected in complete HBSS (Hank’s Buffered Salt Solution (HBSS) supplemented with 2.5 mM Hepes (pH 7.4), 30 mM D-glucose, 1 mM CaCl2, 1 mM MgSO4, and 4 mM NaHCO3). Dissociation was performed in papain (Worthington) supplemented with DNAse I (100 µg/µl) for 15 min at 37 °C. We then performed three washes and triturated manually in plating medium. Cells were plated on cover slips or glass bottoms dishes coated with poly-D-lysine and cultured for 18–21 days in vitro (DIV). During the first 5 days, cells were cultured in neurobasal medium supplemented with 2.5% fetal bovine serum (FBS), B27, and L-glutamine. Half of the medium was replaced every 5 days with the same plating medium without FBS.

Expression of proteins in cell lines and western blotting

Neuro-2a cell line was obtained from ATCC (CCL-131) and cultured according to the recommended protocol. Cells were transfected using Jet-Prime (Polyplus Transfection) according to the manufacturer protocol. After 48 h of transfection cells were scraped and lysed in lysis buffer (20 mM HEPES, 150 mM NaCl, 1 mM EDTA, 1 mM EGTA, 1% Triton-X100, cOmplete Protease Inhibitor Cocktail (Roche), Benzonase (EMD Millipore)) for 30 min at 4 °C. Samples were prepared in Laemmli buffer (Bio-Rad) with 10% 2-Mercaptoethanol and boiled at 95 °C for 5 min. Proteins were separated using SDS-PAGE and transferred to a polyvinylidene difluoride (PVDF) membrane (Immobilon-FL, EMD Millipore). Western blotting was performed using rabbit-anti-SRGAP2 N-terminal (raised against residues 193–205, 1:10006), and goat-anti-rabbit IgG conjugated to IRDye 600RD (1:20000, Li-Cor). Immunoblots were visualized on an Odyssey CLx Imaging System (Li-Cor).

Fixed and live-cell imaging of dissociated neurons

For quantification of SRGA2A levels in soma and dendritic spines, neurons were fixed with 4% paraformaldehyde (Electron Microscopy Sciences) in PBS for 10 min at room temperature. Immunocytochemistry was performed for conditions where we used mTagBFP-HA-tagged SRGAP2A. Cells were incubated for 30 min in 0.2% Triton-X100 in PBS containing 5% goat serum to permeabilize and block nonspecific staining. Cells were incubated overnight at 4 °C with mouse anti-HA primary antibody (1:1000, Anti-HA.11, Biolegend), washed 3x in PBS and incubated with Alexa-conjugated anti-mouse secondary antibody (1:500, Alexa Fluor 647, Invitrogen) for 1 h at room temperature. Coverslips were mounted on slides with Fluoromount-G aqueous mounting medium (ThermoFisher Scientific). Images were acquired on a Nikon A1 confocal microscope using a 60x (soma) or 100x (spines) objective. Exposure time and laser settings were optimized for bright cells to avoid pixel saturation and microscope settings were kept identical for each imaging session.

Live cell imaging was performed using a 40x objective on a Nikon A1 confocal microscope together with a stage incubation chamber (Tokai Hit) to maintain 37 °C/5% CO2 culture conditions. Before imaging, medium was exchanged to recording medium consisting of pre-warmed cHBSS (Hank’s Buffered Salt Solution (HBSS) supplemented with 2.5 mM Hepes (pH 7.4), 30 mM D-glucose, 1 mM CaCl2, 1 mM MgSO4, and 4 mM NaHCO3) with either 10 µM MG-132 or 200 µg/mL Cycloheximide (Sigma-Aldrich). Z-stacks that spanned the entire soma for each neuron were collected every 30 min for the indicated number of hours.

Analysis of images was performed using Nikon NIS-Elements (Nikon Corporation, Melville, NY). Background correction was performed by drawing an ROI at an area without fluorescent labeling and images were subsequently aligned to correct for sample drift. A z-stack maximum projection was generated for each neuron and using the filler signal ROIs were drawn over the soma, excluding the nucleus, or dendritic spines to measure fluorescent intensity.

In utero cortical electroporation and in vivo analysis

In utero cortical electroporation was performed on timed pregnant CD1 females (Crl:CD1(ICR)). Endotoxin-free DNA containing 1 µg/µl of each plasmid was injected into the ventricles of E15.5 embryos. Electroporation was performed using a square wave electroporator (ECM 830, BTX) by applying 5 pulses of 42 V for 50 ms with 500 ms intervals. When animals reached the indicated age they were anaesthetized with isoflurane and intracardiac perfusion with 4% paraformaldehyde (Electron Microscopy Sciences) in PBS was performed. Brains were subsequently isolated and incubated overnight in 4% paraformaldehyde/PBS at 4 °C. Coronal brain sections were prepared by slicing brains at 100 µm using a vibrating microtome (Leica VT1200S). Sections were mounted on glass slides in Fluoromount-G aqueous mounting medium (ThermoFisher Scientific).

Imaging of dendritic spines was performed on optically isolated neurons using a Nikon A1 confocal microscope and 100x objective. Using Nikon NIS-Elements software maximum intensity projections were generated of oblique dendrites originating from the main apical trunk. Spine density and head size was quantified by using the Venus filler signal to draw ROIs around the spine head.

Statistical analysis

A minimum of three independent experiments were performed unless otherwise stated. For in vivo spine analysis, data was obtained from at least six animals from a minimum of two independent litters. Statistical analysis was performed using Prism (Graphpad Software). Normality was checked using D’Agostino-Pearson omnibus normality test. A non-parametric test (Mann-Whitney) was used when distribution deviated significantly from normality. A test was considered significant when p < 0.05.