Linear systems analysis of the fMRI signal

A Randomization-Based, Model-Free Approach to Functional Neuroimaging: A Proof of Concept

Functional neuroimaging analysis takes noisy multidimensional measurements as input and produces statistical inferences regarding the functional properties of brain regions as output. Such inferences are most commonly model-based, in that they assume a model of how neural activity translates to the measured signal (blood oxygenation level-dependent signal in the case of functional MRI). The use of models increases statistical sensitivity and makes it possible to ask fine-grained theoretical questions. However, this comes at the cost of making theoretical assumptions about the underlying data-generating process. An advantage of model-free approaches is that they can be used in cases where model assumptions are known not to hold. To this end, we introduce a randomization-based, model-free approach to functional neuroimaging. TWISTER randomization makes it possible to infer functional selectivity from correlations between experimental runs. We provide a proof of concept in the form of a visuomotor mapping experiment and discuss the possible strengths and limitations of this new approach in light of our empirical results.

Toward a more informative representation of the fetal-neonatal brain connectome using variational autoencoder

Recent advances in functional magnetic resonance imaging (fMRI) have helped elucidate previously inaccessible trajectories of early-life prenatal and neonatal brain development. To date, the interpretation of fetal-neonatal fMRI data has relied on linear analytic models, akin to adult neuroimaging data. However, unlike the adult brain, the fetal and newborn brain develops extraordinarily rapidly, far outpacing any other brain development period across the life span. Consequently, conventional linear computational models may not adequately capture these accelerated and complex neurodevelopmental trajectories during this critical period of brain development along the prenatal-neonatal continuum. To obtain a nuanced understanding of fetal-neonatal brain development, including nonlinear growth, for the first time, we developed quantitative, systems-wide representations of brain activity in a large sample (>500) of fetuses, preterm, and full-term neonates using an unsupervised deep generative model called variational autoencoder (VAE), a model previously shown to be superior to linear models in representing complex resting-state data in healthy adults. Here, we demonstrated that nonlinear brain features, that is, latent variables, derived with the VAE pretrained on rsfMRI of human adults, carried important individual neural signatures, leading to improved representation of prenatal-neonatal brain maturational patterns and more accurate and stable age prediction in the neonate cohort compared to linear models. Using the VAE decoder, we also revealed distinct functional brain networks spanning the sensory and default mode networks. Using the VAE, we are able to reliably capture and quantify complex, nonlinear fetal-neonatal functional neural connectivity. This will lay the critical foundation for detailed mapping of healthy and aberrant functional brain signatures that have their origins in fetal life.

Reconstructing sources location of visual color cortex by the task-irrelevant visual stimuli through machine learning decoding

Visual color sensing is generated by electrical discharges from endocranial neuronal sources that penetrate the skull and reach to the cerebral cortex. However, the space location of the source generated by this neural mechanism remains elusive. In this paper, we emulate the generation of visual color signal by task-irrelevant stimuli to activate brain neurons, where its consequences over the cerebral cortex is experimentally tracked. We first document the changes to brain color sensing using electroencephalography (EEG), and find that the sensing classification accuracy of primary visual cortex (V1) regions was positively correlated with the space correlation of visual evoked potential (VEP) power distribution under machine learning decoding. We then explore the decoded results to trace the brain activity neural source location of EEG inversion problem and assess its reconstructive possibility. We show that visual color EEG in V1 can reconstruct endocranial neuronal source location, through the machine learning decoding of channel location.

What has vision science taught us about functional MRI?

In the domain of human neuroimaging, much attention has been paid to the question of whether and how the development of functional magnetic resonance imaging (fMRI) has advanced our scientific knowledge of the human brain. However, the opposite question is also important; how has our knowledge of the brain advanced our understanding of fMRI? Here, we discuss how and why scientific knowledge about the human and animal visual system has been used to answer fundamental questions about fMRI as a brain measurement tool and how these answers have contributed to scientific discoveries beyond vision science.

Detection of functional activity in brain white matter using fiber architecture informed synchrony mapping

A general linear model is widely used for analyzing fMRI data, in which the blood oxygenation-level dependent (BOLD) signals in gray matter (GM) evoked in response to neural stimulation are modeled by convolving the time course of the expected neural activity with a canonical hemodynamic response function (HRF) obtained a priori. The maps of brain activity produced reflect the magnitude of local BOLD responses. However, detecting BOLD signals in white matter (WM) is more challenging as the BOLD signals are weaker and the HRF is different, and may vary more across the brain. Here we propose a model-free approach to detect changes in BOLD signals in WM by measuring task-evoked increases of BOLD signal synchrony in WM fibers. The proposed approach relies on a simple assumption that, in response to a functional task, BOLD signals in relevant fibers are modulated by stimulus-evoked neural activity and thereby show greater synchrony than when measured in a resting state, even if their magnitudes do not change substantially. This approach is implemented in two technical stages. First, for each voxel a fiber-architecture-informed spatial window is created with orientation distribution functions constructed from diffusion imaging data. This provides the basis for defining neighborhoods in WM that share similar local fiber architectures. Second, a modified principal component analysis (PCA) is used to estimate the synchrony of BOLD signals in each spatial window. The proposed approach is validated using a 3T fMRI dataset from the Human Connectome Project (HCP) at a group level. The results demonstrate that neural activity can be reliably detected as increases in fMRI signal synchrony within WM fibers that are engaged in a task with high sensitivities and reproducibility.

Temporal stability of the hemodynamic response function across the majority of human cerebral cortex

The hemodynamic response function (HRF) measured with functional magnetic resonance imaging is generated by vascular and metabolic responses evoked by brief (<4 s) stimuli. It is known that the human HRF varies across cortex, between subjects, with stimulus paradigms, and even between different measurements in the same cortical location. However, our results demonstrate that strong HRFs are remarkably repeatable across sessions separated by time intervals up to 3 months. In this study, a multisensory stimulus was used to activate and measure the HRF across the majority of cortex (>70%, with lesser reliability observed in some areas of prefrontal cortex). HRFs were measured with high spatial resolution (2‐mm voxels) in central gray matter to minimize variations caused by partial‐volume effects. HRF amplitudes and temporal dynamics were highly repeatable across four sessions in 20 subjects. Positive and negative HRFs were consistently observed across sessions and subjects. Negative HRFs were generally weaker and, thus, more variable than positive HRFs. Statistical measurements showed that across‐session variability is highly correlated to the variability across events within a session; these measurements also indicated a normal distribution of variability across cortex. The overall repeatability of the HRFs over long time scales generally supports the long‐term use of event‐related functional magnetic resonance imaging protocols.

Lag-Optimized Blood Oxygenation Level Dependent Cerebrovascular Reactivity Estimates Derived From Breathing Task Data Have a Stronger Relationship With Baseline Cerebral Blood Flow

