Simultaneous estimation of population receptive field and hemodynamic parameters from single point BOLD responses using Metropolis-Hastings sampling

Bayesian connective field modeling using a Markov Chain Monte Carlo approach

The majority of neurons in the human brain process signals from neurons elsewhere in the brain. Connective Field (CF) modelling is a biologically-grounded method to describe this essential aspect of the brain's circuitry. It allows characterizing the response of a population of neurons in terms of the activity in another part of the brain. CF modelling translates the concept of the receptive field (RF) into the domain of connectivity by assessing, at the voxel level, the spatial dependency between signals in distinct cortical visual field areas. Thus, the approach enables to characterize the functional cortical circuitry of the human cortex. While already very useful, the present CF modelling approach has some intrinsic limitations due to the fact that it only estimates the model's explained variance and not the probability distribution associated with the estimated parameters. If we could resolve this, CF modelling would lend itself much better for statistical comparisons at the level of single voxels and individuals. This is important when trying to gain a detailed understanding of the neurobiology and pathophysiology of the visual cortex, notably in rare cases. To enable this, we present a Bayesian approach to CF modeling (bCF). Using a Markov Chain Monte Carlo (MCMC) procedure, it estimates the posterior probability distribution underlying the CF parameters. Based on this, bCF quantifies, at the voxel level, the uncertainty associated with each parameter estimate. This information can be used in various ways to increase confidence in the CF model predictions. We applied bCF to BOLD responses recorded in the early human visual cortex using 3T fMRI. We estimated both the CF parameters and their associated uncertainties and show they are only weakly correlated. Moreover, we show how bCF facilitates the use of effect size (beta) as a data-driven parameter that can be used to select the most reliable voxels for further analysis. Finally, to further illustrate the functionality gained by bCF, we apply it to perform a voxel-level comparison of a single, circular symmetric, Gaussian versus a Difference-of-Gaussian model. We conclude that our bCF framework provides a comprehensive tool to study human functional cortical circuitry in health and disease.

Foveal pRF properties in the visual cortex depend on the extent of stimulated visual field

Previous studies demonstrated that alterations in functional MRI derived receptive field (pRF) properties in cortical projection zones of retinal lesions can erroneously be mistaken for cortical large-scale reorganization in response to visual system pathologies. We tested, whether such confounds are also evident in the normal cortical projection zone of the fovea for simulated peripheral visual field defects. We applied fMRI-based visual field mapping of the central visual field at 3 T in eight controls to compare the pRF properties of the central visual field of a reference condition (stimulus radius: 14°) and two conditions with simulated peripheral visual field defect, i.e., with a peripheral gray mask, stimulating only the central 7° or 4° radius. We quantified, for the cortical representation of the actually stimulated visual field, the changes in the position and size of the pRFs associated with reduced peripheral stimulation using conventional and advanced pRF modeling. We found foveal pRF-positions (≤3°) to be significantly shifted towards the periphery (p<0.05, corrected). These pRF-shifts were largest for the 4° condition [visual area (mean eccentricity shift): V1 (0.9°), V2 (0.9°), V3 (1.0°)], but also evident for the 7° condition [V1 (0.5°), V2 (0.5°), V3 (0.9°)]. Further, an overall enlargement of pRF-sizes was observed. These findings indicate the dependence of foveal pRF parameters on the spatial extent of the stimulated visual field and are likely associated with methodological biases and/or physiological mechanisms. Consequently, our results imply that, previously reported similar findings in patients with actual peripheral scotomas need to be interpreted with caution and indicate the need for adequate control conditions in investigations of visual cortex reorganization.

Studying Cortical Plasticity in Ophthalmic and Neurological Disorders: From Stimulus-Driven to Cortical Circuitry Modeling Approaches

Unsolved questions in computational visual neuroscience research are whether and how neurons and their connecting cortical networks can adapt when normal vision is compromised by a neurodevelopmental disorder or damage to the visual system. This question on neuroplasticity is particularly relevant in the context of rehabilitation therapies that attempt to overcome limitations or damage, through either perceptual training or retinal and cortical implants. Studies on cortical neuroplasticity have generally made the assumption that neuronal population properties and the resulting visual field maps are stable in healthy observers. Consequently, differences in the estimates of these properties between patients and healthy observers have been taken as a straightforward indication for neuroplasticity. However, recent studies imply that the modeled neuronal properties and the cortical visual maps vary substantially within healthy participants, e.g., in response to specific stimuli or under the influence of cognitive factors such as attention. Although notable advances have been made to improve the reliability of stimulus-driven approaches, the reliance on the visual input remains a challenge for the interpretability of the obtained results. Therefore, we argue that there is an important role in the study of cortical neuroplasticity for approaches that assess intracortical signal processing and circuitry models that can link visual cortex anatomy, function, and dynamics.

Micro-probing enables fine-grained mapping of neuronal populations using fMRI

