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

Test-retest reliability of peak location in the sensorimotor network of resting state fMRI for potential rTMS targets

Most stroke repetitive transcranial magnetic stimulation (rTMS) studies have used hand motor hotspots as rTMS stimulation targets; in addition, recent studies demonstrated that functional magnetic resonance imaging (fMRI) task activation could be used to determine suitable targets due to its ability to reveal individualized precise and stronger functional connectivity with motor-related brain regions. However, rTMS is unlikely to elicit motor evoked potentials in the affected hemisphere, nor would activity be detected when stroke patients with severe hemiplegia perform an fMRI motor task using the affected limbs. The current study proposed that the peak voxel in the resting-state fMRI (RS-fMRI) motor network determined by independent component analysis (ICA) could be a potential stimulation target. Twenty-one healthy young subjects underwent RS-fMRI at three visits (V1 and V2 on a GE MR750 scanner and V3 on a Siemens Prisma) under eyes-open (EO) and eyes-closed (EC) conditions. Single-subject ICA with different total number of components (20, 30, and 40) were evaluated, and then the locations of peak voxels on the left and right sides of the sensorimotor network (SMN) were identified. While most ICA RS-fMRI studies have been carried out on the group level, that is, Group-ICA, the current study performed individual ICA because only the individual analysis could guide the individual target of rTMS. The intra- (test-retest) and inter-scanner reliabilities of the peak location were calculated. The use of 40 components resulted in the highest test-retest reliability of the peak location in both the left and right SMN compared with that determined when 20 and 30 components were used for both EC and EO conditions. ICA with 40 components might be another way to define a potential target in the SMN for poststroke rTMS treatment.

Cardiovascular Pulsatility Increases in Visual Cortex Before Blood Oxygen Level Dependent Response During Stimulus

The physiological pulsations that drive tissue fluid homeostasis are not well characterized during brain activation. Therefore, we used fast magnetic resonance encephalography (MREG) fMRI to measure full band (0–5 Hz) blood oxygen level-dependent (BOLD_FB) signals during a dynamic visual task in 23 subjects. This revealed brain activity in the very low frequency (BOLD_VLF) as well as in cardiac and respiratory bands. The cardiovascular hemodynamic envelope (CHe) signal correlated significantly with the visual BOLD_VLF response, considered as an independent signal source in the V1-V2 visual cortices. The CHe preceded the canonical BOLD_VLF response by an average of 1.3 (± 2.2) s. Physiologically, the observed CHe signal could mark increased regional cardiovascular pulsatility following vasodilation.

Advances in resting state fMRI acquisitions for functional connectomics

Resting state functional magnetic resonance imaging (rs-fMRI) is based on spontaneous fluctuations in the blood oxygen level dependent (BOLD) signal, which occur simultaneously in different brain regions, without the subject performing an explicit task. The low-frequency oscillations of the rs-fMRI signal demonstrate an intrinsic spatiotemporal organization in the brain (brain networks) that may relate to the underlying neural activity. In this review article, we briefly describe the current acquisition techniques for rs-fMRI data, from the most common approaches for resting state acquisition strategies, to more recent investigations with dedicated hardware and ultra-high fields. Specific sequences that allow very fast acquisitions, or multiple echoes, are discussed next. We then consider how acquisition methods weighted towards specific parts of the BOLD signal, like the Cerebral Blood Flow (CBF) or Volume (CBV), can provide more spatially specific network information. These approaches are being developed alongside the commonly used BOLD-weighted acquisitions. Finally, specific applications of rs-fMRI to challenging regions such as the laminae in the neocortex, and the networks within the large areas of subcortical white matter regions are discussed. We finish the review with recommendations for acquisition strategies for a range of typical applications of resting state fMRI.

Topography-Related EEG-fMRI in Surgically Confirmed Epileptic Foci: A Comparison to Spike-Related EEG-fMRI in Clinical Practice

