Acetylcholine Dynamically Controls Spatial Integration in Marmoset Primary Visual Cortex

M. J. Roberts, W. Zinke, K. Guo, R. Robertson, J. S. McDonald, A. Thiele


Recent in vitro studies have shown that acetylcholine (ACh) selectively reduces the efficacy of lateral cortical connections via a muscarinic mechanism, while boosting the efficacy of thalamocortical/feed-forward connections via a nicotinic mechanism. This suggests that high levels of ACh should reduce center-surround interactions of neurons in primary visual cortex, making cells more reliant on feed-forward information. In line with this hypothesis, we show that local iontophoretic application of ACh in primate primary visual cortex reduced the extent of spatial integration, assessed by recording a neurons' length tuning. When ACh was externally applied, neurons' preferred length shifted toward shorter bars, showing reduced impact of the extra-classical receptive field. We fitted a difference and a ratio of Gaussian model to these data to determine the underlying mechanisms of this dynamic change of spatial integration. These models assume overlapping summation and suppression areas with different widths and gains to be responsible for spatial integration and size tuning. ACh significantly reduced the extent of the summation area, but had no significant effect on the extent of the suppression area. In line with previous studies, we also show that applying ACh enhanced the response in the majority of cells, especially in the later (sustained) part of the response. These findings are similar to effects of attention on neuronal activity. The natural release of ACh is strongly linked with states of arousal and attention. Our results may therefore be relevant to the neurobiological mechanism of attention.


Visual cortical neurons respond to stimuli appearing within a restricted area of visual space known as the cell's classical receptive field (CRF); however, visual perception relies on integrating information from across the visual field. Neurons gather information via intracortical connections from distant cells within a cortical area (Brown et al. 2003; Das and Gilbert 1999; Stettler et al. 2002) and from higher areas (Bullier et al. 2001; Moore and Armstrong 2003). These connections provide a neuron with information from a larger region of visual space, described as the cell's “nonclassical receptive field” (nCRF). Stimuli within this region do not elicit action potentials, but can modulate responses to stimuli appearing in the cell's CRF (Albright and Stoner 2002; Born 2000; Cavanaugh et al. 2002; Das and Gilbert 1999; Jones et al. 2002; Kapadia et al. 1999 2000; Knierim and van Essen 1992; Zipser et al. 1996). The nCRF can include a facilitatory inner ring, functionally continuous with the CRF, but that is unable to generate spikes, although it can change the subthreshold membrane potential (Anderson et al. 2001). The much larger part of the nCRF is suppressive in nature (DeAngelis et al. 1992, 1994). A significant part of the CRF response is due to thalamocortical/feed-forward synapses (Angelucci et al. 2002). These synapses carry information into a cortical area from the thalamus or from areas lower in cortical hierarchy to areas higher in hierarchy level (Van Essen et al. 1992). Intracortical (and feedback) synapses recombine information within the cortex and mediate interactions across visual space, they provide contextual nCRF modulation (Das and Gilbert 1999), and they can be facilitatory or suppressive in nature. Contextual modulation may be instrumental in assigning “meaning” to a visual scene (Albright and Stoner 2002). In cluttered environments, contextual information (the clutter) can hinder the detection of small objects, and under natural vision, it may thus be necessary to dynamically adjust the flow of feed-forward and lateral/feedback information; for example, to preferentially process information from a small area of visual space that is of behavioral relevance (Chelazzi 1995; Luck et al. 1997; Reynolds and Desimone 1999). Such dynamical shifting could be done by selectively controlling the synaptic efficacy of feed-forward and intracortical synapses. Recent in vitro studies suggest that acetylcholine (ACh) might be involved in this dynamic shift (Hasselmo 1995; Hasselmo and Bower 1992, 1993; Kimura 2000; Kimura et al. 1999; Linster and Hasselmo 2001). ACh selectively suppresses the efficacy of intracortical synapses by activating muscarinic receptors (Kimura and Baughman 1997). By activating nicotinic receptors, primarily located on thalamocortical/feed-forward synapses (Lavine et al. 1997; Prusky et al. 1987; Sahin et al. 1992), ACh boosts the efficacy of the feed-forward/thalamocortical input that provides information about the CRF (Gil et al. 1997). ACh applied in vivo should thus result in reduced impact of stimuli presented in the nCRF while increasing the effect of stimuli placed within the CRF.

To test this proposal, we investigated length tuning in V1 under conditions of externally applied and not externally applied ACh. Most neurons in V1 respond maximally to bars that extend beyond their CRF (DeAngelis et al. 1994; Sceniak et al. 2001), showing spatial summation from the nCRF. We show that application of ACh caused a shift of a neuron's preferred length toward shorter bars and a decrease in its summation area, supporting the hypothesis that ACh re-balances lateral/feedback and feed-forward connections in favor of feed-forward activation.


All experiments were carried out in accordance with the European Communities Council Directive 1986 (86/609/EEC), the National Institutes of Health Guidelines for Care and Use of Animals for Experimental Procedures, the Society for Neurosciences Policies on the Use of Animals and Humans in Neuroscience Research, and the UK Animals Scientific Procedures Act.

Electrophysiological recordings

