Enhancement of temporal resolution and BOLD sensitivity in real-time fMRI using multi-slab echo-volumar imaging

Characterization of pulsations in the brain and cerebrospinal fluid using ultra-high field magnetic resonance imaging

The development of innovative non-invasive neuroimaging methods and biomarkers is critical for studying brain disease. Imaging of cerebrospinal fluid (CSF) pulsatility may inform the brain fluid dynamics involved in clearance of cerebral metabolic waste. In this work, we developed a methodology to characterize the frequency and spatial localization of whole brain CSF pulsations in humans. Using 7 Tesla (T) human magnetic resonance imaging (MRI) and ultrafast echo-planar imaging (EPI), in-vivo images were obtained to capture pulsations of the CSF signal. Physiological data were simultaneously collected and compared with the 7 T MR data. The primary components of signal pulsations were identified using spectral analysis, with the most evident frequency bands identified around 0.3, 1.2, and 2.4 Hz. These pulsations were mapped spatially and temporally onto the MR image domain and temporally onto the physiological measures of electrocardiogram and respiration. We identified peaks in CSF pulsations that were distinct from peaks in grey matter and white matter regions. This methodology may provide novel in vivo biomarkers of disrupted brain fluid dynamics.

Examining temporal features of BOLD-based cerebrovascular reactivity in clinical populations

Background

Conventional cerebrovascular reactivity (CVR) estimation has demonstrated that many brain diseases and/or conditions are associated with altered CVR. Despite the clinical potential of CVR, characterization of temporal features of a CVR challenge remains uncommon. This work is motivated by the need to develop CVR parameters that characterize individual temporal features of a CVR challenge.

Methods

Data were collected from 54 adults and recruited based on these criteria: (1) Alzheimer's disease diagnosis or subcortical Vascular Cognitive Impairment, (2) sleep apnea, and (3) subjective cognitive impairment concerns. We investigated signal changes in blood oxygenation level dependent (BOLD) contrast images with respect to hypercapnic and normocapnic CVR transition periods during a gas manipulation paradigm. We developed a model-free, non-parametric CVR metric after considering a range of responses through simulations to characterize BOLD signal changes that occur when transitioning from normocapnia to hypercapnia. The non-parametric CVR measure was used to examine regional differences across the insula, hippocampus, thalamus, and centrum semiovale. We also examined the BOLD signal transition from hypercapnia back to normocapnia.

Results

We found a linear association between isolated temporal features of successive CO₂ challenges. Our study concluded that the transition rate from hypercapnia to normocapnia was significantly associated with the second CVR response across all regions of interest (p < 0.001), and this association was highest in the hippocampus (R² = 0.57, p < 0.0125).

Conclusion

This study demonstrates that it is feasible to examine individual responses associated with normocapnic and hypercapnic transition periods of a BOLD-based CVR experiment. Studying these features can provide insight on between-subject differences in CVR.

Hi-Fi fMRI: High-resolution, fast-sampled and sub-second whole-brain functional MRI at 3T in humans

Functional magnetic resonance imaging (fMRI) is a methodological cornerstone of neuroscience. Most studies measure blood-oxygen-level-dependent (BOLD) signal using echo-planar imaging (EPI), Cartesian sampling, and image reconstruction with a one-to-one correspondence between the number of acquired volumes and reconstructed images. However, EPI schemes are subject to trade-offs between spatial and temporal resolutions. We overcome these limitations by measuring BOLD with a gradient recalled echo (GRE) with 3D radial-spiral phyllotaxis trajectory at a high sampling rate (28.24ms) on standard 3T field-strength. The framework enables the reconstruction of 3D signal time courses with whole-brain coverage at simultaneously higher spatial (1mm ³ ) and temporal (up to 250ms) resolutions, as compared to optimized EPI schemes. Additionally, artifacts are corrected before image reconstruction; the desired temporal resolution is chosen after scanning and without assumptions on the shape of the hemodynamic response. By showing activation in the calcarine sulcus of 20 participants performing an ON-OFF visual paradigm, we demonstrate the reliability of our method for cognitive neuroscience research.

Measuring brain beats: Cardiac-aligned fast functional magnetic resonance imaging signals

Blood and cerebrospinal fluid (CSF) pulse and flow throughout the brain, driven by the cardiac cycle. These fluid dynamics, which are essential to healthy brain function, are characterized by several noninvasive magnetic resonance imaging (MRI) methods. Recent developments in fast MRI, specifically simultaneous multislice acquisition methods, provide a new opportunity to rapidly and broadly assess cardiac-driven flow, including CSF spaces, surface vessels and parenchymal vessels. We use these techniques to assess blood and CSF flow dynamics in brief (3.5 min) scans on a conventional 3 T MRI scanner in five subjects. Cardiac pulses are measured with a photoplethysmography (PPG) on the index finger, along with functional MRI (fMRI) signals in the brain. We, retrospectively, align the fMRI signals to the heartbeat. Highly reliable cardiac-gated fMRI temporal signals are observed in CSF and blood on the timescale of one heartbeat (test-retest reliability within subjects R²  > 50%). In blood vessels, a local minimum is observed following systole. In CSF spaces, the ventricles and subarachnoid spaces have a local maximum following systole instead. Slower resting-state scans with slice timing, retrospectively, aligned to the cardiac pulse, reveal similar cardiac-gated responses. The cardiac-gated measurements estimate the amplitude and phase of fMRI pulsations in the CSF relative to those in the arteries, an estimate of the local intracranial impedance. Cardiac aligned fMRI signals can provide new insights about fluid dynamics or diagnostics for diseases where these dynamics are important.

Real-time and Recursive Estimators for Functional MRI Quality Assessment

Real-time quality assessment (rtQA) of functional magnetic resonance imaging (fMRI) based on blood oxygen level-dependent (BOLD) signal changes is critical for neuroimaging research and clinical applications. The losses of BOLD sensitivity because of different types of technical and physiological noise remain major sources of fMRI artifacts. Due to difficulty of subjective visual perception of image distortions during data acquisitions, a comprehensive automatic rtQA is needed. To facilitate rapid rtQA of fMRI data, we applied real-time and recursive quality assessment methods to whole-brain fMRI volumes, as well as time-series of target brain areas and resting-state networks. We estimated recursive temporal signal-to-noise ratio (rtSNR) and contrast-to-noise ratio (rtCNR), and real-time head motion parameters by a framewise rigid-body transformation (translations and rotations) using the conventional current to template volume registration. In addition, we derived real-time framewise (FD) and micro (MD) displacements based on head motion parameters and evaluated the temporal derivative of root mean squared variance over voxels (DVARS). For monitoring time-series of target regions and networks, we estimated the number of spikes and amount of filtered noise by means of a modified Kalman filter. Finally, we applied the incremental general linear modeling (GLM) to evaluate real-time contributions of nuisance regressors (linear trend and head motion). Proposed rtQA was demonstrated in real-time fMRI neurofeedback runs without and with excessive head motion and real-time simulations of neurofeedback and resting-state fMRI data. The rtQA was implemented as an extension of the open-source OpenNFT software written in Python, MATLAB and C++ for neurofeedback, task-based, and resting-state paradigms. We also developed a general Python library to unify real-time fMRI data processing and neurofeedback applications. Flexible estimation and visualization of rtQA facilitates efficient rtQA of fMRI data and helps the robustness of fMRI acquisitions by means of substantiating decisions about the necessity of the interruption and re-start of the experiment and increasing the confidence in neural estimates.

Ultrahigh Resolution fMRI at 7T Using Radial-Cartesian TURBINE Sampling

Purpose

We investigate the use of TURBINE, a 3D radial‐Cartesian acquisition scheme in which EPI planes are rotated about the phase‐encoding axis to acquire a cylindrical k‐space for high‐fidelity ultrahigh isotropic resolution fMRI at 7 Tesla with minimal distortion and blurring.

Methods

An improved, completely self‐navigated version of the TURBINE sampling scheme was designed for fMRI at 7 Telsa. To demonstrate the image quality and spatial specificity of the acquisition, thin‐slab visual and motor BOLD fMRI at 0.67 mm isotropic resolution (16 mm slab, TRvol = 2.32 s), and 0.8 × 0.8 × 2.0 mm (whole‐brain, TRvol = 2.4 s) data were acquired. To prioritize the high spatial fidelity, we employed a temporally regularized reconstruction to improve sensitivity without any spatial bias.

Results

TURBINE images provide high structural fidelity with almost no distortion, dropout, or T₂* blurring for the thin‐slab acquisitions compared to conventional 3D EPI owing to the radial sampling in‐plane and the short echo train used. This results in activation that can be localized to pre‐ and postcentral gyri in a motor task, for example, with excellent correspondence to brain structure measured by a T₁‐MPRAGE. The benefits of TURBINE (low distortion, dropout, blurring) are reduced for the whole‐brain acquisition due to the longer EPI train. We demonstrate robust BOLD activation at 0.67 mm isotropic resolution (thin‐slab) and also anisotropic 0.8 × 0.8 × 2.0 mm (whole‐brain) acquisitions.

Conclusion

TURBINE is a promising acquisition approach for high‐resolution, minimally distorted fMRI at 7 Tesla and could be particularly useful for fMRI in areas of high B₀ inhomogeneity.

