Mapping the MRI voxel volume in which thermal noise matches physiological noise--implications for fMRI

Physiological Contribution in Spontaneous Oscillations: An Approximate Quality-Assurance Index for Resting-State fMRI Signals

Resting-state fMRI (rs-fMRI) is receiving substantial attention for its sensitivity to functional abnormality in the brain networks of people with psychiatric and neurological disorders. However, because of the variety of rs-fMRI processing methods, the necessity of rs-fMRI quality assurance is increasing. Conventionally, the temporal signal-to-noise ratio (tSNR) is generally adopted for quality examination, but the tSNR does not guarantee reliable functional connectivity (FC) outcomes. Theoretically, intrinsic FC is supposed to reflect the spontaneous synchronization of neuronal basis, rather than that from thermal noise or non-neuronal physiological noise. Therefore, we proposed a new quality-assurance index for rs-fMRI to estimate the physiological contributions in spontaneous oscillations (PICSO). The PICSO index was designed as a voxel-wise measure for facilitating practical applications to all existing rs-fMRI data sets on the basis of two assumptions: Gaussian distributions in temporal fluctuations and ultra-slow changes of neural-based physiological fluctuations. To thoroughly validate the sensitivity of the proposed PICSO index to FC, we calibrated the preprocessing steps according to phantom data and verified the relationship between the PICSO and factors that are considered to affect FC in healthy participants (n = 12). Our results demonstrated that FC showed a significantly positive correlation with the PICSO. Moreover, for generating robust FC outcomes, directly acquiring data at a relatively large voxel size was more effective than performing smoothness on high-resolution data sets. In conclusion, compared with tSNR, the PICSO index is more sensitive to the resulting FC, providing a practical quality-assurance indicator for all existing rs-fMRI data sets.

Automatic EEG-assisted retrospective motion correction for fMRI (aE-REMCOR)

Head motions during functional magnetic resonance imaging (fMRI) impair fMRI data quality and introduce systematic artifacts that can affect interpretation of fMRI results. Electroencephalography (EEG) recordings performed simultaneously with fMRI provide high-temporal-resolution information about ongoing brain activity as well as head movements. Recently, an EEG-assisted retrospective motion correction (E-REMCOR) method was introduced. E-REMCOR utilizes EEG motion artifacts to correct the effects of head movements in simultaneously acquired fMRI data on a slice-by-slice basis. While E-REMCOR is an efficient motion correction approach, it involves an independent component analysis (ICA) of the EEG data and identification of motion-related ICs. Here we report an automated implementation of E-REMCOR, referred to as aE-REMCOR, which we developed to facilitate the application of E-REMCOR in large-scale EEG-fMRI studies. The aE-REMCOR algorithm, implemented in MATLAB, enables an automated preprocessing of the EEG data, an ICA decomposition, and, importantly, an automatic identification of motion-related ICs. aE-REMCOR has been used to perform retrospective motion correction for 305 fMRI datasets from 16 subjects, who participated in EEG-fMRI experiments conducted on a 3T MRI scanner. Performance of aE-REMCOR has been evaluated based on improvement in temporal signal-to-noise ratio (TSNR) of the fMRI data, as well as correction efficiency defined in terms of spike reduction in fMRI motion parameters. The results show that aE-REMCOR is capable of substantially reducing head motion artifacts in fMRI data. In particular, when there are significant rapid head movements during the scan, a large TSNR improvement and high correction efficiency can be achieved. Depending on a subject's motion, an average TSNR improvement over the brain upon the application of aE-REMCOR can be as high as 27%, with top ten percent of the TSNR improvement values exceeding 55%. The average correction efficiency over the 305 fMRI scans is 18% and the largest achieved efficiency is 71%. The utility of aE-REMCOR on the resting state fMRI connectivity of the default mode network is also examined. The motion-induced position-dependent error in the DMN connectivity analysis is shown to be reduced when aE-REMCOR is utilized. These results demonstrate that aE-REMCOR can be conveniently and efficiently used to improve fMRI motion correction in large clinical EEG-fMRI studies.

Using High Spatial Resolution to Improve BOLD fMRI Detection at 3T

For different functional magnetic resonance imaging experiments using blood oxygenation level-dependent (BOLD) contrast, the acquisition of T₂*-weighted scans at a high spatial resolution may be advantageous in terms of time-course signal-to-noise ratio and of BOLD sensitivity when the regions are prone to susceptibility artifacts. In this study, we explore this solution by examining how spatial resolution influences activations elicited when appetizing food pictures are viewed. Twenty subjects were imaged at 3 T with two different voxel volumes, 3.4 μl and 27 μl. Despite the diminution of brain coverage, we found that high-resolution acquisition led to a better detection of activations. Though known to suffer to different degrees from susceptibility artifacts, the activations detected by high spatial resolution were notably consistent with those reported in published activation likelihood estimation meta-analyses, corresponding to taste-responsive regions. Furthermore, these regions were found activated bilaterally, in contrast with previous findings. Both the reduction of partial volume effect, which improves BOLD contrast, and the mitigation of susceptibility artifact, which boosts the signal to noise ratio in certain regions, explained the better detection noted with high resolution. The present study provides further evidences that high spatial resolution is a valuable solution for human BOLD fMRI, especially for studying food-related stimuli.

Three-dimensional echo planar imaging with controlled aliasing: A sequence for high temporal resolution functional MRI

PURPOSE

In this work, we combine three-dimensional echo planar imaging (3D-EPI) with controlled aliasing to substantially increase temporal resolution in whole-brain functional MRI (fMRI) while minimizing geometry-factor (g-factor) losses.

THEORY AND METHODS

The study was performed on a 7 Tesla scanner equipped with a 32-channel receive coil. The proposed 3D-EPI-CAIPI sequence was evaluated for: (i) image quality, compared with a conventionally undersampled parallel imaging acquisition; (ii) temporal resolution, the ability to sample and remove physiological signal fluctuations from the fMRI signal of interest and (iii) the ability to distinguish small changes in hemodynamic responses in an event-related fMRI experiment.

RESULTS

Whole-brain fMRI data with a voxel size of 2 × 2 × 2 mm(3) and temporal resolution of 371 ms could be acquired with acceptable image quality. Ten-fold parallel imaging accelerated 3D-EPI-CAIPI data were shown to lower the maximum g-factor losses up to 62% with respect to a 10-fold accelerated 3D-EPI dataset. FMRI with 400 ms temporal resolution allowed the detection of time-to-peak variations in functional responses due to multisensory facilitation in temporal, occipital and frontal cortices.

CONCLUSION

3D-EPI-CAIPI allows increased temporal resolution that enables better characterization of physiological noise, thus improving sensitivity to signal changes of interest. Furthermore, subtle changes in hemodynamic response dynamics can be studied in shorter scan times by avoiding the need for jittering. Magn Reson Med 75:2350-2361, 2016. © 2015 Wiley Periodicals, Inc.

Incorporating relaxivities to more accurately reconstruct MR images

PURPOSE

To develop a mathematical model that incorporates the magnetic resonance relaxivities into the image reconstruction process in a single step.

MATERIALS AND METHODS

In magnetic resonance imaging, the complex-valued measurements of the acquired signal at each point in frequency space are expressed as a Fourier transformation of the proton spin density weighted by Fourier encoding anomalies: T2(⁎), T1, and a phase determined by magnetic field inhomogeneity (∆B) according to the MR signal equation. Such anomalies alter the expected symmetry and the signal strength of the k-space observations, resulting in images distorted by image warping, blurring, and loss in image intensity. Although T1 on tissue relaxation time provides valuable quantitative information on tissue characteristics, the T1 recovery term is typically neglected by assuming a long repetition time. In this study, the linear framework presented in the work of Rowe et al., 2007, and of Nencka et al., 2009 is extended to develop a Fourier reconstruction operation in terms of a real-valued isomorphism that incorporates the effects of T2(⁎), ∆B, and T1. This framework provides a way to precisely quantify the statistical properties of the corrected image-space data by offering a linear relationship between the observed frequency space measurements and reconstructed corrected image-space measurements. The model is illustrated both on theoretical data generated by considering T2(⁎), T1, and/or ∆B effects, and on experimentally acquired fMRI data by focusing on the incorporation of T1. A comparison is also made between the activation statistics computed from the reconstructed data with and without the incorporation of T1 effects.

