Identifying Neural Drivers with Functional MRI: An Electrophysiological Validation

Whether functional magnetic resonance imaging (fMRI) allows the identification of neural drivers remains an open question of particular importance to refine physiological and neuropsychological models of the brain, and/or to understand neurophysiopathology. Here, in a rat model of absence epilepsy showing spontaneous spike-and-wave discharges originating from the first somatosensory cortex (S1BF), we performed simultaneous electroencephalographic (EEG) and fMRI measurements, and subsequent intracerebral EEG (iEEG) recordings in regions strongly activated in fMRI (S1BF, thalamus, and striatum). fMRI connectivity was determined from fMRI time series directly and from hidden state variables using a measure of Granger causality and Dynamic Causal Modelling that relates synaptic activity to fMRI. fMRI connectivity was compared to directed functional coupling estimated from iEEG using asymmetry in generalised synchronisation metrics. The neural driver of spike-and-wave discharges was estimated in S1BF from iEEG, and from fMRI only when hemodynamic effects were explicitly removed. Functional connectivity analysis applied directly on fMRI signals failed because hemodynamics varied between regions, rendering temporal precedence irrelevant. This paper provides the first experimental substantiation of the theoretical possibility to improve interregional coupling estimation from hidden neural states of fMRI. As such, it has important implications for future studies on brain connectivity using functional neuroimaging.     

Network Analysis of the Default Mode Network Using Functional Connectivity MRI in Temporal Lobe Epilepsy

Functional connectivity MRI (fcMRI) is an fMRI method that examines the connectivity of different brain areas based on the correlation of BOLD signal fluctuations over time. Temporal Lobe Epilepsy (TLE) is the most common type of adult epilepsy and involves multiple brain networks. The default mode network (DMN) is involved in conscious, resting state cognition and is thought to be affected in TLE where seizures cause impairment of consciousness. The DMN in epilepsy was examined using seed based fcMRI. The anterior and posterior hubs of the DMN were used as seeds in this analysis. The results show a disconnection between the anterior and posterior hubs of the DMN in TLE during the basal state. In addition, increased DMN connectivity to other brain regions in left TLE along with decreased connectivity in right TLE is revealed. The analysis demonstrates how seed-based fcMRI can be used to probe cerebral networks in brain disorders such as TLE.     

Long-Term Effects of Temporal Lobe Epilepsy on Local Neural Networks: A Graph Theoretical Analysis of Corticography Recordings

Purpose

Pharmaco-resistant temporal lobe epilepsy (TLE) is often treated with surgical intervention at some point. As epilepsy surgery is considered a last resort by most physicians, a long history of epileptic seizures prior to surgery is not uncommon. Little is known about the effects of ongoing TLE on neural functioning. A better understanding of these effects might influence the moment of surgical intervention. Functional connectivity (interaction between spatially distributed brain areas) and network structure (integration and segregation of information processing) are thought to be essential for optimal brain functioning. We report on the impact of TLE duration on temporal lobe functional connectivity and network characteristics.

Methods

Functional connectivity of the temporal lobe at the time of surgery was assessed by means of interictal electrocorticography (ECoG) recordings of 27 TLE patients by using the phase lag index (PLI). Graphs (abstract network representations) were reconstructed from the PLI matrix and characterized by the clustering coefficient C (local clustering), the path length L (overall network interconnectedness), and the “small world index” S (network configuration).

Results

Functional connectivity (average PLI), clustering coefficients, and the small world index were negatively correlated with TLE duration in the broad frequency band (0.5–48 Hz).

Discussion

Temporal lobe functional connectivity is lower in patients with longer TLE history, and longer TLE duration is correlated with more random network configuration. Our findings suggest that the neural networks of TLE patients become more pathological over time, possibly due to temporal lobe changes associated with long-standing lesional epilepsy.     

Increased Intrinsic Connectivity of the Default Mode Network in Temporal Lobe Epilepsy: Evidence from Resting-State MEG Recordings

