Mapping sequences can bias population receptive field estimates

Improving the Reliability and Accuracy of Population Receptive Field Measures Using a Logarithmically Warped Stimulus

The population receptive field method, which measures the region in visual space that elicits a BOLD signal in a voxel in retinotopic cortex, is a powerful tool for investigating the functional organization of human visual cortex with fMRI (Dumoulin & Wandell, 2008). However, recent work has shown that population receptive field (pRF) estimates for early retinotopic visual areas can be biased and unreliable, especially for voxels representing the fovea. Here, we show that a 'log-bar' stimulus that is logarithmically warped along the eccentricity dimension produces more reliable estimates of pRF size and location than the traditional moving bar stimulus. The log-bar stimulus was better able to identify pRFs near the foveal representation, and pRFs were smaller in size, consistent with simulation estimates of receptive field sizes in the fovea.

Comparison of Stimulus Types for Retinotopic Cortical Mapping of Macular Disease

Purpose

Retinotopic maps acquired using functional magnetic resonance imaging (fMRI) provide a valuable adjunct in the assessment of macular function at the level of the visual cortex. The present study quantitatively assessed the performance of different visual stimulation approaches for mapping visual field coverage.

Methods

Twelve patients with geographic atrophy (GA) secondary to age-related macular degeneration (AMD) were examined using high-resolution ultra-high field fMRI (Siemens Magnetom 7T) and microperimetry (MP; Nidek MP-3). The population receptive field (pRF)-based coverage maps obtained with two different stimulus techniques (moving bars, and rotating wedges and expanding rings) were compared with the results of MP. Correspondence between MP and pRF mapping was quantified by calculating the simple matching coefficient (SMC).

Results

Stimulus choice is shown to bias the spatial distribution of pRF centers and eccentricity values with pRF sizes obtained from wedge/ring or bar stimulation showing systematic differences. Wedge/ring stimulation results show a higher number of pRF centers in foveal areas and strongly reduced pRF sizes compared to bar stimulation runs. A statistical comparison shows significantly higher pRF center numbers in the foveal 2.5 degrees region of the visual field for wedge/ring compared to bar stimuli. However, these differences do not significantly influence SMC values when compared to MP (bar <2.5 degrees: 0.88 ± 0.13; bar >2.5 degrees: 0.88 ± 0.11; wedge/ring <2.5 degrees: 0.89 ± 0.12 wedge/ring; >2.5 degrees: 0.86 ± 0.10) for the peripheral visual field.

Conclusions

Both visual stimulation designs examined can be applied successfully in patients with GA. Although the two designs show systematic differences in the distribution of pRF center locations, this variability has minimal impact on the SMC when compared to the MP outcome.

Pitfalls in Post Hoc Analyses of Population Receptive Field Data

Data binning involves grouping observations into bins and calculating bin-wise summary statistics. It can cope with overplotting and noise, making it a versatile tool for comparing many observations. However, data binning goes awry if the same observations are used for binning (selection) and contrasting (selective analysis). This creates circularity, biasing noise components and resulting in artifactual changes in the form of regression towards the mean. Importantly, these artifactual changes are a statistical necessity. Here, we use (null) simulations and empirical repeat data to expose this flaw in the scope of post hoc analyses of population receptive field data. In doing so, we reveal that the type of data analysis, data properties, and circular data cleaning are factors shaping the appearance of such artifactual changes. We furthermore highlight that circular data cleaning and circular sorting of change scores are selection practices that result in artifactual changes even without circular data binning. These pitfalls might have led to erroneous claims about changes in population receptive fields in previous work and can be mitigated by using independent data for selection purposes. Our evaluations highlight the urgency for us researchers to make the validation of analysis pipelines standard practice.

Highly accurate retinotopic maps of the physiological blind spot in human visual cortex

The physiological blind spot is a naturally occurring scotoma corresponding with the optic disc in the retina of each eye. Even during monocular viewing, observers are usually oblivious to the scotoma, in part because the visual system extrapolates information from the surrounding area. Unfortunately, studying this visual field region with neuroimaging has proven difficult, as it occupies only a small part of retinotopic cortex. Here, we used functional magnetic resonance imaging and a novel data-driven method for mapping the retinotopic organization in and around the blind spot representation in V1. Our approach allowed for highly accurate reconstructions of the extent of an observer's blind spot, and out-performed conventional model-based analyses. This method opens exciting opportunities to study the plasticity of receptive fields after visual field loss, and our data add to evidence suggesting that the neural circuitry responsible for impressions of perceptual completion across the physiological blind spot most likely involves regions of extrastriate cortex-beyond V1.

