Iliski, a software for robust calculation of transfer functions

Deconvolution of the Functional Ultrasound Response in the Mouse Visual Pathway Using Block-Term Decomposition

Functional ultrasound (fUS) indirectly measures brain activity by detecting changes in cerebral blood volume following neural activation. Conventional approaches model such functional neuroimaging data as the convolution between an impulse response, known as the hemodynamic response function (HRF), and a binarized representation of the input signal based on the stimulus onsets, the so-called experimental paradigm (EP). However, the EP may not characterize the whole complexity of the activity-inducing signals that evoke the hemodynamic changes. Furthermore, the HRF is known to vary across brain areas and stimuli. To achieve an adaptable framework that can capture such dynamics of the brain function, we model the multivariate fUS time-series as convolutive mixtures and apply block-term decomposition on a set of lagged fUS autocorrelation matrices, revealing both the region-specific HRFs and the source signals that induce the hemodynamic responses. We test our approach on two mouse-based fUS experiments. In the first experiment, we present a single type of visual stimulus to the mouse, and deconvolve the fUS signal measured within the mouse brain's lateral geniculate nucleus, superior colliculus and visual cortex. We show that the proposed method is able to recover back the time instants at which the stimulus was displayed, and we validate the estimated region-specific HRFs based on prior studies. In the second experiment, we alter the location of the visual stimulus displayed to the mouse, and aim at differentiating the various stimulus locations over time by identifying them as separate sources.

Diversity of neurovascular coupling dynamics along vascular arbors in layer II/III somatosensory cortex

The spatial-temporal sequence of cerebral blood flow (CBF), cerebral blood volume (CBV) and blood velocity changes triggered by neuronal activation is critical for understanding functional brain imaging. This sequence follows a stereotypic pattern of changes across different zones of the vasculature in the olfactory bulb, the first relay of olfaction. However, in the cerebral cortex, where most human brain mapping studies are performed, the timing of activity evoked vascular events remains controversial. Here we utilized a single whisker stimulation model to map out functional hyperemia along vascular arbours from layer II/III to the surface of primary somatosensory cortex, in anesthetized and awake Thy1-GCaMP6 mice. We demonstrate that sensory stimulation triggers an increase in blood velocity within the mid-capillary bed and a dilation of upstream large capillaries, and the penetrating and pial arterioles. We report that under physiological stimulation, response onset times are highly variable across compartments of different vascular arbours. Furthermore, generating transfer functions (TFs) between neuronal Ca2+ and vascular dynamics across different brain states demonstrates that anesthesia decelerates neurovascular coupling (NVC). This spatial-temporal pattern of vascular events demonstrates functional diversity not only between different brain regions but also at the level of different vascular arbours within supragranular layers of the cerebral cortex.