The electrophysiological signature of resting state oscillatory functional connectivity within the default mode network (DMN) during spike-free periods in temporal lobe epilepsy (TLE) remains unclear. Using magnetoencephalographic (MEG) recordings, this study investigated how the connectivity within the DMN was altered in TLE, and we examined the effect of lateralized TLE on functional connectivity. Sixteen medically intractable TLE patients and 22 controls participated in this study. Whole-scalp 306-channel MEG epochs without interictal spikes generated from both MEG and EEG data were analyzed using a minimum norm estimate (MNE) and source-based imaginary coherence analysis. With this processing, we obtained the cortical activation and functional connectivity within the DMN. The functional connectivity was increased between DMN and the right medial temporal (MT) region at the delta band and between DMN and the bilateral anterior cingulate cortex (ACC) regions at the theta band. The functional change was associated with the lateralization of TLE. The right TLE showed enhanced DMN connectivity with the right MT while the left TLE demonstrated increased DMN connectivity with the bilateral MT. There was no lateralization effect of TLE upon the DMN connectivity with ACC. These findings suggest that the resting-state functional connectivity within the DMN is reinforced in temporal lobe epilepsy during spike-free periods. Future studies are needed to examine if the altered functional connectivity can be used as a biomarker for treatment responses, cognitive dysfunction and prognosis in patients with TLE.     

Resting State Brain Function Analysis Using Concurrent BOLD in ASL Perfusion fMRI

The past decade has seen astounding discoveries about resting-state brain activity patterns in normal brain as well as their alterations in brain diseases. While the vast majority of resting-state studies are based on the blood-oxygen-level-dependent (BOLD) functional MRI (fMRI), arterial spin labeling (ASL) perfusion fMRI can simultaneously capture BOLD and cerebral blood flow (CBF) signals, providing a unique opportunity for assessing resting brain functions with concurrent BOLD (ccBOLD) and CBF signals. Before taking that benefit, it is necessary to validate the utility of ccBOLD signal for resting-state analysis using conventional BOLD (cvBOLD) signal acquired without ASL modulations. To address this technical issue, resting cvBOLD and ASL perfusion MRI were acquired from a large cohort (n = 89) of healthy subjects. Four widely used resting-state brain function analyses were conducted and compared between the two types of BOLD signal, including the posterior cingulate cortex (PCC) seed-based functional connectivity (FC) analysis, independent component analysis (ICA), analysis of amplitude of low frequency fluctuation (ALFF), and analysis of regional homogeneity (ReHo). Consistent default mode network (DMN) as well as other resting-state networks (RSNs) were observed from cvBOLD and ccBOLD using PCC-FC analysis and ICA. ALFF from both modalities were the same for most of brain regions but were different in peripheral regions suffering from the susceptibility gradients induced signal drop. ReHo showed difference in many brain regions, likely reflecting the SNR and resolution differences between the two BOLD modalities. The DMN and auditory networks showed highest CBF values among all RSNs. These results demonstrated the feasibility of ASL perfusion MRI for assessing resting brain functions using its concurrent BOLD in addition to CBF signal, which provides a potentially useful way to maximize the utility of ASL perfusion MRI.     

Measuring Granger Causality between cortical regions from voxelwise fMRI BOLD signals with LASSO

Functional brain network studies using the Blood Oxygen-Level Dependent (BOLD) signal from functional Magnetic Resonance Imaging (fMRI) are becoming increasingly prevalent in research on the neural basis of human cognition. An important problem in functional brain network analysis is to understand directed functional interactions between brain regions during cognitive performance. This problem has important implications for understanding top-down influences from frontal and parietal control regions to visual occipital cortex in visuospatial attention, the goal motivating the present study. A common approach to measuring directed functional interactions between two brain regions is to first create nodal signals by averaging the BOLD signals of all the voxels in each region, and to then measure directed functional interactions between the nodal signals. Another approach, that avoids averaging, is to measure directed functional interactions between all pairwise combinations of voxels in the two regions. Here we employ an alternative approach that avoids the drawbacks of both averaging and pairwise voxel measures. In this approach, we first use the Least Absolute Shrinkage Selection Operator (LASSO) to pre-select voxels for analysis, then compute a Multivariate Vector AutoRegressive (MVAR) model from the time series of the selected voxels, and finally compute summary Granger Causality (GC) statistics from the model to represent directed interregional interactions. We demonstrate the effectiveness of this approach on both simulated and empirical fMRI data. We also show that averaging regional BOLD activity to create a nodal signal may lead to biased GC estimation of directed interregional interactions. The approach presented here makes it feasible to compute GC between brain regions without the need for averaging. Our results suggest that in the analysis of functional brain networks, careful consideration must be given to the way that network nodes and edges are defined because those definitions may have important implications for the validity of the analysis.     