We recorded extracellular responses of V1 neurons from four adult anesthetized and paralyzed marmosets (Callithrix jacchus, 400–480 g). We attempted to record exclusively from single neurons; however, occasionally the isolation may have been such that multiunit activity (2–3 cells) was recorded. Based on assessment of the autocorrelation, we recorded from 56 single units and 10 multiunits, as evidenced by the absence, or small numbers, of spikes in the bins 1–3 ms following the trigger spike. Anesthesia was induced by intramuscular injection of Saffan (alphadalone/alphaxalone acetate, 1.5 ml/kg) and maintained by continuous intravenous injection of Propofol (Diprivan; 0.8–1.5 ml/kg/h). Analgesia was ensured by continuous injection of Alfentanil (156 μg/kg/h). Paralysis was induced and maintained by intravenous injection of vancuronium (Norcuron, 100 μg/kg/h). Animals were artificially ventilated at a rate of 30–70 strokes/min (3.5–5.5 ml/stroke). End-tidal CO2 was constantly monitored and maintained between 3.5 and 4.5%. In addition, arterial and venous blood pressure and electrocardiogram were continuously monitored and recorded. Prior to paralysis, adequate depth of anesthesia was ensured by repeatedly checking for absence of toe pinch withdrawal reflexes. Level of anesthesia following paralysis was monitored by means of heart rate and/or blood pressure changes following toe pinches. Animals received antibiotic injections every 12 h (Cephuroxide, 125 mg/kg). Eyes were protected with contact lenses and regularly irrigated with saline. Atropine eye drops were regularly applied to induce and maintain mydriasis and cycloplegia.

ACh was applied iontophoretically via a barrel pipette onto which the recording electrode was mounted. The distance between the ACh pipette (5BBL W/FIL 1.2 mm, World Precision Instruments) and electrode (FHC, 1–2 MΩ) tips was 25–50 μm. Pipette impedance was 10–30 MΩ. ACh concentration was 0.8 M (pH 4.5). We applied retention currents of –10 to –5 nA, and ejection currents varied between +1 and +100 nA. We generally tried to adjust the ejection current depending on the strength of the effects of application. To do so we started with a relatively low ejection current (10–20 nA) and monitored the effect on neuronal activity over time, while visual stimuli were presented. If the ejection current did not result in activity changes, we increased the current to 50–60 nA and repeated the procedure. If no results were obtained by this application current, we further increased the current to 80–100 nA. If no ACh effects were obtained using this current strength, we advanced our electrode to the next cell. If no effects were recorded for three consecutive cells, we retracted and replaced the electrode/pipette. In some occasions, ejection currents of 10 nA caused enormous ACh effects, such that the cell increased firing rates dramatically or ceased to fire entirely. If such behavior was encountered, we reduced the ejection current to 1–5 nA. We repeatedly ensured that neuronal activity changes were not due to the currents applied by keeping the overall current constant with the aid of compensation pipettes filled with 0.9% saline. To avoid ACh getting sucked into the compensation pipette under these circumstances and being ejected during the retention phase, we set the compensation currents such that the overall current flow was identical during ACh retention and ejection, while at the same time a positive current was always applied to the saline pipette.

Stimuli and protocol

Stimuli were displayed on a 20-in analogue CRT monitor (75 Hz, 1,600 × 1,200 pixels) positioned 57 cm from the animal. They were presented on a gray background (24.6 cd/m2). Stimuli were brighter or darker than the background, depending on cell preference (70% Michelson contrast). The RF borders were mapped by moving a bar of adjustable size and orientation across the screen [minimum response field (mRF)] (Barlow et al. 1967; Maffei and Fiorentini 1976). mRF locations were within the central 10° for all neurons reported herein. mRF diameter ranged from 0.3 to 1.5°. After determining the cell's preferred orientation (at a resolution of 22.5°), bars of varying length at the preferred orientation were presented centered over the mRF. Bar length was adjusted in seven steps by multiples of the mRF diameter (0.5, 0.75, 1, 1.5, 2, 3, and 5 times mRF diameter in animals 1 and 2; 0.2, 0.4, 0.8, 1.6, 3.2, 6.4, and 12.8 times mRF diameter in animals 3 and 4). Bar width was fixed at 0.15° for RFs > 0.75°, and 0.05° for RFs ≤ 0.75° diameter. To prevent adaptation of neural responses due to presentation of just one orientation, an additional four conditions were included in which bars orthogonal to the preferred orientation were presented (1.5, 2, 3, and 5 times mRF diameter in animals 1 and 2 and 1.6, 3.2, 6.4, 12.8 times mRF diameter in animals 3 and 4). We tested whether adaptation occurred by calculating a linear regression for each condition and recording (separately for initial, application, and control recording). We did find a very small trend for negative regression slopes, i.e., decreasing activity as the experiment progressed (median slope: −0.0623, 25 percentile: −0.204, 75 percentile: 0.144). This means the median firing rate decreased very little over consecutive trial (<2% over 10 trials if the starting firing rate was 50 spikes/s), and if this decrease was due to adaptation, it was unlikely to have influenced our general conclusions, because it was present for the initial recording, application and control, while there was no consistent decrease from the initial recording to the control recording. Additionally we inserted 3- to 5-min waiting times between recordings (initial recording, ACh application, recovery), which are also likely to be sufficient to eliminate adaptation between recordings. Stimuli were presented interleaved ≥15 times each. The presentation time of the stimuli was 500 ms with 500-ms pre- and 200-ms poststimulus time. Stimuli were presented and spike timings were collected with a sampling resolution of 1 ms under the control of Remote Cortex 5.95 (Laboratory of Neuropsychology, NIMH, Bethesda, MD). Recordings were performed during monocular stimulation. Spontaneous activity was calculated from the prestimulus time separately for ACh applied and not applied. All stimulus-driven activity presented here was corrected for spontaneous activity.

