Spatial frequency representation in V2 and V4 of macaque monkey
Abstract
Spatial frequency (SF) is an important attribute in the visual scene and is a defining feature of visual processing channels. However, there remain many unsolved questions about how extrastriate areas in primate visual cortex code this fundamental information. Here, using intrinsic signal optical imaging in visual areas of V2 and V4 of macaque monkeys, we quantify the relationship between SF maps and (1) visual topography and (2) color and orientation maps. We find that in orientation regions, low to high SF is mapped orthogonally to orientation; in color regions, which are reported to contain orthogonal axes of color and lightness, low SFs tend to be represented more frequently than high SFs. This supports a population-based SF fluctuation related to the ‘color/orientation’ organizations. We propose a generalized hypercolumn model across cortical areas, comprised of two orthogonal parameters with additional parameters.
Editor's evaluation
This article makes an important contribution to the study of early visual representation in primates by showing that intermediate cortical areas V2 and V4, as well as primary cortical area V1 (previously shown), contain orthogonal maps of orientation and spatial frequency. The authors provide convincing evidence of this fundamental principle of functional mapping across the two-dimensional cortical surface that ensures and optimizes the complete representation of these two coding dimensions.
https://doi.org/10.7554/eLife.81794.sa0Introduction
Spatial frequency (SF) selectivity is a fundamental feature encoded in the visual system. Previous studies have shown that the organization of SF selectivity is related to orientation and color maps in the primary visual cortex (V1) and have a high degree of periodicity in both cats (Hübener et al., 1997; Issa et al., 2000; Shoham et al., 1997; Tootell et al., 1981) and monkeys (Silverman et al., 1989). Studies have consistently shown an orthogonal mapping of SF and orientation, suggesting an efficient arrangement that provides each orientation access to a wide range of SFs (Issa et al., 2000; Nauhaus et al., 2012; Nauhaus et al., 2016; Xu et al., 2007). In contrast, color representation in V1 (the color is represented in patches commonly referred to as 'blobs' in V1) is generally associated with a range of lower SFs (Silverman et al., 1989; Tootell et al., 1988). Thus, in addition to a gradient from high to low of SF representation across eccentricities (Foster et al., 1985; De Valois et al., 1982; Lu et al., 2018; Yu et al., 2010), SF organization is further specified within local distinct functional regions. This systematic architecture in V1 suggests that SF may be a fundamental feature of the cortical ‘hypercolumn’ (c.f., Silverman et al., 1989: organized cortical modules; Hübener et al., 1997: ‘mosaics’ of functional domains for the different properties; Swindale et al., 2000: uniform coverage of cortical maps).
Whether there are systematic associations between SF and other parameters in extrastriate areas, such as V2 and V4, is not known. The traditional view of the V2 hypercolumn comprises the alternating thin-pale-thick-pale stripe cycle (Horton, 1984; Livingstone and Hubel, 1984; Roe and Ts’o, 1995; Tootell et al., 1983). Within thin stripes, surface properties, typically associated with low SF preferences, such as hue maps (Xiao et al., 2003), ‘brightness’ maps (Roe et al., 2005), and ON/FF maps (Wang et al., 2007) are represented. Within the thick and pale stripes are higher-order orientation maps such as those defined by illusory contours (Ramsden et al., 2001), motion direction maps (Lu et al., 2010), and maps for motion-defined edges (Chen et al., 2016), as well as stereo-defined near-to-far disparity maps (Chen et al., 2008). Neuronal response for features such as texture have also been described (Freeman et al., 2013), but functional organization has not yet been investigated. There is little systematic data relating SF representation in V2 to functional stripes (cf., Gegenfurtner et al., 1996; Levitt et al., 1994; Tootell and Hamilton, 1989) and, despite previous attempts, few studies have demonstrated functional mapping of stripes based on SF alone (Lu and Roe, 2007; Lu et al., 2018).
In V4, surface and shape information are organized into, for lack of better terminology, ‘color’ and ‘orientation’ bands. Within the color bands, maps for hue and for luminance have been described (Tanigawa et al., 2010; Liu et al., 2020; Li et al., 2022). Within orientation bands, there are maps for contrast-defined contours (Hu et al., 2020; Li et al., 2013; Lu et al., 2018; Tang et al., 2020; Tanigawa et al., 2010), disparity-defined contours (Fang et al., 2019), as well as maps for curvature degree and curvature orientation (Hu et al., 2020; Ponce et al., 2017). Despite our growing understanding of functional organization in V4, how SF preference maps (first reported in Lu et al., 2018) relate to other feature maps in V4 remains unknown.
As part of our investigation into ‘hypercolumn’ organization in extrastriate cortical areas, we propose a general hypercolumn layout for V2 and V4 that includes SF (cf., Roe et al., 2009; Ts’o et al., 2009). Based on previous results reported in V1 (Foster et al., 1985; Lu et al., 2018; Nauhaus et al., 2012; Nauhaus et al., 2016; Silverman et al., 1989; Tootell et al., 1988; Yu et al., 2010), we predict that in each area (1) the range of SF associations shift with the topographic location of the ‘hypercolumn’ (Figure 1A), (2) orientation-selective regions (Figure 1B, blue) have a range of low (light gray) to high (dark gray) SF representation; iso-SF contours (purple dashed lines) map orthogonally to iso-orientation contours (green dashed lines), and (3) color-selective regions (Figure 1B, orange) exhibit an association with a range of low SFs (light gray). To address this proposal, we imaged V1, V2, and V4 of macaque monkey via intrinsic signal optical imaging (ISOI) with large cortical fields of view that contained sufficient territory to allow comparisons of functional organization at a range of eccentricities (Figure 2). Quantification of the relationship between SF maps and color and orientation maps in V2 and V4 revealed organizations that generally support our proposal for a hypercolumn architecture.
Results
Overall SF preference across imaged visual cortical areas
Imaging a large field of view of the cortex makes it possible to directly compare the response differences between different visual areas (Figure 2). Figure 3 shows the blood vessel maps (A, E) and corresponding SF preference maps (B, F) for two separate cases wherein we imaged regions spanning V1, V2, and V4. The SF preference maps (Figure 3B and F) were generated by calculating the response amplitude of each pixel in these areas for stimuli with six different SFs, each at orientations of 45° and 135°. We found that in both cases, most of the imaged V1 region favored high SF stimuli, while most of the imaged V4 region favored low SF stimuli. To examine the proportion of cortical area dedicated to a single SF preference in each area (V1, V2, and V4), we calculated for each SF a coverage ratio (the area of SF-preferential response divided by the whole area, i.e., V1, V2, or V4, within the field of view). For V1, the coverage ratio peaks at high SF (light gray bars in Figure 3C and G), while for V4, the coverage ratio peaks at low SF (black bars in Figure 3C and G), consistent with previous findings that V4 prefers lower SF than V1 (Lu et al., 2018). The preferred SFs (Figure 3D and H) significantly decrease from V1 to V4 (see Table 1, Kolmogorov–Smirnov test, V1 vs. V4, V2 vs. V4, V1 vs. V2, p<0.001). Our maps, in contrast to earlier studies (Foster et al., 1985; Levitt et al., 1994; Lu and Roe, 2007; Lu et al., 2018; Silverman et al., 1989), directly show the overall SF preference in cortical space.
To confirm that the use of just two orientations, 45° and 135°, detects a complete picture of SF preference in the imaged area, we compared the SF preference results acquired by different orientations (45° + 135°, 45°, 135°) (see Figure 3—figure supplement 1) and did two-way ANOVA analysis (nine values went into the two-way ANOVA, including coverage ratios of high SF regions or low SF regions in V1, V2, and V4 acquired from different orientation comparisons, 45°, 135°, 45 + 135°). We found that, for different orientations, there were no significant differences (two-way ANOVA, p>0.05) in the coverage ratios of high SF preference (SF = 4 cycles/deg) and low SF preference (SF = 0.25 cycles/deg) regions in the imaged area. In contrast, significant differences were found for different visual areas (two-way ANOVA, p<0.05).
Foveal to parafoveal shift of SF representation in V4
The large-scale imaging allowed us to capture highly structured maps of functional domains (e.g., orientation and color domains in Figure 4A and B, respectively) and reveal changes to those maps for different SF conditions. Below, we present our results on orientation maps and then color maps.
In the imaging results shown in Figure 4D (left panel), we acquired differential orientation maps using six different SFs. We found that the functional map at a specific SF is not always easily ascertained. Gratings with low SFs (e.g., 0.25 cycles/deg, 0.5 cycles/deg) evoked clear selective responses (Figure 4D, thresholded, right panel) in all imaged V4 regions, whereas gratings with high SFs only evoked responses in certain regions, mostly in the lateral part of the cortex (toward foveal representation, red star).
To quantify this, for each tested SF, we analyzed the spatial distribution of the orientation-selective domains (see Figure 4D right panel, only pixels with significantly differential response to 45° vs. 135° were included). For each orientation map obtained at a given SF, we defined the ‘selective activation center’ as the geometric centroid of all significantly activated pixels (Figure 4C, two-tailed t-test, p<0.01). As SF increases, the location of the selective activation center shifts from medial to lateral across the cortex (Figure 4C, blue dots). We confirmed that this remains true for each of the other cases (Figure 4—figure supplement 1, four cases from four hemispheres of three monkeys). Thus, the lack of response in the medial region to high SFs leads to a spatial medial-to-lateral shift of the selective activation center with increasing SF (Figure 4—figure supplement 1B).
We also examined color maps. Many previous V4 imaging studies employed grating stimuli to obtain color maps. However, most of these studies (Li et al., 2014; Tanigawa et al., 2010) have not addressed how SF affects color selective response. To test this, we recorded cortical responses to red/green isoluminance sinusoidal gratings with six different SFs and compared these responses with those to achromatic stimuli of the same SF (Figure 4E). Color domains that showed significantly stronger responses (two-tailed t-test, p<0.01, N = 30) to red/green stimuli were marked in black (Figure 4E, right panel). Similar to orientation domains, under low SF conditions (<1 cycle/deg), color domains were detected in the medial region of the cortex. At high SF, the color-selective response is no longer easily discernible in the medial region (SF = 3, 4 cycles/deg, Figure 4E, right panel). But in the lateral cortical region, color-selective response was detected regardless of SF. Thus, the activation center of color-selective response also shifted from medial to lateral with increasing SF (Figure 4C, orange dots). In a second case, we obtained similar results (Figure 4—figure supplement 2A) and found the activation center shifted laterally as SF increased (Figure 4—figure supplement 2B). We found that, when only a few SFs are tested, the coverage ratio of the orientation and color domains is underestimated (see white numbers in upper corners in Figure 4A, B, D and E), underscoring the importance of testing a wide range of SFs.
In addition to extracting the activation centers for different SFs in V4, we also examined the spatial distributions for each SF separately. We plotted the proportion of the activated pixels at different distances from lateral (distance = 0 mm) to medial under different SF conditions (see Figure 4F and G). As reported in previous V4 studies (Tanigawa et al., 2010; Li et al., 2013), orientation domains and color domains tend to separate in space, forming different functional bands. As shown in Figure 4F and G, at distances of 2, 6, 10, and 14 mm, the percentages of activated orientation regions decrease while the percentage of activated color regions increases. We found that for the lateral orientation and color bands (see the two bands < 6 mm), the percentage values were higher for high SFs (blue/cyan lines vs. red/orange lines). In comparison, for the medial orientation and color bands (see the two bands >6 mm), the percentage values were higher for low SFs (red/orange lines vs. blue/cyan lines).
We found surprisingly that the relative proportions of high to low SF preferences differ for color and orientation. As shown in Figure 5, these proportions differ markedly in foveal vs. parafoveal locations. For each SF, we divided the number of activated pixels in the lateral and the medial parts (see Figure 5—figure supplement 1) by the total number of the activated pixels preferring this SF (pixels selectively activated with a single SF). We found that the proportions of activated pixels in color and orientation domains in the lateral and medial parts of V4 change with SFs in distinct ways. In parafoveal cortex (medial), for both orientation and color domains, the proportion of pixels tends to decrease with increasing SF; in foveal cortex (lateral), for both orientation and color domains, the proportion of pixels increase as SF increases.
The exposed areas of V1 are in the lateral region of the hemisphere, corresponding to the eccentricity of 0–2°; here low SFs barely evoke measurable selective orientation responses (Figure 4—figure supplement 1A, left column, cases 1–3). The orientation map in the lateral region is only apparent at high SFs (Figure 4—figure supplement 1A, >1 cycle/deg). Although the optimal SF differs between V1 and V4, for both areas, the selective activation center moves from medial to lateral as SF increases (Figure 4—figure supplement 1, cases 1–4, dots of different shading correspond to the activation centers at different SFs). In monkey visual cortex, the foveal regions of V1, V2, and V4 are located in the lateral part of the cortex (Gattass et al., 1981; Gattass et al., 1988 see Figure 2D), in agreement with the findings from electrophysiology studies that foveal regions favor high SF (Perry and Cowey, 1985; Schein, 1988; Wässle et al., 1989; Wilder et al., 1996).
To further confirm the spatial relationship of functional maps acquired at different SFs, we calculated the cross-correlation values of these functional maps acquired at different SF conditions. Each response map was divided into two halves: foveal region, the left half close to V4 foveal region; and parafoveal region, the right half away from foveal region. The correlation values of the two halves were calculated separately (Figure 4—figure supplement 3A and B, left matrix: foveal region; right matrix: parafoveal region). For orientation maps (Figure 4—figure supplement 3A), high correlation values (>0.5) appeared in the comparisons among high SF conditions in foveal region. In contrast, high correlation values appeared in the comparisons among low SF conditions in parafoveal region. For color functional maps (Figure 4—figure supplement 3B), higher correlation values were measured under low SF conditions in both foveal and parafoveal regions, while high correlation values appeared only in the foveal region under high SF conditions.
These differences reflect distinct capabilities of different cortical regions (fovea vs. parafovea). Independent of the type of visual information presented (orientation or color), parafoveal regions tend to process the visual input containing relatively low SF components, while visual input containing high SF components is better processed in foveal regions. In V4, the foveal region is capable of processing visual information with a broad range of SF sensitivity.
Relationship between SF and orientation maps in V4 and V2
Having obtained orientation and SF preference maps from the same cortical region, it becomes possible to analyze the spatial relationships between these maps. To identify the regions which are highly selective for orientation, for each pixel, we calculated normalized orientation values (from 0 untuned to 1 highly tuned) and set a threshold of 0.5 or larger (see orientation selectivity maps in Figure 6B and Figure 6—figure supplement 2B; and the selectivity thresholded orientation preference map in Figure 6C). We determined the iso-orientation and iso-SF contours based on the smoothed orientation preference map (see Figure 6A) and SF preference map (see example in Figure 6—figure supplement 1, 18 iso-orientation gradient contours and 5 iso-SF gradient contours). In V2 and V4, the iso-orientation and iso-SF contours that intersected at large angles were found most frequently (see Figure 6G, K and O, Figure 6—figure supplement 2G and K). To demonstrate that the intersection angles are more frequently detected at a large angle, we divided the detected intersection angles into three groups (small: 0–30°; medium: 30–60°; large: 60–90°) and compared the percentage difference among these groups. The results indicate that there are more (percentage value) 60–90° intersection angles than other kinds of intersection angles in the orientation-selective regions (see white patches in Figure 6F, J and N). The percentage of the large angle group is significantly higher than the small (Wilcoxon rank-sum test, p=1.60 × 10–5, n = 15 from five regions, two V2 regions, and three V4 regions) and medium groups (Wilcoxon rank-sum test, p=1.73 × 10–4). In addition, we compare the distribution between groups with strong orientation selectivity (e.g., Figure 6D–O) and weak orientation selectivity (see Figure 6P–U). The percentage difference in large intersection angle (60–90°) is also significant (strong orientation selectivity group, n = 15 from five regions, two V2 regions, and three V4 regions; weak orientation selectivity group, labeled as '60–90° Ref.' in Figure 6V, n = 6 from two V4 regions; Wilcoxon rank-sum test, p=0.0057, see Figure 6V).
SF bias in V4 color domains
Since the first report on V4 SF preference domains (Lu et al., 2018), the relationship between V4 SF domains and V4 color domains remains unclear. It should be noted that SF preference domains (i.e., preference for low vs. high, Figure 7B and C) are distinct from SF preference maps (i.e., six-value SF maps, Figure 3B and F), as the domains are determined by statistical analysis (two-tailed t-test, p<0.01) and effectively distinguish high SF from low SF. Here, by comparing cortical responses recorded using high SF stimuli at two orientations: 45°, 135°, and low SF stimuli (at the same two orientations), we generated differential SF maps (see Figure 7A). A distinct segregation of light and dark regions is visible in this functional map. The black patches are the regions that prefer high SF to low SF, while the white patches have the opposite preference.
We compared the spatial relationship of color domains with SF domains (Figure 7A–D). Figure 7B shows the low SF preference domains color coded in gray relative to the color domains (orange patches). Similarly, high SF preference domains are color coded in gray relative to the same domains in Figure 7C. In this case, low SF preference domains tend to co-localize with V4 color domains to a greater extent than high SF preference domains (21.3% vs 8.3% of all pixels in color domains).
Additional examples are shown in Figure 7—figure supplement 1 (three cases from three different hemispheres). In these four cases, we found that the total amount of cortex able to discriminate high SF vs. low SF (high SF plus low SF pixels) is only a small portion of the total V4 area (coverage ratio of high SF preference domains in V4: 3.9% ± 1.5%, mean ± SD; coverage ratio of low SF preference domains in V4: 6.0% ± 3.7%, mean ± SD.; coverage ratio: the area of SF domains divided by the imaged area of V4). Although the difference is not statistically significant (four cases, Wilcoxon test, p>0.05), in all cases, we found a great tendency for low SF domains to overlap with color domains (the percentage of high SF preference regions in color domains: 3.9% ± 3.3%; the percentage of low SF preference regions in color domains: 13.3% ± 8.4%, Figure 7D).
Change of SF preference in V2
Previous studies in V2 have shown that relying on SF to distinguish different types of stripes (thin, pale, thick) is difficult (Levitt et al., 1994; Lu et al., 2018; Tootell and Hamilton, 1989). However, one study did report a stripe specific SF selectivity (Lu and Roe, 2007). We hypothesize that this controversy is due to the inability of cytochrome oxidase staining to reliably identify stripe type, which functional imaging can securely address. To further explore this, we first examined whether there were measurable SF preference differences in V2. As indicated in the map, the SF preference changed across V2 (same as V2 in Figure 3B, the region between lunate sulcus and V1/V2 border indicated by red arrow). We adopted the same method as Lu (Lu et al., 2018) to obtain a differential SF map (Figure 7F). The low SF preference regions obtained using these two methods (white patches, in the differential SF map, see Figure 7F; red/orange regions, in the colored SF preference map, see Figure 7I and J) are well correlated with color preference domains (red dashed circles in Figure 7F and black patches in Figure 7G and H). In addition, we acquired differential SF maps by subtractions between different SF pairs (see Figure 7—figure supplement 2). For subtraction between high SF (>1 cycle/deg) and low SF (0.5 cycles/deg) conditions, high SF preference domains can be detected in the foveal region (red star, see black patches indicated by red arrows in Figure 7—figure supplement 2A and B). However, for the subtraction between medium SF (1 cycle/deg) and low SF (0.25 cycles/deg) conditions, the relatively high SF preference domains can only be detected in the parafoveal region (see black patches indicated by blue arrows in Figure 7—figure supplement 2C). These results indicate from parafoveal to foveal region, the preferred SF tends to increase gradually.
To provide an additional illustration of how SF preference varies along V2, we selected V2 regions that exhibited two separated color domains (Figure 7E, regions 1 and 2) and compared the changes of SF preference against color selectivity along the V2 long axis (parallel to the black arrows in Figure 7G and H). Both SF preference (red lines) and color selectivity (black lines) vary along the V1/V2 border (Figure 7G–L). Color selectivity was found in or near low SF preference regions (see also Figure 7F, red dashed circles and white patches in V2); however, not all low SF preference regions exhibited strong color selectivity. Thus, SF preference differences vary uniquely within V2 (see Figure 7M). The differences in SF preference between the regions on the two sides of a thin stripe (see Figure 7K and L, thin stripe location: around the peaks of the black lines) may correlate to the reported two types (medial vs. lateral) of pale stripes in V2 (Felleman et al., 2015).
Discussion
Using optical imaging of intrinsic signals in a large area of visual cortex, we were able to simultaneously characterize the functional architecture underlying color, orientation, and SF domains in areas V2 and V4. We found (1) Topography: With respect to topography, the SF population response of V4 orientation and color domains shift systematically with topographic location. To be specific, the geometric centroids of the selective response shift toward foveolar (lateral) parts of the cortex as SF increases. This finding quantifies, at a population level, the gradation of SF preference across the visual field. Interestingly, distinct from V1 fovea, which was selectively responsive to orientations of higher SFs (see Figure 4—figure supplement 1A), the foveal region of V4 was responsive to a broad range of SF from high to low, suggesting greater SF integration at higher cortical levels, even in foveal regions. (2) Orthogonal primary axes: Within V4 and V2, similar to what was previously shown in V1, the gradients of orientation and SF maps tend to be orthogonal to one another. (3) Low SF bias in V4 and V2 color domains: We find there is a tendency for low SF preference domains to overlap with color domains in V4 and V2.
General organization rules of cortical space
As mentioned in many previous studies, orthogonal crossings between different cortical maps facilitate the maximum combination of response properties in a local area with columnar organization. This kind of spatial relationship has now been reported in different animals (cat, monkey) and between different functional maps including orientation vs. ocular dominance, orientation vs. SF in cat area 17 (Hübener et al., 1997); orientation vs. ocular dominance, orientation vs. SF in monkey V1 (Bartfeld and Grinvald, 1992; Nauhaus et al., 2012; Obermayer and Blasdel, 1993); hue vs. lightness in macaque V1 (Li et al., 2022); orientation vs. disparity in macaque V2 (Chen et al., 2008; Ts’o et al., 2009) and, in this study, orientation vs. SF in macaque V2 and V4. These examples suggest that to effectively use cortical space, this orthogonality is established by key visual attributes.
We suggest that this orthogonality may be a common principle and reveals key parameters specific to each cortical area. Thus, for object structure, orientation and SF are two key parameters; for color, the key parameters are hue and luminance (Li et al., 2022). These combinations ensure a complete representation of basic shape and surface information at each cortical level. Thus, cortical mosaics contain distinct regions of orthogonal feature parameters, as observed in the color and orientation stripes of V2 and the color and orientation bands of V4, and may present for other parameters (e.g., face space, object space) in higher cortical areas.
Population-selective responses across the cortex
The ISOI method enables us to study the population-selective responses across a large area of the visual cortex. Consistent with previous findings (Desimone and Schein, 1987; Lu et al., 2018), our results indicate that the majority of V4 regions favor low SF relative to V1 (Figure 3). For both V4 orientation- and color-selective response, there is a tendency for representation to shift to higher SFs from parafovea to fovea (Figure 4 with three supplements, and Figure 5 with one supplement), supporting the inverse relationship between SF preference and retinal eccentricity (Desimone and Schein, 1987; Lu et al., 2018). However, it should be noted that in foveal V4 region, even at low SF, robust orientation- and color-selective responses are detected. This points to an important difference between V1 and V4. That is, foveal representation in V4 may be better organized for processing complex images (e.g., natural scenes) with multiple SF components.
SF preference organization in V2
We explored whether SF is spatially organized relative to the stripes in V2. We found that similar to color-selective response, SF preference changes within the exposed V2 area, forming different SF preference patches, which supports a general functional layout for SF coding in the visual system (preference for low SF in V1 blob, V2/V4 low SF domains; preference for high SF in V1 interblob, V2/V4 high SF domains). Based on the resolution of ISOI, we cannot achieve cellular level resolution. We cannot do further analysis in regions containing neurons with complex selectivity (e.g., orientation pinwheel). These tiny structures are best studied with other methods (e.g., two-photon imaging). More evidence and new techniques (e.g., ultra-high-field 7T fMRI) could also be introduced to test whether the SF preference in the entire V2 changes periodically as other features (e.g., color, orientation, direction, disparity).
Why does SF preference in V2 change in this way? As suggested by modeling V2 retinotopic maps of tree shrews (Sedigh-Sarvestani et al., 2021), in elongated cortical areas such as V2 there tend to be periodic changes in response features across the cortical surface. Another putative functional implication of this periodic distribution is to aid in the integration across features spaces (wiring minimization) via horizontal connections in V2 (Chklovskii et al., 2002; Cowey, 1979; Hubel and Wiesel, 1977; Mitchison, 1991; Koulakov and Chklovskii, 2001).
Thoughts about the hypercolumn
Based on the above findings, we suggest a nested hierarchy of organizations (see Figure 1). At the scale of visual field representation, there is a broad and downward shifting range of SFs from center to periphery. Within this global SF map lie hypercolumns of repeated orientation and color representation, each of which contains two orthogonally arranged primary parameters. In the ‘orientation’ regions, SF is systematically and orthogonally mapped in relation to the orientation map (this study); in the ‘color’ regions, hue and lightness are orthogonally mapped (Li et al., 2022). Analogous to how the primary parameter spaces mapped in each cortical area change from one area to another (e.g., V1: ocular dominance, orientation, color; V2: higher-order orientation, hue, disparity; V4: curvature, 3D shape from shading, hue and lightness, Hu et al., 2020; Srinath et al., 2021; face areas: face maps, Kanwisher et al., 1997, Chang and Tsao, 2017; object areas: object maps, Bao et al., 2020), we hypothesize that subregions of a cortical hypercolumn also organize for different parameters. Thus, while SF is an important primary axis in orientation regions, in color regions which by nature are associated with low SFs, the rationale for a systematic SF map is weakened. One could view the color regions as evolutionary ‘add-ons,’ which became tacked on to the low end of the SF continuum. It should be noted that the functional maps or regions (e.g., orientation, color) in V2 and V4 are not simple repeats of the functional representation in V1, although here they have the same name (orientation, color). These maps or regions are dedicated to processing different visual information in different visual areas. We suggest the architecture of SF representation, which describes distinct SF representations within orientation and color regions, further extends and supports the view of continued parallel streaming of feature-specific pathways.
Materials and methods
Data was acquired from five hemispheres of three adult macaque monkeys (one male and two female, Macaca mulatta). All procedures were performed in accordance with the National Institutes of Health Guidelines and were approved by the Zhejiang University Institutional Animal Care and Use Committee (approval no. ZJU20200022 and ZJU20200023).
Animal preparation
Request a detailed protocolChronic optical windows were implanted in contact with the cortex above areas V1, V2, and V4, containing lunate sulcus and superior temporal sulcus as described previously (Hu et al., 2020, also see Figure 2B). The eccentricity of the visual field corresponding to exposed V4 was 0–5° and for V1/V2 was 0–2°. Following the craniotomy surgery, optical images were acquired in order to generate basic functional maps. Monkeys were artificially ventilated and anesthetized with propofol (induction 5–10 mg/kg, maintenance 5–10 mg/kg/hr, i.v.) and isoflurane (0.2–1.5%). Anesthetic depth was assessed continuously by monitoring heart rate, end-tidal CO2, blood oximetry, and EEG. Rectal temperature was maintained at 38°C. Animals were paralyzed (vecuronium bromide, induction 0.25 mg/kg, maintenance 0.05–0.1 mg/kg/hr, i.v.) and respirated. Pupils were dilated (atropine sulfate 1%) and eyes fitted with contact lenses of appropriate curvature to focus on a stimulus screen 57 cm from the eyes.
Visual stimuli for optical imaging
Request a detailed protocolVisual stimuli were created using ViSaGe (Cambridge Research Systems Ltd.) and displayed on a calibrated 27-inch monitor (Philips 272G5D) operating at 60 Hz refresh rate. The luminance for white stimuli was 206.52 cd/m2 and black was 0.50 cd/m2. Full-screen visual grating stimuli were used to locate different functional domains. To acquire color maps (Figure 4E, Figure 4—figure supplement 2, Figure 7B, C and E, Figure 7—figure supplement 1), red/green isoluminance (red: CIExyY, 0.662, 0.328, 40; green: CIExyY, 0.320, 0.613, 40) and black-white sine-wave drifting grating stimuli, as shown in Figure 2A, were presented at two different orientations (e.g., 45° and 135°) with various SFs. To acquire orientation maps (differential orientation maps: Figure 4D and Figure 4—figure supplement 1; orientation preference maps: Figure 6A, Figure 6—figure supplement 1A, Figure 6—figure supplement 2A; orientation selectivity map: Figure 6B, Figure 6—figure supplement 2B; selectivity thresholded orientation preference map: Figure 6C, Figure 6—figure supplement 2C) and SF maps (differential SF maps: Figure 7A and F, Figure 7—figure supplement 2; SF preference maps: Figure 3B and F, Figure 3—figure supplement 1A–C, G–I, Figure 6—figure supplement 1B, Figure 7I and J), gratings with different orientations (0°, 45°, 90°, 135°) and different SFs (0.25, 0.5, 1, 1.5, 3, 4 cycles/deg) were presented. The temporal frequency of the gratings was fixed to 4 Hz, and the corresponding drifting speeds of these SF conditions (0.25, 0.5, 1, 1.5, 3, 4 cycle/deg) are 16, 8, 4, 2.7, 1.3, 1 deg/s. The different directions of motion were randomly interleaved.
Optical imaging
Request a detailed protocolThe brain was imaged through a glass window mounted in contact with cortex. Images of cortical reflectance changes (intrinsic hemodynamic signals) corresponding to local cortical activity were acquired (Imager 3001, Optical Imaging Inc, German Town, NY) with 630 nm illumination. Image size was 1080 × 1308 pixels representing 14.4 × 17.4 or 8.7 × 10.5 mm field of view. Visual stimuli were presented in a random order. Each stimulus was presented for 3.5–4.5 s. Frames were acquired at 4 Hz for 4–5 s synchronized to respiration. Visual stimuli were presented 0.5 s after beginning image acquisition. The imaging data were stored in a block fashion. Each block contained the imaging data recorded from the stimulus conditions (presented one time). Each stimulus was presented at least 25 times.
Data analysis
Generation of functional maps
Request a detailed protocolWith the following formula,, we assessed the response differences between two comparison groups. and are the mean dR/R values (, R9-end is the averaged response from frames 9 to the last frame, R1-3 is the averaged response from frames 1–3) in the two compared conditions for pixel i, N is the number of trials, and Si is the standard deviation of . Single condition maps were obtained by comparing the images acquired during stimulus and during a blank.
Color maps were obtained by comparing red/green and white/black grating images, differential SF maps were obtained by comparing high (2–6 cycles/deg) and low (0.25–0.5 cycles/deg) SF images, and differential orientation maps were obtained by comparing two orthogonal orientation images (45° vs. 135°). Maps were low-pass filtered (Gaussian filter, ~30–80 μm diameter) and low-frequency noise was reduced by convolving a given map with a~1–2 mm diameter Gaussian filter and subtracting from the original map. Within a single experimental session, the same filtering parameters were always used to ensure that this filtering procedure did not influence the observed differences. The border between V1 and V2 was determined based on color map: in V1 color response has a blob-like distribution, whereas in V2 color response has a stripe-like distribution.
To generate SF preference maps (e.g., in Figure 3B and F), for each pixel we compared its activation under different single SF conditions. The preferred SF of each pixel was defined as the SF at which the strongest activation signal (amplitude averaged from frame 10 to frame 20 in each condition) for that pixel was observed. The comparison includes two orientations (45° and 135°); for each orientation, six different SFs, 0.25, 0.5, 1, 1.5, 3, and 4 cycles/deg, were presented. Each pixel in a given SF preference map was assigned a unique color to represent the preferred SF. Orientation preference maps (Figure 6A, Figure 6—figure supplement 2A) were calculated based on single orientation condition maps (four orientations, 0°, 45°, 90°, 135°), and each pixel was assigned a unique color to represent the preferred orientation (Bosking et al., 1997).
Locating the positions of selective activation and determining the activation center
Request a detailed protocolFunctional domains were identified by selecting the pixels with a significant difference in dR/R (two-tailed t-test, p<0.01) under two comparison conditions (see Table 2).
For a given activated region, the activation center was defined as the geometric centroid of all significantly activated pixels within the region; for example, the orientation-selective activation center in Figure 4D is the centroid of all the pixels of 45° and 135° orientation domains under one SF condition and the color activation center in Figure 4E is the centroid of all the pixels of color domains under one SF condition. The overlap between different functional domains was also calculated based on these thus-defined functional domains (see Figure 7B–D, Figure 7—figure supplement 1).
Calculating the correlation of pairs of maps
Request a detailed protocolTo quantify the correlation between two functional maps, we isolated the significant responses in the imaged area of V4 (regions that were significantly activated by the visual stimuli, two-tailed t-test, p<0.01) and calculated the cross-correlation values between the maps (Figure 4D and E, right panels) acquired under different SF conditions. To compare the difference between foveal and parafoveal regions, we divided the imaged V4 regions into two halves: the left half of the region was designated as foveal, the right half was designated as parafoveal, and the correlation value for each half was calculated separately.
Comparing the spatial relationship between SF and orientation maps
Request a detailed protocolIso-orientation contours (18 contours, 5°, 15°, 25°, 35°, 45°, 55°, 65°, 75°, 85°, 95°, 105°, 115°, 125°, 135°, 145°, 155°, 165°, and 175°) and iso-SF contours (five contours, 0.25, 0.5, 1, 2, and 3 cycles/deg) were drawn based on these smoothed maps using the MATLAB ‘contour’ function. We quantified the orientation selectivity of each pixel by calculating the vector sum of the responses to the four tested orientations (0°, 45°, 90°, and 135°). The length of each vector was normalized to a range of 0–1 by dividing the largest vector length. We calculated the difference between the two gradients at each intersection within the strong orientation-selective regions (normalized orientation selectivity > 0.5, e.g., Figure 6, V2 and V4 parts) or weak orientation-selective regions (normalized orientation selectivity < 0.5, e.g., Figure 6, V4 refs) to determine the spatial relationship between SF and orientation maps.
Characterizing the layout of SF preference in V2
Request a detailed protocolAs reported in previous studies, V2 color-selective response changes periodically along the long axis of V2 (Levitt et al., 1994; Roe and Ts’o, 1995). To characterize the periodic change of SF preference in V2, we chose a region of V2 with clearly identifiable periodic changes in color response (at least two well-separated color domains) for further analysis. We slightly rotated the selected V2 region to align the V1/V2 border horizontally in the cropped small map (see Figure 7G and H). For each of these small maps, we quantified and normalized color selectivity and SF preference for all pixels in the map. The average value for the pixels along each vertical line at different distances from left (distance = 0 mm) to right (distance = 4.5 mm) was then computed and plotted.
Calculating the coverage ratio of selective activation
Request a detailed protocolWe first determined the activated orientation and color domains to calculate the proportion of activated regions in a given area. The coverage ratio (e.g., the values at the corners in Figure 4D and E) was calculated by dividing the number of the selectively activated pixels recorded at each SF by the total number of pixels in the given area. To evaluate the weights of the activated pixels at different distances from lateral (distance = 0 mm) to medial (Figure 4F and G), we calculated by dividing the number of the selectively activated pixels along each vertical line at different distances by the total number of activated pixels with the given SF. To evaluate the response weights of V4 lateral and medial parts in different SFs (Figure 5B and C) at each single SF condition, we calculated the ratio of the activated pixels in lateral and medial parts to all the activated pixels for a given SF.
After merging all the orientation- and color-selective pixels at different SF conditions (at least three SFs, low: 0.25–0.5 cycles/deg; medium: 1–2 cycles/deg; and high: 3–6 cycles/deg), we obtained the regions representing nearly the entirety of orientation/color domains (see example in Figure 4A and B) and calculated the coverage ratio of the functional domains within the area (values in Figure 4A and B).
Data availability
All data generated or analysed during this study are included in the manuscript and supporting file. All data used in the figures have been deposited at Open Science Framework (https://osf.io/agkr7/).
-
Open Science FrameworkID agkr7. Spatial frequency representation in V2 and V4 of macaque monkey.
References
-
Orientation selectivity and the arrangement of horizontal connections in tree shrew striate cortexThe Journal of Neuroscience 17:2112–2127.https://doi.org/10.1523/JNEUROSCI.17-06-02112.1997
-
An orientation map for motion boundaries in macaque V2Cerebral Cortex 26:279–287.https://doi.org/10.1093/cercor/bhu235
-
Cortical maps and visual perception: the grindley memorial lectureThe Quarterly Journal of Experimental Psychology 31:1–17.https://doi.org/10.1080/14640747908400703
-
Visual properties of neurons in area V4 of the macaque: sensitivity to stimulus formJournal of Neurophysiology 57:835–868.https://doi.org/10.1152/jn.1987.57.3.835
-
An orientation map for disparity-defined edges in area V4Cerebral Cortex 29:666–679.https://doi.org/10.1093/cercor/bhx348
-
The representation of orientation in macaque V2: four stripes not threeCerebral Cortex 25:2354–2369.https://doi.org/10.1093/cercor/bhu033
-
A functional and perceptual signature of the second visual area in primatesNature Neuroscience 16:974–981.https://doi.org/10.1038/nn.3402
-
Visual topography of V2 in the macaqueThe Journal of Comparative Neurology 201:519–539.https://doi.org/10.1002/cne.902010405
-
Visuotopic organization and extent of V3 and V4 of the macaqueThe Journal of Neuroscience 8:1831–1845.https://doi.org/10.1523/JNEUROSCI.08-06-01831.1988
-
Processing of color, form, and motion in macaque area V2Visual Neuroscience 13:161–172.https://doi.org/10.1017/s0952523800007203
-
Cytochrome oxidase patches: a new cytoarchitectonic feature of monkey visual cortexPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 304:199–253.https://doi.org/10.1098/rstb.1984.0021
-
Ferrier lecture: functional architecture of macaque monkey visual cortexProceedings of the Royal Society of London. Series B, Biological Sciences 198:1–59.https://doi.org/10.1098/rspb.1977.0085
-
Spatial relationships among three columnar systems in cat area 17The Journal of Neuroscience 17:9270–9284.https://doi.org/10.1523/JNEUROSCI.17-23-09270.1997
-
Spatial frequency maps in cat visual cortexThe Journal of Neuroscience 20:8504–8514.https://doi.org/10.1523/JNEUROSCI.20-22-08504.2000
-
The fusiform face area: a module in human extrastriate cortex specialized for face perceptionThe Journal of Neuroscience 17:4302–4311.https://doi.org/10.1523/JNEUROSCI.17-11-04302.1997
-
Receptive fields and functional architecture of macaque V2Journal of Neurophysiology 71:2517–2542.https://doi.org/10.1152/jn.1994.71.6.2517
-
Perceptual color map in macaque visual area V4The Journal of Neuroscience 34:202–217.https://doi.org/10.1523/JNEUROSCI.4549-12.2014
-
Anatomy and physiology of a color system in the primate visual cortexThe Journal of Neuroscience 4:309–356.https://doi.org/10.1523/JNEUROSCI.04-01-00309.1984
-
Optical imaging of contrast response in macaque monkey V1 and V2Cerebral Cortex 17:2675–2695.https://doi.org/10.1093/cercor/bhl177
-
Neuronal branching patterns and the economy of cortical wiringProceedings. Biological Sciences 245:151–158.https://doi.org/10.1098/rspb.1991.0102
-
Geometry of orientation and ocular dominance columns in monkey striate cortexThe Journal of Neuroscience 13:4114–4129.https://doi.org/10.1523/JNEUROSCI.13-10-04114.1993
-
End-stopping predicts curvature tuning along the ventral streamThe Journal of Neuroscience 37:648–659.https://doi.org/10.1523/JNEUROSCI.2507-16.2016
-
Visual topography in primate V2: multiple representation across functional stripesThe Journal of Neuroscience 15:3689–3715.https://doi.org/10.1523/JNEUROSCI.15-05-03689.1995
-
BookVisual system: functional architecture of area V2In: Squire LR, editors. Encyclopedia of Neuroscience. Elsevier. pp. 331–349.https://doi.org/10.1016/B978-008045046-9.00215-1
-
Anatomy of macaque fovea and spatial densities of neurons in foveal representationThe Journal of Comparative Neurology 269:479–505.https://doi.org/10.1002/cne.902690403
-
Visual cortex maps are optimized for uniform coverageNature Neuroscience 3:822–826.https://doi.org/10.1038/77731
-
Functional organization for color and orientation in macaque V4Nature Neuroscience 13:1542–1548.https://doi.org/10.1038/nn.2676
-
Functional anatomy of macaque striate cortexV. Spatial Frequency. Journal of Neuroscience 8:1610–1624.https://doi.org/10.1523/JNEUROSCI.08-05-01610.1988
-
Functional anatomy of the second visual area (V2) in the macaqueThe Journal of Neuroscience 9:2620–2644.https://doi.org/10.1523/JNEUROSCI.09-08-02620.1989
-
Whither the hypercolumn?The Journal of Physiology 587:2791–2805.https://doi.org/10.1113/jphysiol.2009.171082
-
How do functional maps in primary visual cortex vary with eccentricity?The Journal of Comparative Neurology 501:741–755.https://doi.org/10.1002/cne.21277
-
Spatial and temporal frequency tuning in striate cortex: functional uniformity and specializations related to receptive field eccentricityThe European Journal of Neuroscience 31:1043–1062.https://doi.org/10.1111/j.1460-9568.2010.07118.x
Article and author information
Author details
Funding
China Brain Initiative (2021ZD0200401)
- Anna Wang Roe
National Key Research and Development Program of China (2018YFA0701400)
- Anna Wang Roe
National Natural Science Foundation of China (31627802)
- Anna Wang Roe
National Natural Science Foundation of China (U20A20221)
- Anna Wang Roe
National Natural Science Foundation of China (81961128029)
- Anna Wang Roe
The MOE Frontier Science Center for Brain Science & Brain-Machine Integration
- Anna Wang Roe
National Natural Science Foundation of China (32100802)
- Jia Ming Hu
China Postdoctoral Science Foundation (2020M681829)
- Jia Ming Hu
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Yin Liu and Meizhen Qian for help with the animal experiments. This research was conducted at Zhejiang University and was supported by the China Brain Initiative (grant no. 2021ZD0200401 to AWR), National Key R&D Program of China (grant no. 2018YFA0701400 to AWR), the National Natural Science Foundation of China (grant nos. 31627802, U20A20221, and 81961128029 to AWR; grant no. 32100802 to JMH), China Postdoctoral Science Foundation (grant no. 2020M681829 to JMH), and the MOE Frontier Science Center for Brain Science & Brain-Machine Integration, Zhejiang University.
Ethics
All procedures were performed in accordance with the National Institutes of Health Guidelines and were approved by the Zhejiang University Institutional Animal Care and Use Committee (Permit Number: ZJU20200022 and ZJU20200023).
Copyright
© 2023, Zhang et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 880
- views
-
- 173
- downloads
-
- 10
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Neural diversity can expand the encoding capacity of a circuitry. A striking example of diverse structure and function is presented by the afferent synapses between inner hair cells (IHCs) and spiral ganglion neurons (SGNs) in the cochlea. Presynaptic active zones at the pillar IHC side activate at lower IHC potentials than those of the modiolar side that have more presynaptic Ca2+ channels. The postsynaptic SGNs differ in their spontaneous firing rates, sound thresholds, and operating ranges. While a causal relationship between synaptic heterogeneity and neural response diversity seems likely, experimental evidence linking synaptic and SGN physiology has remained difficult to obtain. Here, we aimed at bridging this gap by ex vivo paired recordings of murine IHCs and postsynaptic SGN boutons with stimuli and conditions aimed to mimic those of in vivo SGN characterization. Synapses with high spontaneous rate of release (SR) were found predominantly on the pillar side of the IHC. These high SR synapses had larger and more temporally compact spontaneous EPSCs, lower voltage thresholds, tighter coupling of Ca2+ channels and vesicular release sites, shorter response latencies, and higher initial release rates. This study indicates that synaptic heterogeneity in IHCs directly contributes to the diversity of spontaneous and sound-evoked firing of SGNs.
-
- Neuroscience
Chronic back pain (CBP) is a global health concern with significant societal and economic burden. While various predictors of back pain chronicity have been proposed, including demographic and psychosocial factors, neuroimaging studies have pointed to brain characteristics as predictors of CBP. However, large-scale, multisite validation of these predictors is currently lacking. In two independent longitudinal studies, we examined white matter diffusion imaging data and pain characteristics in patients with subacute back pain (SBP) over 6- and 12-month periods. Diffusion data from individuals with CBP and healthy controls (HC) were analyzed for comparison. Whole-brain tract-based spatial statistics analyses revealed that a cluster in the right superior longitudinal fasciculus (SLF) tract had larger fractional anisotropy (FA) values in patients who recovered (SBPr) compared to those with persistent pain (SBPp), and predicted changes in pain severity. The SLF FA values accurately classified patients at baseline and follow-up in a third publicly available dataset (Area under the Receiver Operating Curve ~0.70). Notably, patients who recovered had FA values larger than those of HC suggesting a potential role of SLF integrity in resilience to CBP. Structural connectivity-based models also classified SBPp and SBPr patients from the three data sets (validation accuracy 67%). Our results validate the right SLF as a robust predictor of CBP development, with potential for clinical translation. Cognitive and behavioral processes dependent on the right SLF, such as proprioception and visuospatial attention, should be analyzed in subacute stages as they could prove important for back pain chronicity.