Concordance of Epileptic Networks Associated with Epileptic Spikes Measured by High-Density EEG and Fast fMRI

Objective

The present study aims to investigate whether a newly developed fast fMRI called MREG (magnetic resonance encephalography) measures metabolic changes related to interictal epileptic discharges (IED). For this purpose BOLD changes are correlated with the IED distribution and variability.

Methods

Patients with focal epilepsy underwent EEG-MREG using a 64 channel cap. IED voltage maps were generated using 32 and 64 channels and compared regarding their correspondence to the BOLD response. The extents of IEDs (defined as number of channels with >50% of maximum IED negativity) were correlated with the extents of positive and negative BOLD responses. Differences in inter-spike variability were investigated between interictal epileptic discharges (IED) sets with and without concordant positive or negative BOLD responses.

Results

17 patients showed 32 separate IED types. In 50% of IED types the BOLD changes could be confirmed by another independent imaging method. The IED extent significantly correlated with the positive BOLD extent (p = 0.04). In 6 patients the 64-channel EEG voltage maps better reflected the positive or negative BOLD response than the 32-channel EEG; in all others no difference was seen. Inter-spike variability was significantly lower in IED sets with than without concordant positive or negative BOLD responses (with p = 0.04).

Significance

Higher density EEG and fast fMRI seem to improve the value of EEG-fMRI in epilepsy. The correlation of positive BOLD and IED extent could suggest that widespread BOLD responses reflect the IED network. Inter-spike variability influences the likelihood to find IED concordant positive or negative BOLD responses, which is why single IED analysis may be promising.     

Voxel-Wise Motion Artifacts in Population-Level Whole-Brain Connectivity Analysis of Resting-State fMRI

Functional Magnetic Resonance Imaging (fMRI) based brain connectivity analysis maps the functional networks of the brain by estimating the degree of synchronous neuronal activity between brain regions. Recent studies have demonstrated that “resting-state” fMRI-based brain connectivity conclusions may be erroneous when motion artifacts have a differential effect on fMRI BOLD signals for between group comparisons. A potential explanation could be that in-scanner displacement, due to rotational components, is not spatially constant in the whole brain. However, this localized nature of motion artifacts is poorly understood and is rarely considered in brain connectivity studies. In this study, we initially demonstrate the local correspondence between head displacement and the changes in the resting-state fMRI BOLD signal. Than, we investigate how connectivity strength is affected by the population-level variation in the spatial pattern of regional displacement. We introduce Regional Displacement Interaction (RDI), a new covariate parameter set for second-level connectivity analysis and demonstrate its effectiveness in reducing motion related confounds in comparisons of groups with different voxel-vise displacement pattern and preprocessed using various nuisance regression methods. The effect of using RDI as second-level covariate is than demonstrated in autism-related group comparisons. The relationship between the proposed method and some of the prevailing subject-level nuisance regression techniques is evaluated. Our results show that, depending on experimental design, treating in-scanner head motion as a global confound may not be appropriate. The degree of displacement is highly variable among various brain regions, both within and between subjects. These regional differences bias correlation-based measures of brain connectivity. The inclusion of the proposed second-level covariate into the analysis successfully reduces artifactual motion-related group differences and preserves real neuronal differences, as demonstrated by the autism-related comparisons.     

Epileptic Discharges Affect the Default Mode Network – fMRI and Intracerebral EEG Evidence

