Abstract
Cortical maps are remarkably precise, with organized arrays of thalamocortical afferents (TCAs) that project into distinct neuronal modules. Here, we present evidence for the involvement of efficient neurotransmitter release in mouse cortical barrel map development using barrelless mice, a loss-of-function mutant of calcium/calmodulin-activated adenylyl cyclase I (AC1), and mice with a mutation in Rab3-interacting molecule 1α (RIM1α), an active zone protein that regulates neurotransmitter release. We demonstrate that release efficacy is substantially decreased in barrelless TCAs. We identify RIMs as important phosphorylation targets for AC1 in the presynaptic terminal. We further show that RIM1α mutant mice have reduced TCA neurotransmitter release efficacy and barrel map deficits, although not as severe as those found in barrelless mice. This supports the role of RIM proteins in mediating, in part, AC1 signaling in barrel map development. Finally, we present a model to show how inadequacies in presynaptic function can interfere with activity-dependent processes in neuronal circuit formation. These results demonstrate how efficient synaptic transmission mediated by AC1 function contributes to the development of cortical barrel maps.
Introduction
The mouse primary somatosensory cortex is characterized by a row-and-arc pattern of “barrels” in cortical layer IV, replicating the arrangement of whiskers on the snout (Woolsey and Van der Loos, 1970; Killackey and Leshin, 1975). Each cylindrical barrel is composed of two parts: a presynaptic component with aggregates of thalamocortical afferents (TCAs) that relay sensory information from a single whisker; and a postsynaptic component with discrete rings of layer IV neurons that project their dendrites inward to synapse on one cluster of TCAs.
The precise mapping of facial whiskers to cortical barrels requires the coordinated development of TCAs and the dendrites of postsynaptic cortical neurons. Several lines of evidence suggest that input from the sensory periphery is instructive in this process. First, in mutant mice with an atypical number of whiskers, there is close correspondence between the barrel pattern in the cortex and the pattern of whiskers on the snout (Van der Loos and Dorfl, 1978; Van der Loos et al., 1984). Second, embryonic visual cortex transplanted to a location corresponding to the primary somatosensory cortex forms barrels (Schlaggar and O’Leary, 1991). Third, knocking out the serotonin transporter, which is expressed in TCAs but not in cortical neurons, leads to a disruption of the barrel map (Persico et al., 2001). Finally, presynaptic 5-HT1B receptors on TCAs play an essential role in the barrel phenotype of monoamine oxidase A mutants (Salichon et al., 2001; Rebsam et al., 2002). Despite the evidence for an instructive role of the sensory periphery in barrel map development, an integral role for the cortex in this process has also been well documented. For example, introducing an ectopic source of the morphogen fibroblast growth factor 8 (FGF8) into posterior cortex produces a duplicate barrel map (Fukuchi-Shimogori and Grove, 2001), and the cortex-specific deletion of the NMDA receptor subunit NR1 severely disrupts barrel map development (Iwasato et al., 2000; Datwani et al., 2002; Lee et al., 2005). Thus, it is likely that interplay between cortical neurons and TCAs regulates the emergence of a barrel pattern.
Barrelless (brl) (Adcy1brl) mice have a loss-of-function mutation in the adenylyl cyclase I (AC1) gene (Abdel-Majid et al., 1998). In barrelless, TCAs reach their cortical target but do not cluster into barrels (Van der Loos et al., 1986; Welker et al., 1996), despite normal subcortical barrel structures in the somatosensory thalamus and brainstem. In vivo single-cell recordings from the presumptive barrel field of barrelless mice show that layer IV neurons have somatotopically appropriate receptive fields but respond inappropriately to stimulation of more than one facial whisker. We described previously how deficiencies in postsynaptic AMPA receptor (AMPAR) trafficking and long-term synaptic plasticity at thalamocortical (TC) synapses in barrelless mice likely contribute to their phenotype (Lu et al., 2003). However, barrel defects in barrelless mice are quite severe, with no evidence of a barrel pattern in TCAs or cortical layer IV neurons. This is unexpected if the locus of the deficit is purely postsynaptic, as other barrel mutant mice with presumed postsynaptic deficits, such as cortex-specific NMDA receptor subunit 1 (Cx-NR1) mutant mice (Iwasato et al., 2000) and phospholipase C-β1 (PLC-β1) mutant mice (Erzurumlu and Kind, 2001; Hannan et al., 2001), show much more severe abnormalities in cortical layer IV neuron barrel patterning than in TCA clustering. The complete absence of a barrel pattern in barrelless mice suggests functional deficits may also exist in the presynaptic compartment of TC synapses in barrelless mice.
We present electrophysiological evidence for neurotransmitter release defects in barrelless mutant mice. We then use biochemical techniques to demonstrate that reduced AC1-dependent phosphorylation of Rab3-interacting molecules (RIMs), active zone proteins critical in regulating neurotransmitter release, contributes to these deficits in barrelless mice. We further show that RIM1α mutant mice have TCA neurotransmitter release defects and impaired barrel maps, confirming that the efficient release of neurotransmitter contributes to barrel map development. Finally, we present a computational model that demonstrates how inadequacies in presynaptic function can disrupt activity-dependent segregation of sensory afferents into an organized topographic map.
Materials and Methods
Electrophysiology.
Acute thalamocortical brain slices, standard artificial CSF (ACSF), and whole-cell “internal” solutions were prepared as described previously (Lu et al., 2001). The ACSF was saturated with 95% O2 and 5% CO2 and consists of the following (in mm): 124 NaCl, 5 KCl, 1.25 NaH2PO4, 1.3 MgSO4, 2 CaCl2, 26 NaHCO3, and 11 glucose, pH 7.2, 300–310 mOsm. Experiments were performed at room temperature (∼26°C). The whole-cell recording solution contained the following (in mm): 99 cesium gluconate, 17.5 CsCl, 8 NaCl, 10 HEPES, 0.2 EGTA, 4 Mg-ATP, 0.3 GTP, 7 phosphocreatine, and 10 BAPTA, pH 7.2, 290–300 mOsm. BAPTA was included to prevent inadvertent potentiation in the postsynaptic compartment. Whole-cell voltage-clamp recordings were used to record from layer IV neurons in somatosensory cortex using a MultiClamp 700A (Molecular Devices, Union City, CA). Stimuli (30–300 μA; 0.1 ms) were applied through bipolar sharpened and insulated stainless steel microelectrodes (FHC, Bowdoinham, ME) placed into the ventrobasal thalamus. Data were collected with 10 kHz sampling with 5 kHz filter and analyzed on-line using a computer-driven acquisition system (PCS with National Instrument AD boards; National Instruments, Austin, TX) and software that was written within the Igor programming environment. Input and series resistances were measured continuously to evaluate and monitor the health of the cell. Responses were considered to be monosynaptic only if they exhibited short and constant latencies that did not change with increasing stimulus intensity. To eliminate inhibitory currents, 100 μm picrotoxin or 10 μm GABAzin (Tocris Cookson, Ellisville, MO) was included in the ACSF. To measure NMDA currents, the cell was held at +40 mV and 10 μm 2,3-dihydro-6-nitro-7-sulfamoyl-benzo(f)quinoxaline (NBQX; Tocris Cookson) was used to block AMPARs. For (+)-5-methyl-10,11-dihydro-5H-dibenzo [a,d] cyclohepten-5,10-imine maleate (MK801) experiments, NMDA current was isolated and only cells responding with a significant amount of NMDA current (>15 pA) were used. After acquiring a stable baseline response (15–20 min), the stimulation was stopped and ACSF containing (+)-MK801 (40 μm; Tocris Cookson) was perfused onto the slice (1 ml/min flow rate). After 10 min (to allow enough time for MK801 to equilibrate in the recording chamber), we resumed stimulation and recorded until no NMDA current could be detected for >5 min. The time course of MK801-dependent blockade was analyzed by normalizing all of the responses to the amplitude of the first NMDA response after MK801 was applied, and the blockade rates were estimated by fitting an exponential function to the decay.
For calcium-dependence experiments, NMDA-EPSCs were isolated and evoked at 0.2–0.33 Hz for a period of 15–20 min at 2 mm [Ca2+] at first to achieve a stable baseline for NMDA-EPSCs, and then the external [Ca2+] was sequentially decreased to 1.75–0.5 mm [Ca2+] while the [Mg2+] was adjusted to maintain appropriate cation concentrations. The combinations of the CaCl2 and MgSO4 concentration were as follows: 2 mm Ca2+/1.3 mm Mg2+, 1.75 mm Ca2+/3 mm Mg2+, 1.5 mm Ca2+/3 mm Mg2+, 1.25 mm Ca2+/4 mm Mg2+, 1 mm Ca2+/5 mm Mg2+, 0.75 mm Ca2+/5 mm Mg2+, and 0.5 mm Ca2+/5 mm Mg2+.
Data were accepted and analyzed only when neurons have an input resistance >300 MΩ, and the synaptic response was stable with no sign of drift for at least 10 min before any treatment. EPSC amplitudes were defined as the mean current during a fixed 3–4 ms window at the peak of the EPSC minus the mean current during a similar window immediately before the stimulus artifact. All summary data are presented as mean ± SEM. The Student’s t test was used to determine statistical significance, except where otherwise specified. A two-way ANOVA was used to compare the PP ratios and calcium-dependent release. Differences were considered significant if p < 0.05.
Animals and genotyping.
Barrelless mice were discovered as a spontaneous mutation in a line from the ICR stock at Université de Lausanne (Van der Loos et al., 1986). Barrelless mice of the 10th backcross generation of the incipient C57BL/6-Adcy1brl congenic inbred strain were used in these studies. For all of the electrophysiological experiments in this study, control (ctrl) animals were the wild-type (wt) and heterozygous (het) littermates of barrelless mice with the same genetic background. No significant differences between wt and het animals were detected, and thus for most of the data comparisons, the data points from wt and het animals were lumped together as the control group. The loss-of-function mutants of RIM1α and synapsin I/II were generated in Dr. T. C. Südhof’s laboratory (Rosahl et al., 1995; Castillo et al., 2002). Genotypes of these mice were determined as described previously (Rosahl et al., 1995; Castillo et al., 2002; Lu et al., 2003). Animals were treated in compliance with the United States Department of Health and Human Services and Baylor College of Medicine guidelines.
Immunoblotting.
Somatosensory cortices were isolated from postnatal day 4 (P4) to P11 mice and homogenized in homogenization buffer consisting of 0.32 m sucrose, 10 mm HEPES, pH 7.4, with protease and phosphatase inhibitor mixture (Sigma, St. Louis, MO). Total tissue homogenates were prepared as described previously (Lu et al., 2003). Twenty to 30 μg of total protein were loaded onto 5–20% gradient SDS-PAGE gels and electrophoretically transferred to nitrocellulose membranes (Bio-Rad, Hercules, CA). Primary antibody against the protein in question was applied, and the signals were detected with HRP-conjugated secondary antibodies. The immunoreactivity was quantified with ECL (Amersham Biosciences, Piscataway, NJ) and densitometric quantification with NIH Image software (Scion, Frederick, MD). A paired t test was used for most of the immunoblotting studies.
Histology.
The cytochrome oxidase (CO) and Nissl staining of tangential sections of somatosensory cortex were conducted as described previously (Lu et al., 2001). Briefly, a fixed, flattened cortex from 4% paraformaldehyde-perfused brains (P7 to adult) were sectioned tangentially (parallel to layer IV) with a Leica (Nussloch, Germany) Vibratome to 100 or 50 μm thickness and placed in PBS. For CO staining, 100 μm tangential slices were placed in CO reaction solution (5 mg of diaminobenzidine, 5 mg of cytochrome c, 0.4 g of sucrose in 10 ml of PBS, pH 7.4) for 4–18 h at room temperature. After staining, slices were rinsed with PBS and coverslipped with Fluoromount-G (EMS, Hatfield, PA). For Nissl staining, serial 50-μm-thick tangential sections were mounted on Superfrost-plus slides (Fisher Scientific, Hampton, NH). After drying, the entire series of sections throughout the cortical depth were stained with cresyl violet solutions. Sections were rinsed, dehydrated in graded alcohols to xylene, and then coverslipped with xylene-based mounting medium (Cytoseal; Richard-Allen Scientific, Kalamazoo, MI).
Antibodies.
Primary antibodies were as follows: synapsin (Synaptic Systems, Göttingen, Germany), synapsin phosphorylation site-specific antibodies (Jovanovic et al., 2001), RIMs (BD Biosciences, San Diego, CA), p-RIM(PKA) (Lonart et al., 2003), p-rabphilin and rabphilin (Lonart et al., 2001), and β-actin (Sigma). Secondary antibodies were as follows: goat anti-mouse IgG-Cy3, I-125 mouse IgG, and rabbit IgG were obtained from Amersham Biosciences; rabbit IgG Alexa 488 was obtained from Invitrogen (Eugene, OR), and mouse IgG-HRP and rabbit-IgG-HRP were obtained from Pierce (Rockford, IL).
Model.
We simulate a single postsynaptic neuron that receives input from 30 thalamic afferents from two whiskers (see Fig. 7A): the primary whisker (PW) and an adjacent whisker (AW). Barrel formation in this case corresponds to the evolution of the simulation such that afferents from the AW are eliminated. Because we are interested in discovering a situation in which this segregation fails, we initially bias the competition in favor of the PW with 20 initial contacts to the 10 contacts from the AW, so that segregation failing in this case is likely representative.
Although early postnatal thalamic activity in vivo has not been directly measured, much is known about adult thalamic activity (Simons and Carvell, 1989). Because thalamic development generally precedes cortical development, we assume thalamic activity during cortical barrel development is similar to the adult and is comprised of both isolated action potentials and bursts stimulated by vibrissa movement. Thalamic afferents fire isolated spikes with uniform probability at a rate of 8 Hz. Bursts involving all afferents of one whisker occur randomly at a rate of 0.05 Hz. During a burst, afferents fire a spike within 5 ms of burst onset and then randomly at a rate of 36 Hz for 100 ms.
To simulate synaptic depression, we use synapses with several active zones (NAZ = 10). Each active zone has a single vesicle, which it releases with probability Pr during a presynaptic action potential. Vesicles are replenished a time TR = 100 ms after release, and a time TB is chosen so that a maximum of one vesicle per active zone releases during a burst. With these assumptions, the probability of release is simply related to the paired-pulse (PP) ratio (PPR): PPR = 1-Pr.
We assume that thalamic activity causes some form of feedforward inhibition, likely mediated through intracortical pathways, to limit cortical neuronal response to thalamic input to one or two action potentials. In the adult, feedforward inhibition is prominent during both whisker movement and direct electrical stimulation in the thalamus and results in single-spike cortical response despite bursting thalamic input (Pinto et al., 2003). For simplicity, we assume that feedforward IPSPs occur regardless of whether a given active zone releases, with a 10 ms delay. EPSPs and IPSPs are assumed to be uniform over their duration, with EPSPs lasting TE = 10 ms and IPSPs lasting TI = 50 ms. The size of EPSPs in a given time step of length Δt is Wsyn(i) × (Δt/TE), where Wsyn(i) is the synaptic weight of the ith active zone, and the size of IPSPs is given by winh × (Δt/TI), where winh is the strength of inhibitory synapses. Unlike the excitatory weight Wsyn(i), the inhibitory weight winh is not plastic and is chosen such that the strength of inhibitory inputs at the beginning of the simulation is equal to 0.5 times the strength of excitatory inputs: winh = 0.5PrNAZW0, where W0 is the initial weight of excitatory synapses. Excitatory weights are chosen to be uniform across all synapses at the beginning of the simulation.
The cortical neuron is modeled as a simple integrate-and-fire neuron, the membrane potential of which at a time t is given by V(t) = V(t − Δt) exp(−Δt/τM) + (synaptic input), where τM = 100 ms is the membrane time constant. When the membrane potential exceeds the threshold θ, the neuron fires an action potential, and V is reset to zero. The threshold θ is chosen for a given simulation so that the postsynaptic firing rate is equal to the rate of all presynaptic bursts (0.1 Hz).
We use a spike time-dependent learning rule (see Fig. 7B) like that measured at cortical layer IV to layer II/III synapses of somatosensory cortex (Feldman, 2000). If a presynaptic vesicle release precedes a postsynaptic spike by a time t, the synaptic weight is strengthened by A+exp(−t/τ+), and if a vesicle release follows a postsynaptic spike, the associated synaptic weight is weakened by A−exp(−t/τ−). Here, τ+ = 10 ms, τ− = 50 ms, and A+/A− = 0.52. The effect of the learning rule over time is mediated by the learning rate A+, which ranged from 0.0001 to 0.001 with no change in the behavior of the simulation. Synaptic weights cannot fall below 0, nor can they exceed a maximum value given by θ/(2PrNAZ), such that, at their strongest, an average of two simultaneously spiking afferents are needed to evoke a postsynaptic spike. For the simulations shown in Figure 7F, PPR and the strength of LTD (A−) were held constant and the long-term potentiation (LTP)/long-term depression (LTD) ratio was varied as shown.
The simulation was run at a time step of Δt = 1 ms, which was sufficiently small relative to the other simulation time scales (the learning rule and burst dynamics). The simulation ran until either the AW weights went to 0 (segregation) or the layer IV neuron stopped firing. For sufficiently high release probability (Pr > 0.4, or PPR < 0.6), PW afferents strengthened and AW afferents weakened, meaning the simulation segregated (see Fig. 7C–E). Segregation occurs as a result of a balance of potentiation and depression such that PW afferents preferentially strengthen from bursting activity and AW afferents weaken from spontaneous activity. For coincidence to lead to synaptic strengthening, it is necessary that Pr A − Pr (1 − Pr) B1 − Pr (1 − Pr)2 B2 − Pr(1 − Pr)3 B3 > 0, where A, B1, B2, and B3 are approximate values from the spike time-dependent learning rule at the time of the first, second, and third spikes, respectively, and Pr is the probability of release. This condition is not met in the simulation for Pr < 0.4, meaning that bursts no longer lead to synaptic strengthening, upsetting the balance between potentiation and depression, and resulting in weakening of both PW and AW afferents until they can no longer drive cortical spikes.
The length of time that it takes to reach either outcome depends on the learning rate, but higher or lower values of the learning rate did not change the outcome of the simulation. As a result, the simulation time is an arbitrary number and is not labeled on the axes of Figure 7. The simulation was written in ANSI c++, and the code is available on request.
The behavior of this simulation depends on several parameter choices that are not constrained by experiment. As a result, we verified that its qualitative behavior could be obtained for a range of such unconstrained parameter values. At the same time, we found that some general aspects of the model structure were required for our findings. In particular, it is necessary to assume some form of cortical inhibition, which prevents the cortical neuron from bursting during thalamic bursts, consistent with observations of layer IV in vivo activity at the earliest ages studied (Foeller and Feldman, 2004). This inhibition may be GABAergic, although it may well be mediated through other mechanisms early in development such as presynaptic inhibition (Rhoades et al., 1994). Without inhibition, cortical neuron spiking will always follow strong presynaptic input, meaning that bursts will always lead to strengthening, and decreasing Pr will not prevent segregation.
Results
Impaired short-term synaptic plasticity in barrelless TCAs
We first used a paired-pulse protocol to characterize neurotransmitter release properties of TC synapses in developing barrelless and control neurons. Paired-pulse (PP) ratios are a standard measure of presynaptic function that is typically inversely related to the probability (Pr) of neurotransmitter release (Zucker and Regehr, 2002). PP ratios are determined by measuring the relative response amplitude to two closely spaced stimuli. Synapses with high Pr regularly show paired-pulse depression (PPD), whereas synapses with low Pr show paired-pulse facilitation. TC synapses are glutamatergic and contain a mixture of AMPARs, kainate receptors, and NMDA receptors (NMDARs) (Kidd and Isaac, 1999; Lu et al., 2001). We initially used NMDAR-mediated EPSCs to examine PP ratios at TC synapses with a range of interstimulus intervals (ISIs; 50, 75, 100, and 500 ms). In control cells, paired stimuli produced strong PPD (Fig. 1A,C), consistent with the high Pr found at TC synapses (Gil et al., 1999; Yanagisawa et al., 2004). Significantly less depression was evoked in barrelless TC synapses at all ISIs tested (Fig. 1B,C) (p < 0.001 for difference between brl and ctrl across the range of ISIs; two-way ANOVA; p < 0.001 for difference between brl and ctrl at 50 and 75 ms ISIs and p < 0.05 for difference between brl and ctrl at 100 and 500 ms ISIs using a pairwise multiple comparison or Tukey’s test).
We also measured AMPA-PP ratios in control and barrelless TC synapses (Fig. 1D,E). As with NMDA-PP ratios, paired stimuli induced PPD in both control and barrelless neurons. AMPA-PP ratios were more variable in barrelless than in control neurons, probably because AMPA currents themselves are small in barrelless mutants (Lu et al., 2003), increasing their variability (pharmacologically reducing the amplitude of AMPA currents in control neurons leads to similar variability) (supplemental Fig. 1, available at www.jneurosci.org as supplemental material). Nonetheless, PP depression was significantly less in barrelless TC synapses compared with control across the range of ISIs (Fig. 1D,E) (p < 0.01 for difference between brl and ctrl across the range of ISIs; two-way ANOVA; no statistically significant difference for any one ISI using a pairwise multiple comparison or Tukey’s test).
We noted, as have others (Yanagisawa et al., 2004), that PP ratios for AMPAR-EPSCs were larger than for NMDAR-EPSCs at TC synapses. This suggests that postsynaptic mechanisms, such as NMDAR saturation (Chen et al., 2002), may contribute to short-term synaptic plasticity. To examine this, NMDA-PP ratios were measured in the presence of a low-affinity antagonist to NMDARs, l-AP-5 (800 μm). The substantial increase in NMDA-PP ratio with l-AP-5 demonstrates that NMDAR saturation normally contributes to short-term plasticity (Fig. 1F) (supplemental Fig. 2, available at www.jneurosci.org as supplemental material). However, even in the presence of l-AP-5, barrelless PP ratios remain significantly higher than in control neurons (p < 0.02), demonstrating that differences in NMDAR saturation are not causing the observed differences in short-term plasticity between barrelless and control neurons. There are many mechanisms that could contribute to the observed difference in PP ratios between barrelless and control TC synapses, including differences in the calcium-binding affinity of calcium sensors, calcium buffers, or the distance between calcium channels and neurotransmitter release machinery. We used high Ca2+ (4 mm) Ringer’s solution to maximize Pr at TC synapses and to saturate calcium sensors. Barrelless PP ratios measured in high Ca2+ Ringer’s solution are still significantly elevated relative to control neurons (p = 0.02), despite the slight reduction in overall PP ratio. The insignificant reduction in PP ratio in high Ca2+ Ringer’s solution in both control (p = 0.09) and barrelless neurons (p = 0.25) also suggests that Pr for TCAs is close to maximal in normal (2 mm) Ca2+ Ringer’s solution. This implies that the short-term plasticity deficits observed in barrelless TCAs are not attributable to intrinsic defects in neurotransmitter release machinery in the absence of AC1 function. Together, the increased PP ratios in barrelless, detected with both AMPA- and NMDA-EPSCs, suggests that short-term synaptic plasticity at barrelless TCAs is defective and the probability of neurotransmitter release (Pr) is reduced.
Slow block of NMDA-EPSCs by MK801 in barrelless TC synapses
The rate at which the amplitude of NMDAR-mediated EPSCs are attenuated by MK801 during repetitive stimulation is an indirect indication of presynaptic release probability (Hessler et al., 1993; Rosenmund et al., 1993; Manabe and Nicoll, 1994). This method depends on the fact that NMDA receptors must first be opened before they can be blocked by MK801, and the blockade is essentially irreversible. NMDA-EPSC amplitudes decrease progressively with repeated stimuli in the presence of MK801, and the rate of NMDA channel blockade varies inversely with the probability of vesicular release. NMDA-EPSCs decay quickly in the presence of MK801 in control neurons (Fig. 2A,C), indicating that Pr is high at TC synapses. Blockade by MK801 in barrelless neurons is much slower, and substantial residual response was detected even after 40 stimuli (Fig. 2B,C). Fitting the decay of the NMDA current in MK801 by an exponential reflects this difference, with a time constant of 19.37 ± 3.4 stimuli in barrelless neurons (n = 8) versus 4.64 ± 0.8 stimuli in control neurons (n = 14; p < 0.001; Mann–Whitney rank sum test). This difference in MK801 blockade is not a result of aberrant NMDAR function, because both the subunit composition of NMDARs (Lu et al., 2003) and the decay kinetics of NMDA currents (supplemental Fig. 3, available at www.jneurosci.org as supplemental material) are normal in barrelless synapses. The difference in the rate of NMDA-EPSC blockade by MK801 between control and barrelless neurons as well as the difference in paired-pulse ratio in barrelless neurons both suggest that Pr is lower in barrelless TC synapses.
Calcium-dependent release is different in barrelless TC synapses
Another measure of presynaptic function is the sensitivity of release to intracellular calcium concentration ([Ca2+]) (Wolfel and Schneggenburger, 2003). To test whether a change in Ca2+ sensitivity can account for, at least in part, the reduced Pr in barrelless TC synapses, we examined the response in P4–P6 control and barrelless neurons to a change in extracellular [Ca2+] in the range of 0.5–2 mm. At these relatively low concentrations, the intracellular [Ca2+] that drives vesicle fusion is likely to have a linear relationship with the extracellular [Ca2+] (Dodge and Rahamimoff, 1967; Reid et al., 1998). We measured NMDA-EPSC amplitudes as a function of extracellular [Ca2+] and quantitatively compared responses in control and barrelless neurons (Fig. 3). As expected, evoked EPSCs were significantly reduced when the external [Ca2+] was decreased (Fig. 3A,B). However, the relationship between EPSC amplitude and [Ca2+] was very different in control and barrelless synapses (p < 0.001 by ANOVA) (Fig. 3C), with barrelless neurons requiring higher [Ca2+] to reach comparable evoked amplitudes. Generally, the relationship between EPSC amplitude and extracellular [Ca2+] is different in low Pr synapses and high Pr terminals (Rozov et al., 2001), so the reduced sensitivity to [Ca2+] in barrelless TC synapses likely reflects reduced Pr. These data are consistent with our findings from the paired-pulse and MK801 experiments, all of which indicate functional deficits in the neurotransmitter release machinery in barrelless TCAs.
Synapsin and RIM proteins are putative AC1 targets in barrel cortex
Functional interactions between synaptic vesicle proteins and active zone proteins in the neurotransmitter release apparatus are modulated by neural activity through cAMP-dependent protein kinase A (PKA) phosphorylation (Lonart, 2002). A prominent PKA target in the presynaptic release apparatus is the RIM family of active zone proteins. RIM1α is directly phosphorylated by PKA, and the N-terminal PKA site (Ser-413) on RIM1α is required for cerebellar parallel fiber LTP (Lonart et al., 2003). We examined the degree of PKA phosphorylation of RIM using the N-terminal PKA site-specific antibody [PKA-phosphorylated RIM (P-RIM)] (Fig. 4A,B). Tissue homogenates were prepared from the somatosensory cortex of P4–P10 wild-type and barrelless mice. A decrease in the normalized P-RIM/RIM level was already evident in P4–P6 barrelless somatosensory cortex (p < 0.01) and was also obvious in P7–P10 barrelless mice (p < 0.01). There was no detectable difference in total RIM levels. We also examined the degree of PKA phosphorylation on rabphilin, another presynaptic PKA target (Lonart and Sudhof, 1998, 2001; Schluter et al., 1999). No difference was detected in the degree of PKA phosphorylation on rabphilin or in the total amount of rabphilin (Fig. 4A,C).
Synapsins are integral synaptic vesicle membrane proteins that are also major phosphoproteins in the neurotransmitter release machinery (Czernik et al., 1987; Greengard, 1987; Hilfiker et al., 1999; Hosaka et al., 1999; Jovanovic et al., 2001; Wang and Sudhof, 2003). Synapsins are phosphorylated by several kinases, including PKA and MAPK (mitogen-activated protein kinase). PKA phosphorylation on synapsins is implicated in the regulation of neurotransmitter release through synaptic vesicle mobilization (Chi et al., 2001). There are three known synapsin isoforms (Feng et al., 2002), and studies with synapsin I/II double knock-out (KO) mice suggest these proteins are required to maintain a normal number of synaptic vesicles and to regulate short-term synaptic plasticity (Rosahl et al., 1993, 1995; Li et al., 1995). Using a PKA site-specific antibody that recognizes PKA-phosphorylated synapsin (Jovanovic et al., 2001), we found that the fraction of PKA-phosphorylated synapsin relative to total synapsin is significantly reduced (p < 0.01) in barrelless somatosensory cortex (Fig. 4D,E). We also examined the quantity of MAPK-phosphorylated synapsin with two MAPK site-specific antibodies. No detectable difference was observed for the abundance of MAPK-phosphorylated synapsin or total synapsin between wild-type and barrelless synapses (Fig. 4D,F,G). In summary, we identified synapsins and RIMs, proteins involved in the regulation of neurotransmitter release, as putative AC1 targets in mouse barrel cortex.
Barrel map deficits in RIM1α KO mice
To clarify the role of RIMs and synapsins in barrel development, we examined the organization of barrel maps in RIM1α KO and synapsin I/II double KO mice with CO and Nissl histochemistry (Fig. 5), standard procedures for revealing whisker-related patterns at the subcortical and cortical levels of the rodent trigeminal pathway (Killackey and Leshin, 1975). In young (P8–P21) RIM1α KO mice, CO barrels corresponding to the small tactile hairs on the snout were absent or indistinct (Fig. 5B) (n = 10), whereas CO barrels for the large whiskers remained intact. Clear CO patterns for all of the barrels were observed in paired littermate controls (Fig. 5A) (n = 10). In adult RIM1α KO mice (3–6 months of age; n = 4), small CO barrel boundaries were again indistinct (Fig. 5D,F), in contrast to littermate controls (Fig. 5C) (n = 4) and synapsin I/II double KO mice (Fig. 5E) (n = 2). We also examined layer IV barrel cytoarchitecture with Nissl-stained serial sections from RIM1α KO mice and paired littermate controls. Distinct barrel walls throughout the barrel field were always evident in control mice (Fig. 5I) (n = 8). In contrast, a rather uniform distribution of layer IV neurons was observed in the somatosensory cortices of RIM1α KO mice (Fig. 5J) (n = 8), except for a few large barrels with a rudimentary barrel pattern. No obvious abnormality in the pattern of thalamic barreloids was found in RIM1α KO mice (Fig. 5G,H) (n = 3 pairs). A representative series of Nissl-stained sections through layer IV from a pair of control (Fig. 6A) and RIM1α KO mice (Fig. 6B) at P24 confirms the barrel cytoarchitectural deficits in RIM1α KO mice. These data strongly support the notion that RIM1α is an AC1 target in cortical barrel map formation.
In hippocampal neurons, RIM1α mutants show an increased PP ratio and slower block by MK801, reflecting functional deficits in presynaptic release efficacy (Schoch et al., 2002; Calakos et al., 2004). We therefore investigated whether RIM1α mediates, at least in part, AC1 modulation of presynaptic release efficacy at TC synapses during barrel map development. Indeed, PP ratios for both NMDAR-EPSCs and AMPAR-EPSCs were significantly higher in RIM1α KO TC synapses than littermate control cells for 50, 75, and 100 ms ISIs (Fig. 6C,D). The degree of depression triggered by 500 ms ISIs in RIM1α KO synapses was not different from control neurons for both NMDA- and AMPA-EPSCs. This is similar to what was observed in hippocampal synapses (Schoch et al., 2002; Calakos et al., 2004) but different from what was observed at barrelless TC synapses. This suggests that PPD triggered by long ISIs at TC synapses is mediated through an AC1-dependent pathway that is independent of RIM1α. The barrel cytoarchitecture and TCA neurotransmitter release deficits in RIM1α mutants together strongly suggest that efficient presynaptic function plays an important role in cortical barrel map development.
Decreased presynaptic release probability can disrupt barrel segregation
To demonstrate how the observed deficits in synaptic function in barrelless mice could disrupt barrel development, we constructed a model for neural activity-dependent (Hebbian) barrel map formation and tested the effects of reducing neurotransmitter release efficacy on barrel segregation (Fig. 7) (see Materials and Methods, Model). The activity-dependent development of neural circuits is thought to be driven by synaptic learning rules that translate patterns of neural activity in presynaptic and postsynaptic neurons into changes in synaptic strength (Crair and Malenka, 1995; Katz and Shatz, 1996). We used a spike time-dependent learning rule (Fig. 7B) (Abbott and Nelson, 2000; Dan and Poo, 2004) similar to that known to drive synaptic development in layer IV to layer II/III synapses in somatosensory cortex (Feldman, 2000) and retinotopic map development at the retinotectal synapse (Zhang et al., 1998).
With a spike timing-dependent learning rule (Fig. 7B), strengthening of the synapse occurs when an EPSC precedes a cortical action potential (spike), and weakening occurs when an EPSC follows the cortical spike. Because thalamic afferents burst while cortical neurons typically respond with isolated spikes (Simons and Carvell, 1989; Pinto et al., 2003), afferent spikes can occur before or after an evoked cortical spike, in proportions dictated by the probability of release (Fig. 7B). For high probability of release synapses, afferent spikes will typically lead cortical spikes, resulting in more LTP than LTD. Because the initially strong PW synapses evoke more cortical spikes than the AW, they are further strengthened. The stronger PW synapses also cause more cortical spikes to be uncorrelated with AW activity, resulting in more AW activity that falls into the larger LTD window. Thus, for high release probability (low PPR) (Fig. 7C), segregation occurs because of “potentiation” from afferent activity that drives cortical spikes, and “depression” from afferent activity uncorrelated with cortical spiking: PW synapses strengthen because there is more LTP than LTD, and AW afferents weaken because the ratio balances in favor of LTD.
This LTP/LTD balance is achieved only if afferent whisker stimulation, which produces bursts of thalamic activity in vivo, does not induce bursts of cortical activity; otherwise, LTP will always dominate because afferent spikes will precede most cortical spikes. It is therefore essential to include an inhibitory component in the model, which dampens cortical response to whisker stimulation to one or two spikes, consistent with observations of layer IV in vivo activity at the earliest ages studied (Foeller and Feldman, 2004). This inhibition may occur through feedforward GABAergic circuits (Foeller and Feldman, 2004), or via other developmentally regulated mechanisms, such as modulation of neurotransmitter release by monoamines (Rhoades et al., 1994).
Reducing the probability of vesicular release can also upset the balance between LTP and LTD, because afferent activity more often fails to elicit vesicular release before the postsynaptic spike (Fig. 7B) (PPR, 0.7 instead of 0.5). As a result, the afferent bursting that evokes a cortical spike will result in less strengthening, whereas the amount of weakening from uncorrelated afferent activity stays constant. This causes the weakening of both PW and AW afferents until they can no longer drive cortical spikes (Fig. 7D,E) (PP ratios >0.6). The failure to segregate for high PP ratios predicts the final state of TC synapses observed in barrelless, in which afferents do not segregate into barrels, and AMPAR-mediated currents are small. Thus, presynaptic neurotransmitter release deficits will favor LTD over LTP, disrupt segregation, and lead to weaker TC synapses.
The presynaptic deficits described here are consistent with a more general situation in which the strength of LTP is decreased relative to LTD. Directly changing the ratio of LTP relative to LTD in the model also prevents segregation (Fig. 7F). Our model does not specifically distinguish between the role of the presynaptic deficit described here in barrelless mice and the previously observed postsynaptic deficits in AMPAR trafficking (Lu et al., 2003), which also might affect the overall strength of LTP relative to LTD. Rather, it illustrates how interfering with either may disrupt thalamic afferent segregation and barrel formation.
Discussion
We presented evidence that barrelless (AC1) mutant mice have impaired neurotransmitter release machinery that contributes to their barrelless phenotype. With biochemical techniques, we identified RIM proteins as likely targets for AC1/PKA signaling. We further showed that RIM1α mutant mice have TCA neurotransmitter release deficits and impaired barrel patterning, which confirms that interfering with presynaptic release negatively impacts neural circuit development. Finally, we used a computational model to demonstrate how inadequacies in presynaptic function will obstruct the activity-dependent process responsible for the segregation of sensory afferents. Together, our data suggest that efficient neurotransmitter release at thalamocortical synapses, which is a general feature of afferents in all sensory pathways, is important in the development of precise cortical maps.
AC1 regulates presynaptic release in TC synapses through RIMs
We found a significant reduction in the PKA phosphorylation level of RIMs in barrelless somatosensory cortex, suggesting that AC1 regulates presynaptic vesicular release through the phosphorylation of these proteins. Like barrelless mutants, we also found that RIM1α mutant TC synapses have decreased paired-pulse depression, and others have noted in the hippocampus of RIM1α mutants similar deficits in short-term plasticity and differences in the rate of synaptic blockade by MK801 (Schoch et al., 2002; Calakos et al., 2004). The barrel pattern defects in RIM1α knock-out mice also strongly suggest that the AC1-mediated regulation of presynaptic function responsible for barrel development acts at least in part through RIMs. However, there are a number of ways independent of RIM proteins in which AC1 may also modulate presynaptic release. For example, PKA phosphorylation of synapsins, which is also reduced in barrelless mice, may modulate neurotransmitter release through synaptic vesicle mobilization (Chi et al., 2001), although synapsin mutants have normal barrels. Other direct PKA targets that may be responsible for the decrease in presynaptic release efficacy observed in barrelless mutants include cysteine string protein (Evans et al., 2001), SNAP25 (synaptosome-associated protein 25) (Risinger and Bennett, 1999; Hepp et al., 2002; Nagy et al., 2004), snapin (Chheda et al., 2001), and syntaphilin (Lao et al., 2000; Boczan et al., 2002). The changes in release sensitivity to [Ca2+] we observed in barrelless mutants may also be mediated by a downstream target other then RIM, because [Ca2+] sensitivity is unchanged in RIM1α mutant mice (Calakos et al., 2004). It is also possible that AC1 modulates neurotransmitter release independent of PKA, because RIM1α binds to cAMP-GEFII (guanine nucleotide exchange factor II) (Ozaki et al., 2000), a direct target of cAMP that is responsible for cAMP-dependent, but PKA-independent exocytosis. A detailed understanding of the signaling cascades triggered by AC1 that regulate neurotransmitter release will require alternative experimental approaches, possibly autaptic neuronal cultures or some other reduced preparation with greater power to assay the molecular basis for defects in neurotransmitter release. Regardless, our data clearly demonstrate a presynaptic role for AC1 in regulating the efficiency of vesicular release, and show that the AC1 phosphorylation target RIM1α contributes to barrel map development.
Barrel map deficits in RIM1α mutants
There are a handful of mutants that display barrel map deficits (Erzurumlu and Kind, 2001; Hannan et al., 2001; Lopez-Bendito and Molnar, 2003). Barrelless (AC1 KO) mice rank among the strongest phenotypes; with no hint of barrel patterning in either the thalamocortical afferent arbors or the postsynaptic neuronal cytoarchitecture. It is more common, as in RIM1α knock-out mice, for the cortical neuronal cytoarchitecture to be absent or dramatically disturbed while thalamocortical afferent patterning is mostly intact (Cx-NR1 and PLCβ1 mutants are similar). The small anterior barrels are also typically more disturbed than the larger barrels in the posterior medial barrel subfield (Erzurumlu and Kind, 2001; Hannan et al., 2001). The barrel deficits observed in RIM1α knock-outs are therefore not as severe as in some mutants but are also not atypical. Even perturbations that are thought to specifically affect thalamic afferent terminals and not cortical neurons cause more dramatic effects on the small barrels and postsynaptic barrel cytoarchitecture than presynaptic thalamocortical afferent clustering of the large posterior medial barrels (Salichon et al., 2001). Nonetheless, the greater effect of mutating RIM1α on postsynaptic cortical barrel cytoarchitecture than on presynaptic TCA clustering could be attributable to a purely cortical role of RIM1α on barrel development in which presynaptic function between cortical neurons is disturbed in RIM1α mutants. The more profound defects in barrel map development and neurotransmitter release observed in barrelless compared with RIM1α mutants suggests that AC1 signaling is mediated through presynaptic pathways beyond RIM1α, perhaps other RIMs, or combines with established postsynaptic pathways mediated by AC1 (Lu et al., 2003) to produce a more dramatic barrel phenotype.
Cortical role in barrel map formation
Some have argued that factors intrinsic to the cortex instruct sensory map formation (Fukuchi-Shimogori and Grove, 2001), whereas others contend that barrel map features depend on, at least in part, signals communicated by thalamic axons (Lopez-Bendito and Molnar, 2003). The segregation of TCAs in Cx-NR1 KO mice is fuzzy, with TCA arbors that are more diffuse than normal (Iwasato et al., 2000; Lee et al., 2005). This observation suggests that cortical NR1 is required for the maintenance and refinement of TCA arbors, and that communication between TCAs and layer IV neurons is crucial for the development of a barrel pattern. This situation is likely to be the case for barrelless as well. Previously, we showed that NMDAR-dependent LTP at thalamocortical synapses requires postsynaptic PKA activity and is difficult to induce in barrelless mice, probably because of an inability to properly regulate AMPAR trafficking (Lu et al., 2003). This LTP defect, as is captured in our model (Fig. 7F), can mask activity-dependent signals responsible for the regulated addition or elimination of particular synaptic targets on the postsynaptic neuron. In a mutation that only affects postsynaptic neurons, one might expect that the asymmetric orienting of dendrites from layer IV cortical neurons into barrel hollows is more severely affected than the clustering of presynaptic afferents into a barrel pattern. For instance, Cx-NR1 mutant mice completely lack postsynaptic layer IV barrel cytoarchitecture, but TCAs still aggregate to form a rudimentary barrel pattern (Iwasato et al., 2000). The complete absence of both presynaptic and postsynaptic barrel patterning in barrelless (AC1) mutants suggests that presynaptic AC1 function plays an important role in the barrelless phenotype. This issue will not be fully resolved, however, until thalamus-specific or cortex-specific AC1 mutants are developed and tested.
How do afferents instruct their postsynaptic partners to form barrels?
TCAs that innervate the barrel field initially (approximately P2–P3) have relatively diffuse projections, and the progressive addition of branches and synapses at the cortical location corresponding to the mapped facial whisker and the commensurate elimination of inappropriate branches (Senft and Woolsey, 1991; Agmon et al., 1995; Rebsam et al., 2002) in an NMDA receptor-dependent manner (Lee et al., 2005) generate an afferent barrel pattern. The extensive tangential coverage of arbors from single barrelless thalamocortical afferents (Welker et al., 1996) is similar to the enlarged retinal ganglion cell (RGC) axon projection in the superior colliculus of barrelless mutants (Ravary et al., 2003; Plas et al., 2004), and RGC arbors after retinal activity blockade in a variety of species (O’Leary et al., 1986; Sretavan et al., 1988). The pattern of RGC activity is crucial, not just its presence, because enlarged projections are found when correlated retinal waves are eliminated but uncorrelated spontaneous RGC activity persists (McLaughlin et al., 2003; Chandrasekaran et al., 2005). In the frog optic tectum, visual activity governs the development of axonal projections by altering the retraction dynamics of axon terminals (Ruthazer et al., 2003). The establishment and maintenance of sensory maps are also perturbed when the excitability of sensory neurons is diminished (Yu et al., 2004; Hua et al., 2005), and neurons that cannot release neurotransmitter fail to recruit postsynaptic AMPARs (Harms et al., 2005) or maintain their central connections. The barrelless phenotype is likely attributable to similar barriers in this activity-dependent process.
How do deficits in presynaptic function interfere with barrel map development? Our model demonstrates that interfering with neurotransmitter release disrupts the balance between spike timing-dependent LTP and LTD and prevents thalamocortical afferents with correlated postsynaptic activity from strengthening, which is necessary for segregation to occur (Fig. 7). Thus, sufficient disruption of presynaptic function will interfere with Hebbian circuit development, as will direct intervention with postsynaptic plasticity mechanisms. Interestingly, neurons that cannot release neurotransmitter, in a competitive scenario, are able to form synapses but fail to recruit postsynaptic AMPARs (Harms et al., 2005). It is a common feature of many sensory pathways for thalamocortical afferents to have relatively high release probability and potent short-term depression (Stratford et al., 1996; Gil et al., 1999; Amitai, 2001; Iwasaki and Takahashi, 2001; Chung et al., 2002; Murphy et al., 2004). We propose that high release probability at TC synapses permits postsynaptic partners to use a correlation-based mechanism to mediate circuit development. Presumably, the activity-dependent development of other sensory maps, including ocular dominance columns and retinotopic maps, relies on afferents having high release probability and are similarly sensitive to perturbations that interfere with presynaptic function.
Footnotes
H.-C.L. was supported by National Institute of Neurological Disorders Grant NS48884, American Heart Association Grant 0465158Y, and the Gordon and Mary Cain Foundation. D.A.B. was supported by a Charles A. King Trust postdoctoral fellowship (Bank of America CoTrustee, Boston, MA). R.J. was supported by the W. M. Keck Foundation. M.C.C. was supported by National Institute of Mental Health Grant MH62639, the Merck and Klingenstein Foundations, and Mental Retardation Research Center at Baylor College of Medicine Grant HD24064. We thank Drs. R. Heidelberger, C. Rosenmund, and T. C. Südhof and members of the Crair laboratory for comments and discussion on this manuscript. Anti-phosphorylated RIMs and synapsin antibodies were generous gifts provided by Drs. G. Lonart and P. Greengaurd. The Developmental Studies Hybridoma Bank (University of Iowa) provided the anti-SV2 antibody.
- Correspondence should be addressed to either of the following: Hui-Chen Lu, The Cain Foundation Laboratories, Department of Pediatrics, Division of Neurology/Developmental Neuroscience, Baylor College of Medicine, 6621 Fannin Street, MC3-6365, Houston, TX 77030, Email: hclu{at}bcm.tmc.edu or Michael C. Crair, Department of Neuroscience, Baylor College of Medicine, One Baylor Plaza, S-603, Houston, TX 77030, mcrair{at}bcm.tmc.edu