Length tuning was initially measured with no external ACh applied (control condition), at least once with ACh applied, and subsequently, in at least one repeated control condition. Neuronal activity was compared across control (ACh not applied) conditions to ensure that full recovery following ACh application occurred (P > 0.05, 2-way Kruskal Wallis ANOVA). Cells that showed significant differences across control conditions (P < 0.05) were excluded from further analysis. For the remaining cells, data from control conditions were combined and compared with data obtained during ACh application. Only cells showing a significant effect of ACh during either spontaneous or stimulus driven activity were included in further analysis. For a cell to be included, it was sufficient that either the spontaneous activity was significantly affected by ACh or the response to presentation of any of the preferred orientation bars was significantly affected by ACh application (P < 0.05, 2-way Kruskal-Wallis ANOVA).

Difference of Gaussian model

Length tuning data were fitted with a difference of Gaussian model (DOG model) (Sceniak et al. 2001). In this model, the narrower Gaussian represents the RF excitatory center, while the broader Gaussian represents the inhibitory surround. Each Gaussian is described by a strength (gain) and a space constant, determining its height and width, respectively. This function captures the shape of measured length tuning curves and it allows the relative contribution and size of excitation (summation) and inhibition (suppression) areas to be separated. The fitted function is of the form Math(1) Where R corresponds to the model's response to a bar of length y, Ke corresponds to the excitatory component amplitude (the summation gain), a corresponds to the size constant of excitatory area (referred to as “summation area”), Ki corresponds to the inhibitory component amplitude (the suppression gain), and b is the size constant of the inhibitory area (referred to as “suppressive area”). Fits of the summation area and suppression area were constrained such that the suppression area was larger than the summation area, and we chose the maximum size to be 30 times mRF. We have tested different constraints, i.e., allowing these areas to become substantially larger (no constraints, size diameters 12.8 and 20 times mRF) with the same general outcome as described in results.

Ratio of Gaussian model

An alternative description of spatial integration can be given by a ratio of Gaussian Model (ROG) (Cavanaugh et al. 2002). It is principally similar to the difference of Gaussian model, but assumes that the influence of the suppression area is a normalization Math(2) Again, R corresponds to the model's response to a bar of length y, Ke corresponds to the excitatory component amplitude, a corresponds to the size constant of excitatory area, Ki corresponds to the inhibitory component amplitude, and b is the size constant of the inhibitory area.

Fit optimization

Fits were optimized to minimize the χ2 error (Press et al. 2002). In short, we searched for the model (m) parameters (a) that minimized the error function Math(3) where r is the response to stimulus s, and var(r) is the response variance. Since response variance is generally proportional to the firing rate (Carandini et al. 1997), we separately estimated the variance fitting a polynomial function Math(4) to our variance versus mean firing rate data that were obtained with ACh applied and not applied. The estimates for the variables a and b were different for the two conditions (ACh not applied: var = 5.9271 × mean1.1403; ACH applied: var = 14.0572 × mean0.9317). To avoid giving too much importance to data points with low firing rate and variance, we set all variance values <1 to be equal to 1. To increase the probability that our fitting routine yielded small error values (and thus good fits), we initially fitted our data with a set of 24 different starting positions for the different parameters. The starting parameters that resulted in the smallest χ2 errors were used for the final optimization. Empirical evidence showed that starting parameters needed to be different for the DOG and the ROG models to produce adequate fits with small χ2 errors, and were thus different for the two fitting functions.

We assessed the goodness of each fit by calculating the normalized χ2 error (χ2N) between that data and the model predictions (Cavanaugh et al. 2002) and calculated the associated P value (Press et al. 2002) to determine whether fits were acceptable. We only included cells where P values were >0.05. Additionally we calculated the percentage of variance accounted for by the model (Carandini et al. 1997). The preferred length of a neuron corresponds to the length at which the peak of the fitted function occurred.

Determination of significance of differences for individual cells

We applied a bootstrap procedure to determine the reliability of the effects of length tuning shifts, changes in gain, and changes of summation/inhibition areas. To this end, we selected for each stimulus condition a set of 15–45 trials (depending on the number of repetitions measured for the respective cell) at random (random with replacement) and performed the model fitting based on these selections of trials. This bootstrapping procedure was performed 100 times for each cell when ACh was applied and when it was not applied, thus resulting in 100 different preferred lengths, gain, and summation/inhibition area estimates with ACh applied and 100 estimates with ACh not applied for each cell recorded. We used a signed-rank test to determine whether these estimates were significantly different for the two conditions.

Tonic index

To determine whether ACh changed the response profile of a cell, we calculated the tonic index (TI) as the firing rate during late response period (Rlate, 250–500 ms after stimulus onset), divided by the firing rate during the early response period (Rearly, 30–250 ms after stimulus onset) Math(5) This index determined whether cells responded more or less tonically (sustained) in the presence of externally applied ACh.


