Population receptive fields in nonhuman primates from whole-brain fMRI and large-scale neurophysiology in visual cortex
Abstract
Population receptive field (pRF) modeling is a popular fMRI method to map the retinotopic organization of the human brain. While fMRI-based pRF maps are qualitatively similar to invasively recorded single-cell receptive fields in animals, it remains unclear what neuronal signal they represent. We addressed this question in awake nonhuman primates comparing whole-brain fMRI and large-scale neurophysiological recordings in areas V1 and V4 of the visual cortex. We examined the fits of several pRF models based on the fMRI blood-oxygen-level-dependent (BOLD) signal, multi-unit spiking activity (MUA), and local field potential (LFP) power in different frequency bands. We found that pRFs derived from BOLD-fMRI were most similar to MUA-pRFs in V1 and V4, while pRFs based on LFP gamma power also gave a good approximation. fMRI-based pRFs thus reliably reflect neuronal receptive field properties in the primate brain. In addition to our results in V1 and V4, the whole-brain fMRI measurements revealed retinotopic tuning in many other cortical and subcortical areas with a consistent increase in pRF size with increasing eccentricity, as well as a retinotopically specific deactivation of default mode network nodes similar to previous observations in humans.
Editor's evaluation
This study is a detailed, systematic comparison of visually evoked population receptive fields (pRFs) measured non-invasively with MRI and invasively with electrophysiology in the same primate species. The authors show that MRI pRFs provide a good estimate of receptive fields based on multi-unit spiking activity in early visual areas. These results make an important contribution to our understanding of human imaging data in research and in the clinic.
https://doi.org/10.7554/eLife.67304.sa0Introduction
The concept of a receptive field (RF) is crucial for our understanding of the mechanisms underlying perception, cognition, and action. RFs (Hartline, 1938; Sherrington, 1906) typically describe stimulus locations that evoke or modulate neuronal responses, but they can be generalized to different stimulus features such as color or spatial frequency. RFs are usually measured by determining the neuronal firing rate elicited by visual stimuli (Hubel and Wiesel, 1998; Hubel and Wiesel, 1968; Hubel and Wiesel, 1959), but they can also be defined based on other neuronal signals such as subthreshold activity (Priebe, 2008), properties of the local field potential (LFP; Victor et al., 1994), or calcium levels that can, for instance, be measured with fluorescent calcium indicators (van Beest et al., 2021; Bonin et al., 2011).
Noninvasive methods lack the spatial resolution to measure the RF properties of single neurons, but they can characterize the RF properties of the aggregate neural signals being measured. The retinotopic organization of the human brain has now been characterized with functional magnetic resonance imaging for decades (Wandell et al., 2007; Wandell and Winawer, 2011). Early studies used phase-encoding with ‘rotating wedge’ and ‘expanding or contracting ring’ stimuli to identify RF position (Engel, 2012; Engel et al., 1994; Sereno et al., 1995), while later studies increasingly used the ‘population receptive field’ (pRF) method that estimates RF size in addition to position (Dumoulin and Wandell, 2008; Wandell et al., 2007; Wandell and Winawer, 2015; Wandell and Winawer, 2011). The method is popular and has been used to map a range of visual and cognitive functions (Binda et al., 2018; Ekman et al., 2020; Harvey et al., 2020; Harvey et al., 2015; He et al., 2019; Hughes et al., 2019; Mo et al., 2018; Poltoratski et al., 2019; Poltoratski and Tong, 2020; Puckett et al., 2020; Shao et al., 2013; Shen et al., 2020; Silson et al., 2018; Stoll et al., 2020; Thomas et al., 2015; Welbourne et al., 2018; Zuiderbaan et al., 2017), dysfunctions (Ahmadi et al., 2020; Alvarez et al., 2020; de Best et al., 2019; Dumoulin and Knapen, 2018; Green et al., 2019; Schwarzkopf et al., 2014), mechanisms of brain development (Dekker et al., 2019), cortical evolution (Keliris et al., 2019; Kolster et al., 2014; Zhu and Vanduffel, 2019), and information transfer across different brain areas (Haak et al., 2013).
The term ‘pRF’ highlights the analogy to neuronal RFs. It assumes that the blood-oxygen-level-dependent (BOLD) signal measured with fMRI reflects the aggregate response of a large population of neurons within a voxel. Indeed, pRFs of the human visual cortex are qualitatively similar to the RFs of single neurons or multi-unit activity (MUA) in animals (Dumoulin and Wandell, 2008). However, most of the previous comparisons were between species, and between studies that used different techniques to measure pRFs/RFs (Barlow et al., 1966; Hubel and Wiesel, 1968). Exceptions are pRF studies based on intracranial recordings in human patients (Harvey et al., 2013; Winawer et al., 2013), which used a limited number of surface electrodes to measure ECoG or intracranial EEG, but not spiking activity. The pRFs derived from the LFP exhibited similar properties to pRFs derived from BOLD signals, including similar spatial summation characteristics (Winawer et al., 2013). Another study compared BOLD-based pRFs in monkeys to the RF properties of single units in published work (Kolster et al., 2014). Furthermore, Keliris et al., 2019 found that single-unit RFs in one of their monkeys were smaller than pRFs measured with BOLD and proposed another method to estimate RF sizes. Their study included MUA but not the LFP. A systematic within-species comparison of pRFs derived from BOLD, MUA, and LFP has however never been carried out.
Here, we fill that gap with extensive pRF modeling based on BOLD, MUA, and LFP signals in macaque monkeys. The question that neuronal signal forms the basis of the fMRI-BOLD signal has far-reaching consequences for the interpretation of human neuroimaging results and is therefore a topic of ongoing debate and rigorous investigation (Arthurs and Boniface, 2002; Bartels et al., 2008; Boynton, 2011; Drew, 2019; Ekstrom, 2010; Goense and Logothetis, 2008; Logothetis, 2010; Logothetis, 2003; Logothetis et al., 2001; Logothetis and Wandell, 2004; Maier et al., 2008; Schölvinck et al., 2010; Sirotin and Das, 2009; Winawer et al., 2013; Winder et al., 2017). Some studies reported that properties of the BOLD signal resemble features of both neuronal spiking and the LFP (Mukamel et al., 2005; Nir et al., 2007; Rees et al., 2000), others that they resemble the LFP but not spiking (Bartolo et al., 2011; Maier et al., 2008; Niessing et al., 2005; Viswanathan and Freeman, 2007), and yet others that they resemble spiking rather than the LFP (Lima et al., 2014). Here, we examine the degree to which pRFs based on the BOLD signal resemble pRFs based on spiking activity and distinct frequency bands of the LFP (Buzsáki, 2006; Buzsáki and Draguhn, 2004; Einevoll et al., 2013; van Kerkoerle et al., 2014). In nonhuman primates, we measured BOLD-pRFs using whole-brain fMRI and determined neuronal pRFs with large-scale neurophysiological recordings in V1 and V4 (Figure 1). Besides showing the presence of retinotopic information throughout the brain based on the fMRI data, we could directly compare V1 and V4 pRFs obtained with fMRI and electrophysiology. This intraspecies comparison provides new insight into the neurophysiological basis of the BOLD-defined pRFs and offers a benchmark for visual field maps obtained with fMRI.
The original method of estimating pRFs from BOLD responses (Dumoulin and Wandell, 2008) uses a forward modeling approach to fit the location and size of a symmetrical two-dimensional Gaussian to the BOLD responses. This approach is often used to predict neuronal activity elicited by moving bar-shaped stimuli. The RF model minimizes the difference between measured and predicted responses by multiplying the pRF with the stimulus and convolving the result with a hemodynamic response function (HRF), which accounts for the time course of neurovascular coupling (Figure 1). Later refinements implemented a difference-of-Gaussians pRF profile (DoG) to account for center-surround interactions (Zuiderbaan et al., 2012; Figure 1C), with substantial improvements in early visual cortex. Another refinement is the introduction of a static nonlinearity that models nonlinear spatial summation across RFs (Britten and Heuer, 1999; Kay et al., 2013; Oleksiak et al., 2011; Winawer et al., 2013). In such a model, the best parameters indicate subadditive spatial summation in all visual areas (Kay et al., 2013). This means that if stimulus S1 elicits a response R1 and a nonoverlapping stimulus S2 elicits response R2, the response to the combined stimulus, S1 + S2, is smaller than the sum, R1 + R2. For this reason, the nonlinear spatial summation model has also been called the ‘compressive spatial summation’ (CSS) model. A third extension has been the modeling of negative pRFs. Standard approaches tend to only include voxels that show increases in the BOLD signal in response to a stimulus. The inclusion of ‘negative’ pRFs with decreased BOLD activity has revealed the retinotopic organization of a number of areas in the so-called default mode network (DMN) (Szinte and Knapen, 2020). In our analysis of pRFs based on BOLD, MUA, and LFPs, we explored several pRF models, allowing us to investigate the potential presence of nonlinear spatial summation and negative pRFs.
Results
Four macaque monkeys (Macaca mulatta) participated in this study. They were rewarded with fluid for maintaining their gaze inside a 1.5° window around a fixation point that was presented at the center of a frontoparallel screen. While they fixated, a 2° wide bar containing full-contrast moving checkerboards traversed the screen in eight different directions (Figure 1). Two animals performed this task in a 3T horizontal bore MRI scanner. Two other monkeys were each implanted with 1024 electrodes in the visual cortex (V1, V4). They performed the same task while neuronal activity (MUA and LFP) was recorded simultaneously from all electrodes. In the MRI setup, the stimulus covered 16° of the visual field, which was the maximum possible with the monitor located just outside of the scanner bore. The bar traveled across this screen in 20 steps of 2.5 s (1 TR). In the electrophysiology setup, the monitor was closer to the animal, allowing a visual field coverage of 28°. The stimulus bar moved across this aperture in 30 steps of 500 ms. For both the MRI and electrophysiology recordings, we only included data from epochs when the animals maintained fixation for >80% of the time.
After preprocessing (see Materials and methods), we independently fit four pRF models to the average BOLD time courses. These models were (1) a linear pRF model constrained to have positive responses (P-LIN) (Dumoulin and Wandell, 2008), (2) an unconstrained version of the linear pRF model that can also model negative responses (U-LIN), (3) a DoG pRF model (Zuiderbaan et al., 2012), and (4) a nonlinear CSS pRF model (Kay et al., 2013; Figure 1). We used the fitting method to determine pRF size, shape, and location. A cross-validated goodness of fit (R2) was determined by fitting the model to one half of the data and calculating R2 using the other half of the data. Cross-validation allows the comparison of fit quality between models with different numbers of parameters.
pRFs measured with BOLD-fMRI
All models provided good fits to the BOLD time courses in a range of cortical and subcortical areas known to be involved in visual processing. For both monkeys (M1 and M2), we found robust retinotopic information in occipital, temporal, and parietal cortex (Figure 2). pRFs in all these areas were in the contralateral visual field and retinotopic maps were consistent with previous reports (Figure 2B), some of which were more extensive (Arcaro et al., 2011; Arcaro and Livingstone, 2017; Brewer et al., 2002; Janssens et al., 2014; Kolster et al., 2014; Rima et al., 2020; Zhu and Vanduffel, 2019). Weaker and sparser retinotopic information was also observed in the frontal cortex, for example, around the arcuate sulcus (area 8, including the frontal eye fields [FEF]) and in the ventrolateral prefrontal cortex (VLPFC). Throughout this study, we will use a voxel inclusion criterion of R2 > 5% unless otherwise noted. While R2 was generally much higher in visual areas (Figure 2A, bottom panel), this relatively low threshold also reveals retinotopic information in more frontal areas and some subcortical regions (where the signal picked up by surface coils has a lower signal-to-noise ratio [SNR]). Figure 2C shows the number of voxels within a range of areas for which the models explained more than 5% of the variance (Figure 2—figure supplement 1 shows the proportion per region of interest [ROI]). The functional parcellation of visual areas based on field sign inversions around horizontal and vertical meridians lined up well with a probabilistic atlas, co-registered to the individual animal’s anatomy (D99, Reveley et al., 2017).
Subcortically, we could segregate the lateral geniculate nucleus (LGN), pulvinar and some striatal regions from their surrounding areas on the basis of a higher R2. In both monkeys, the LGN of both hemispheres contained clear retinotopic maps (Figure 3A). Retinotopic information was also evident in the bilateral pulvinar of M1, but some of the pRFs in the pulvinar of M2 were large and crossed the vertical meridian, resulting in noisy polar angle maps (Figure 3B). Striatal retinotopy was less pronounced in M1 and even more variable in terms of polar angle maps (Figure 3—figure supplement 1).
We calculated cross-validated R2 values to compare the four pRF models: P-LIN, U-LIN, DoG, and CSS (Figure 1C). A comparison of model performance pooled over subjects and voxels confirmed that there were significant differences across models, with the CSS model outperforming the P-LIN model (Kruskal–Wallis test on all four models, H = 21.33, df = 3, p<0.0001; post-hoc Tukey’s HSD multiple comparisons of mean rank, R2CSS > R2P-LIN, p<0.0001). Indeed, the CSS model fit was better than that of P-LIN in all ROIs in M1 and in 38 out of 39 ROIs in M2 (Wilcoxon signed-rank, p<0.05; no difference in premotor area F7 in M2) (Figure 4). Since the CSS model also provided the best fits for the neurophysiological signals (described later), we report results from this model and only extend the analysis to other models where this is useful (e.g., in the case of negative gain values). The advantage of the CSS model over the P-LIN model generally increased in higher visual areas (Figure 4—figure supplement 1).
The static nonlinearity parameter, or ‘pRF exponent’ that models the nonlinearity of spatial summation (Materials and methods: Equation 4), was in the range of 0.2–0.4 and significantly below 1 in all areas (Wilcoxon signed-rank, one-tailed, all ROIs with more than four voxels R2 > 5%, p<0.001). These values of the pRF exponent represent subadditive (compressive) spatial summation, in accordance with observations in the human visual cortex (Kay et al., 2013; Winawer et al., 2013). The pRF exponent in early visual cortex is comparable to previously reported values for human V1. The exponent value in higher visual areas was similar to that in early visual cortex of the two monkeys, and higher than previously observed in human extrastriate cortex, which suggests that spatial suppression is less pronounced in higher areas of the monkey visual cortex than in higher areas of the human visual cortex.
Both the U-LIN model and the DoG model also performed better than the standard P-LIN model (Kruskal–Wallis, Tukey’s HSD, both ps<0.0001). The DoG had slightly better fits across all pooled voxels than the U-LIN model (Kruskal–Wallis, Tukey’s HSD, p<0.0001). The advantage of the DoG model over the P-LIN model was most pronounced in V1 and decreased in higher visual areas.
Negative pRFs from suppressed BOLD responses
There was a subset of voxels with negative BOLD responses for which both the U-LIN and DoG models provided much better pRF fits than the P-LIN and CSS models (arrows in Figure 4A, Figure 4—figure supplement 1C). We inspected the voxels for which the R2 in the U-LIN/DOG models was at least 5% higher than in the P-LIN model (gray, dashed triangles in Figure 4A). The U-LIN model estimated a negative gain for these voxels (median gain = –0.31, Wilcoxon signed-rank, one-tailed, z = –43.9, p<0.0001) (Figure 4B) and the DoG model returned a high level of inhibition (median normalized suppressive amplitude = 1.14, interquartile range [IQR], 0.98–1.29).
There were two categories of voxels with negative responses (Figure 4—figure supplement 2). For the first category of voxels, the P-LIN model estimated pRFs outside the boundaries of the stimulated visual field. This result suggests that the negative response represents surround suppression that is particularly strong around the fovea (Sereno et al., 1995; Shmuel et al., 2006; Smith et al., 2004). The retinotopy of these voxels is consistent with this explanation. In V1, for instance, the voxels were on the medial side of the occipital pole, which represents the peripheral visual field (Figure 4C). The second category of voxels with negative pRFs appeared to be different. Here, P-LIN and CSS models could not fit any pRF, suggesting purely negative BOLD responses. These voxels were primarily located in the medial occipital parietal cortex (mOP), at the superior border of the superior temporal sulcus in the lateral occipital parietal cortex (lOP) and around the lateral sulcus (LatS), which includes parts of the insula, cingulate, parietal, and premotor cortices (Figure 4D). These areas have all previously been identified as being part of the monkey DMN (Mantini et al., 2011). This finding aligns with recent research in humans that revealed similar negatively tuned pRFs in corresponding nodes of the human DMN (Szinte and Knapen, 2020).
pRF size as a function of eccentricity
As expected, pRF sizes were larger at higher eccentricities. This relationship was evident in all areas with larger numbers of well-fit voxels. pRFs were also larger in higher areas, which exhibited a steeper slope of the eccentricity-size relationship (Figure 5, Figure 5—figure supplements 1 and 2). The differences between the slopes in V1 and V2 are smaller than expected based on previous electrophysiological studies, but this is not uncommon with fMRI (Kay et al., 2013). In one animal, we also unexpectedly observed retinotopy in a number of higher areas, such as the anterior cingulate cortex (Figure 5—figure supplements 1 and 2). This brain region has been studied predominantly in the context of decision-making (Amiez et al., 2006; Fouragnan et al., 2019), but it does have resting state correlations with V1 (Griffis et al., 2017). Our design lacked the power for a more detailed investigation of this retinotopic organization, but this result may inspire future work on brain-wide retinotopic tuning.
Multi-unit spiking activity RFs
We determined the RFs of MUA recorded with chronically implanted electrode arrays (Utah arrays) in areas V1 and V4 in two additional monkeys (M3 and M4) that did not participate in the fMRI experiments. We used a 1024-channel cortical implant, consisting of a titanium pedestal that was connected to 16 Utah arrays (Rousche and Normann, 1998). Each Utah array had 8 × 8 shanks with a length of 1.5 mm. In both monkeys, 14 arrays were placed in V1 and 2 in V4 of the left hemisphere (Figure 6). The stimulus was similar to that used in the fMRI experiments with some small differences due to constraints of the two setups (in the electrophysiology setup, the stimulus covered a larger portion of the visual field because the screen was closer to the animal) and the intrinsic nature of the recorded signals (stimulus steps were faster in the electrophysiology experiments because the electrophysiology signals are much faster than the BOLD response). We fit the four pRF models to the MUA responses and to the LFP power in five distinct frequency bands. We compared the MUA pRFs to a more conventional MUA RF-mapping method (cRFs). For this cRF method, we selected channels with an SNR larger than 3 (i.e., visual responses that were more than three times larger than the standard deviation of the spontaneous activity) and derived the RF borders from the onset and offset of the neuronal activity elicited by a smoothly moving light bar (see Materials and methods). Whenever both methods were able to estimate a pRF and cRF (R2 > 25% for the pRF method, SNR > 3 and R2 > 25% for the cRF method), the estimated locations were highly similar (median distance between pRF and cRF center, V1: 0.34, IQR 0.18–0.49 dva; V4: 0.90, IQR 0.40–1.41 dva). Compared to the P-LIN pRF model, the moving bar method estimated smaller cRFs (median size difference pRFsz-cRFsz: 0.50, IQR: 0.08–0.92) (Figure 7). The CSS model, however, returned pRF size estimates that were very similar to the cRF sizes or even a little bit smaller (median size difference pRFsz-cRFsz: –0.13; IQR: –0.41–0.14), suggesting that nonlinear spatial summation might indeed be better at capturing the RF properties of a small population of neurons than linear summation.
The pRF models that were used in this study are all based on circular (symmetric) RFs. A recent study (Silson et al., 2018) suggested that pRFs in human early visual cortex might be elliptical rather than circular, although this suggestion goes against previous work (Greene et al., 2014; Merkel et al., 2018; Zeidman et al., 2018). A later study demonstrated that the elliptical fits were an artifact of the software that had been used in the analysis (Lerma-Usabiaga et al., 2021). The cRF method separately estimates the width and height of the RF, and can thus be used to calculate a simplified RF aspect ratio to investigate RF symmetry. While this measure differs from RF ellipticity (a 45° tilted ellipse has the same aspect ratio as a circle), it does provide some insight into the symmetry of the MUA RFs. We did observe a few cRFs with aspect ratios (σlarge/σsmall) that were larger than 2 (M3: 18/753; M4: 10/527; together 2.2% of all cRFS), but the vast majority of cRFs in both animals had aspect ratios close to 1 (M3 median: 1.12, IQR: 1.04–1.20; M4 median 1.13, IQR: 1.03–1.22) indicating near-symmetric RFs.
We obtained excellent fits to the MUA for all pRF models (see Figure 1E for an example fit). These pRFs covered a large proportion of the lower-right visual field (Figure 6, Figure 6—figure supplement 1). As expected, the pRFs from electrodes of the same arrays (shown in the same color in Figure 6) were clustered in space, and their locations were in accordance with established retinotopy (e.g., Hubel and Wiesel, 1974). The average R2 (over all electrodes with R2 > 0) was 64% in V1 (M3: 54%; M4: 73%) and 53% in V4 (M3: 37%; M4: 68%), which is substantially higher than the average R2 of 11% in both V1 and V4 for the MRI data (all voxels with R2 > 0; V1 M1: 14%, M2: 8%; V4 M1: 14%, M2: 8%).
Cross-validated comparisons revealed significant differences between the four models (Kruskal–Wallis test on all four models: HV1 = 204, dfV1 = 3, pV1 < 0.0001; HV4 = 13.4, dfV4 = 3, pV4 < 0.01) (Figure 8; similar patterns were present in each individual animal). Post-hoc pairwise comparisons (Tukey’s HSD) revealed that the CSS and DoG models provided a better fit than the linear models, although for V4 the advantage of the CSS model over the P-LIN model was only significant when electrodes with a poor fit (R2 < 25%) were excluded (V1, all electrodes: CSS vs. P-LIN, p<0.001, DoG vs. P-LIN, p<0.001; V4, all electrodes: CSS vs. P-LIN, p=0.16, DoG vs. P-LIN, p<0.001; V1, electrodes with R2 > 25%: CSS vs. P-LIN, p<0.001, DoG vs. P-LIN, p<0.001; V4, electrodes with R2 > 25%: CSS vs. P-LIN, p<0.02, DoG vs. P-LIN, p<0.01) (Figure 8). The improved fit of the DoG model was caused by the suppressive surround (median normalized suppressive amplitude = 0.71, IQR 0.56–0.86).
The pRF exponent of the CSS model was 0.38 ± 0.23 in V1 (mean ± SD; 1358 electrodes with R2 > 25%), which was significantly smaller than 1 (Wilcoxon signed-rank, one-tailed, z = –31.59, p<0.0001) and similar to the MRI-based values, which had a mean of 0.34 ± 0.19 (6341 voxels). Likewise, in V4, MUA pRF exponent values were significantly smaller than 1 (0.33 ± 0.19; z = 11.13, p<0.0001; n = 165), and comparable to MRI-based values in the same area (0.30 ± 0.16; n = 2324). The similarity in the values of the pRF exponent indicates that subadditive spatial summation is a prominent feature in these areas.
The size of the estimated MUA pRFs increased with eccentricity (Figure 9). However, the pRF sizes for two of the V1 arrays in both monkeys were approximately three times smaller than expected when compared to the data from the other arrays. These outlying arrays were located on the posterior-medial side of the surface of V1 in both monkeys where the gray matter is relatively thin (Figure 9A). We therefore suspect that the 1.5-mm-long shanks of the Utah arrays were pushed into the white matter, where they picked up activity of thalamic afferents (i.e., the geniculostriate pathway). We therefore excluded these electrodes from the analyses of pRF sizes. The remaining MUA pRF sizes and eccentricities were highly similar to RFs reported in previous electrophysiology studies at similar eccentricities (Gattass et al., 1987; Gattass et al., 1981; Van Essen et al., 1984; Victor et al., 1994). Next, we compared pRFs between the electrophysiological signals and the fMRI-BOLD signal.
Local field potential pRFs
The LFP was split into five frequency bands: θ (4–8 Hz), α (8–16 Hz), β (16–30 Hz), γlow (30–60 Hz), and γhigh (60–120 Hz). We fit the four models to estimate pRFs for each frequency band. The results for γlow and γhigh resembled those for the MUA with high R2 values for a large proportion of the electrodes (especially in V1). The CSS model again outperformed the other models, and we did not observe negative pRFs in these V1 and V4 regions (Figure 8—figure supplements 1 and 2). Electrodes with good MUA-pRF fits usually also had good LFP-pRFs, but the opposite was not always true (Figure 8—figure supplement 3). The pRFs in lower frequency bands differed. Whereas θ generally yielded low R2 values, α and β yielded good fits for a substantial number of electrodes. Interestingly, there were two classes of electrodes in these frequency bands. For the first class, the power increased with visual stimulation and CSS model fits were best. In contrast, the second class of electrodes had negative pRFs, that is, the stimulus suppressed power (Figure 8—figure supplement 1). Given the few electrodes with good low-frequency LFP-pRFs in V4, we focused our analysis on the positive and negative responses on V1 (the split into positive/negative pRF was based on the parameters of the U-LIN model; Figure 10, Figure 10—figure supplement 1).
We observed a number of remarkable differences between the positive and negative pRFs. First, negative pRFs generally had lower eccentricities than positive pRFs (Figure 10B, Figure 10—figure supplement 1B). Second, negative pRFs were larger than positive pRFs (Figure 10C, Figure 10—figure supplement 1C). Third, there was a systematic difference in the distance from the γ-pRF of the same electrode. Specifically, we calculated a ‘separation index’ by dividing the distance between the centers of the low- and high-frequency LFP-pRFs by their summed size estimates (SI = Distance/(σlf + σhf)). A separation index of less than 1 indicates pRF overlap. For both α and β power, the negative pRFs were farther from the γ-pRFs than the positive pRFs (Figure 10D, Figure 10—figure supplement 1D) and generally shifted in the direction of fixation (Figure 10E, Figure 10—figure supplements 1E and 2). These results suggest that the positive pRFs represent visually driven activity, whereas the negative pRFs signal a form of suppression that is strongest at smaller eccentricities, close to the fixation point. One possible explanation is that the monkeys directed attention to the fixation point, which may have caused a ring of suppression closely surrounding it. Another possibility is that the negative pRFs might be a consequence of small eye movements around the fixation point. The size of α-pRFs did not depend on eccentricity (Figure 10F), whereas the size of β-pRFs increased with eccentricity, although this relation was very weak for positive pRFs (Figure 10—figure supplement 1F).
The pRFs derived from all electrophysiology signals on the same electrode had similar locations, with less than 1 dva between their centers on average (CSS model, Figure 11A). We next analyzed RF sizes, normalizing the size estimates to the MUA-cRF. All LFP-pRFs estimates were larger than MUA-pRFs, and lower frequency LFP components yielded larger pRFs than the higher frequencies (Figure 11B; we observed the same pattern present in each animal). The pRF exponent was well below 1 for all LFP components (V1 and V4 electrodes with R2 > 25%; Wilcoxon signed-rank, one-tailed <1, p<0.001) and smaller at lower frequencies, indicative of stronger CSS. We also compared the exponents to those of the MRI-pRFs. The exponent of the MRI was not significantly different from that of γlow in V1, and both γlow and γhigh in V4 (p>0.05, Figure 11C). Differences between the MRI exponent and those of the other LFP bands were significant (V1: Kruskal–Wallis, H = 505.49, df = 6, p<0.0001; V4: Kruskal–Wallis, H = 21.65, df = 3, p<0.001; Tukey’s HSD for multiple comparisons).
Comparison of pRF eccentricity-size relationship between fMRI and electrophysiology signals
We next compared the eccentricity-size relationship of the electrophysiological signals to that of the BOLD-fMRI pRFs using linear mixed models (LMMs) to evaluate the pRF estimates of the CSS model, separately for V1 and V4 (electrodes with R2 > 50% and R2 > 5% for MRI). We included voxels for which the pRF was in the lower-right visual quadrant where we had electrode coverage. In V1, positive eccentricity-size relationships existed for MRI, MUA, β, γlow, and γhigh (Figure 12A). We compared these signals with a single LMM, revealing an interaction between signal type and eccentricity (F = 14.44, df = 4, p<0.0001), which indicated that the slopes differed. We further used pairwise LMMs to compare the slope of the MRI-pRFs to the electrophysiological signals with a significant eccentricity-size relationship. The fMRI eccentricity-size relationship was similar to that of MUA (F = 0.38, df = 1, p=0.54), whereas it was significantly different from all LFP signals (β: F = 7.97, df = 1, p<0.01; γlow, F = 20.02, df = 1, p<0.001; γhigh, F = 20.05, df = 1, p<0.001) (Figure 12C). We repeated this analysis in V4 (Figure 12B) where we had fewer electrodes. We obtained positive eccentricity-size slopes for MUA, γlow, and γhigh but did not obtain fits of sufficient quality for the θ, α, and β bands. In contrast to V1, there was no clear difference across these signals (F = 1.11, df = 3, p=0.24). Pairwise comparisons with the MRI size-eccentricity relationship did not reveal a difference between those for γlow (F = 0.36, df = 1, p=0.55), γhigh (F = 0.28, df = 1, p=0.60), or MUA (F = 2.40, df = 1, p=0.12). Because the visual field coverage in V4 differed between the two animals that were used for electrophysiology, we repeated the analysis in both individuals and observed similar results (not shown).
To further investigate the robustness of these results, we repeated the analysis (1) with the inclusion of either all V1 and V4 voxels or a subset of the voxels in approximately the same anatomical location as the electrode arrays and (2) by varying the R2-based inclusion criteria for electrodes and voxels (Figure 12—figure supplement 1). MUA-pRFs in V1 and V4 were generally similar to the BOLD-pRFs, although in V4 γlow and γhigh also approximated the fMRI results in some of the comparisons.
Discussion
The current study conducted a systematic comparison of pRFs obtained with fMRI and electrophysiological recordings in V1 and V4 of awake behaving macaque monkeys. Within the same species, we fit several pRF models to seven different signal types (BOLD, MUA, and the power in five LFP frequency bands) to gain insight into the neuronal basis of MRI-BOLD-pRF measurements (Wandell and Winawer, 2015). Our results demonstrate retinotopic tuning in many brain regions and the presence of negative pRFs in areas of the DMN. We found that subadditive spatial summation is a prominent feature of many measures of brain activity. Furthermore, the results establish a clear relationship between BOLD-pRFs and electrophysiologically determined pRFs, with MUA-pRFs being most similar to BOLD-pRFs (as discussed below).
Cortical and subcortical retinotopic tuning of the BOLD signal
Retinotopic information was present in occipital, temporal, parietal, and frontal cortex as well as in a subcortical areas. We could differentiate the LGN and pulvinar from their surrounding areas based on their higher R2. Retinotopic maps in the LGN were roughly in line with previously published retinotopic organization of the macaque LGN (Erwin et al., 1999), but they comprised few voxels and the variability across animals prohibited detailed inferences. In the pulvinar, the organization of the retinotopic information was even more variable across animals. In humans, subcortical retinotopic maps have been observed in LGN, pulvinar, superior colliculus (SC), thalamic reticular nucleus (TRN), and substantia nigra (SN) (Cotton and Smith, 2007; DeSimone et al., 2015; Schneider et al., 2004). Here, we also observed voxels with pRFs located lateral and medial to the pulvinar, but more targeted investigations are necessary for detailed individual-subject segmentation of the thalamus (DeSimone et al., 2015; Tani et al., 2011).
We obtained good pRF fits in the occipital, temporal, and parietal cortical areas, regions for which previous studies demonstrated retinotopic organization with phase-encoded retinotopic mapping (Arcaro et al., 2011; Janssens et al., 2014; Kolster et al., 2014; Kolster et al., 2010; Kolster et al., 2009; Patel et al., 2010). Despite the fact that checkerboard stimuli may not be ideal for the activation of frontal areas, which are better driven by more complex stimuli (Janssens et al., 2014; Saygin and Sereno, 2008), we also observed retinotopy in several frontal areas, including the insula, cingulate cortex, FEFs (area 8), orbitofrontal cortex, ventromedial prefrontal cortex, and dorsolateral prefrontal cortex. These areas are thought to be involved in visual processing and visual attention. Retinotopic maps have previously been reported in the cerebellum in humans (van Es et al., 2019). We did not detect such maps, likely because of a lower SNR in the cerebellum due to lower field strength, a less optimal coil placement, and the sphinx position of the monkeys.
Negative pRFs
We observed negative pRFs in the MRI data and in the low-frequency LFP, but they were of a different nature. We identified two classes of negative pRFs in the MRI data. Negative pRFs in the visual cortex were usually accompanied by positive responses at peripheral visual field locations, close the boundaries of the visual display. Similar negative visual BOLD responses have been reported in human visual cortex (Smith et al., 2004) and are presumably caused by surround suppression (Allman et al., 1985; Cavanaugh et al., 2002; Hubel and Wiesel, 1962; Knierim and van Essen, 1992). We found a second class of negative pRFs around the LatS, in the mOP, and at the superior border of the superior temporal sulcus in the IOP, all areas previously implicated in the monkey’s DMN (Mantini et al., 2011). A similar retinotopy of negative responses was recently found in the human DMN (Szinte and Knapen, 2020). The retinotopy of nodes of the DMN in humans and monkeys implies a sensory-based organization of this part of the DMN, which could play a role in visual cognition (Arsenault et al., 2018).
Negative pRFs were also observed for a subset of the electrodes in the α and β range of the LFP, although pRF fits were generally of lower quality than those for the γ-LFP or MUA signals. It is conceivable that these negative pRFS reflect shifts of attention, which are known to modulate α and β power (Griffiths et al., 2019; Siegel et al., 2008; Womelsdorf and Fries, 2007; Worden et al., 2000). In accordance with this idea, negative α/β pRFs were larger than their positive counterparts and shifted towards the fixation point. Future studies are needed for a closer examination of this phenomenon.
The relation between BOLD-based pRFs and electrophysiologically determined pRFs
We compared BOLD-based pRFs to electrophysiological measures in terms of their spatial summation characteristics, size, and the relationship between eccentricity and size (Amano et al., 2009; Dumoulin and Wandell, 2008; Felleman and Van Essen, 1987; Gattass et al., 2005; Kay et al., 2013; Larsson and Heeger, 2006; Van Essen et al., 1984; Victor et al., 1994). Eye movements can have an effect on the pRF size estimates, and it is therefore important to note that (1) we only included data from recordings where the animals maintained a high fixation performance (IQR-span of the horizontal and vertical eye position, M1: 0.23, 0.34 dva; M2: 0.36, 0.49 dva; M3: 0.18, 0.38 dva; M4: 0.37, 0.60 dva), and (2) we averaged across multiple stimulus presentations to obtain robust response profiles. Hence, variations in eye position can have had only minor effects on the present results.
Compressive spatial summation
The CSS pRF model includes nonlinear spatial summation (Kay et al., 2013). Previous studies observed subadditive spatial summation, also known as spatial compression, throughout the human visual cortex in BOLD-based pRFs, with stronger compression in higher visual areas than in V1 (Kay et al., 2013; Winawer et al., 2013). In the present study, spatial compression was present in subcortical and cortical areas with good BOLD-pRF fits, indicating that it is a widespread, if not universal, characteristic of the primate visual system. The strength of compression did not differ much across areas, and it was similar to that in human V1. For the MUA data, the CSS model estimated pRF sizes that were very similar to RF size estimates derived from conventional methods that use moving luminance bars, whereas the P-LIN model systematically returned larger pRF estimates. This difference might indicate that spatial compression indeed better captures the neuronal RF properties, at least in V1 and V4.
In human iEEG recordings, spatial compression in the broadband iEEG signal was reported to be similar in strength to that of the BOLD signal (Winawer et al., 2013). We also observed CSS in both the MUA-pRFs and in all frequency bands of the LFP, and the pRF exponent tended to be smaller for lower-frequency components, indicative of stronger compression. The lower frequencies also exhibited stronger suppression than the BOLD signal. In contrast, the pRF exponent of MUA and γlow was similar to that of BOLD-fMRI in both V1 and V4, but in V1 the exponent for γhigh was larger. Hence, MUA and γlow were the two electrophysiological markers for which spatial compression resembled that of BOLD-fMRI most.
pRF size and the eccentricity-size relationship
As expected, pRFs were larger at higher eccentricities, with a larger slope in higher visual cortical areas. The pRF sizes and the eccentricity-size relationships that we found in visual cortex were in line with previous results from human and monkey neuroimaging studies (Amano et al., 2009; DeSimone et al., 2015; Dumoulin and Wandell, 2008; Kay et al., 2013; Keliris et al., 2019; Kolster et al., 2014; Kolster et al., 2010; Welbourne et al., 2018; Zhu and Vanduffel, 2019; Zuiderbaan et al., 2012). Furthermore, the results are in keeping with previous electrophysiological recordings of single units and MUA in early visual areas such as V1 (Gattass et al., 1987; Gattass et al., 1981; Van Essen et al., 1984; Victor et al., 1994), V2 (Burkhalter and Van Essen, 1986; Gattass et al., 1981; Rosa et al., 1988), V3 (Burkhalter and Van Essen, 1986; Felleman and Van Essen, 1987; Newsome et al., 1986; Rosa et al., 2000), and V4 (Gattass et al., 1988). Finally, they also match the results of human electrophysiology with subdural electrodes (Harvey et al., 2013; Yoshor et al., 2007).
The spatial scope of fMRI-BOLD, MUA, and LFP
The spatial scale of recorded neural signals depends both on the nature of the signal and the recording method. For extracellular electrophysiology recordings, the impedance and size of an electrode determine its sensitivity for single-unit spiking activity or MUA. This relationship is complex (Viswam et al., 2019; Ward et al., 2009), but high-impedance (~1 MΩ) electrodes are generally better suited for the detection of spiking activity, with smaller contact sites and higher impedance thought to sample from smaller neural populations and being more likely to pick up single-neuron activity. The spatial range of the MUA signal is ~140 µm (Buzsáki and Draguhn, 2004). Electrode impedance, size, and shape are less important for LFP recordings, at least within commonly used ranges (Nelson and Pouget, 2010; Viswam et al., 2019). The spatial extent of the LFP signal is a topic of ongoing debate (Kajikawa and Schroeder, 2011) with some authors estimating it as low as 120–250 µm in visual cortex (Katzner et al., 2009; Xing et al., 2009), while others suggest it may stretch up to several millimeters (Berens et al., 2008; Katzner et al., 2009; Kreiman et al., 2006), and may even be detected more than a centimeter away from the source (Kajikawa and Schroeder, 2011). The hypothesis that low-frequency components of the LFP have a larger spatial reach than high-frequency components is supported by computational modeling (Leski et al., 2013), although it has also been challenged (Dubey and Ray, 2016). Yet, the exact origin of the LFP and its relationship to spiking remain to be completely understood (Einevoll et al., 2013).
Given the larger spatial spread of the LFP compared to MUA, one would expect pRFs based on LFPs to be larger than those based on MUA. This is indeed what we found. The size of LFP-pRFs furthermore depended on the frequency component of the LFP, with lower frequencies yielding larger pRFs, suggesting a frequency-dependent spread of visual information in the LFP signal. In fMRI, there is a trade-off between spatial and temporal resolution. While higher magnetic field strengths and specialized acquisition methods are continuously increasing spatiotemporal resolution, studies with voxel sizes on the order of 1–2 mm isotropic and a repetition time of 2–3 s at 3T field strengths are still common. Due to the complex relationship between the hemodynamic BOLD signal and the underlying neural activity, it is difficult to predict its stimulus sensitivity from neuroimaging parameters such as voxel size. Nonetheless, a comparison of the eccentricity-size relationship between MRI and electrophysiology in V1 revealed that the BOLD-based pRFs resembled MUA in V1 and V4. In V4, however, the BOLD slope was also very similar to that of the LFP gamma power (Figure 12—figure supplement 1). One possible reason for this difference between areas is that V4 is smaller and more heterogeneous (Kolster et al., 2014; Zhu and Vanduffel, 2019) than V1. A V4 voxel; therefore, samples from a neuronal population with more heterogeneous spatial tuning than an equally sized V1 voxel may therefore reflect the activity of a larger population of neurons, which is better approximated by gamma power. Other factors that might play a role are potential differences in the quality of the recorded signal across areas or partial volume effects that are likely to be more prevalent in V4 than in V1.
Conclusions
Our comparison of fMRI with large-scale neurophysiological recordings in visual cortex revealed that pRFs derived from the BOLD signal resemble MUA RFs. Subadditive spatial summation is a general feature of many brain areas and occurs for BOLD, MUA, and LFP. We observed negative pRFs in the monkey DMN and as part of center-surround organization of pRFs in early visual areas. The spatial compression and the eccentricity-size relationship of BOLD resembles that of MUA, but also bears a resemblance to gamma power. We conclude that BOLD-pRFs accurately represent the spatial tuning of the underlying neuronal populations.
Materials and methods
Subject details
Request a detailed protocolFour male macaques (M. mulatta; 7–12 kg, 5–8 years old) participated in this study. Animal care and experimental procedures were in accordance with the ILAR’s Guide for the Care and Use of Laboratory Animals, the European legislation (Directive 2010/63/EU), and approved by the Institutional Animal Care and Use Committee of the Royal Netherlands Academy of Arts and Sciences and the Central Authority for Scientific Procedures on Animals (CCD) in the Netherlands (license numbers AVD8010020173789 and AVD8010020171046). The animals were socially housed in an enriched specialized primate facility with natural daylight, controlled temperature and humidity, and fed with standard primate chow, supplemented with raisins, fresh fruits, and vegetables. Their access to fluid was controlled, according to a carefully designed regime for fluid uptake. During weekdays, the animals received diluted fruit juice in the experimental setup. We ensured that the animals drank sufficient fluid in the setup and received extra fluid after experimental sessions if needed. On the weekends, animals received at least 700 ml of water in the home cage. The animals were regularly checked by veterinary staff and animal caretakers, and their weight and general appearance were recorded in an electronic logbook on a daily basis during fluid-control periods.
Surgical procedures
Request a detailed protocolTwo animals (M1 and M2) participated in the MRI experiments and were implanted with an MRI-compatible plastic (PEEK) head-post, fixed to the skull with ceramic bone screws and acrylic (Papageorgiou et al., 2014; Vanduffel et al., 2001). Anesthetics, analgesics, and monitoring procedures were similar to previous surgical procedures in our laboratory and are described in detail elsewhere (Klink et al., 2017; Poort et al., 2012; Supèr and Roelfsema, 2005). Two other animals (M3 and M4) participated in the electrophysiology experiments. They were implanted with a custom 3D-printed titanium head-post that was designed in-house based on a CT scan of the skull, aligned to a T1-weighted anatomical MRI scan of the brain (Chen et al., 2020; Chen et al., 2017). The titanium head-post was attached to the skull with titanium bone screws, and the skin was closed around the implant without the use of any acrylic. In a second surgery, each animal was additionally implanted with a total of 1024 electrodes in 16 Utah electrode arrays (Blackrock Microsystems) in their visual cortices (14 arrays in V1, 2 arrays in V4; Figure 1B). Each array contained an 8-by-8 grid of 64 iridium oxide electrodes with a length of 1.5 mm spaced at a distance of 400 μm from each other. Pre-implantation electrode impedances ranged from 6 to 12 kΩ. A custom-designed 1024-channel pedestal was attached to the skull with titanium bone screws, and the skin was closed around it. More details on the surgical procedures have been published elsewhere (Chen et al., 2020; Chen et al., 2017).
Visual stimuli and procedures
Request a detailed protocolIn the MRI experiment, animals were head-fixed, sitting in the sphinx position (Papageorgiou et al., 2014; Vanduffel et al., 2001), and viewing a 32″ screen (1920 × 1080 pixels, 100 Hz) (Cambridge Research Systems) at the end of the bore, 130 cm away. pRFs were measured using conventional moving bar stimuli that traversed the screen in eight different directions behind a large virtual circular aperture (Figure 1). The borders of this virtual aperture were invisible because both the foreground and background had the same gray level (22.3 cd/m2) (Dumoulin and Wandell, 2008). In the MRI experiments, the bar sweep spanned 16° (diameter) in 20 steps (Dumoulin and Wandell, 2008). The moving bars were 2° wide and contained a checkerboard pattern (100% contrast; 0.5° checkers; luminance of white checkers: 106.8 cd/m2; luminance of black checkers: 0.2 cd/m2) that moved parallel to the bar’s orientation. Each bar position was on the screen for 2.5 s (1 TR), making one full bar sweep last 50 s. Bar sweep series (all directions presented once) were preceded and followed by 37.5 s (15 TRs) of neutral gray background. Each horizontal or vertical bar sweep was followed by a neutral gray background period of 25 s. The order of the bar sweep directions was 270°, 315°, 180°, 225°, 90°, 135°, 0°, 45° on most runs, but for one animal we inverted the directions to 90°, 135°, 0°, 45°, 270°, 315°, 180°, 225° on some runs to compensate for the animal’s tendency to fall asleep near the end of runs. We included data from 8 scanning sessions for monkey M1 (34 runs, 268 bar sweeps) and 10 sessions for monkey M2 (66 runs, 406 bar sweeps). During stimulus presentation, the animals received fluid rewards (Crist Instruments, Hagerstown, MD) for maintaining fixation within a circular fixation window with a diameter of 2°, centered on a 0.15° red fixation dot, surrounded by a 0.75° (diameter) aperture of neutral gray background color. In the electrophysiology experiments, the stimulus and task were very similar, but bar sweeps now spanned 28°, which was possible because the animals were closer to the monitor (luminance values of this monitor were black: 0 cd/m2; white: 92.1 cd/m2, neutral gray: 14.8 cd/m2). Bars traveled along this path in 30 steps of 500 ms, and the neutral gray luminance intervals were reduced to 2.5 s due to the much faster neuronal responses (compared to the BOLD signal). The fixation window in the electrophysiology experiments was slightly smaller with a diameter of 1.5°. In the MRI experiment, eye position and pupil diameter were tracked with an MRI-compatible infrared eye-tracking system at 120 Hz (ISCAN ETL-200). Hand positions were also monitored using fiber optic amplifiers (Omron E3X-NH) and optic fibers. To reduce body movement-related imaging artifacts, the animals were trained to maintain their hands inside a response box by making reward delivery contingent on both eye and hand position. In the electrophysiology experiments, animals were head-fixed in a conventional vertical primate chair and viewed a 21″ CRT monitor (1024 × 768, 85 Hz) at a distance of 64 cm while their eye position and pupil diameters were tracked at 230 Hz using an infrared eye tracker (TREC ET-49B, Thomas Recording GmbH).
MRI acquisition
Request a detailed protocolMRI was performed in a standard Philips Ingenia 3.0T horizontal bore full-body scanner (Spinoza Center for Neuroimaging, Amsterdam, the Netherlands). We used a custom-built eight-channel phased array receive coil system (Ekstrom et al., 2008; Kolster et al., 2009) (KU Leuven) and the scanner’s full-body transmit coil. Functional images were obtained using a gradient-echo T2* echo-planar sequence (44 horizontal slices, in-plane 72 × 68 matrix, TR = 2500 ms, TE = 20 ms, flip angle = 77.2°, 1.25 × 1.25 × 1.25 mm isotropic voxels, SENSE-factor of 2 in the AP direction, and phase-encoding in the AP direction).
fMRI preprocessing
Request a detailed protocolAll fMRI data were preprocessed with a custom-written Nipype pipeline that we have made available online (RRID:021813; https://github.com/visionandcognition/NHP-BIDS). In short, MRI scans were exported from the scanner as DICOM images and converted to NIFTI files with the dcm2niix tool (Li et al., 2016). The volumes were then reoriented to correct for the animal being in the sphinx position and resampled to 1 mm3 isotropic voxels. The resulting images were realigned using a nonrigid slice-by-slice registration algorithm based on AFNI tools (Cox, 1996) followed by an FSL-based motion correction procedure MCFLIRT (Jenkinson et al., 2002). Functional volumes were linearly aligned to the individual high-resolution anatomical volumes, which were in turn nonlinearly registered to the NMT standard space (Seidlitz et al., 2018). Preprocessed data were further processed with a combination of custom-written MATLAB (MathWorks, Natick, MA) and shell scripts (https://gin.g-node.org/ChrisKlink/NHP-PRF). BOLD time courses for each voxel were normalized to percentage signal change and averaged across runs (or parts of runs) for which fixation was maintained at 80% of the time or more. We averaged odd and even runs separately to allow for a cross-validation approach in the evaluation of the pRF model fits. Anatomical ROIs were defined based on a probabilistic atlas (Reveley et al., 2017; Seidlitz et al., 2018) and refined using individual retinotopic maps.
Post-fit comparisons across pRF models, HRFs, and ROIs were performed in MATLAB based on the volumetric results. For visualization of the fMRI data, volumetric results were also projected to the individual cortical surfaces. To create these surfaces, we averaged multiple anatomical scans (T1-weighted, 3D-FFE, TE = 6 ms, TR = 13 ms, TI = 900 ms, flip angle = 8°, 100 horizontal slices, in-plane 224 × 224 matrix, 0.6 × 0.6 × 0.6 mm isotropic voxels, and phase-encoding in the AP direction) and processed the result with customized tools based on FreeSurfer (Fischl, 2012) and Pycortex (Gao et al., 2015) that were adjusted to handle our NHP data. These tools and their documentation can be found at https://github.com/VisionandCognition/NHP-Freesurfer; (Klink, 2021; copy archived at swh:1:rev:8d4b89337b865fb194e196cf1b2af4967e14d607) and https://github.com/VisionandCognition/NHP-pycortex, respectively (Messinger et al., 2021).
Electrophysiology acquisition
Request a detailed protocolNeuronal activity was acquired from 1024 channels simultaneously at a 30 kHz sampling rate. The 1024-channel pedestal was connected to eight 128-channel CerePlex M head-stages through an electronic interface board. Each head-stage processed signals from two 64-channel electrode arrays with a 0.3–7500 Hz analog filter at unity gain (i.e., no amplification). After analog-to-digital conversion, the signal from each head-stage was sent to a 128-channel Digital Hub (Blackrock Microsystems) where it was converted into an optical output signal and sent to a 128-channel Neural Signal Processor (NSP, Blackrock Microsystems) for storage and further processing. The eight NSPs were controlled with eight simultaneously running instances of the Blackrock Central Software Suite (Blackrock Microsystems) distributed over two computers (four instances each) (Chen et al., 2020).
Electrophysiology data preprocessing
Request a detailed protocolThe neuronal signal that was acquired using different software instances was temporally aligned using common TTL pulses sent by the stimulus computer. The data were then separated in (1) envelope MUA and (2) broadband LFP. MUA represents the spiking activity of a local population of neurons around the electrode (Cohen and Maunsell, 2009; Palmer et al., 2007; Supèr and Roelfsema, 2005). To extract MUA, we amplified the raw neuronal signal, band-pass filtered it between 500 Hz and 9 kHz, full-wave rectified it, and applied a low-pass filter of 200 Hz. The resulting time series were downsampled to 1 kHz. We subtracted the baseline MUA activity in a 1000 ms prestimulus time window. Baseline-corrected MUA responses were then averaged, first across runs and then within a 50–500 ms time window for each stimulus position. The broadband LFP signal was generated by low-pass filtering the raw signal at 150 Hz and downsampling it to 500 Hz. The LFP signal was further processed with a multi-taper method using the Chronux toolbox (Bokil et al., 2010). Power spectra were calculated in a 500 ms moving window (step size 50 ms), using a time bandwidth product of five and nine tapers. LFP power was averaged within five distinct frequency bands: 4–8 Hz (theta), 8–16 Hz (alpha), 16–30 Hz (beta), 30–60 Hz (low gamma), and 60–120 Hz (high gamma). Baseline power in a 1000 ms prestimulus period was subtracted and the power for each recording site was averaged across runs, within a 50–500 ms time window during each stimulus position.
pRF models and fitting procedure
Request a detailed protocolWe fit four pRF models to all the data (voxels and electrode channels) using a customized version of the analyzePRF toolbox (Kay et al., 2013) for MATLAB. In the fitting procedure, the stimuli were spatially downsampled to a resolution of 10 pixels per dva and converted to ‘effective stimuli,’ consisting of binary representations that encode stimulus position. Response predictions were calculated as the product of the effective stimulus and the pRF shape (Equation 1). This method also yields a prediction for pRFs that are only partially stimulated, and the best-fitting pRF can have a center location outside of the stimulus aperture. Because pRF estimates are generally more reliable for pRF that are strongly driven by a visual stimulus, we include indications of the directly stimulated visual field in our figures. The four pRF models differed in the pRF shape. Linear models (Equations 1–3) describe a single isotropic 2D Gaussian-shaped pRF and assume linear spatial summation across the visual field (Dumoulin and Wandell, 2008). We implemented two linear model versions. For the first model, responses were constrained to be positively related to the visual stimuli (P-LIN). A second version lacked this constraint and also allowed negative responses, that is, stimulus-driven activity reductions (U-LIN). Negative BOLD responses have been demonstrated in some brain areas (Shmuel et al., 2006; Szinte and Knapen, 2020). The nonlinear spatial summation model (Kay et al., 2013) expands the linear model by capturing nonlinear summation of signals across the visual field. It has previously been shown that the value of the exponent is generally smaller than 1 in human visual cortex, indicating subadditive spatial summation or CSS (Kay et al., 2013). This model is therefore generally referred to as the CSS model. Because the pRF size and static nonlinearity interact in the nonlinear model, the pRF size is defined as the standard deviation of the predicted Gaussian response profile to a point stimulus for all models (Equation 3; Kay et al., 2013). The mathematical descriptions of the linear and nonlinear pRF models are in Equations 1–3, where Resppred indicates the predicted response, g is a gain factor to scale the response, S(x,y) is the effective stimulus, G(x,y) is the Gaussian pRF profile, and n is the exponent that determines the static spatial nonlinearity. For the P-LIN model, the gain g was constrained to positive values, while for the U-LIN model, gain values could be negative as well. Negative gain factors imply stimulus-induced reductions of activity. In both linear models, the exponent n was fixed to be 1. In the definition of the Gaussian, (x0,y0) defines the center and σ the standard deviation of the pRF.
Linear and nonlinear spatial summation pRF models:
The DoG model uses a standard linear Gaussian (G1 in Equation 4) to describe the excitatory center of a pRF and subtracts a second Gaussian profile (G2) to model an inhibitory surround component (Zuiderbaan et al., 2012). This second Gaussian is by definition broader than the one describing the center. The size σ2 and amplitude a of the surround Gaussian are additional parameters (Equations 4–6).
DoG pRF model:
The models assume a near-instantaneous link between the stimulus and the response dynamics, which holds for the electrophysiological signals. For these signals, the stimulus changes position every 500 ms, much slower than the onset or decay of activity. The BOLD response is much slower than the speed with which the stimulus traverses the screen in the fMRI experiments (2500 ms per position). We therefore convolved the predicted response with an HRF at a resolution of 1.25 s per sample (twice the acquisition rate of TR = 2.5 s). We used both a canonical human HRF and a standard monkey HRF that we derived from separate scanning sessions (Figure 1—figure supplement 1). In short, we presented the animals with brief (0.1 s) full-contrast and full-screen checkerboard stimuli. We then used FMRIB’s Linear Optimal Basis Sets (FLOBS) (Woolrich et al., 2004) toolkit from the FSL software package to estimate the relative contributions of a set of basis functions for those voxels in the primary visual cortex that were activated by the stimulus. We then calculated a single weighted average HRF function based on these basis functions and used it as the standard monkey HRF. The monkey HRF was narrower than the canonical human HRF. It had a faster time-to-peak (4.2 s vs. 4.8 s) and peak-to-fall time (6.2 s vs. 12.6 s). The fits with the monkey HRF were slightly better, especially in the lower visual areas (Figure 1—figure supplement 1). However, size and location estimates were highly similar for the two HRFs for all models, and we only report results from model fits based on the monkey HRF.
Model fitting was performed on a cluster computer (LISA, SURFsara) using nonlinear optimization (MATLAB Optimization Toolbox). The accuracy of the different model fits was quantified as the cross-validated percentage of variance (R2) explained of the BOLD response (Equation 6: DATA) by the model prediction (Equation 6: MODEL). For cross-validation, we divided the data into two nonoverlapping sets (odd and even runs) and tested the prediction of a model that was fit one data set against the time series of the other data set and vice versa. This yielded two R2 values per voxel or electrode that were averaged. The cross-validated comparison of model performance is valid for models with different numbers of parameters and prevents overfitting. Fit results are available as voxel-based maps warped to the NMT template space on Neurovault.org (Fox et al., 2021; Gorgolewski et al., 2015) at https://identifiers.org/neurovault.collection:8082.
Model accuracy:
Comparison of pRF (and cRF) estimates
Request a detailed protocolAfter fitting the pRF models to all voxels and recording sites, we compared the pRF estimates both within and across recording modalities. We pooled voxels and recording sites across subjects and used nonparametric statistical tests (Wilcoxon signed-rank, Wilcoxon rank-sum, or Kruskal–Wallis) with post-hoc Tukey’s HSD tests to correct for multiple comparisons. For the fMRI data, we compared R2 across models and HRFs. We constructed retinotopic maps using the best pRF model and HRF and investigated the relationship between pRF eccentricity and size for a subset of ROIs with good model fits (R2 > 5%). For the electrophysiological data, we compared model accuracy and pRF estimates across MUA and LFP components to unravel to which extent retinotopic information is available in the different neuronal signals. We also compared the pRF estimates to a more conventional RF mapping technique for MUA based on responses to thin, moving bar stimuli (cRF). For recording sites with an SNR larger than 3 (i.e., visual responses that were more than three times larger than the standard deviation of the spontaneous activity), we fitted a Gaussian to the averaged MUA traces and determined the onset and offset of the visual response as the mean of this Gaussian plus or minus its standard deviation (SD). Horizontal and vertical RF boundaries were then derived from the onset and offset times for stimuli moving in opposite directions (Supèr and Roelfsema, 2005). RF centers were defined as the midpoint between the horizontal and vertical borders. For comparison with the pRFs, we calculated RF sizes as half the diagonal of the rectangular area between the horizontal and vertical cRF borders (). This measure approximates the RF radius based on the SD of a Gaussian response profile and can be directly compared to the sigma estimated by pRF models (Figure 7). It is smaller than cRF sizes typically reported in electrophysiological studies because neurons are activated by stimuli farther than 1 SD from their RF center (our lab usually defines cRF diameter as 3.3 * SD).
To compare pRFs based on fMRI-BOLD and electrophysiology, we combined data from individual animals to create one pool of BOLD-based voxel pRFs and six pools of electrophysiology-based electrode pRFs (MUA, and the different frequency bands of the LFP) for V1 and V4 data. We compared the relationship between RF eccentricity and size with a set of LMMs. We first tested for a correlation between eccentricity and size for the different signal types. Signal types with a significant correlation were subsequently tested together in a single LMM to determine whether the eccentricity-size relationship differed between signal types (interaction SIGNAL × ECC). Finally, we compared electrophysiological signals with the MRI results. For this analysis, we only selected V1 and V4 voxels with pRFs falling within the eccentricity range of the electrode arrays and voxels and electrodes with a fit accuracy above a predetermined threshold (fMRI threshold: R2 > 5%, electrophysiology threshold: R2 > 50%). To test the robustness of the results for this analysis, we repeated it with different combinations of criteria for initial voxel inclusion (i.e., including only voxels that roughly correspond to the location of the electrode arrays, or including all V1 and V4 voxels) and fit quality (fMRI threshold: R2 > 5% or 10%, electrophysiology threshold: R2 > 25% or 50%; Figure 12—figure supplement 1).
Data availability
All data and code are available on GIN: https://doi.gin.g-node.org/10.12751/g-node.2j01af. Unthresholded fMRI model fitting results are available on Neurovault: https://identifiers.org/neurovault.collection:8082.
-
NeuroVaultID 8082. Visual population receptive fields - macaque fMRI.
-
GINPopulation receptive fields in non-human primates from whole-brain fMRI and large-scale neurophysiology in visual cortex.https://doi.org/10.12751/g-node.2j01af
References
-
Visual Field Maps, Population Receptive Field Sizes, and Visual Field Coverage in the Human MT+ ComplexJournal of Neurophysiology 102:2704–2718.https://doi.org/10.1152/jn.00102.2009
-
Reward Encoding in the Monkey Anterior Cingulate CortexCerebral Cortex 16:1040–1055.https://doi.org/10.1093/cercor/bhj046
-
Visuotopic Organization of Macaque Posterior Parietal Cortex: A Functional Magnetic Resonance Imaging StudyThe Journal of Neuroscience 31:2064–2078.https://doi.org/10.1523/JNEUROSCI.3334-10.2011
-
Retinotopic Organization of Scene Areas in Macaque Inferior Temporal CortexThe Journal of Neuroscience 37:7373–7389.https://doi.org/10.1523/JNEUROSCI.0569-17.2017
-
Attention Shifts Recruit the Monkey Default Mode NetworkThe Journal of Neuroscience 38:1202–1217.https://doi.org/10.1523/JNEUROSCI.1111-17.2017
-
How well do we understand the neural origins of the fMRI BOLD signal?Trends in Neurosciences 25:27–31.https://doi.org/10.1016/s0166-2236(00)01995-0
-
fMRI and its interpretations: an illustration on directional selectivity in area V5/MTTrends in Neurosciences 31:444–453.https://doi.org/10.1016/j.tins.2008.06.004
-
Feature selectivity of the gamma-band of the local field potential in primate primary visual cortexFrontiers in Neuroscience 2:199–207.https://doi.org/10.3389/neuro.01.037.2008
-
Chronux: A platform for analyzing neural signalsJournal of Neuroscience Methods 192:146–151.https://doi.org/10.1016/j.jneumeth.2010.06.020
-
Local diversity and fine-scale organization of receptive fields in mouse visual cortexJournal of Neuroscience 31:18506–18521.https://doi.org/10.1523/JNEUROSCI.2974-11.2011
-
Visual Areas in Macaque Cortex Measured Using Functional Magnetic Resonance ImagingThe Journal of Neuroscience 22:10416–10426.
-
Spatial summation in the receptive fields of MT neuronsThe Journal of Neuroscience 19:5074–5084.https://doi.org/10.1523/JNEUROSCI.19-12-05074.1999
-
Processing of color, form and disparity information in visual areas VP and V2 of ventral extrastriate cortex in the macaque monkeyThe Journal of Neuroscience 6:2327–2351.
-
Neuronal Oscillations in Cortical NetworksScience 304:1926–1929.https://doi.org/10.1126/science.1099745
-
BookRhythms of the BrainOxford University Press.https://doi.org/10.1093/acprof:oso/9780195301069.001.0001
-
Nature and interaction of signals from the receptive field center and surround in macaque V1 neuronsJournal of Neurophysiology 88:2530–2546.https://doi.org/10.1152/jn.00692.2001
-
3D printing and modelling of customized implants and surgical guides for non-human primatesJournal of Neuroscience Methods 286:38–55.https://doi.org/10.1016/j.jneumeth.2017.05.013
-
Attention improves performance primarily by reducing interneuronal correlationsNature Neuroscience 12:1594–1600.https://doi.org/10.1038/nn.2439
-
Contralateral Visual Hemifield Representations in the Human Pulvinar NucleusJournal of Neurophysiology 98:1600–1609.https://doi.org/10.1152/jn.00419.2007
-
AFNI: Software for Analysis and Visualization of Functional Magnetic Resonance NeuroimagesComputers and Biomedical Research, an International Journal 29:162–173.https://doi.org/10.1006/cbmr.1996.0014
-
Population receptive field tuning properties of visual cortex during childhoodDevelopmental Cognitive Neuroscience 37:100614.https://doi.org/10.1016/j.dcn.2019.01.001
-
Population Receptive Field Estimation Reveals New Retinotopic Maps in Human SubcortexThe Journal of Neuroscience 35:9836–9847.https://doi.org/10.1523/JNEUROSCI.3840-14.2015
-
Vascular and neural basis of the BOLD signalCurrent Opinion in Neurobiology 58:61–69.https://doi.org/10.1016/j.conb.2019.06.004
-
Spatial spread of local field potential is band-pass in the primary visual cortexJournal of Neurophysiology 116:1986–1999.https://doi.org/10.1152/jn.00443.2016
-
How Visual Cortical Organization Is Altered by Ophthalmologic and Neurologic DisordersAnnual Review of Vision Science 4:357–379.https://doi.org/10.1146/annurev-vision-091517-033948
-
Modelling and analysis of local field potentials for studying the function of cortical circuitsNature Reviews. Neuroscience 14:770–785.https://doi.org/10.1038/nrn3599
-
Object selection by automatic spreading of top-down attentional signals in V1The Journal of Neuroscience 40:9250–9259.https://doi.org/10.1523/jneurosci.0438-20.2020
-
Receptive field properties of neurons in area V3 of macaque monkey extrastriate cortexJournal of Neurophysiology 57:889–920.https://doi.org/10.1152/jn.1987.57.4.889
-
Pycortex: an interactive surface visualizer for fMRIFrontiers in Neuroinformatics 9:23.https://doi.org/10.3389/fninf.2015.00023
-
Visual topography of V2 in the macaqueThe Journal of Comparative Neurology 201:519–539.https://doi.org/10.1002/cne.902010405
-
Visual topography of V1 in the Cebus monkeyThe Journal of Comparative Neurology 259:529–548.https://doi.org/10.1002/cne.902590404
-
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
-
Cortical visual areas in monkeys: location, topography, connections, columns, plasticity and cortical dynamicsPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 360:709–731.https://doi.org/10.1098/rstb.2005.1629
-
Neurophysiology of the BOLD fMRI signal in awake monkeysCurrent Biology 18:631–640.https://doi.org/10.1016/j.cub.2008.03.054
-
X-Chromosome Insufficiency Alters Receptive Fields across the Human Early Visual CortexThe Journal of Neuroscience 39:8079–8088.https://doi.org/10.1523/JNEUROSCI.2745-18.2019
-
The response of single optic nerve fibers of the vertebrate eye to illumination of the retinaAmerican Journal of Physiology-Legacy Content 121:400–415.https://doi.org/10.1152/ajplegacy.1938.121.2.400
-
Functional properties of monkey caudate neurons. II. Visual and auditory responsesJournal of Neurophysiology 61:799–813.https://doi.org/10.1152/jn.1989.61.4.799
-
Functional properties of monkey caudate neurons. III. Activities related to expectation of target and rewardJournal of Neurophysiology 61:814–832.https://doi.org/10.1152/jn.1989.61.4.814
-
Receptive fields of single neurones in the cat’s striate cortexThe Journal of Physiology 148:574–591.https://doi.org/10.1113/jphysiol.1959.sp006308
-
Receptive fields, binocular interaction and functional architecture in the cat’s visual cortexThe Journal of Physiology 160:106–154.https://doi.org/10.1113/jphysiol.1962.sp006837
-
Receptive fields and functional architecture of monkey striate cortexThe Journal of Physiology 195:215–243.https://doi.org/10.1113/jphysiol.1968.sp008455
-
Uniformity of monkey striate cortex: a parallel relationship between field size, scatter, and magnification factorThe Journal of Comparative Neurology 158:295–305.https://doi.org/10.1002/cne.901580305
-
Compressive spatial summation in human visual cortexJournal of Neurophysiology 110:481–494.https://doi.org/10.1152/jn.00105.2013
-
SoftwareNHP_Freesurfer, version swh:1:rev:8d4b89337b865fb194e196cf1b2af4967e14d607Software Heritage.
-
Neuronal responses to static texture patterns in area V1 of the alert macaque monkeyJournal of Neurophysiology 67:961–980.https://doi.org/10.1152/jn.1992.67.4.961
-
Visual field map clusters in macaque extrastriate visual cortexThe Journal of Neuroscience 29:7031–7039.https://doi.org/10.1523/JNEUROSCI.0518-09.2009
-
The retinotopic organization of the human middle temporal area MT/V5 and its cortical neighborsThe Journal of Neuroscience 30:9801–9820.https://doi.org/10.1523/JNEUROSCI.2069-10.2010
-
Two retinotopic visual areas in human lateral occipital cortexThe Journal of Neuroscience 26:13128–13142.https://doi.org/10.1523/JNEUROSCI.1657-06.2006
-
Population Receptive Field Shapes in Early Visual Cortex Are Nearly CircularThe Journal of Neuroscience 41:2420–2427.https://doi.org/10.1523/JNEUROSCI.3052-20.2021
-
Frequency Dependence of Signal Power and Spatial Reach of the Local Field PotentialPLOS Computational Biology 9:e1003137.https://doi.org/10.1371/journal.pcbi.1003137
-
The first step for neuroimaging data analysis: DICOM to NIfTI conversionJournal of Neuroscience Methods 264:47–56.https://doi.org/10.1016/j.jneumeth.2016.03.001
-
Stimulus-Related Neuroimaging in Task-Engaged Subjects Is Best Predicted by Concurrent SpikingThe Journal of Neuroscience 34:13878–13891.https://doi.org/10.1523/JNEUROSCI.1595-14.2014
-
The underpinnings of the BOLD functional magnetic resonance imaging signalThe Journal of Neuroscience 23:3963–3971.https://doi.org/10.1523/JNEUROSCI.23-10-03963.2003
-
Interpreting the BOLD SignalAnnual Review of Physiology 66:735–769.https://doi.org/10.1146/annurev.physiol.66.082602.092845
-
Neurovascular uncoupling: Much ado about nothingFrontiers in Neuroenergetics 2:2.https://doi.org/10.3389/fnene.2010.00002
-
Divergence of fMRI and neural signals in V1 during perceptual suppression in the awake monkeyNature Neuroscience 11:1193–1200.https://doi.org/10.1038/nn.2173
-
Default mode of brain function in monkeysThe Journal of Neuroscience 31:12954–12962.https://doi.org/10.1523/JNEUROSCI.2318-11.2011
-
Attention priority map of face images in human early visual cortexThe Journal of Neuroscience 38:149–157.https://doi.org/10.1523/JNEUROSCI.1206-17.2017
-
Do electrode properties create a problem in interpreting local field potential recordings?Journal of Neurophysiology 103:2315–2317.https://doi.org/10.1152/jn.00157.2010
-
Ventral posterior visual area of the macaque: visual topography and areal boundariesThe Journal of Comparative Neurology 252:139–153.https://doi.org/10.1002/cne.902520202
-
Spatial summation in macaque parietal area 7a follows a winner-take-all ruleJournal of Neurophysiology 105:1150–1158.https://doi.org/10.1152/jn.00907.2010
-
Linking neuronal and behavioral performance in a reaction-time visual detection taskThe Journal of Neuroscience 27:8122–8137.https://doi.org/10.1523/JNEUROSCI.1940-07.2007
-
BookAdvanced Brain Neuroimaging Topics in Health and DiseaseMethods and Applications.https://doi.org/10.5772/58256
-
Resolving the spatial profile of figure enhancement in human V1 through population receptive field modelingThe Journal of Neuroscience 40:3292–3303.https://doi.org/10.1523/JNEUROSCI.2377-19.2020
-
The relationship between subthreshold and suprathreshold ocular dominance in cat primary visual cortexThe Journal of Neuroscience 28:8553–8559.https://doi.org/10.1523/JNEUROSCI.2182-08.2008
-
Three-Dimensional Digital Template Atlas of the Macaque BrainCerebral Cortex 27:4463–4477.https://doi.org/10.1093/cercor/bhw248
-
Wide-field retinotopy reveals a new visuotopic cluster in macaque posterior parietal cortexBrain Structure & Function 225:2447–2461.https://doi.org/10.1007/s00429-020-02134-2
-
Responses of striatal neurons in the behaving monkey. 1. Head of the caudate nucleusBehavioural Brain Research 7:179–210.https://doi.org/10.1016/0166-4328(83)90191-2
-
Representation of the visual field in the second visual area in the Cebus monkeyThe Journal of Comparative Neurology 275:326–345.https://doi.org/10.1002/cne.902750303
-
“Third tier” ventral extrastriate cortex in the New World monkey, Cebus apellaExperimental Brain Research 132:287–305.https://doi.org/10.1007/s002210000344
-
Chronic recording capability of the Utah Intracortical Electrode Array in cat sensory cortexJournal of Neuroscience Methods 82:1–15.https://doi.org/10.1016/s0165-0270(98)00031-4
-
Larger Extrastriate Population Receptive Fields in Autism Spectrum DisordersThe Journal of Neuroscience 34:2713–2724.https://doi.org/10.1523/JNEUROSCI.4416-13.2014
-
Visual cortex organisation in a macaque monkey with macular degenerationThe European Journal of Neuroscience 38:3456–3464.https://doi.org/10.1111/ejn.12349
-
Apparent motion induces activity suppression in early visual cortex and impairs visual detectionThe Journal of Neuroscience 40:5471–5479.https://doi.org/10.1523/JNEUROSCI.0563-20.2020
-
Observations on the scratch-reflex in the spinal dogThe Journal of Physiology 34:1–50.https://doi.org/10.1113/jphysiol.1906.sp001139
-
Differential Sampling of Visual Space in Ventral and Dorsal Early Visual CortexThe Journal of Neuroscience 38:2294–2303.https://doi.org/10.1523/JNEUROSCI.2717-17.2018
-
Negative BOLD in the visual cortex: Evidence against blood stealingHuman Brain Mapping 21:213–220.https://doi.org/10.1002/hbm.20017
-
Chronic multiunit recordings in behaving animals: advantages and limitationsProgress in Brain Research 147:263–282.https://doi.org/10.1016/S0079-6123(04)47020-4
-
Visual Organization of the Default NetworkCerebral Cortex 30:3518–3527.https://doi.org/10.1093/cercor/bhz323
-
Mouse visual cortex contains a region of enhanced spatial resolutionNature Communications 12:4029.https://doi.org/10.1038/s41467-021-24311-5
-
Topographic Maps of Visual Space in the Human CerebellumCurrent Biology 29:1689–1694.https://doi.org/10.1016/j.cub.2019.04.012
-
Population encoding of spatial frequency, orientation, and color in macaque V1Journal of Neurophysiology 72:2151–2166.https://doi.org/10.1152/jn.1994.72.5.2151
-
Neurometabolic coupling in cerebral cortex reflects synaptic more than spiking activityNature Neuroscience 10:1308–1312.https://doi.org/10.1038/nn1977
-
Imaging retinotopic maps in the human brainVision Research 51:718–737.https://doi.org/10.1016/j.visres.2010.08.004
-
Computational neuroimaging and population receptive fieldsTrends in Cognitive Sciences 19:349–357.https://doi.org/10.1016/j.tics.2015.03.009
-
The role of neuronal synchronization in selective attentionCurrent Opinion in Neurobiology 17:154–160.https://doi.org/10.1016/j.conb.2007.02.002
-
Spatial spread of the local field potential and its laminar variation in visual cortexThe Journal of Neuroscience 29:11540–11549.https://doi.org/10.1523/JNEUROSCI.2573-09.2009
-
Receptive Fields in Human Visual Cortex Mapped with Surface ElectrodesCerebral Cortex 17:2293–2302.https://doi.org/10.1093/cercor/bhl138
Article and author information
Author details
Funding
Nederlandse Organisatie voor Wetenschappelijk Onderzoek (VENI 451.13.023)
- P Christiaan Klink
Nederlandse Organisatie voor Wetenschappelijk Onderzoek (STW-Perspectief P15-42 "NESTOR")
- Xing Chen
- Pieter Roelfsema
FP7 Ideas: European Research Council (ERC 339490 "Cortic_al_gorithms")
- Pieter R Roelfsema
Human Brain Project (Agreements 720270 and 785907 "Human Brain Project SGA1 and SGA2")
- Pieter R Roelfsema
- Wim Vanduffel
Nederlandse Organisatie voor Wetenschappelijk Onderzoek (Crossover Program 17619 "INTENSE")
- Pieter R Roelfsema
- P Christiaan Klink
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Jonathan Williford for his contributions to the fMRI preprocessing pipeline; Pieter Buur, Wietske van der Zwaag, Diederick Stoffers, and the Laboratory of Neuro- and Psychophysiology of KU Leuven for technical assistance in setting up the nonhuman primate MR infrastructure; Kor Brandsma, Anneke Ditewig, and Lex Beekman for animal care and biotechnical assistance; Feng Wang for help with electrophysiology data collection; Chris van der Togt for help with data management; and Tomas Knapen and Serge Dumoulin for fruitful discussion and comments on an earlier version of the manuscript. This work was supported by NWO (Crossover Program 17619 'INTENSE'; STW-Perspectief P15-42 'NESTOR'; VENI 451.13.023), the European Union FP7 (ERC 339490 'Cortic_al_gorithms'), the Human Brain Project (agreements 720270 and 785907, 'Human Brain Project SGA1 and SGA2'), and the Friends Foundation of the Netherlands Institute for Neuroscience.
Ethics
Animal care and experimental procedures were in accordance with the ILAR's Guide for the Care and Use of Laboratory Animals, the European legislation (Directive 2010/63/EU) and approved by the institutional animal care and use committee of the Royal Netherlands Academy of Arts and Sciences and the Central Authority for Scientific Procedures on Animals (CCD) in the Netherlands (License numbers AVD8010020173789 and AVD8010020171046).
Copyright
© 2021, Klink 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
-
- 2,814
- views
-
- 390
- downloads
-
- 39
- 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
Specialized chemosensory signals elicit innate social behaviors in individuals of several vertebrate species, a process that is mediated via the accessory olfactory system (AOS). The AOS comprising the peripheral sensory vomeronasal organ has evolved elaborate molecular and cellular mechanisms to detect chemo signals. To gain insight into the cell types, developmental gene expression patterns, and functional differences amongst neurons, we performed single-cell transcriptomics of the mouse vomeronasal sensory epithelium. Our analysis reveals diverse cell types with gene expression patterns specific to each, which we made available as a searchable web resource accessed from https://www.scvnoexplorer.com. Pseudo-time developmental analysis indicates that neurons originating from common progenitors diverge in their gene expression during maturation with transient and persistent transcription factor expression at critical branch points. Comparative analysis across two of the major neuronal subtypes that express divergent GPCR families and the G-protein subunits Gnai2 or Gnao1, reveals significantly higher expression of endoplasmic reticulum (ER) associated genes within Gnao1 neurons. In addition, differences in ER content and prevalence of cubic membrane ER ultrastructure revealed by electron microscopy, indicate fundamental differences in ER function.
-
- Developmental Biology
- Neuroscience
Williams syndrome (WS; OMIM#194050) is a rare disorder, which is caused by the microdeletion of one copy of 25–27 genes, and WS patients display diverse neuronal deficits. Although remarkable progresses have been achieved, the mechanisms for these distinct deficits are still largely unknown. Here, we have shown that neural progenitor cells (NPCs) in WS forebrain organoids display abnormal proliferation and differentiation capabilities, and synapse formation. Genes with altered expression are related to neuronal development and neurogenesis. Single cell RNA-seq (scRNA-seq) data analysis revealed 13 clusters in healthy control and WS organoids. WS organoids show an aberrant generation of excitatory neurons. Mechanistically, the expression of transthyretin (TTR) are remarkably decreased in WS forebrain organoids. We have found that GTF2IRD1 encoded by one WS associated gene GTF2IRD1 binds to TTR promoter regions and regulates the expression of TTR. In addition, exogenous TTR can activate ERK signaling and rescue neurogenic deficits of WS forebrain organoids. Gtf2ird1-deficient mice display similar neurodevelopmental deficits as observed in WS organoids. Collectively, our study reveals critical function of GTF2IRD1 in regulating neurodevelopment of WS forebrain organoids and mice through regulating TTR-ERK pathway.