Click here for author‐reader discussions

Monofractal analysis of functional magnetic resonance imaging: An introductory review

The following review will aid readers in providing an overview of scale-free dynamics and monofractal analysis, as well as its applications and potential in functional magnetic resonance imaging (fMRI) neuroscience and clinical research. Like natural phenomena such as the growth of a tree or crashing ocean waves, the brain expresses scale-invariant, or fractal, patterns in neural signals that can be measured. While neural phenomena may represent both monofractal and multifractal processes and can be quantified with many different interrelated parameters, this review will focus on monofractal analysis using the Hurst exponent (H). Monofractal analysis of fMRI data is an advanced analysis technique that measures the complexity of brain signaling by quantifying its degree of scale-invariance. As such, the H value of the blood oxygenation level-dependent (BOLD) signal specifies how the degree of correlation in the signal may mediate brain functions. This review presents a brief overview of the theory of fMRI monofractal analysis followed by notable findings in the field. Through highlighting the advantages and challenges of the technique, the article provides insight into how to best conduct fMRI fractal analysis and properly interpret the findings with physiological relevance. Furthermore, we identify the future directions necessary for its progression towards impactful functional neuroscience discoveries and widespread clinical use. Ultimately, this presenting review aims to build a foundation of knowledge among readers to facilitate greater understanding, discussion, and use of this unique yet powerful imaging analysis technique.

MRI magnetic compatible electrical neural interface: From materials to application

Neural electrical interfaces are important tools for local neural stimulation and recording, which potentially have wide application in the diagnosis and treatment of neural diseases, as well as in the transmission of neural activity for brain-computer interface (BCI) systems. At the same time, magnetic resonance imaging (MRI) is one of the effective and non-invasive techniques for recording whole-brain signals, providing details of brain structures and also activation pattern maps. Simultaneous recording of extracellular neural signals and MRI combines two expressions of the same neural activity and is believed to be of great importance for the understanding of brain function. However, this combination makes requests on the magnetic and electronic performance of neural interface devices. MRI-compatibility refers here to a technological approach to simultaneous MRI and electrode recording or stimulation without artifacts in imaging. Trade-offs between materials magnetic susceptibility selection and electrical function should be considered. Herein, prominent trends in selecting materials of suitable magnetic properties are analyzed and material design, function and application of neural interfaces are outlined together with the remaining challenge to fabricate MRI-compatible neural interface.

Three-dimensional reduced field-of-view imaging (3D-rFOVI)

PURPOSE

This study aimed at developing a 3D reduced field-of-view imaging (3D-rFOVI) technique using a 2D radiofrequency (RF) pulse, and demonstrating its ability to achieve isotropic high spatial resolution and reduced image distortion in echo planar imaging (EPI).

METHODS

The proposed 3D-rFOVI technique takes advantage of a 2D RF pulse to excite a slab along the conventional slice-selection direction (i.e., z-direction) while limiting the spatial extent along the phase-encoded direction (i.e., y-direction) within the slab. The slab is phase-encoded in both through-slab and in-slab phase-encoded directions. The 3D-rFOVI technique was implemented at 3T in gradient-echo and spin-echo EPI pulse sequences for functional MRI (fMRI) and diffusion-weighted imaging (DWI), respectively. 3D-rFOVI experiments were performed on a phantom and human brain to illustrate image distortion reduction, as well as isotropic high spatial resolution, in comparison with 3D full-FOV imaging.

RESULTS

In both the phantom and the human brain, image voxel dislocation was substantially reduced by 3D-rFOVI when compared with full-FOV imaging. In the fMRI experiment with visual stimulation, 3D isotropic spatial resolution of (2 × 2 × 2 mm³ ) was achieved with an adequate signal-to-noise ratio (81.5) and blood oxygen level-dependent (BOLD) contrast (2.5%). In the DWI experiment, diffusion-weighted brain images with an isotropic resolution of (1 × 1 × 1 mm³ ) was obtained without appreciable image distortion.

CONCLUSION

This study indicates that 3D-rFOVI is a viable approach to 3D neuroimaging over a zoomed region.

Slab boundary artifact correction in multislab imaging using convolutional-neural-network-enabled inversion for slab profile encoding

PURPOSE

This study aims to propose a novel algorithm for slab boundary artifact correction in both single-band multislab imaging and simultaneous multislab (SMSlab) imaging.

THEORY AND METHODS

In image domain, the formation of slab boundary artifacts can be regarded as modulating the artifact-free images using the slab profiles and introducing aliasing along the slice direction. Slab boundary artifact correction is the inverse problem of this process. An iterative algorithm based on convolutional neural networks (CNNs) is proposed to solve the problem, termed CNN-enabled inversion for slab profile encoding (CPEN). Diffusion-weighted SMSlab images and reference images without slab boundary artifacts were acquired in 7 healthy subjects for training. Images of 5 healthy subjects were acquired for testing, including single-band multislab and SMSlab images with 1.3-mm or 1-mm isotropic resolution. CNN-enabled inversion for slab profile encoding was compared with a previously reported method (i.e., nonlinear inversion for slab profile encoding [NPEN]).

RESULTS

CNN-enabled inversion for slab profile encoding reduces the slab boundary artifacts in both single-band multislab and SMSlab images. It also suppresses the slab boundary artifacts in the diffusion metric maps. Compared with NPEN, CPEN shows fewer residual artifacts in different acquisition protocols and more significant improvements in quantitative assessment, and it also accelerates the computation by more than 35 times.

CONCLUSION

CNN-enabled inversion for slab profile encoding can reduce the slab boundary artifacts in multislab acquisitions. It shows better slab boundary artifact correction capacity, higher robustness, and computation efficiency when compared with NPEN. It has the potential to improve the accuracy of multislab acquisitions in high-resolution DWI and functional MRI.

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

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

On the acquisition of the water signal during water suppression: High-speed MR spectroscopic imaging with water referencing and concurrent functional MRI

This study evaluated the utility of concurrent water signal acquisition as part of the water suppression in MR spectroscopic imaging (MRSI), to allow simultaneous water referencing for metabolite quantification, and to concurrently acquire functional MRI (fMRI) data. We integrated a spatial-spectral binomial water excitation RF pulse and a short spatial-spectral echo-planar readout into the water suppression module of 2D and 3D proton-echo-planar-spectroscopic-imaging (PEPSI) with a voxel size as small as 4 x 4 x 6 mm³ . Metabolite quantification in reference to tissue water was validated in healthy controls for different prelocalization methods (spin-echo, PRESS and semi-LASER) and the clinical feasibility of a 3-minute 3D semi-Laser PEPSI scan (TR/TE: 1250/32 ms) with water referencing in patients with brain tumors was demonstrated. Spectral quality, SNR, Cramer-Rao-lower-bounds and water suppression efficiency were comparable with conventional PEPSI. Metabolite concentration values in reference to tissue water, using custom LCModel-based spectral fitting with relaxation correction, were in the range of previous studies and independent of the prelocalization method used. Next, we added a phase-encoding undersampled echo-volumar imaging (EVI) module during water suppression to concurrently acquire metabolite maps with water referencing and fMRI data during task execution and resting state in healthy controls. Integration of multimodal signal acquisition prolongated minimum TR by less than 50 ms on average. Visual and motor activation in concurrent fMRI/MRSI (TR: 1250-1500 ms, voxel size: 4 x 4 x 6 mm³ ) was readily detectable in single-task blocks with percent signal change comparable with conventional fMRI. Resting-state connectivity in sensory and motor networks was detectable in 4 minutes. This hybrid water suppression approach for multimodal imaging has the potential to significantly reduce scan time and extend neuroscience research and clinical applications through concurrent quantitative MRSI and fMRI acquisitions.

Trading off spatio-temporal properties in 3D high-speed fMRI using interleaved stack-of-spirals trajectories

PURPOSE

Highly undersampled acquisitions have been proposed to push the limits of temporal resolution in functional MRI. This contribution is aimed at identifying parameter sets that let the user trade-off between ultra-high temporal resolution and spatial signal quality by varying the sampling densities. The proposed method maintains the synergies of a temporal resolution that enables direct filtering of physiological artifacts for highest statistical power, and 3D read-outs with optimal use of encoding capabilities of multi-coil arrays for efficient sampling and high signal-to-noise ratio (SNR).

METHODS

One- to four-shot interleaved spherical stack-of-spiral trajectories with repetition times from 96 to 352 ms at a nominal resolution of 3 mm using different sampling densities were compared for image quality and temporal SNR (tSNR). The one- and three-shot trajectories were employed in a resting state study for functional characterization.

RESULTS

Compared to a previously described single-shot trajectory, denser sampled trajectories of the same type are shown to be less prone to blurring and off-resonance vulnerability that appear in addition to the variable density artifacts of the point spread function. While the multi-shot trajectories lead to a decrease in tSNR efficiency, the high SNR due to the 3D read-out, combined with notable increases in image quality, leads to superior overall results of the three-shot interleaved stack of spirals. A resting state analysis of 15 subjects shows significantly improved functional sensitivity in areas of high off-resonance gradients.

CONCLUSION