We recorded the length tuning of a total of 120 neurons. Of these, 66 units showed a significant effect of externally applied ACh on firing rate (P < 0.05, 2-way Kruskal Wallis ANOVA) and a return to baseline following recovery. How many of the remaining 54 neurons were unaffected by ACh is impossible to say, because the pipette may either have leaked occasionally, causing high ambient ACh level throughout individual recordings, it may have been blocked despite current flow across the tip, or neurons may have shown long-term ACh effects that cannot be distinguished from drifts in firing rate. Due to these difficulties, we did not attempt to quantify the number of neurons that were susceptible to ACh application. Although we cannot be sure whether endogenous ACh was released when ACh was not applied or small amounts leaked from the pipette, we will refer to the condition of externally applied ACh as “ACh present” and to the condition when ACh was not externally applied as “ACh absent” throughout the remainder of the text. Presence of ACh caused a facilitation of responses in 41/66 (62.1%) units and suppression in 25/66 units (37.9%) compared with absence of ACh. This number is comparable with reports from earlier studies (Metherate et al. 1988; Murphy and Sillito 1991; Sato et al. 1987; Sillito and Kemp 1983; Sillito and Murphy 1987).

To determine the preferred length, the summation area, the inhibition area, and the excitatory and inhibitory gains we fitted our data with a DOG and with a ROG model (see methods). These models capture length tuning properties and provide independent estimates of the relative strength and size of the summation and suppression areas (Cavanaugh et al. 2002; Sceniak et al. 2001). Both models provided good fits to the data (the median variance accounted for was >90% for all conditions). The ROG model yielded marginally better fits as evident from the variance accounted for and from the χ2 values. The median variance across stimuli for which the DOG model accounted was 90.3% (25–75 percentile: 78.2–96.7%) for data fits with ACh absent, and it was 92.7% (25–75 percentile: 76.7–98.0%) for data fits with ACh present. For the ROG model, these values were 92.05% (25–75 percentile: 78.4–96.7%) and 93.49% (25–75 percentile: 81.2–98.2%), respectively. The median χ2 value for the DOG model with ACh absent was 0.136 (25–75 percentile: 0.061–0.234) and with ACh present it was 0.080 (25–75 percentile: 0.034–0.156). For the ROG model, the median χ2 value with ACh absent was 0.117 (25–75 percentile: 0.053–0.212) and with ACh present it was 0.072 (25–75 percentile: 0.024–0.153). Based on these marginal differences between the quality of fits for the DOG and ROG model, we do not attempt to determine which of the two is a better descriptor of V1 receptive fields. The important finding is that the two models yielded similar results regarding the changes of spatial integration under conditions of ACh absent and present as detailed below.

Length tuning

For the majority of cells (69.7%), the highest activity with ACh absent occurred at lengths greater than the diameter of the mRF, showing that the area surrounding the mRF facilitated the response to long bars of the preferred orientation. With ACh present, the preferred length tended to shift toward shorter bars (Figs. 1 and 2, left column; see also Table 1), suggesting that ACh reduced modulation from outside the mRF. The preferred length with ACh absent and present was determined from the curve fitted to the data (corresponding to the peak location). With ACh present, the preferred length was significantly reduced across the population [Fig. 2, left column, P < 0.001, signed-rank test (SRT); Table 1 provides more details regarding these changes]. The effect was found both for cells that were facilitated (P = 0.008, SRT) as well as cells that were suppressed (P = 0.011, SRT) by ACh. As we find a reduction of preferred length in facilitated and suppressed cells, the results cannot be due to response saturation. The reduction of preferred length was largely due to a decrease of the summation area (Fig. 2, column B; Table 1). The suppression area (Fig. 2C; Table 1) was not significantly affected by ACh application. The strength of the excitatory gain was significantly increased in the DOG model (P = 0.014, Fig. 2D, top row), but not in the ROG model (P > 0.05). The strength of the inhibitory gain was not significantly affected (P > 0.05; Fig. 2E; Table 1). Thus the consistent finding with both models was a decreased preferred length, and a decreased summation area with ACh present.

FIG. 1.

Single cell examples. Mean stimulus driven response to seven stimuli of varying length in the presence (black) and absence (gray) of acetylcholine (ACh). Smooth lines show fitted difference of Gaussians (DOG) models. Vertical arrows mark preferred length, taken as the peak of the fitted curve. Preferred length was shortened with ACh present. Error bars show SE. Histograms at the base of graphs show distribution of preferred length determined by the bootstrap method with ACh present (black histogram) and absent (gray histogram). x-axis: bar length in multiples of classical receptive field (minimum response field mapping).

FIG. 2.

Preferred length and fitting parameters extracted from DOG (top row) and the ratio of Gaussian model (bottom row). A: preferred lengths of neurons (n = 66 cells) with ACH absent (x-axes) and ACh present (y-axes). Cells that shifted their tuning toward shorter bars during ACh application appear below the diagonal. B: spatial summation area with ACh absent and present. Summation area decreased in the presence of ACh. C: suppression area with ACh absent and present. Width of suppression area was not affected by ACh application D: gain of summation area with ACh absent and present. E: gain of suppression area with ACh absent and present. Black dots, neurons for which the parameter of interest changed significantly as a function of ACh application; open squares, neurons for which no significant change was found (bootstrap method); symbols (dots, squares), horizontal, and vertical bars, median and 25–75 CIs for the parameter of interest (bootstrap method).

