Bimodal modulation and continuous stimulation in optical imaging to map direction selectivity

Refractory density model of cortical direction selectivity: Lagged-nonlagged, transient-sustained, and On-Off thalamic neuron-based mechanisms and intracortical amplification

A biophysically detailed description of the mechanisms of the primary vision is still being developed. We have incorporated a simplified, filter-based description of retino-thalamic visual signal processing into the detailed, conductance-based refractory density description of the neuronal population activity of the primary visual cortex. We compared four mechanisms of the direction selectivity (DS), three of them being based on asymmetrical projections of different types of thalamic neurons to the cortex, distinguishing between (i) lagged and nonlagged, (ii) transient and sustained, and (iii) On and Off neurons. The fourth mechanism implies a lack of subcortical bias and is an epiphenomenon of intracortical interactions between orientation columns. The simulations of the cortical response to moving gratings have verified that first three mechanisms provide DS to an extent compared with experimental data and that the biophysical model realistically reproduces characteristics of the visual cortex activity, such as membrane potential, firing rate, and synaptic conductances. The proposed model reveals the difference between the mechanisms of both the intact and the silenced cortex, favoring the second mechanism. In the fourth case, DS is weaker but significant; it completely vanishes in the silenced cortex.DS in the On-Off mechanism derives from the nonlinear interactions within the orientation map. Results of simulations can help to identify a prevailing mechanism of DS in V1. This is a step towards a comprehensive biophysical modeling of the primary visual system in the frameworks of the population rate coding concept.

A major mechanism that underlies tuning of cortical neurons to the direction of a moving stimulus is still debated. Considering the visual cortex structured with orientation-selective columns, we have realized and compared in our biophysically detailed mathematical model four hypothetical mechanisms of the direction selectivity (DS) known from experiments. The present model accomplishes our previous model that was tuned to experimental data on excitability in slices and reproduces orientation tuning effects in vivo. In simulations, we have found that the convergence of inputs from so-called transient and sustained (or lagged and nonlagged) thalamic neurons in the cortex provides an initial bias for DS, whereas cortical interactions amplify the tuning. In the absence of any bias, DS emerges as an epiphenomenon of the orientation map. In the case of a biased convergence of On- and Off- thalamic inputs, DS emerges with the help of the intracortical interactions on the orientation map, also. Thus, we have proposed a comprehensive description of the primary vision and revealed characteristic features of different mechanisms of DS in the visual cortex with columnar structure.

The hemodynamic signal as a first-order low-pass temporal filter: Evidence and implications for neuroimaging studies

Neuronal activation triggers local changes in blood flow and hemoglobin oxygenation. These hemodynamic signals can be recorded through functional magnetic resonance imaging or intrinsic optical imaging, and allows inferring neural activity in response to stimuli. These techniques are widely used to uncover functional brain architectures. However, their accuracy suffers from distortions inherent to hemodynamic responses and noise. The analysis of these signals currently relies on models of impulse hemodynamic responses to brief stimuli. Here, in order to infer precise functional architectures, we focused on integrated signals associated to the dynamic response of functional maps. To this end, we recorded orientation and direction maps in cat primary visual cortex and compared two protocols: the conventional episodic stimulation technique and a continuous, periodic stimulation paradigm. Conventional methods show that the dynamics of activation and deactivation of the functional maps follows a linear first-order differential equation representing a low-pass filter. Comparison with the periodic stimulation methods confirmed this observation: the phase shifts and magnitude attenuations extracted at various frequencies were consistent with a low-pass filter with a 5s time constant. This dynamics presumably reflects the variations in deoxyhemoglobin mediated by arterial dilations. This dynamics open new avenues in the analysis of neuroimaging data that differs from common methods based on the hemodynamic response function. In particular, we demonstrate that inverting this first-order low-pass filter minimized the distortions of the signal and enabled a much faster and accurate reconstruction of functional maps.

Mesoscale Mapping of Mouse Cortex Reveals Frequency-Dependent Cycling between Distinct Macroscale Functional Modules

Connectivity mapping based on resting-state activity in mice has revealed functional motifs of correlated activity. However, the rules by which motifs organize into larger functional modules that lead to hemisphere wide spatial-temporal activity sequences is not clear. We explore cortical activity parcellation in head-fixed, quiet awake GCaMP6 mice from both sexes by using mesoscopic calcium imaging. Spectral decomposition of spontaneous cortical activity revealed the presence of two dominant frequency modes (<1 and ∼3 Hz), each of them associated with a unique spatial signature of cortical macro-parcellation not predicted by classical cytoarchitectonic definitions of cortical areas. Based on assessment of 0.1-1 Hz activity, we define two macro-organizing principles: the first being a rotating polymodal-association pinwheel structure around which activity flows sequentially from visual to barrel then to hindlimb somatosensory; the second principle is correlated activity symmetry planes that exist on many levels within a single domain such as intrahemispheric reflections of sensory and motor cortices. In contrast, higher frequency activity >1 Hz yielded two larger clusters of coactivated areas with an enlarged default mode network-like posterior region. We suggest that the apparent constrained structure for intra-areal cortical activity flow could be exploited in future efforts to normalize activity in diseases of the nervous system.SIGNIFICANCE STATEMENT Increasingly, functional connectivity mapping of spontaneous activity is being used to reveal the organization of the brain. However, because the brain operates across multiple space and time domains a more detailed understanding of this organization is necessary. We used in vivo wide-field calcium imaging of the indicator GCaMP6 in head-fixed, awake mice to characterize the organization of spontaneous cortical activity at different spatiotemporal scales. Correlation analysis defines the presence of two to three superclusters of activity that span traditionally defined functional territories and were frequency dependent. This work helps define the rules for how different cortical areas interact in time and space. We provide a framework necessary for future studies that explore functional reorganization of brain circuits in disease models.

