Reliably detecting unexpected sounds is important for environmental awareness and survival. By selectively reducing responses to frequently, but not rarely, occurring sounds, auditory cortical neurons are thought to enhance the brain's ability to detect unexpected events through stimulus-specific adaptation (SSA). The majority of neurons in the primary auditory cortex exhibit SSA, yet little is known about the underlying cortical circuits. We found that two types of cortical interneurons differentially amplify SSA in putative excitatory neurons. Parvalbumin-positive interneurons (PVs) amplify SSA by providing non-specific inhibition: optogenetic suppression of PVs led to an equal increase in responses to frequent and rare tones. In contrast, somatostatin-positive interneurons (SOMs) selectively reduce excitatory responses to frequent tones: suppression of SOMs led to an increase in responses to frequent, but not to rare tones. A mutually coupled excitatory-inhibitory network model accounts for distinct mechanisms by which cortical inhibitory neurons enhance the brain's sensitivity to unexpected sounds.https://doi.org/10.7554/eLife.09868.001
In everyday life, we are often exposed to a mix of different sounds. An essential task for our brain is to separate the important sounds from the unimportant ones. For example, stepping out onto a busy street, you may at first be very aware of the noise of traffic. Later, you may start to ignore the din and instead only notice sounds that break the monotony: a honking car horn or maybe a stranger's voice. This is because the neurons in the auditory pathway respond differently to common and rare sounds. In particular, excitatory neurons in the region termed the ‘auditory cortex’ send fewer nerve impulses in response to frequent sounds, but respond vigorously to rare sounds. This phenomenon is called ‘stimulus-specific adaptation’, but it is not known exactly which neurons in this brain region enable this process to occur.
Now, Natan et al. have combined different cutting-edge neuroscience techniques to identify the circuit of brain cells that drives this stimulus specific adaptation. A technique called optogenetics was used to effectively ‘turn off’ each of two kinds of inhibitory neuron in the auditory cortex of mice, by exposing the brain to colored light from a laser.
Natan et al. found that both kinds of inhibitory neuron amplified stimulus-specific adaptation, but via different mechanisms. One of these neuron types, called ‘parvalbumin-positive interneurons’, exerted a general effect on excitatory neurons and suppressed responses to both frequent and rare sounds As the responses to rare sounds started off greater than the responses to frequent sounds, suppressing both by an equal amount actually led to an increase in the relative difference between them. On the other hand, the second kind of inhibitory neuron, called ‘somatostatin-positive interneurons’, only reduced the excitatory neurons' responses to frequent sounds; these neurons had no effect on responses to rare noises.
Future studies will test how specific adaptation in different contexts can help us to behaviorally detect rare sounds while ignoring common ones, and search for the circuits beyond the auditory cortex that support hearing in complex sound environments.https://doi.org/10.7554/eLife.09868.002
Across sensory modalities, cortical neurons exhibit adaptation, attenuating their responses to redundant stimuli (Das and Gilbert, 1999; Ulanovsky et al., 2003; Garcia-Lazaro et al., 2007; Asari and Zador, 2009; Khatri et al., 2009). Adaptation to stimulus context is thought to increase efficiency of sensory coding under the constraints of limited resources (Barlow, 1961). Yet, the neuronal-circuit mechanisms that facilitate adaptation in the cortex remain poorly understood. In the primary auditory cortex (A1), the vast majority of neurons exhibit stimulus-specific adaptation (SSA, Figure 1). When presented with a sequence of two tones, one of which occurs frequently (termed ‘standard’) and another rarely (termed ‘deviant’), the neuron's response to the standard tone becomes weaker, but the response to the deviant tone remains strong (Ulanovsky et al., 2003; Szymanski et al., 2009; Farley et al., 2010; Fishman and Steinschneider, 2012). Whereas SSA has also been found in sub-cortical structures e.g. in the auditory midbrain (Malmierca et al., 2009; Zhao et al., 2011; Thomas et al., 2012) and the auditory thalamus (Kraus et al., 1994; Anderson et al., 2009; Antunes et al., 2010; Bauerle et al., 2011), it is weak in the lemniscal areas of the auditory pathway, which project to A1, and stronger in those non-lemniscal areas which receive feedback from A1 (Ulanovsky et al., 2004, Perez-Gonzalez et al., 2005, Duque et al., 2012). Therefore, cortical circuits are proposed to contribute to and amplify SSA in A1 (Ulanovsky et al., 2003, Szymanski et al., 2009, Bauerle et al., 2011, Fishman and Steinschneider, 2012, Escera and Malmierca, 2014), through a combination of plastic modulation of thalamocortical inputs and intra-cortical inhibitory circuits, which would allow for selective suppression of neuronal responses to specific stimuli (Nelken, 2014). Our study tests whether and how inhibitory neurons contribute to cortical SSA.
Auditory cortex, like other sensory cortices, contains morphologically and physiologically diverse inhibitory interneurons, which form dense interconnected networks with excitatory neurons (DeFelipe, 2002; Douglas and Martin, 2004). While different interneuron types have been hypothesized to carry out specialized complementary functions in sensory processing (DeFelipe, 2002; Markram et al., 2004; Isaacson and Scanziani, 2011; Kepecs and Fishell, 2014; Marlin et al., 2015), their function in driving changes in dynamic auditory processing has not been previously established. We hypothesized that the two most common types of interneurons in the cortex, parvalbumin- (PVs) and somatostatin-positive interneurons (SOMs) (Xu et al., 2010; Rudy et al., 2011), facilitate SSA in excitatory neurons of A1 in a complementary fashion. PVs, a subset of which receive direct thalamic inputs (Staiger et al., 1996), may amplify SSA in excitatory neurons by providing a constant inhibitory drive; equally strong inhibitory drive would attenuate the weak response to standard tones relatively more than the strong response to deviant tones, leading to a greater differential between standard vs deviant tone spiking response. SOMs, which target distal dendrites of pyramidal cells (McGarry et al., 2010; Gentet et al., 2012), have excitatory synapses that exhibit facilitation upon repetitive stimulation (Reyes et al., 1998; Silberberg and Markram, 2007). Therefore, inputs from SOMs may exert a stimulus-specific increase in suppression of excitatory neurons that is selective to the standard tone and does not generalize to the deviant tone. Alternatively, they may contribute to selective adaptation in excitatory neurons through differential post-synaptic integration.
To tease apart the function of different inhibitory types in SSA, we tested whether optogenetic suppression of either PV or SOM interneurons during sound presentation reduced SSA in putative excitatory neurons in the auditory cortex (Hamilton et al., 2013; Pi et al., 2013; Weible et al., 2014). We found that both types of interneurons contribute to SSA in the cortex, with PVs providing constant inhibition, and SOMs increasing their effect with repeated tones.
We recorded spiking activity of neurons as well as local field potentials (LFPs) in A1 in head-fixed mice under light isoflurane anesthesia. SSA was measured from the firing rate (FR) of neurons in response to tones presented as a series of ‘oddball’ stimuli. Each oddball stimulus consisted of a sequence of tone pips at one of two frequencies (tones A and B). In each oddball stimulus, one tone was presented as the rare (deviant) tone, while the other was presented as the frequent (standard) tone (A to B ratio of 90:10 or 10:90, Figure 1A). A third stimulus was also presented (equal stimulus), with tones A and B being presented equally often (50:50). The frequencies of tone A and B were selected at 0.39 octave intervals, narrower than the typical tuning bandwidth of A1 neurons (Hackett et al., 2011; Guo et al., 2012; Kanold et al., 2014), such that they activated the majority of recorded neurons on each session (Figure 1B).
As expected, for a representative neuron recorded in A1, the mean FR in response to a tone was lower when the tone was presented as the standard than as the deviant (Figure 1C), exhibiting SSA. To quantify the level of adaptation for each neuron, we computed the index of the change in FR to the same tone when it was presented as the deviant vs the standard (SSA index). SSA index is 1 when adaptation is complete (i.e., no response to the standard, and significant response to the deviant), and 0 when there is no adaptation (i.e., the response to the standard and deviant is equal). Almost all neurons recorded in A1 exhibited significant SSA (Figure 1D, standard tone-evoked FR significantly lower than the deviant tone-evoked FR in N = 138 out of 147 neurons, Wilcoxon rank sum test p < 0.05).
We first tested whether SSA is present in inputs from the thalamus. Current source density (CSD) analysis has been extensively used to quantify inputs from the thalamus (Metherate and Cruikshank, 1999; Kaur et al., 2005; Szymanski et al., 2009; Happel et al., 2014). We used a linear probe to record LFPs using electrodes spaced 50 microns apart inserted perpendicularly to brain surface in the primary auditory cortex. The multi-electrode probe is 775-µm long, spanning layers 1–6 of mouse A1. CSD is computed as the second spatial derivative of the LFPs across the depth of the cortex (Figure 1E, Figure 1—figure supplement 1A, 20 sessions, 15 mice). Typically, in response to tones, CSD exhibits a negative basin, termed sink, within a short delay of tone onset, localized to electrodes in thalamo-recipient layer (Figure 1F, Figure 1—figure supplement 1B) (Kaur et al., 2005; Szymanski et al., 2009). The amplitude of current in the sink was taken as a measure of the combined strength of post-synaptic inputs onto layer 4 neurons, which should reflect the strength of the thalamic inputs to the cortex (Metherate and Cruikshank, 1999; Kaur et al., 2005; Szymanski et al., 2009; Happel et al., 2014).
We compared the amplitude of the CSD sink for each tone when presented as a deviant or standard, and computed their ratio (Figure 1F). The sink amplitude was lower for the standard as compared to the deviant tones (Figure 1F,G), suggesting that excitatory signals produced by thalamo-cortical inputs exhibit SSA, consistent with previous findings (Szymanski et al., 2009). This finding supports the ‘adaptation in narrowly tuned inputs’ model, which postulates that SSA in broadly tuned neurons in A1 reflects adaptation in either thalamocortical inputs or at the stage of integration of thalamocortical inputs, specific to inputs tuned to the standard tone (Mill et al., 2011; Taaseh et al., 2011; Nelken, 2014). Importantly, across sessions, the SSA index of the granular layer CSD sinks was significantly lower than that of either the non-thalamo-recipient layers (Δ = −28%, p-value from one-sided test after correction (p1) = 6e−4, z = −3.4, Bonferroni corrected for two tests (C = 2)) or the SSA index of the mean spiking activity of A1 neurons (Δ = 23%, p1 = 0.029, z = −2.1, C = 2) in each session (N = 20 sessions in 15 mice, Figure 1G), suggesting that additional intra-cortical mechanisms may contribute to SSA in the cortex.
We next tested whether cortical inhibitory interneurons may contribute to SSA. Since different inhibitory neuronal subtypes can differentially affect sensory responses of putative excitatory neurons (Lee et al., 2012; Wilson et al., 2012; Cottam et al., 2013), we separately tested the role of PVs and SOMs. We used targeted viral delivery in the auditory cortex of mice to drive Archaerhodopsin (Arch) expression, which hyperpolarizes neurons when stimulated by light, in either PVs or SOMs (Chow et al., 2010). A modified adeno-associated virus (AAV) encoding anti-sense code for Arch and a fluorescent reporter, under the FLEX cassette, was injected into PV-Cre or SOM-Cre mice (Boyden et al., 2005; Sohal et al., 2009; Cardin et al., 2010; Zhang et al., 2010; Deisseroth, 2011) (Figure 2A). 2 weeks following virus injection, Arch was expressed selectively in PVs or SOMs in auditory cortex at expected levels (Kvitsiani et al., 2013) (Figure 2B,C PV-Cre: N = 250 neurons in 4 mice, specificity = 92 ± 1%, efficiency = 73 ± 5%. SOM-Cre: N = 149 neurons in 5 mice, specificity = 95 ± 2%, efficiency = 86 ± 5%). To activate Arch, a light guide was positioned to cast 180 mW/mm 532-nm light onto A1 surface, perpendiular to cortical layers. In vitro intracellular recordings from optically identified PVs or SOMs (Figure 2—figure supplement 1, Figure 2—figure supplement 2) demonstrate that light cast over the auditory cortex in vitro drives a strong suppressive current (Figure 2D, Figure 2—figure supplements 1C,D, 2C,D) and hyperpolarizes the membrane potential in these neurons (Figure 2—figure supplements 1B, 2B). Assuming a 100-fold attenuation of light over 1 mm of brain tissue (Aravanis et al., 2007), the estimated irradiance in the deepest cortical layer (1.8 mW/mm2) was strong enough to induce hyperpolarizing current in neurons in vitro (Figure 2D). In vivo, in both PV-Cre and SOM-Cre mice, illuminating the auditory cortex suppressed spiking activity in a small subset of recorded neurons (Figure 2E,F, left, putative inhibitory neurons) and increased activity in a great majority of recorded neurons (Figure 2E,F right, putative excitatory neurons). Shining light over A1 increased spontaneous neuronal activity in the majority of the recorded neurons in both PV-Cre mice (N = 115 neurons, 102 increased, 0 decreased, in 10 mice) (Figure 2G) and SOM-Cre mice (N = 104 neurons, 61 increased, 3 decreased, in 9 mice) (Figure 2H). These measurements demonstrate that casting light over A1 selectively and effectively suppresses the activity of either PVs or SOMs.
To test the function of PVs and SOMs in SSA, their activity was suppressed during every fifth tone of the oddball stimulus by illuminating A1 (Figure 3A). To directly test the effect of interneuron suppression, we computed the SSA index separately on light-on and light-off trials for neurons responsive to both tones A and B (SSA was found in 63 out of 67 tone-responsive neurons in PV-Cre mice, 42 out of 43 tone-responsive neurons in SOM-Cre mice). Photosuppression of either PVs or SOMs affected the responses of neurons to the tones (Figure 3B,C), resulting in a significant reduction in SSA index across the population (Figure 3E,F, PV-Cre: Δ = −41%, p1 = 1e−12, t(66) = 8.6. SOM-Cre: Δ = −25%, p1 = 2e−6, t(42) = 5.4). Photo-manipulation-affected responses only to the tone during which it was presented, but not to subsequent tones (Figure 3—figure supplement 1). Additionally, photo-manipulation was limited to cortex since it did not affect thalamo-recipient layer CSD tone responses and SSA (Figure 3—figure supplement 2). In a control group of PV-Cre or SOM-Cre mice (6 mice), we injected a modified AAV, which encoded anti-sense fluorescent reporter alone under the FLEX cassette, and computed the effect of casting light on SSA (SSA was found in 33 out of 37 tone-responsive neurons in control mice). In this control group, SSA was not affected by light (Figure 3D,G, p > 0.05, t(36) = −2.0), confirming that Arch expression was required for the effect of the light. Therefore, the effects of interneurons are specific to intra-cortical mechanisms. These results demonstrate that both types of interneurons contribute to the reduction of the response of the neuron to the stimulus during SSA.
A decrease in the SSA index may be due to several factors: (1) an increase in response to the standard only, (2) a decrease in response to the deviant, or (3) an increase in response both to the standard and the deviant, but with a relatively greater increase for the standard. Therefore, we next investigated the effect of interneuron photosuppression on FR of putative excitatory neurons evoked by the standard and deviant tones separately. The effects of PVs and SOMs diverged; in addition to increasing spontaneous activity (Δ = 185%, p-value from one-sided t-test after correction (p2) = 3e−11, t(159) = −7.2), suppressing PVs led to increased FR to both the standard (Δ = 102%, p2 = 3e−11, t(159) = −7.2) and deviant (Δ = 56%, p2 = 9e−12, t(159) = −7.4) tones (N = 160, Figure 4A–C, Figure 4—figure supplement 1A). 83% of neurons exhibited greater FR to the standard and 46% to the deviant during PV suppression (Figure 4—figure supplement 1B). The difference in FR due to suppression of PVs was not significantly different between the standard and deviant tones (p2 > 0.05, t(159) = −0.1, C = 2) but both were greater than the difference in the spontaneous FR (standard: Δ = 25%, p2 = 0.001, t(159) = −3.6, C = 2. Deviant: Δ = 26%, p2 = 0.039, t(159) = −2.4, C = 2), indicating that the change in tone-evoked FR was similar regardless of tone probability (Figure 4B, bottom panel). Because an equal increase in the FR produces a weaker relative effect on the response to the deviant (which is higher than to the standard), PV inactivation decreases SSA index (Figure 3E).
By contrast, suppressing SOMs led to an increase in FR for spontaneous activity (Δ = 46%, p2 = 2e−9, t(113) = −6.5) and during the standard (Δ = 29%, p2 = 2e−8, t(113) = −6.1) but not deviant (p2 > 0.05, t(113) = −0.8) tone (N = 114, Figure 4D–F, Figure 4—figure supplement 1C). 52% of neurons exhibited greater FR to the standard and only 11% to the deviant during PV suppression (Figure 4—figure supplement 1D). The increase in FR for spontaneous activity was not different than that during the standard tone (p2 > 0.05, t(113) = 0.2, C = 2) and the differences in FR due to suppression of SOMs were stronger for spontaneous activity and the standard tone than the deviant tone (spontaneous: Δ = 390%, p2 = 0.004, t(113) = 3.1. Standard: Δ = 378%, p2 = 0.005, t(113) = 3.1) (Figure 4E, bottom panel), thereby accounting for the change in SSA with SOM inactivation (Figure 3F). Responses to the equal stimulus evoked consistent, yet weaker effects (Figure 4—figure supplement 2).
PVs and SOMs differ in their density among different layers of the cortex and in laminar sources and targets of their inputs and outputs (Markram et al., 2004; Xu and Callaway, 2009; Fino et al., 2013). The effects of PV and SOM suppression on SSA had differential laminar distribution (Figure 4—figure supplement 3). The effect of PVs on SSA was equally strong in the supra-granular and infra-granular layers, but stronger in the granular layer, that is, the thalamo-recipient layer. This differential effect is consistent with the relative proportion of cortical interneurons that are PVs, which is higher in granular than either in infra- or supra-granular layers (Lee et al., 2010; Xu et al., 2010; Ouellet and de Villers-Sidani, 2014). In contrast, suppressing SOMs reduced SSA in the granular and infra-granular, but not supra-granular layers. The relative proportion of cortical interneurons that are SOMs is greatest in the granular and infra-granular layers, but still present in supra-granular layers (Lee et al., 2010; Xu et al., 2010; Ouellet and de Villers-Sidani, 2014). As some SOMs predominantly target the distal dendrites of pyramidal neurons (Markram et al., 2004), the effect of suppressing SOMs in supra-granular layers may be evident in recordings of pyramidal neurons with cell bodies in deeper layers, supporting our results. In addition, cortical extracellular recordings may be biased toward neurons in granular and infra-granular layers, precluding adequate sampling of activity in superficial layers. In controls, we did not observe a difference in the effect of light on SSA across layers, demonstrating that the differences are not due to differential artifact of light stimulation.
Our results indicate that both PVs and SOMs affect SSA, but in different ways: (1) the increase in the FR of putative excitatory neurons due to PV suppression is constant, either during presentation of the standard or the deviant, and greater than changes in spontaneous activity. Thus, PVs amplify SSA in excitatory neurons by exerting a relatively stronger inhibitory drive for the standard than for the deviant. (2) Suppression of SOMs leads to increased putative excitatory neuron activity only during the spontaneous firing or the presentation of the standard, but not for the deviant. This suggests that the strength of SOM-mediated inhibitory drive is not significant in response to the deviant but increases with repeated presentations of the standard.
In neurons exhibiting SSA, responses to the deviant are stronger than to the standard. This difference might lead to a ‘ceiling’ effect, reducing the effect of PV photosuppression on FR to the deviant, but not standard (Olsen et al., 2012). However, restricting the analysis to two subpopulations of neurons, which have matched mean and standard deviation of FR to the standard vs the deviant tones (Ulanovsky et al., 2004; Rust and Dicarlo, 2010), preserved the observed effects of photosuppression (Figure 4—figure supplement 4). Suppressing PVs led to an equal increase in FR to both the standard and the deviant tone (N = 54—standard: Δ = 62%, p2 = 6e−8, t(53) = 6.3. Deviant: Δ = 55%, p2 = 3e−5, t(53) = 4.5. Standard vs deviant: p2 > 0.05, t(53) = 0.5). In contrast, suppressing SOMs led to a significant increase in FR to the standard, but no change in FR to the deviant (N = 44—standard: Δ = 30%, p2 = 7e−6, t(43) = 5.1. Deviant: p2 > 0.19, t(43) = 1.3. Standard vs deviant: Δ = 382%, p2 = 6e−4, t(43) = 3.7).
For neurons that responded more strongly to one of the tones (‘strong’ vs ‘weak’ tone), a ceiling effect would predict that the effect of interneuron suppression would be stronger for the weak than the strong tone. However, PV and SOM suppression exhibited a similar effect on responses to the strong and the weak tones in neurons that exhibited differential responses to two tones (Figure 4—figure supplements 5, 6). Suppressing PVs led to similar increases in tone-evoked FR between weak and strong tones for both deviant (N = 51, p2 > 0.05, t(50) = 1.0) and standard tones (p2 > 0.05, t(50) = −1.9). Suppressing SOMs also led to similar differential effects between strong and weak tones; standard tone-evoked FR increased equally (N = 34, p2 > 0.05, t(33) = 1.1) and deviant tone-evoked FR was equally unchanged (p2 = 0.05, t(33) = −0.1). Combined, these analyses demonstrate that the effect of PV photosuppression on SSA cannot be explained by the ceiling effect for either PVs or SOMs.
Although Arch drove strong currents in both SOM and PV neurons (Figure 2D, Figure 2—figure supplements 1, 2), there might be a difference in expression level or efficacy of Arch between SOM-Cre and PV-Cre mice, leading to a stronger effect of photosuppression in PV-Cre than in SOM-Cre mice on tone-evoked FRs (Figure 4B,E). Alternatively, the difference might be attributable to the morphological or functional differences between SOMs and PVs. To address this confound, we selected tone responses that exhibited matched difference in standard tone-evoked FR between light-on and light-off trials (N = 66, Figure 4—figure supplement 7). Within these matched subpopulations, PV and SOM photosuppression exhibited differential effects similar to those of the whole population. The change in FR due to PV suppression was not significantly different between responses to the standard and deviant (p2 > 0.05, t(65) = −0.3, C = 3). By contrast, the change in deviant tone-evoked FR due to SOM suppression was significantly weaker than that for the standard tone (Δ = −78%, p2 = 0.003, t(3.5), C = 3). By the design of the analysis, the effect of PV or SOM suppression on standard tone-evoked FR was nearly identical (p1 > 0.05, t(65) = −0.1, C = 3). However, the change in deviant tone-evoked FR was greater for PV photosuppression than SOM photosuppression (Δ = 404%, p1 = 0.029, t(65) = 2.4, C = 3). Since the observed differential effects of PV and SOM suppression persisted in subsets of neurons that were matched for photosuppression-induced change in standard tone-evoked FR, these differences are unlikely due to differential expression or efficacy of Arch in the PV-Cre and SOM-Cre mice, but rather reflect functional differences between the two types of interneurons.
Within the oddball sequence, after the presentation of the deviant tone, SSA takes several repeats of the standard tone to reach an adapted state (Ulanovsky et al., 2004). Consistent with previous findings (Ulanovsky et al., 2004), presentation of the deviant tone temporarily reduced SSA without photosuppression (Figure 5A–C, dark color bars); following the deviant tone (T0), the first two standard tones (T1 and T2) evoked elevated FRs compared to the fourth standard tone (T4) (PV-Cre, Figure 5B—N = 148, T1: Δ = 60%, p2 = 3e−8, t(146) = 6.3, C = 11, T2: Δ = 26%, p2 = 0.043, t(146) = 2.9, C = 11. SOM-Cre, Figure 5C—N = 102, T1: Δ = 72%, p2 = 1e−5, t(101) = 5.2, C = 11, T2: Δ = 31%, p2 = 0.013, t(101) = 3.3, C = 11). The third standard tone (T3) and the tone prior to the deviant tone (T−1) evoked responses similar to T4 (PV-Cre, Figure 5B—T−1 and T3: p2 > 0.05, t(146) < 2.5, C = 11. SOM-Cre, Figure 5C—T−1 and T3: p2 > 0.05, t(101) < 2.9, C = 11). Neurons in which response to T0 did not produce spikes were excluded. Suppressing PVs led to a significant and equal increase in FR to four consecutive presentations of the standard following the deviant (Figure 5B, left, for each tone, T−1 through T4, with light-on compared to T4 with light-off: Δ > 132%, p2 < 2e−9, t(146) ≥ 6.8, C = 11. Figure 5B, right, change in FR between light-on and light-off responses to each T−1 through T3 as compared to T4: p > 0.05, t(146) < 1.8, C = 5). In contrast with PVs, suppressing SOMs led to a progressively increasing effect on FR to consecutive presentations of the standard tone following the deviant (Figure 5C, left, for each standard tone, T−1 through T4, with light-on compared to T4 with light-off: Δ > 64%, p2 < 9e−4, t(101) ≥ 4.1, C = 11. Figure 5C, right, difference between FR change in T1 and T4 with light-on: p = 0.008, t(101) = −3.2, C = 5. Repeated measures ANOVA with tone number (T1 through T4) as a factor: F(3, 300) = 4.30, p = 0.0054). These results are consistent with the interpretation that the inhibitory drive from PVs is constant throughout the stimulus regardless of tone history, whereas the effect of SOM modulation increases with repeated presentations of the standard tone.
The time course of the effect of interneuron photosuppression on FR of the putative excitatory neurons at the beginning of each oddball sequence exhibited similar differences between PVs and SOMs. After the onset of each oddball sequence, SSA develops over the course of several standard tone presentations (Ulanovsky, 2004). As expected, on light-off trials, FR decreased in response to the standard tone over the first 20 repetitions of the tone (Figure 5—figure supplement 1). For PV-Cre mice, the difference in FR to the standard tone between light-on and light-off trials did not change over this time and stayed positive for the remainder of the oddball stimulus (Figure 5—figure supplement 1A). Over the first 20 trials, FR adapted with a similar time course for both the light-on and light-off trials, so the change due to PV photosuppression in FR to standard stayed constant (Figure 5—figure supplement 1A,B). In contrast, for SOM-Cre mice, FR on light-on trials increased over the first 40 trials, whereas on light-off trials, it decreased (Figure 5—figure supplement 1A). As a result, the difference due to photo-manipulation in FR to the standard tone increased over the first 40 trials and then stayed consistently positive throughout the stimulus presentation (Figure 5—figure supplement 1C). These results demonstrate that the PV-mediated effect on putative excitatory neuronal responses did not change with repeated presentations of the standard tone, whereas the SOM-mediated effect increased with the repeated stimulus.
In order to understand how PVs and SOMs exert differential control of SSA in putative excitatory neurons, we used optogenetic tagging to identify the specific interneurons and to quantify whether PVs and SOMs exhibited SSA (Lima et al., 2009). Through targeted viral delivery to AC, we drove Channelrhodopsin-2 (ChR2) expression, which depolarizes neurons when stimulated by light, in either PVs or SOMs (Chow et al., 2010) (Figure 6A,D, Figure 6—figure supplement 1A). A modified AAV encoding anti-sense code for ChR2 and a fluorescent reporter, under the FLEX cassette, was injected into PV-Cre or SOM-Cre mice (Boyden et al., 2005; Sohal et al., 2009; Cardin et al., 2010; Zhang et al., 2010; Deisseroth, 2011) and resulted in specific expression of ChR2, localized to PVs or SOMs (Figure 6—figure supplement 1B, c PV-Cre; N = 183 neurons in 3 mice, specificity = 67 ± 1%, efficiency = 76 ± 5%. SOM-Cre: N = 202 neurons in 4 mice, specificity = 90 ± 3%, efficiency = 81 ± 4%). Neurons were identified as PVs or SOMs if they responded to brief (5 ms) flashes of light with spikes within 1.5–4.5 ms of laser pulse onset (Figure 6A,D).
Both PVs and SOMs exhibited SSA, evidenced by a significant reduction in standard tone-evoked FR compared to the deviant tone response (Figure 6B,C,E,F, PV: N = 16, Δ = −32%, p2 = 0.023, z = −2.5, C = 2. SOM: N = 28, Δ = −41%, p2 = 0.002, z = −3.3, C = 2. Signed-rank test). The SSA index was not significantly different between PVs and SOMs (Figure 6G, neurons responsive to both tones A and B—PV: N = 5, SOM: N = 12. PV and SOM: p2 > 0.05, C = 2. Rank sum test) and both were similar to the mean SSA index in putative excitatory neurons (Figure 6G—Exc: N = 67. Exc vs PV: p > 0.05, z = 0.7, C = 2. Exc vs SOM: p > 0.05, z = 0.4, C = 2). PVs and SOMs exhibited some differences in relative response changes between the deviant, the standard, and the equal tones (Figure 6, Figure 6—figure supplement 2B,D); PVs' response to the equal tones did not decrease significantly as compared to deviant tones (N = 16 p2 > 0.05, z = −1.7, C = 2), whereas SOMs adapted in their response to equal tones (Δ = −36%, p2 = 0.049, z = −2.3, C = 2), and then further to standard tones (N = 28, Δ = −49%, p2 = 0.022, z = −2.6, C = 2). These results suggest that SOMs may adapt at a faster time scale than PVs with repeated presentation of tones.
Our results of recordings from PVs and SOMs present a surprising finding that PVs and SOMs adapt in response to repeated tones, countering our initial hypothesis that SOMs saturate in responses to the deviant, or facilitate with repeated presentation of a tone. How can an adapting interneuron contribute to added adaptation in excitatory neurons? To address this question, we next developed a model of coupled excitatory–inhibitory neuronal populations. Excitatory and inhibitory neurons form tight mutually coupled networks in A1, and we hypothesized that through differential post-synaptic integration by excitatory neurons, interneurons can amplify adaptation in excitatory neurons.
As a proof-of-principle that would account for our findings that PVs and SOMs exhibit similar magnitude of SSA, yet have a differential effect on SSA in putative excitatory neurons, we constructed a simplified model of mutually coupled inhibitory–excitatory neuronal populations. We tested how responses of the model putative excitatory neurons are affected by manipulation of activity of PVs or SOMs (Figure 7A). Thalamocortical tone-evoked inputs were modeled including an adaptation term and resulted in reduced responses of excitatory, PV, and SOM populations to repeated tones (Figure 7—figure supplement 1A,B). The model replicated the differential effects of manipulation of PV and SOM activity on responses to standard and deviant tones in putative excitatory neurons (Figure 7B–E): when PVs were suppressed optogenetically, the responses to both the standard and the deviant tones increased (Figure 7B,C). By contrast, when SOMs were suppressed, although the spontaneous FR and standard tone-evoked FR were elevated, the responses to the deviant tone remained constant, whereas the responses to the standard tone increased (Figure 7D,E). SOMs have been shown to inhibit PVs (Cottam et al., 2013; Pfeffer et al., 2013; Sturgill and Isaacson, 2015). Including inhibition between SOMs and PVs did not affect the model outcome, with suppression of PVs resulting in suppression of excitatory responses to both the standard and the deviant, and suppression of SOMs driving specific suppression of excitatory responses to the standard, but not the deviant (Figure 7—figure supplement 2).
An explanation for the difference of the effects of PVs and SOMs can be provided by examining the combined transfer function between pre-synaptic inputs and post-synaptic activity of excitatory neurons separately for PVs and SOM suppression (Figure 7A, insets): light-driven modulation of PV activity has the same effect on excitatory neuron responses at spontaneous, standard tone-evoked, and deviant tone-evoked activity (Figure 7A, left inset). Spontaneous, standard, and deviant input levels all fall within the linear portion of the transfer function between inputs and change in the excitatory neuron activity. On the other hand, for SOMs, modulation of their activity in the deviant tone-evoked regime drives small to no changes in excitatory neuronal activity, whereas modulation of SOM activity in the spontaneous and standard tone-evoked regime drives significant changes in excitatory neuronal activity (Figure 7A, right inset). The deviant tone-evoked activity falls on the saturating part of the input–output transfer function, whereas the standard tone-evoked and spontaneous inputs fall on the linear part of the transfer function. Then, shifts in SOM inputs due to photosuppression evoke small changes during deviant tone responses, but larger changes during either standard or spontaneous activity. Either PV or SOM manipulation would result in reduction of combined SSA of excitatory neurons.
The majority of neurons in the auditory cortex selectively reduce their responses to frequent, but not rare sounds, exhibiting SSA. However, the cortical mechanisms involved in the production and stimulus-specificity of SSA within the auditory cortex are not well understood. Here, we found that, in addition to adaptation at the level of thalamocortical inputs, two distinct types of interneurons, PVs and SOMs, differentially contributed to SSA in the primary auditory cortex. Optogenetic suppression of either PVs or SOMs led to a reduction in SSA in putative excitatory neurons (Figure 3). Suppression of PVs led to an equal increase in the FR of the putative excitatory neurons in response to the standard tone and the deviant tone (Figure 4). By contrast, suppression of SOMs significantly increased the response to the standard tone but lacked a significant effect on the response to the deviant tone (Figure 4). This series of findings expands on the ‘adaptation in narrowly tuned units’ model, which proposes that repeated presentation of the standard stimulus drives adaptation within more narrowly tuned inputs, such as thalamocortical inputs (Mill et al., 2011; Taaseh et al., 2011; Nelken, 2014). Our data indicate dual effects of cortical inhibition on SSA: (1) PVs contribute to SSA by providing a constant level of inhibition, resulting in a relatively higher inhibitory drive during the presentation of the standard, as compared to the deviant. Taking into account the non-linear synaptic input to FR output function of a typical pyramidal neuron, the constant inhibition amplifies the effect of thalamocortical depression in suppressing the response of the neuron to repeated stimulus (Figure 7A). (2) The selective increase of the inhibitory drive from SOMs for standard stimulus as compared to the deviant stimulus responses might be explained by a shift in the non-linear transfer function between inputs to SOMs and their outputs to excitatory neurons, possibly due to facilitation of SOM-to-excitatory neuron synapses (Beierlein et al., 2003; Silberberg and Markram, 2007) (Figure 7A).
Surprisingly, we found that, despite the differential effect of PV and SOM suppression on tone-evoked responses in putative excitatory neurons, both PVs and SOMs exhibit SSA. This finding is consistent with previous results that found that thalamocortical synapses onto inhibitory neurons and synapses from inhibitory neurons to excitatory cells can be depressing (Tan et al., 2008; Ma et al., 2012). How does suppression of these interneurons result in differential reduction in SSA in excitatory neurons? Our model provides an intuition for this effect: the mutually coupled excitatory–inhibitory network model demonstrates that the observed differential effects of PV and SOM suppression may be due to their differential action on excitatory neuronal responses in the unadapted and adapted state (Figure 7). Tone-evoked responses of PVs would fall on the linear portion of the transfer function between PV activity and excitatory neuron depolarization, while the same tones maximally affect inputs from SOMs onto excitatory neurons, with SSA shifting the inputs to the linear, more sensitive range of inputs from SOMs. Thus, SSA may serve an additional function: to adjust the responses of neurons in a range that is more sensitive to small changes in the inputs from both excitatory and inhibitory neuronal populations. More generally, the simulation demonstrates that a circuit element, such as PVs or SOMs, that itself adapts may further amplify adaptation in the excitatory neurons.
To estimate the differential contribution of PVs or SOMs inputs to the excitatory neurons, we measured the difference in the FR of neurons due to optogenetic partial suppression of their firing. This measurement provides an estimate of the change in the FR of the putative excitatory neurons with the change in combined inputs to the inhibitory neurons, thereby allowing estimation of the synaptic transfer function (Figure 7A, insets). A simple biologically plausible network incorporating these transfer functions can reproduce the observed responses (Figure 7). There are several caveats to this interpretation. First, the FR may not linearly translate onto synaptic input strength because of the spiking non-linear rectification between the inputs and outputs of the putative excitatory neuron: a small change in FR in the low-FR regime might correspond to a greater change in the synaptic drive than a similar-sized change in FR in the high-FR regime. However, our findings would still hold were this the case: in examining the effect of SOM suppression on response to the deviant, the actual difference in the synaptic drive between the deviant and the standard would then be even greater than observed. At the other end of the non-linearity, the analysis of neuronal responses sorted based on their FR to the standard tone and the deviant tone revealed that the ‘ceiling effect’ would not contribute to a decreased effect of photostimulation on the response to the deviant in SOM-Cre mice (Figure 4—figure supplements 4–7). Second, PVs and SOMs may inhibit not only the excitatory neurons, but also each other. SOMs make synapses onto PVs (Isaacson and Scanziani, 2011; Ma et al., 2012; Cottam et al., 2013), thereby potentially suppressing them with repeated presentation of the standard. Therefore, when SOMs are suppressed, some PVs may be disinhibited and provide a stronger suppression of excitatory neurons. The null effect on responses to the deviant tone during SOM suppression could result from a combination of increase in inhibition from disinhibited PVs in addition to reduced inhibition of SOMs onto excitatory neurons. Including inhibition from SOMs to PVs in the proof-of-principle model supported experimental findings (Figure 7—figure supplement 2). Third, other interneuron types, such as vasopressin-positive interneuron may be involved in the circuit (Pi et al., 2013), and the changes that we observe may reflect several inhibitory stages of processing.
One must be cautious in translating the data from our experiments as a strict description of neuronal activity in awake animals, as our results were based on recordings from mice under light isoflurane anesthesia. Other forms of anesthesia, such as pentobarbital-based (Cheung et al., 2001; Gaese and Ostwald, 2001), ketamine (Otazu et al., 2009) and high concentrations of isoflurane (Cheung et al., 2001; Ter-Mikaelian et al., 2007), can affect multiple aspects of sound-evoked responses in the auditory cortex. Nonetheless, our results are likely to extend for awake mice, since isoflurane anesthesia-induced effects on neuronal activity decrease as the concentration of isoflurane is reduced to the levels used in our recordings (Land et al., 2012). In addition, all recordings and manipulations were performed under identical anesthetic conditions, and our conclusions are based on the relative comparison of the effects of suppressing PVs and SOMs, which are expected to hold under awake conditions (Centanni et al., 2013).
While not demonstrated directly, SSA has been linked to detection of deviant sounds (Ulanovsky et al., 2003), which may be facilitated by a relatively enhanced neuronal response to a change in the ongoing sound (Nelken et al., 2003; Winkler et al., 2009; Grimm and Escera, 2012). By suppressing the responses to a frequently presented tone, the responses of neurons to a rare stimulus become relatively enhanced. However, whether and how modulating SSA in the auditory cortex affects auditory behavior has not yet been tested. Inhibitory interneurons may prove to have a complementary role in shaping auditory perception in addition to receptive field reorganization driven by synaptic plasticity (Froemke et al., 2013). The use of optogenetic methods to test the function of inhibitory interneurons in SSA overcomes the limitations of lesion or pharmacological studies (Elliott and Trahiotis, 1972; Duque et al., 2014), which only allow for prolonged, non-selective inactivation (Moore et al., 2001). By combining optogenetic manipulation of interneuron activity with behavioral measurements, future experiments will explore whether interneuron-mediated SSA indeed affects the auditory behavior of the subject, such as enhanced ability to detect unexpected events.
All experiments were performed in adult male mice (Jackson Laboratories, Bar Harbor, ME, United States; age, 12–15 weeks; weight, 22–32 g; PV-Cre mice, strain: B6;129P2-Pvalbtm1(cre)Arbr/J; SOM-Cre: Ssttm2.1(cre)Zjh/J) housed at 28°C on a 12-hr light:dark cycle with water and food provided ad libitum. In PV-Cre mice, Cre recombinase (Cre) is expressed in parvalbumin-positive interneurons; in SOM-Cre mice, Cre is expressed in somatostatin-positive interneurons (Taniguchi et al., 2011). This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to a protocol approved by the Institutional Animal Care and Use Committee of the University of Pennsylvania (Protocol Number: 803266). All surgery was performed under isoflurane anesthesia, and every effort was made to minimize suffering.
Modified AAVs were obtained from Penn Vector Core (Philadelphia, PA, United States). Modified AAV encoding Arch under FLEX promoter was used for selective suppression of PVs or SOMs (catalog number AV-9-PV2432, AAV9.CBA.Flex.Arch-GFP.WPRE.SV40, Addgene22222, serotype 2/9) (Chow et al., 2010). Modified AAV encoding GFP alone under FLEX promoter was used as a control for the specific action of Arch on the neuronal populations (catalog number AV-9-ALL854, AAV9.CAG.Flex.eGFP.WPRE.bGH, Allen Institute 854, serotype 2/9). Modified AAV encoding ChR2 under FLEX promoter was used for selective excitation of PVs or SOMs (catalog number AV-9-18917P, AAV9.CAGGS.Flex.ChR2-tdTomato.WPRE.SV40, Addgene18917, serotype 2/9).
2–3 weeks prior to the start of experimental recordings, a 0.5-mm diameter craniotomy was drilled over primary auditory cortex (2.6 mm caudal and 4.1 mm lateral from bregma) under aseptic conditions while the mouse was anesthetized with isoflurane. A 750 nl bolus of AAV in water was injected into A1 (1 mm ventral from pia mater) using a stereotaxic syringe pump (Pump 11 Elite Nanomite, Havard Apparatus, Holliston, MA, United States). The craniotomy was covered with bone wax and a small custom head-post was secured to the skull with dental acrylic.
All recordings were carried out inside a double-walled acoustic isolation booth (Industrial Acoustics, Bronx, NY, United States). Electrodes were targeted to A1 on the basis of stereotaxic coordinates and in relation to blood vessels. In electrophysiological recordings, the location was confirmed by examining the click and tone-pip responses of the recorded units for characteristic responses of neurons in core auditory areas, as described previously by our group in the rat (Carruthers et al., 2013) and by other groups in the mouse (Linden and Schreiner, 2003; Guo et al., 2012; Marlin et al., 2015). While the electrodes were targeted to A1, some recordings may include data from the anterior auditory field, adjacent to A1 (Linden et al., 2003). Mice were placed in the recording chamber, anesthetized with isoflurane, and the headpost secured to a custom base, immobilizing the head. After drilling a craniotomy and creating a durotomy exposing auditory cortex, a silicon multi-channel probe (A1x32-Poly2-5mm-50s-177 [Poly-2] or A1x32-tri-5mm-91-121-A32 [Triode], Neuronexus Ann Arbor, MI, United States) was slowly lowered to between 750 µm and 1 mm into the cortex, perpendicular to the cortical surface and used to record electrical activity. Raw signals from 32 channels were bandpass filtered at 600–6000 Hz and thresholded for spike analysis, or at 10–300 Hz for LFP and CSD analysis (Poly-2 probe only), digitized at 32 kHz and stored for offline analysis (Neuralynx, Bozeman, MT, United States). Common-mode noise was removed by referencing a probe inserted in the brain outside the auditory cortex. On the Poly-2 probe, two rows of 16 electrodes each on a single shank were arranged such that each electrode site was 50 µm away from all three closest neighbors. This arrangement allowed us to record densely across depth, that is, one electrode for every 25 µm in depth. On the triode, electrodes were arranged in groups of three equidistant sites, forming an equilateral triangle (25-µm separation). The triodes were separated vertically by 91-µm center-to-center distance, spanning 1 mm, with two additional single sites, one on each end.
Spike sorting was performed using commercial software (Offline Sorter, Plexon, Dallas, TX, United States) (Carruthers et al., 2013). In order to improve isolation of single units from recordings using low-impedance probes, spiking activity was sorted across three (triode, 25-µm separation) or four (poly-2, 50-µm separation) adjacent electrode sites (Niell and Stryker, 2008; Olsen et al., 2012). We used a stringent set of criteria to isolate single units from multiunit clusters (Otazu et al., 2009; Bizley et al., 2010; Brasselet et al., 2012; Durand et al., 2012; Carruthers et al., 2013; Picard et al., 2014; Carruthers et al., 2015). Single-unit clusters contained <1% of spikes within a 1.0-ms interspike interval, and the spike waveforms across 3 or 4 channels had to form a visually identifiable distinct cluster in a projection onto a three-dimensional subspace. Putative excitatory neurons were identified based on their expected response patterns to sounds and the lack of significant suppression of the spontaneous FR due to light (Lima et al., 2009; Moore and Wehr, 2013). While this subpopulation may still contain inhibitory neurons, only 2% of all recorded neurons were significantly photosuppressed at baseline (one-sided paired t-test, significance taken at p < 0.05). The low impedance of the extracellular probes precluded us from conducting a more detailed analysis of cortical subpopulations based on the spike waveform (Bartho et al., 2004; Moore and Wehr, 2013).
Stimuli were delivered via a magnetic speaker (Tucker-David Technologies, Alachua, FL, United States), directed toward the mouse's head. Speakers were calibrated prior to the experiments to ±3 dB over frequencies between 1 and 40 kHz, by placing a microphone (Brüel and Kjaer, Denmark) in the location of the ear contralateral to the recorded A1 hemisphere, recording speaker output and filtering stimuli to compensate for acoustic aberrations (Carruthers et al., 2013). First, to measure tuning, a train of 50 pure tones of frequencies spaced logarithmically between 1 and 80 kHz, at 65-dB sound pressure level (SPL) relative to 20 µPa, in pseudorandom order, was presented 20 times. Each tone was 100-ms long, with an inter-stimulus interval (ISI) of 300 ms. Frequency response functions were calculated online for several multiunits, and two frequencies (separated by 0.39 octaves), which elicited spiking responses of similar strength, were selected as tone A and B. Next, a series of stimuli composed of tones A and B were presented in interleaved blocks, repeated four times. Each oddball stimulus consisted of a train of 653 A and B tones (100-ms long, 300-ms ISI, 65-dB SPL). In oddball stimulus 1, 90% of the tones were A (standard), while 10% of the tones were B (deviant). We used a frozen sequence of standard and deviant tones in pseudorandom order and counterbalanced with respect to the number of standard tones preceding each deviant. In oddball 2, the probabilities of tones A and B were reversed so that tone B was the standard and A the deviant. In the equal probability stimulus, A and B each comprised 50% of tones.
An optic fiber was use to direct 532-nm laser light (Shanghai Laser & Optics Century, China). After positioning the silicon probe, an optic fiber was placed over the surface of auditory cortex. To limit Becquerel effect artifacts due to light-striking electrodes, we positioned the optical fiber parallel to the silicon probe (Han et al., 2009; Kvitsiani et al., 2013). During every fifth tone of the oddball and equal probability stimuli, light was cast over A1 to suppress interneurons. The light onset was 100 ms prior to tone onset, and lasted for 250 ms. At 180 mW/mm2, light pulses were intense enough to significantly modulate multiunit activity throughout all cortical layers. The effect of optical stimulation was not significant for responses to subsequent tones (Figure 3—figure supplement 1).
Brains were post-fixed in paraformaldehyde (4%, PFA) and cryoprotected in 30% sucrose. Coronal sections (40 µm) were cut using a cryostat (CM1860, Leica, Allendale, NJ, United States), washed in PBS containing 0.1% Triton X-100 (PBST; three washes, 5 min), incubated at room temperature in blocking solution (for PV, 10% normal goat serum and 5% bovine serum albumin in PBST; for SOM, 10% normal goat serum with 0.1% sodium azide and 2% cold water fish gelatin in PBS; 3 hr), and then incubated in primary antibody diluted in blocking solution overnight at 4°C. The following primary antibodies were used: anti-PV (PV 25 rabbit polyclonal, 1:500, Swant, Switzerland) or anti-SOM (AB5494 rabbit polyclonal, 1:200, Millipore, Billerica, MA, United States). After incubation, sections were washed in blocking solution (three washes, 5 min), incubated for 2 hr at room temperature with secondary antibodies (Alexa 594 goat anti-rabbit IgG; for PV 1:1000 and SOM 1:400), and then washed in PBS (three washes, 5 min each). Sections were mounted using Fluoromount-G (Southern Biotech, Birmingham, AL, United States) and confocal images were acquired (Leica SP5). Cells were identified in independent fluorescent channels and subsequently scored for co-localization by hand using ImageJ's cell counter plug-in. Transfection efficiency is the percent of antibody-labeled neurons, which are co-labeled with GFP. Transfection specificity is the percent of GFP-expressing neurons, which are co-labeled with the antibody.
Acute brain slices were prepared from mice using standard techniques essentially as previously described (Goldberg et al., 2011). Mice were anesthetized via inhaled isoflurane and then transcardially perfused with 10 ml of oxygenated, ice-cold artificial cerebrospinal fluid (ACSF) at a rate of 5 ml/min, that contained, in mM: 87 NaCl, 75 sucrose, 2.5 KCl, 1.25 NaH2PO4, 26 NaHCO3, 10 glucose, 0.5 CaCl2, 4 MgSO4. Slices (300-µm thick) were cut on a Leica VT1200S and incubated in cutting solution in a holding chamber at 32°C for approximately 30 min followed by continued incubation at room temperature prior to electrophysiological recording, at which point slices were transferred to a submersion-type recording chamber attached to the microscope stage. ACSF used for recording contained, in mM: 125 NaCl, 2.5 KCl, 1.25 NaH2PO4, 26 NaHCO3, 10 glucose, 2 CaCl2, and 1 MgSO4. The solution was continuously bubbled with 95% O2 and 5% CO2 throughout cutting, slice incubation, and recording, so as to maintain a pH of approximately 7.4.
Cells were identified via GFP expression under epifluorescence microscopy and subsequently visualized using a 40×, 0.8 NA water-immersion objective (Olympus, Center Valley, PA, United States) on an Olympus BX-61 upright microscope equipped with infrared differential interference contrast optics. Recordings were performed using the whole-cell patch clamp technique. Access resistance (Ra) was <25 MΩ upon break-in; data obtained from a given cell were rejected if Ra changed by >20% during the course of the experiment. Internal solution contained, in mM: potassium gluconate, 130; potassium chloride, 6.3; EGTA, 0.5; MgCl2, 1.0; HEPES, 10; Mg-ATP, 4; Na-GTP, 0.3; biocytin, 0.1%. Osmolarity was adjusted to 285–290 mOsm using 30% sucrose. Voltage was recorded using a MultiClamp 700B amplifier (Molecular Devices, Union City, CA, United States), lowpass filtered at 10 kHz, digitized at 16-bit resolution (Digidata 1550, Axon Instruments, Sunnyvale, CA, United States), and sampled at 20 kHz. pCLAMP 10 software was used for data acquisition, and analysis was performed using the Clampfit module of pCLAMP.
Cells were illuminated with a 561-nm solid state laser (Coherent, Santa Clara, CA, United States) routed to the standard X-Y galvanometer of a two-photon microscope (Bruker Corporation, Billerica, MA, United States) via a single-mode fiber. Illuminance at the specimen was estimated using a 10-µm pinhole aperture (Edmund Optics, Barrington, NJ, United States) and a photodiode power sensor (Thorlabs, Newton, NJ, United States).
For each putative excitatory neuron, the spontaneous FR and tone-evoked FRs were measured as the mean FR over 50 ms pre- and post-tone onset, respectively. For each identified interneuron, FRs were measured 100 ms pre- and post-tone onset. FR was measured separately for each tone, A and B, as standard, deviant, and equal probabilities, and for light-off and light-on trials. FR normalization was carried out separately for each tone, A and B, for each neuron by dividing the response under all conditions by the maximum FR (across 5-ms bins) of the deviant tone, light-off condition. Performing this normalization by dividing response in all conditions by the mean, rather than maximum FR of the deviant tone, light-off condition did not alter significant results (Figure 4—figure supplement 8). For all FR analyses, each neuron's responses to tones A and B were treated separately, and each was only included if the light-off deviant tone-evoked FR was significantly greater than the spontaneous FR (Wilcoxon signed-rank test p < 0.05). Further, tone responses were only included in analysis if the neuronal FR during each oddball stimulus exceeded 0.02 Hz, and the neuron was significantly tuned to the tone. Tuning was considered significant if the spike count in response to a tone (A or B) was significantly higher than the pool of spike counts across all tones outside one octave band centered on tones A and B (N = 42, t-test, p1 < 0.05). Population responses in each condition were measured as the mean and standard error of FRs across tone responses in each experimental group.
For each neuron, SSA index is a measure of the strength of SSA based on its mean FR with respect to tone probability. FRs to tones A and B were summed according to their standard or deviant probability within each oddball stimulus (Ulanovsky et al., 2003). Thus, SSA index was computed as:
where S and D indicate the mean FR for standard and deviant trials, respectively, and their subscripts indicate the tone frequency condition. SSA index was computed separately for light-off and light-on conditions. Population SSA indices were measured as the mean and standard error of SSA indices across all neurons of each population. Criteria for inclusion in the analysis were the same as in tone response FR analysis described above, with the added criterion that the deviant tone-evoked FR must be greater than spontaneous FR for each of tones A and B (Wilcoxon signed-rank test p < 0.05).
To calculate the CSD, the net current density moving through cortical tissue at 32 positions along the cortical axis was calculated based on LFPs of responses to tones recorded on each electrode, by using the second order central finite difference to calculate the second spatial derivative across the LFPs over the vertically arranged electrodes (Szymanski et al., 2009). Across the CSD profile, the deepest current sink corresponds to the thalamo-recipient granular layer (Kaur et al., 2005; Szymanski et al., 2009) allowing us to reconstruct the laminar location of recorded neurons. Neurons recorded on electrodes falling within the deepest sink were assigned to the granular layer, while those superior and inferior were assigned to the supra-granular and infra-granular layers (Figure 1E,F, Figure 4—figure supplement 3). The tone-evoked amplitude of the CSD was measured by first calculating root mean square of each channel during the first 50-ms post-tone onset, and then calculating the mean across all electrodes determined to fall within either the deepest short latency sink (granular layer) or pooled across all electrodes either above (supra-granular layer) or below (infra-granular layer). For each session, the granular layer CSD amplitude for all tone conditions was normalized across conditions by the deviant tone, light-off condition, and the mean across sessions was statistically analyzed. The SSA index was calculated as described in SSA index on the basis of the amplitude.
For all statistical tests in which N ≥ 30, we applied the Student's t-test (Matlab, Mathworks, Natick, MA, United States) unless specified otherwise, and reported the p-value, degrees of freedom, and t-statistic. For all tests with N < 30, sample variance was tested for normality using the Komogorov–Smirnov test. If any group's variance was non-normal, we applied a non-parametric test, for example, Wilcoxon sign rank or rank sum test (Matlab), and provided the z-statistic for any group with a normal distribution. For all tests, Bonferroni correction was applied for multiple comparisons and reported as ‘C = X’ where X is the factor by which the p-value was adjusted. Statistical tests were single-tailed if there was a reasonable prior expectation about the direction of the difference between samples. p1 refers to one-sided, and p2 refers to two-sided statistics set. In all figures, single, double, and triple stars indicate p < 0.05, 0.01, and 0.001, respectively. Error bars in all figures represent the standard error of the mean, unless otherwise noted.
We constructed models of the excitatory–inhibitory neuronal circuit to understand the coupling of excitatory interneurons with PV and SOM interneurons. We constructed firing-rate models based on Wilson–Cowan dynamics (Staiger et al., 1996; Xu et al., 2010, Rudy et al., 2011). The parameters were chosen in order to achieve a match to experimental data. The mean activity level of each population was modeled as:
where E(t) is the activity of the excitatory population; I(t) is the activity of the inhibitory population; S(x) is the transfer function between the combined ‘synaptic’ input and the neuronal FR. S(x) is linear with respect to intermediate inputs, but imposes a minimum and maximum activation limits. Sinh(x) is the transfer function between the inhibitory FR and excitatory post-synaptic current; jEI and jIE are excitatory–inhibitory and inhibitory–excitatory synaptic weights (0.2 and −1.0 for PVs, 0.05 and −0.39 for SOMs, respectively); jETone(t) and jITone(t) are tone-evoked input currents to excitatory and inhibitory neurons, respectively, modeled as 50-ms long exponentially decaying inputs of maximum amplitude 3 (delayed by 7 ms for SOMs, which do not receive direct thalamic inputs, relative to PVs, which receive direct thalamic inputs); τE (10 ms) and τI (10 ms) are synaptic time constants for excitatory and inhibitory neurons; k and r represent the maximum and minimum FR of neurons, respectively, (k = 15, r = 1); jinh(t) is the negative input to inhibitory neurons due to Arch. The optogenetic modulation was modeled as a unitary 250-ms pulse. To capture the differences in inputs due to repeated tone exposures, we modeled thalamic inputs reflecting the tone inputs with synaptic depression. We modeled the conductance of the thalamic projections, gInp, as changing according to the equation:
where g0 is the maximum conductance (g0 = 1), r is the gating coefficient representing tone-evoked thalamic input, Tg is the time scale for replenishment (Tg = 3 s), Tr is the time scale for depletion (Tr = 80 ms). We took r to be a step function with an exponential decay (with 40-ms time constant and amplitude of 3). The full input to auditory cortical neural populations is then equal to gInpr. In a train of four tones, the first tone-evoked response was taken as the deviant tone, and the fourth tone as the standard tone.
For the inhibitory-to-excitatory inputs, we used a sigmoidal transfer function and showed the existence of parameter regimes consistent with our results. For PVs, we used a sigmoid of the form:
where p = 0.3 and θ = 9. This gives a facilitating response at low input levels and a linear response at high input levels. For SOMs, we used a hyperbolic tangent that provided a saturating non-linearity:
where s = 3. For visualization, the baseline FR of neurons was removed and the peak response to a ‘deviant’ tone without optogenetic manipulation normalized to 1.
We also constructed a model with additional coupling between the PV and SOM interneurons using a generalization of the above dynamics, which may be written as:
where Ni is the FR of the ith population (EXC, PV, SOM), Ti = 10 ms is the time constant for each population, k = 15, r = 1, S has different maximum and minimum values for each population (xmin,E = −1, xmax,E = 1.75, xmin,PV = −0.5, xmax,PV = 4, xmin,SOM = 0, xmax,SOM = 3). SE(x) = x, and SSOM and SPV use the definitions above. jE,E = jSOM,SOM = jPV,PV = jPV,SOM = 0, jE,SOM = 0.25, jSOM,E = −0.25, jE,PV = 0.4, jPV,E = −1, and jSOM,PV = −0.1. jext,PV = 1.5, jext,SOM = 1. Tone inputs are the same as described above.
Stimulus-specific adaptation occurs in the auditory thalamusThe Journal of Neuroscience 29:7359–7363.https://doi.org/10.1523/JNEUROSCI.0793-09.2009
An optical neural interface: in vivo control of rodent motor cortex with integrated fiberoptic and optogenetic technologyJournal of Neural Engineering 4:S143–S156.https://doi.org/10.1088/1741-2560/4/3/S02
Long-lasting context dependence constrains neural encoding models in rodent auditory cortexJournal of Neurophysiology 102:2638–2656.https://doi.org/10.1152/jn.00577.2009
Possible principles underlying the transformation of sensory messagesPossible principles underlying the transformation of sensory messages.
Two dynamically distinct inhibitory networks in layer 4 of the neocortexJournal of Neurophysiology 90:2987–3000.https://doi.org/10.1152/jn.00283.2003
Neural ensemble codes for stimulus periodicity in auditory cortexThe Journal of Neuroscience 30:5078–5091.https://doi.org/10.1523/JNEUROSCI.5475-09.2010
Millisecond-timescale, genetically targeted optical control of neural activityNature Neuroscience 8:1263–1268.https://doi.org/10.1038/nn1525
Encoding of ultrasonic vocalizations in the auditory cortexJournal of Neurophysiology 109:1912–1927.https://doi.org/10.1152/jn.00483.2012
Emergence of invariant representation of vocalizations in the auditory cortexJournal of Neurophysiology, jn 00095 02015, 10.1152/jn.00095.2015.
Target-specific effects of somatostatin-expressing interneurons on neocortical visual processingThe Journal of Neuroscience 33:19567–19578.https://doi.org/10.1523/JNEUROSCI.2624-13.2013
Cortical interneurons: from Cajal to 2001Progress in Brain Research 136:215–238.
Neuronal circuits of the neocortexAnnual Review of Neuroscience 27:419–451.https://doi.org/10.1146/annurev.neuro.27.070203.144152
The logic of inhibitory connectivity in the neocortexThe Neuroscientist 19:228–237.https://doi.org/10.1177/1073858412456743
Long-term modification of cortical synapses improves sensory perceptionNature Neuroscience 16:79–88.https://doi.org/10.1038/nn.3274
Anesthesia changes frequency tuning of neurons in the rat primary auditory cortexJournal of Neurophysiology 86:1062–1066.
Shifting and scaling adaptation to dynamic stimuli in somatosensory cortexThe European Journal of Neuroscience 26:2359–2368.https://doi.org/10.1111/j.1460-9568.2007.05847.x
Auditory deviance detection revisited: evidence for a hierarchical novelty systemInternational Journal of Psychophysiology 85:88–92.https://doi.org/10.1016/j.ijpsycho.2011.05.012
Robustness of cortical topography across fields, laminae, anesthetic states, and neurophysiological signal typesThe Journal of Neuroscience 32:9159–9172.https://doi.org/10.1523/JNEUROSCI.0065-12.2012
Linking topography to tonotopy in the mouse auditory thalamocortical circuitThe Journal of Neuroscience 31:2983–2995.https://doi.org/10.1523/JNEUROSCI.5333-10.2011
Local versus global scales of organization in auditory cortexTrends in Neurosciences 37:502–510.https://doi.org/10.1016/j.tins.2014.06.003
Nonprimary auditory thalamic representation of acoustic changeJournal of Neurophysiology 72:1270–1277.
The largest group of superficial neocortical GABAergic interneurons expresses ionotropic serotonin receptorsThe Journal of Neuroscience 30:16796–16808.https://doi.org/10.1523/JNEUROSCI.1869-10.2010
Spectrotemporal structure of receptive fields in areas AI and AAF of mouse auditory cortexJournal of Neurophysiology 90:2660–2675.https://doi.org/10.1152/jn.00751.2002
Stimulus-specific adaptation in the inferior colliculus of the anesthetized ratThe Journal of Neuroscience 29:5483–5493.https://doi.org/10.1523/JNEUROSCI.4153-08.2009
Thalamocortical inputs trigger a propagating envelope of gamma-band activity in auditory cortex in vitroExperimental Brain Research 126:160–174.https://doi.org/10.1007/s002210050726
A neurocomputational model of stimulus-specific adaptation to oddball and Markov sequencesPLOS Computational Biology 7:e1002117.https://doi.org/10.1371/journal.pcbi.1002117
Parvalbumin-expressing inhibitory interneurons in auditory cortex are well-tuned for frequencyThe Journal of Neuroscience 33:13713–13723.https://doi.org/10.1523/JNEUROSCI.0663-13.2013
Primary auditory cortex of cats: feature detection or something else?Biological Cybernetics 89:397–406.https://doi.org/10.1007/s00422-003-0445-3
Stimulus-specific adaptation and deviance detection in the auditory system: experiments and modelsBiological Cybernetics, 10.1007/s00422-014-0585-7.
Highly selective receptive fields in mouse visual cortexThe Journal of Neuroscience 28:7520–7536.https://doi.org/10.1523/JNEUROSCI.0623-08.2008
Engaging in an auditory task suppresses responses in auditory cortexNature Neuroscience 12:646–654.https://doi.org/10.1038/nn.2306
Novelty detector neurons in the mammalian auditory midbrainThe European Journal of Neuroscience 22:2879–2885.
Aberrant development and plasticity of excitatory visual cortical networks in the absence of cpg15The Journal of Neuroscience 34:3517–3522.https://doi.org/10.1523/JNEUROSCI.2955-13.2014
Target-cell-specific facilitation and depression in neocortical circuitsNature Neuroscience 1:279–285.https://doi.org/10.1038/1092
Three groups of interneurons account for nearly 100% of neocortical GABAergic neuronsDevelopmental Neurobiology 71:45–61.https://doi.org/10.1002/dneu.20853
Selectivity and tolerance ("invariance") both increase as visual information propagates from cortical area V4 to ITThe Journal of Neuroscience 30:12978–12995.https://doi.org/10.1523/JNEUROSCI.0179-10.2010
Current source density profiles of stimulus-specific adaptation in rat auditory cortexJournal of Neurophysiology 102:1483–1490.https://doi.org/10.1152/jn.00240.2009
Robust but delayed thalamocortical activation of dendritic-targeting inhibitory interneuronsProceedings of the National Academy of Sciences of USA 105:2187–2192.https://doi.org/10.1073/pnas.0710628105
Transformation of temporal properties between auditory midbrain and cortex in the awake Mongolian gerbilThe Journal of Neuroscience 27:6091–6102.https://doi.org/10.1523/JNEUROSCI.4848-06.2007
Processing of low-probability sounds by cortical neuronsNature Neuroscience 6:391–398.https://doi.org/10.1038/nn1032
Multiple time scales of adaptation in auditory cortex neuronsThe Journal of Neuroscience 24:10440–10453.https://doi.org/10.1523/JNEUROSCI.1905-04.2004
Modeling the auditory scene: predictive regularity representations and perceptual objectsTrends in Cognitive Sciences 13:532–540.https://doi.org/10.1016/j.tics.2009.09.003
Laminar specificity of functional input to distinct types of inhibitory cortical neuronsThe Journal of Neuroscience 29:70–85.https://doi.org/10.1523/JNEUROSCI.4104-08.2009
Immunochemical characterization of inhibitory mouse cortical neurons: three chemically distinct classes of inhibitory cellsThe Journal of Comparative Neurology 518:389–404.https://doi.org/10.1002/cne.22229
Andrew J KingReviewing Editor; University of Oxford, United Kingdom
eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.
Thank you for submitting your work entitled “Complementary control of sensory adaptation by two types of cortical interneurons” for peer review at eLife. Your submission has been favorably evaluated by Gary Westbrook (Senior editor), and three reviewers, one of whom is a member of the Board of Reviewing Editors.
The following individuals responsible for the peer review of your submission have agreed to reveal their identity: Andrew King (Reviewing editor); Jennifer Linden (peer reviewer); Ben Willmore (peer reviewer).
The reviewers have discussed the reviews with one another and the Reviewing editor has drafted this decision to help you prepare a revised submission.
This paper investigates the contribution of two types of inhibitory interneuron to the generation of stimulus-specific adaptation (SSA) by neurons in the mouse auditory cortex. SSA is exhibited by neurons in non-lemniscal regions of the inferior colliculus and medial geniculate body, but appears to be a largely emergent property in the primary auditory cortex. Consequently, investigating the cortical circuitry that gives rise to this phenomenon is an important and largely unexplored question.
The authors show that optogenetic suppression of parvalbumin-positive (PV) cells produces a non-specific reduction of inhibition, which augments cortical responses to both deviant and standard tones by a similar amount (but produces an apparent reduction in the conventional “SSA index”, because the increase in response is proportionally greater for standard than for deviant tones). More interestingly, this work also shows that optogenetic suppression of somatostatin-positive (SOM) cortical interneurons specifically augments cortical responses to standard over deviant tones, and therefore reduces SSA. These results therefore suggest that PV cells contribute to SSA only through non-specific inhibition of responses to both deviant and standard tones, but that SOM cells selectively reduce cortical responses to frequently repeated standard tones and specifically amplify SSA.
The authors provide a large number of important and worthwhile controls to rule out alternative explanations for their findings, and also include a mathematical model to support their hypothesis that differences in the impact of PV and SOM inhibition on SSA could arise if the two interneuron types operate in different ranges of an essentially sigmoidal synaptic transfer function. The reviewers agree that this study provides a valuable insight into the mechanisms of SSA, but have raised a number of concerns about the analysis, modeling and presentation of the results.
1) Low impedance linear probes with recording sites separated by 50 µm were used. In a number of places, the authors refer to single unit activity. It is unlikely that the authors would have been able to isolate single units with low-impedance Neuronexus probes. At very least, more information about spike sorting and sorting quality, as well as the percentage of recordings that were actually multi-unit clusters need to be provided. The results of the present experiments are likely to be robust even if the proportion of multi-unit recordings is high, but this needs to be made clear in the manuscript (for example, for interpretation of the data in Figures 4 and 6).
2) Most of the figures show normalized firing rate measures, and most of the reported analyses were conducted on normalized firing rates. The Methods section indicates that firing rates were normalized for each neuron by dividing by the *maximum* response evoked by a deviant tone. Response measures based on a maximum (or minimum) are extremely sensitive to noise. It is very important to confirm that none of the key results would have been different had the authors used a more noise-robust normalization strategy, such as normalizing by the average response to deviant tones rather than the maximum.
3) The statistics in this paper are poorly presented, which undermines confidence in the results. Details of the statistics used are missing throughout, and need to be presented in a conventional format. The paragraph on “Statistical tests” in the Methods makes a general statement, but does not clarify what was done in individual cases. Multiple t-tests appear to be have been over used. In most places it would have been better to use a non-parametric test, such as the Wilcoxon tests that are referred to in 3 places. Where the use of multiple t-tests can be justified, it is essential to state when one- or two-tailed tests were used, and to provide effect sizes and indicate when p-values were corrected for multiple comparisons. Faith in the authors' statistical analysis is not helped by the fact that greater-than (>) signs appear to have been used where less-than (<) signs were intended in many cases (e.g. p>0.05).
4) The known impact of isofluorane and other volatile anesthetics on temporal aspects of neuronal responses – as demonstrated in rodent barrel cortex – could be a significant issue for SSA measurements Published data in cats suggests that isoflurane suppresses responses, increases latencies and increases thresholds (Cheung et al. 2001). Whilst acknowledging that the effects of modulating the activity of two different inhibitory interneuron populations were studied under the same anesthetic conditions, this issue should be mentioned in the Discussion.
5) It appears that the parameters of the mathematical model (e.g. for synaptic transfer functions for simulated PV and SOM cells) were not taken from experimental data, but rather were chosen to ensure that model replicated as closely as possible the impact of PV and SOM manipulations on SSA. There is a reference in the third paragraph of the Discussion, to direct estimation of the synaptic transfer function for PV and SOM cells, but it is not evident what data were used for this estimation (and the referenced insets in Figure 7A have no quantitative axis labels). Justification for the parameters and sigmoid functions used needs to be provided. Also the model assigns equivalent connectivity patterns between PV, SOM and excitatory neurons. Although both types of interneuron receive local excitatory inputs, only PV interneurons receive thalamocortical and feedforward corticocortical inputs and PV neurons provide largely unidirectional inhibition of SOM neurons. This should be taken into account in the model. At the very least, the authors should discuss this issue as a major limitation of the model (along with the fact that there is no attempt made to model differences in the location of PV versus SOM synaptic connections of the proximal versus distal parts of the excitatory cell's dendritic tree). While it is interesting that differences in the operating ranges for PV versus SOM synaptic transfer functions might in principle explain the observed differences in their impact on SSA, the model does not “reveal” mechanisms by which the different interneuron types exert their effects – and could not do so unless the underlying parameter choices were independently validated experimentally. We are not recommending that the modelling should be removed completely, but this section of the paper is currently rather weak.
6) The way the paper is written implies that the origin of SSA in the cortex is explained here (e.g. first sentence of Discussion). In fact, however, it seems that the explanation presented here depends on SSA already being present in excitatory, SOM and PV neurons. This paper is therefore really more about how PV/SOM neurons amplify SSA effects, rather than the origin of SSA itself. This needs to be clarified throughout. In particular, the PV interneurons appear to provide a general amplification of inhibition, rather than having any role that is specific to SSA.
7) The authors show using the expression of chanelrhodopsin in cre-PV and cre-SOM mice that both types of interneuron exhibit SSA. This is an important result, but are there any differences in the SSA exhibited by these interneurons and by cortical excitatory neurons? The hypothesis put forward in the Introduction that SOM neurons show facilitation to repetitive stimulation – and are therefore likely to have a greater suppressive effect on the response of excitatory neurons to the frequently presented standard tones – implies that this effect would be mediated at the level of the synaptic inputs to the SOM neurons. Consequently, SOM neurons should respond differently from PV neurons to repeated stimulation. Moreover, if, as suggested in the subsection “PVs and SOMs differentially suppress putative excitatory neuron responses to standard and deviant tones”, SOM neurons are firing at saturation to deviant tones, the authors should have the data to show this. It would strengthen the paper if the authors could show that SOM neurons eventually facilitate with repetitive stimulation, and that SOM cell responses are saturated for deviant tones.https://doi.org/10.7554/eLife.09868.028
- Maria Neimark Geffen
- Maria Neimark Geffen
- Maria Neimark Geffen
- Maria Neimark Geffen
- Mark Aizenberg
- Maria Neimark Geffen
- Ethan M Goldberg
- Ryan G Natan
- John J Briguglio
- John J Briguglio
- Maria Neimark Geffen
- Ethan M Goldberg
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
This work was supported by NIH R03DC013660, NIH R01DC014479, Klingenstein Award in Neuroscience, Human Frontier in Science Foundation Young Investigator Award, and the Pennsylvania Lions Club Hearing Research Fellowship to MNG and NARSAD Young Investigator Award to MA. RGN is supported by NIH NIMH T32MH017168. JJB is supported by NSF PHY1058202 and US-Israel BSF 2011058. EMG is partially supported by NINDS K12 NS049453 and the Burroughs Wellcome Fund Career Award for Medical Scientists. MNG is the recipient of the Burroughs Wellcome Fund Award at the Scientific Interface. EMG is the recipient of the Burroughs Wellcome Fund Career Award for Medical Scientists. The authors thank Jennifer Blackwell, Lisa Liu, Anh Nguyen, Daniel Feingold, Danielle Mohabir, Joshua Margolis, and Andrew Chen for technical support with experiments; Vijay Balasubramanian, Yale Cohen, Steven Eliades, Diego Contreras, Andrea Hasenstaub, and the members of the Geffen laboratory for insightful discussions and helpful comments; and Karl Deisseroth for providing the MTA for the use of some of the viruses.
Animal experimentation: All experimental procedures are in accordance with NIH guidelines and approved by the IACUC at University of Pennsylvania (protocol number 803266).
- Andrew J King, Reviewing Editor, University of Oxford, United Kingdom
© 2015, Natan et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.