EEG-fMRI has gained increasing importance in epilepsy pre-surgical diagnosis. However, 40-70% of EEG-fMRI recordings in patients lack interictal epileptiform discharges (IEDs) during the scan, which could be overcome by detecting matching topography maps. We tried to validate this method in clinical settings taking various electroclinical factors into consideration. Eleven patients who had undergone EEG-fMRI during pre-surgical evaluation for drug-resistant epilepsy and who had had clinical long-term video-EEG were studied. Spike-related blood oxygen level-dependent (BOLD) maps were created using IEDs occurring during the EEG-fMRI scan. Separate maps were then generated from IEDs marked on the clinical long-term EEG recordings, which were averaged to produce topographical IED maps and correlated with the EEGs recorded inside the scanner yielding a correlation coefficient time course. Epileptogenic zones were defined by an expert panel during pre-surgical evaluation and validated by an epilepsy surgery resulting in a good outcome. Both techniques' performance was evaluated according to factors including arousal during IED recording, IED topography and lateralization, lesion type, and localization. Topography-related EEG-fMRI yielded more specific results compared to the spike-related method. Superficial lesion location and ipsilateral IED seem to result in a higher concordance of BOLD maps. The polarity of BOLD responses may be lesion-dependent, and both positive and negative BOLD changes may be associated with the irritative zone. Topography-related EEG-fMRI may show improved specificity especially for superficial lesions producing ipsilateral spikes. This method can be used as an alternative either in the absence of spikes during the simultaneous EEG-fMRI acquisition or to sharpen a diffusely activated BOLD-map.

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.

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.

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.

Single-subject manual independent component analysis and resting state fMRI connectivity outcomes in patients with juvenile absence epilepsy

The quality of fMRI data impacts functional connectivity measures and consequently, the decisions that clinicians and researchers make regarding functional connectivity interpretation. The present study used resting state fMRI to investigate resting state network connectivity in a sample of patients with Juvenile Absence Epilepsy. Single-subject manual independent component analysis was used in two levels, whereby all noise components were removed, and cerebrospinal fluid pulsation components only were isolated and removed. Improved temporal signal to noise ratios and functional connectivity metrics were observed in each of the cleaning levels for both epilepsy and control cohorts. Results showed full, single-subject manual independent component analysis reduced the number of functional connectivity correlations and increased the strength of these correlations. Similar effects were also observed for the cerebrospinal fluid pulsation only cleaned data relative to the uncleaned, and fully cleaned data. Single-subject manual independent component analysis coupled with short TR multiband acquisition can significantly improve the validity of findings derived from fMRI data sets.

Direct modelling of gradient artifacts for EEG-fMRI denoising and motion tracking

OBJECTIVE

Simultaneous electroencephalography and functional magnetic resonance imaging recording (EEG-fMRI) has been widely used in neuroscientific and clinical research. The artifacts in the recorded EEG resulting from rapidly switching magnetic field gradients are usually corrected by average-artifact subtraction (AAS) due to their repetitive nature. But the performance of AAS is often disrupted by altered artifact waveforms across epochs, notably due to head motion.

APPROACH

Here, a method is proposed to make use of the known MR sequence gradient waveforms for a direct modelling of gradient artifacts. After accounting for filtering effects on the gradient artifacts, a continuous modulation of the gradient waveforms superimposed on the EEG signal is obtained.

MAIN RESULTS

Although a moving AAS template can adjust to slow drifts in gradient artifact variation, it fails to adapt to abrupt motion, resulting in residual noise. We demonstrate how this modelling approach can reduce motion-affected gradient artifacts without distorting the underlying neuronal signals. Moreover, the method provides useful head motion information highly correlated with motion tracked by an optical camera.

SIGNIFICANCE

Our work provides a novel way to improve gradient artifact removal in EEG-fMRI, and shows a potential to detect head motion without requiring additional hardware.

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.

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.

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.