RESULT

Accounting for T1 effects in image reconstruction is shown to recover image contrast that exists prior to T1 equilibrium. The incorporation of T1 is also shown to induce negligible correlation in reconstructed images and preserve functional activations.

CONCLUSION

With the use of the proposed method, the effects of T2(⁎) and ∆B can be corrected, and T1 can be incorporated into the time series image-space data during image reconstruction in a single step. Incorporation of T1 provides improved tissue segmentation over the course of time series and therefore can improve the precision of motion correction and image registration.

Does functional MRI detect activation in white matter? A review of emerging evidence, issues, and future directions

Functional magnetic resonance imaging (fMRI) is a non-invasive technique that allows for visualization of activated brain regions. Until recently, fMRI studies have focused on gray matter. There are two main reasons white matter fMRI remains controversial: (1) the blood oxygen level dependent (BOLD) fMRI signal depends on cerebral blood flow and volume, which are lower in white matter than gray matter and (2) fMRI signal has been associated with post-synaptic potentials (mainly localized in gray matter) as opposed to action potentials (the primary type of neural activity in white matter). Despite these observations, there is no direct evidence against measuring fMRI activation in white matter and reports of fMRI activation in white matter continue to increase. The questions underlying white matter fMRI activation are important. White matter fMRI activation has the potential to greatly expand the breadth of brain connectivity research, as well as improve the assessment and diagnosis of white matter and connectivity disorders. The current review provides an overview of the motivation to investigate white matter fMRI activation, as well as the published evidence of this phenomenon. We speculate on possible neurophysiologic bases of white matter fMRI signals, and discuss potential explanations for why reports of white matter fMRI activation are relatively scarce. We end with a discussion of future basic and clinical research directions in the study of white matter fMRI.

Voluntary enhancement of neural signatures of affiliative emotion using FMRI neurofeedback

In Ridley Scott's film "Blade Runner", empathy-detection devices are employed to measure affiliative emotions. Despite recent neurocomputational advances, it is unknown whether brain signatures of affiliative emotions, such as tenderness/affection, can be decoded and voluntarily modulated. Here, we employed multivariate voxel pattern analysis and real-time fMRI to address this question. We found that participants were able to use visual feedback based on decoded fMRI patterns as a neurofeedback signal to increase brain activation characteristic of tenderness/affection relative to pride, an equally complex control emotion. Such improvement was not observed in a control group performing the same fMRI task without neurofeedback. Furthermore, the neurofeedback-driven enhancement of tenderness/affection-related distributed patterns was associated with local fMRI responses in the septohypothalamic area and frontopolar cortex, regions previously implicated in affiliative emotion. This demonstrates that humans can voluntarily enhance brain signatures of tenderness/affection, unlocking new possibilities for promoting prosocial emotions and countering antisocial behavior.

You and your kin: Neural signatures of family-based group perception in the subgenual cortex

Attachment to one's kin as an in-group emerges from a fundamental human motivation and is vital for human survival. Despite important recent advances in the field of social neuroscience, the neural mechanisms underlying family-related in-group perception remain obscure. To examine the neural basis of perceiving family-related in-group boundaries in response to written kinship scenarios, we used functional magnetic resonance imaging in 27 healthy adults and obtained self-report ratings of family-related entitativity, which measures to what degree participants perceive their family as a coherent and distinct group in society. We expected that activity in the subgenual cingulate cortex and septo-hypothalamic region would track individual differences in entitativity. Perceiving one's family as a distinct and cohesive group (high entitativity) was associated with increased subgenual cortex response to kinship scenarios. The subgenual cingulate cortex may represent a key link between kin-related emotional attachment and group perception, providing a neurobiological basis for group belongingness.

Abnormal striatal BOLD responses to reward anticipation and reward delivery in ADHD

Altered reward processing has been proposed to contribute to the symptoms of attention deficit hyperactivity disorder (ADHD). The neurobiological mechanism underlying this alteration remains unclear. We hypothesize that the transfer of dopamine release from reward to reward-predicting cues, as normally observed in animal studies, may be deficient in ADHD. Functional magnetic resonance imaging (fMRI) was used to investigate striatal responses to reward-predicting cues and reward delivery in a classical conditioning paradigm. Data from 14 high-functioning and stimulant-naïve young adults with elevated lifetime symptoms of ADHD (8 males, 6 females) and 15 well-matched controls (8 males, 7 females) were included in the analyses. During reward anticipation, increased blood-oxygen-level-dependent (BOLD) responses in the right ventral and left dorsal striatum were observed in controls, but not in the ADHD group. The opposite pattern was observed in response to reward delivery; the ADHD group demonstrated significantly greater BOLD responses in the ventral striatum bilaterally and the left dorsal striatum relative to controls. In the ADHD group, the number of current hyperactivity/impulsivity symptoms was inversely related to ventral striatal responses during reward anticipation and positively associated with responses to reward. The BOLD response patterns observed in the striatum are consistent with impaired predictive dopamine signaling in ADHD, which may explain altered reward-contingent behaviors and symptoms of ADHD.

Enhanced phase regression with Savitzky-Golay filtering for high-resolution BOLD fMRI

Phase regression exploits the temporal evolution of phase in individual voxels to suppress blood oxygenation level dependent (BOLD) signal fluctuations caused by larger vessels and draining veins while preserving signal changes from microvascular effects. However, this process does not perform well when phase time series have low signal-to-noise ratios because of high levels of physiological noise. We demonstrate that Savitzky-Golay filters may be used to recover the underlying change in phase and completely restore the efficacy of phase regression. We do not make a priori assumptions regarding phase evolution and perform a data-driven exploration of parameter space to select the Savitzky-Golay filter parameters that minimize temporal variance in each voxel after phase regression. This approach is shown to work well on data acquired with single-shot and multi-shot pulse sequences, and should therefore be useful for both human and animal gradient-echo fMRI at high spatial resolutions at high fields. The ability to improve the spatial specificity of BOLD activation may be especially advantageous for clinical applications of fMRI that rely upon the accuracy of individual subject's activation maps to assist with presurgical planning and clinical decision-making. Enhanced phase regression with Savitzky-Golay filtering may also find other uses in analyses of resting state functional connectivity.

Functional dissociation of ventral frontal and dorsomedial default mode network components during resting state and emotional autobiographical recall

Humans spend a substantial share of their lives mind‐wandering. This spontaneous thinking activity usually comprises autobiographical recall, emotional, and self‐referential components. While neuroimaging studies have demonstrated that a specific brain “default mode network” (DMN) is consistently engaged by the “resting state” of the mind, the relative contribution of key cognitive components to DMN activity is still poorly understood. Here we used fMRI to investigate whether activity in neural components of the DMN can be differentially explained by active recall of relevant emotional autobiographical memories as compared with the resting state. Our study design combined emotional autobiographical memory, neutral memory and resting state conditions, separated by a serial subtraction control task. Shared patterns of activation in the DMN were observed in both emotional autobiographical and resting conditions, when compared with serial subtraction. Directly contrasting autobiographical and resting conditions demonstrated a striking dissociation within the DMN in that emotional autobiographical retrieval led to stronger activation of the dorsomedial core regions (medial prefrontal cortex, posterior cingulate cortex), whereas the resting state condition engaged a ventral frontal network (ventral striatum, subgenual and ventral anterior cingulate cortices) in addition to the IPL. Our results reveal an as yet unreported dissociation within the DMN. Whereas the dorsomedial component can be explained by emotional autobiographical memory, the ventral frontal one is predominantly associated with the resting state proper, possibly underlying fundamental motivational mechanisms engaged during spontaneous unconstrained ideation. Hum Brain Mapp 35:3302–3313, 2014. © 2013 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc..