View this table:

Median and 25/75 percentiles for the preferred length (peak in multiples of mRF diameter), excitatory gain (Ke, arbitrary units), inhibitory gain (Ki, arbitrary units), summation area (in multiples of mRF diameter), and inhibition area (in multiples of mRF diameter) as a function of ACh application for the difference and the ratio of Gaussian model

In a recent paper, Cavanaugh et al. (2002) argued that stimulation-induced changes of spatial integration were not due to changes in the size of the summation or suppression area, but rather due to changes in the gain of these mechanisms. They employed the ROG model. To determine whether our data could be described by changes in response gain alone (gain model; Cavanaugh et al. 2002), we fitted the data such that the space constants were forced to take the same values with ACh absent and present, while excitatory gain and inhibitory gains could change. Not surprisingly, these fits were worse, but after accounting for the difference in degrees of freedom resulted in similar χ2N (median ROG Gain model χ2N: 0.035, 25 percentile: 0.0197, 75 percentile: 0.0686, median ROG model χ2N: 0.038, 25 percentile: 0.0198, 75 percentile: 0.0689). For the DOG model, the effects were similar (median DOG Gain model χ2N: 0.044, 25 percentile: 0.0237, 75 percentile: 0.0945, median DOG model χ2N: 0.038, 25 percentile: 0.0213, 75 percentile: 0.0761). These fits indicate that the changes in length tuning we observed during ACh application are unlikely to have occurred by changes in gain alone, particularly because the change in summation area showed the most consistent effect (evident by smaller P values and by changing into the same direction; Table 1). We take this as evidence that the major effect of ACh application was to change the size of the summation area, rather than (or in addition to) the gains of facilitation or suppression.

Effect of ACh on the time course and evolution of the response

Figure 3 shows the population response for facilitated (Fig. 3A) and suppressed (Fig. 3B) cells as a function of bar length and ACh application. It shows the data for cells measured at bar lengths ranging from 0.2 to 12.8 times mRF diameter (data for cells measured at 0.5–5 times mRF diameter gave similar results). With ACh present, cells facilitated by ACh showed higher activities mostly during the later part of the response (Fig. 3A) compared with ACh absent, while cells suppressed by ACh showed a reduced response when ACh was present mostly in the initial response period (Fig. 3B). We calculated the preferred population length as a function of time from response onset by fitting the DOG model to the population responses for both cell groups. We took the population response that occurred in 5-ms bins starting 35 ms after stimulus onset (smoothed with a half Gaussian, σ = 10 ms) for the different bar lengths as input data for the fitting routine. Figure 4 shows the preferred length and the spatial summation area as a function of time and ACh present/absent for facilitated (Fig. 4A) and suppressed cells (Fig. 4B). Cells facilitated by ACh preferred shorter bars from shortly after stimulus onset; however, the difference in preferred length became particularly pronounced from ∼150 ms after stimulus onset. These differences seem to occur somewhat later for cells that were suppressed by ACh application. The summation area showed a very similar behavior to the preferred length for both cell groups, while the suppression area was not systematically affected by ACh application (data not shown). It may seem puzzling that facilitated cells preferred shorter bars from ∼150 ms after response onset, although the ACh induced changes occurred only after the initial peak, while suppressed cells showed the effects later, where these suppressed cells showed ACh induced suppression mostly (but not exclusively) during the early response phase. From Fig. 3, it is apparent that this suppression was strongest at intermediate bar lengths, and these bar lengths also elicited the strongest responses during the early response phase in the absence of ACh. Thus the ACh-induced reduction during the early response phase is to some extend scaled with the response strengths itself, and therefore has resulted in similar preferred bar lengths. This was different for the late response phase. Here the ACh induced reduction was almost the same for medium and long bars, while there was no decrease (or even a small increase) in firing rate for small bars. This resulted therefore in a shift of preferred length toward shorted bars during the late response phase for cells suppressed by ACh.

FIG. 3.

Population activity as a function of bar length and ACh application. Population activity for cells facilitated with and without spontaneous activity subtracted (left 2 columns) and suppressed (right columns) by ACh with ACh absent (solid gray line) and ACh present (dashed black line) for bar length 0.2–12.8 times the minimum response field (mRF) diameter (top to bottom). Neurons facilitated by ACh showed increased activity from response onset particularly when bars were shorter than the neuron's mRF diameter (0.2 and 0.4 mRF diam), but the largest differences occurred during later response periods for all bar lengths. Effect of ACh on suppressed neurons was strongest when bars larger than the mRF were presented.

FIG. 4.

Preferred length and summation area as a function of time for cells facilitated by ACh and cells suppressed by ACh. Preferred length and summation area with ACh absent (gray line, gray area shows SE) and present (dashed black line, dashed flanking lines show SE) for cells facilitated (left) and suppressed (right) by ACh. With ACh present, preferred length and summation area were reduced from ∼150–200 ms after stimulus onset in cells facilitated by ACh, whereas this difference occurred somewhat later (230–280 ms) for cells suppressed by ACh, and it did not become significant until ∼320 ms after response onset in the latter cell group. Histograms at the base show time-resolved P values (1/P value) of differences between the ACh present/absent condition.