Functional neuroimaging studies of epilepsy patients often show, at the time of epileptic activity, deactivation in default mode network (DMN) regions, which is hypothesized to reflect altered consciousness. We aimed to study the metabolic and electrophysiological correlates of these changes in the DMN regions. We studied six epilepsy patients that underwent scalp EEG-fMRI and later stereotaxic intracerebral EEG (SEEG) sampling regions of DMN (posterior cingulate cortex, Pre-cuneus, inferior parietal lobule, medial prefrontal cortex and dorsolateral frontal cortex) as well as non-DMN regions. SEEG recordings were subject to frequency analyses comparing sections with interictal epileptic discharges (IED) to IED-free baselines in the IED-generating region, DMN and non-DMN regions. EEG-fMRI and SEEG were obtained at rest. During IEDs, EEG-fMRI demonstrated deactivation in various DMN nodes in 5 of 6 patients, most frequently the pre-cuneus and inferior parietal lobule, and less frequently the other DMN nodes. SEEG analyses demonstrated decrease in gamma power (50–150 Hz), and increase in the power of lower frequencies (<30 Hz) at times of IEDs, in at least one DMN node in all patients. These changes were not apparent in the non-DMN regions. We demonstrate that, at the time of IEDs, DMN regions decrease their metabolic demand and undergo an EEG change consisting of decreased gamma and increased lower frequencies. These findings, specific to DMN regions, confirm in a pathological condition a direct relationship between DMN BOLD activity and EEG activity. They indicate that epileptic activity affects the DMN, and therefore may momentarily reduce the consciousness level and cognitive reserve.     

Resting State BOLD Functional Connectivity at 3T: Spin Echo versus Gradient Echo EPI

Previous evidence showed that, due to refocusing of static dephasing effects around large vessels, spin-echo (SE) BOLD signals offer an increased linearity and promptness with respect to gradient-echo (GE) acquisition, even at low field. These characteristics suggest that, despite the reduced sensitivity, SE fMRI might also provide a potential benefit when investigating spontaneous fluctuations of brain activity. However, there are no reports on the application of spin-echo fMRI for connectivity studies at low field. In this study we compared resting state functional connectivity as measured with GE and SE EPI sequences at 3T. Main results showed that, within subject, the GE sensitivity is overall larger with respect to that of SE, but to a less extent than previously reported for activation studies. Noteworthy, the reduced sensitivity of SE was counterbalanced by a reduced inter-subject variability, resulting in comparable group statistical connectivity maps for the two sequences. Furthermore, the SE method performed better in the ventral portion of the default mode network, a region affected by signal dropout in standard GE acquisition. Future studies should clarify if these features of the SE BOLD signal can be beneficial to distinguish subtle variations of functional connectivity across different populations and/or treatments when vascular confounds or regions affected by signal dropout can be a critical issue.     

Simultaneous BOLD fMRI and fiber-optic calcium recording in rat neocortex

Functional magnetic resonance imaging (fMRI) based on blood oxygen level-dependent (BOLD) contrast is widely used for probing brain activity, but its relationship to underlying neural activity remains elusive. Here, we combined fMRI with fiber-optic recordings of fluorescent calcium indicator signals to investigate this relationship in rat somatosensory cortex. Electrical forepaw stimulation (1-10 Hz) evoked fast calcium signals of neuronal origin that showed frequency-dependent adaptation. Additionally, slower calcium signals occurred in astrocyte networks, as verified by astrocyte-specific staining and two-photon microscopy. Without apparent glia activation, we could predict BOLD responses well from simultaneously recorded fiber-optic signals, assuming an impulse response function and taking into account neuronal adaptation. In cases with glia activation, we uncovered additional prolonged BOLD signal components. Our findings highlight the complexity of fMRI BOLD signals, involving both neuronal and glial activity. Combined fMRI and fiber-optic recordings should help to clarify cellular mechanisms underlying BOLD signals.     

Protective Effect of Resveratrol against Kainate-induced Temporal Lobe Epilepsy in Rats