Mild variable-density sampling leads to excellent tSNR behavior and no increased off-resonance vulnerability, and is suggested unless maximum temporal resolution is sought.

Vascular origins of low-frequency oscillations in the cerebrospinal fluid signal in resting-state fMRI: Interpretation using photoplethysmography

In vivo mapping of cerebrovascular oscillations in the 0.05–0.15 Hz remains difficult. Oscillations in the cerebrospinal fluid (CSF) represent a possible avenue for noninvasively tracking these oscillations using resting‐state functional MRI (rs‐fMRI), and have been used to correct for vascular oscillations in rs‐fMRI functional connectivity. However, the relationship between low‐frequency CSF and vascular oscillations remains unclear. In this study, we investigate this relationship using fast simultaneous rs‐fMRI and photoplethysmogram (PPG), examining the 0.1 Hz PPG signal, heart‐rate variability (HRV), pulse‐intensity ratio (PIR), and the second derivative of the PPG (SDPPG). The main findings of this study are: (a) signals in different CSF regions are not equivalent in their associations with vascular and tissue rs‐fMRI signals; (b) the PPG signal is maximally coherent with the arterial and CSF signals at the cardiac frequency, but coherent with brain tissue at ~0.2 Hz; (c) PIR is maximally coherent with the CSF signal near 0.03 Hz; and (d) PPG‐related vascular oscillations only contribute to ~15% of the CSF (and arterial) signal in rs‐fMRI. These findings caution against averaging all CSF regions when extracting physiological nuisance regressors in rs‐fMRI applications, and indicate the drivers of the CSF signal are more than simply cardiac. Our study is an initial attempt at the refinement and standardization of how the CSF signal in rs‐fMRI can be used and interpreted. It also paves the way for using rs‐fMRI in the CSF as a potential tool for tracking cerebrovascular health through, for instance, the potential relationship between PIR and the CSF signal.

Improving the sensitivity of spin-echo fMRI at 3T by highly accelerated acquisitions

PURPOSE

Spin-echo (SE) functional MRI (fMRI) can be highly advantageous compared to gradient-echo (GE) fMRI with respect to magnetic field-inhomogeneity artifacts. However, at 3T, the majority of blood oxygenation level-dependent (BOLD) fMRI experiments are performed using T 2 ∗ -weighted GE sequences because of their superior sensitivity compared to SE-fMRI. The presented SE implementation of a highly accelerated GE pulse sequence therefore aims to improve the sensitivity of SE-fMRI while profiting from a reduction of susceptibility-induced signal dropout.

METHODS

Spin-echo MR encephalography (SE-MREG) is compared with the more conventionally used spin-echo echo-planar imaging (SE-EPI) and spin-echo simultaneous multislice (SE-SMS) at 3T in terms of capability to detect neuronal activations and resting-state functional connectivity. For activation analysis, healthy subjects underwent consecutive SE-MREG (pulse repetition time [TR] = 0.25 seconds), SE-SMS (TR = 1.3 seconds), and SE-EPI (TR = 4.4 seconds) scans in pseudorandomized order applied to a visual block design paradigm for generation of t-statistics maps. For the investigation of functional connectivity, additional resting-state data were acquired for 5 minutes and a seed-based correlation analysis using Stanford's FIND (Functional Imaging in Neuropsychiatric Disorders) atlas was performed.

RESULTS

The increased sampling rate of SE-MREG relative to SE-SMS and SE-EPI improves the sensitivity to detect BOLD activation by 33% and 54%, respectively, and increases the capability to extract resting-state networks. Compared with a brain region that is not affected by magnetic field inhomogeneities, SE-MREG shows 2.5 times higher relative signal strength than GE-MREG in mesial temporal structures.

CONCLUSION

SE-MREG offers a viable possibility for whole-brain fMRI with consideration of brain regions that are affected by strong susceptibility-induced magnetic field gradients.

15 Years MR-encephalography

Objective

This review article gives an account of the development of the MR-encephalography (MREG) method, which started as a mere ‘Gedankenexperiment’ in 2005 and gradually developed into a method for ultrafast measurement of physiological activities in the brain. After going through different approaches covering k-space with radial, rosette, and concentric shell trajectories we have settled on a stack-of-spiral trajectory, which allows full brain coverage with (nominal) 3 mm isotropic resolution in 100 ms. The very high acceleration factor is facilitated by the near-isotropic k-space coverage, which allows high acceleration in all three spatial dimensions.

Methods

The methodological section covers the basic sequence design as well as recent advances in image reconstruction including the targeted reconstruction, which allows real-time feedback applications, and—most recently—the time-domain principal component reconstruction (tPCR), which applies a principal component analysis of the acquired time domain data as a sparsifying transformation to improve reconstruction speed as well as quality.

Applications

Although the BOLD-response is rather slow, the high speed acquisition of MREG allows separation of BOLD-effects from cardiac and breathing related pulsatility. The increased sensitivity enables direct detection of the dynamic variability of resting state networks as well as localization of single interictal events in epilepsy patients. A separate and highly intriguing application is aimed at the investigation of the glymphatic system by assessment of the spatiotemporal patterns of cardiac and breathing related pulsatility.

Discussion

MREG has been developed to push the speed limits of fMRI. Compared to multiband-EPI this allows considerably faster acquisition at the cost of reduced image quality and spatial resolution.

Real-Time Resting-State Functional Magnetic Resonance Imaging Using Averaged Sliding Windows with Partial Correlations and Regression of Confounding Signals

Background/Introduction: There is considerable interest in using real-time functional magnetic resonance imaging (fMRI) for monitoring functional connectivity dynamics. To date, the majority of real-time resting-state fMRI studies have examined a limited number of brain regions. This is, in part, due to the computational demands of traditional seed- and independent component analysis-based methods, in particular when using increasingly available high-speed fMRI methods. Methods: This study describes a computationally efficient, real-time, seed-based, resting-state fMRI analysis pipeline using moving averaged sliding-windows (ASW) with partial correlations and regression of motion parameters and signals from white matter and cerebrospinal fluid. Results: Analytical and numerical analyses of ASW correlation and sliding-window regression as a function of window width show selectable bandpass filter characteristics and effective suppression of artifactual correlations resulting from signal drifts and transients. The analysis pipeline is compatible with multislab echo-volumar imaging and simultaneous multislice echo-planar imaging with repetition times as short as 136 msec. High-speed, resting-state fMRI data in healthy controls demonstrate the effectiveness of this approach for minimizing artifactual correlations in white and gray matter, which was comparable to conventional regression across the entire scan. Integrating sliding-window averaging (width: W₁) within a second-level sliding-window (width: W₂) enabled monitoring of intra- and internetwork correlation dynamics of up to 12 resting-state networks with bandpass filter characteristics determined by the first-level sliding-window and temporal resolution W₁ + W₂. Conclusions: The computational performance and confound tolerance make this seed-based, resting-state fMRI approach suitable for real-time monitoring of data quality and resting-state connectivity dynamics in neuroscience and clinical research studies. Impact statement Using averaged sliding-windows for seed-based correlation and regression of confounding signals provides a powerful model-free approach to increase tolerance to artifactual signal transients in resting-state analysis. The algorithmic efficiency of this sliding-window approach enables real-time, seed-based, resting-state functional magnetic resonance imaging (fMRI) of multiple networks with computation of connectivity matrices and online monitoring of data quality. Integration of a second-level sliding-window enables mapping of resting-state connectivity dynamics. Sensitivity and tolerance to confounding signals compare favorably with conventional correlation and confound regression across the entire scan. This methodological advance has the potential to enhance the clinical utility of resting-state fMRI and facilitate neuroscience applications.

Quality and denoising in real-time functional magnetic resonance imaging neurofeedback: A methods review

