Abstract
Pyramidal neurons in layers 2/3 and 5 of primary somatosensory cortex (S1) exhibit somewhat modest synaptic plasticity after whisker input deprivation. Whether neurons involved at earlier steps of sensory processing show more or less plasticity has not yet been examined. Here, we used longitudinal in vivo two-photon microscopy to investigate dendritic spine dynamics in apical tufts of GFP-expressing layer 4 (L4) pyramidal neurons of the vibrissal (barrel) S1 after unilateral whisker trimming. First, we characterize the molecular, anatomical, and electrophysiological properties of identified L4 neurons in Ebf2-Cre transgenic mice. Next, we show that input deprivation results in a substantial (∼50%) increase in the rate of dendritic spine loss, acutely (4–8 d) after whisker trimming. This robust synaptic plasticity in L4 suggests that primary thalamic recipient pyramidal neurons in S1 may be particularly sensitive to changes in sensory experience. Ebf2-Cre mice thus provide a useful tool for future assessment of initial steps of sensory processing in S1.
Introduction
As recipients of most synaptic inputs in the neocortex, dendritic spines are central to neuronal function. Using time-lapse in vivo two-photon microscopy of fluorescently labeled neurons, it has been possible to study dynamic aspects of synapse formation and elimination (Holtmaat and Svoboda, 2009), which correlate with synaptic adaptations taking place during normal brain development (Holtmaat et al., 2005; Zuo et al., 2005a; Cruz-Martín et al., 2010), in response to changes in sensory experience (Holtmaat and Svoboda, 2009) or as a result of learning (Xu et al., 2009; Yang et al., 2014). The importance of spine dynamics has been exemplified in the context of experience-dependent plasticity in the vibrissal primary somatosensory cortex (S1), where sensory input deprivation can alter the rate of spine turnover of layer (L)5B pyramidal neurons or the likelihood that they are stabilized (Zuo et al., 2005b; Holtmaat et al., 2006; Wilbrecht et al., 2010; Schubert et al., 2013).
The cortical representation of the rodent vibrissae is somatotopically organized into discrete columnar units in S1 known as “barrels” (Feldmeyer et al., 2013). Complete whisker trimming in adult mice leads to reduced spine elimination 1 month later, but not in the first 2 weeks after trimming (Zuo et al., 2005b). In contrast, “chessboard” trimming (which promotes competition between neighboring barrels) results in a greater probability of stabilization of new persistent (>8 d) spines in L5B neurons at the border between spared and deprived barrels (Holtmaat et al., 2006; Wilbrecht et al., 2010). Intriguingly, the more drastic paradigm of whisker follicle ablation does not change either the density or the dynamics of spines from L5B neurons, but leads to an increase in the density of new spines of L2/3 pyramidal neurons (Schubert et al., 2013).
A limitation of these studies has been the focus on L5B (occasionally L2/3) pyramidal neurons because those are the cells labeled in Thy1 GFP-M and YFP-H mice (Feng et al., 2000). In contrast, spines of L4 neurons, which are the primary recipients of thalamic sensory inputs, have not yet been studied in vivo. Here, we tested the hypothesis that L4 pyramidal neuron spines exhibit experience-dependent plasticity, because of the privileged role of L4 in the cortical sensory processing hierarchy (Feldmeyer et al., 2013), albeit direct inputs from thalamus to L5 neurons have also been described by Constantinople and Bruno (2013). Therefore, we examined the effect of input deprivation on apical dendritic spines of L4 pyramidal neurons (readily accessible for in vivo imaging), which are thought to both receive direct thalamic sensory information and integrate information from multiple inputs (Staiger et al., 2004; Schoonover et al., 2014). Using longitudinal two-photon microscopy in the barrel cortex of Ebf2-Cre mice, we find that sensory deprivation by contralateral whisker trimming leads to a rapid loss of spines (within 4–8 d) in L4 pyramids.
Materials and Methods
Animals and rAAV constructs.
We used Ebf2-Cre BAC-transgenic mice (MGI:4421668) of either sex expressing Cre recombinase under the Ebf2 promoter and genotyped for Cre (Chiara et al., 2012). All the procedures described in this study were approved by the University of California Chancellor's Animal Research Committee and the UCLA EH&S Biosafety Division. rAAVs were purchased from the University of Pennsylvania Vector Core: rAAV2/1.CAG.FLEX.EGFP.WPRE.bGH (“rAAV-EGFP”; Allen Institute 854) and rAAV2/1.EF1a.DIO.hChR2(H134R)-EYFP.WPRE.hGH (“rAAV-ChR2”; Addgene 20298).
Whisker trimming.
Unilateral trimming of the contralateral vibrissae was performed every 1–2 d for 16 d using a whisker pocket trimmer (Wahl) without anesthesia.
In utero rAAV injections.
In utero injections were described previously (Cruz-Martín et al., 2010). Approximately 0.3 μl of rAAV-solution (titer = e12 − e13 with 0.1% Fast Green) was pressure-injected into the right lateral cerebral ventricle at E15.
Histology.
Mouse perfusion and tissue preparation was performed as described previously (Chiara et al., 2012). Briefly, sections (30–40 μm) were blocked in PBS containing 5% Goat serum (Thermo Fisher Scientific PCN5000) and 0.1% Triton X-100 (Sigma-Aldrich), and then incubated with primary antibodies overnight at RT at the following dilutions: chicken anti-GFP (1:500, Millipore Bioscience Research Reagents AB16901), guinea pig anti-VGlut2 (1:2000, Millipore Bioscience Research Reagents, AB2251), rabbit anti-Cux1 (1:500, Santa-Cruz Biotechnology, sc-13024), rat anti-Ctip2 (1:750, Abcam AB18465), and rabbit anti-GABA (1:500, Sigma-Aldrich A2052). For secondary antibodies, we used donkey anti-chicken-FITC (Jackson Laboratories), AlexaFluor 568 donkey anti-rabbit, AlexaFluor 568 goat anti-guinea pig, and Cy5-goat anti-rat (all from Life Technologies). Vectashield (Vector Laboratories) was used for mounting.
Morphological analysis and cell quantification.
Neuron reconstruction and morphological analysis was performed in Neurolucida (MicroBrightField; n = 48/11 neurons/mice). For unbiased cluster analysis, we used k means in MATLAB using three parameters: maximum branch order, total length of apical dendrite, and depth at first bifurcation. Cell density quantifications (Fig. 1) were done in ImageJ (n = 5 ROIs 0.2 × 0.8 mm, 5 sections from 3 mice) comprising the entire cortical depth at S1; normalized values were averaged and plotted in GraphPad Prism and expressed as mean ± SEM.
Electrophysiology and optogenetics.
Acute brain slices (350 μm thick) from adult Ebf2-Cre;rAAV-ChR2 mice were obtained using vibratome (VT1000, Leica), placed on a chamber maintained at 32°–35°C and submerged in ACSF containing the following (in mm): 119 NaCl, 2.5 KCl, 1 MgSO4, 1.25 NaH2PO4, 2 CaCl2, 26 NaHCO3, and 25 glucose, perfused at a rate of 2–4 ml/min, and bubbled with 95% O2 and 5% CO2. L4 and L2/3 neurons were identified using EYFP fluorescence and differential interference contrast (DIC) optics, respectively, using a 40×/0.9 NA water-immersion objective (Olympus). Glass microelectrodes (4–6 MΩ) for patch-clamp recordings were filled with internal solution containing the following (in mm): 130 K-gluconate, 10 KCl, 10 HEPES, 10 Na phosphocreatine, 4 Mg-ATP, 0.4 GTP, 2 NaCl, and 0.02 AlexaFluor 594. Recordings were performed using whole-cell technique in current-clamp configuration with a patch-clamp amplifier (Multiclamp, Molecular Devices). Input resistance (Rin) was calculated at −70 pA; Vm and Rin were expressed as median ± SD (n = 3 neurons, 3 mice). For optogenetics, short pulses (3 ms) of blue light (473 nm) were delivered through the objective using a custom-made LED stimulation system.
Cranial window surgery.
Chronic glass-covered cranial windows were implanted in adult mice as described previously (Mostany and Portera-Cailliau, 2008; Holtmaat et al., 2009).
IOS.
Intrinsic optical signal (IOS) imaging was performed to map the barrel field using contralateral stimulation of a bundle of whiskers and custom-written MATLAB routines, as described previously (Johnston et al., 2013).
High-resolution in vivo two-photon imaging of dendritic structure.
Imaging was done under light isoflurane anesthesia (1–1.5%) with a custom-built two-photon microscope, using a Ti:Sapph laser (910 nm, Chameleon Ultra II, Coherent), as previously described (Mostany et al., 2013), using ScanImage software (Pologruto et al., 2003). A preliminary imaging session at low-magnification was performed to identify potential candidates for dendritic imaging. High-magnification images (512 × 512 pixels, 0.152 μm/pixel, 1.5 μm z-steps) were obtained for analysis of dendritic spines in apical tufts of L4 pyramidal neurons. For display purposes, we used “best” projections (xyz) of dendritic stacks, as described previously (Mostany et al., 2013).
Spine analysis.
Data on the density and dynamics of spines was obtained using the Spine Analysis program included in ScanImage software. We analyzed dendritic spines in n = 35 dendritic segments from 29 neurons in nine mice (4 M/5 F; age = 5.5 ± 0.5 months). Total number of spines n = 511 over 0.95 mm of total dendritic length. Dynamics of spines are expressed as a fraction of the total spine number. To avoid unnecessary statistical contrasts, imaging sessions were grouped (in days) as follows: pre/basal (−8, −4, and 0), early (+4, +8), and late (+12, +16) after treatment (beginning of whisker trimming), and we compared post-treatment groups with baseline. We performed a linear mixed model with a fixed treatment effect and a random effect with “dendrite” nested within “mouse,” which takes into account that observations from the same mouse are correlated. Multiple comparisons were controlled using the false discovery rate at 5%. Survival fraction was calculated as described previously (Mostany et al., 2013), with groups compared using one-way-ANOVA with Tukey's post hoc test (95% confidence interval).
Results
L4 excitatory neurons are labeled in Ebf2-Cre mice
We identified expression in L4 neurons of adult S1 cortex when we injected E15 Ebf2-Cre mouse embryos (Chiara et al., 2012) with a conditional rAAV-flex-EGFP (rAAV-EGFP) construct (see Materials and Methods; Fig. 1a). A similar pattern of expression in L4 was seen when we crossed Ebf2-Cre mice to fluorescent reporter R26-stop-EYFP mice to generate Ebf2-Cre;R26R-YFP offspring (data not shown), though isolated expression was also found in infragranular layers (especially L6), in the subplate (SP), and in surviving Cajal-Retzius neurons (Chiara et al., 2012). Compared with Ebf2-Cre;R26R-YFP, rAAV-EGFP injections in Ebf2-Cre mice produced minimal expression in L6/SP, presumably due to selective rAAV tropism (Fig. 1a,b). We found that the enrichment of Ebf2-Cre expression (referred to as Ebf2+) in L4 begins early in postnatal cortical development [postnatal day (P)10–P11] and does not colocalize with GABA (data not shown).
We further identified Ebf2+ cells as excitatory neurons of L4 with additional experiments. First, we injected Ebf2-Cre embryos at E15 with a conditional rAAV-DIO-ChR2-EYFP (rAAV-ChR2) vector, which revealed barrel-like clusters of expression in S1. Of note, although the vast majority of Ebf2+ neurons exhibited a morphology typical of L4 spiny stellate neurons, a subset of Ebf2+ neurons were identified as L4 pyramidal (Ebf2+p) neurons, as evidenced by the presence of an apical dendrite extending toward the pial surface (Fig. 1c′,d′). Second, immunostaining against Cux1 and Ctip2 (Fig. 1c,d), markers of excitatory neurons of L2–4 and L5–6, respectively (Li et al., 2013; Pouchelon et al., 2014) showed that Ebf2 expression completely overlapped with the lower border of Cux1 expression (Fig. 1c,c′) but was excluded from the territory of Ctip2 expression (Fig. 1d). Importantly, Ebf2+p neurons clearly expressed Cux1 at the soma but were devoid of Ctip2 expression (Fig. 1c′,d′). Third, Vglut2 labeling of thalamocortical axons at P11 (Li et al., 2013; Fig. 1e,f) and cytochrome oxidase staining of adult mice (data not shown) both showed that the vast majority of Ebf2+ neurons reside within barrels.
We next performed whole-cell patch-clamp recordings of visually identified Ebf2+p neurons in acute brain slices from Ebf2-Cre;rAAV-EGFP mice. Ebf2+p neurons in L4 had a resting membrane potential (Vm) of 65 ± 3.7 mV and input resistance (Rin) of 144.3 ± 16.5. A depolarizing current injection pulse of 100 pA triggered a train of action potentials and spike adaptation characteristic of regular spiking pyramidal neurons (Fig. 1g). Optogenetic stimulation of Ebf2-Cre;rAAV-ChR2 expressing L4 neurons (including pyramidal cells), generated reliable action potential responses to brief (3 ms) pulses of blue light at rates of 10 Hz (Fig. 1g). Optogenetic stimulation of Ebf2+ neurons was sufficient to synaptically trigger EPSP responses in individually recorded L2/3 pyramidal neurons in a DNQX-dependent manner (Fig. 1h).
Anatomical features of Ebf2+ L4 pyramidal neurons based on the morphology of their apical dendritic tufts
Although excitatory neurons of the mouse somatosensory L4 are mainly represented by spiny stellate neurons, pyramidal and star pyramidal neurons have also been described (Lübke et al., 2000; Staiger et al., 2004; Schoonover et al., 2014). The sparse labeling of Ebf2+p obtained with our approach allowed us to perform in vivo two-photon imaging of their apical dendritic tufts through a cranial window implanted over barrel cortex. We used Neurolucida to trace dendritic arbors in the image stacks and analyzed the morphological characteristics of 48 Ebf2+p neurons (Fig. 2a–d; see Materials and Methods). In 19 of these neurons (40%) we successfully traced the primary apical dendrite all the way back to its cell body at an average depth of 399.8 ± 15.5 μm (which is similar to the average depth of spiny stellate neurons; 377.7 ± 11.9 μm, n = 54), corresponding to the expected depth for L4 in adult mice (310–520 μm below the pia (Jia et al., 2014)). The average total length of the entire dendritic tuft (measured from the first branching point of the main apical branch) was 730.6 ± 84.0 μm; the average depth below the pial surface of the first branching point was 92.0 ± 7.4 μm; the average number of branches per neuron was 9.75 ± 0.89; and the average maximum branch order was 4.2 ± 0.2 (Fig. 2e–h). We used a k means cluster analysis to sort Ebf2+p neurons into simple and complex categories based on morphological characteristics (Fig. 2c); most (∼94%) Ebf2+p neurons were “simple” (possessing simple apical dendritic arbors), while a minority (3/48) had more complex apical dendritic arbors that bifurcated deeper below the pia and had cell bodies at depths >500 μm below the pia (data not shown).
Experience-dependent plasticity of apical dendritic spines in L4 Ebf2+ pyramidal neurons in barrel cortex
To investigate whether L4 Ebf2+p neurons exhibit dynamic changes in response to whisker input deprivation, we longitudinally imaged apical dendritic branches parallel to the pial surface in adult Ebf2-Cre;rAAV-EYFP mice. We confirmed that neurons we imaged were in barrel cortex with IOS imaging (Fig. 3a; see Materials and Methods). We chose a standard protocol for imaging spines every 4 d (Holtmaat et al., 2006; Mostany et al., 2010, 2013; Wilbrecht et al., 2010; Johnston et al., 2013; Kuhlman et al., 2014), which accounts well for the existence of populations of spines with different lifetimes (Mostany et al., 2013). Spine dynamics were analyzed before and after whisker trimming at two time-points: early-post (4–8 d) and late-post (12–16 d; Fig. 3a,b).
After sensory input deprivation (unilateral trimming of all whiskers on the contralateral snout every 2 d), the density of spines did not change significantly (0.36 ± 0.03 spines/μm at “pre” conditions vs 0.37 ± 0.03 spines/μm at early-post, and 0.036 ± 0.03 spines/μm at late-post; p = 0.99; n = 35 dendritic segments, 29 neurons, 9 mice; Fig. 3c). Similarly, the fraction of spines gained after each imaging session was not significantly altered (p = 0.50; 0.13 ± 0.01 pre vs 0.11 ± 0.01 early, and 0.11 ± 0.02 spines/μm late; Fig. 3d). In contrast, the fraction of lost spines exhibited a significant ∼50% increase shortly after input deprivation (0.10 ± 0.02 at basal pre vs 0.15 ± 0.02 at 4–8 d at early-post; p = 0.029 pre vs early-post; Fig. 3e). The fraction of lost spines then returned to baseline levels 12–16 d after the onset of trimming (0.10 ± 0.01 at late-post; p = 0.92, pre vs late-post; Fig. 3e). We also analyzed the survival fraction of spines present at −8, 0, and +8 d that remains at two subsequent imaging sessions, and found no significant differences between groups (pre-trimming baseline: 0.86 ± 0.02; early-post trimming: 0.80 ± 0.02; late-post trimming: 0.86 ± 0.02; p = 0.11, ANOVA; Fig. 3g). Despite the increase in spine loss, overall spine turnover rate (defined as the fraction of spines that appear and disappear over time) did not change significantly after whisker trimming (0.12 ± 0.01 at pre vs 0.13 ± 0.01 at early-post and 0.10 ± 0.01 at late-post; p = 0.16; Fig. 3f).
Discussion
Using chronic in vivo two-photon imaging at the barrel cortex, we investigated adult experience-dependent plasticity in dendritic spines of L4 pyramidal neurons, known recipients of thalamic innervation and ideal candidates for input integration due to the columnar and translaminar nature of their dendritic projections (Lübke et al., 2000; Staiger et al., 2004; Schoonover et al., 2014). We identified robust synaptic plasticity with a dramatic (50%) increase in spine elimination immediately after whisker trimming, followed by a return to normal spine dynamics at 2 weeks of continuous whisker deprivation. To our knowledge, this is the first demonstration of experience-dependent changes in spine dynamics in L4 neurons in vivo. Thus, complete unilateral input deprivation leads to a loss of pre-existing spines in L4 pyramidal neurons in the contralateral cortex, presumably in response to loss of afferent sensory input.
In the absence of compensatory changes in the fraction of spines gained (Fig. 3d), we might have expected to find a slight decrease in spine density (on the order of 0.02/μm) but this was within our margin of error (±0.03/μm) when estimating spine density (Fig. 3c). Notably, previous in vivo imaging studies of experience-dependent spine plasticity in L5B and L2/3 pyramidal neurons after input deprivation (Fig. 4) also found stable spine densities (Trachtenberg et al., 2002; Zuo et al., 2005b; Holtmaat et al., 2006; Wilbrecht et al., 2010; Schubert et al., 2013). The increased rate of spine elimination was transient and could not be detected after 2 weeks of continuous whisker trimming, which fits with prior observations in L5B, where whisker input deprivation led to delayed spine elimination >2 months (Zuo et al., 2005b), but not at 2 weeks after deprivation. Perhaps not surprisingly, spine plasticity in L4 neurons precedes that in L5B and correlates with the reduction of thalamocortical input density in L4 upon unilateral trimming (>3 d; Oberlaender et al., 2012). Of course, differences in cell types (L4 vs L5B) and/or trimming paradigms (complete unilateral vs chessboard) may also account for differences in results between different studies, as previously argued (Chen et al., 2014).
This study was possible because we identified a BAC transgenic mouse line (Ebf2-Cre) that drives strong Cre expression in L4 excitatory neurons (both spiny stellates and pyramidal neurons). A subset of these L4 neurons demonstrated a pyramidal morphology (Fig. 2). Cluster analysis of quantified morphological parameters of apical dendrites from L4 pyramidal neuron revealed that Ebf2+ L4 pyramidal neurons are typically simple-tufted and comparable to what others have referred to as “class 3” narrow L5 pyramidal neurons in primary visual cortex (Tsiola et al., 2003). A small fraction (6%) of Ebf2+p cells was similar to L5B complex pyramidal neurons labeled in Thy1 mice (Holtmaat et al., 2006; Mostany et al., 2013).
Sensory processing in the barrel cortex involves the orchestrated activity of different cell types. Both L4 and infragranular layers receive parallel independent thalamic inputs (Constantinople and Bruno, 2013) and L4 neurons are readily activated in vivo by such inputs (Jia et al., 2014; Schoonover et al., 2014). Moreover, L4 (but not L5) neurons are affected by impairments in thalamocortical drive (Li et al., 2013). Whether the experience-dependent spine plasticity that we observed in Ebf2+p neurons is a direct consequence of thalamocortical input loss [presumably from posteromedial nucleus (POm); Meyer et al., 2010] or a correlate of more complex information processing in the cortex, will require further studies. What is clear is that, despite accounting for a smaller fraction of the total input to L4, POm synapses (paralemniscal pathway; Fig. 4) may be the main drivers of activity onto apical dendrites of L4 neurons (Feldmeyer et al., 2013; Pouchelon et al., 2014) and thus have a substantial impact on the circuit.
Footnotes
This study was supported by UCLA's CTSI Grant Number UL1TR000124 (NIH/NCATS), by Grants 1R01MH083785 and 1R21MH100614 from the National Institute for Mental Health (C.P.-C.), a NARSAD Young Investigator Grant from the Brain and Behavior Research Foundation (A.M.), an Undergraduate Research Fellowship Program funded by The Milton Gottlieb Scholarship (S.K.), and the Italian Telethon Foundation (G.G.C.). We thank Drs Carlos Cepeda and Joyce Wondolowsky for helpful discussions regarding the electrophysiology, Dr Karel Svoboda and Tim O'Connor for the spine analysis software, Dr Daniel Fiole for assistance with MATLAB, Dr Bennett Novitch, Dr Caroline Pearson, Dr Barbara Rust, and Aaron Lulla for assistance with immunohistochemistry, Drs Matthew Shtrahman and Raul Serrano for help with optogenetics hardware, and Sitaram Vangala, Tristan Grogan (UCLA-CTSI), Jeffrey Gornbein (UCLA-Biomathematics), and Dr. Anubhuthi Goel for help with statistics.
The authors declare no competing financial interests.
- Correspondence should be addressed to either Dr Amaya Miquelajauregui or Dr Carlos Portera-Cailliau, Departments of Neurology and Neurobiology, David Geffen School of Medicine at UCLA, Reed Neurological Research Center, Room A-145, 710 Westwood Plaza, Los Angeles, CA 90095. amaya.miq{at}gmail.com or cpcailliau{at}mednet.ucla.edu