Cerebrovascular reactivity (CVR), an important indicator of cerebrovascular health, is commonly studied with the Blood Oxygenation Level Dependent functional MRI (BOLD-fMRI) response to a vasoactive stimulus. Theoretical and empirical evidence suggests that baseline cerebral blood flow (CBF) modulates BOLD signal amplitude and may influence BOLD-CVR estimates. We address how acquisition and modeling choices affect the relationship between baseline cerebral blood flow (bCBF) and BOLD-CVR: whether BOLD-CVR is modeled with the inclusion of a breathing task, and whether BOLD-CVR amplitudes are optimized for hemodynamic lag effects. We assessed between-subject correlations of average GM values and within-subject spatial correlations across cortical regions. Our results suggest that a breathing task addition to a resting-state acquisition, alongside lag-optimization within BOLD-CVR modeling, can improve BOLD-CVR correlations with bCBF, both between- and within-subjects, likely because these CVR estimates are more physiologically accurate. We report positive correlations between bCBF and BOLD-CVR, both between- and within-subjects. The physiological explanation of this positive correlation is unclear; research with larger samples and tightly controlled vasoactive stimuli is needed. Insights into what drives variability in BOLD-CVR measurements and related measurements of cerebrovascular function are particularly relevant when interpreting results in populations with altered vascular and/or metabolic baselines or impaired cerebrovascular reserve.

Temporal dynamics of the neural representation of hue and luminance polarity

Hue and luminance contrast are basic visual features. Here we use multivariate analyses of magnetoencephalography data to investigate the timing of the neural computations that extract them, and whether they depend on common neural circuits. We show that hue and luminance-contrast polarity can be decoded from MEG data and, with lower accuracy, both features can be decoded across changes in the other feature. These results are consistent with the existence of both common and separable neural mechanisms. The decoding time course is earlier and more temporally precise for luminance polarity than hue, a result that does not depend on task, suggesting that luminance contrast is an updating signal that separates visual events. Meanwhile, cross-temporal generalization is slightly greater for representations of hue compared to luminance polarity, providing a neural correlate of the preeminence of hue in perceptual grouping and memory. Finally, decoding of luminance polarity varies depending on the hues used to obtain training and testing data. The pattern of results is consistent with observations that luminance contrast is mediated by both L-M and S cone sub-cortical mechanisms.

Robust priors for regularized regression

Induction benefits from useful priors. Penalized regression approaches, like ridge regression, shrink weights toward zero but zero association is usually not a sensible prior. Inspired by simple and robust decision heuristics humans use, we constructed non-zero priors for penalized regression models that provide robust and interpretable solutions across several tasks. Our approach enables estimates from a constrained model to serve as a prior for a more general model, yielding a principled way to interpolate between models of differing complexity. We successfully applied this approach to a number of decision and classification problems, as well as analyzing simulated brain imaging data. Models with robust priors had excellent worst-case performance. Solutions followed from the form of the heuristic that was used to derive the prior. These new algorithms can serve applications in data analysis and machine learning, as well as help in understanding how people transition from novice to expert performance.

A tutorial on generalized eigendecomposition for denoising, contrast enhancement, and dimension reduction in multichannel electrophysiology

The goal of this paper is to present a theoretical and practical introduction to generalized eigendecomposition (GED), which is a robust and flexible framework used for dimension reduction and source separation in multichannel signal processing. In cognitive electrophysiology, GED is used to create spatial filters that maximize a researcher-specified contrast. For example, one may wish to exploit an assumption that different sources have different frequency content, or that sources vary in magnitude across experimental conditions. GED is fast and easy to compute, performs well in simulated and real data, and is easily adaptable to a variety of specific research goals. This paper introduces GED in a way that ties together myriad individual publications and applications of GED in electrophysiology, and provides sample MATLAB and Python code that can be tested and adapted. Practical considerations and issues that often arise in applications are discussed.

Imaging faster neural dynamics with fast fMRI: A need for updated models of the hemodynamic response

Fast fMRI enables the detection of neural dynamics over timescales of hundreds of milliseconds, suggesting it may provide a new avenue for studying subsecond neural processes in the human brain. The magnitudes of these fast fMRI dynamics are far greater than predicted by canonical models of the hemodynamic response. Several studies have established nonlinear properties of the hemodynamic response that have significant implications for fast fMRI. We first review nonlinear properties of the hemodynamic response function that may underlie fast fMRI signals. We then illustrate the breakdown of canonical hemodynamic response models in the context of fast neural dynamics. We will then argue that the canonical hemodynamic response function is not likely to reflect the BOLD response to neuronal activity driven by sparse or naturalistic stimuli or perhaps to spontaneous neuronal fluctuations in the resting state. These properties suggest that fast fMRI is capable of tracking surprisingly fast neuronal dynamics, and we discuss the neuroscientific questions that could be addressed using this approach.

rsHRF: A toolbox for resting-state HRF estimation and deconvolution

The hemodynamic response function (HRF) greatly influences the intra- and inter-subject variability of brain activation and connectivity, and might confound the estimation of temporal precedence in connectivity analyses, making its estimation necessary for a correct interpretation of neuroimaging studies. Additionally, the HRF shape itself is a useful local measure. However, most algorithms for HRF estimation are specific for task-related fMRI data, and only a few can be directly applied to resting-state protocols. Here we introduce rsHRF, a Matlab and Python toolbox that implements HRF estimation and deconvolution from the resting-state BOLD signal. We first provide an overview of the main algorithm, practical implementations, and then demonstrate the feasibility and usefulness of rsHRF by validation experiments with a publicly available resting-state fMRI dataset. We also provide tools for statistical analyses and visualization. We believe that this toolbox may significantly contribute to a better analysis and understanding of the components and variability of BOLD signals.

Voxel-Wise Linearity Analysis of Increments and Decrements in BOLD Responses in Human Visual Cortex Using a Contrast Adaptation Paradigm

The linearity of BOLD responses is a fundamental presumption in most analysis procedures for BOLD fMRI studies. Previous studies have examined the linearity of BOLD signal increments, but less is known about the linearity of BOLD signal decrements. The present study assessed the linearity of both BOLD signal increments and decrements in the human primary visual cortex using a contrast adaptation paradigm. Results showed that both BOLD signal increments and decrements kept linearity to long stimuli (e.g., 3 s, 6 s), yet, deviated from linearity to transient stimuli (e.g., 1 s). Furthermore, a voxel-wise analysis showed that the deviation patterns were different for BOLD signal increments and decrements: while the BOLD signal increments demonstrated a consistent overestimation pattern, the patterns for BOLD signal decrements varied from overestimation to underestimation. Our results suggested that corrections to deviations from linearity of transient responses should consider the different effects of BOLD signal increments and decrements.

Statistical power or more precise insights into neuro-temporal dynamics? Assessing the benefits of rapid temporal sampling in fMRI

Functional magnetic resonance imaging (fMRI), a non-invasive and widely used human neuroimaging method, is most known for its spatial precision. However, there is a growing interest in its temporal sensitivity. This is despite the temporal blurring of neuronal events by the blood oxygen level dependent (BOLD) signal, the peak of which lags neuronal firing by 4-6 seconds. Given this, the goal of this review is to answer a seemingly simple question - "What are the benefits of increased temporal sampling for fMRI?". To answer this, we have combined fMRI data collected at multiple temporal scales, from 323 to 1000 milliseconds, with a review of both historical and contemporary temporal literature. After a brief discussion of technological developments that have rekindled interest in temporal research, we next consider the potential statistical and methodological benefits. Most importantly, we explore how fast fMRI can uncover previously unobserved neuro-temporal dynamics - effects that are entirely missed when sampling at conventional 1 to 2 second rates. With the intrinsic link between space and time in fMRI, this temporal renaissance also delivers improvements in spatial precision. Far from producing only statistical gains, the array of benefits suggest that the continued temporal work is worth the effort.

Physiological noise modeling in fMRI based on the pulsatile component of photoplethysmograph