Neurofeedback training using real-time functional magnetic resonance imaging (rtfMRI-NF) allows subjects voluntary control of localised and distributed brain activity. It has sparked increased interest as a promising non-invasive treatment option in neuropsychiatric and neurocognitive disorders, although its efficacy and clinical significance are yet to be determined. In this work, we present the first extensive review of acquisition, processing and quality control methods available to improve the quality of the neurofeedback signal. Furthermore, we investigate the state of denoising and quality control practices in 128 recently published rtfMRI-NF studies. We found: (a) that less than a third of the studies reported implementing standard real-time fMRI denoising steps, (b) significant room for improvement with regards to methods reporting and (c) the need for methodological studies quantifying and comparing the contribution of denoising steps to the neurofeedback signal quality. Advances in rtfMRI-NF research depend on reproducibility of methods and results. Notably, a systematic effort is needed to build up evidence that disentangles the various mechanisms influencing neurofeedback effects. To this end, we recommend that future rtfMRI-NF studies: (a) report implementation of a set of standard real-time fMRI denoising steps according to a proposed COBIDAS-style checklist (https://osf.io/kjwhf/), (b) ensure the quality of the neurofeedback signal by calculating and reporting community-informed quality metrics and applying offline control checks and (c) strive to adopt transparent principles in the form of methods and data sharing and support of open-source rtfMRI-NF software. Code and data for reproducibility, as well as an interactive environment to explore the study data, can be accessed at https://github.com/jsheunis/quality-and-denoising-in-rtfmri-nf.

Time-domain principal component reconstruction (tPCR): A more efficient and stable iterative reconstruction framework for non-Cartesian functional MRI

PURPOSE

To improve the reconstruction efficiency (i.e., computational load) and stability of iterative reconstruction for non-Cartesian fMRI when using high undersampling rates and/or in the presence of strong off-resonance effects.

THEORY AND METHODS

The magnetic resonance encephalography (MREG) sequence with 3D non-Cartesian trajectory and 0.1s repetition time (TR) was applied to acquire fMRI datasets. Different from a conventional time-point-by-time-point sequential reconstruction (SR), the proposed time-domain principal component reconstruction (tPCR) performs three steps: (1) decomposing the k-t-space fMRI datasets into time-domain principal component space using singular value decomposition, (2) reconstructing each principal component with redistributed computation power according to their weights, and (3) combining the reconstructed principal components back to image-t-space. The comparison of reconstruction accuracy was performed by simulation experiments and then verified in real fMRI data.

RESULTS

The simulation experiments showed that the proposed tPCR was able to significantly reduce reconstruction errors, and subsequent functional activation errors, relative to SR at identical computational cost. Alternatively, at fixed reconstruction accuracy, computation time was greatly reduced. The improved performance was particularly obvious for L1-norm nonlinear reconstructions relative to L2-norm linear reconstructions and robust to different regularization strength, undersampling rates, and off-resonance effects intensity. By examining activation maps, tPCR was also found to give similar improvements in real fMRI experiments.

CONCLUSION

The proposed proof-of-concept tPCR framework could improve (1) the reconstruction efficiency of iterative reconstruction, and (2) the reconstruction stability especially for nonlinear reconstructions. As a practical consideration, the improved reconstruction speed promotes the application of highly undersampled non-Cartesian fast fMRI.

Real-time presurgical resting-state fMRI in patients with brain tumors: Quality control and comparison with task-fMRI and intraoperative mapping

Resting-state functional magnetic resonance imaging (rsfMRI) is a promising task-free functional imaging approach, which may complement or replace task-based fMRI (tfMRI) in patients who have difficulties performing required tasks. However, rsfMRI is highly sensitive to head movement and physiological noise, and validation relative to tfMRI and intraoperative electrocortical mapping is still necessary. In this study, we investigate (a) the feasibility of real-time rsfMRI for presurgical mapping of eloquent networks with monitoring of data quality in patients with brain tumors and (b) rsfMRI localization of eloquent cortex compared with tfMRI and intraoperative electrocortical stimulation (ECS) in retrospective analysis. Five brain tumor patients were studied with rsfMRI and tfMRI on a clinical 3T scanner using MultiBand(8)-echo planar imaging (EPI) with repetition time: 400 ms. Moving-averaged sliding-window correlation analysis with regression of motion parameters and signals from white matter and cerebrospinal fluid was used to map sensorimotor and language resting-state networks. Data quality monitoring enabled rapid optimization of scan protocols, early identification of task noncompliance, and head movement-related false-positive connectivity to determine scan continuation or repetition. Sensorimotor and language resting-state networks were identifiable within 1 min of scan time. The Euclidean distance between ECS and rsfMRI connectivity and task-activation in motor cortex, Broca's, and Wernicke's areas was 5-10 mm, with the exception of discordant rsfMRI and ECS localization of Wernicke's area in one patient due to possible cortical reorganization and/or altered neurovascular coupling. This study demonstrates the potential of real-time high-speed rsfMRI for presurgical mapping of eloquent cortex with real-time data quality control, and clinically acceptable concordance of rsfMRI with tfMRI and ECS localization.

A 3D k-space Fourier encoding and reconstruction framework for simultaneous multi-slab acquisition

PURPOSE

To propose a novel 3D k-space Fourier encoding and reconstruction framework for simultaneous multi-slab (SMSlab) acquisition and demonstrate its efficacy in high-resolution imaging.

METHODS

First, it is illustrated in theory how the inter-slab gap interferes with the formation of the SMSlab 3D k-space. Then, joint RF and gradient encoding are applied to remove the inter-slab gap interference and form a SMSlab 3D k-space. In vivo experiments are performed to validate the proposed theory. Acceleration in the proposed SMSlab 3D k-space is also evaluated.

RESULTS

High-resolution (1.0 mm isotropic) images can be reconstructed using the proposed SMSlab 3D framework. Controlled aliasing in parallel imaging sampling and 2D GRAPPA reconstruction can also be applied in the SMSlab 3D k-space. Compared with conventional multi-slab acquisition, SMSlab exhibits better SNR maintainability (such as lower g-factors), especially at high acceleration factors.

CONCLUSION

It is demonstrated that the joint application of RF and gradient encoding enables SMSlab within a 3D Fourier encoding framework. Images with high isotropic resolution can be reconstructed, and further acceleration is also applicable. The proposed SMSlab 3D k-space can be valuable for both high-resolution and high-efficiency diffusion and functional MRI.

The potential of MR-Encephalography for BCI/Neurofeedback applications with high temporal resolution

Real-time functional magnetic resonance imaging (rt-fMRI) enables the update of various brain-activity measures during an ongoing experiment as soon as a new brain volume is acquired. However, the recorded Blood-oxygen-level dependent (BOLD) signal also contains physiological artifacts such as breathing and heartbeat, which potentially cause misleading false positive effects especially problematic in brain-computer interface (BCI) and neurofeedback (NF) setups. The low temporal resolution of echo planar imaging (EPI) sequences (which is in the range of seconds) prevents a proper separation of these artifacts from the BOLD signal. MR-Encephalography (MREG) has been shown to provide the high temporal resolution required to unalias and correct for physiological fluctuations and leads to increased specificity and sensitivity for mapping task-based activation and functional connectivity as well as for detecting dynamic changes in connectivity over time. By comparing a simultaneous multislice echo planar imaging (SMS-EPI) sequence and an MREG sequence using the same nominal spatial resolution in an offline analysis for three different experimental fMRI paradigms (perception of house and face stimuli, motor imagery, Stroop task), the potential of this novel technique for future BCI and NF applications was investigated. First, adapted general linear model pre-whitening which accounts for the high temporal resolution in MREG was implemented to calculate proper statistical results and be able to compare these with the SMS-EPI sequence. Furthermore, the respiration- and cardiac pulsation-related signals were successfully separated from the MREG signal using independent component analysis which were then included as regressors for a GLM analysis. Only the MREG sequence allowed to clearly separate cardiac pulsation and respiration components from the signal time course. It could be shown that these components highly correlate with the recorded respiration and cardiac pulsation signals using a respiratory belt and fingertip pulse plethysmograph. Temporal signal-to-noise ratios of SMS-EPI and MREG were comparable. Functional connectivity analysis using partial correlation showed a reduced standard error in MREG compared to SMS-EPI. Also, direct time course comparisons by down-sampling the MREG signal to the SMS-EPI temporal resolution showed lower variance in MREG. In general, we show that the higher temporal resolution is beneficial for fMRI time course modeling and this aspect can be exploited in offline application but also, is especially attractive, for real-time BCI and NF applications.

Olfactory fMRI: Implications of Stimulation Length and Repetition Time

Studying olfaction with functional magnetic resonance imaging (fMRI) poses various methodological challenges. This study aimed to investigate the effects of stimulation length and repetition time (TR) on the activation pattern of 4 olfactory brain regions: the anterior and the posterior piriform cortex, the orbitofrontal cortex, and the insula. Twenty-two healthy participants with normal olfaction were examined with fMRI, with 2 stimulation lengths (6 s and 15 s) and 2 TRs (0.901 s and 1.34 s). Data were analyzed using General Linear Model (GLM), Tensorial Independent Component Analysis (TICA), and by plotting the event-related time course of brain activation in the 4 olfactory regions of interest. The statistical analysis of the time courses revealed that short TR was associated with more pronounced signal increase and short stimulation was associated with shorter time to peak signal. Additionally, both long stimulation and short TR were associated with oscillatory time courses, whereas both short stimulation and short TR resulted in more typical time courses. GLM analysis showed that the combination of short stimulation and short TR could result in visually larger activation within these olfactory areas. TICA validated that the tested paradigm was spatially and temporally associated with a functionally connected network that included all 4 olfactory regions. In conclusion, the combination of short stimulation and short TR is associated with higher signal increase and shorter time to peak, making it more amenable to standard GLM-type analyses than long stimulation and long TR, and it should, thus, be preferable for olfactory fMRI.

Ultrafast fMRI of the rodent brain using simultaneous multi-slice EPI

Increasing spatial and temporal resolutions of functional MRI (fMRI) measurement has been shown to benefit the study of neural dynamics and functional interaction. However, acceleration of rodent brain fMRI using parallel and simultaneous multi-slice imaging techniques is hampered by the lack of high-density phased-array coils for the small brain. To overcome this limitation, we adapted phase-offset multiplanar and blipped-controlled aliasing echo planar imaging (EPI) to enable simultaneous multi-slice fMRI of the mouse brain using a single loop coil on a 9.4T scanner. Four slice bands of 0.3 × 0.3 × 0.5 mm³ resolution can be simultaneously acquired to cover the whole brain at a temporal resolution of 300 ms or the whole cerebrum in 150 ms. Instead of losing signal-to-noise ratio (SNR), both spatial and temporal SNR can be increased due to the increased k-space sampling compared to a standard single-band EPI. Task fMRI using a visual stimulation shows close to 80% increase of z-score and 4 times increase of activated area in the visual cortex using the multiband EPI due to the highly increased temporal samples. Resting-state fMRI shows reliable detection of bilateral connectivity by both single-band and multiband EPI, but no significant difference was found. Without the need of a dedicated hardware, we have demonstrated a practical method that can enable unparallelly fast whole-brain fMRI for preclinical studies. This technique can be used to increase sensitivity, distinguish transient response or acquire high spatiotemporal resolution fMRI.

Spectral entropy indicates electrophysiological and hemodynamic changes in drug-resistant epilepsy - A multimodal MREG study

Objective

Epilepsy causes measurable irregularity over a range of brain signal frequencies, as well as autonomic nervous system functions that modulate heart and respiratory rate variability. Imaging dynamic neuronal signals utilizing simultaneously acquired ultra-fast 10 Hz magnetic resonance encephalography (MREG), direct current electroencephalography (DC-EEG), and near-infrared spectroscopy (NIRS) can provide a more comprehensive picture of human brain function. Spectral entropy (SE) is a nonlinear method to summarize signal power irregularity over measured frequencies. SE was used as a joint measure to study whether spectral signal irregularity over a range of brain signal frequencies based on synchronous multimodal brain signals could provide new insights in the neural underpinnings of epileptiform activity.

Methods

Ten patients with focal drug-resistant epilepsy (DRE) and ten healthy controls (HC) were scanned with 10 Hz MREG sequence in combination with EEG, NIRS (measuring oxygenated, deoxygenated, and total hemoglobin: HbO, Hb, and HbT, respectively), and cardiorespiratory signals. After pre-processing, voxelwise SE_MREG was estimated from MREG data. Different neurophysiological and physiological subfrequency band signals were further estimated from MREG, DC-EEG, and NIRS: fullband (0–5 Hz, FB), near FB (0.08–5 Hz, NFB), brain pulsations in very-low (0.009–0.08 Hz, VLFP), respiratory (0.12–0.4 Hz, RFP), and cardiac (0.7–1.6 Hz, CFP) frequency bands. Global dynamic fluctuations in MREG and NIRS were analyzed in windows of 2 min with 50% overlap.

Results

Right thalamus, cingulate gyrus, inferior frontal gyrus, and frontal pole showed significantly higher SE_MREG in DRE patients compared to HC. In DRE patients, SE of cortical Hb was significantly reduced in FB (p = .045), NFB (p = .017), and CFP (p = .038), while both HbO and HbT were significantly reduced in RFP (p = .038, p = .045, respectively). Dynamic SE of HbT was reduced in DRE patients in RFP during minutes 2 to 6. Fitting to the frontal MREG and NIRS results, DRE patients showed a significant increase in SE_EEG in FB in fronto-central and parieto-occipital regions, in VLFP in parieto-central region, accompanied with a significant decrease in RFP in frontal pole and parietal and occipital (O2, Oz) regions.

Conclusion

This is the first study to show altered spectral entropy from synchronous MREG, EEG, and NIRS in DRE patients. Higher SE_MREG in DRE patients in anterior cingulate gyrus together with SE_EEG and SE_NIRS results in 0.12–0.4 Hz can be linked to altered parasympathetic function and respiratory pulsations in the brain. Higher SE_MREG in thalamus in DRE patients is connected to disturbances in anatomical and functional connections in epilepsy. Findings suggest that spectral irregularity of both electrophysiological and hemodynamic signals are altered in specific way depending on the physiological frequency range.

•

Simultaneous imaging methods indicate spectral irregularity in neurovascular and electrophysiological brain pulsations in DRE.

•

Altered spectral entropy in EEG, NIRS and BOLD indicate dysfunctional brain pulsations in respiratory frequency in epilepsy.

•

Spectral irregularity (0-5 Hz) of BOLD in right thalamus supports previous structural and functional findings in epilepsy.

On the analysis of rapidly sampled fMRI data

Recent advances in parallel imaging and simultaneous multi-slice techniques have permitted whole-brain fMRI acquisitions at sub-second sampling intervals, without significantly sacrificing the spatial coverage and resolution. Apart from probing brain function at finer temporal scales, faster sampling rates may potentially lead to enhanced functional sensitivity, owing possibly to both cleaner neural representations (due to less aliased physiological noise) and additional statistical benefits (due to more degrees of freedom for a fixed scan duration). Accompanying these intriguing aspects of fast acquisitions, however, confusion has also arisen regarding (1) how to preprocess/analyze these fast fMRI data, and (2) what exactly is the extent of benefits with fast acquisitions, i.e., how fast is fast enough for a specific research aim? The first question is motivated by the altered spectral distribution and noise characteristics at short sampling intervals, while the second question seeks to reconcile the complicated trade-offs between the functional contrast-to-noise ratio and the effective degrees of freedom. Although there have been recent efforts to empirically approach different aspects of these two questions, in this work we discuss, from a theoretical perspective accompanied by some illustrative, proof-of-concept experimental in vivo human fMRI data, a few considerations that are rarely mentioned, yet are important for both preprocessing and optimizing statistical inferences for studies that employ acquisitions with sub-second sampling intervals. Several summary recommendations include concerns regarding advisability of relying on low-pass filtering to de-noise physiological contributions, employment of statistical models with sufficient complexity to account for the substantially increased serial correlation, and cautions regarding using rapid sampling to enhance functional sensitivity given that different analysis models may associate with distinct trade-offs between contrast-to-noise ratios and the effective degrees of freedom. As an example, we demonstrate that as TR shortens, the intrinsic differences in how noise is accommodated in general linear models and Pearson correlation analyses (assuming Gaussian distributed stochastic signals and noise) can result in quite different outcomes, either gaining or losing statistical power.

No time for drifting: Comparing performance and applicability of signal detrending algorithms for real-time fMRI

As a consequence of recent technological advances in the field of functional magnetic resonance imaging (fMRI), results can now be made available in real-time. This allows for novel applications such as online quality assurance of the acquisition, intra-operative fMRI, brain-computer-interfaces, and neurofeedback. To that aim, signal processing algorithms for real-time fMRI must reliably correct signal contaminations due to physiological noise, head motion, and scanner drift. The aim of this study was to compare performance of the commonly used online detrending algorithms exponential moving average (EMA), incremental general linear model (iGLM) and sliding window iGLM (iGLM^window). For comparison, we also included offline detrending algorithms (i.e., MATLAB's and SPM8's native detrending functions). Additionally, we optimized the EMA control parameter, by assessing the algorithm's performance on a simulated data set with an exhaustive set of realistic experimental design parameters. First, we optimized the free parameters of the online and offline detrending algorithms. Next, using simulated data, we systematically compared the performance of the algorithms with respect to varying levels of Gaussian and colored noise, linear and non-linear drifts, spikes, and step function artifacts. Additionally, using in vivo data from an actual rt-fMRI experiment, we validated our results in a post hoc offline comparison of the different detrending algorithms. Quantitative measures show that all algorithms perform well, even though they are differently affected by the different artifact types. The iGLM approach outperforms the other online algorithms and achieves online detrending performance that is as good as that of offline procedures. These results may guide developers and users of real-time fMRI analyses tools to best account for the problem of signal drifts in real-time fMRI.

•

fMRI time series are almost always affected by signal drifts.

•

Signal drifts can be reduced in real-time using different correction methods.

•

Robustness of these methods was tested in the presence of different artifacts.

•

Incremental GLM (iGLM, iGLM^window) were found optimal in most cases.

An approach to directly link ICA and seed-based functional connectivity: Application to schizophrenia

Independent component analysis (ICA) and seed-based analyses are widely used techniques for studying intrinsic neuronal activity in task-based or resting scans. In this work, we show there is a direct link between the two, and show that there are some important differences between the two approaches in terms of what information they capture. We developed an enhanced connectivity-matrix independent component analysis (cmICA) for calculating whole brain voxel maps of functional connectivity, which reduces the computational complexity of voxel-based connectivity analysis on performing many temporal correlations. We also show there is a mathematical equivalency between parcellations on voxel-to-voxel functional connectivity and simplified cmICA. Next, we used this cost-efficient data-driven method to examine the resting state fMRI connectivity in schizophrenia patients (SZ) and healthy controls (HC) on a whole brain scale and further quantified the relationship between brain functional connectivity and cognitive performances measured by the Measurement and Treatment Research to Improve Cognition in Schizophrenia (MATRICS) battery. Current results suggest that SZ exhibit a wide-range abnormality, primarily a decrease, in functional connectivity both between networks and within different network hubs. Specific functional connectivity decreases were associated with MATRICS performance deficits. In addition, we found that resting state functional connectivity decreases was extensively associated with aging regardless of groups. In contrast, there was no relationship between positive and negative symptoms in the patients and functional connectivity. In sum, we have developed a novel mathematical relationship between ICA and seed-based connectivity that reduces computational complexity, which has broad applicability, and showed a specific application of this approach to characterize connectivity changes associated with cognitive scores in SZ.

Targeted partial reconstruction for real-time fMRI with arbitrary trajectories

PURPOSE

A partial image reconstruction formalism is introduced for the targeted extraction of real-time feedback from arbitrary trajectories when full image reconstruction in real time is computationally too demanding.

METHODS

Explicit calculation and storage of linear combinations of lines of the reconstruction matrix by an incomplete basis change in spatial coordinates lead to translation of the expensive full reconstruction from a frame-wise application to a region of interest (ROI)-wise application. This step is independent from signal data and can be executed before the experiment. Subsequently, the results of the sum over fully reconstructed voxels can be evaluated directly. Data from a high-speed fMRI acquisition was used to investigate the targeted partial reconstruction of a functional ROI atlas, incorporating an intravolume dephasing correction. The same data and ROIs were used for a comparison of the time series obtained with those obtained from already existing methods for compartment-wise reconstruction. To examine real-time feasibility, the reconstruction was implemented and tested for online reconstruction performance.

RESULTS

The reconstruction yields results that are virtually identical to the standard reconstruction (i.e., the magnitude sums over the ROIs), with negligible discrepancies even after termination of the conjugate gradient algorithm at a feasible number of iterations. Notably, more discrepancies arise with existing compartment-wise reconstructions. The online real-time implementation evaluated 1 ROI within 2.8 ms in the case of a highly parallel 3D whole brain acquisition.

CONCLUSION

The high reconstruction fidelity and speed are satisfying for the exemplary application of real-time functional feedback using a highly parallel 3D whole brain acquisition.

Sparse Estimation of Resting-State Effective Connectivity From fMRI Cross-Spectra

In functional magnetic resonance imaging (fMRI), functional connectivity is conventionally characterized by correlations between fMRI time series, which are intrinsically undirected measures of connectivity. Yet, some information about the directionality of network connections can nevertheless be extracted from the matrix of pairwise temporal correlations between all considered time series, when expressed in the frequency-domain as a cross-spectral density matrix. Using a sparsity prior, it then becomes possible to determine a unique directed network topology that best explains the observed undirected correlations, without having to rely on temporal precedence relationships that may not be valid in fMRI. Applying this method on simulated data with 100 nodes yielded excellent retrieval of the underlying directed networks under a wide variety of conditions. Importantly, the method did not depend on temporal precedence to establish directionality, thus reducing susceptibility to hemodynamic variability. The computational efficiency of the algorithm was sufficient to enable whole-brain estimations, thus circumventing the problem of missing nodes that otherwise occurs in partial-brain analyses. Applying the method to real resting-state fMRI data acquired with a high temporal resolution, the inferred networks showed good consistency with structural connectivity obtained from diffusion tractography in the same subjects. Interestingly, this agreement could also be seen when considering high-frequency rather than low-frequency connectivity (average correlation: r = 0.26 for f < 0.3 Hz, r = 0.43 for 0.3 < f < 5 Hz). Moreover, this concordance was significantly better (p < 0.05) than for networks obtained with conventional functional connectivity based on correlations (average correlation r = 0.18). The presented methodology thus appears to be well-suited for fMRI, particularly given its lack of explicit dependence on temporal lag structure, and is readily applicable to whole-brain effective connectivity estimation.

Fast Spatial Resolution Analysis of Quadratic Penalized Least-Squares Image Reconstruction With Separate Real and Imaginary Roughness Penalty: Application to fMRI

Penalized least-squares iterative image reconstruction algorithms used for spatial resolution-limited imaging, such as functional magnetic resonance imaging (fMRI), commonly use a quadratic roughness penalty to regularize the reconstructed images. When used for complex-valued images, the conventional roughness penalty regularizes the real and imaginary parts equally. However, these imaging methods sometimes benefit from separate penalties for each part. The spatial smoothness from the roughness penalty on the reconstructed image is dictated by the regularization parameter(s). One method to set the parameter to a desired smoothness level is to evaluate the full width at half maximum of the reconstruction method's local impulse response. Previous work has shown that when using the conventional quadratic roughness penalty, one can approximate the local impulse response using an FFT-based calculation. However, that acceleration method cannot be applied directly for separate real and imaginary regularization. This paper proposes a fast and stable calculation for this case that also uses FFT-based calculations to approximate the local impulse responses of the real and imaginary parts. This approach is demonstrated with a quadratic image reconstruction of fMRI data that uses separate roughness penalties for the real and imaginary parts.

Simultaneous multi-slice inverse imaging of the human brain

Ultrafast functional magnetic resonance imaging (fMRI) can measure blood oxygen level dependent (BOLD) signals with high sensitivity and specificity. Here we propose a novel method: simultaneous multi-slice inverse imaging (SMS-InI) — a combination of simultaneous multi-slice excitation, simultaneous echo refocusing (SER), blipped controlled aliasing in parallel imaging echo-planar imaging (EPI), and regularized image reconstruction. Using a 32-channel head coil array on a 3 T scanner, SMS-InI achieves nominal isotropic 5-mm spatial resolution and 10 Hz sampling rate at the whole-brain level. Compared with traditional inverse imaging, we found that SMS-InI has higher spatial resolution with lower signal leakage and higher time-domain signal-to-noise ratio with the optimized regularization parameter in the reconstruction. SMS-InI achieved higher effective resolution and higher detection power in detecting visual cortex activity than InI. SMS-InI also detected subcortical fMRI signals with the similar sensitivity and localization accuracy like EPI. The spatiotemporal resolution of SMS-InI was used to reveal that presenting visual stimuli with 0.2 s latency between left and right visual hemifield led to 0.2 s relative hemodynamic response latency between the left and right visual cortices. Together, these results indicate that SMS-InI is a useful tool in measuring cortical and subcortical hemodynamic responses with high spatiotemporal resolution.

Novel imaging techniques to study the functional organization of the human brain

Despite more than a century of investigation into the cortical organization of motor function, the existence of motor somatotopy is still debated. We review functional magnetic resonance imaging (fMRI) studies examining motor somatotopy in the cerebral cortex. In spite of a substantial overlap of representations corresponding to different body parts, especially in non-primary motor cortices, geographic approaches are capable of revealing somatotopic ordering. From the iconic homunculus in the contralateral primary cortex to the subtleties of ipsilateral somatotopy and its relations with lateralization, we outline potential reasons for the lack of segregation between motor representations. Among these are the difficulties in distinguishing activity that arises from multiple muscular effectors, the need for flexible motor control and coordination of complex movements through functional integration and artefacts in fMRI. Methodological advances with regard to the optimization of experimental design and fMRI acquisition protocols as well as improvements in spatial registration of images and indices aiming at the quantification of the degree of segregation between different functional representations are inspected. Additionally, we give some hints as to how the functional organization of motor function might be related to various anatomical landmarks in brain morphometry.

Improving temporal resolution in fMRI using a 3D spiral acquisition and low rank plus sparse (L+S) reconstruction

Rapid whole-brain dynamic Magnetic Resonance Imaging (MRI) is of particular interest in Blood Oxygen Level Dependent (BOLD) functional MRI (fMRI). Faster acquisitions with higher temporal sampling of the BOLD time-course provide several advantages including increased sensitivity in detecting functional activation, the possibility of filtering out physiological noise for improving temporal SNR, and freezing out head motion. Generally, faster acquisitions require undersampling of the data which results in aliasing artifacts in the object domain. A recently developed low-rank (L) plus sparse (S) matrix decomposition model (L+S) is one of the methods that has been introduced to reconstruct images from undersampled dynamic MRI data. The L+S approach assumes that the dynamic MRI data, represented as a space-time matrix M, is a linear superposition of L and S components, where L represents highly spatially and temporally correlated elements, such as the image background, while S captures dynamic information that is sparse in an appropriate transform domain. This suggests that L+S might be suited for undersampled task or slow event-related fMRI acquisitions because the periodic nature of the BOLD signal is sparse in the temporal Fourier transform domain and slowly varying low-rank brain background signals, such as physiological noise and drift, will be predominantly low-rank. In this work, as a proof of concept, we exploit the L+S method for accelerating block-design fMRI using a 3D stack of spirals (SoS) acquisition where undersampling is performed in the k_z-t domain. We examined the feasibility of the L+S method to accurately separate temporally correlated brain background information in the L component while capturing periodic BOLD signals in the S component. We present results acquired in control human volunteers at 3T for both retrospective and prospectively acquired fMRI data for a visual activation block-design task. We show that a SoS fMRI acquisition with an acceleration of four and L+S reconstruction can achieve a brain coverage of 40 slices at 2mm isotropic resolution and 64 x 64 matrix size every 500ms.

Safety and EEG data quality of concurrent high-density EEG and high-speed fMRI at 3 Tesla

PURPOSE

Concurrent EEG and fMRI is increasingly used to characterize the spatial-temporal dynamics of brain activity. However, most studies to date have been limited to conventional echo-planar imaging (EPI). There is considerable interest in integrating recently developed high-speed fMRI methods with high-density EEG to increase temporal resolution and sensitivity for task-based and resting state fMRI, and for detecting interictal spikes in epilepsy. In the present study using concurrent high-density EEG and recently developed high-speed fMRI methods, we investigate safety of radiofrequency (RF) related heating, the effect of EEG on cortical signal-to-noise ratio (SNR) in fMRI, and assess EEG data quality.

MATERIALS AND METHODS

The study compared EPI, multi-echo EPI, multi-band EPI and multi-slab echo-volumar imaging pulse sequences, using clinical 3 Tesla MR scanners from two different vendors that were equipped with 64- and 256-channel MR-compatible EEG systems, respectively, and receive only array head coils. Data were collected in 11 healthy controls (3 males, age range 18-70 years) and 13 patients with epilepsy (8 males, age range 21-67 years). Three of the healthy controls were scanned with the 256-channel EEG system, the other subjects were scanned with the 64-channel EEG system. Scalp surface temperature, SNR in occipital cortex and head movement were measured with and without the EEG cap. The degree of artifacts and the ability to identify background activity was assessed by visual analysis by a trained expert in the 64 channel EEG data (7 healthy controls, 13 patients).

RESULTS

RF induced heating at the surface of the EEG electrodes during a 30-minute scan period with stable temperature prior to scanning did not exceed 1.0° C with either EEG system and any of the pulse sequences used in this study. There was no significant decrease in cortical SNR due to the presence of the EEG cap (p > 0.05). No significant differences in the visually analyzed EEG data quality were found between EEG recorded during high-speed fMRI and during conventional EPI (p = 0.78). Residual ballistocardiographic artifacts resulted in 58% of EEG data being rated as poor quality.

CONCLUSION

This study demonstrates that high-density EEG can be safely implemented in conjunction with high-speed fMRI and that high-speed fMRI does not adversely affect EEG data quality. However, the deterioration of the EEG quality due to residual ballistocardiographic artifacts remains a significant constraint for routine clinical applications of concurrent EEG-fMRI.

On the detection of high frequency correlations in resting state fMRI

Current studies of resting-state connectivity rely on coherent signal fluctuations at frequencies below 0.1 Hz, however, recent studies using high-speed fMRI have shown that fluctuations above 0.5 Hz may exist. This study replicates the feasibility of measuring high frequency (HF) correlations in six healthy controls and a patient with a brain tumor while analyzing non-physiological signal sources via simulation. Resting-state data were acquired using a high-speed multi-slab echo-volumar imaging pulse sequence with 136 ms temporal resolution. Bandpass frequency filtering in combination with sliding window seed-based connectivity analysis using running mean of the correlation maps was employed to map HF correlations up to 3.7 Hz. Computer simulations of Rician noise and the underlying point spread function were analyzed to estimate baseline spatial autocorrelation levels in four major networks (auditory, sensorimotor, visual, and default-mode). Using seed regions based on Brodmann areas, the auditory and default-mode networks were observed to have significant frequency band dependent HF correlations above baseline spatial autocorrelation levels. Correlations in the sensorimotor network were at trend level. The auditory network was still observed using a unilateral single voxel seed. In the patient, HF auditory correlations showed a spatial displacement near the tumor consistent with the displacement seen at low frequencies. In conclusion, our data suggest that HF connectivity in the human brain may be observable with high-speed fMRI, however, the detection sensitivity may depend on the network observed, data acquisition technique, and analysis method.

Nuisance Regression of High-Frequency Functional Magnetic Resonance Imaging Data: Denoising Can Be Noisy

Recently, emerging studies have demonstrated the existence of brain resting-state spontaneous activity at frequencies higher than the conventional 0.1 Hz. A few groups utilizing accelerated acquisitions have reported persisting signals beyond 1 Hz, which seems too high to be accommodated by the sluggish hemodynamic process underpinning blood oxygen level-dependent contrasts (the upper limit of the canonical model is ∼0.3 Hz). It is thus questionable whether the observed high-frequency (HF) functional connectivity originates from alternative mechanisms (e.g., inflow effects, proton density changes in or near activated neural tissue) or rather is artificially introduced by improper preprocessing operations. In this study, we examined the influence of a common preprocessing step-whole-band linear nuisance regression (WB-LNR)-on resting-state functional connectivity (RSFC) and demonstrated through both simulation and analysis of real dataset that WB-LNR can introduce spurious network structures into the HF bands of functional magnetic resonance imaging (fMRI) signals. Findings of present study call into question whether published observations on HF-RSFC are partly attributable to improper data preprocessing instead of actual neural activities.

Increased fMRI Sensitivity at Equal Data Burden Using Averaged Shifted Echo Acquisition

There is growing evidence as to the benefits of collecting BOLD fMRI data with increased sampling rates. However, many of the newly developed acquisition techniques developed to collect BOLD data with ultra-short TRs require hardware, software, and non-standard analytic pipelines that may not be accessible to all researchers. We propose to incorporate the method of shifted echo into a standard multi-slice, gradient echo EPI sequence to achieve a higher sampling rate with a TR of <1 s with acceptable spatial resolution. We further propose to incorporate temporal averaging of consecutively acquired EPI volumes to both ameliorate the reduced temporal signal-to-noise inherent in ultra-fast EPI sequences and reduce the data burden. BOLD data were collected from 11 healthy subjects performing a simple, event-related visual-motor task with four different EPI sequences: (1) reference EPI sequence with TR = 1440 ms, (2) shifted echo EPI sequence with TR = 700 ms, (3) shifted echo EPI sequence with every two consecutively acquired EPI volumes averaged and effective TR = 1400 ms, and (4) shifted echo EPI sequence with every four consecutively acquired EPI volumes averaged and effective TR = 2800 ms. Both the temporally averaged sequences exhibited increased temporal signal-to-noise over the shifted echo EPI sequence. The shifted echo sequence with every two EPI volumes averaged also had significantly increased BOLD signal change compared with the other three sequences, while the shifted echo sequence with every four EPI volumes averaged had significantly decreased BOLD signal change compared with the other three sequences. The results indicated that incorporating the method of shifted echo into a standard multi-slice EPI sequence is a viable method for achieving increased sampling rate for collecting event-related BOLD data. Further, consecutively averaging every two consecutively acquired EPI volumes significantly increased the measured BOLD signal change and the subsequently calculated activation map statistics.

Enhanced subject-specific resting-state network detection and extraction with fast fMRI

Resting-state networks have become an important tool for the study of brain function. An ultra-fast imaging technique that allows to measure brain function, called Magnetic Resonance Encephalography (MREG), achieves an order of magnitude higher temporal resolution than standard echo-planar imaging (EPI). This new sequence helps to correct physiological artifacts and improves the sensitivity of the fMRI analysis. In this study, EPI is compared with MREG in terms of capability to extract resting-state networks. Healthy controls underwent two consecutive resting-state scans, one with EPI and the other with MREG. Subject-level independent component analyses (ICA) were performed separately for each of the two datasets. Using Stanford FIND atlas parcels as network templates, the presence of ICA maps corresponding to each network was quantified in each subject. The number of detected individual networks was significantly higher in the MREG data set than for EPI. Moreover, using short time segments of MREG data, such as 50 seconds, one can still detect and track consistent networks. Fast fMRI thus results in an increased capability to extract distinct functional regions at the individual subject level for the same scan times, and also allow the extraction of consistent networks within shorter time intervals than when using EPI, which is notably relevant for the analysis of dynamic functional connectivity fluctuations. Hum Brain Mapp 38:817-830, 2017. © 2016 Wiley Periodicals, Inc.

Marker-based ballistocardiographic artifact correction improves spike identification in EEG-fMRI of focal epilepsy patients

OBJECTIVES

Ballistocardiographic (BCG) artifacts resemble interictal epileptic discharges (IEDs) and can lead to incorrect IED identification in EEG-fMRI. This study investigates IEDs marked in EEGs corrected using information from a moiré phase tracking (MPT) marker.

METHODS

EEG-fMRI from 18 patients was processed with conventional methods for BCG removal, while 9 patients used a MPT marker. IEDs were marked first without ECG information. In a second review, suspicious IEDs synchronous with the BCG were discarded. After each review, an event-related fMRI analysis was performed on the marked IEDs.

RESULTS

No difference was found in the proportion of suspicious IEDs in the 2 patient groups. However, the distribution of IED timings was significantly related to the cardiac cycle in 11 of 18 patients recorded without MPT marker, but only 2 of 9 patients with marker. In patients recorded without marker, failing to discard suspicious IEDs led to more inaccurate fMRI maps and more distant activations.

CONCLUSIONS

BCG artifact correction based on MPT recordings allowed a more straightforward identification of IEDs that did not require ECG information in the large majority of patients.

SIGNIFICANCE

Marker-based ballistocardiographic artifact correction greatly facilitates the study of the generators of interictal discharges with EEG-fMRI.

Ultra-fast magnetic resonance encephalography of physiological brain activity - Glymphatic pulsation mechanisms?

The theory on the glymphatic convection mechanism of cerebrospinal fluid holds that cardiac pulsations in part pump cerebrospinal fluid from the peri-arterial spaces through the extracellular tissue into the peri-venous spaces facilitated by aquaporin water channels. Since cardiac pulses cannot be the sole mechanism of glymphatic propulsion, we searched for additional cerebrospinal fluid pulsations in the human brain with ultra-fast magnetic resonance encephalography. We detected three types of physiological mechanisms affecting cerebral cerebrospinal fluid pulsations: cardiac, respiratory, and very low frequency pulsations. The cardiac pulsations induce a negative magnetic resonance encephalography signal change in peri-arterial regions that extends centrifugally and covers the brain in ≈1 Hz cycles. The respiratory ≈0.3 Hz pulsations are centripetal periodical pulses that occur dominantly in peri-venous areas. The third type of pulsation was very low frequency (VLF 0.001-0.023 Hz) and low frequency (LF 0.023-0.73 Hz) waves that both propagate with unique spatiotemporal patterns. Our findings using critically sampled magnetic resonance encephalography open a new view into cerebral fluid dynamics. Since glymphatic system failure may precede protein accumulations in diseases such as Alzheimer's dementia, this methodological advance offers a novel approach to image brain fluid dynamics that potentially can enable early detection and intervention in neurodegenerative diseases.

Exploitation of temporal redundancy in compressed sensing reconstruction of fMRI studies with a prior-based algorithm (PICCS)

PURPOSE

Compressed sensing is a technique used to accelerate magnetic resonance imaging (MRI) acquisition without compromising image quality. While it has proven particularly useful in dynamic imaging procedures such as cardiac cine, very few authors have applied it to functional magnetic resonance imaging (fMRI). The purpose of the present study was to check whether the prior image constrained compressed sensing (PICCS) algorithm, which is based on an available prior image, can improve the statistical maps in fMRI better than other strategies that also exploit temporal redundancy.

METHODS

PICCS was compared to spatiotemporal total variation (TTV) and k-t FASTER, since they have already demonstrated high performance and robustness in other MRI applications, such as cardiac cine MRI and resting state fMRI, respectively. The prior image for PICCS was the average of all undersampled data. Both PICCS and TTV were solved using the split Bregman formulation. K-t FASTER algorithm relies on matrix completion to reconstruct the undersampled k-spaces. The three algorithms were evaluated using two datasets with high and low signal-to-noise ratio (SNR)-BOLD contrast-acquired in a 7 T preclinical MRI scanner and retrospectively undersampled at various rates (i.e., acceleration factors). The authors evaluated their performance in terms of the sensitivity/specificity of BOLD detection through receiver operating characteristic curves and by visual inspection of the statistical maps.

RESULTS

With high SNR studies, PICCS performed similarly to the state-of-the-art algorithms TTV and k-t FASTER and provided consistent BOLD signal at the ROI. In scenarios with low SNR and high acceleration factors, PICCS still provided consistent maps and higher sensitivity/specificity than TTV, whereas k-t FASTER failed to provide significant maps.

CONCLUSIONS

The authors performed a comparison between three reconstructions (PICCS, TTV, and k-t FASTER) that exploit temporal redundancy in fMRI. The prior-based algorithm, PICCS, preserved BOLD activation and sensitivity/specificity better than TTV and k-t FASTER in noisy scenarios. The PICCS algorithm can potentially reach an acceleration factor of ×8 and still provide BOLD contrast in the ROI with an area under the curve over 0.99.

Using fMRI to study reward processing in humans: past, present, and future

Functional magnetic resonance imaging (fMRI) is a noninvasive tool used to probe cognitive and affective processes. Although fMRI provides indirect measures of neural activity, the advent of fMRI has allowed for1) the corroboration of significant animal findings in the human brain, and2) the expansion of models to include more common human attributes that inform behavior. In this review, we briefly consider the neural basis of the blood oxygenation level dependent signal to set up a discussion of how fMRI studies have applied it in examining cognitive models in humans and the promise of using fMRI to advance such models. Specifically, we illustrate the contribution that fMRI has made to the study of reward processing, focusing on the role of the striatum in encoding reward-related learning signals that drive anticipatory and consummatory behaviors. For instance, we discuss how fMRI can be used to link neural signals (e.g., striatal responses to rewards) to individual differences in behavior and traits. While this functional segregation approach has been constructive to our understanding of reward-related functions, many fMRI studies have also benefitted from a functional integration approach that takes into account how interconnected regions (e.g., corticostriatal circuits) contribute to reward processing. We contend that future work using fMRI will profit from using a multimodal approach, such as combining fMRI with noninvasive brain stimulation tools (e.g., transcranial electrical stimulation), that can identify causal mechanisms underlying reward processing. Consequently, advancements in implementing fMRI will promise new translational opportunities to inform our understanding of psychopathologies.