Interestingly, preferred length and summation areas show quite a degree of temporal dynamics when ACh was present and when it was absent. Whereas preferred length and summation area were fairly large during the initial response period, they decreased substantially during the period of ∼130–200 ms after stimulus onset and increased again thereafter.

Response profile

Responses with ACh present were more sustained, i.e., the difference in firing strength between the transient and sustained part of the response was decreased. This occurred for cells whose stimulus driven response was increased by ACh (Fig. 3A), as well as those whose response was decreased by ACh (Fig. 3B). However, there were notable differences between these changes. While cells facilitated by ACh showed an increase mostly during the sustained response, cells suppressed by ACh mostly showed a reduction of the transient response. We quantified the effect by calculating the TI (see methods): the average TI with ACh present and absent as a function of bar length is plotted in Fig. 5. ACh presence significantly increased TI, i.e., cell responses became more sustained. This effect was significant for cells facilitated by ACh (cells measured at 0.2–12.8 mRF diameter: P < 0.001, n = 19; cells measured at 0.5–5 times mRF: P < 0.001, n = 22; 2 factor ANOVA, factor 1: treatment, factor 2: bar length), but the trend also occurred for cells suppressed by ACh, although it did not reach significance for cells measured with bar lengths 0.5–5 times mRF, probably due to the small sample size (cells measured at 0.2–12.8 mRF diameter: P = 0.007, n = 18; cells measured at 0.5–5 times mRF: P = 0.086, n = 7; 2 factor ANOVA). The response profile also seemed to depend somewhat on bar length, although bar length affected TI to a lesser extent than ACh application. TI was not significantly affected by bar length for cells facilitated by ACh that were measured at 0.2–12.8 mRF diameter (P = 0.159, 2-factor ANOVA), while TI was significantly affected by bar length for cells measured at 0.5–5 times mRF diameter (P < 0.001, 2-factor ANOVA). Bar length also appeared to affect TI for cells suppressed by ACh (cells measured at 0.2–12.8 mRF diameter: P = 0.008; cells measured at 0.5–5 times mRF: P = 0.059; 2-factor ANOVA). There was no significant interaction between ACh present/absent and bar length for any of these measurements, i.e., the effects of ACh presence/absence on response profile did not depend on what length was present. These results show that ACh application changed a cells' response to be more sustained. Increasing the bar length had a similar, but less profound and less consistent, effect.

FIG. 5.

Strength of sustained response as a function of ACh. Tonic index for cells measured with bar length ranging from 0.2 to 12.8 times mRF (top) and measured with bar length ranging from 0.5 to 5 times mRF (bottom) for cells facilitated by ACh (left) and cells suppressed by ACh (right). Sustained response was on average increased in the presence of ACh (black dashed curves and error bars) compared with its absence (gray solid curves and error bars). Error bars are SE.

Location of recording sites

We did not make electrolytic lesions at the end of each recording track, because that would have destroyed parts of the intrinsic V1 network, which would likely have affected our results regarding the effect of ACh on intra-areal processing. However, we took care to monitor our recording depth precisely. After making a small incision into the dura prior to each track, we positioned our electrode/pipette under microscope guidance such that our zero depth registration corresponded to the location where the pipette tip just touched the cortical surface. We attempted to make penetrations perpendicular to the cortical surface, thereby hoping to be able to reconstruct the depth (and potentially the layers cells were recorded from) with reasonable precision. From these measurements, we reconstructed the recording depth. There was no obvious correlation between recording depth and whether cells were facilitated or suppressed by ACh. Neither was there an obvious correlation between recording depth and an ACh induced decrease of a cell's summation area (and thus length preference).

The above data were obtained in marmoset monkeys (C. jacchus). We also recorded 19 cells with a significant effect of ACh in one anesthetized macaque monkey (Macaca mulatta) under otherwise identical conditions. These cells showed identical trends, i.e., they showed a systematic reduction in preferred length and a reduction of the summation area.


The transmitter ACh has been implicated in a variety of functions, ranging from improved sensory processing, learning, arousal, attention, and even awareness (Davidson et al. 1999; Everitt and Robbins 1997; Robbins 1997; Sarter and Bruno 1997; Sarter et al. 2001). Moreover, deficits in cortical cholinergic functions often cause substantial cognitive deficits (Perry et al. 1999; Sarter and Bruno 1998). Despite its wide implication in various functions, the precise effects of ACh on cortical processing remain unknown. Several recent in vitro studies suggest that a key function of cortical ACh may be to control the flow of neuronal information by selectively suppressing lateral intracortical synapses while leaving thalamocortical/feed-forward synapses unaffected (Hasselmo and Bower 1992; Kimura et al. 1999) or even increasing their efficacy (Gil et al. 1997; Hsieh et al. 2000). In the visual cortex, these two types of synapse separately influence nCRF modulation and CRF activation. Thus high levels of ACh should attenuate nCRF interactions while potentially facilitating the CRF response. In line with this proposal, we show that, during ACh application, cells shift their length preference toward shorter bars, showing reduced summation from outside the classical receptive field. Fitting the length tuning data with a DOG or ROG model showed that ACh reduced the size of the summation area. Despite reduced spatial summation, most cells responded more strongly in the presence of ACh, suggesting that CRF activation was also boosted. We did not find significant changes in the size of the suppression area under conditions of increased ACh. This potentially reflects the fact that inhibitory synapses are less affected by ACh than excitatory synapses (Kimura and Baughman 1997).