The influence of spatial resolution and smoothing on the detectability of resting-state and task fMRI

Functional MRI blood oxygen level-dependent (BOLD) signal changes can be subtle, motivating the use of imaging parameters and processing strategies that maximize the temporal signal-to-noise ratio (tSNR) and thus the detection power of neuronal activity-induced fluctuations. Previous studies have shown that acquiring data at higher spatial resolutions results in greater percent BOLD signal changes, and furthermore that spatially smoothing higher resolution fMRI data improves tSNR beyond that of data originally acquired at a lower resolution. However, higher resolution images come at the cost of increased acquisition time, and the number of image volumes also influences detectability. The goal of our study is to determine how the detection power of neuronally induced BOLD fluctuations acquired at higher spatial resolutions and then spatially smoothed compares to data acquired at the lower resolutions with the same imaging duration. The number of time points acquired during a given amount of imaging time is a practical consideration given the limited ability of certain populations to lie still in the MRI scanner. We compare acquisitions at three different in-plane spatial resolutions (3.50×3.50mm(2), 2.33×2.33mm(2), 1.75×1.75mm(2)) in terms of their tSNR, contrast-to-noise ratio, and the power to detect both task-related activation and resting-state functional connectivity. The impact of SENSE acceleration, which speeds up acquisition time increasing the number of images collected, is also evaluated. Our results show that after spatially smoothing the data to the same intrinsic resolution, lower resolution acquisitions have a slightly higher detection power of task-activation in some, but not all, brain areas. There were no significant differences in functional connectivity as a function of resolution after smoothing. Similarly, the reduced tSNR of fMRI data acquired with a SENSE factor of 2 is offset by the greater number of images acquired, resulting in few significant differences in detection power of either functional activation or connectivity after spatial smoothing.

Whole-head rapid fMRI acquisition using echo-shifted magnetic resonance inverse imaging

The acquisition time of BOLD contrast functional MRI (fMRI) data with whole-brain coverage typically requires a sampling rate of one volume in 1-3s. Although the volumetric sampling time of a few seconds is adequate for measuring the sluggish hemodynamic response (HDR) to neuronal activation, faster sampling of fMRI might allow for monitoring of rapid physiological fluctuations and detection of subtle neuronal activation timing information embedded in BOLD signals. Previous studies utilizing a highly accelerated volumetric MR inverse imaging (InI) technique have provided a sampling rate of one volume per 100 ms with 5mm spatial resolution. Here, we propose a novel modification of this technique, the echo-shifted InI, which allows TE to be longer than TR, to measure BOLD fMRI at an even faster sampling rate of one volume per 25 ms with whole-brain coverage. Compared with conventional EPI, echo-shifted InI provided an 80-fold speedup with similar spatial resolution and less than 2-fold temporal SNR loss. The capability of echo-shifted InI to detect HDR timing differences was tested empirically. At the group level (n=6), echo-spaced InI was able to detect statistically significant HDR timing differences of as low as 50 ms in visual stimulus presentation. At the level of individual subjects, significant differences in HDR timing were detected for 400 ms stimulus-onset differences. Our results also show that the temporal resolution of 25 ms is necessary for maintaining the temporal detecting capability at this level. With the capabilities of being able to distinguish the timing differences in the millisecond scale, echo-shifted InI could be a useful fMRI tool for obtaining temporal information at a time scale closer to that of neuronal dynamics.

Reliable quantification of BOLD fMRI cerebrovascular reactivity despite poor breath-hold performance

Cerebrovascular reactivity (CVR) can be mapped using BOLD fMRI to provide a clinical insight into vascular health that can be used to diagnose cerebrovascular disease. Breath-holds are a readily accessible method for producing the required arterial CO2 increases but their implementation into clinical studies is limited by concerns that patients will demonstrate highly variable performance of breath-hold challenges. This study assesses the repeatability of CVR measurements despite poor task performance, to determine if and how robust results could be achieved with breath-holds in patients. Twelve healthy volunteers were scanned at 3 T. Six functional scans were acquired, each consisting of 6 breath-hold challenges (10, 15, or 20 s duration) interleaved with periods of paced breathing. These scans simulated the varying breath-hold consistency and ability levels that may occur in patient data. Uniform ramps, time-scaled ramps, and end-tidal CO2 data were used as regressors in a general linear model in order to measure CVR at the grey matter, regional, and voxelwise level. The intraclass correlation coefficient (ICC) quantified the repeatability of the CVR measurement for each breath-hold regressor type and scale of interest across the variable task performances. The ramp regressors did not fully account for variability in breath-hold performance and did not achieve acceptable repeatability (ICC<0.4) in several regions analysed. In contrast, the end-tidal CO2 regressors resulted in "excellent" repeatability (ICC=0.82) in the average grey matter data, and resulted in acceptable repeatability in all smaller regions tested (ICC>0.4). Further analysis of intra-subject CVR variability across the brain (ICCspatial and voxelwise correlation) supported the use of end-tidal CO2 data to extract robust whole-brain CVR maps, despite variability in breath-hold performance. We conclude that the incorporation of end-tidal CO2 monitoring into scanning enables robust, repeatable measurement of CVR that makes breath-hold challenges suitable for routine clinical practice.

Human Connectome Project informatics: quality control, database services, and data visualization

The Human Connectome Project (HCP) has developed protocols, standard operating and quality control procedures, and a suite of informatics tools to enable high throughput data collection, data sharing, automated data processing and analysis, and data mining and visualization. Quality control procedures include methods to maintain data collection consistency over time, to measure head motion, and to establish quantitative modality-specific overall quality assessments. Database services developed as customizations of the XNAT imaging informatics platform support both internal daily operations and open access data sharing. The Connectome Workbench visualization environment enables user interaction with HCP data and is increasingly integrated with the HCP's database services. Here we describe the current state of these procedures and tools and their application in the ongoing HCP study.

Correlated slow fluctuations in respiration, EEG, and BOLD fMRI

Low-frequency temporal fluctuations of physiological signals (<0.1 Hz), such as the respiration and cardiac pulse rate, occur naturally during rest and have been shown to be correlated with blood-oxygenation-level-dependent (BOLD) signal fluctuation. Such physiological signal modulations have been considered as sources of noise and their effects on BOLD signal are commonly removed in functional magnetic resonance imaging (fMRI) studies. However, possible neural correlates of the physiological fluctuations have not been considered nor examined in detail. In the present study we investigated this possibility by simultaneously acquiring electroencephalogram (EEG) with BOLD fMRI data, respiratory and cardiac waveforms in healthy human subjects at eyes-closed and eyes-open resting. We quantified the concurrent changes of the EEG power in the alpha frequency band, the respiration volume, and the cardiac pulse rate, then assessed the temporal correlations between alpha EEG power and physiological signal fluctuations. In addition, time-shifted time courses of alpha EEG power or physiological data were included as regressors to examine their correlations with the whole-brain BOLD fMRI signals. We observed a significant correlation between alpha EEG global field power and respiration, particularly at eyes-closed resting condition. Similar spatial patterns were observed between the correlation maps of BOLD with alpha EEG power and respiration, with negative correlations coinciding in the visual cortex, superior/middle temporal gyrus, inferior frontal gyrus, and inferior parietal lobule and positive correlations in the thalamus and caudate. Regressing out the physiological variations in the BOLD signal resulted in reduced correlation between BOLD and alpha EEG power. These results suggest a mutual link of neuronal origin between alpha EEG power, respiration, and BOLD signals.