Wavelet-based clustering of resting state MRI data in the rat

While functional connectivity has typically been calculated over the entire length of the scan (5-10min), interest has been growing in dynamic analysis methods that can detect changes in connectivity on the order of cognitive processes (seconds). Previous work with sliding window correlation has shown that changes in functional connectivity can be observed on these time scales in the awake human and in anesthetized animals. This exciting advance creates a need for improved approaches to characterize dynamic functional networks in the brain. Previous studies were performed using sliding window analysis on regions of interest defined based on anatomy or obtained from traditional steady-state analysis methods. The parcellation of the brain may therefore be suboptimal, and the characteristics of the time-varying connectivity between regions are dependent upon the length of the sliding window chosen. This manuscript describes an algorithm based on wavelet decomposition that allows data-driven clustering of voxels into functional regions based on temporal and spectral properties. Previous work has shown that different networks have characteristic frequency fingerprints, and the use of wavelets ensures that both the frequency and the timing of the BOLD fluctuations are considered during the clustering process. The method was applied to resting state data acquired from anesthetized rats, and the resulting clusters agreed well with known anatomical areas. Clusters were highly reproducible across subjects. Wavelet cross-correlation values between clusters from a single scan were significantly higher than the values from randomly matched clusters that shared no temporal information, indicating that wavelet-based analysis is sensitive to the relationship between areas.

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.