3D functional ultrasound imaging of the cerebral visual system in rodents

3D functional imaging of the whole brain activity during visual task is a challenging task in rodents due to the complex tri-dimensional shape of involved brain regions and the fine spatial and temporal resolutions required to reveal the visual tract. By coupling functional ultrasound (fUS) imaging with a translational motorized stage and an episodic visual stimulation device, we managed to accurately map and to recover the activity of the visual cortices, the Superior Colliculus (SC) and the Lateral Geniculate Nuclei (LGN) in 3D. Cerebral Blood Volume (CBV) responses during visual stimuli were found to be highly correlated with the visual stimulus time profile in visual cortices (r=0.6), SC (r=0.7) and LGN (r=0.7). These responses were found dependent on flickering frequency and contrast, and optimal stimulus parameters for largest CBV increases were obtained. In particular, increasing the flickering frequency higher than 7 Hz revealed a decrease of visual cortices response while the SC response was preserved. Finally, cross-correlation between CBV signals exhibited significant delays (d=0.35 s +/−0.1 s) between blood volume response in SC and visual cortices in response to our visual stimulus. These results emphasize the interest of fUS imaging as a whole brain neuroimaging modality for brain vision studies in rodent models.

Pinwheel-dipole configuration in cat early visual cortex

In the early visual cortex, information is processed within functional maps whose layouts are thought to underlie visual perception. However, the precise organization of these functional maps as well as their interrelationships remain unsettled. Here, we show that spatial frequency representation in cat early visual cortex exhibits singularities around which the map organizes like an electric dipole potential. These singularities are precisely co-located with singularities of the orientation map: the pinwheel centers. To show this, we used high resolution intrinsic optical imaging in cat areas 17 and 18. First, we show that a majority of pinwheel centers exhibit in their neighborhood both semi-global maximum and minimum in the spatial frequency map (i.e. extreme values of the spatial frequency in a hypercolumn). This contradicts pioneering studies suggesting that pinwheel centers are placed at the locus of a single spatial frequency extremum. Based on an analogy with electromagnetism, we proposed a mathematical model for a dipolar structure, accurately fitting optical imaging data. We conclude that a majority of orientation pinwheel centers form spatial frequency dipoles in cat early visual cortex. Given the functional specificities of neurons at singularities in the visual cortex, it is argued that the dipolar organization of spatial frequency around pinwheel centers could be fundamental for visual processing.

Mesoscale transcranial spontaneous activity mapping in GCaMP3 transgenic mice reveals extensive reciprocal connections between areas of somatomotor cortex

Transgenic mice expressing genetically encoded activity indicators are an attractive means of mapping mesoscopic regional functional cortical connectivity given widespread stable and cell-specific expression compatible with chronic recordings. Cortical functional connectivity was evaluated using wide-field imaging in lightly anesthetized Emx1-creXRosa26-GCaMP3 mice expressing calcium sensor in cortical neurons. Challenges exist because green fluorescence signals overlap with endogenous activity-dependent autofluorescence and are affected by changes in blood volume and oxygenation. Under the conditions used for imaging and analysis (0.1-1 Hz frequency band), autofluorescence and hemodynamic effects contributed 3% and 8% of the SD of spontaneous activity-dependent GCaMP3 fluorescence when signals were recorded through intact bone. To evaluate the accuracy and sensitivity of this approach, the topology of functional connections between somatomotor cortex (primary S1 and secondary S2 somatosensory, and primary motor cortex M1) was estimated. During sequences of spontaneous activity, calcium signals recorded at each location of area S1 were correlated with activity in contralateral area S1, ipsilateral area S2, and bilateral areas M1. Reciprocal results were observed when "seed pixels" were placed in S2 and M1. Coactivation of areas implies functional connections but could also be attributed to both regions receiving common upstream drive. These apparent connections revealed during spontaneous activity coactivation by GCaMP3 were confirmed by intracortical microstimulation but were more difficult to detect using intrinsic signals from reflected red light. We anticipate GCAMP wide-field imaging will enable longitudinal studies during plasticity paradigms or after models of CNS disease, such as stroke, where the weighting within these connectivity maps may be altered.

Vobi One: a data processing software package for functional optical imaging