Resveratrol (Res) is a phytoalexin produced naturally by several plants, which has multi functional effects such as neuroprotection, anti-inflammatory, and anti-cancer. The present study was to evaluate a possible anti-epileptic effect of Res against kainate-induced temporal lobe epilepsy (TLE) in rat. We performed behavior monitoring, intracranial electroencepholography (IEEG) recording, histological analysis, and Western blotting to evaluate the anti-epilepsy effect of Res in kainate-induced epileptic rats. Res decreased the frequency of spontaneous seizures and inhibited the epileptiform discharges. Moreover, Res could protect neurons against kainate-induced neuronal cell death in CA1 and CA3a regions and depressed mossy fiber sprouting, which are general histological characteristics both in TLE patients and animal models. Western blot revealed that the expression level of kainate receptors (KARs) in hippocampus was reduced in Res-administrated rats compared to that in epileptic ones. These results suggest that Res is a potent anti-epilepsy agent, which protects against epileptogenesis and progression of the kainate-induced TLE animal. The authors Z. Wu and Q. Xu contributed equally to this work.     

Self-regulation of local brain activity using real-time functional magnetic resonance imaging (fMRI)

Functional magnetic resonance imaging (fMRI) measures the blood oxygen level-dependent (BOLD) signal related to neuronal activity. So far, this technique has been limited by time-consuming data analysis impeding on-line analysis. In particular, no brain-computer interface (BCI) was available which provided on-line feedback to learn physiological self-regulation of the BOLD signal. Recently, studies have shown that fMRI feedback is feasible and facilitates voluntary control of brain activity. Here we review these studies to make the fMRI feedback methodology accessible to a broader scientific community such as researchers concerned with functional brain imaging and the neurobiology of learning. Methodological and conceptual limitations were substantially reduced by artefact control, sensitivity improvements, real-time algorithms, and adapted experimental designs. Physiological self-regulation of the local BOLD response is a new paradigm for cognitive neuroscience to study brain plasticity and the functional relevance of regulated brain areas by modification of behaviour. Voluntary control of abnormal activity in circumscribed brain areas may even be applied as psychophysiological treatment.     

The neural basis of the blood-oxygen-level-dependent functional magnetic resonance imaging signal

Magnetic resonance imaging (MRI) has rapidly become an important tool in clinical medicine and biological research. Its functional variant (functional magnetic resonance imaging; fMRI) is currently the most widely used method for brain mapping and studying the neural basis of human cognition. While the method is widespread, there is insufficient knowledge of the physiological basis of the fMRI signal to interpret the data confidently with respect to neural activity. This paper reviews the basic principles of MRI and fMRI, and subsequently discusses in some detail the relationship between the blood-oxygen-level-dependent (BOLD) fMRI signal and the neural activity elicited during sensory stimulation. To examine this relationship, we conducted the first simultaneous intracortical recordings of neural signals and BOLD responses. Depending on the temporal characteristics of the stimulus, a moderate to strong correlation was found between the neural activity measured with microelectrodes and the BOLD signal averaged over a small area around the microelectrode tips. However, the BOLD signal had significantly higher variability than the neural activity, indicating that human fMRI combined with traditional statistical methods underestimates the reliability of the neuronal activity. To understand the relative contribution of several types of neuronal signals to the haemodynamic response, we compared local field potentials (LFPs), single- and multi-unit activity (MUA) with high spatio-temporal fMRI responses recorded simultaneously in monkey visual cortex. At recording sites characterized by transient responses, only the LFP signal was significantly correlated with the haemodynamic response. Furthermore, the LFPs had the largest magnitude signal and linear systems analysis showed that the LFPs were better than the MUAs at predicting the fMRI responses. These findings, together with an analysis of the neural signals, indicate that the BOLD signal primarily measures the input and processing of neuronal information within a region and not the output signal transmitted to other brain regions.     

Functional Connectivity Estimated from Intracranial EEG Predicts Surgical Outcome in Intractable Temporal Lobe Epilepsy