Cortical depth-dependent modeling of visual hemodynamic responses

A physiologically based three-dimensional (3D) hemodynamic model is developed to predict the experimentally observed blood oxygen level dependent (BOLD) responses versus the cortical depth induced by visual stimuli. Prior 2D approximations are relaxed in order to analyze 3D blood flow dynamics as a function of cortical depth. Comparison of the predictions with experimental data for evoked stimuli demonstrates that the full 3D model performs at least as well as previous approaches while remaining parsimonious. In particular, the 3D model requires significantly fewer assumptions and model parameters than previous models such that there is no longer need to define depth-specific parameter values for spatial spreading, peak amplitude, and hemodynamic velocity.

Visual Field Reconstruction in Hemianopia Using fMRI Based Mapping Techniques

PURPOSE

A stroke that includes the primary visual cortex unilaterally leads to a loss of visual field (VF) representation in the hemifield contralateral to the damage. While behavioral procedures for measuring the VF, such as perimetry, may indicate that a patient cannot see in a particular area, detailed psychophysical testing often detects the ability to perform detection or discrimination of visual stimuli ("blindsight"). The aim of this study was to determine whether functional magnetic resonance imaging (fMRI) could be used to determine whether perimetrically blind regions of the VF were still represented in VF maps reconstructed on the basis of visually evoked neural activity.

METHODS

Thirteen patients with hemianopia and nine control participants were scanned using 3T MRI while presented with visual stimulation. Two runs of a dynamic "wedge and ring" mapping stimulus, totaling approximately 10 min, were performed while participants fixated centrally. Two different analysis approaches were taken: the conventional population receptive field (pRF) analysis and micro-probing (MP). The latter is a variant of the former that makes fewer assumptions when modeling the visually evoked neural activity. Both methods were used to reconstruct the VF by projecting modeled activity back onto the VF. Following a normalization step, these "coverage maps" can be compared to the VF sensitivity plots obtained using perimetry.

RESULTS

While both fMRI-based approaches revealed regions of neural activity within the perimetrically "blind" sections of the VF, the MP approach uncovered more voxels in the lesioned hemisphere in which a modest degree of visual sensitivity was retained. Furthermore, MP-based analysis indicated that both early (V1/V2) and extrastriate visual areas contributed equally to the retained sensitivity in both patients and controls.

CONCLUSION

In hemianopic patients, fMRI-based approaches for reconstructing the VF can pick up activity in perimetrically blind regions of the VF. Such regions of the VF may be particularly amenable for rehabilitation to regain visual function. Compared to conventional pRF modeling, MP reveals more voxels with retained visual sensitivity, suggesting it is a more sensitive approach for VF reconstruction.

Inferior Occipital Gyrus Is Organized along Common Gradients of Spatial and Face-Part Selectivity

The ventral visual stream of the human brain is subdivided into patches with categorical stimulus preferences, like faces or scenes. However, the functional organization within these areas is less clear. Here, we used functional magnetic resonance imaging and vertex-wise tuning models to independently probe spatial and face-part preferences in the inferior occipital gyrus (IOG) of healthy adult males and females. The majority of responses were well explained by Gaussian population tuning curves for both retinotopic location and the preferred relative position within a face. Parameter maps revealed a common gradient of spatial and face-part selectivity, with the width of tuning curves drastically increasing from posterior to anterior IOG. Tuning peaks clustered more idiosyncratically but were also correlated across maps of visual and face space. Preferences for the upper visual field went along with significantly increased coverage of the upper half of the face, matching recently discovered biases in human perception. Our findings reveal a broad range of neural face-part selectivity in IOG, ranging from narrow to "holistic." IOG is functionally organized along this gradient, which in turn is correlated with retinotopy.SIGNIFICANCE STATEMENT Brain imaging has revealed a lot about the large-scale organization of the human brain and visual system. For example, occipital cortex contains map-like representations of the visual field, while neurons in ventral areas cluster into patches with categorical preferences, like faces or scenes. Much less is known about the functional organization within these areas. Here, we focused on a well established face-preferring area-the inferior occipital gyrus (IOG). A novel neuroimaging paradigm allowed us to map the retinotopic and face-part tuning of many recording sites in IOG independently. We found a steep posterior-anterior gradient of decreasing face-part selectivity, which correlated with retinotopy. This suggests the functional role of ventral areas is not uniform and may follow retinotopic "protomaps."