Optical imaging is the only technique that allows to record the activity of a neuronal population at the mesoscopic scale. A large region of the cortex (10-20 mm diameter) is directly imaged with a CCD camera while the animal performs a behavioral task, producing spatio-temporal data with an unprecedented combination of spatial and temporal resolutions (respectively, tens of micrometers and milliseconds). However, researchers who have developed and used this technique have relied on heterogeneous software and methods to analyze their data. In this paper, we introduce Vobi One, a software package entirely dedicated to the processing of functional optical imaging data. It has been designed to facilitate the processing of data and the comparison of different analysis methods. Moreover, it should help bring good analysis practices to the community because it relies on a database and a standard format for data handling and it provides tools that allow producing reproducible research. Vobi One is an extension of the BrainVISA software platform, entirely written with the Python programming language, open source and freely available for download at https://trac.int.univ-amu.fr/vobi_one.

Organization and origin of spatial frequency maps in cat visual cortex

It remains controversial whether and how spatial frequency (SF) is represented tangentially in cat visual cortex. Several models were proposed, but there is no consensus. Worse still, some data indicate that the SF organization previously revealed by optical imaging techniques simply reflects non-stimulus-specific responses. Instead, stimulus-specific responses arise from the homogeneous distribution of geniculo-cortical afferents representing X and Y pathways. To clarify this, we developed a new imaging method allowing rapid stimulation with a wide range of SFs covering more than 6 octaves with only 0.2 octave resolution. A benefit of this method is to avoid error of high-pass filtering methods which systematically under-represent dominant selectivity features near pinwheel centers. We show unequivocally that SF is organized into maps in cat area 17 (A17) and area 18 (A18). The SF organization in each area displays a global anteroposterior SF gradient and local patches. Its layout is constrained to that of the orientation map, and it is suggested that both maps share a common functional architecture. A17 and A18 are bound at the transition zone by another SF gradient involving the geniculo-cortical and the callosal pathways. A model based on principal component analysis shows that SF maps integrate three different SF-dependent channels. Two of these reflect the segregated excitatory input from X and Y geniculate cells to A17 and A18. The third one conveys a specific combination of excitatory and suppressive inputs to the visual cortex. In a manner coherent with anatomical and electrophysiological data, it is interpreted as originating from a subtype of Y geniculate cells.

Quantitative investigation of the effect of the extra-cerebral vasculature in diffuse optical imaging: a simulation study

Diffuse optical imaging (DOI) is a non invasive technique allowing the recovery of hemodynamic changes in the brain. Due to the diffusive nature of photon propagation in turbid media and the fact that cerebral tissues are located around 1.5 cm under the adult human scalp, DOI measurements are subject to partial volume errors. DOI measurements are also sensitive to large pial vessels because oxygenated and deoxygenated hemoglobin are the dominant chromophores in the near infrared window. In this study, the effect of the extra-cerebral vasculature in proximity of the sagittal sinus was investigated for its impact on DOI measurements simulated over the human adult visual cortex. Numerical Monte Carlo simulations were performed on two specific models of the human head derived from magnetic resonance imaging (MRI) scans. The first model included the extra-cerebral vasculature in which constant hemoglobin concentrations were assumed while the second did not. The screening effect of the vasculature was quantified by comparing recovered hemoglobin changes from each model for different optical arrays and regions of activation. A correction factor accounting for the difference between the recovered and the simulated hemoglobin changes was computed in each case. The results show that changes in hemoglobin concentration are better estimated when the extra-cerebral vasculature is modeled and the correction factors obtained in this case were at least 1.4-fold lower. The effect of the vasculature was also examined in a high-density diffuse optical tomography configuration. In this case, the difference between changes in hemoglobin concentration recovered with each model was reduced down to 10%.

Evaluation of receptive field size from higher harmonics in visuotopic mapping using continuous stimulation optical imaging

The extraction and detection of specific responses from a large amount of background noise has been the subject of a considerable body of research in brain functional imaging, and more specifically in optical intrinsic signal imaging. Recent work by Kalatsky and Stryker (2003) showed that by combining different conditions and using periodic stimuli, recording times can be reduced. Spectral decomposition is then used to provide amplitude and phase information locked to the stimulus. A drawback of the above method is that by focusing only on a single harmonic, response information is limited. The shape of the hemodynamic response function (HRF) and the temporal variations in the neural responses cannot be assessed. In this work it is argued that additional information about neural responses can be gathered by using higher harmonics. Moving bars were used to generate visuotopic maps on large portions of the cat visual cortex. Up to four simultaneously bars moving repetitively across the visual field at different frequencies were used to sample the HRF in the Fourier domain. The HRF profile obtained with continuous stimulation was spatially homogeneous throughout the cortex and similar to the HRF profile obtained using episodic stimulation. Furthermore, by modeling the optical response as a convolution between HRF and neuronal responses, the ratio of the second harmonic to the first provided an estimation of the receptive field size. This was further validated by measuring spatial frequency selectivity. Therefore, the use of higher harmonics opens new avenues to estimate receptive field size from temporal signals.