This project aimed to determine if a correlation-based measure of functional connectivity can identify epileptogenic zones from intracranial EEG signals, as well as to investigate the prognostic significance of such a measure on seizure outcome following temporal lobe lobectomy. To this end, we retrospectively analyzed 23 adult patients with intractable temporal lobe epilepsy (TLE) who underwent an invasive stereo-EEG (SEEG) evaluation between January 2009 year and January 2012. A follow-up of at least one year was required. The primary outcome measure was complete seizure-freedom at last follow-up. Functional connectivity between two areas in the temporal lobe that were sampled by two SEEG electrode contacts was defined as Pearson’s correlation coefficient of interictal activity between those areas. SEEG signals were filtered between 5 and 50 Hz prior to computing this correlation. The mean and standard deviation of the off diagonal elements in the connectivity matrix were also calculated. Analysis of the mean and standard deviation of the functional connections for each patient reveals that 90% of the patients who had weak and homogenous connections were seizure free one year after temporal lobectomy, whereas 85% of the patients who had stronger and more heterogeneous connections within the temporal lobe had recurrence of seizures. This suggests that temporal lobectomy is ineffective in preventing seizure recurrence for patients in whom the temporal lobe is characterized by weakly connected, homogenous networks. This pilot study shows promising potential of a simple measure of functional brain connectivity to identify epileptogenicity and predict the outcome of epilepsy surgery.     

Disrupted Ipsilateral Network Connectivity in Temporal Lobe Epilepsy

ObjectiveThe current practice under which patients with refractory epilepsy are surgically treated is based mainly on the identification of specific cortical areas, mainly the epileptogenic zone, which is believed to be responsible for generation of seizures. A better understanding of the whole epileptic network and its components and properties is required before more effective and less invasive therapies can be developed. The aim of the present study was to partially characterize the evolution of the functional network during the preictal-ictal transition in partial seizures in patients with temporal lobe epilepsy (TLE).MethodsScalp and foramen ovale (FOE) recordings from twenty-two TLE patients were analyzed under the complex network perspective. The density of links, average path length, average clustering coefficient, and modularity were calculated during the preictal and the ictal stages. Both linear–Pearson correlation–and non-linear–phase synchronization–measures were used as proxies of functional connectivity between the electrode locations areas. The transition from one stage to the other was evaluated in the whole network and in the mesial sub-networks. The results were compared with a voltage-dependent measure, namely, the spectral entropy.ResultsChanges in the global functional network during the transition from the preictal to the ictal stage show, in the linear case, that in sixteen cases (72.7%) the density of the links increased during the seizure, with a decrease in the average path length in fifteen cases (68.1%). There was also a preictal and ictal imbalance in functional connectivity during both stages (77.2% to 86.3%). The SE dropped during the seizure in 95.4% of the cases, but did not show any tendency towards lateralization. When using the nonlinear measure of functional connectivity, the phase synchronization, similar results were obtained.ConclusionsIn TLE patients, the transition to the ictal stage is accompanied by increasing global synchronization and a more ordered spectral content of the signals, indicated by lower spectral entropy. The interictal connectivity imbalance (lower ipsilateral connectivity) is sustained during the seizure, irrespective of any appreciable imbalance in the spectral entropy of the mesial recordings.     

Reproducibility of swallow-induced cortical BOLD positive and negative fMRI activity

Functional MRI (fMRI) studies have demonstrated that a number of brain regions (cingulate, insula, prefrontal, and sensory/motor cortices) display blood oxygen level-dependent (BOLD) positive activity during swallow. Negative BOLD activations and reproducibility of these activations have not been systematically studied. The aim of our study was to investigate the reproducibility of swallow-related cortical positive and negative BOLD activity across different fMRI sessions. We studied 16 healthy volunteers utilizing an fMRI event-related analysis. Individual analysis using a general linear model was used to remove undesirable signal changes correlated with motion, white matter, and cerebrospinal fluid. The group analysis used a mixed-effects multilevel model to identify active cortical regions. The volume and magnitude of a BOLD signal within each cluster was compared between the two study sessions. All subjects showed significant clustered BOLD activity within the known areas of cortical swallowing network across both sessions. The cross-correlation coefficient of percent fMRI signal change and the number of activated voxels across both positive and negative BOLD networks were similar between the two studies (r ≥ 0.87, P < 0.0001). Swallow-associated negative BOLD activity was comparable to the well-defined "default-mode" network, and positive BOLD activity had noticeable overlap with the previously described "task-positive" network. Swallow activates two parallel cortical networks. These include a positive and a negative BOLD network, respectively, correlated and anticorrelated with swallow stimulus. Group cortical activity maps, as well as extent and amplitude of activity induced by volitional swallowing in the cortical swallowing network, are reproducible between study sessions.     