In our study, we have defined the spatial diameter of the receptive field as the mRF (Barlow et al. 1967; Blasdel and Fitzpatrick 1984); the area within which presentation of a small bar elicits an extracellular response. The term CRF is often used as a synonym (Knierim and van Essen 1992). Surrounding the CRF/mRF is an area that can modulate the response to stimuli presented within the CRF/mRF, and this modulation can be facilitatory or inhibitory. The excitatory and inhibitory parts of the receptive field are usually assumed to extend over the CRF and the nCRF. The excitatory area is continuous in the sense that within its spatial sampling range inputs elicit excitatory postsynaptic potentials (EPSPs), but the number of EPSPs elicited and/or efficacy of changing the membrane potential near the axon hillock decreases as the stimulus is moved away from the center of the receptive field (Orban et al. 1979b; Sillito 1977). The CRF can be well modeled with a Gaussian envelope (Jones and Palmer 1987). According to the DOG model (DeAngelis et al. 1992), the CRF is surrounded by an excitatory fringe, which is continuous with the CRF and contributes to spatial summation (Cavanaugh et al. 2002; Sceniak et al. 1999). DeAngelis et al. (1992) argue that a cell's excitatory receptive field is most appropriately “described as the region within which a stimulus can either elicit an excitatory response or add to the response elicited by another stimulus.” By having used the mRF technique to determine the RF size, we are likely to have only elicited subthreshold responses in the insensitive excitatory RF fringes, and we therefore saw changes in suprathreshold responses when presenting longer bars. As a result, the preferred bar length was almost always larger than the mRF (the median preferred length in the absence of ACh was 1.33 time the mRF size). The suppressive area seems to be organized slightly differently and subdivided into two systems. One is acting in the RF center and is hardly orientation tuned (DeAngelis et al. 1992), the other is responsible for end and side inhibition, is orientation tuned, and originates from outside the excitatory RF, although it can overlap with it (DeAngelis et al. 1994; Orban et al. 1979a). End and side inhibition seem to be mediated through intracortical inhibitory interactions between binocular neurons (DeAngelis et al. 1994). It is often assumed that the long horizontal axonal projections within V1, particularly within layers 2 and 3, mediate the influences (excitation and suppression) from beyond the CRF (Gilbert et al. 1996; Hupe et al. 2001a; Stettler et al. 2002), whereas others have argued that modulation from outside the CRF is due to feedback projections from higher areas (Bair et al. 2003; Knierim and van Essen 1992; Lamme 1995; but see also Hupe et al. 2001a). Despite these arguments, the exact source of the signals mediating CRF and nCRF contributions is not well known. The excitatory area is likely to receive a mixture of feed-forward connections from the thalamus (Angelucci et al. 2002), vertical connections from within V1 (Callaway 1996), lateral connections from within area V1 (Gilbert et al. 1996; Hupe et al. 2001a; Stettler et al. 2002), and feedback projections from higher areas (Bullier et al. 2001; Hupe et al. 1998, 2001b). Inhibitory surround mechanisms are likely to originate to some extend within V1 (Gilbert et al. 1996; Hupe et al. 2001a; Stettler et al. 2002), but also from feedback projections (Bair et al. 2003). Although the source of inputs to the excitatory fringe surrounding the CRF is not precisely known, our finding of a reduced influence of this fringe (a reduced summation area) when ACh was applied suggests that the inputs to this fringe are predominantly intracortical synapses, which are suppressed by a muscarinic mechanism on ACh application (Kimura and Baughman 1997). Interestingly, we did not find an influence of ACh on the size or gain of the inhibitory area. If the surround suppression was indeed largely mediated by feedback projections (Bair et al. 2003), this could suggest that feedback connections are affected by ACh in a different manner than lateral connections. Future intracellular in vivo studies are necessary to determine the sources of these interactions and how they are affected by neuromodulators.

Previous experiments have shown that the relative strength and size of the summation and suppression areas can change according to stimulus parameters and context (Kapadia et al. 1999, 2000; Sceniak et al. 1999, 2001). The mechanism behind this change is still debated. While Sceniak et al. (2001) suggest that different stimulus configurations can cause a change in the width of the summation area, Cavanaugh et al. (2002) argue that these changes are better explained by gain changes alone. Contrary to the reports of Cavanaugh et al., our data were better explained by a change in spatial summation area, rather than a gain change. This is probably due to different mechanisms involved in the two phenomena. While changes in stimulus contrast are likely to involve contrast normalization (Heeger 1992; Heeger et al. 1996; Simoncelli and Heeger 1998), possibly mediated by GABAergic mechanisms (Heeger 1992; Heeger et al. 1996; Simoncelli and Heeger 1998; Thiele et al. 2004) or synaptic depression (Abbott et al. 1997; Carandini et al. 2002), ACh acts through a variety of presynaptic and postsynaptic mechanisms and receptors (Alkondon et al. 2000; Gebhard et al. 1993; Gil et al. 1997; Iannazzo and Majewski 2000; Kimura 2000; Kimura and Baughman 1997; Kimura et al. 1999; Levin and Simon 1998; McCormick and Prince 1985). Thus the different findings are not necessarily contradictory.