The blood oxygenation level-dependent (BOLD) contrast mechanism allows the noninvasive monitoring of changes in deoxyhemoglobin content. As such, it is commonly used in functional magnetic resonance imaging (fMRI) to study brain activity since levels of deoxyhemoglobin are indirectly related to local neuronal activity through neurovascular coupling mechanisms. However, the BOLD signal is severely affected by physiological processes as well as motion. Due to this, several noise correction techniques have been developed to correct for the associated confounds. The present study focuses on cardiac pulsatility fMRI confounds, aiming to refine model-based techniques that utilize the photoplethysmograph (PPG) signal. Specifically, we propose a new technique based on convolution filtering, termed cardiac pulsatility model (CPM) and compare its performance with the cardiac-related RETROICOR (Card-RETROICOR), which is a technique commonly used to model fMRI fluctuations due to cardiac pulsatility. Further, we investigate whether variations in the amplitude of the PPG pulses (PPG-Amp) covary with variations in amplitude of pulse-related fMRI fluctuations, as well as with the systemic low frequency oscillations (SLFOs) component of the fMRI global signal (GS - defined as the mean signal across all gray matter voxels). Capitalizing on 3T fMRI data from the Human Connectome Project, CPM was found to explain a significantly larger fraction of the fMRI signal variance compared to Card-RETROICOR, particularly for subjects with larger heart rate variability during the scan. The amplitude of the fMRI pulse-related fluctuations did not covary with PPG-Amp; however, PPG-Amp explained significant variance in the GS that was not attributed to variations in heart rate or breathing patterns. Our results suggest that the proposed approach can model high-frequency fluctuations due to pulsation as well as low-frequency physiological fluctuations more accurately compared to model-based techniques commonly employed in fMRI studies.

Origins of 1/f-like tissue oxygenation fluctuations in the murine cortex

The concentration of oxygen in the brain spontaneously fluctuates, and the distribution of power in these fluctuations has a 1/f-like spectra, where the power present at low frequencies of the power spectrum is orders of magnitude higher than at higher frequencies. Though these oscillations have been interpreted as being driven by neural activity, the origin of these 1/f-like oscillations is not well understood. Here, to gain insight of the origin of the 1/f-like oxygen fluctuations, we investigated the dynamics of tissue oxygenation and neural activity in awake behaving mice. We found that oxygen signal recorded from the cortex of mice had 1/f-like spectra. However, band-limited power in the local field potential did not show corresponding 1/f-like fluctuations. When local neural activity was suppressed, the 1/f-like fluctuations in oxygen concentration persisted. Two-photon measurements of erythrocyte spacing fluctuations and mathematical modeling show that stochastic fluctuations in erythrocyte flow could underlie 1/f-like dynamics in oxygenation. These results suggest that the discrete nature of erythrocytes and their irregular flow, rather than fluctuations in neural activity, could drive 1/f-like fluctuations in tissue oxygenation.

Iliski, a software for robust calculation of transfer functions