The Spectral Diversity of Resting-State Fluctuations in the Human Brain

In order to assess whole-brain resting-state fluctuations at a wide range of frequencies, resting-state fMRI data of 20 healthy subjects were acquired using a multiband EPI sequence with a low TR (354 ms) and compared to 20 resting-state datasets from standard, high-TR (1800 ms) EPI scans. The spatial distribution of fluctuations in various frequency ranges are analyzed along with the spectra of the time-series in voxels from different regions of interest. Functional connectivity specific to different frequency ranges (<0.1 Hz; 0.1–0.25 Hz; 0.25–0.75 Hz; 0.75–1.4 Hz) was computed for both the low-TR and (for the two lower-frequency ranges) the high-TR datasets using bandpass filters. In the low-TR data, cortical regions exhibited highest contribution of low-frequency fluctuations and the most marked low-frequency peak in the spectrum, while the time courses in subcortical grey matter regions as well as the insula were strongly contaminated by high-frequency signals. White matter and CSF regions had highest contribution of high-frequency fluctuations and a mostly flat power spectrum. In the high-TR data, the basic patterns of the low-TR data can be recognized, but the high-frequency proportions of the signal fluctuations are folded into the low frequency range, thus obfuscating the low-frequency dynamics. Regions with higher proportion of high-frequency oscillations in the low-TR data showed flatter power spectra in the high-TR data due to aliasing of the high-frequency signal components, leading to loss of specificity in the signal from these regions in high-TR data. Functional connectivity analyses showed that there are correlations between resting-state signal fluctuations of distant brain regions even at high frequencies, which can be measured using low-TR fMRI. On the other hand, in the high-TR data, loss of specificity of measured fluctuations leads to lower sensitivity in detecting functional connectivity. This underlines the advantages of low-TR EPI sequences for resting-state and potentially also task-related fMRI experiments.     

Biophysical models of fMRI responses

Functional magnetic resonance imaging (fMRI) is used to investigate where the neural implementation of specific cognitive processes occurs. The standard approach uses linear convolution models that relate experimentally designed inputs, through a haemodynamic response function, to observed blood oxygen level dependent (BOLD) signals. Such models are, however, blind to the causal mechanisms that underlie observed BOLD responses. Recent developments have focused on how BOLD responses are generated and include biophysical input-state-output models with neural and haemodynamic state equations and models of functional integration that explain local dynamics through interactions with remote areas. Forward models with parameters at the neural level, such as dynamic causal modelling, combine both approaches, modelling the whole causal chain from external stimuli, via induced neural dynamics, to observed BOLD responses.     

Assessing the influence of different ROI selection strategies on functional connectivity analyses of fMRI data acquired during steady-state conditions

In blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI), assessing functional connectivity between and within brain networks from datasets acquired during steady-state conditions has become increasingly common. However, in contrast to connectivity analyses based on task-evoked signal changes, selecting the optimal spatial location of the regions of interest (ROIs) whose timecourses will be extracted and used in subsequent analyses is not straightforward. Moreover, it is also unknown how different choices of the precise anatomical locations within given brain regions influence the estimates of functional connectivity under steady-state conditions. The objective of the present study was to assess the variability in estimates of functional connectivity induced by different anatomical choices of ROI locations for a given brain network. We here targeted the default mode network (DMN) sampled during both resting-state and a continuous verbal 2-back working memory task to compare four different methods to extract ROIs in terms of ROI features (spatial overlap, spatial functional heterogeneity), signal features (signal distribution, mean, variance, correlation) as well as strength of functional connectivity as a function of condition. We show that, while different ROI selection methods produced quantitatively different results, all tested ROI selection methods agreed on the final conclusion that functional connectivity within the DMN decreased during the continuous working memory task compared to rest.