Combining stimulus types for improved coverage in population receptive field mapping

Retinotopy experiments using population receptive field (pRF) mapping are ideal for assigning regions in the visual field to cortical brain areas. While various designs for visual stimulation were suggested in the literature, all have specific shortcomings regarding visual field coverage. Here we acquired high-resolution 7 Tesla fMRI data to compare pRF-based coverage maps obtained with the two most commonly used stimulus variants: moving bars; rotating wedges and expanding rings. We find that stimulus selection biases the spatial distribution of pRF centres. In addition, eccentricity values and pRF sizes obtained from wedge/ring or bar stimulation runs show systematic differences. Wedge/ring stimulation results show lower eccentricity values and strongly reduced pRF sizes compared to bar stimulation runs. Statistical comparison shows significantly higher pRF centre numbers in the foveal 2° region of the visual field for wedge/ring compared to bar stimuli. We suggest and evaluate approaches for combining pRF data from different visual stimulus patterns to obtain improved mapping results.

A validation framework for neuroimaging software: The case of population receptive fields

Neuroimaging software methods are complex, making it a near certainty that some implementations will contain errors. Modern computational techniques (i.e., public code and data repositories, continuous integration, containerization) enable the reproducibility of the analyses and reduce coding errors, but they do not guarantee the scientific validity of the results. It is difficult, nay impossible, for researchers to check the accuracy of software by reading the source code; ground truth test datasets are needed. Computational reproducibility means providing software so that for the same input anyone obtains the same result, right or wrong. Computational validity means obtaining the right result for the ground-truth test data. We describe a framework for validating and sharing software implementations, and we illustrate its usage with an example application: population receptive field (pRF) methods for functional MRI data. The framework is composed of three main components implemented with containerization methods to guarantee computational reproducibility. In our example pRF application, those components are: (1) synthesis of fMRI time series from ground-truth pRF parameters, (2) implementation of four public pRF analysis tools and standardization of inputs and outputs, and (3) report creation to compare the results with the ground truth parameters. The framework was useful in identifying realistic conditions that lead to imperfect parameter recovery in all four pRF implementations, that would remain undetected using classic validation methods. We provide means to mitigate these problems in future experiments. A computational validation framework supports scientific rigor and creativity, as opposed to the oft-repeated suggestion that investigators rely upon a few agreed upon packages. We hope that the framework will be helpful to validate other critical neuroimaging algorithms, as having a validation framework helps (1) developers to build new software, (2) research scientists to verify the software's accuracy, and (3) reviewers to evaluate the methods used in publications and grants.

Comparison of human population receptive field estimates between scanners and the effect of temporal filtering

Background: Population receptive field (pRF) analysis with functional magnetic resonance imaging (fMRI) is an increasingly popular method for mapping visual field representations and estimating the spatial selectivity of voxels in human visual cortex. However, the multitude of experimental setups and processing methods used makes comparisons of results between studies difficult. Methods: Here, we compared pRF maps acquired in the same three individuals using comparable scanning parameters on a 1.5 and a 3 Tesla scanner located in two different countries. We also tested the effect of low-pass filtering of the time series on pRF estimates. Results: As expected, the signal-to-noise ratio for the 3 Tesla data was superior; critically, however, estimates of pRF size and cortical magnification did not reveal any systematic differences between the sites. Unsurprisingly, low-pass filtering enhanced goodness-of-fit, presumably by removing high-frequency noise. However, there was no substantial increase in the number of voxels containing meaningful retinotopic signals after low-pass filtering. Importantly, filtering also increased estimates of pRF size in the early visual areas which could substantially skew interpretations of spatial tuning properties. Conclusion: Our results therefore suggest that pRF estimates are generally comparable between scanners of different field strengths, but temporal filtering should be used with caution.