Resting-state fMRI confounds and cleanup

The goal of resting-state functional magnetic resonance imaging (fMRI) is to investigate the brain's functional connections by using the temporal similarity between blood oxygenation level dependent (BOLD) signals in different regions of the brain "at rest" as an indicator of synchronous neural activity. Since this measure relies on the temporal correlation of fMRI signal changes between different parts of the brain, any non-neural activity-related process that affects the signals will influence the measure of functional connectivity, yielding spurious results. To understand the sources of these resting-state fMRI confounds, this article describes the origins of the BOLD signal in terms of MR physics and cerebral physiology. Potential confounds arising from motion, cardiac and respiratory cycles, arterial CO₂ concentration, blood pressure/cerebral autoregulation, and vasomotion are discussed. Two classes of techniques to remove confounds from resting-state BOLD time series are reviewed: 1) those utilising external recordings of physiology and 2) data-based cleanup methods that only use the resting-state fMRI data itself. Further methods that remove noise from functional connectivity measures at a group level are also discussed. For successful interpretation of resting-state fMRI comparisons and results, noise cleanup is an often over-looked but essential step in the analysis pipeline.

Sensitivity to white matter FMRI activation increases with field strength

Functional magnetic resonance imaging (fMRI) activation in white matter is controversial. Given that many of the studies that report fMRI activation in white matter used high field MRI systems, we investigated the field strength dependence of sensitivity to white matter fMRI activation. In addition, we evaluated the temporal signal to noise ratio (tSNR) of the different tissue types as a function of field strength. Data were acquired during a motor task (finger tapping) at 1.5 T and 4 T. Group and individual level activation results were considered in both the sensorimotor cortex and the posterior limb of the internal capsule. We found that sensitivity increases associated with field strength were greater for white matter than gray matter. The analysis of tSNR suggested that white matter might be less susceptible to increases in physiological noise related to increased field strength. We therefore conclude that high field MRI may be particularly advantageous for fMRI studies aimed at investigating activation in both gray and white matter.

Subject specific BOLD fMRI respiratory and cardiac response functions obtained from global signal

Subtle changes in either breathing pattern or cardiac pulse rate alter blood oxygen level dependent functional magnetic resonance imaging signal (BOLD fMRI). This is problematic because such fluctuations could possibly not be related to underlying neuronal activations of interest but instead the source of physiological noise. Several methods have been proposed to eliminate physiological noise in BOLD fMRI data. One such method is to derive a template based on average multi-subject data for respiratory response function (RRF) and cardiac response function (CRF) by simultaneously utilizing an external recording of cardiac and respiratory waveforms with the fMRI. Standard templates can then be used to model, map, and remove respiration and cardiac fluctuations from fMRI data. Utilizing these does not, however, account for intra-subject variations in physiological response. Thus, performing a more individualized approach for single subject physiological noise correction becomes more desirable, especially for clinical purposes. Here we propose a novel approach that employs subject-specific RRF and CRF response functions obtained from the whole brain or brain tissue-specific global signals (GS). Averaging multiple voxels in global signal computation ensures physiological noise dominance over thermal and system noise in even high-spatial-resolution fMRI data, making the GS suitable for deriving robust estimations of both RRF and CRF for individual subjects. Using these individualized response functions instead of standard templates based on multi-subject averages judiciously removes physiological noise from the data, assuming that there is minimal neuronal contribution in the derived individualized filters. Subject-specific physiological response functions obtained from the GS better maps individuals' physiological characteristics.

Modelling temporal stability of EPI time series using magnitude images acquired with multi-channel receiver coils

In 2001, Krueger and Glover introduced a model describing the temporal SNR (tSNR) of an EPI time series as a function of image SNR (SNR₀). This model has been used to study physiological noise in fMRI, to optimize fMRI acquisition parameters, and to estimate maximum attainable tSNR for a given set of MR image acquisition and processing parameters. In its current form, this noise model requires the accurate estimation of image SNR. For multi-channel receiver coils, this is not straightforward because it requires export and reconstruction of large amounts of k-space raw data and detailed, custom-made image reconstruction methods. Here we present a simple extension to the model that allows characterization of the temporal noise properties of EPI time series acquired with multi-channel receiver coils, and reconstructed with standard root-sum-of-squares combination, without the need for raw data or custom-made image reconstruction. The proposed extended model includes an additional parameter κ which reflects the impact of noise correlations between receiver channels on the data and scales an apparent image SNR (SNR′₀) measured directly from root-sum-of-squares reconstructed magnitude images so that κ = SNR′₀/SNR₀ (under the condition of SNR₀>50 and number of channels ≤32). Using Monte Carlo simulations we show that the extended model parameters can be estimated with high accuracy. The estimation of the parameter κ was validated using an independent measure of the actual SNR₀ for non-accelerated phantom data acquired at 3T with a 32-channel receiver coil. We also demonstrate that compared to the original model the extended model results in an improved fit to human task-free non-accelerated fMRI data acquired at 7T with a 24-channel receiver coil. In particular, the extended model improves the prediction of low to medium tSNR values and so can play an important role in the optimization of high-resolution fMRI experiments at lower SNR levels.

A neural signature of affiliative emotion in the human septohypothalamic area

Comparative studies have established that a number of structures within the rostromedial basal forebrain are critical for affiliative behaviors and social attachment. Lesion and neuroimaging studies concur with the importance of these regions for attachment and the experience of affiliation in humans as well. Yet it remains obscure whether the neural bases of affiliative experiences can be differentiated from the emotional valence with which they are inextricably associated at the experiential level. Here we show, using functional MRI, that kinship-related social scenarios evocative of affiliative emotion induce septal-preoptic-anterior hypothalamic activity that cannot be explained by positive or negative emotional valence alone. Our findings suggest that a phylogenetically conserved ensemble of basal forebrain structures, especially the septohypothalamic area, may play a key role in enabling human affiliative emotion. Our finding of a neural signature of human affiliative experience bears direct implications for the neurobiological mechanisms underpinning impaired affiliative experiences and behaviors in neuropsychiatric conditions.

Effects of image contrast on functional MRI image registration

Lack of tissue contrast and existing inhomogeneous bias fields from multi-channel coils have the potential to degrade the output of registration algorithms; and consequently degrade group analysis and any attempt to accurately localize brain function. Non-invasive ways to improve tissue contrast in fMRI images include the use of low flip angles (FAs) well below the Ernst angle and longer repetition times (TR). Techniques to correct intensity inhomogeneity are also available in most mainstream fMRI data analysis packages; but are not used as part of the pre-processing pipeline in many studies. In this work, we use a combination of real data and simulations to show that simple-to-implement acquisition/pre-processing techniques can significantly improve the outcome of both functional-to-functional and anatomical-to-functional image registrations. We also emphasize the need of tissue contrast on EPI images to be able to appropriately evaluate the quality of the alignment. In particular, we show that the use of low FAs (e.g., θ≤40°), when physiological noise considerations permit such an approach, significantly improves accuracy, consistency and stability of registration for data acquired at relatively short TRs (TR≤2s). Moreover, we also show that the application of bias correction techniques significantly improves alignment both for array-coil data (known to contain high intensity inhomogeneity) as well as birdcage-coil data. Finally, improvements in alignment derived from the use of the first infinite-TR volumes (ITVs) as targets for registration are also demonstrated. For the purpose of quantitatively evaluating the different scenarios, two novel metrics were developed: Mean Voxel Distance (MVD) to evaluate registration consistency, and Deviation of Mean Voxel Distance (dMVD) to evaluate registration stability across successive alignment attempts.