Evaluation of 2D multiband EPI imaging for high-resolution, whole-brain, task-based fMRI studies at 3T: Sensitivity and slice leakage artifacts

Functional magnetic resonance imaging (fMRI) studies that require high-resolution whole-brain coverage have long scan times that are primarily driven by the large number of thin slices acquired. Two-dimensional multiband echo-planar imaging (EPI) sequences accelerate the data acquisition along the slice direction and therefore represent an attractive approach to such studies by improving the temporal resolution without sacrificing spatial resolution. In this work, a 2D multiband EPI sequence was optimized for 1.5 mm isotropic whole-brain acquisitions at 3 T with 10 healthy volunteers imaged while performing simultaneous visual and motor tasks. The performance of the sequence was evaluated in terms of BOLD sensitivity and false-positive activation at multiband (MB) factors of 1, 2, 4, and 6, combined with in-plane GRAPPA acceleration of 2 × (GRAPPA 2), and the two reconstruction approaches of Slice-GRAPPA and Split Slice-GRAPPA. Sensitivity results demonstrate significant gains in temporal signal-to-noise ratio (tSNR) and t-score statistics for MB 2, 4, and 6 compared to MB 1. The MB factor for optimal sensitivity varied depending on anatomical location and reconstruction method. When using Slice-GRAPPA reconstruction, evidence of false-positive activation due to signal leakage between simultaneously excited slices was seen in one instance, 35 instances, and 70 instances over the ten volunteers for the respective accelerations of MB 2 × GRAPPA 2, MB 4 × GRAPPA 2, and MB 6 × GRAPPA 2. The use of Split Slice-GRAPPA reconstruction suppressed the prevalence of false positives significantly, to 1 instance, 5 instances, and 5 instances for the same respective acceleration factors. Imaging protocols using an acceleration factor of MB 2 × GRAPPA 2 can be confidently used for high-resolution whole-brain imaging to improve BOLD sensitivity with very low probability for false-positive activation due to slice leakage. Imaging protocols using higher acceleration factors (MB 3 or MB 4 × GRAPPA 2) can likely provide even greater gains in sensitivity but should be carefully optimized to minimize the possibility of false activations.

