Abstract
In the primary visual cortex (V1), orientation-selective neurons can be categorized into simple and complex cells primarily based on their receptive field (RF) structures. In mouse V1, although previous studies have examined the excitatory/inhibitory interplay underlying orientation selectivity (OS) of simple cells, the synaptic bases for that of complex cells have remained obscure. Here, by combining in vivo loose-patch and whole-cell recordings, we found that complex cells, identified by their overlapping on/off subfields, had significantly weaker OS than simple cells at both spiking and subthreshold membrane potential response levels. Voltage-clamp recordings further revealed that although excitatory inputs to complex and simple cells exhibited a similar degree of OS, inhibition in complex cells was more narrowly tuned than excitation, whereas in simple cells inhibition was more broadly tuned than excitation. The differential inhibitory tuning can primarily account for the difference in OS between complex and simple cells. Interestingly, the differential synaptic tuning correlated well with the spatial organization of synaptic input: the inhibitory visual RF in complex cells was more elongated in shape than its excitatory counterpart and also was more elongated than that in simple cells. Together, our results demonstrate that OS of complex and simple cells is differentially shaped by cortical inhibition based on its orientation tuning profile relative to excitation, which is contributed at least partially by the spatial organization of RFs of presynaptic inhibitory neurons.
SIGNIFICANCE STATEMENT Simple and complex cells, two classes of principal neurons in the primary visual cortex (V1), are generally thought to be equally selective for orientation. In mouse V1, we report that complex cells, identified by their overlapping on/off subfields, has significantly weaker orientation selectivity (OS) than simple cells. This can be primarily attributed to the differential tuning selectivity of inhibitory synaptic input: inhibition in complex cells is more narrowly tuned than excitation, whereas in simple cells inhibition is more broadly tuned than excitation. In addition, there is a good correlation between inhibitory tuning selectivity and the spatial organization of inhibitory inputs. These complex and simple cells with differential degree of OS may provide functionally distinct signals to different downstream targets.
- excitation-inhibition balance
- in vivo whole-cell recording
- orientation tuning
- receptive field
- synaptic input
Introduction
Orientation selectivity (OS) of neuronal responses is considered to be fundamental for visual perception of contours. In the primary visual cortex (V1), orientation-selective principal neurons are categorized into two distinct classes, simple and complex cells, based on their spike responses to either flashing or drifting stimuli (Hubel and Wiesel, 1962; Campbell et al., 1968; De Valois et al., 1982; Skottun et al., 1991; Niell and Stryker, 2008). The two cell types can be primarily distinguished by their different receptive field (RF) structures: simple cells have spatially segregated on and off subfields, while complex cells display overlapping on and off subfields (Hubel and Wiesel, 1962; Heggelund, 1986). Although simple and complex cells are generally considered to be equally selective for stimulus orientation, there have been results from several studies in cats and monkeys suggesting that complex cells are somewhat less selectively tuned than simple cells (Henry et al., 1974; Rose and Blakemore, 1974; Watkins and Berkley, 1974; Ikeda and Wright, 1975; Schiller et al., 1976; De Valois et al., 1982; Ringach et al., 2002).
The mechanisms for the potential differential degree of OS between simple and complex cells have not been explored previously. In the hierarchical model for visual processing (Hubel and Wiesel, 1962), it is thought that complex cells receive converging inputs from simple cells displaying similar orientation preferences, thus inheriting OS from the group of presynaptic neurons. It is certainly possible that the presynaptic simple cells do not perfectly register in orientation tuning profile, and that the convergence of inputs from them results in an averaging/smoothing effect, leading to the reduced tuning selectivity of the postsynaptic complex cell. This mechanism may be reflected by more weakly tuned excitatory input in complex than simple cells. On the other hand, in our previous study of simple cells in mouse V1, we have demonstrated that their orientation tuning is critically shaped by the interplay between moderately tuned excitation and even more broadly tuned inhibition as compared with excitation (Liu et al., 2011). The latter appears to play an essential role in sharpening OS of simple cells (Liu et al., 2011). Thus, an alternative mechanism could be that a differential excitatory/inhibitory interplay results in relatively weak selectivity of complex cells. To further understand the synaptic bases for OS in mouse V1, we examined complex cells primarily in layer 2/3 by combining in vivo cell-attached and whole-cell voltage-clamp recordings. Our results indicated that the relatively weak OS of complex cells compared with simple cells could be primarily attributed to an inhibitory synaptic mechanism: although there was no significant difference in excitation between simple and complex cells, complex cells received more narrowly tuned inhibitory inputs compared with their excitation, in contrast to the more broadly tuned inhibition in simple cells. Finally, the differential orientation tuning of inhibitory input correlated well with the spatial organization of its visual RF.
Materials and Methods
Animal preparation.
All experimental procedures used in this study were approved by the Animal Care and Use Committee of the University of Southern California. Female adult (12–16 weeks) C57BL/6 mice were used. Animals were sedated with an intramuscular injection of chlorprothixene (4 mg/ml) and anesthetized with urethane (1.2 g/kg, i.p., at 20% w/v in saline). Surgical procedure was performed to expose the visual cortex as previously described (Niell and Stryker, 2008; Liu et al., 2009, 2010). Eyelids were sutured during the surgical procedure. After the surgery, the right eyelid was reopened. The eye was rinsed with saline and a thin layer of silicone oil (30,000 centistokes) was applied to prevent drying while allowing clear optical transmission. Multiunit recordings were performed to determine the retinotopic map and location of V1, as previously described (Liu et al., 2009, 2010). The eye movement and RF drift of single units were negligible within the recording time windows (Mangini and Pearlman, 1980; Liu et al., 2009).
In vivo electrophysiology.
Cell-attached and whole-cell recordings were performed with an Axopatch 200B (Molecular Devices) according to previous studies (Moore and Nelson, 1998; Zhang et al., 2003; Liu et al., 2010). The patch pipette had a tip opening of ∼2 μm (4–6 MΩ). For most whole-cell recordings, we used a Cs+-based intrapipette solution containing the following (in mm): 125 Cs-gluconate, 5 TEA-Cl, 4 MgATP, 0.3 GTP, 8 phosphocreatine, 10 HEPES, 10 EGTA, 2 CsCl, 1 QX-314, 0.75 MK-801, pH 7.25. A K+-based intrapipette solution containing the following (in mm): 130 K-gluconate, 2 KCl, 1 CaCl2, 4 MgATP, 0.3 GTP, 8 phosphocreatine, 10 HEPES, 11 EGTA, pH 7.25, was used in some current-clamp recordings to measure membrane potential responses, resting membrane potential,and spike threshold. The pipette capacitance and whole-cell capacitance were compensated completely. The series resistance (25–50 MΩ) was compensated by 50–60% (at 100 μs lag) to achieve an effective series resistance of 10–25 MΩ. An 11 mV junction potential was corrected. Signals were filtered at 2 kHz for voltage-clamp recording, 5 kHz for current-clamp recording, and sampled at 10 kHz. The evoked excitatory and inhibitory synaptic currents were resolved by clamping the cell at −70 and 0 mV, respectively. As discussed before, our blind whole-cell recording method highly biases sampling toward pyramidal neurons (Liu et al., 2009, 2010, 2011). For cell-attached recordings, glass electrodes containing the ACSF were used. A 100–250 MΩ seal was formed on the patched cell. The pipette capacitance was completely compensated. Spikes were recorded in voltage-clamp mode, with a commend potential applied to achieve a zero baseline current. The spike signal was filtered at 10 kHz and sampled at 20 kHz. All neurons recorded in this study were located at a depth of 150–400 μm below the pia according to the microdrive reading, corresponding to layer 2/3 and a small top part of layer 4.
Visual stimulation.
Software for data acquisition and visual stimulation was custom-developed in LabVIEW (National Instruments) and MATLAB (MathWorks) respectively. Visual stimuli were provided by a 34.5 × 25.9 cm monitor (refresh rate 120 Hz, mean luminance ∼41.1 cd/m2) placed 0.25 m away from the right eye, which is equivalent to infinity (Liu et al., 2010). The center of monitor was placed at 45° azimuth, 0° elevation, and it covered ±35° horizontally and ±27° vertically of the mouse visual field, corresponding to the monocular zone of V1. To map the RF, either flashing bright (57.5 cd/m2) and dark (24.7 cd/m2) squares (11 × 11, each square was 5° × 5°) or flashing bars (15 bars, each bar was 4° × 60°) of the preferred orientation were applied (over the gray background, contrast = 95%) at different positions in a pseudorandom sequence. RFs were mapped for 5–10 repetitions. The stimulus duration was 200 ms and interstimulus interval was 240 ms. To measure orientation tuning, drifting sinusoidal gratings (2 Hz, 0.04 cycle/°, contrast 95%) or drifting bar (4° width, 60° length, 50°/s speed) of 12 directions (30° step) were applied. For drifting gratings, stationary grating of one orientation was first presented on the full screen for 1.8 s before it drifted for 1.5 s. The grating stopped drifting for 500 ms before another grating pattern appeared. Ideally, OS would be measured at the best spatial frequency. However, because we could only maintain a good recording within a limited time window and to apply all different combinations of orientation and spatial frequency took enormous time, we chose to measure OS at a spatial frequency of 0.04 cpd, as the largest fraction of mouse V1 cells prefers this spatial frequency and the average spatial frequency in layer 2/3 is 0.04 cpd (Niell and Stryker, 2008). It should be noted that this limitation does not affect our conclusions since measurements with single moving bars generated similar results (see below). For measuring response modulation (i.e., F1/F0 ratio), drifting sinusoidal gratings of the preferred direction (with temporal frequency of 2 Hz) were presented for 10 cycles, at various spatial frequencies (0.01, 0.02, 0.04, 0.08, 0.16, and 0.32 cycle/°). In total, ∼50–100 cycles were collected. F1/F0 ratio was calculated at the preferred spatial frequency. Drifting bar was moved across the screen at an interstimulus interval of 3 s. The 12 patterns were presented in a pseudorandom sequence, and were repeated 5–10 times. Spontaneous activity was recorded when a uniform gray background was presented.
Data analysis.
Spikes were sorted offline. Spikes evoked by flashing stimuli were counted within a 70–270 ms time window after the stimulus onset. Spikes evoked by drifting stimuli were counted within a 70–1570 ms window after the start of drifting. The spontaneous firing rate was subtracted from the stimulus-evoked spike rate. The cells included in the analysis should first show evoked spike responses to both flash stimuli and drifting bar/grating stimuli. Evoked responses are defined by firing rates exceeding the average spontaneous rate by three SDs of baseline fluctuations. Within this group of cells, cells that have well defined on and off subfields are included for analysis (81%). Subfields are identified as spatially contiguous pixels with evoked responses. Firing rates of other pixels in the RF were set as zero. The cells analyzed in this study had one on and one off subfield.
To quantify the spatial overlap between two dimensional on and off subfields, we first fit the hard boundary of the subfield with an ellipse (Li et al., 2013), and quantified the goodness of fitting by calculating the adjust r square. Most of our fittings (91% of all subfields) have a value larger than 0.9, indicating that the ellipses did well in fitting the subfields. Then we calculated an overlap index (Liu et al., 2009; Li et al., 2013): OI = (w1 + w2 − d)/(w1 + w2 + d), where OI is the overlap index, d is the distance between the centers of two ellipses, and w1 and w2 are the half-length of the line segment connecting the centers of two ellipses and intercepted by each ellipse. OI = 1 gives a complete overlap, and OI ≤ 0 gives a complete separation. To calculate OI for One-dimensional (1-D) subfields, the hard boundary of the subfields were determined in a similar manner as in the two-dimensional RFs. The envelope for the peak response amplitudes (i.e., the spatial tuning curve) was generated using MATLAB software Envelope 1.1 (Sun et al., 2010), where d was defined then as the distance between the peaks of on and off subfields, and w1 and w2 are the distance between the peak and the hard boundary of the subfield on the inner side. OI was calculated using the same equation as described above.
We also quantified the spatial overlap between on and off subfields using Pearson's correlation coefficient (rRF) according to the following question (Priebe et al., 2004; Mata and Ringach, 2005; Liu et al., 2009):
Ron,i and Roff,i and are individual responses to an on or off stimulus respectively. R̄on and R̄off are the mean on and off responses within the subfields.
To quantify the strength of OS, orientation selectivity index (OSI) and tuning width were calculated. The responses to two directions at each orientation were averaged to obtain the orientation tuning curve between 0 and 180°, which was then fit with a Gaussian function: R(θ) = A × e−0.5×(θ−θpref)2/σ2 + B. θpref is the preferred orientation determined by the vector sum of responses across orientations. OSI was quantified as follows: (Rpref − Rorth)/(Rpref − Rorth) = A/(A + 2B), where Rpref is the response level at the preferred orientation θpref, and Rorth is that at the angle of θpref + 90°. Tuning width was defined as the half-width at half-height of the fitted Gaussian curve above the offset B.
To calculate the F1/F0 ratio, first the poststimulus spike time histogram (PSTH) was generated from all the cycles for responses over multiple repetitions. F1 response was calculated from the PSTH as the amplitude of the best-fitting sinusoid at the modulation frequency (Mata and Ringach, 2005). F0 response was the mean evoked spike rate during the drifting grating stimulus.
In current-clamp recordings, subthreshold membrane potential responses were analyzed after removing spikes with a median filter at 8 ms (Carandini and Ferster, 2000). Putative complex cells were identified by OI of spike response ≥0.33 or OI of membrane potential response >0.7.
In voltage-clamp recordings, average excitatory and inhibitory response traces were first smoothed by averaging within a sliding 40 ms window (Li et al., 2012a,b). Excitatory and inhibitory synaptic conductance were derived according to the following equation (Wehr and Zador, 2003; Liu et al., 2010):
I(t) is the amplitude of current at any time point, Gr and Er are the resting leak conductance and resting membrane potential respectively, Ge and Gi are the excitatory and inhibitory synaptic conductance, respectively, Vm(t) is the membrane voltage, and Ee (0 mV) and Ei (−70 mV) are the reversal potentials. Vm(t) is corrected by Vm(t) = Vh − Rs × I(t), where Rs was the effective series resistance and Vh is the applied holding voltage. Measurement of currents at two different voltages yielded a system of two equations that could be solved for Ge and Gi at any particular t. We used the bootstrapping method (bootstrap sampling = 1000 times) to evaluate the variability of the derived conductance (Li et al., 2012a; Ma et al., 2013).
Neuron modeling.
A conductance-based integrate-and-fire neuron model was used to simulate the membrane potential response (Liu et al., 2010; Li et al., 2012a,b, 2014):
Vm(t) is the membrane potential at time t, C is the whole-cell capacitance, Gr is the resting leak conductance, and Er is the resting membrane potential (−65 mV). Ee and Ei were set at 0 mV and −70 mV, respectively. Gr was calculated based on the equation Gr = C × Gm/Cm, where Gm, the specific membrane conductance is 2e − 5 S/cm2, and Cm, the specific membrane capacitance is 1e − 6 F/cm2.
Synaptic conductance was simulated with a skew normal distribution function:
where ϕ and ∅ are the standard normal probability density function and the cumulative distribution function, respectively, and ξ determines the location. The scale factor (ω) was set at 145 ms and the shape factor (α) at 1.5 for both excitatory and inhibitory conductance. Inhibition was set 25 ms delay relative to excitation (Li et al., 2012a,b). Spike responses were simulated using a power law (Carandini and Ferster, 2000; Priebe and Ferster, 2008; Liu et al., 2011): R(Vm) = k ⌊ Vm − Vth ⌋+p. R is the firing rate, k is the gain factor, p is the exponent, and “+” indicates rectification.
Results
OS of complex cells is weaker than simple cells
To examine whether complex and simple cells in mouse V1 differ in the degree of orientation selectivity, we performed cell-attached loose-patch recordings from excitatory neurons mainly in layer 2/3 (see Materials and Methods). We mapped the spiking visual RF of each recorded cell to determine the simple/complex cell type, and to reveal the relationship between OS and the RF structure. Figure 1A shows the spatial patterns of spike responses of a cell to a set of bright (on) and dark (off) squares flashed for 200 ms. The largely overlapped on and off response regions (i.e., on and off subfields) indicated that the cell was a complex cell. Figure 1B shows the cell's responses to drifting sinusoidal gratings of different orientations and directions. The cell responded to gratings of all orientations, but responded maximally to gratings at a certain orientation angle. Figure 1C,D, shows the responses of another cell to flash and drifting stimuli. The on and off subfields of the cell were clearly segregated (Fig. 1C), indicating that the cell could be functionally identified as a simple cell. In response to drifting sinusoidal gratings, it exhibited apparently temporally modulated spiking responses which are sharply tuned for orientation (Fig. 1D).
Complex cells are less orientation selective than simple cells. A, Spike on/off subfields of an example complex cell examined by cell-attached recording. Top, Arrays of PSTHs for spike responses of the neuron to a set of bright (on) or dark (off) flash squares. PSTHs were arranged in 2-D space according to the locations of corresponding stimuli. Each pixel represents visual space of 5°. Red and green ovals depict the elliptic fitting of the on and off subfields, respectively. Spontaneous firing rate is 1.2 Hz. Scale bars, 20 Hz, 300 ms. Bottom, Color maps for spike on (red) and off (green) responses. The brightness of color represents the average evoked firing rate. The maps were smoothed by bilinear interpolation. The white line depicts the RF major axis, defined by the long axis of the oval. B, Orientation tuning of the same cell. Left, PSTHs for its spike responses to drifting gratings at 12 directions. Red arrow indicates the drifting direction. Scale bars, 20 Hz, 250 ms. Right, Orientation tuning curve of evoked spike numbers. The responses to gratings of opposite directions were averaged for each orientation. Red curve indicates the Gaussian fit. Blue dash line indicates the level of spontaneous spiking rate. C, D, An example simple cell. Data are displayed in a similar manner as in A and B. Spontaneous firing rate is 0.6 Hz. Scale bars: C, 21 Hz, 300 ms; D, 21 Hz, 300 ms. E, Distribution of overlap indices in the recorded cell population. The distribution is significantly bimodal (p < 0.05, Hartigan's dip test). Dash line marks OI = 0.33 (which separates simple from complex cells). F, Depth distribution of the recorded cells. G, Plot of F1/F0 ratio against OI. Blue dash line is the best-fit linear regression line. Correlation coefficient r = −0.74. Black dash lines mark OI = 0.33 and F1/F0 = 1, respectively. H, Plot of OSI against OI. OSI was measured from responses to either drifting gratings (open triangle) or drifting bars (gray dots). Blue dash line is the best-fit linear regression line (r = −0.63). Black dash line marks OI = 0.33. Inset, Comparison of mean OSIs of complex (black, n = 57) and simple (gray, n = 42) cells defined by OI values. ***p < 0.001, two-tailed t test. I, Plot of OSI against rRF measured between on and off subregions (r = −0.59). Black dash line marks rRF = 0, which separates complex from simple cells. Inset, Comparison of mean OSIs of complex (n = 54) and simple (n = 45) cells defined by rRF values. ***p < 0.001, two-tailed t test. J, Plot of tuning width against OI (r = 0.54). Inset, Comparison of mean tuning widths of complex and simple cells defined by OI values. ***p < 0.001, two-tailed t test. K, Plot of tuning width against rRF (r = 0.51). Inset, Comparison of mean tuning widths of complex and simple cells defined by rRF values. ***p < 0.001, two-tailed t test. L, Plot of angle of the major axis of on subfield versus that of off subfield for complex cells (n = 43). Red dash line is the identity line. Note that only complex cells with OSI ≥0.3 were selected. M, Plot of angle of the preferred orientation (measured with drifting gratings or bars) and that of RF major axis (measured with flash squares). Red dash line is the identity line.
We classified cell types based on the spatial overlap between on and off subfields, which was quantified with both an OI and a rRF (see Materials and Methods). Consistent with our previous report (Liu et al., 2010), OIs of our recorded excitatory cells exhibited a bimodal distribution (Fig. 1E), indicating that they could be categorized into two classes. We used OI = 0.33 as a criteria to separate complex cells from simple cells, following our previous study (Liu et al., 2009). In a total of 128 cells, 70 cells were classified as complex cells and 58 were classified as simple cells. These simple and complex cells slightly differ in depth distribution, in that complex cells were more frequently encountered than simple cells in shallower cortical locations (p < 0.01, two-tailed t test; Fig. 1F). To measure OS, drifting sinusoidal gratings were applied to 72 and drifting bars were applied to 27 of these cells. While all tested simple cells showed strong temporal response modulation to sinusoidal gratings, as reflected by a F1/F0 ratio larger than 1, only 16 of 41 complex cells exhibited a F1/F0 ratio smaller than 1 (Fig. 1G). These complex cells would then be classified as simple cells if based on F1/F0 ratio, another criterion for separating simple and complex cells (Skottun et al., 1991; Carandini et al., 1997; Anderson et al., 2000; Mata and Ringach, 2005). Therefore, the cell classifications based on the two criteria do not strictly agree with each other. There was a negative correlation between F1/F0 ratio and OI (Fig. 1G; correlation coefficient (r) = −0.74), reminiscent of a previous observation in macaque V1 (Mata and Ringach, 2005). We next quantified OS with both an OSI and orientation tuning width (see Materials and Methods), and examined the relationship between OS and the RF structure. We found that OSI was negatively correlated with both OI (Fig. 1H) and rRF (Fig. 1I), whereas tuning width was positively correlated with OI and rRF (Fig. 1J,K). Complex cells had significantly weaker OS (Fig. 1H,I) and broader tuning widths than simple cells (Fig. 1J, K). Although most of the simple cells exhibited strong OS (with OSI >0.5), there were some complex cells that were not orientation selective (defined by OSI <0.3).
For 43 orientation-selective complex cells (with OSI ≥0.3), we further examined the relationship between OS and the RF structure. We noticed that subfields of complex cells exhibited a slightly elongated shape (Fig. 1A). By fitting the subfield with an ellipse, we determined the major axis for the subfield (see Materials and Methods). The major axis for the on subfield was essentially identical to that for the off subfield (Fig. 1L). The RF major axis also agreed with the preferred orientation of the cell (Fig. 1M). Based on these results, in the following intracellular recordings we applied flashing bars at the cell's preferred orientation to map the RF and determine the cell type. Furthermore, we did not find any significant difference in either evoked (32.0 ± 24.2 Hz for simple cells; 37.2 ± 33.4 Hz for complex cells; p = 0.80, two-tailed t test) or spontaneous firing rate (2.3 ± 2.8 Hz for simple cells; 2.4 ± 4.0 Hz for complex cells; p = 0.96, two-tailed t test) between simple and complex cells. This indicates that complex and simple cells in our study have similar signal-to-noise ratios in their responses.
Subthreshold membrane potential responses of complex cells
We next performed in vivo whole-cell recordings with a K+-gluconate based internal solution to examine OS of subthreshold membrane potential (Vm) responses. We first measured the orientation tuning of Vm responses under current-clamp by applying drifting sinusoidal gratings at 12 directions. The spike responses of the cell were also recorded at the same time. One-dimensional RF structure was then measured by applying flashing bars of preferred orientation at different locations, and the cell type was determined based on the overlap between spike on and off responses in the 1-D space. The example cell in Figure 2A was identified as a complex cell, based on its largely overlapped on and off spike subfields. The cell's spiking response was moderately orientation tuned, whereas its subthreshold Vm response was much more broadly tuned (Fig. 2B, 2C). Nonetheless, the spike and Vm responses exhibited a similar preferred orientation (Fig. 2C). An example simple cell was displayed in a similar manner (Fig. 2D–F).
Orientation tuning of subthreshold membrane potential responses in complex and simple cells. A, PSTHs for spike responses (Sp, top) and average subthreshold membrane potential responses (Vm, bottom; after filtering out spikes) to 15 flash bars at different locations, plotted for an example complex cell examined by whole-cell current-clamp recording. Red dash lines indicate the level of resting membrane potential. Scale bars, 30 Hz (spike)/22 mV (Vm), 300 ms. B, PSTHs for spike responses (left) and average subthreshold Vm responses (right) to drifting gratings, plotted for the same cell in A. Scale bars, 40 Hz/21 mV, 200 ms. C, Orientation tuning curves of spike (top) and Vm (bottom) responses for the same cell. Red dash curves depict Gaussian fits. Vm response was measured as the peak depolarization relative to the resting membrane potential in the cycle-averaged waveform. D–F, An example simple cell. Data are displayed in a similar manner as in A-C. Scale bars: D, 40 Hz/24 mV, 300 ms; E, 60 Hz/36 mV, 200 ms. G, Plot of OI measured from spiking responses versus OI measured from Vm responses. Horizontal dash line marks OI_spike = 0.33. Vertical dash line marks OI_Vm = 0.7, which separates complex from simple cells. H, Plot of OSI measured from spike responses versus OSI measured from Vm responses. Red dash line is the identity line. Blue dash line is the best-fit linear regression line (r = 0.81). Horizontal dash line marks OSI_spike = 0.3. Vertical dash line marks OSI_Vm = 0.06, which separates orientation selective from nonselective neurons. I, Plot of angle of the preferred orientation measured from spike responses versus that from Vm responses, for complex cells with OSI_spike ≥0.3 only. Red dash line is the identity line. J, Distribution OSIs measured from Vm responses. There is a significant difference between simple and complex cells. p < 0.05, two-tailed t test. K, Plot of firing rate versus Vm response, fit with a power-law function, for an example complex cell. The power value is indicated. L, Plot of firing rate versus Vm response for an example simple cell. M, Comparison of powers between complex (n = 24) and simple (n = 25) cells. p = 0.28, two tailed t test. Bar = SEM. N, Comparison of spike threshold and resting membrane potential between complex (n = 24) and simple (n = 25) cells. p = 0.43, 0.47, respectively, two tailed t test. Error bar indicates SEM.
In 49 similarly recorded cells, 24 were identified as complex cells (spike OI >0.33) and 25 as simple cells (spike OI <0.33; Fig. 2G). We found a strong positive correlation between OI values calculated from spike and Vm responses (Fig. 2G, r = 0.91), indicating that OI calculated from the Vm response can be used as a measure to predict simple or complex cell type. Here we set OI_Vm = 0.7 as a separation point between simple and complex cells (Fig. 2G). As for orientation tuning, for both the complex and simple cells, OSI calculated from the spike response was much larger than that from the Vm response, consistent with the notion that the spike thresholding mechanism in general greatly enhances neuronal selectivity (Carandini and Ferster, 2000; Priebe and Ferster, 2008; Liu et al., 2011). In addition, for both the complex and simple cells, there was a strong positive correlation between OSIs calculated from spike and subthreshold Vm responses (r = 0.82 for complex cells; r = 0.83 for simple cells; r = 0.81 for all the cells; Fig. 2H), suggesting that OSI_Vm can be used as a measure to predict the degree of OS of the cell. If OSI_spike = 0.3 is set as a criterion for defining orientation selective neurons, it would correspond approximately to OSI_Vm = 0.06. The latter was then used as a criterion in the following intracellular recordings to identify putative orientation selective cells. Similar as the loose-patch recording results, the complex and simple cell groups differed significantly in the mean value of OSI_spike (p < 0.001, two-tailed t test), and a small subset of the complex cells were not orientation selective (Fig. 2H). For the orientation selective complex cells, their spike and subthreshold Vm responses displayed a similar preferred orientation (Fig. 2I), similar as simple cells (Liu et al., 2011). The Vm responses of complex cells were more broadly tuned than those of simple cells, although there was a large overlap between the distributions of OSI_Vm values of these two cell groups (Fig. 2J; p < 0.05, t test). To determine whether the transformation from membrane potential to spike would contribute to the differential OS between complex and simple cells, we fit the relationship between membrane potential and spike responses with a power-law function and calculated the power (Priebe et al., 2004; Priebe and Ferster, 2008). As shown by the example complex and simple cells, the data were well described by the function (Fig. 2K,L). No significant difference was found for the power between complex and simple cells (5.36 ± 0.86 for complex; 5.66 ± 0.71 for simple; p = 0.28, two-tailed t test Fig. 2M). Because we did not find significant difference in the level of spike threshold (21.8 ± 5.8 mV for complex; 21.4 ± 3.9 mV for simple; p = 0.43, two-tailed t test) or that of resting membrane potential (−61.5 ± 6.2 mV for complex; −61.2 ± 8.9 mV for simple; p = 0.47, two-tailed t test; Fig. 2N), the differential degree of OS between complex and simple cells could be primarily attributed to the differential tuning of their Vm responses.
Inhibition is more narrowly tuned than excitation in complex cells
We further dissected the synaptic mechanisms underlying the weaker orientation tuning of complex cells by performing whole-cell voltage-clamp recordings with a cesium-based internal solution (see Materials and Methods). Before switching recording to voltage-clamp, we mapped orientation tuning of Vm response with drifting bars as well as on and off subfields with flashing bars of the preferred orientation (see Materials and Methods). OI_Vm < 0.7 was used as a criterion to identify putative complex cells (Fig. 2G) and only cells with OSI_Vm > 0.06 (putative orientation selective neurons, 14 of 16 recorded neurons) were further analyzed. Excitatory and inhibitory synaptic currents evoked by drifting bars were then recorded under voltage-clamp, by clamping the cell's membrane potential at −70 and 0 mV, respectively. The cell shown in Figure 3, A–C, was a putative complex cell, since it displayed largely overlapped on and off subfields (Fig. 3C). Its Vm responses (traces not shown) exhibited an OSI of 0.18, indicating that it was a presumptive orientation selective neuron. We calculated evoked synaptic conductance (Fig. 3B, top inset) from the recorded currents (Fig. 3A), and found that the inhibitory input displayed narrower orientation tuning than the excitatory input (Fig. 3B) in this complex cell, as measured by the peak evoked synaptic conductance. On the other hand, in a putative simple cell, we found that the inhibitory input was more broadly tuned than the excitatory input (Fig. 3D–F).
Synaptic inputs to complex and simple cells. A, Average excitatory (Exc) and inhibitory (Inh) currents in an example complex cell evoked by drifting bars at 12 directions (5 trials). Red arrow indicates the moving direction. Scale, 150 pA (Exc)/210 pA (Inh), 300 ms. B, Orientation tuning curves for peak evoked excitatory (red) and inhibitory (blue) conductance for the same cell shown in A. Top, Inset, Traces of derived synaptic conductance at the preferred direction. Shading indicates 95% confidence level. C, 1-D RF of Vm responses (recorded in current-clamp) to flashing bars of preferred orientation, for the same cell in A. Each bar is 4° in width. Dash and solid curves in the bottom depict the spatial tuning (i.e., the envelope of peak response amplitudes) of on and off responses respectively. Scale, 8 mV, 200 ms. OI value is indicated. D–F, An example simple cell. Data are presented in the same as in A–C. Scales: D, 200 pA (Exc)/600 pA (Inh), 300 ms; F, 10 mV, 200 ms.
The average normalized tuning curves for a total of 14 similarly identified orientation selective complex cells showed that inhibition was more narrowly tuned than excitation (Fig. 4A), with a narrower tuning width compared with excitation (Fig. 4A, inset). In contrast, in a total of 12 simple cells, inhibition was more broadly tuned than excitation (Fig. 4B), with a broader tuning width compared with excitation (Fig. 4B, inset). These trends were observed in all individual complex and simple cells, as demonstrated by the plot of ratio of OSI for inhibition over that for excitation in the same cell versus OI measured from Vm responses (Fig. 4C). For complex cells, the mean OSI ± SD was 0.23 ± 0.08 for excitation and 0.33 ± 0.11 for inhibition (Fig. 4D; p < 0.001, two-tailed paired t test). For simple cells, the mean OSI was 0.26 ± 0.09 for excitation and 0.12 ± 0.10 for inhibition (Fig. 4D; p < 0.001, two-tailed paired t test). We did not find a significant difference in the tuning level of excitation between complex and simple cells (Fig. 4D; p = 0.37, two-tailed t test), as the normalized amplitude of excitation at the orthogonal orientation was 0.62 ± 0.14 for simple cells and 0.70 ± 0.22 for complex cells (Fig. 4A, B; p = 0.33, two-tailed t test). However, there was a significant difference in the tuning level of inhibition between these two cell classes (Fig. 4D; p < 0.001, two-tailed t test). In addition, for both orientation selective complex and simple cells, we found that inhibition and excitation had essentially identical orientation preferences (Fig. 4E), an observation in agreement with some measurements in the cat (Anderson et al., 2000; Mariño et al., 2005). Specifically, the excitation/inhibition (E/I) ratio at the orthogonal orientation was significantly larger in complex than simple cells (p < 0.001, two-tailed t test), whereas that at the preferred orientation did not differ between the two cell classes (p = 0.36, two-tailed t test; Fig. 4F).
Differential orientation tuning of excitatory and inhibitory inputs in complex and simple cells. A, Average normalized orientation tuning curves for excitatory (red) and inhibitory (blue) inputs to complex cells (n = 14). Error bar indicates SD. Inset, comparison of mean tuning widths (σ) of complex and simple cells. **p < 0.01, two-tailed paired t test. Error bar indicates SEM. B, Average normalized orientation tuning curves for excitatory and inhibitory inputs to simple cells (n = 12). Inset, Comparison of tuning widths. *p < 0.05, two-tailed paired t test. C, Ratio of OSI of inhibitory input over that of the excitatory input in the same cell plotted against OI measured from Vm responses. Blue dash line is the best-fit linear regression line (r = 0.82). Horizontal dash line marks OSI_inh = OSI_exc. Vertical dash line marks OI_Vm = 0.7, which separates complex from simple cells. D, Average OSIs for excitation and inhibition in complex and simple cell groups. ***p < 0.001, two-tailed paired t test (within group) or t test (between groups). Error bar indicates SEM. E, Preferred orientation of inhibition versus that of excitation. Red dash line is the identity line. Each data point represents one cell. F, E/I ratios at the preferred and orthogonal orientations. Blue symbols display mean ± SD. ***p < 0.001, two-tailed paired t test or t test. Data points for the same cell are connected by a line.
Narrowly tuned inhibition contributes to the weaker OS of complex cells
We have shown that a prominent difference between complex and simple cells lies in the relative tuning relationship between excitation and inhibition. In complex cells, E/I ratio was higher for the orthogonal orientation than the preferred orientation, while in simple cells E/I ratio was lower for the orthogonal orientation than the preferred (Fig. 4F). The relatively weaker inhibition at the orthogonal orientation may reduce the effectiveness of inhibitory sharpening of orientation tuning of output responses. Consistent with this idea, we found that the Vm response in complex cells showed significantly weaker orientation tuning than their excitatory input (OSI = 0.15 ± 0.08 for Vm vs 0.20 ± 0.11 for excitation, p < 0.01, two-tailed paired t test; Fig. 5A), whereas in simple cells OSI of Vm response was similar to that of excitation (OSI = 0.25 ± 0.10 for Vm vs 0.26 ± 0.09 for excitation, p = 0.62, two-tailed paired t test; Fig. 5A).
Inhibition can determine the differential orientation tuning in complex and simple cells. A, Comparison of OSIs of excitatory input and Vm response for complex and simple cells. Data points for the same cell are connected by a line. **p < 0.01, two-tailed paired t test. *p < 0.05, two-tailed t test. B, Comparison of peak evoked synaptic conductance at the preferred orientation between complex (C) and simple (S) cells. There is no significant difference between C and S for excitation or inhibition (p = 0.63, 0.21, respectively, two-tailed t test). C, Temporal profiles of simulated evoked excitatory (red) and inhibitory (blue) synaptic conductanc (top), and the derived Vm response using the neuron model (bottom). The resting membrane potential was set at −65 mV. D, Simulated orientation tuning curves for synaptic conductance. Black curve is for excitation (fixed at OSI = 0.25). Colored curves are for inhibition at various OSIs. E, Normalized tuning curves for derived Vm responses resulting from integration of excitation and inhibition of different tuning relationships (colors correspond to different conditions in D). F, OSIs of Vm (solid square) and spike (open triangle) responses at varying degrees of inhibitory tuning.
To have a more thorough understanding of how the relative tuning relationship between excitation and inhibition shapes OS of output responses, we performed simulations in a conductance-based neuron model (see Materials and Methods). Because we did not find significant differences in synaptic strengths between simple and complex cells at the preferred orientation (Fig. 5B), we fixed the amplitudes of excitation and inhibition at the preferred orientation, as well as the tuning profile of excitation, while allowing the tuning profile of inhibition to vary. We simulated excitatory and inhibitory conductance evoked by the preferred directional bar using mathematical functions fitting the average experimental data (Fig. 5C, top). By integrating excitation and inhibition in the neuron model, the expected Vm response could be derived (Fig. 5C, bottom). The amplitudes of excitation and inhibition at each of other stimulus orientations were than scaled according to their respective orientation tuning profiles. The OSI for excitation was set at 0.25 (Fig. 5D, black curve). The OSI for inhibition varied from 0.1 to 0.4 at a step of 0.1 (Fig. 5D, color curves). The normalized tuning curves of Vm responses derived from these simulated conductance showed that as the inhibitory tuning became more and more narrowed, the tuning of Vm response became broader and broader (Fig. 5E). As OSI of inhibition increased from 0 (i.e., untuned) to 0.4, OSI of Vm response reduced from 0.22 to 0.10 (Fig. 5F, solid square). This would result in an even faster decrease in OSI of spike response (Fig. 5F, open triangle), as demonstrated by our simulated spiking responses (see Materials and Methods). These modeling results indicate that inhibition being more narrowly tuned than excitation can largely contribute to the weaker OS of complex cells relative to simple cells.
Synaptic RFs of complex cells
Previous studies have shown that OS of simple cells can be predicted, at least partially, from their RF organization (Jones and Palmer, 1987; Gardner et al., 1999; Lampl et al., 2001; Liu et al., 2011). OS of complex cells however is generally thought not correlated with their RF architecture. In our loose-patch recording data, we noticed that the spike subfields were often slightly elongated, and the axis for the elongation was in accordance with the preferred orientation of the cell (Fig. 1M). Here, we further looked into this issue at the synaptic input level, by mapping 2-D on and off subfields of excitatory and inhibitory inputs (Fig. 6A,B). Consistent with the spiking response results, synaptic subfields (both excitatory and inhibitory) in complex cells appeared elongated (Fig. 6A). To measure the degree of elongation, we fit the synaptic subfields with ellipses (see Materials and Methods). An aspect ratio, defined as the ratio of the subfield long axis over the short axis, was calculated for the subfield of a given contrast (on or off). As shown by the plot of aspect ratio for inhibitory subfield versus that for its excitatory counterpart, all of the subfields had an aspect ratio well over 1 (Fig. 6C), indicating that they were indeed elongated. In addition, the majority of data points were above the identity line, indicating that in individual complex cells, the inhibitory subfield tended to be more elongated than the corresponding excitatory subfield (Fig. 6C; p < 0.01, two-tailed paired t test). Conversely, simple cells tended to show a more elongated excitatory than inhibitory subfield (Fig. 6D; p < 0.05, two-tailed paired t test). Comparison of distributions of aspect ratios (Fig. 6E) indicated that inhibitory subfields as a population were significantly more elongated in complex than simple cells (p < 0.001, two-tailed t test), whereas excitatory subfields were only slightly more elongated in complex than simple cells (p < 0.05, two-tailed t test). For individual cells, the ratio between aspect ratios of inhibitory and excitatory subfields was significantly >1 in complex cells (p < 0.001, two-tailed t test), and also was significantly higher in complex than simple cells (p < 0.001, two-tailed t test; Fig. 6F). These results demonstrate that in complex cells, the inhibitory subfield tended to be more elongated than the excitatory subfield responding to the same stimulus contrast, whereas in simple cells inhibitory subfield tended to be slightly less elongated. These observations correlate well with the finding that orientation tuning of inhibition is narrower than excitation in complex cells, and that it is narrower in complex than simple cells.
Geometry of synaptic subfields in complex and simple cells. A, Average excitatory and inhibitory responses to flash on and off square stimuli in an example complex cell. Square size is 5°. Ovals depict the elliptic fitting of the synaptic subfields. Scale, 300 pA (Exc)/560 pA (Inh), 300 ms. B, Synaptic RFs for an example simple cell. Scale, 250 pA (Exc)/510 pA (Inh), 300 ms. C, Aspect ratio for inhibitory subfield versus that for the excitatory subfield in the same cell and of the same contrast, plotted for individual complex cells (n = 38 pairs of subfields). Magenta symbol displays mean ± SD. Red dash line is the identity line. D, Aspect ratio for inhibitory subfield versus that for excitatory subfield, plotted for individual simple cells (n = 24 pairs of subfields). E, Distribution of aspect ratios for excitatory (top) and inhibitory (bottom) subfields in complex (black) and simple (gray) cell groups. Arrows indicate the mean values. F, Ratios of aspect ratio for inhibitory subfield over that for the excitatory subfield in the same cell and of the same contrast, in complex and simple cell groups. Blue symbols display mean ± SD. ***p < 0.001, two-tailed t test.
Discussion
In this study, we have measured the orientation tuning of complex cells in mostly the superficial layers of mouse V1 at three different levels: spiking response, membrane potential response, and synaptic conductance. We found that complex cells are less selectively tuned than simple cells at both the spiking and membrane potential response levels. By dissecting excitatory and inhibitory synaptic inputs with voltage-clamp recordings, we found that the difference in the degree of OS between complex and simple cells can be primarily attributed to an inhibitory synaptic mechanism: although the selectivity level of excitation is similar between complex and simple cells, inhibition is more selectively tuned in complex than simple cells. The two types of neuron thus differ in the tuning relationship between excitation and inhibition: in complex cells inhibition is more narrowly tuned than excitation, whereas in simple cells inhibition is more broadly tuned than excitation. Our modeling results further suggested that the differential excitatory/inhibitory tuning relationship itself can lead to differential degrees of OS.
Simple and complex cells
In cat V1, the thalamorecipient layers (layer 4 and upper layer 6) are dominated by simple cells, whereas regions beyond thalamocortical stages of processing such as superficial layers contain mainly complex cells (Gilbert, 1977; Martinez et al., 2005; Hirsch and Martinez, 2006). These laminar distributions of the cell types are consistent with the hierarchical model in which simple cells are constructed by convergent thalamic inputs, and complex cells are built by sampling simple-cell inputs. In mouse V1, we have previously shown that layer 4 is dominated by monocontrast cells (a special type of simple cell that displays only one spike subfield, either on or off; Liu et al., 2009; Y. T. Li et al., 2014), whereas layer 2/3 contains a large number of conventional simple and complex cells (defined by segregated and overlapping on and off spike subfields, respectively; Fig. 1F). At the synaptic level, the separation of on and off subfields of these cells exhibited a graded variation (Fig. 2G). These results suggest that simple and complex cells are not necessarily formed at different stages of hierarchical processing, but they can be formed at the same stage of processing within a common circuit (Alonso and Martinez, 1998; Priebe et al., 2004; Liu et al., 2010), likely by pooling monocontrast layer 4 inputs and intralaminar inputs.
It is generally thought that simple and complex cells are about equally selective for orientation. Although there are reports about differential orientation tuning of simple and complex cells, the difference is relatively small. For example, the median tuning bandwidths of simple and complex cells are 40° and 44° respectively in macaque V1 (De Valois et al., 1982). In a more recent study in macaque V1, no significant difference in tuning bandwidth has been found between simple and complex cells (Ringach et al., 2002). However, when OS is quantified with a global measure (circular variance), simple cells do exhibit higher OS than complex cells (Ringach et al., 2002). An earlier single-unit study in mouse V1 has reported that superficial neurons are highly selective for orientation without distinguishing between simple and complex cells (Niell and Stryker, 2008). In the current study, using a conventional index to measure OS, we found that complex cells are significantly less orientation tuned than simple cells. The functional significance of this differential level of OS is unclear. One possibility is that these V1 cells can provide functionally distinct inputs to different downstream areas, accounting for the functional specialization of these targets (Andermann et al., 2011; Marshel et al., 2011; Glickfeld et al., 2013).
Heterogeneous excitatory/inhibitory tuning relationship
Tuning properties of sensory cortical neurons are predominantly shaped by excitatory and inhibitory inputs they receive (Isaacson and Scanziani, 2011; Zhang et al., 2011). Our voltage-clamp recording data indicate that excitatory inputs to simple and complex cells do not differ in the degree of orientation tuning, which is consistent with the notion that these two types of cell may receive excitation from a common population of presynaptic neurons. However, they do differ in the tuning of inhibitory input, as well as the relative tuning relationship between excitation and inhibition. Our data in the mouse clearly indicate that excitation and inhibition have similar orientation preferences (Fig. 4E; Liu et al., 2011), although in the cat there is still debate as to whether this is the case (Anderson et al., 2000; Monier et al., 2003; Mariño et al., 2005). Under such circumstance, because inhibition is narrower than excitation in complex cells but broader than excitation in simple cells, if we average tuning curves from all the cells we would obtain a result suggesting that excitation and inhibition have similar tuning (Fig. 4C); a conclusion that was made by a previous study (Tan et al., 2011). Therefore, behind a big picture of cotuned or balanced excitation and inhibition (Mariño et al., 2005; Tan et al., 2011), the exact or detailed tuning relationship between excitation and inhibition can be in fact heterogeneous.
Inhibitory influence on tuning of output responses
Because complex and simple cells do not differ in excitatory input tuning, in synaptic amplitudes at the optimal stimulus orientation, or in intrinsic properties such as the effective spike threshold, their differential selectivity levels must be attributed mainly to the differential inhibitory tuning. Previously we have shown in simple cells that inhibition, while having a similar orientation preference as excitation, plays a critical role in sharpening OS. The sharpening is achieved by not only lowering the level of membrane potential responses in general, but also expanding the input dynamic range and rendering the input–output relationship more linear (Liu et al., 2011). Here, we found that in complex cells, inhibition also has the same orientation preference as excitation. However, inhibition is more sharply tuned than excitation. Such scenario determines that in complex cells, inhibition would have a much weaker or even no sharpening effect compared with simple cells. Indeed, in our neuron modeling when we manipulated the tuning level of inhibition from more broadly tuned to more sharply tuned, the tuning of resulting membrane potential responses became gradually broadened (Fig. 5F). Our current-clamp recording results confirmed that the tuning of subthreshold membrane potential response is broader in complex than simple cells (Fig. 2J), and that it is broader compared with excitatory input in complex cells (Fig. 5A). Recently, using dynamic clamp recordings, we also demonstrated in cortical cells that by broadening the sensory tuning of inhibitory input, the tuning of spiking response could be narrowed (Li et al., 2012a; L.Y. Li et al., 2014). In other words, the tuning of spiking response will be broadened if the tuning of inhibitory input is narrowed. Modulating the inhibitory tuning to shape OS seems a dominant mechanism that also accounts for the contrast invariance of orientation tuning of excitatory cells (Li et al., 2012b), and the developmental sharpening of orientation tuning of layer 4 excitatory neurons (Li et al., 2012a).
RF architecture and OS
A previous study has demonstrated that OS of simple cells can be predicted from the geometry of their RF maps at the membrane potential response level (Lampl et al., 2001). In principle, the spatial organization of synaptic RFs can also contribute to the orientation tuning of complex cells. Here, we found that excitatory subfields in complex cells are slightly but significantly more elongated in complex than simple cells (Fig. 6E). However, the on and off excitatory subfields in simple cells are spatially more separated compared with complex cells, and they are often displayed along the minor axis of the subfields. As demonstrated previously (Lien and Scanziani, 2013), such displacement of on and off excitatory inputs alone will result in orientation tuning. Together, the subfield shape and on/off displacement may render the summed excitatory input to simple cells to have a similar degree of orientation tuning as complex cells.
Potential origin of inhibitory input tuning
That inhibitory tuning is sharper in complex than simple cells can be explained by two possible mechanisms. First, the presynaptic neurons targeting a complex cell are more sharply tuned than those targeting a simple cell. Inhibitory neurons vary broadly in their tuning selectivity. For example, somatostatin-positive neurons are found to be much more sharply tuned than parvalbumin (PV) positive neurons (Ma et al., 2010). Even in the population of PV neurons which are generally considered broadly tuned (Kerlin et al., 2010; Ma et al., 2010; Hofer et al., 2011), some PV cells can be much more sharply tuned than the others (Runyan and Sur, 2013). In addition, other types of inhibitory neuron (Gonchar et al., 2007; Xu et al., 2010; Rudy et al., 2011) may also be selective for orientation, which remain to be investigated.
Second, even if presynaptic inhibitory neurons are all broadly tuned, an elongated spatial arrangement of their visual RFs can result in sharply tuned summed inhibitory input, in a similar manner as the convergence of untuned thalamic inputs on a simple cell in cat V1 (Reid and Alonso, 1995; Ferster et al., 1996; Alonso and Martinez, 1998; Lampl et al., 2001). That inhibitory subfields are more elongated in complex than simple cells (Fig. 6E) correlates well with the observation that the inhibitory tuning is sharper in complex than simple cells. Although the difference in inhibitory subfield geometry between complex and simple cells appears somewhat moderate, it can certainly contribute to a difference in inhibitory tuning strength. Together with our previous understanding that inhibitory on and off subfields in simple cells are overlapped similarly as in complex cells (Liu et al., 2010), our data suggest that the way inhibitory neuron RFs are spatially arranged is a likely mechanism underlying the more sharply tuned inhibition in complex cells.
Recent slice recording data demonstrate that inhibitory neurons in general connect with nearby excitatory cells rather densely and promiscuously (Fino and Yuste, 2011; Packer and Yuste, 2011; Pfeffer et al., 2013). These results cast doubts on whether the summed inhibitory input can be sharply tuned. However, that inhibition has the same orientation preference as excitation in complex cells argues against that inhibitory–excitatory neuron connections are entirely nonspecific (Yoshimura and Callaway, 2005; Ohki and Reid, 2007). These findings then raise new questions: why do complex cells have a specific linear arrangement of their inhibitory inputs? What kind of synaptic plasticity mechanism supports the formation of such circuit configuration? Future investigations are needed to answer these intriguing questions.
Footnotes
This work was supported by grants from the National Institutes of Health (EY019049) and Kirchgessner Foundation to H.W.T., and National Institutes of Health Grant DC008983 to L.I.Z.
The authors declare no competing financial interests.
- Correspondence should be addressed to Huizhong W. Tao, University of Southern California, 1501 San Pablo Street, ZNI 339, Los Angeles, CA 90033. htao{at}usc.edu