Layer-specific BOLD activation in awake monkey V1 revealed by ultra-high spatial resolution functional magnetic resonance imaging

The laminar structure of the cortex has previously been explored both in non-human primates and human subjects using high-resolution functional magnetic resonance imaging (fMRI). However, whether the spatial specificity of the blood-oxygenation-level-dependent (BOLD) fMRI is sufficiently high to reveal lamina specific organization in the cortex reliably is still unclear. In this study we demonstrate for the first time the detection of such layer-specific activation in awake monkeys at the spatial resolution of 200 × 200 × 1000 μm(3) in a vertical 4.7 T scanner. Results collected in trained monkeys are high in contrast-to-noise ratio and low in motion artifacts. Isolation of laminar activation was aided by choosing the optimal slice orientation and thickness using a novel pial vein pattern analysis derived from optical imaging. We found that the percent change of GE-BOLD signal is the highest at a depth corresponding to layer IV. Changes in the middle layers (layer IV) were 30% greater than changes in the top layers (layers I-III), and 32% greater than the bottom layers (layers V/VI). The laminar distribution of BOLD signal correlates well with neural activity reported in the literature. Our results suggest that the high intrinsic spatial resolution of GE-BOLD signal is sufficient for mapping sub-millimeter functional structures in awake monkeys. This degree of spatial specificity will be useful for mapping both laminar activations and columnar structures in the cerebral cortex.

Signal fluctuations in fMRI data acquired with 2D-EPI and 3D-EPI at 7 Tesla

Segmented three-dimensional echo planar imaging (3D-EPI) provides higher image signal-to-noise ratio (SNR) than standard single-shot two-dimensional echo planar imaging (2D-EPI), but is more sensitive to physiological noise. The aim of this study was to compare physiological noise removal efficiency in single-shot 2D-EPI and segmented 3D-EPI acquired at 7 Tesla. Two approaches were investigated based either on physiological regressors (PR) derived from cardiac and respiratory phases, or on principal component analysis (PCA) using additional resting-state data. Results show that, prior to physiological noise removal, 2D-EPI data had higher temporal SNR (tSNR), while spatial SNR was higher in 3D-EPI. Blood oxygen level dependent (BOLD) sensitivity was similar for both methods. The PR-based approach allowed characterization of relative contributions from different noise sources, confirming significant increases in physiological noise from 2D to 3D prior to correction. Both physiological noise removal approaches produced significant increases in tSNR and BOLD sensitivity, and these increases were larger for 3D-EPI, resulting in higher BOLD sensitivity in the 3D-EPI than in the 2D-EPI data. The PCA-based approach was the most effective correction method, yielding higher tSNR values for 3D-EPI than for 2D-EPI postcorrection.

Ultrafast inverse imaging techniques for fMRI

Inverse imaging (InI) supercharges the sampling rate of traditional functional MRI 10-100 fold at a cost of a moderate reduction in spatial resolution. The technique is inspired by similarities between multi-sensor magnetoencephalography (MEG) and highly parallel radio-frequency (RF) MRI detector arrays. Using presently available 32-channel head coils at 3T, InI can be sampled at 10 Hz and provides about 5-mm cortical spatial resolution with whole-brain coverage. Here we discuss the present applications of InI, as well as potential future challenges and opportunities in further improving its spatiotemporal resolution and sensitivity. InI may become a helpful tool for clinicians and neuroscientists for revealing the complex dynamics of brain functions during task-related and resting states.

EEG-assisted retrospective motion correction for fMRI: E-REMCOR

We propose a method for retrospective motion correction of fMRI data in simultaneous EEG-fMRI that employs the EEG array as a sensitive motion detector. EEG motion artifacts are used to generate motion regressors describing rotational head movements with millisecond temporal resolution. These regressors are utilized for slice-specific motion correction of unprocessed fMRI data. Performance of the method is demonstrated by correction of fMRI data from five patients with major depressive disorder, who exhibited head movements by 1-3mm during a resting EEG-fMRI run. The fMRI datasets, corrected using eight to ten EEG-based motion regressors, show significant improvements in temporal SNR (TSNR) of fMRI time series, particularly in the frontal brain regions and near the surface of the brain. The TSNR improvements are as high as 50% for large brain areas in single-subject analysis and as high as 25% when the results are averaged across the subjects. Simultaneous application of the EEG-based motion correction and physiological noise correction by means of RETROICOR leads to average TSNR enhancements as high as 35% for extended brain regions. These TSNR improvements are largely preserved after the subsequent fMRI volume registration and regression of fMRI motion parameters. The proposed EEG-assisted method of retrospective fMRI motion correction (referred to as E-REMCOR) can be applied to improve quality of fMRI data with severe motion artifacts and to reduce spurious correlations between the EEG and fMRI data caused by head movements. It does not require any specialized equipment beyond the standard EEG-fMRI instrumentation and can be applied retrospectively to any existing EEG-fMRI data set.

Multishot versus single-shot pulse sequences in very high field fMRI: a comparison using retinotopic mapping

High-resolution functional MRI is a leading application for very high field (7 Tesla) human MR imaging. Though higher field strengths promise improvements in signal-to-noise ratios (SNR) and BOLD contrast relative to fMRI at 3 Tesla, these benefits may be partially offset by accompanying increases in geometric distortion and other off-resonance effects. Such effects may be especially pronounced with the single-shot EPI pulse sequences typically used for fMRI at standard field strengths. As an alternative, one might consider multishot pulse sequences, which may lead to somewhat lower temporal SNR than standard EPI, but which are also often substantially less susceptible to off-resonance effects. Here we consider retinotopic mapping of human visual cortex as a practical test case by which to compare examples of these sequence types for high-resolution fMRI at 7 Tesla. We performed polar angle retinotopic mapping at each of 3 isotropic resolutions (2.0, 1.7, and 1.1 mm) using both accelerated single-shot 2D EPI and accelerated multishot 3D gradient-echo pulse sequences. We found that single-shot EPI indeed led to greater temporal SNR and contrast-to-noise ratios (CNR) than the multishot sequences. However, additional distortion correction in postprocessing was required in order to fully realize these advantages, particularly at higher resolutions. The retinotopic maps produced by both sequence types were qualitatively comparable, and showed equivalent test/retest reliability. Thus, when surface-based analyses are planned, or in other circumstances where geometric distortion is of particular concern, multishot pulse sequences could provide a viable alternative to single-shot EPI.

Spatiotemporal dynamics of the brain at rest--exploring EEG microstates as electrophysiological signatures of BOLD resting state networks

Neuroimaging research suggests that the resting cerebral physiology is characterized by complex patterns of neuronal activity in widely distributed functional networks. As studied using functional magnetic resonance imaging (fMRI) of the blood-oxygenation-level dependent (BOLD) signal, the resting brain activity is associated with slowly fluctuating hemodynamic signals (~10s). More recently, multimodal functional imaging studies involving simultaneous acquisition of BOLD-fMRI and electroencephalography (EEG) data have suggested that the relatively slow hemodynamic fluctuations of some resting state networks (RSNs) evinced in the BOLD data are related to much faster (~100 ms) transient brain states reflected in EEG signals, that are referred to as "microstates". To further elucidate the relationship between microstates and RSNs, we developed a fully data-driven approach that combines information from simultaneously recorded, high-density EEG and BOLD-fMRI data. Using independent component analysis (ICA) of the combined EEG and fMRI data, we identified thirteen microstates and ten RSNs that are organized independently in their temporal and spatial characteristics, respectively. We hypothesized that the intrinsic brain networks that are active at rest would be reflected in both the EEG data and the fMRI data. To test this hypothesis, the rapid fluctuations associated with each microstate were correlated with the BOLD-fMRI signal associated with each RSN. We found that each RSN was characterized further by a specific electrophysiological signature involving from one to a combination of several microstates. Moreover, by comparing the time course of EEG microstates to that of the whole-brain BOLD signal, on a multi-subject group level, we unraveled for the first time a set of microstate-associated networks that correspond to a range of previously described RSNs, including visual, sensorimotor, auditory, attention, frontal, visceromotor and default mode networks. These results extend our understanding of the electrophysiological signature of BOLD RSNs and demonstrate the intrinsic connection between the fast neuronal activity and slow hemodynamic fluctuations.