The characterization of receptive field (RF) properties is fundamental to understanding the neural basis of sensory and cognitive behaviour. The combination of non-invasive imaging, such as fMRI, with biologically inspired neural modelling has enabled the estimation of population RFs directly in humans. However, current approaches require making numerous a priori assumptions, so these cannot reveal unpredicted properties, such as fragmented RFs or subpopulations. This is a critical limitation in studies on adaptation, pathology or reorganization. Here, we introduce micro-probing (MP), a technique for fine-grained and largely assumption free characterization of multiple pRFs within a voxel. It overcomes many limitations of current approaches by enabling detection of unexpected RF shapes, properties and subpopulations, by enhancing the spatial detail with which we analyze the data. MP is based on tiny, fixed-size, Gaussian models that efficiently sample the entire visual space and create fine-grained probe maps. Subsequently, we derived population receptive fields (pRFs) from these maps. We demonstrate the scope of our method through simulations and by mapping the visual fields of healthy participants and of a patient group with highly abnormal RFs due to a congenital pathway disorder. Without using specific stimuli or adapted models, MP mapped the bilateral pRFs characteristic of observers with albinism. In healthy observers, MP revealed that voxels may capture the activity of multiple subpopulations RFs that sample distinct regions of the visual field. Thus, MP provides a versatile framework to visualize, analyze and model, without restrictions, the diverse RFs of cortical subpopulations in health and disease.

Population receptive field and connectivity properties of the early visual cortex in human albinism

In albinism, the pathological decussation of the temporal retinal afferents at the optic chiasm leads to superimposed representations of opposing hemifields in the visual cortex. Here, we assessed the equivalence of the two representations and the cortico-cortical connectivity of the early visual areas. Applying fMRI-based population receptive field (pRF)-mapping (both hemifield and bilateral mapping) and connective field (CF)-modeling, we investigated the early visual cortex in 6 albinotic participants and 4 controls. In albinism, superimposed retinotopic representations of the contra- and ipsilateral visual hemifield were observed on the hemisphere contralateral to the stimulated eye. This was confirmed by the observation of bilateral pRFs during bilateral mapping. Hemifield mapping revealed similar pRF-sizes for both hemifield representations throughout V1 to V3. The typical increase of V1-sampling extent for V3 compared to V2 was not found for the albinotic participants. The similarity of the pRF-sizes for opposing visual hemifield representations highlights the equivalence of the two maps in the early visual cortex. The altered V1-sampling extent in V3 might indicate the adaptation of cortico-cortical connections to visual pathway abnormalities in albinism. These findings thus suggest that conservative developmental mechanisms are complemented by alterations of the extrastriate cortico-cortical connectivity.

Methods for computing the maximum performance of computational models of fMRI responses

Computational neuroimaging methods aim to predict brain responses (measured e.g. with functional magnetic resonance imaging [fMRI]) on the basis of stimulus features obtained through computational models. The accuracy of such prediction is used as an indicator of how well the model describes the computations underlying the brain function that is being considered. However, the prediction accuracy is bounded by the proportion of the variance of the brain response which is related to the measurement noise and not to the stimuli (or cognitive functions). This bound to the performance of a computational model has been referred to as the noise ceiling. In previous fMRI applications two methods have been proposed to estimate the noise ceiling based on either a split-half procedure or Monte Carlo simulations. These methods make different assumptions over the nature of the effects underlying the data, and, importantly, their relation has not been clarified yet. Here, we derive an analytical form for the noise ceiling that does not require computationally expensive simulations or a splitting procedure that reduce the amount of data. The validity of this analytical definition is proved in simulations, we show that the analytical solution results in the same estimate of the noise ceiling as the Monte Carlo method. Considering different simulated noise structure, we evaluate different estimators of the variance of the responses and their impact on the estimation of the noise ceiling. We furthermore evaluate the interplay between regularization (often used to estimate model fits to the data when the number of computational features in the model is large) and model complexity on the performance with respect to the noise ceiling. Our results indicate that when considering the variance of the responses across runs, computing the noise ceiling analytically results in similar estimates as the split half estimator and approaches the true noise ceiling under a variety of simulated noise scenarios. Finally, the methods are tested on real fMRI data acquired at 7 Tesla.

Encoding computational models in brain responses measured with fMRI allows testing the algorithmic representations carried out by the neural population within voxels. The accuracy of a model in predicting new responses is used as a measure of the brain validity of the computational model being tested, but the result of this analysis is determined not only by how precisely the model describes the responses but also by the quality of the data. In this article, we evaluate existing approaches to estimate the best possible accuracy that any computational model can achieve conditioned to the amount of measurement noise that is present in the experimental data (i.e. the noise ceiling). Additionally we introduce a close form estimation of the noise ceiling that does not require computationally or data expensive procedures. All the methods are compared using simulated and real fMRI data. We draw conclusions over the impact of regularization procedures and make practical recommendations on how to report the results of computational models in neuroimaging.