Understanding the relationships between biological processes is paramount to unravel pathophysiological mechanisms. These relationships can be modeled with Transfer Functions (TFs), with no need of a priori hypotheses as to the shape of the transfer function. Here we present Iliski, a software dedicated to TFs computation between two signals. It includes different pre-treatment routines and TF computation processes: deconvolution, deterministic and non-deterministic optimization algorithms that are adapted to disparate datasets. We apply Iliski to data on neurovascular coupling, an ensemble of cellular mechanisms that link neuronal activity to local changes of blood flow, highlighting the software benefits and caveats in the computation and evaluation of TFs. We also propose a workflow that will help users to choose the best computation according to the dataset. Iliski is available under the open-source license CC BY 4.0 on GitHub (https://github.com/alike-aydin/Iliski) and can be used on the most common operating systems, either within the MATLAB environment, or as a standalone application.

Iliski is a software helping the user to find the relationship between two sets of data, namely transfer functions. Although transfer functions are widely used in many scientific fields to link two signals, their computation can be tricky due to data features such as multisource noise, or to specific shape requirements imposed by the nature of the signals, e.g. in biological data. Iliski offers a user-friendly graphical interface to ease the computation of transfer functions for both experienced and users with no coding skills. It proposes several signal pre-processing methods and allows rapid testing of different computing approaches, either based on deconvolution or on optimization of multi-parametric functions. This article, combined with a User Manual, provides a detailed description of Iliski functionalities and a thorough description of the advantages and drawbacks of each computing method using experimental biological data. In the era of Big Data, scientists strive to find new models for patho-physiological mechanisms, and Iliski fulfils the requirements of rigorous, flexible, and fast data driven hypothesis testing.

The neural-vascular basis of age-related processing speed decline

Most studies examining neurocognitive aging are based on the blood-oxygen level-dependent signal obtained during functional magnetic resonance imaging (fMRI). The physiological basis of this signal is neural-vascular coupling, the process by which neurons signal cerebrovasculature to dilate in response to an increase in active neural metabolism due to stimulation. These fMRI studies of aging rely on the hemodynamic equivalence assumption that this process is not disrupted by physiologic deterioration associated with aging. Studies of neural-vascular coupling challenge this assumption and show that neural-vascular coupling is closely related to cognition. In this review, we put forward a theory of processing speed decline in aging and how it is related to age-related neural-vascular coupling changes based on the results of studies elucidating the relationships between cognition, cerebrovascular dynamics, and aging.

Hemodynamic matrix factorization for functional magnetic resonance imaging

The General Linear Model (GLM) used in task-fMRI relates activated brain areas to extrinsic task conditions. The translation of resulting neural activation into a hemodynamic response is commonly approximated with a linear convolution model using a hemodynamic response function (HRF). There are two major limitations in GLM analysis. Firstly, the GLM assumes that neural activation is either on or off and matches the exact stimulus duration in the corresponding task timings. Secondly, brain networks observed in resting-state fMRI experiments present also during task experiments, but the GLM approach models these task-unrelated brain activity as noise. A novel kernel matrix factorization approach, called hemodynamic matrix factorization (HMF), is therefore proposed that addresses both limitations by assuming that task-related and task-unrelated brain activity can be modeled with the same convolution model as in GLM analysis. By contrast to the GLM, the proposed HMF is a blind source separation (BSS) technique, which decomposes fMRI data into modes. Each mode comprises of a neural activation time course and a spatial mapping. Two versions of HMF are proposed in which the neural activation time course of each mode is convolved with either the canonical HRF or predetermined subject-specific HRFs. Firstly, HMF with the canonical HRF is applied to two open-source cohorts. These cohorts comprise of several task experiments including motor, incidental memory, spatial coherence discrimination, verbal discrimination task and a very short localization task, engaging multiple parts of the eloquent cortex. HMF modes were obtained whose neural activation time course followed original task timings and whose corresponding spatial map matched cortical areas known to be involved in the respective task processing. Secondly, the alignment of these neural activation time courses to task timings were further improved by replacing the canonical HRF with subject-specific HRFs during HMF mode computation. In addition to task-related modes, HMF also produced seemingly task-unrelated modes whose spatial maps matched known resting-state networks. The validity of a fMRI task experiment relies on the assumption that the exposure to a stimulus for a given time causes an imminent increase in neural activation of equal duration. The proposed HMF is an attempt to falsify this assumption and allows to identify subject task participation that does not comply with the experiment instructions.

Validating Linear Systems Analysis for Laminar fMRI: Temporal Additivity for Stimulus Duration Manipulations

Advancements in ultra-high field (7 T and higher) magnetic resonance imaging (MRI) scanners have made it possible to investigate both the structure and function of the human brain at a sub-millimeter scale. As neuronal feedforward and feedback information arrives in different layers, sub-millimeter functional MRI has the potential to uncover information processing between cortical micro-circuits across cortical depth, i.e. laminar fMRI. For nearly all conventional fMRI analyses, the main assumption is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal adheres to the principles of linear systems theory. For laminar fMRI, however, directional blood pooling across cortical depth stemming from the anatomy of the cortical vasculature, potentially violates these linear system assumptions, thereby complicating analysis and interpretation. Here we assess whether the temporal additivity requirement of linear systems theory holds for laminar fMRI. We measured responses elicited by viewing stimuli presented for different durations and evaluated how well the responses to shorter durations predicted those elicited by longer durations. We find that BOLD response predictions are consistently good predictors for observed responses, across all cortical depths, and in all measured visual field maps (V1, V2, and V3). Our results suggest that the temporal additivity assumption for linear systems theory holds for laminar fMRI. We thus show that the temporal additivity assumption holds across cortical depth for sub-millimeter gradient-echo BOLD fMRI in early visual cortex.

Comparing hand movement rate dependence of cerebral blood volume and BOLD responses at 7T

Functional magnetic resonance imaging (fMRI) based on the Blood Oxygenation Level Dependent (BOLD) contrast takes advantage of the coupling between neuronal activity and the hemodynamics to allow a non-invasive localisation of the neuronal activity. In general, fMRI experiments assume a linear relationship between neuronal activation and the observed hemodynamics. However, the relationship between BOLD responses, neuronal activity, and behaviour are often nonlinear. In addition, the nonlinearity between BOLD responses and behaviour may be related to neuronal process rather than a neurovascular uncoupling. Further, part of the nonlinearity may be driven by vascular nonlinearity effects in particular from large vessel contributions. fMRI based on cerebral blood volume (CBV), promises a higher microvascular specificity, potentially without vascular nonlinearity effects and reduced contamination of the large draining vessels compared to BOLD. In this study, we aimed to investigate differences in BOLD and VASO-CBV signal changes during a hand movement task over a broad range of movement rates. We used a double readout 3D-EPI sequence at 7T to simultaneously measure VASO-CBV and BOLD responses in the sensorimotor cortex. The measured BOLD and VASO-CBV responses increased very similarly in a nonlinear fashion, plateauing for movement rates larger than 1 Hz. Our findings show a tight relationship between BOLD and VASO-CBV responses, indicating that the overall interplay of CBV and BOLD responses are similar for the assessed range of movement rates. These results suggest that the observed nonlinearity of neuronal origin is already present in VASO-CBV measurements, and consequently shows relatively unchanged BOLD responses.

Physiological Gaussian process priors for the hemodynamics in fMRI analysis

BACKGROUND

Inference from fMRI data faces the challenge that the hemodynamic system that relates neural activity to the observed BOLD fMRI signal is unknown.

NEW METHOD

We propose a new Bayesian model for task fMRI data with the following features: (i) joint estimation of brain activity and the underlying hemodynamics, (ii) the hemodynamics is modeled nonparametrically with a Gaussian process (GP) prior guided by physiological information and (iii) the predicted BOLD is not necessarily generated by a linear time-invariant (LTI) system. We place a GP prior directly on the predicted BOLD response, rather than on the hemodynamic response function as in previous literature. This allows us to incorporate physiological information via the GP prior mean in a flexible way, and simultaneously gives us the nonparametric flexibility of the GP.

RESULTS

Results on simulated data show that the proposed model is able to discriminate between active and non-active voxels also when the GP prior deviates from the true hemodynamics. Our model finds time varying dynamics when applied to real fMRI data.

COMPARISON WITH EXISTING METHOD(S)

The proposed model is better at detecting activity in simulated data than standard models, without inflating the false positive rate. When applied to real fMRI data, our GP model in several cases finds brain activity where previously proposed LTI models does not.

CONCLUSIONS

We have proposed a new non-linear model for the hemodynamics in task fMRI, that is able to detect active voxels, and gives the opportunity to ask new kinds of questions related to hemodynamics.

So young, yet so mature? Electrophysiological and vascular correlates of phonotactic processing in 18-month-olds

The present study investigated neural correlates of implicit phonotactic processing in 18-month-old children that just reached an important step in language development: the vocabulary spurt. Pseudowords, either phonotactically legal or illegal with respect to their native language, were acoustically presented to monolingually German raised infants. Neural activity was simultaneously assessed by means of electroencephalography (EEG) and functional near-infrared spectroscopy (fNIRS). The former method excellently tracks fast processing mechanisms, whereas the latter reveals brain areas recruited. Results of the present study indicate that 18-month-olds recognize the linguistic properties of their native language based on phonotactics. This manifested in an increased N400 for legal compared to illegal pseudowords in the EEG conforming to adult-like mechanisms. Unfortunately, fNIRS findings did not support this discrimination ability. Possible methodological and brain maturational reasons might explain this null finding. This study provides evidence for the advantage of a multi-methodological approach in order to get a clear picture on neural language development.

Uncovering electrophysiological and vascular signatures of implicit emotional prosody

The capability of differentiating between various emotional states in speech displays a crucial prerequisite for successful social interactions. The aim of the present study was to investigate neural processes underlying this differentiating ability by applying a simultaneous neuroscientific approach in order to gain both electrophysiological (via electroencephalography, EEG) and vascular (via functional near-infrared-spectroscopy, fNIRS) responses. Pseudowords conforming to angry, happy, and neutral prosody were presented acoustically to participants using a passive listening paradigm in order to capture implicit mechanisms of emotional prosody processing. Event-related brain potentials (ERPs) revealed a larger P200 and an increased late positive potential (LPP) for happy prosody as well as larger negativities for angry and neutral prosody compared to happy prosody around 500 ms. FNIRS results showed increased activations for angry prosody at right fronto-temporal areas. Correlation between negativity in the EEG and activation in fNIRS for angry prosody suggests analogous underlying processes resembling a negativity bias. Overall, results indicate that mechanisms of emotional and phonological encoding (P200), emotional evaluation (increased negativities) as well as emotional arousal and relevance (LPP) are present during implicit processing of emotional prosody.

Linear systems analysis for laminar fMRI: Evaluating BOLD amplitude scaling for luminance contrast manipulations

A fundamental assumption of nearly all functional magnetic resonance imaging (fMRI) analyses is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal can be described as following linear systems theory. With the advent of ultra-high field (7T and higher) MRI scanners, it has become possible to perform sub-millimeter resolution fMRI in humans. A novel and promising application of sub-millimeter fMRI is measuring responses across cortical depth, i.e. laminar imaging. However, the cortical vasculature and associated directional blood pooling towards the pial surface strongly influence the cortical depth-dependent BOLD signal, particularly for gradient-echo BOLD. This directional pooling may potentially affect BOLD linearity across cortical depth. Here we assess whether the amplitude scaling assumption for linear systems theory holds across cortical depth. For this, we use stimuli with different luminance contrasts to elicit different BOLD response amplitudes. We find that BOLD amplitude across cortical depth scales with luminance contrast, and that this scaling is identical across cortical depth. Although nonlinearities may be present for different stimulus configurations and acquisition protocols, our results suggest that the amplitude scaling assumption for linear systems theory across cortical depth holds for luminance contrast manipulations in sub-millimeter laminar BOLD fMRI.

How the Brain Understands Spoken and Sung Sentences

The present study investigates whether meaning is similarly extracted from spoken and sung sentences. For this purpose, subjects listened to semantically correct and incorrect sentences while performing a correctness judgement task. In order to examine underlying neural mechanisms, a multi-methodological approach was chosen combining two neuroscientific methods with behavioral data. In particular, fast dynamic changes reflected in the semantically associated N400 component of the electroencephalography (EEG) were simultaneously assessed with the topographically more fine-grained vascular signals acquired by the functional near-infrared spectroscopy (fNIRS). EEG results revealed a larger N400 for incorrect compared to correct sentences in both spoken and sung sentences. However, the N400 was delayed for sung sentences, potentially due to the longer sentence duration. fNIRS results revealed larger activations for spoken compared to sung sentences irrespective of semantic correctness at predominantly left-hemispheric areas, potentially suggesting a greater familiarity with spoken material. Furthermore, the fNIRS revealed a widespread activation for correct compared to incorrect sentences irrespective of modality, potentially indicating a successful processing of sentence meaning. The combined results indicate similar semantic processing in speech and song.

Vascular and neural basis of the BOLD signal

Neural activity in the brain is usually coupled to increases in local cerebral blood flow, leading to the increase in oxygenation that generates the BOLD fMRI signal. Recent work has begun to elucidate the vascular and neural mechanisms underlying the BOLD signal. The dilatory response is distributed throughout the vascular network. Arteries actively dilate within a second following neural activity increases, while venous distensions are passive and have a time course that last tens of seconds. Vasodilation, and thus local blood flow, is controlled by the activity of both neurons and astrocytes via multiple different pathways. The relationship between sensory-driven neural activity and the vascular dynamics in sensory areas are well-captured with a linear convolution model. However, depending on the behavioral state or brain region, the coupling between neural activity and hemodynamic signals can be weak or even inverted.

Mesoscopic and microscopic imaging of sensory responses in the same animal

Imaging based on blood flow dynamics is widely used to study sensory processing. Here we investigated the extent to which local neuronal and capillary responses (two-photon microscopy) are correlated to mesoscopic responses detected with fast ultrasound (fUS) and BOLD-fMRI. Using a specialized chronic olfactory bulb preparation, we report that sequential imaging of the same mouse allows quantitative comparison of odour responses, imaged at both microscopic and mesoscopic scales. Under these conditions, functional hyperaemia occurred at the threshold of neuronal activation and fUS-CBV signals could be detected at the level of single voxels with activation maps varying according to blood velocity. Both neuronal and vascular responses increase non-linearly as a function of odour concentration, whereas both microscopic and mesoscopic vascular responses are linearly correlated to local neuronal calcium. These data establish strengths and limits of mesoscopic imaging techniques to report neural activity.

Neuronal activity leads to a local increase in blood flow and volume, a process termed hyperaemia. Here, the authors employ multiple imaging approaches of neuronal and vascular activity at varying resolution to delineate the spatiotemporal dynamics of neurovascular coupling evoked by odours in the olfactory bulb.

Bilingual and Monolingual First Language Acquisition Experience Differentially Shapes Children's Property Term Learning: Evidence from Behavioral and Neurophysiological Measures

Studies of novel noun learning show bilingual children rely less on the Mutual Exclusivity Constraint (MEC) for word learning than monolinguals. Shifting the focus to learning novel property terms (adjectives), the present study compared 3.5- and five-year-old bilingual and monolingual preschoolers' adherence to the MEC. We found no bilingual-monolingual differences on a behavioral forced-choice task for the 3.5-year-olds, but five-year-old monolinguals adhered more to the MEC than bilinguals did. Older bilinguals adhered less to the MEC than younger ones, while there was no difference in MEC adherence between the younger and older monolinguals. In the 5-year-olds, we additionally acquired neurophysiological data using functional near-infrared spectroscopy (fNIRS) to allow for a first explorative look at potential neuronal underpinnings. The data show that, compared to bilinguals, monolinguals reveal higher activation over three brain regions (right frontal, left temporo-parietal, and left prefrontal) that may be involved in exploiting the MEC, building on conflict detection, inhibition, solution of a disjunction, and working memory processes. Taken together, our behavioral and neurophysiological findings reveal different paths towards novel property term learning depending on children's language acquisition context.

Optimising neonatal fMRI data analysis: Design and validation of an extended dHCP preprocessing pipeline to characterise noxious-evoked brain activity in infants

The infant brain is unlike the adult brain, with considerable differences in morphological, neurodynamic, and haemodynamic features. As the majority of current MRI analysis tools were designed for use in adults, a primary objective of the Developing Human Connectome Project (dHCP) is to develop optimised methodological pipelines for the analysis of neonatal structural, resting state, and diffusion MRI data. Here, in an independent neonatal dataset we have extended and optimised the dHCP fMRI preprocessing pipeline for the analysis of stimulus-response fMRI data. We describe and validate this extended dHCP fMRI preprocessing pipeline to analyse changes in brain activity evoked following an acute noxious stimulus applied to the infant's foot. We compare the results obtained from this extended dHCP pipeline to results obtained from a typical FSL FEAT-based analysis pipeline, evaluating the pipelines' outputs using a wide range of tests. We demonstrate that a substantial increase in spatial specificity and sensitivity to signal can be attained with a bespoke neonatal preprocessing pipeline through optimised motion and distortion correction, ICA-based denoising, and haemodynamic modelling. The improved sensitivity and specificity, made possible with this extended dHCP pipeline, will be paramount in making further progress in our understanding of the development of sensory processing in the infant brain.

Unraveling the spatiotemporal brain dynamics during a simulated reach-to-eat task

The reach-to-eat task involves a sequence of action components including looking, reaching, grasping, and feeding. While cortical representations of individual action components have been mapped in human functional magnetic resonance imaging (fMRI) studies, little is known about the continuous spatiotemporal dynamics among these representations during the reach-to-eat task. In a periodic event-related fMRI experiment, subjects were scanned while they reached toward a food image, grasped the virtual food, and brought it to their mouth within each 16-s cycle. Fourier-based analysis of fMRI time series revealed periodic signals and noise distributed across the brain. Independent component analysis was used to remove periodic or aperiodic motion artifacts. Time-frequency analysis was used to analyze the temporal characteristics of periodic signals in each voxel. Circular statistics was then used to estimate mean phase angles of periodic signals and select voxels based on the distribution of phase angles. By sorting mean phase angles across regions, we were able to show the real-time spatiotemporal brain dynamics as continuous traveling waves over the cortical surface. The activation sequence consisted of approximately the following stages: (1) stimulus related activations in occipital and temporal cortices; (2) movement planning related activations in dorsal premotor and superior parietal cortices; (3) reaching related activations in primary sensorimotor cortex and supplementary motor area; (4) grasping related activations in postcentral gyrus and sulcus; (5) feeding related activations in orofacial areas. These results suggest that phase-encoded design and analysis can be used to unravel sequential activations among brain regions during a simulated reach-to-eat task.

Task-Evoked Negative BOLD Response in the Default Mode Network Does Not Alter Its Functional Connectivity

While functional connectivity networks are often extracted from resting-state fMRI scans, they have been shown to be active during task performance as well. However, the effect of an in-scanner task on functional connectivity networks is not completely understood. While there is evidence that task-evoked positive BOLD response can alter functional connectivity networks, particularly in the primary sensorimotor cortices, the effect of task-evoked negative BOLD response on the functional connectivity of the Default mode network (DMN) is somewhat ambiguous. In this study, we aim to investigate whether task performance, which is associated with negative BOLD response in the DMN regions, alters the time-course of functional connectivity in the same regions obtained by independent component analysis (ICA). ICA has been used to effectively extract functional connectivity networks during task performance and resting-state. We first demonstrate that performing a simple visual-motor task alters the temporal time-course of the network extracted from the primary visual cortex. Then we show that despite detecting a robust task-evoked negative BOLD response in the DMN regions, a simple visual-motor task does not alter the functional connectivity of the DMN regions. Our findings suggest that different mechanisms may underlie the relationship between task-related activation/deactivation networks and the overlapping functional connectivity networks in the human large-scale brain networks.

Large-scale dynamic modeling of task-fMRI signals via subspace system identification

OBJECTIVE

We analyze task-based fMRI time series to produce large-scale dynamical models that are capable of approximating the observed signal with good accuracy.

APPROACH

We extend subspace system identification methods for deterministic and stochastic state-space models with external inputs. The dynamic behavior of the generated models is characterized using control-theoretic analysis tools. To validate their effectiveness, we perform a probabilistic inversion of the identified input-output relationships via joint state-input maximum likelihood estimation. Our experimental setup explores a large dataset generated using state-of-the-art acquisition and pre-processing methods from the Human Connectome Project.

MAIN RESULTS

We analyze both anatomically parcellated and spatially dense time series, and propose an efficient algorithm to address the high-dimensional optimization problem resulting from the latter. Our results enable the quantification of input-output transfer functions between each task condition and each region of the cortex, as exemplified by a motor task. Further, the identified models produce impulse response functions between task conditions and cortical regions that are compatible with typical hemodynamic response functions. We then extend subspace methods to account for multi-subject experimental configurations, identifying models that capture common dynamical characteristics across subjects. Finally, we show that system inversion via maximum-likelihood allows the time-of-occurrence of the task stimuli to be estimated from the observed outputs.

SIGNIFICANCE

The ability to produce dynamical input-output models might have an impact in the expanding field of neurofeedback. In particular, the models we produce allow the partial quantification of the effect of external task-related inputs on the metabolic response of the brain, conditioned on its current state. Such a notion provides a basis for leveraging control-theoretic approaches to neuromodulation and self-regulation in therapeutic applications.

A quantitative framework for motion visibility in human cortex

Despite the central use of motion visibility to reveal the neural basis of perception, perceptual decision making, and sensory inference there exists no comprehensive quantitative framework establishing how motion visibility parameters modulate human cortical response. Random-dot motion stimuli can be made less visible by reducing image contrast or motion coherence, or by shortening the stimulus duration. Because each of these manipulations modulates the strength of sensory neural responses they have all been extensively used to reveal cognitive and other nonsensory phenomena such as the influence of priors, attention, and choice-history biases. However, each of these manipulations is thought to influence response in different ways across different cortical regions and a comprehensive study is required to interpret this literature. Here, human participants observed random-dot stimuli varying across a large range of contrast, coherence, and stimulus durations as we measured blood-oxygen-level dependent responses. We developed a framework for modeling these responses that quantifies their functional form and sensitivity across areas. Our framework demonstrates the sensitivity of all visual areas to each parameter, with early visual areas V1-V4 showing more parametric sensitivity to changes in contrast and V3A and the human middle temporal area to coherence. Our results suggest that while motion contrast, coherence, and duration share cortical representation, they are encoded with distinct functional forms and sensitivity. Thus, our quantitative framework serves as a reference for interpretation of the vast perceptual literature manipulating these parameters and shows that different manipulations of visibility will have different effects across human visual cortex and need to be interpreted accordingly. NEW & NOTEWORTHY Manipulations of motion visibility have served as a key tool for understanding the neural basis for visual perception. Here we measured human cortical response to changes in visibility across a comprehensive range of motion visibility parameters and modeled these with a quantitative framework. Our quantitative framework can be used as a reference for linking human cortical response to perception and underscores that different manipulations of motion visibility can have greatly different effects on cortical representation.

Task-dependent enhancement of facial expression and identity representations in human cortex

What cortical mechanisms allow humans to easily discern the expression or identity of a face? Subjects detected changes in expression or identity of a stream of dynamic faces while we measured BOLD responses from topographically and functionally defined areas throughout the visual hierarchy. Responses in dorsal areas increased during the expression task, whereas responses in ventral areas increased during the identity task, consistent with previous studies. Similar to ventral areas, early visual areas showed increased activity during the identity task. If visual responses are weighted by perceptual mechanisms according to their magnitude, these increased responses would lead to improved attentional selection of the task-appropriate facial aspect. Alternatively, increased responses could be a signature of a sensitivity enhancement mechanism that improves representations of the attended facial aspect. Consistent with the latter sensitivity enhancement mechanism, attending to expression led to enhanced decoding of exemplars of expression both in early visual and dorsal areas relative to attending identity. Similarly, decoding identity exemplars when attending to identity was improved in dorsal and ventral areas. We conclude that attending to expression or identity of dynamic faces is associated with increased selectivity in representations consistent with sensitivity enhancement.

Encoding model of temporal processing in human visual cortex

How is temporal information processed in human visual cortex? Visual input is relayed to V1 through segregated transient and sustained channels in the retina and lateral geniculate nucleus (LGN). However, there is intense debate as to how sustained and transient temporal channels contribute to visual processing beyond V1. The prevailing view associates transient processing predominately with motion-sensitive regions and sustained processing with ventral stream regions, while the opposing view suggests that both temporal channels contribute to neural processing beyond V1. Using fMRI, we measured cortical responses to time-varying stimuli and then implemented a two temporal channel-encoding model to evaluate the contributions of each channel. Different from the general linear model of fMRI that predicts responses directly from the stimulus, the encoding approach first models neural responses to the stimulus from which fMRI responses are derived. This encoding approach not only predicts cortical responses to time-varying stimuli from milliseconds to seconds but also, reveals differential contributions of temporal channels across visual cortex. Consistent with the prevailing view, motion-sensitive regions and adjacent lateral occipitotemporal regions are dominated by transient responses. However, ventral occipitotemporal regions are driven by both sustained and transient channels, with transient responses exceeding the sustained. These findings propose a rethinking of temporal processing in the ventral stream and suggest that transient processing may contribute to rapid extraction of the content of the visual input. Importantly, our encoding approach has vast implications, because it can be applied with fMRI to decipher neural computations in millisecond resolution in any part of the brain.

Compressive Temporal Summation in Human Visual Cortex

Combining sensory inputs over space and time is fundamental to vision. Population receptive field models have been successful in characterizing spatial encoding throughout the human visual pathways. A parallel question, how visual areas in the human brain process information distributed over time, has received less attention. One challenge is that the most widely used neuroimaging method, fMRI, has coarse temporal resolution compared with the time-scale of neural dynamics. Here, via carefully controlled temporally modulated stimuli, we show that information about temporal processing can be readily derived from fMRI signal amplitudes in male and female subjects. We find that all visual areas exhibit subadditive summation, whereby responses to longer stimuli are less than the linear prediction from briefer stimuli. We also find fMRI evidence that the neural response to two stimuli is reduced for brief interstimulus intervals (indicating adaptation). These effects are more pronounced in visual areas anterior to V1-V3. Finally, we develop a general model that shows how these effects can be captured with two simple operations: temporal summation followed by a compressive nonlinearity. This model operates for arbitrary temporal stimulation patterns and provides a simple and interpretable set of computations that can be used to characterize neural response properties across the visual hierarchy. Importantly, compressive temporal summation directly parallels earlier findings of compressive spatial summation in visual cortex describing responses to stimuli distributed across space. This indicates that, for space and time, cortex uses a similar processing strategy to achieve higher-level and increasingly invariant representations of the visual world.SIGNIFICANCE STATEMENT Combining sensory inputs over time is fundamental to seeing. Two important temporal phenomena are summation, the accumulation of sensory inputs over time, and adaptation, a response reduction for repeated or sustained stimuli. We investigated these phenomena in the human visual system using fMRI. We built predictive models that operate on arbitrary temporal patterns of stimulation using two simple computations: temporal summation followed by a compressive nonlinearity. Our new temporal compressive summation model captures (1) subadditive temporal summation, and (2) adaptation. We show that the model accounts for systematic differences in these phenomena across visual areas. Finally, we show that for space and time, the visual system uses a similar strategy to achieve increasingly invariant representations of the visual world.

Inverted Encoding Models of Human Population Response Conflate Noise and Neural Tuning Width

Channel-encoding models offer the ability to bridge different scales of neuronal measurement by interpreting population responses, typically measured with BOLD imaging in humans, as linear sums of groups of neurons (channels) tuned for visual stimulus properties. Inverting these models to form predicted channel responses from population measurements in humans seemingly offers the potential to infer neuronal tuning properties. Here, we test the ability to make inferences about neural tuning width from inverted encoding models. We examined contrast invariance of orientation selectivity in human V1 (both sexes) and found that inverting the encoding model resulted in channel response functions that became broader with lower contrast, thus apparently violating contrast invariance. Simulations showed that this broadening could be explained by contrast-invariant single-unit tuning with the measured decrease in response amplitude at lower contrast. The decrease in response lowers the signal-to-noise ratio of population responses that results in poorer population representation of orientation. Simulations further showed that increasing signal to noise makes channel response functions less sensitive to underlying neural tuning width, and in the limit of zero noise will reconstruct the channel function assumed by the model regardless of the bandwidth of single units. We conclude that our data are consistent with contrast-invariant orientation tuning in human V1. More generally, our results demonstrate that population selectivity measures obtained by encoding models can deviate substantially from the behavior of single units because they conflate neural tuning width and noise and are therefore better used to estimate the uncertainty of decoded stimulus properties.SIGNIFICANCE STATEMENT It is widely recognized that perceptual experience arises from large populations of neurons, rather than a few single units. Yet, much theory and experiment have examined links between single units and perception. Encoding models offer a way to bridge this gap by explicitly interpreting population activity as the aggregate response of many single neurons with known tuning properties. Here we use this approach to examine contrast-invariant orientation tuning of human V1. We show with experiment and modeling that due to lower signal to noise, contrast-invariant orientation tuning of single units manifests in population response functions that broaden at lower contrast, rather than remain contrast-invariant. These results highlight the need for explicit quantitative modeling when making a reverse inference from population response profiles to single-unit responses.

BOLD neurovascular coupling does not change significantly with normal aging

Studies of cognitive function that compare the blood oxygenation level dependent (BOLD) signal across age groups often require the assumption that neurovascular coupling does not change with age. Tests of this assumption have produced mixed results regarding the strength of the coupling and its relative time course. Using deconvolution, we found that age does not have a significant effect on the time course of the hemodynamic impulse response function or on the slope of the BOLD versus stimulus duration relationship. These results suggest that in cognitive studies of healthy aging, group differences in BOLD activation are likely due to age-related changes in cognitive-neural interactions and information processing rather than to impairments in neurovascular coupling. Hum Brain Mapp 38:3538-3551, 2017. © 2017 Wiley Periodicals, Inc.

At risk of being risky: The relationship between "brain age" under emotional states and risk preference

Developmental differences regarding decision making are often reported in the absence of emotional stimuli and without context, failing to explain why some individuals are more likely to have a greater inclination toward risk. The current study (N=212; 10-25y) examined the influence of emotional context on underlying functional brain connectivity over development and its impact on risk preference. Using functional imaging data in a neutral brain-state we first identify the "brain age" of a given individual then validate it with an independent measure of cortical thickness. We then show, on average, that "brain age" across the group during the teen years has the propensity to look younger in emotional contexts. Further, we show this phenotype (i.e. a younger brain age in emotional contexts) relates to a group mean difference in risk perception - a pattern exemplified greatest in young-adults (ages 18-21). The results are suggestive of a specified functional brain phenotype that relates to being at "risk to be risky."

Validation of periodic fMRI signals in response to wearable tactile stimulation

To map cortical representations of the body, we recently developed a wearable technology for automatic tactile stimulation in human functional magnetic resonance imaging (fMRI) experiments. In a two-condition block design experiment, air puffs were delivered to the face and hands periodically. Surface-based regions of interest (S-ROIs) were initially identified by thresholding a linear statistical measure of signal-to-noise ratio of periodic response. Across subjects, S-ROIs were found in the frontal, primary sensorimotor, posterior parietal, insular, temporal, cingulate, and occipital cortices. To validate and differentiate these S-ROIs, we develop a measure of temporal stability of response based on the assumption that a periodic stimulation evokes stable (low-variance) periodic fMRI signals throughout the entire scan. Toward this end, we apply time-frequency analysis to fMRI time series and use circular statistics to characterize the distribution of phase angles for data selection. We then assess the temporal variability of a periodic signal by measuring the path length of its trajectory in the complex plane. Both within and outside the primary sensorimotor cortex, S-ROIs with high temporal variability and deviant phase angles are rejected. A surface-based probabilistic group-average map is constructed for spatial screening of S-ROIs with low to moderate temporal variability in non-sensorimotor regions. Areas commonly activated across subjects are also summarized in the group-average map. In summary, this study demonstrates that analyzing temporal characteristics of the entire fMRI time series is essential for second-level selection and interpretation of S-ROIs initially defined by an overall linear statistical measure.

Cortical interactions during the resolution of information processing demands in autism spectrum disorders

INTRODUCTION

Our flexible and adaptive interactions with the environment are guided by our individual representation of the physical world, estimated through sensation and evaluation of available information against prior knowledge. When linking sensory evidence with higher-level expectations for action, the central nervous system (CNS) in typically developing (TD) individuals relies in part on distributed and interacting cortical regions to communicate neuronal signals flexibly across the brain. Increasing evidence suggests that the balance between levels of signal and noise during information processing may be disrupted in individuals with Autism Spectrum Disorders (ASD).

METHODS

Participants with and without ASD performed a visuospatial interference task while undergoing functional Magnetic Resonance Imaging (fMRI). We empirically estimated parameters characterizing participants' latencies and their subtle fluctuations (noise accumulation) over the 16-min scan. We modeled hemodynamic activation and used seed-based analyses of neural coupling to study dysfunction in interference-specific connectivity in a subset of ASD participants who were nonparametrically matched to TD participants on age, male-to-female ratio, and magnitude of movement during the scan.

RESULTS

Stochastic patterns of response fluctuations reveal significantly higher noise-to-signal levels and a more random and noisy structure in ASD versus TD participants, and in particular ASD adults who have the greatest clinical autistic deficits. While individuals with ASD show an overall weaker modulation of interference-specific functional connectivity relative to TD individuals, in particular between the seeds of Anterior Cingulate Cortex (ACC) and Inferior Parietal Sulcus (IPS) and the rest of the brain, we found that in ASD, higher uncertainty during the task is linked to increased interference-specific coupling between bilateral anterior insula and prefrontal cortex.

CONCLUSIONS

Subtle and informative differences in the structure of experiencing information exist between ASD and TD individuals. Our findings reveal in ASD an atypical capacity to apply previously perceived information in a manner optimal for adaptive functioning, plausibly revealing suboptimal message-passing across the CNS.

Modeling Neurovascular Coupling from Clustered Parameter Sets for Multimodal EEG-NIRS

Despite significant improvements in neuroimaging technologies and analysis methods, the fundamental relationship between local changes in cerebral hemodynamics and the underlying neural activity remains largely unknown. In this study, a data driven approach is proposed for modeling this neurovascular coupling relationship from simultaneously acquired electroencephalographic (EEG) and near-infrared spectroscopic (NIRS) data. The approach uses gamma transfer functions to map EEG spectral envelopes that reflect time-varying power variations in neural rhythms to hemodynamics measured with NIRS during median nerve stimulation. The approach is evaluated first with simulated EEG-NIRS data and then by applying the method to experimental EEG-NIRS data measured from 3 human subjects. Results from the experimental data indicate that the neurovascular coupling relationship can be modeled using multiple sets of gamma transfer functions. By applying cluster analysis, statistically significant parameter sets were found to predict NIRS hemodynamics from EEG spectral envelopes. All subjects were found to have significant clustered parameters (P < 0.05) for EEG-NIRS data fitted using gamma transfer functions. These results suggest that the use of gamma transfer functions followed by cluster analysis of the resulting parameter sets may provide insights into neurovascular coupling in human neuroimaging data.

An investigation of positive and inverted hemodynamic response functions across multiple visual areas

Recent studies have demonstrated significant regional variability in the hemodynamic response function (HRF), highlighting the difficulty of correctly interpreting functional MRI (fMRI) data without proper modeling of the HRF. The focus of this study was to investigate the HRF variability within visual cortex. The HRF was estimated for a number of cortical visual areas by deconvolution of fMRI blood oxygenation level dependent (BOLD) responses to brief, large-field visual stimulation. Significant HRF variation was found across visual areas V1, V2, V3, V4, VO-1,2, V3AB, IPS-0,1,2,3, LO-1,2, and TO-1,2. Additionally, a subpopulation of voxels was identified that exhibited an impulse response waveform that was similar, but not identical, to an inverted version of the commonly described and modeled positive HRF. These voxels were found within the retinotopic confines of the stimulus and were intermixed with those showing positive responses. The spatial distribution and variability of these HRFs suggest a vascular origin for the inverted waveforms. We suggest that the polarity of the HRF is a separate factor that is independent of the suppressive or activating nature of the underlying neuronal activity. Correctly modeling the polarity of the HRF allows one to recover an estimate of the underlying neuronal activity rather than discard the responses from these voxels on the assumption that they are artifactual. We demonstrate this approach on phase-encoded retinotopic mapping data as an example of the benefits of accurately modeling the HRF during the analysis of fMRI data.

Optimizing the general linear model for functional near-infrared spectroscopy: an adaptive hemodynamic response function approach

An increasing number of functional near-infrared spectroscopy (fNIRS) studies utilize a general linear model (GLM) approach, which serves as a standard statistical method for functional magnetic resonance imaging (fMRI) data analysis. While fMRI solely measures the blood oxygen level dependent (BOLD) signal, fNIRS measures the changes of oxy-hemoglobin (oxy-Hb) and deoxy-hemoglobin (deoxy-Hb) signals at a temporal resolution severalfold higher. This suggests the necessity of adjusting the temporal parameters of a GLM for fNIRS signals. Thus, we devised a GLM-based method utilizing an adaptive hemodynamic response function (HRF). We sought the optimum temporal parameters to best explain the observed time series data during verbal fluency and naming tasks. The peak delay of the HRF was systematically changed to achieve the best-fit model for the observed oxy- and deoxy-Hb time series data. The optimized peak delay showed different values for each Hb signal and task. When the optimized peak delays were adopted, the deoxy-Hb data yielded comparable activations with similar statistical power and spatial patterns to oxy-Hb data. The adaptive HRF method could suitably explain the behaviors of both Hb parameters during tasks with the different cognitive loads during a time course, and thus would serve as an objective method to fully utilize the temporal structures of all fNIRS data.

A mixed L2 norm regularized HRF estimation method for rapid event-related fMRI experiments

Brain state decoding or "mind reading" via multivoxel pattern analysis (MVPA) has become a popular focus of functional magnetic resonance imaging (fMRI) studies. In brain decoding, stimulus presentation rate is increased as fast as possible to collect many training samples and obtain an effective and reliable classifier or computational model. However, for extremely rapid event-related experiments, the blood-oxygen-level-dependent (BOLD) signals evoked by adjacent trials are heavily overlapped in the time domain. Thus, identifying trial-specific BOLD responses is difficult. In addition, voxel-specific hemodynamic response function (HRF), which is useful in MVPA, should be used in estimation to decrease the loss of weak information across voxels and obtain fine-grained spatial information. Regularization methods have been widely used to increase the efficiency of HRF estimates. In this study, we propose a regularization framework called mixed L2 norm regularization. This framework involves Tikhonov regularization and an additional L2 norm regularization term to calculate reliable HRF estimates. This technique improves the accuracy of HRF estimates and significantly increases the classification accuracy of the brain decoding task when applied to a rapid event-related four-category object classification experiment. At last, some essential issues such as the impact of low-frequency fluctuation (LFF) and the influence of smoothing are discussed for rapid event-related experiments.

Convolution models for induced electromagnetic responses

In Kilner et al. [Kilner, J.M., Kiebel, S.J., Friston, K.J., 2005. Applications of random field theory to electrophysiology. Neurosci. Lett. 374, 174-178.] we described a fairly general analysis of induced responses-in electromagnetic brain signals-using the summary statistic approach and statistical parametric mapping. This involves localising induced responses-in peristimulus time and frequency-by testing for effects in time-frequency images that summarise the response of each subject to each trial type. Conventionally, these time-frequency summaries are estimated using post-hoc averaging of epoched data. However, post-hoc averaging of this sort fails when the induced responses overlap or when there are multiple response components that have variable timing within each trial (for example stimulus and response components associated with different reaction times). In these situations, it is advantageous to estimate response components using a convolution model of the sort that is standard in the analysis of fMRI time series. In this paper, we describe one such approach, based upon ordinary least squares deconvolution of induced responses to input functions encoding the onset of different components within each trial. There are a number of fundamental advantages to this approach: for example; (i) one can disambiguate induced responses to stimulus onsets and variably timed responses; (ii) one can test for the modulation of induced responses-over peristimulus time and frequency-by parametric experimental factors and (iii) one can gracefully handle confounds-such as slow drifts in power-by including them in the model. In what follows, we consider optimal forms for convolution models of induced responses, in terms of impulse response basis function sets and illustrate the utility of deconvolution estimators using simulated and real MEG data.