Variation of noise in multi-run functional MRI using generalized autocalibrating partially parallel acquisition (GRAPPA)

PURPOSE

To investigate the noise variation in multi-run functional MRI (fMRI) scans using generalized autocalibrating partially parallel acquisition (GRAPPA), with a focus on the cause of this variation.

MATERIALS AND METHODS

A phantom was continuously scanned for 10 runs using echo-planar imaging (EPI) combined with GRAPPA to simulate a multi-run fMRI exam. The variation of noise between runs was examined for different GRAPPA acceleration factors. The noise variation was also evaluated in a real fMRI experiment with human subjects at an acceleration factor of two. The cause of noise variation was explored by offline reconstruction using different GRAPPA weights and numerical simulation of GRAPPA reference scans.

RESULTS

It was found that the noise distribution in the image is stable within a run but may vary randomly from run to run. The variation of noise was also observed in fMRI experiments with human subjects. The variation can be significantly reduced if all the images from individual runs are reconstructed using the same reference scan data.

CONCLUSION

Both phantom experiments and human data showed that the noise pattern may change in different fMRI runs. The variation is mainly caused by the random noise in separate reference scans for GRAPPA in each run.

Phase stability in fMRI time series: effect of noise regression, off-resonance correction and spatial filtering techniques

Although the majority of fMRI studies exploit magnitude changes only, there is an increasing interest regarding the potential additive information conveyed by the phase signal. This integrated part of the complex number furnished by the MR scanners can also be used for exploring direct detection of neuronal activity and for thermography. Few studies have explicitly addressed the issue of the available signal stability in the context of phase time-series, and therefore we explored the spatial pattern of frequency specific phase fluctuations, and evaluated the effect of physiological noise components (heart beat and respiration) on the phase signal. Three categories of retrospective noise reduction techniques were explored and the temporal signal stability was evaluated in terms of a physiologic noise model, for seven fMRI measurement protocols in eight healthy subjects at 3T, for segmented CSF, gray and white matter voxels. We confirmed that for most processing methods, an efficient use of the phase information is hampered by the fact that noise from physiological and instrumental sources contributes significantly more to the phase than to the magnitude instability. Noise regression based on the phase evolution of the central k-space point, RETROICOR, or an orthonormalized combination of these were able to reduce their impact, but without bringing phase stability down to levels expected from the magnitude signal. Similar results were obtained after targeted removal of scan-to-scan variations in the bulk magnetic field by the dynamic off-resonance in k-space (DORK) method and by the temporal off-resonance alignment of single-echo time series technique (TOAST). We found that spatial high-pass filtering was necessary, and in vivo a Gaussian filter width of 20mm was sufficient to suppress physiological noise and bring the phase fluctuations to magnitude levels. Stronger filters brought the fluctuations down to levels dictated by thermal noise contributions, and for 62.5mm(3) voxels the phase stability was as low as 5 mrad (0.27°). In conditions of low SNR(o) and high temporal sampling rate (short TR); we achieved an upper bound for the phase instabilities at 0.0017 ppm, which is close to the dHb contribution to the GM/WM phase contrast.

High-field FMRI for human applications: an overview of spatial resolution and signal specificity

In the last decade, dozens of 7 Tesla scanners have been purchased or installed around the world, while 3 Tesla systems have become a standard. This increased interest in higher field strengths is driven by a demonstrated advantage of high fields for available signal-to-noise ratio (SNR) in the magnetic resonance signal. Functional imaging studies have additional advantages of increases in both the contrast and the spatial specificity of the susceptibility based BOLD signal. One use of this resultant increase in the contrast to noise ratio (CNR) for functional MRI studies at high field is increased image resolution. However, there are many factors to consider in predicting exactly what kind of resolution gains might be made at high fields, and what the opportunity costs might be. The first part of this article discusses both hardware and image quality considerations for higher resolution functional imaging. The second part draws distinctions between image resolution, spatial specificity, and functional specificity of the fMRI signals that can be acquired at high fields, suggesting practical limitations for attainable resolutions of fMRI experiments at a given field, given the current state of the art in imaging techniques. Finally, practical resolution limitations and pulse sequence options for studies in human subjects are considered.

Physiological noise reduction using volumetric functional magnetic resonance inverse imaging

Physiological noise arising from a variety of sources can significantly degrade the detection of task-related activity in BOLD-contrast fMRI experiments. If whole head spatial coverage is desired, effective suppression of oscillatory physiological noise from cardiac and respiratory fluctuations is quite difficult without external monitoring, since traditional EPI acquisition methods cannot sample the signal rapidly enough to satisfy the Nyquist sampling theorem, leading to temporal aliasing of noise. Using a combination of high speed magnetic resonance inverse imaging (InI) and digital filtering, we demonstrate that it is possible to suppress cardiac and respiratory noise without auxiliary monitoring, while achieving whole head spatial coverage and reasonable spatial resolution. Our systematic study of the effects of different moving average (MA) digital filters demonstrates that a MA filter with a 2 s window can effectively reduce the variance in the hemodynamic baseline signal, thereby achieving 57%-58% improvements in peak z-statistic values compared to unfiltered InI or spatially smoothed EPI data (FWHM = 8.6 mm). In conclusion, the high temporal sampling rates achievable with InI permit significant reductions in physiological noise using standard temporal filtering techniques that result in significant improvements in hemodynamic response estimation.

Cardiac-induced physiological noise in 3D gradient echo brain imaging: effect of k-space sampling scheme

The physiological noise in 3D image acquisition is shown to depend strongly on the sampling scheme. Five sampling schemes are considered: Linear, Centric, Segmented, Random and Tuned. Tuned acquisition means that data acquisition at k-space positions k and -k are separated with a specific time interval. We model physiological noise as a periodic temporal oscillation with arbitrary spatial amplitude in the physical object and develop a general framework to describe how this is rendered in the reconstructed image. Reconstructed noise can be decomposed in one component that is in phase with the signal (parallel) and one that is 90° out of phase (orthogonal). Only the former has a significant influence on the magnitude of the signal. The study focuses on fMRI using 3D EPI. Each k-space plane is acquired in a single shot in a time much shorter than the period of the physiological noise. The above mentioned sampling schemes are applied in the slow k-space direction and noise propagates almost exclusively in this direction. The problem then, is effectively one-dimensional. Numerical simulations and analytical expressions are presented. 3D noise measurements and 2D measurements with high temporal resolution are conducted. The measurements are performed under breath-hold to isolate the effect of cardiac-induced pulsatile motion. We compare the time-course stability of the sampling schemes and the extent to which noise propagates from a localized source into other parts of the imaging volume. Tuned and Linear acquisitions perform better than Centric, Segmented and Random.

Temporal SNR characteristics in segmented 3D-EPI at 7T

