Neocortical circuitry can alter throughout life with experience. However, the contributions of changes in synaptic strength and modifications in neuronal wiring to experience-dependent plasticity in mature animals remain unclear. We trimmed whiskers of rats and made electrophysiological recordings after whisker cortical maps have developed. Measurements of miniature EPSPs suggested that synaptic inputs to layer 2/3 pyramidal neurons were altered at the junction of deprived and spared cortex in primary somatosensory cortex. Whole-cell recordings were made from pairs of synaptically connected pyramidal neurons to investigate possible changes in local excitatory connections between layer 2/3 pyramidal neurons. The neurons were filled with fluorescent dyes during recording and reconstructed in three dimensions using confocal microscopy and image deconvolution to identify putative synapses. We show that sensory deprivation induces a striking reduction in connectivity between layer 2/3 pyramidal neurons in deprived cortex without large-scale, compensatory increases in the strength of remaining local excitatory connections. A markedly different situation occurs in spared cortex. Connection strength is potentiated, but local excitatory connectivity and synapse number per connection are unchanged. Our data suggest that alterations in local excitatory circuitry enhance the expansion of spared representations into deprived cortex. Moreover, our findings offer one explanation for how the responses of spared and deprived cortex to sensory deprivation can be dissociated in developed animals.
The neocortex is not hardwired but reorganizes with persistent changes in the sensory periphery during development (Hensch, 2004; Feldman and Brecht, 2005) and adulthood (Fox and Wong, 2005). Cortical map reorganization induced by brief periods of sensory deprivation is a popular experimental paradigm with which to study experience-dependent plasticity. Plasticity is commonly greatest in layer 2/3 (L2/3) at the junction of cortex that has lost its sensory input (deprived) and cortex that has retained its sensory input (spared) (Fox, 1992; Gilbert and Wiesel, 1992; Diamond et al., 1994; Glazewski and Fox, 1996). Distinct mechanisms with different temporal profiles drive changes in sensory responses to deprived and spared inputs (Fox, 2002). In somatosensory cortex, responses to deprived inputs depress rapidly (over a few days) after whisker deprivation during adolescence, but this effect becomes less prominent when deprivation is started in adult animals (Glazewski and Fox, 1996; Wallace and Fox, 1999). In contrast, responses to spared sensory inputs are enhanced. This enhancement develops more slowly than the depression of responses in deprived cortex of adolescent animals but can be elicited in adulthood (Glazewski et al., 1996).
The relative contributions of changes in synaptic strength and modifications in neuronal wiring to either enhancement or depression of sensory responses in neocortex have not been fully elucidated (Chklovskii et al., 2004). Structural changes in axons (Darian-Smith and Gilbert, 1994; Kossut and Juliano, 1999) and dendritic spine turnover (Trachtenberg et al., 2002; Zuo et al., 2005b) have been reported after sensory deprivation, but the new synaptic partners acquired after axonal growth or new spine formation are not always identified.
In vitro studies of the effects of experience-dependent plasticity on synaptic strength in neocortex have commonly focused on vertical intracortical pathways or long-distance (hundreds of micrometers to millimeters) horizontal connections that cross cortical columns (Rioult-Pedotti et al., 1998; Finnerty et al., 1999; Allen et al., 2003; Bender et al., 2006; Clem and Barth, 2006). Changes in the strength of short-range (less than a few hundred micrometers) horizontal connections within cortical columns have attracted less attention. Pyramidal neurons in L2/3 and L5 do not connect randomly with neighboring pyramidal neurons but form local excitatory networks (Song et al., 2005; Yoshimura et al., 2005). These networks may be an important substrate for sensory processing that occurs on a subcolumnar basis (Brumberg et al., 1999) [e.g., angular tuning responses of whiskers (Andermann and Moore, 2006)]. The output from local networks in L2/3 provides much of the excitatory drive for the lateral, intracortical spread of activity in vivo (Petersen et al., 2003). Hence, local excitatory connections are in a position to play a key role in regulating the size of cortical maps. Here we report that alterations in sensory input, which start after cortical maps have developed, induce a marked reduction in local connectivity between L2/3 pyramidal neurons in deprived cortex of developed animals without large-scale modifications in connection strength. In contrast, connectivity in local excitatory networks within spared cortex is unchanged, but connection strength is potentiated.
Materials and Methods
All procedures were performed in accordance with the United Kingdom Animals (Scientific Procedures) Act 1986. We trimmed bilaterally the dorsal three (A–C) or the ventral two (D–E) rows of rats' whiskers flush with the skin every day. Whisker trimming or sham-trimming started at postnatal day 19 (P19), which is after the L2/3 whisker map has developed (Stern et al., 2001), and continued for 11–41 d.
Brain slice preparation and recording.
On the day of electrophysiological recording, all whiskers were trimmed acutely by an assistant to blind experimenters to the rat's sensory history. The rat's brain was cooled rapidly by transcardial perfusion with ice-cold dissection buffer comprising the following (in mm): 108 choline-Cl, 3 KCl, 26 NaHCO3, 1.25 NaHPO4, 25 d-glucose, 3 Na pyruvate, 2 CaCl2, and 1 MgCl2 (bubbled with 95%O2/5%CO2). Brain slices 300–400 μm thick were cut at 65° to the midline (Fig. 1A), incubated in artificial CSF containing 120 mm NaCl, 3 mm KCl, 23 mm NaHCO3, 1.25 mm NaHPO4, 10 mm d-glucose, 2 mm CaCl2, and 1 mm MgCl2 bubbled with 95%O2/5%CO2 at 35°C for 45 min and then maintained at room temperature until transferred to the recording bath, which was heated to 36–37°C. Bright-field illumination was used to identify brain slices containing whisker barrels (Fig. 1A). Neurons in layer 2/3 of the C or D barrel columns were visualized using infrared differential interference contrast. Recording pipettes (4–7 MΩ) contained the following: 130 mm KMeSO4, 8 mm NaCl, 2 mm KH2PO4, 2 mm d-glucose, 10 mm HEPES, 1 mm Alexa Fluor 488 (AF488) or 1 mm Alexa Fluor 568 (AF568) (Invitrogen, Carlsbad, CA), and 3 mg/ml biocytin.
The loose-seal current injection protocol was designed to minimize the risk of missing small-amplitude or facilitating connections. It comprised four trials of five pulses (5 ms, 1–10 nA, 20 Hz). Trials were repeated every 5 s. If there was a definite or possible response in the “postsynaptic” neuron, then the “presynaptic” neuron that had the loose seal made on it was repatched with a second pipette containing the internal solution described above and whole-cell recording was established. A pair of neurons was termed “connected” when both cells were recorded in whole-cell mode and action potentials in one neuron evoked EPSPs in the other neuron (see Figs. 1B, 5D). Connections were denoted “reciprocal” when the connection was bidirectional. Statistical analysis of connectivity only uses connections confirmed with whole-cell recordings. It does not incorporate the numbers of reciprocal connections. Connectivity data from initial nonblind experiments were pooled with data from recordings performed blind to the sensory history of the rats because there was no significant difference between the results for the two groups (data not shown).
The presynaptic neuron was stimulated to fire a train of eight action potentials by repetitive injection of brief (2–3 ms), depolarizing current pulses and the postsynaptic potential was recorded (Multiclamp 700A; Molecular Devices, Sunnyvale, CA). The amplitude of the depolarizing current pulses was adjusted to ensure that only one action potential was fired by the presynaptic neuron per current pulse. Action potentials were filtered at 10 kHz and postsynaptic responses at 3 kHz. All signals were digitized at 20 kHz. We fitted a single-exponential curve to a 3–10 ms period before EPSP onset and a second function, a + b × t × e−c × t + d × t × e−g × t, to the EPSP over a 3–20 ms window encompassing the EPSP peak (Fig. 1B). EPSP amplitude was estimated from the function generated by subtracting the extrapolated baseline curve from the EPSP curve. Steady-state amplitude was defined as the mean of the last four responses in the train. Data acquisition and analysis were performed with custom-written programs in Labview (National Instruments, Austin, TX) and Matlab (MathWorks, Natick, MA). Neuronal location was confirmed after recording by overlaying 5× images of the slice and the fluorescent dye-filled neurons. Septal width was estimated by measuring the distance between the D and E barrels (maximum of 120 μm), which is the largest interbarrel gap with this slice orientation (Fig. 1A). We excluded paired recordings if either neuron in the pair was <60 μm from the center of the supraseptal region between C and D barrel columns because L2/3 pyramidal neurons above L4 septa have different circuitry (Kim and Ebner, 1999; Brecht et al., 2003) and experience-dependent plasticity has different effects on circuitry in barrel columns and the septal projection to L2/3 during development (Shepherd et al., 2003).
Whole-cell recordings of miniature EPSPs (mEPSPs) were made at resting membrane potential with 1 μm tetrodotoxin and 100 μm picrotoxin (Tocris Bioscience, Bristol, UK) in the recording bath. One neuron per slice was studied for the longer-duration deprivation protocol (22–34 d). Neurons in the C, D, and E barrel columns of each slice were recorded for the shorter-duration deprivation protocol (13–20 d). The data acquisition parameters were the same as used for unitary EPSPs. We used shape criteria derived from measurements of evoked unitary EPSPs (see above, Paired recordings) (Fig. 1B) to identify mEPSPs in our traces. The criteria required a depolarizing potential with amplitude >62.5 μV (2.5 times the root mean square of baseline noise) (Clements and Bekkers, 1997), a time-to-peak <5 ms, and a rate of decay after peak that was between 2% of peak mEPSP amplitude per ms and (50/time to peak)% of peak mEPSP amplitude per millisecond. EPSPs that overlapped (collided) were considered to be separate when the second EPSP started ≥1 ms after the peak of the first EPSP. EPSP amplitudes were measured using custom-written software (Labview) under manual control. A minimum of 50 mEPSPs were measured per cell.
Slices were fixed in 4% paraformaldehyde for 30 min, mounted with Gel/Mount (Biomeda, Foster City, CA) and imaged on a Zeiss (Oberkochen, Germany) LSM510 META confocal microscope with a 63×, 1.2 numerical aperture C-Apochromat water immersion objective. Pinhole diameters (maximum, 1 Airy unit) were adjusted to give identical optical slice thicknesses for AF488 and AF568 tracks. Medium resolution z-stacks (≤120 μm) were acquired initially. Regions in which synapses were suspected were then reimaged at high resolution (voxel, 0.07 × 0.07 × 0.2 μm). Datasets were deconvolved with AutoDeblur (AutoQuant, Troy, NY) and imaged with Imaris (Bitplane, Zurich, Switzerland) to visualize synapses. The optical performance of the confocal microscope (scan head alignment, correction for chromatic aberrations, z-axis drift) combined with deconvolution was tested with 0.1 and 4 μm Tetraspeck fluorescent microspheres (Invitrogen) mounted similarly to the brain slices. The peak fluorescence signal from AF488 and AF568 fluorophores was within 0.1 μm in the image plane and within the same optical section. Criteria for synapses were based on previously reported electron microscopic reconstructions of excitatory synapses onto L2/3 pyramidal neurons (Peters and Kaiserman-Abramof, 1969) convolved with the point spread function of emitted light. Putative synapses were defined to occur when there was >0.40 μm overlap of a presynaptic varicosity with a postsynaptic spine or dendritic shaft in the x–y plane and the overlap extended over a minimum of four contiguous 0.2 μm slices in the confocal z-stack (supplemental video, available at www.jneurosci.org as supplemental material).
All results are given as mean ± SEM unless noted. Mean mEPSP amplitudes underwent a reciprocal transformation, and mean mEPSP frequencies underwent a logarithmic transformation to normalize the data and stabilize their variances before performing a two-way ANOVA. We modeled the dependence of connectivity, mean mEPSP amplitude, and unitary EPSP amplitude on age and deprivation status under the generalized linear model framework (McCullagh and Nelder, 1989) using the “stats” Package implemented in R (R Project for Statistical Computing, http://www.r-project.org/) with models of the following general form: g[E(Y)] = α + β × age + γ × DEP + δ × SP + ε × age × DEP + ζ × age × SP, where E(Y) is the expected value of connectivity, mean mEPSP amplitude, or unitary EPSP amplitude, g is a link function, DEP is a dummy variable (1 = deprived, otherwise 0) for deprived L2/3 neurons, and SP is a dummy variable (1 = spared, otherwise 0) for spared L2/3 pyramidal neurons. For simplicity, interactions between deprivation status and the duration of deprivation were modeled using age, which equals the duration of deprivation plus 19 d. The link function and distribution family, e.g., binomial, Gaussian, are specified for each model in Results.
The distributions of control and spared unitary EPSP amplitudes were well fitted by gamma distributions whose probability density functions had the following general form: g(x; r, λ, μ) = λr × (x − μ)(r−1) × exp(−λ × (x − μ))/Γ(r), where x is the unitary EPSP amplitude, r and λ are constants >0, μ is the offset constant, and Γ is the gamma function. The best fits were found by minimizing the Anderson–Darling statistic (A2), which places greater emphasis on the tails of distributions than the Kolmogorov–Smirnov test (Stephens, 1986), and was implemented in Matlab (MathWorks). The null hypothesis for the two-sample Anderson–Darling test was that control and spared unitary EPSP amplitudes were two random samples drawn from the same parent population, which we assumed also had a gamma distribution. The hypothesis was tested by finding the gamma distribution that maximized the probability that the spared and control unitary EPSP amplitudes were both drawn from it. The probability for each gamma distribution tested required the distribution of the Anderson–Darling statistic, which we found with Monte Carlo simulations that generated distributions with 100,000 data points. To verify our Monte Carlo distribution results, we compared them with published tables (Stephens, 1986) and found that they agreed to three significant figures.
Reduced connectivity in deprived cortex
We altered the sensory input to somatosensory cortex of rats after the whisker cortical map in L2/3 has developed (Stern et al., 2001) by trimming bilaterally the dorsal three (A–C) or the ventral two (D–E) rows of whiskers for 11–37 d (mean, 23 d). We then prepared brain slices that permitted visualization of the whisker barrel rows in somatosensory cortex (Fig. 1A). Sham-trimmed rats served as controls. We looked for evidence of altered connectivity by searching for pairs of synaptically connected L2/3 pyramidal neurons (Fig. 1B,C) that were within the same barrel column and close to the boundary between deprived and spared cortex using the loose-seal technique (Feldmeyer et al., 1999) (see Materials and Methods) (see Figs. 5A,B, 6A).
We performed 2683 loose-seal tests. The frequency of connections in control cortex was 1:35 (32 connections, 1117 loose-seals). A total of 22% (7 of 32) of these connections were reciprocal. Connectivity in spared cortex (frequency of 1:45, 20 connections identified from 902 loose seals; 3 (15%) reciprocal) was not changed compared with controls (χ2 = 0.60; df = 1; p = 0.44). In contrast, connectivity in deprived cortex was dramatically reduced (frequency of 1:83; 8 connections identified from 664 loose seals; 0% reciprocal) compared with controls (χ2 = 4.50; df = 1; p = 0.034). Addition of the glutamatergic antagonists CNQX (20 μm) and AP-5 (50 μm) to the recording bath blocked evoked EPSPs (n = 8 connections) (Fig. 1B). There was no difference between the groups in the age of the animals studied, duration of sensory deprivation, or passive membrane properties (Table 1).
The effect of age and the duration of deprivation (age minus 19 d) on connectivity in spared and deprived cortex over our recording period (P30–P60) was explored under the generalized linear model framework. We first used the following model: logit [E(connectivity)] = α + β × age + γ × DEP + δ × SP, with a logit link function and the binomial family (see Materials and Methods). The results suggested that connectivity between L2/3 pyramidal neurons was reduced in deprived cortex (γ = −0.951 ± 0.399; p = 0.017; z = −2.38) but not in spared cortex (δ = −0.272 ± 0.286; p = 0.342; z = −0.95) and that connectivity was not modified by age in isolation (β = −0.003 ± 0.019; p = 0.887; z = −0.14). We could not determine whether the effect of deprivation increased with the duration of deprivation by incorporating interaction terms into this model because of collinearity. Therefore, we analyzed the connectivity data using a second “interaction” model: logit [E(connectivity)] = α + γ × DEP × age + δ × SP. The interaction model did not result in a significant improvement in fit to the data (null deviance, 123; non-interaction model, residual deviance, 116; interaction model, residual deviance, 116) but was consistent with a progressive decrease in connectivity in deprived cortex with age (γ = −0.024 ± 0.010, p = 0.017, z = −2.38; α = −3.461 ± 0.179, p < 0.001, z = −19.3). The expected frequency of connections in deprived cortex at P30 was 1:66. This decreased to 1:132 at P60. Our findings indicate that sensory deprivation results in a selective loss of local excitatory connectivity in deprived cortex close to spared cortex. The statistical model using DEP as an explanatory variable effectively assumed that the effect of deprivation was independent of age and was as good at explaining the data as the model incorporating DEP × age. Hence, a reasonable interpretation of the statistical models is that a large part of the reduction in connectivity between L2/3 pyramidal neurons in deprived cortex had developed by P30 and that an additional reduction in connectivity probably occurred between P30 and P60.
Miniature EPSPs are larger in spared and deprived cortex
The reduced connectivity that we found indicated that pyramidal neurons in deprived cortex received less synaptic input from neighboring pyramidal neurons. Moreover, we found no compensatory modifications in either the excitability or the passive membrane properties of L2/3 pyramidal neurons (Fig. 1D,E, Table 1). Hence, we surmised that this combination of changes would manifest as a reduction in measures of synaptic input onto deprived pyramidal neurons.
Rats' whiskers were trimmed for 21–34 d, and brain slices were prepared. We blocked action potentials and fast inhibition to isolate mEPSPs (Fig. 2A) (see Materials and Methods), which are thought to comprise postsynaptic responses to release of single vesicles of excitatory neurotransmitter (Fatt and Katz, 1952). mEPSPs were recorded with the same parameters used for paired recordings to facilitate comparison with evoked unitary EPSPs of small amplitude. The mean frequency of mEPSPs in either spared cortex (10 ± 1 Hz; p = 0.72, t test; n = 10 neurons) or deprived cortex (9 ± 1 Hz; p = 0.27, t test; n = 10 neurons) was not different from control cortex (10 ± 1 Hz; n = 10 neurons) (Fig. 2B). The frequency of mEPSPs that we found is similar to values reported previously for recordings at physiological temperatures in younger animals (Hardingham and Larkman, 1998; Simkus and Stricker, 2002). However, the mean amplitude of mEPSPs in both spared cortex (0.18 ± 0.02 mV; p = 0.016, t test; n = 10 neurons) and deprived cortex (0.19 ± 0.02 mV; p = 0.030, t test; n = 10 neurons) was greater than in control cortex (0.13 ± 0.01 mV; n = 10 neurons) (Fig. 2C–F). The empirical distribution function of the mean amplitude of mEPSPs in deprived cortex was displaced to the right relative to controls (p = 0.031, Kolmogorov–Smirnov test; n = 10 neurons per group) (Fig. 2F). The empirical distribution function for mEPSP amplitude in spared cortex showed a similar rightward shift when compared with control cortex, but this did not reach statistical significance (Fig. 2E) (p = 0.11, Kolmogorov–Smirnov test; n = 10 neurons per group). The resting membrane potential and input resistance of neurons used for mEPSP analysis were similar in control, spared, and deprived cortex (data not shown), suggesting that our results were not attributable to differences in EPSP driving force or passive membrane properties of L2/3 pyramidal neurons.
Our recordings did not show the reduction in mEPSP frequency in deprived cortex that we had expected from the loss of local excitatory connectivity. We investigated this further by comparing mEPSPs with evoked unitary EPSPs (Fig. 1Bi). The range of amplitudes spanned by the 0 to 95th centiles of the control mEPSP distribution was 0.06–0.30 mV (Fig. 2C,D). The mean percentage of evoked unitary EPSPs with amplitudes within this range was 41% (n = 562 unitary EPSPs from 16 control connections). Hence, many evoked unitary EPSPs were similar in amplitude to mEPSPs. However, more detailed analysis revealed that the vast majority of evoked unitary EPSPs in L2/3 pyramidal neurons were larger than 0.10 mV (557 of 562 unitary EPSPs from 16 control connections). These findings indicated, first, that the smallest evoked unitary EPSPs were probably quantal responses and, second, that mEPSPs attributable to L2/3 inputs from neighboring pyramidal neurons were almost always larger than 0.10 mV. Many mEPSPs were smaller than this (Fig. 2C,D). Hence, it is possible that we did not detect a decrease in mEPSP frequency in deprived cortex because it was obscured by variability introduced by small-amplitude (<0.10 mV) mEPSPs. Therefore, we rethresholded our data to include only those mEPSPs larger than 0.10 mV in amplitude. We refer to these datasets as L2/3-biased mEPSPs. The mean frequency of L2/3-biased mEPSPs in either spared cortex (7 ± 1 Hz; p = 0.62, Mann–Whitney rank sum test; n = 10 neurons) or deprived cortex (7 ± 1 Hz; p = 0.66, Mann–Whitney rank sum test; n = 10 neurons) was not different from control cortex (7 ± 1 Hz; n = 10 neurons). The mean amplitude of L2/3-biased mEPSPs in spared cortex (0.22 ± 0.01 mV; p = 0.009, t test; n = 10 neurons) was greater than control cortex (0.17 ± 0.01 mV; n = 10 neurons). The mean amplitude of L2/3-biased mEPSPs in deprived cortex (0.20 ± 0.02 mV; p = 0.23, Mann–Whitney rank sum test; n = 10 neurons) was not statistically different from control values. This analysis supports the finding that mEPSP frequency was not reduced in deprived cortex despite loss of synaptic inputs from neighboring L2/3 pyramidal neurons.
The connectivity and mEPSP data together suggested that different processes were occurring in spared and deprived cortex. In spared cortex, the lack of change in local connectivity, unchanged mEPSP frequency, and increased mean mEPSP amplitude are consistent with strengthening of synapses without changes in synapse number. In deprived cortex, the preservation of mEPSP frequency in the presence of reduced local connectivity combined with the increased mean mEPSP amplitude in deprived cortex close to spared cortex pointed to a change in the source of excitatory inputs.
We explored the temporal and spatial development of the changes in mEPSPs in more detail with emphasis on deprived cortex. L2/3 pyramidal neurons in the C, D, and E barrel columns of the same brain slice were recorded after trimming the D and E whisker rows (Fig. 3A). In this experiment, whisker trimming started at P19 as before and continued for 13–20 d, which is shorter than the deprivation period used for the previous mEPSP study. Sham-trimmed rats served as controls. The frequency of mEPSPs in L2/3 of the spared C barrel column, the deprived D barrel column, and the deprived E barrel column were not different from the frequency of mEPSPs recorded from the corresponding barrel columns in control cortex (C-spared, 13 ± 1 Hz; C-control, 13 ± 1 Hz; p = 0.79; D-deprived, 13 ± 1 Hz; D-control, 13 ± 1 Hz; p = 0.99; E-deprived, 12 ± 1 Hz; E-control, 14 ± 1 Hz; p = 0.24, t tests; n = 10 neurons for all groups). The mean amplitude of mEPSPs was larger in L2/3 of the spared C barrel column and the deprived D barrel column but not in the deprived E barrel column (C-spared, 0.14 ± 0.01 mV; C-control, 0.11 ± 0.01 mV; p = 0.036, t test; D-deprived, 0.14 ± 0.01 mV; D-control, 0.11 ± 0.01 mV; p = 0.039, t test; E-deprived, 0.14 ± 0.01 mV; E-control, 0.12 ± 0.01 mV; p = 0.27, Mann–Whitney rank sum test; n = 10 neurons for all groups) (Fig. 3B–D). These data suggested that the synaptic changes in spared and deprived cortex were present after whiskers were trimmed for 13–20 d and that the increase in mEPSP amplitude in deprived cortex was most prominent at the junction of deprived and spared cortex.
Comparison of the results from animals that had been trimmed from P19 for 13–20 d with animals trimmed from P19 for 21–34 d suggested that mEPSP frequency decreased with age (p = 0.001; F(1,64) = 11.8; two-way ANOVA) but was not affected by deprivation status, i.e., spared, deprived, or sham trimmed (p = 0.85; F(2,64) = 0.17; two-way ANOVA). In contrast, mean mEPSP amplitude increased with age (p < 0.001; F(1,64) = 24.6; two-way ANOVA) and deprivation status (p < 0.001; F(2,64) = 8.5; two-way ANOVA) (see Materials and Methods).
We investigated the effect of the age, deprivation status, and duration of deprivation on mean mEPSP amplitude under the generalized linear model framework. We first used the following model: inverse[E(mean mEPSP amplitude)] = α + β × age + γ × DEP + δ × SP, with an inverse link function and the Gaussian family (see Materials and Methods). Mean mEPSP amplitudes were described reasonably well by the model (null deviance, 0.129; residual deviance, 0.088; n = 70 neurons) and suggested that mEPSP amplitude increased with both age and deprivation status (α = 13.1 ± 1.6, p < 0.001, t = 8.3; β = −0.12 ± 0.04, p = 0.001, t = −3.4; γ = −2.17 ± 0.48, p < 0.001, t = −4.1; δ = −1.92 ± 0.55, p < 0.001, t = −3.5). We could not ascertain whether the effect of deprivation progressed during the deprivation protocol by incorporating interaction terms into this model because of collinearity. Therefore, we analyzed the mEPSP data using a second interaction model (inverse link function and the Gaussian family): inverse[E(mean mEPSP amplitude)] = α + β × age + γ × DEP × age + δ × SP × age. The interaction model did not result in a significant improvement in fit to the data (residual deviance, 0.087) but was consistent with a progressive increase in mean mEPSP amplitude in spared and deprived cortex with longer deprivation (γ = −0.053 ± 0.013, p < 0.001, t = −4.24; δ = −0.045 ± 0.013, p < 0.001, t = −3.47; α = 12.0 ± 1.6, p < 0.001, t = 8.0). The effect of deprivation was superimposed on a slow increase in mEPSP amplitude attributable to age alone (β = −0.09 ± 0.04, p = 0.011, t = −2.6).
Our analysis suggested that the increase in mean mEPSP amplitude in L2/3 pyramidal neurons has at least two components. First, there is a normal, slow increase in mEPSP amplitude between P32 and P53. This occurs in both sham-trimmed and whisker-trimmed animals. Superimposed on this normal increase is a second component attributable to deprivation, which further increases mEPSP amplitude at the border between spared and deprived cortex. The statistical models indicated that an effect of deprivation had developed by the time that we started recording and that this probably progressed over the recording period. In our animals, the deprivation component had a larger effect on mEPSP amplitude than the physiological increase with age [Kolmogorov–Smirnov tests, P32–P39 group compared with P40–P53 group: control, p = 0.18 (Fig. 4A); spared, p = 0.001 (Fig. 4B); deprived, p = 0.001 (Fig. 4C)].
Increased unitary EPSP amplitude in spared cortex
Our data suggested strengthening of synapses in spared cortex but could not specify the origin of the strengthened inputs. L2/3 pyramidal neurons receive much of their excitatory input from neighboring pyramidal neurons (Binzegger et al., 2004). Therefore, we investigated whether sensory deprivation modified the EPSP amplitude at unitary connections between pairs of pyramidal neurons in the C or D barrel columns at the junction of deprived and spared cortex (see Figs. 5A,B, 6A). A change in unitary EPSP amplitude could be attributable to alterations in synapse strength or synapse number. Hence, we used the short-term synaptic dynamics at unitary connections as a second probe to assay whether there were functional changes at synapses. The mean baseline unitary EPSP amplitude in control cortex was 0.31 ± 0.07 mV (n = 32 connections) (Fig. 5C). During 20 Hz stimulation, unitary EPSPs facilitated, depressed, or showed a combination of facilitation and depression (Fig. 5D). In spared cortex, baseline unitary EPSP amplitudes (0.53 ± 0.12 mV; n = 21) were larger than controls (p = 0.027, Mann–Whitney rank sum test). Unitary EPSP amplitudes in control and spared cortex were well described by gamma distributions [control(x; 1.24, 4.95, 0.03), p = 0.916; spared(x; 0.96, 2.13, 0.06), p = 0.991; Anderson–Darling test; see Materials and Methods] (Fig. 5E). The empirical distribution plot showed that unitary amplitudes in spared cortex were tilted toward higher values relative to controls (p < 0.046, two-sample Anderson–Darling test; see Materials and Methods). Moreover, unitary EPSPs in spared cortex depressed more during 20 Hz trains compared with controls (spared mean change, −0.16 ± 0.06 mV, n = 21; control mean change, −0.04 ± 0.03 mV, n = 32; p = 0.020, Mann–Whitney rank sum test) (Fig. 5F,G). Despite the increased short-term depression, excitatory drive was greater throughout 20 Hz trains (Fig. 5F). The empirical distribution plot of change in unitary EPSP amplitude at steady state showed a leftwards shift relative to controls (p = 0.036, Kolmogorov–Smirnov test) (Fig. 5G), suggesting that there is greater depression at all unitary connections in spared cortex rather than a selective change in synaptic dynamics at a subpopulation of unitary connections. The ages of sham-trimmed rats and trimmed rats were not different (Table 1). There was no difference in the locations of presynaptic or postsynaptic neurons recorded in spared or control cortex (Fig. 5A,B).
We investigated whether the increase in unitary EPSP amplitude was affected by age and the duration of deprivation (Fig. 5H) with the following model: ln(unitary EPSP amplitude) = α + β × age + γ × age × SP, which used a natural log link function and Gaussian family (see Materials and Methods). Unitary EPSP amplitude in spared cortex increased with the duration of deprivation (δ = 0.019 ± 0.007; t = 2.61; p = 0.012) but was not affected by age in isolation (β = 0.033 ± 0.021; t = 1.61; p = 0.114). We concluded that synapses at local excitatory connections in spared cortex were stronger after deprivation and that this strengthening continued to develop over the deprivation period.
Connection strength and synaptic dynamics do not alter in deprived cortex
Unitary connections in deprived cortex gave very different results from spared cortex. The baseline unitary EPSP amplitude (0.39 ± 0.09 mV; p = 0.220, Mann–Whitney rank sum test; n = 8 connections) (Fig. 6B), the change in EPSP amplitude at steady state (mean change, −0.08 ± 0.03 mV; p = 0.106, Mann–Whitney rank sum test) (Fig. 6C), and its distribution (p = 0.18, Kolmogorov–Smirnov test) (Fig. 6D) were not different from controls. The paucity of connections recorded from deprived cortex means that our results have to be interpreted with caution but indicate that there were no large-scale changes in the strength of the remaining local excitatory connections in deprived cortex.
Synapse number is unchanged at local excitatory connections in spared cortex
At least part of the increase in unitary EPSP amplitude and greater synaptic depression in spared cortex could have been attributable to an increase in synapse number per unitary connection with new synapses showing greater synaptic depression. Therefore, we used laser scanning confocal microscopy to reconstruct pairs of synaptically connected neurons filled with fluorescent dye during recording (Fig. 1C) and counted putative synapses between them (Fig. 7A,B) (supplemental video, available at www.jneurosci.org as supplemental material) (see Materials and Methods).
We found no difference in the mean number of synapses per unitary connection in spared (3.0 ± 0.9; n = 8 connections) and control (3.5 ± 0.6; p = 0.72, t test; n = 11 connections) cortex. The proportions of synapses on apical and basal dendrites formed by neighboring pyramidal neurons in spared cortex were not different from controls (apical:basal synapses, control, 22:16; spared, 11:13; χ2 = 0.44; df = 1; p = 0.51) (Fig. 7C). The locations of presynaptic and postsynaptic neurons used for synapse counts were representative of their respective groups (data not shown). We estimated the error associated with our optical criteria for identifying putative synapses by reconstructing either the unconnected direction from pairs of L2/3 pyramidal neurons that were unidirectionally connected (n = 8) or pairs of unconnected L2/3 pyramidal neurons (n = 2). The number of points of close proximity fulfilling our criteria for putative synapses was 0.17 per L2/3 axon (n = 12) unconnected to an adjacent L2/3 pyramidal neuron. This is ∼5% of the synapse number found for synaptically connected pairs of neurons. The cell bodies of unconnected neurons were separated by the same distance as the cell bodies of connected pairs (intersomatic distance: unconnected pair, 60 ± 5 μm; connected pair, 63 ± 5 μm; p = 0.75, Mann–Whitney rank sum test). Hence, the frequency of false-positive results returned by our optical criteria is within acceptable limits. We concluded that deprivation did not increase synapse number at connections between neighboring L2/3 pyramidal neurons in spared cortex.
We were able to reconstruct two pairs of synaptically connected L2/3 pyramidal neurons in deprived cortex, which formed a total of four putative synapses (Fig. 7C). The small sample size meant that we could not draw any reliable conclusions about structural reorganization in deprived cortex from these data.
Our study shows that whisker deprivation induces site-specific modifications in local excitatory networks in primary somatosensory cortex. We found that there was a dramatic reduction in connectivity between L2/3 pyramidal neurons in deprived cortex, which was not accompanied by a compensatory increase in the strength of remaining local excitatory network connections. A markedly different situation occurred in adjacent spared cortex in which there was strengthening of connections between neighboring pyramidal neurons without a change in synapse number per connection or local excitatory connectivity.
Reduction of local excitatory connectivity in deprived cortex
The precise cellular mechanisms underlying the marked reduction in connectivity that we report in deprived cortex remain to be elucidated. Our analysis indicated that the structural reorganization of local excitatory circuitry was well developed by P30. This made it difficult to be certain whether the reorganization was a global process affecting all connections or whether specific connections were targeted.
The lack of change in mEPSP frequency in deprived cortex suggests that one outcome of the reorganization process may be replacement of synaptic inputs previously formed by neighboring pyramidal neurons in deprived cortex by other excitatory inputs. Our recordings indicated that the increase in mean mEPSP amplitude was not uniform throughout deprived cortex but was greatest at the junction of deprived and spared cortex. Hence, one likely source for new excitatory synapses is horizontal intracortical L2/3 inputs from spared cortex, which have been reported to increase their axonal arborization into deprived cortex (Darian-Smith and Gilbert, 1994; Kossut and Juliano, 1999). These inputs strengthen after whisker deprivation (Finnerty et al., 1999). Therefore, it is possible that formation of new, stronger excitatory synapses by L2/3 inputs from spared cortex with pyramidal neurons in deprived cortex drives the loss of synapses between neighboring pyramidal neurons in deprived cortex. This hypothesis predicts that there should be increased synaptic inputs from L2/3 pyramidal neurons in spared cortex onto L2/3 pyramidal neurons in deprived cortex. We could not test this hypothesis directly because unitary connections between L2/3 pyramidal neurons in C and D barrel columns were difficult to find (no connections in 221 tests) in our brain slice preparation. This is not surprising given the distance (>120 μm) between potential recording sites (Holmgren et al., 2003). Hence, additional work possibly combined with technical advances will be needed to test this prediction.
It has been proposed that synapse elimination in cortex resembles the process at immature neuromuscular junctions in which motor afferents compete to innervate a muscle cell (Lichtman and Colman, 2000). One of the motor inputs strengthens, with concurrent weakening of the other inputs before axon withdrawal of the weakened inputs (Colman et al., 1997). This hypothesis predicts that (1) unitary EPSP amplitudes of local excitatory connections in deprived cortex would decrease and (2) there would be a left shift in the empirical distribution function of unitary EPSP amplitudes compared with control unitary EPSP amplitudes. Our data do not support this hypothesis. However, the low connectivity in deprived cortex means that we found few connections. Hence, this hypothesis cannot be excluded and requires additional investigation with recordings before connections are lost.
The reduced connectivity in deprived cortex that we found is opposite to the increased connectivity in layer 4 of deprived visual cortex induced by monocular deprivation during the early critical period (Maffei et al., 2004). The difference in the results suggests that the effect of sensory deprivation on connectivity depends on several factors, including developmental status.
Strengthening of local excitatory connections in spared cortex
Our data indicated that functional changes (progressive increase in unitary EPSP amplitude and greater synaptic depression) occurred at unitary connections between L2/3 pyramidal neurons in spared cortex without an alteration in the number of putative synaptic contacts. These findings suggest that whisker deprivation in developed animals primarily causes strengthening of excitatory unitary connections in spared cortex by modifying the properties of excitatory synapses at those connections rather than by increasing synapse number. The most parsimonious explanation for our results is that whisker trimming induced functional changes at existing synapses. The changes in unitary EPSP amplitude and synaptic dynamics that we found in spared cortex are similar to those that accompany long-term potentiation (LTP) induced by pairing presynaptic action potentials with postsynaptic depolarization in vitro (Markram and Tsodyks, 1996). This raises the possibility that sensory deprivation strengthens local excitatory connections in spared cortex in vivo by a process that resembles LTP in vitro. However, we cannot exclude the possibility that synapse turnover (Trachtenberg et al., 2002; Zuo et al., 2005a; De Paola et al., 2006; Stettler et al., 2006) results in new, stronger synapses being formed during the deprivation protocol.
Strengthening of local excitatory connections may play a role in cortical map expansion by either increasing firing within local excitatory networks in L2/3 of spared cortex or by altering the pattern of firing within those networks, e.g., by increasing synchronization. This would augment excitatory drive along horizontal intracortical pathways and, thereby, enhance expansion of spared representations into supragranular layers of deprived cortex (Das and Gilbert, 1995; Buonomano and Merzenich, 1998; Finnerty et al., 1999) or modify the representation of spared inputs within spared cortex (Diamond et al., 1994).
The effects of whisker deprivation on neuronal firing vary with the duration and pattern of deprivation (Fox, 2002). A commonly used method to study this involves mechanically displacing one whisker at low frequency. Evoked firing in L4 of a barrel column receiving its principal input from a spared whisker increases after whisker pairing (Armstrong-James et al., 1994) and chessboard deprivation (Wallace and Fox, 1999). Alterations in evoked principal whisker responses in L2/3 of the spared barrel column are not as prominent as those in L4 but include increased firing in a proportion of neurons (Barth et al., 2000; Wallace et al., 2001; Fox, 2002). In contrast, there is a robust increase in firing in L2/3 of deprived barrel columns (Glazewski and Fox, 1996; Wallace and Fox, 1999) and spared barrel columns (Diamond et al., 1994) that are adjacent to the barrel column receiving its principal input from the spared whisker. The latter finding suggests that there may be an interaction between spared whiskers. The magnitude of this effect in our preparation is not clear, but it raises the possibility that results from deprivation paradigms, which spare multiple whiskers, may not generalize to rodents with single-whisker experience.
The effects of deprivation on mEPSP amplitude were superimposed on an age-dependent increase in mEPSP amplitude and decrease in mEPSP frequency. Spine density on neocortical neurons peaks between P15 and P30 before declining slowly (Wise et al., 1979; Markus and Petit, 1987; Zuo et al., 2005a). The decrease in mEPSP frequency that we found with age is consistent with the physiological spine loss that occurs in developed animals. The increase in mEPSP amplitude may reflect either loss of weaker synapses or strengthening of some synapses over our recording period.
The effects of deprivation on mEPSPs that we found in developed animals appear to differ from a study in the developing primary somatosensory cortex, which reported that mEPSC amplitudes in L2/3 pyramidal neurons after D-row deprivation were not different in the deprived D row compared with the spared B row (Bender et al., 2006). However, dissimilarities in the deprivation protocols make direct comparison difficult. In particular, experiments on visual cortex suggest that the timing of the deprivation protocol with respect to maturation of inputs to L4 has a marked effect on the outcome of deprivation-induced changes in mEPSCs recorded from L2/3 pyramidal neurons (Desai et al., 2002).
Local excitatory circuitry and whisker representations
We propose a simple model whereby changes in local excitatory networks could contribute to cortical map reorganization. This model incorporates the modifications in longer-distance horizontal connections reported previously (Finnerty et al., 1999) and complements findings from experiments on vertical intracortical connections (Finnerty et al., 1999; Allen et al., 2003; Bender et al., 2006; Clem and Barth, 2006). The alterations in responses at the junction of deprived and spared cortex induced by sensory deprivation are activity dependent (Wallace et al., 2001). In freely behaving rats, sensory deprivation is followed acutely by increased firing in spared cortex and decreased firing in deprived cortex (Kelly et al., 1999). We propose that the increased firing of pyramidal neurons in L2/3 of spared cortex has two effects. First, it results in strengthening of existing excitatory synapses in spared cortex, which manifests as potentiation of local excitatory connections. Second, increased firing or changes in the pattern of firing in L2/3 of spared cortex will augment the excitatory drive from spared to deprived cortex. We suggest that this promotes loss of excitatory synapses in deprived cortex. Hence, depression of sensory responses in deprived cortex and potentiation of sensory responses in spared cortex can be dissociated by manipulations that block synaptic potentiation (Glazewski et al., 1996) but do not affect the acute changes in firing. The apparently greater dependence of plasticity in deprived cortex on changes in wiring may explain why depression of deprived responses is abolished in adult animals, whereas potentiation of spared responses persists (Fox, 2002).
The work was funded by the Wellcome Trust, King's College Hospital Medical Research Trust (C.E.J.C.), and the Medical Research Council (C.E.J.E.). G.T.F. is a Wellcome Trust Senior Fellow in the Clinical Sciences. We thank Barry Connors, John Parnavelas, and Ben Philpot for helpful comments.
- Correspondence should be addressed to Gerald T. Finnerty at the above address.