A recent paper (Ozeki et al. 2004) showed that applying the GABAa receptor antagonist, bicuculline methiodide (BMI), to cat V1 cells only slightly widened the cells' size tuning, while at the same time surround suppression retained its orientation selectivity. Based on this, they argued that size tuning does not arise from intracortical interactions but is more likely to arise from size-tuned cells in the lateral geniculate nucleus. Our data do not support this proposal, but rather argue for cortical mechanisms of spatial integration. There are a few problems with the BMI arguments. Application of BMI does not affect excitatory connections within the cortex, and thus does not address whether spatial summation is mediated intracortically. Moreover, BMI application only affects GABAa receptors, leaving GABAb receptors unaffected, i.e., GABAergic inhibition was only partly affected by BMI injection.

ACh facilitated the response of the majority of cells, and one might therefore argue that the shift of length preference toward shorter bar length might be due to response saturation, i.e., in the presence of ACh neurons fire at maximum level when relatively short bars are presented, and therefore firing rate cannot increase further when longer bars are presented, although synaptic activity might still increase (which can only be determined through intracellular recordings). Contrary to this argument, we also found a shift in length tuning toward shorter bars in neurons suppressed under conditions of high ACh, arguing against an explanation based on response saturation.

Our finding of changes in the spatial integration of cortical neurons on ACh application is contrary to reports from cat dorsal lateral geniculate nucleus (dLGN). Receptive fields of relay cells in that area mainly show increases in gain, with small increases of receptive field size (summation area) (Fjeld et al. 2002). ACh in the dLGN excites relay cells (Eysel et al. 1986; Sillito et al. 1983), it seems to have an inhibitory effect on inhibitory interneurons (McCormick and Pape 1988), and inhibits perigeniculate neurons, which in turn inhibit relay cells (McCormick and Prince 1987). The network involved is thus different from the cortical network (Gil et al. 1997; Kimura 2000; Kimura and Baughman 1997; Kimura et al. 1999), and therefore these differences of ACh action are not necessarily surprising.

Here we show that cortical RF integration can be dynamically modulated by internal factors, such as the neuromodulator ACh. Since the natural release of ACh is bound to states of arousal and attention (Everitt and Robbins 1997; Sarter and Bruno 1997; Sarter et al. 2001), ACh may be involved in dynamic RF changes associated with spatial attention (Connor et al. 1997; Li et al. 2004; Thiele 2004; Treue and Trujillo 1999). It has been suggested that the function of the nCRF is to allow visual neurons to code natural scenes more efficiently by exploiting redundancies in the scene (Young 2000). Coding more efficiently by this mechanism relies on inference about the visual world (Young 2000) and may therefore come at the cost of an increased error rate (Dayan and Yu 2001; Yu and Dayan 2002). By reducing the power of the nCRF, the presence of ACh may cause cells to process stimuli within its CRF more accurately without moderation from the wider context, thereby increasing local information processing, and potentially reducing errors due to inferential processing. One parallel between the effects of attentional modulation and the results here is that attentional modulation largely skips the onset transient (Fries et al. 2001; Roelfsema et al. 1998), which we also found in facilitated cells. However, it might be argued that the cholinergic system lacks the speed and spatial precision to mediate effects of spatial attention which operates rapidly with high spatial resolution. At least a subset of cholinergic nerve terminals establish classical synapses in the cortex (Turrini et al. 2001), which would allow for fast processing compared with volume transmission. Moreover, a variety of experiments show that spatial/regional specificity of ACh release is higher than originally thought (Carey and Rieck 1987; Fournier et al. 2004; Price and Stern 1983). Additionally it might be argued that the speed and spatial specificity of spatial attention is mediated by an interaction of cholinergic input and feedback (synapto/synaptic) connection. This would, however, imply that ACh affects feedback projections in a different manner from its effects on lateral-intracortical connections, something which future experiments will have to determine.

In this study, we have shown that the application of ACh to cells in primate V1 can significantly change the cell's length tuning. This result is compatible with the hypothesis that the effect of cortical ACh is to control the flow of neuronal information such that the efficacy of information arriving from the senses is boosted relative to information from within the cortex. Various earlier lines of evidence have contributed to this hypothesis; however, this is the first study to test it directly in vivo in the primate. A number of features of our data are strikingly similar to data from studies investigating the effects of spatial attention on neuronal processing. For example, our data resemble data from an attention experiment that showed that neuronal response profiles become more tonic during attentive states (Roelfsema et al. 1998). Perhaps more striking is the finding that spatial attention reduces contextual influences (Ito et al. 1998) in a manner similar to our demonstration that the application of ACh reduces contextual influences. To what extent neural effects of attention may be explained by the action of cortical ACh is currently an ongoing project in our laboratory.


This work was supported by the Wellcome Trust (070380/Z/03/Z) and the Medical Research Council UK (G0100407; G78/7853).


We thank K. Dobkins and two anonymous reviewers for valuable comments on the manuscript.

Present address of J. S. McDonald: Max Planck Institute für Kybernetik, Tübingen, Germany.


  • The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


View Abstract