Three-dimensional segmented echo planar imaging (3D-EPI) is a promising approach for high-resolution functional magnetic resonance imaging, as it provides an increased signal-to-noise ratio (SNR) at similar temporal resolution to traditional multislice 2D-EPI readouts. Recently, the 3D-EPI technique has become more frequently used and it is important to better understand its implications for fMRI. In this study, the temporal SNR characteristics of 3D-EPI with varying numbers of segments are studied. It is shown that, in humans, the temporal variance increases with the number of segments used to form the EPI acquisition and that for segmented acquisitions, the maximum available temporal SNR is reduced compared to single shot acquisitions. This reduction with increased segmentation is not found in phantom data and thus likely due to physiological processes. When operating in the thermal noise dominated regime, fMRI experiments with a motor task revealed that the 3D variant outperforms the 2D-EPI in terms of temporal SNR and sensitivity to detect activated brain regions. Thus, the theoretical SNR advantage of a segmented 3D-EPI sequence for fMRI only exists in a low SNR situation. However, other advantages of 3D-EPI, such as the application of parallel imaging techniques in two dimensions and the low specific absorption rate requirements, may encourage the use of the 3D-EPI sequence for fMRI in situations with higher SNR.

The impact of physiological noise correction on fMRI at 7 T

Cognitive neuroimaging studies typically require fast whole brain image acquisition with maximal sensitivity to small BOLD signal changes. To increase the sensitivity, higher field strengths are often employed, since they provide an increased image signal-to-noise ratio (SNR). However, as image SNR increases, the relative contribution of physiological noise to the total time series noise will be greater compared to that from thermal noise. At 7 T, we studied how the physiological noise contribution can be best reduced for EPI time series acquired at three different spatial resolutions (1.1 mm × 1.1 mm × 1.8 mm, 2 mm × 2 mm × 2 mm and 3 mm × 3 mm × 3 mm). Applying optimal physiological noise correction methods improved temporal SNR (tSNR) and increased the numbers of significantly activated voxels in fMRI visual activation studies for all sets of acquisition parameters. The most dramatic results were achieved for the lowest spatial resolution, an acquisition parameter combination commonly used in cognitive neuroimaging which requires high functional sensitivity and temporal resolution (i.e. 3mm isotropic resolution and whole brain image repetition time of 2s). For this data, physiological noise models based on cardio-respiratory information improved tSNR by approximately 25% in the visual cortex and 35% sub-cortically. When the time series were additionally corrected for the residual effects of head motion after retrospective realignment, the tSNR was increased by around 58% in the visual cortex and 71% sub-cortically, exceeding tSNR ~140. In conclusion, optimal physiological noise correction at 7 T increases tSNR significantly, resulting in the highest tSNR per unit time published so far. This tSNR improvement translates into a significant increase in BOLD sensitivity, facilitating the study of even subtle BOLD responses.

Phasic and sustained fear in humans elicits distinct patterns of brain activity

Aversive events are typically more debilitating when they occur unpredictably than predictably. Studies in humans and animals indicate that predictable and unpredictable aversive events can induce phasic and sustained fear, respectively. Research in rodents suggests that anatomically related but distinct neural circuits may mediate phasic and sustained fear. We explored this issue in humans by examining threat predictability in three virtual reality contexts, one in which electric shocks were predictably signaled by a cue, a second in which shocks occurred unpredictably but never paired with a cue, and a third in which no shocks were delivered. Evidence of threat-induced phasic and sustained fear was presented using fear ratings and skin conductance. Utilizing recent advances in functional magnetic resonance imaging (fMRI), we were able to conduct whole-brain fMRI at relatively high spatial resolution and still have enough sensitivity to detect transient and sustained signal changes in the basal forebrain. We found that both predictable and unpredictable threat evoked transient activity in the dorsal amygdala, but that only unpredictable threat produced sustained activity in a forebrain region corresponding to the bed nucleus of the stria terminalis complex. Consistent with animal models hypothesizing a role for the cortex in generating sustained fear, sustained signal increases to unpredictable threat were also found in anterior insula and a frontoparietal cortical network associated with hypervigilance. In addition, unpredictable threat led to transient activity in the ventral amygdala-hippocampal area and pregenual anterior cingulate cortex, as well as transient activation and subsequent deactivation of subgenual anterior cingulate cortex, limbic structures that have been implicated in the regulation of emotional behavior and stress responses. In line with basic findings in rodents, these results provide evidence that phasic and sustained fear in humans may manifest similar signs of distress, but appear to be associated with different patterns of neural activity in the human basal forebrain.

Physiological noise effects on the flip angle selection in BOLD fMRI

This work addresses the choice of imaging flip angle in blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI). When noise of physiological origin becomes the dominant noise source in fMRI timeseries, it causes a nonlinear dependence of the temporal signal-to-noise ratio (TSNR) versus signal-to-noise ratio (SNR) that can be exploited to perform BOLD fMRI at angles well below the Ernst angle without any detrimental effect on our ability to detect sites of neuronal activation. We show, both experimentally and theoretically, that for situations where available SNR is high and physiological noise dominates over system/thermal noise, although TSNR still reaches it maximum for the Ernst angle, reduction of imaging flip angle well below this angle results in negligible loss in TSNR. Moreover, we provide a way to compute a suggested imaging flip angle, which constitutes a conservative estimate of the minimum flip angle that can be used under given experimental SNR and physiological noise levels. For our experimental conditions, this suggested angle equals 7.63° for the grey matter compartment, while the Ernst angle=77°. Finally, using data from eight subjects with a combined visual-motor task we show that imaging at angles as low as 9° introduces no significant differences in observed hemodynamic response time-course, contrast-to-noise ratio, voxel-wise effect size or statistical maps of activation as compared to imaging at 75° (an angle close to the Ernst angle). These results suggest that using low flip angles in BOLD fMRI experimentation to obtain benefits such as (1) reduction of RF power, (2) limitation of apparent T(1)-related inflow effects, (3) reduction of through-plane motion artifacts, (4) lower levels of physiological noise, and (5) improved tissue contrast is feasible when physiological noise dominates and SNR is high.

Human 7T MRI: First Clinical and Neuroscientific Applications

After many years of development and niche applications in a very limited number of laboratories, human 7T magnetic resonance imaging systems are becoming available to a number of clinical and neuroscientific researchers. The spectrum of available methods and their robustness is increasing and the first studies are underway to evaluate the potential applications and benefits for larger clinical studies or even clinical diagnosis. A number of technical and methodological challenges currently limit the application mainly to examinations of the brain. Some of the current possibilities of ultra-high field systems and examples of first applications in patient and research studies are demonstrated to give the reader an overview.

Functional connectivity in the rat at 11.7T: Impact of physiological noise in resting state fMRI

Resting state functional MRI (rs-fMRI) of the brain has the potential to elicit networks of functional connectivity and to reveal changes thereof in animal models of neurological disorders. In the present study, we investigate the contribution of physiological noise and its impact on assessment of functional connectivity in rs-fMRI of medetomidine sedated, spontaneously breathing rats at ultrahigh field of 11.7 Tesla. We employed gradient echo planar imaging (EPI) with repetition times of 3s and used simultaneous recordings of physiological parameters. A model of linear regression was applied to quantify the amount of BOLD fMRI signal fluctuations attributable to physiological sources. Our results indicate that physiological noise - mainly originating from the respiratory cycle -dominates the rs-fMRI time course in the form of spatially complex correlation patterns. As a consequence, these physiological fluctuations introduce severe artifacts into seed-based correlation maps and lead to misinterpretation of corresponding connectivity measures. We demonstrate that a scheme of motion correction and linear regression can significantly reduce physiological noise in the rs-fMRI time course, remove artifacts, and hence improve the reproducibility of functional connectivity assessment. In conclusion, physiological noise can severely compromise functional connectivity MRI (fcMRI) of the rodent at high fields and must be carefully considered in design and interpretation of future studies. Motion correction should be considered the primary strategy for reduction of apparent motion related to respiratory fluctuations. Combined with subsequent regression of physiological confounders, this strategy has proven successful in reducing physiological noise and related artifacts affecting functional connectivity analysis. The proposed new and rigorous protocol now opens the potential of fcMRI to elicit the role of brain connectivity in pathological processes without concerns of confounding contributions from physiological noise.