•

MB factors 1, 2, 4, and 6 and two reconstructions were evaluated for fMRI performance.

•

MB accelerations 2, 4, and 6 improved BOLD sensitivity metrics over MB 1.

•

False-positive activation due to signal leakage was seen at high accelerations.

•

Use of Split Slice-GRAPPA reconstruction significantly reduces false positives.

The neural basis of time-varying resting-state functional connectivity

Dynamic network analysis based on resting-state magnetic resonance imaging (rsMRI) is a fairly new and potentially powerful tool for neuroscience and clinical research. Dynamic analysis can be sensitive to changes that occur in psychiatric or neurologic disorders and can detect variations related to performance on individual trials in healthy subjects. However, the appearance of time-varying connectivity can also arise in signals that share no temporal information, complicating the interpretation of dynamic functional connectivity studies. Researchers have begun utilizing simultaneous imaging and electrophysiological recording to elucidate the neural basis of the networks and their variability in animals and in humans. In this article, we review findings that link changes in electrically recorded brain states to changes in the networks obtained with rsMRI and discuss some of the challenges inherent in interpretation of these studies. The literature suggests that multiple brain processes may contribute to the dynamics observed, and we speculate that it may be possible to separate particular aspects of the rsMRI signal to enhance sensitivity to certain types of neural activity, providing new tools for basic neuroscience and clinical research.