Selective functional integration between anterior temporal and distinct fronto-mesolimbic regions during guilt and indignation

It has been hypothesized that the experience of different moral sentiments such as guilt and indignation is underpinned by activation in temporal and fronto-mesolimbic regions and that functional integration between these regions is necessary for the differentiated experience of these moral sentiments. A recent fMRI study revealed that the right superior anterior temporal lobe (ATL) was activated irrespective of the context of moral feelings (guilt or indignation). This region has been associated with context-independent conceptual social knowledge which allows us to make fine-grained differentiations between qualities of social behaviours (e.g. "critical" and "faultfinding"). This knowledge is required to make emotional evaluations of social behaviour. In contrast to the context-independent activation of the ATL, there were context-dependent activations within different fronto-mesolimbic regions for guilt and indignation. However, it is unknown whether functional integration occurs between these regions and whether regional patterns of integration are distinctive for the experience of different moral sentiments. Here, we used fMRI and psychophysiological interaction analysis, an established measure of functional integration to investigate this issue. We found selective functional integration between the right superior ATL and a subgenual cingulate region during the experience of guilt and between the right superior ATL and the lateral orbitofrontal cortex for indignation. Our data provide the first evidence for functional integration of conceptual social knowledge representations in the right superior ATL with representations of different feeling contexts in fronto-mesolimbic regions. We speculate that this functional architecture allows for the conceptually differentiated experience of moral sentiments in healthy individuals.

Retinotopic mapping with spin echo BOLD at 7T

For blood oxygenation level-dependent (BOLD) functional MRI experiments, contrast-to-noise ratio (CNR) increases with increasing field strength for both gradient echo (GE) and spin echo (SE) BOLD techniques. However, susceptibility artifacts and nonuniform coil sensitivity profiles complicate large field-of-view fMRI experiments (e.g., experiments covering multiple visual areas instead of focusing on a single cortical region). Here, we use SE BOLD to acquire retinotopic mapping data in early visual areas, testing the feasibility of SE BOLD experiments spanning multiple cortical areas at 7T. We also use a recently developed method for normalizing signal intensity in T(1)-weighted anatomical images to enable automated segmentation of the cortical gray matter for scans acquired at 7T with either surface or volume coils. We find that the CNR of the 7T GE data (average single-voxel, single-scan stimulus coherence: 0.41) is almost twice that of the 3T GE BOLD data (average coherence: 0.25), with the CNR of the SE BOLD data (average coherence: 0.23) comparable to that of the 3T GE data. Repeated measurements in individual subjects find that maps acquired with 1.8-mm resolution at 3T and 7T with GE BOLD and at 7T with SE BOLD show no systematic differences in either the area or the boundary locations for V1, V2 and V3, demonstrating the feasibility of high-resolution SE BOLD experiments with good sensitivity throughout multiple visual areas.

Human fronto-tectal and fronto-striatal-tectal pathways activate differently during anti-saccades

Almost all cortical areas in the vertebrate brain take part in recurrent connections through the subcortical basal ganglia (BG) nuclei, through parallel inhibitory and excitatory loops. It has been suggested that these circuits can modulate our reactions to external events such that appropriate reactions are chosen from many available options, thereby imposing volitional control over behavior. The saccade system is an excellent model system to study cortico-BG interactions. In this study two possible pathways were investigated that might regulate automaticity of eye movements in the human brain; the cortico-tectal pathway, running directly between the frontal eye fields (FEF) and superior colliculus (SC) and the cortico-striatal pathway from the FEF to the SC involving the caudate nucleus (CN) in the BG. In an event-related functional magnetic resonance imaging (fMRI) paradigm participants made pro- and anti-saccades. A diffusion tensor imaging (DTI) scan was made for reconstruction of white matter tracts between the FEF, CN and SC. DTI fiber tracts were used to divide both the left and right FEF into two sub-areas, projecting to either ipsilateral SC or CN. For each of these FEF zones an event-related fMRI timecourse was extracted. In general activity in the FEF was larger for anti-saccades. This increase in activity was lateralized with respect to anti-saccade direction in FEF zones connected to the SC but not for zones only connected to the CN. These findings suggest that activity along the contralateral FEF-SC projection is responsible for directly generating anti-saccades, whereas the pathway through the BG might merely have a gating function withholding or allowing a pro-saccade.

Rapid 3D radial multi-echo functional magnetic resonance imaging

Functional magnetic resonance imaging with readouts at multiple echo times is useful for optimizing sensitivity across a range of tissue T2* values as well as for quantifying T2*. With single-shot acquisitions, both the minimum TE value and the number of TEs which it is possible to collect within a single TR are limited by the long echo-planar imaging readout duration (20-40 ms). In the present work, a multi-shot 3D radial acquisition which allows rapid whole-brain imaging at a range of echo times is proposed. The proposed 3D k-space coverage is implemented via a series of rotations of a single 2D interleaf. Data can be reconstructed at a variety of temporal resolutions from a single dataset, allowing for a flexible tradeoff between temporal resolution and BOLD contrast to noise ratio. It is demonstrated that whole-brain images at 5 echo times (TEs from 10 to 46 ms) can be acquired at a temporal rate as rapid as 400 ms/volume (3.75 mm isotropic resolution). Activation maps for a simultaneous motor/visual task consistent across multiple acceleration factors are obtained. Weighted combination of the echoes results in Z-scores that are significantly (p=0.016) higher than those resulting from any of the individual echo time images.

Flexible, phase-matched, linear receive arrays for high-field MRI in monkeys

High signal-to-noise ratios (SNR) are essential for high-resolution anatomical and functional MRI. Phased arrays are advantageous for this but have the drawback that they often have inflexible and bulky configurations. Particularly in experiments where functional MRI is combined with simultaneous electrophysiology, space constraints can be prohibitive. To this end we developed a highly flexible multiple receive element phased array for use on anesthetized monkeys. The elements are interchangeable and different sizes and combinations of coil elements can be used, for instance, combinations of single and overlapped elements. The preamplifiers including control electronics are detachable and can serve a variety of prefabricated and phase matched arrays of different configurations, allowing the elements to always be placed in close proximity to the area of interest. Optimizing performance of the individual elements ensured high SNR at the cortical surface as well as in deeper laying structures. Performance of a variety of arrangements of gapped linear arrays was evaluated at 4.7 and 7T in high-resolution anatomical and functional MRI.

Face-identity change activation outside the face system: "release from adaptation" may not always indicate neuronal selectivity

Face recognition is a complex cognitive process that requires distinguishable neuronal representations of individual faces. Previous functional magnetic resonance imaging (fMRI) studies using the "fMRI-adaptation" technique have suggested the existence of face-identity representations in face-selective regions, including the fusiform face area (FFA). Here, we present face-identity adaptation findings that are not well explained in terms of face-identity representations. We performed blood-oxygen level-dependent (BOLD) fMRI measurements, while participants viewed familiar faces that were shown repeatedly throughout the experiment. We found decreased activation for repeated faces in face-selective regions, as expected based on previous studies. However, we found similar effects in regions that are not face-selective, including the parahippocampal place area (PPA) and early visual cortex (EVC). These effects were present for exact-image (same view and lighting) as well as different-image (different view and/or lighting) repetition, but more widespread for exact-image repetition. Given the known functional properties of PPA and EVC, it appears unlikely that they contain domain-specific face-identity representations. Alternative interpretations include general attentional effects and carryover of activation from connected regions. These results remind us that fMRI stimulus-change effects can have a range of causes and do not provide conclusive evidence for a neuronal representation of the changed stimulus property.