Visual Stimuli Modulate Local Field Potentials But Drive No High-Frequency Activity in Human Auditory Cortex.

Neuroimaging studies suggest cross-sensory visual influences in human auditory cortices (ACs). Whether these influences reflect active visual processing in human ACs, which drives neuronal firing and concurrent broadband high-frequency activity (BHFA; >70 Hz), or whether they merely modulate sound processing is still debatable. Here, we presented auditory, visual, and audiovisual stimuli to 16 participants (7 women, 9 men) with stereo-EEG depth electrodes implanted near ACs for presurgical monitoring. Anatomically normalized group analyses were facilitated by inverse modeling of intracranial source currents. Analyses of intracranial event-related potentials (iERPs) suggested cross-sensory responses to visual stimuli in ACs, which lagged the earliest auditory responses by several tens of milliseconds. Visual stimuli also modulated the phase of intrinsic low-frequency oscillations and triggered 15-30 Hz event-related desynchronization in ACs. However, BHFA, a putative correlate of neuronal firing, was not significantly increased in ACs after visual stimuli, not even when they coincided with auditory stimuli. Intracranial recordings demonstrate cross-sensory modulations, but no indication of active visual processing in human ACs.

Ballistocardiogram suppression in concurrent EEG-MRI by dynamic modeling of heartbeats.

The ballistocardiogram (BCG), the induced electric potentials by the head motion originating from heartbeats, is a prominent source of noise in electroencephalography (EEG) data during magnetic resonance imaging (MRI). Although methods have been proposed to suppress the BCG artifact, more work considering the variability of cardiac cycles and head motion across time and subjects is needed to provide highly robust correction. Here, a method called "dynamic modeling of heartbeats" (DMH) is proposed to reduce BCG artifacts in EEG data recorded inside an MRI system. The DMH method models BCG artifacts by combining EEG points at time instants with similar dynamics. The modeled BCG artifact is then subtracted from the EEG recording to suppress the BCG artifact. Performance of DMH was tested and specifically compared with the Optimal Basis Set (OBS) method on EEG data recorded inside a 3T MRI system with either no MRI acquisition (Inside-MRI), echo-planar imaging (EPI-EEG), or fast MRI acquisition using simultaneous multi-slice and inverse imaging methods (SMS-InI-EEG). In a steady-state visual evoked response (SSVEP) paradigm, the 15-Hz oscillatory neuronal activity at the visual cortex after DMH processing was about 130% of that achieved by OBS processing for Inside-MRI, SMS-InI-EEG, and EPI-EEG conditions. The DMH method is computationally efficient for suppressing BCG artifacts and in the future may help to improve the quality of EEG data recorded in high-field MRI systems for neuroscientific and clinical applications.

Distributed source modeling of intracranial stereoelectro-encephalographic measurements.

Intracranial stereoelectroencephalography (sEEG) provides unsurpassed sensitivity and specificity for human neurophysiology. However, functional mapping of brain functions has been limited because the implantations have sparse coverage and differ greatly across individuals. Here, we developed a distributed, anatomically realistic sEEG source-modeling approach for within- and between-subject analyses. In addition to intracranial event-related potentials (iERP), we estimated the sources of high broadband gamma activity (HBBG), a putative correlate of local neural firing. Our novel approach accounted for a significant portion of the variance of the sEEG measurements in leave-one-out cross-validation. After logarithmic transformations, the sensitivity and signal-to-noise ratio were linearly inversely related to the minimal distance between the brain location and electrode contacts (slope≈-3.6). The signa-to-noise ratio and sensitivity in the thalamus and brain stem were comparable to those locations at the vicinity of electrode contact implantation. The HGGB source estimates were remarkably consistent with analyses of intracranial-contact data. In conclusion, distributed sEEG source modeling provides a powerful neuroimaging tool, which facilitates anatomically-normalized functional mapping of human brain using both iERP and HBBG data.

Concurrent electrophysiological and hemodynamic measurements of evoked neural oscillations in human visual cortex using sparsely interleaved fast fMRI and EEG.

Electroencephalography (EEG) concurrently collected with functional magnetic resonance imaging (fMRI) is heavily distorted by the repetitive gradient coil switching during the fMRI acquisition. The performance of the typical template-based gradient artifact suppression method can be suboptimal because the artifact changes over time. Gradient artifact residuals also impede the subsequent suppression of ballistocardiography artifacts. Here we propose recording continuous EEG with temporally sparse fast fMRI (fast fMRI-EEG) to minimize the EEG artifacts caused by MRI gradient coil switching without significantly compromising the field-of-view and spatiotemporal resolution of fMRI. Using simultaneous multi-slice inverse imaging to achieve whole-brain fMRI with isotropic 5-mm resolution in 0.1 â€‹s, and performing these acquisitions once every 2 â€‹s, we have 95% of the duty cycle available to record EEG with substantially less gradient artifact. We found that the standard deviation of EEG signals over the entire acquisition period in fast fMRI-EEG was reduced to 54% of that in conventional concurrent echo-planar imaging (EPI) and EEG recordings (EPI-EEG) across participants. When measuring 15-Hz steady-state visual evoked potentials (SSVEPs), the baseline-normalized oscillatory neural response in fast fMRI-EEG was 2.5-fold of that in EPI-EEG. The functional MRI responses associated with the SSVEP delineated by EPI and fast fMRI were similar in the spatial distribution, the elicited waveform, and detection power. Sparsely interleaved fast fMRI-EEG provides high-quality EEG without substantially compromising the quality of fMRI in evoked response measurements, and has the potential utility for applications where the onset of the target stimulus cannot be precisely determined, such as epilepsy.

An orthogonal shim coil for 3T brain imaging.

PURPOSE: We designed and implemented an orthogonal shim array consisting of shim coils with their planes perpendicular to the planes of neighboring RF coils. This shim coil improves the magnetic field homogeneity by minimizing the interference to RF coils. METHODS: Using realistic off-resonance maps of the human brain, we first evaluated the performance of shim coils in different orientations. Based on simulations, we developed a 7-channel orthogonal shim array, whose coil plan was perpendicular to neighboring RF coils, at the forehead. A programmable open-source current driver supplied shim currents. RESULTS: The 7-channel orthogonal shim array caused only marginal SNR loss to the integrated 32-channel RF-shim array. The 7-channel orthogonal shim array itself improved the magnetic field homogeneity by 30% in slice-optimized shimming, comparable to the baseline shimming offered by the scanner's 2nd order spherical harmonic shimming. CONCLUSION: Orthogonal shim coils can improve the field homogeneity while maintaining high image SNR.

Reduction of lipid contamination in MR spectroscopy imaging using signal space projection.

PURPOSE: Lipid contamination can complicate the metabolite quantification in MR spectroscopic imaging (MRSI). In addition to various experimental methods demonstrated to be feasible for lipid suppression, the postprocessing method is beneficial in the flexibility of applications. In this study, the signal space projection (SSP) algorithm is proposed to suppress the lipid signal in the MRSI. METHODS: The performance of lipid suppression using SSP and SSP combined with the Papoulis-Gerchberg (PG) algorithm (PG+SSP) is examined in 2D MRSI data and the results were compared with outer volume saturation (OVS) methods. Up to 10 lipid spatial components were extracted by SSP from lipid signals in the range of 0.8~1.5 ppm. RESULTS: Our results show that most lipid signals were found in the first 4 to 5 components and that lipid signals on the spectra can be suppressed using 4 to 5 components. Metabolites concentrations were quantified using LCModel. Two regions of interest (ROIs) were manually selected on the peripheral and inner brain regions. The quantification of metabolites in terms of fitting reliability (CRLB) and spatial variations within ROIs (SpaVar) is improved using SSP. When 5 to 6 components were used in SSP and PG+SSP, the metabolite concentrations and the associated SpaVar and CRLB are at the same level as those from the OVS. CONCLUSION: We have demonstrated that the SSP method can be used to suppress the lipid signals of MRSI and SSP with 5 to 6 components is suggested to have a similar suppression performance as the OVS method.

Relative latency and temporal variability of hemodynamic responses at the human primary visual cortex.

The blood-oxygen-level-dependent (BOLD) functional MRI (fMRI) signal is a robust surrogate for local neuronal activity. However, it has been shown to vary substantially across subjects, brain regions, and repetitive measurements. This variability represents a limit to the precision of the BOLD response and the ability to reliably discriminate brain hemodynamic responses elicited by external stimuli or behavior that are nearby in time. While the temporal variability of the BOLD signal at human visual cortex has been found in the range of a few hundreds of milliseconds, the spatial distributions of the average and standard deviation of this temporal variability have not been quantitatively characterized. Here we use fMRI measurements with a high sampling rate (10Hz) to map the latency, intra- and inter-subject variability of the evoked BOLD signal in human primary (V1) visual cortices using an event-related fMRI paradigm. The latency relative to the average BOLD signal evoked by 30 stimuli was estimated to be 0.03±0.20s. Within V1, the absolute value of the relative BOLD latency was found correlated to intra- and inter-subject temporal variability. After comparing these measures to retinotopic maps, we found that locations with V1 areas sensitive to smaller eccentricity have later responses and smaller inter-subject variabilities. These correlations were found from data with either short inter-stimulus interval (ISI; average 4s) or long ISI (average 30s). Maps of the relative latency as well as inter-/intra-subject variability were found visually asymmetric between hemispheres. Our results suggest that the latency and variability of regional BOLD signal measured with high spatiotemporal resolution may be used to detect regional differences in hemodynamics to inform fMRI studies. However, the physiological origins of timing index distributions and their hemispheric asymmetry remain to be investigated.

Mitigation of B1+ inhomogeneity using spatially selective excitation with jointly designed quadratic spatial encoding magnetic fields and RF shimming.

PURPOSE: The inhomogeneity of flip angle distribution is a major challenge impeding the application of high-field MRI. We report a method combining spatially selective excitation using generalized spatial encoding magnetic fields (SAGS) with radiofrequency (RF) shimming to achieve homogeneous excitation. This method can be an alternative approach to address the challenge of B1+ inhomogeneity using nonlinear gradients. METHODS: We proposed a two-step algorithm that jointly optimizes the combination of nonlinear spatial encoding magnetic fields and the combination of multiple RF transmitter coils and then optimizes the locations, RF amplitudes, and phases of the spokes. RESULTS: Our results show that jointly designed SAGS and RF shimming can provide a more homogeneous flip angle distribution than using SAGS or RF shimming alone. Compared with RF shimming alone, our approach can reduce the relative standard deviation of flip angle by 56% and 52% using phantom and human head data, respectively. CONCLUSION: The jointly designed SAGS and RF shimming method can be used to achieve homogeneous flip angle distributions when fully parallel RF transmission is not available. Magn Reson Med 78:577-587, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Rotary scanning acquisition in ultra-low-field MRI.

PURPOSE: To develop a method of achieving large field of view (FOV) imaging with a smaller amount of data in ultra-low-field (ULF) MRI. THEORY: In rotary scanning acquisition (RSA), data from the imaging object is acquired at multiple angles by rotating the object or the scanner. RSA is similar to radial-trajectory acquisition but simplifies the measurement and image reconstruction when concomitant fields are nonnegligible. METHODS: RSA was implemented to achieve large FOV with only three localized superconductive quantum interference device (SQUID) sensors at the ULF-MRI field of 50 µT. RESULTS: Simulations suggest benefits of RSA, including reduced concomitant field artifacts, large FOV imaging, and SNR improvement. Experimental data demonstrate the feasibility of reconstructing large FOV images using only three SQUID sensors with 33% of the amount of data collected using a Cartesian trajectory. CONCLUSION: RSA can be useful in low-field, low-weight, or portable MRI to generate large FOV images with only a few sensors. Magn Reson Med 75:2255-2264, 2016. © 2015 Wiley Periodicals, Inc.

Significant feed-forward connectivity revealed by high frequency components of BOLD fMRI signals.

Granger causality analysis has been suggested as a method of estimating causal modulation without specifying the direction of information flow a priori. Using BOLD-contrast functional MRI (fMRI) data, such analysis has been typically implemented in the time domain. In this study, we used magnetic resonance inverse imaging, a method of fast fMRI enabled by massively parallel detection allowing up to 10 Hz sampling rate, to investigate the causal modulation at different frequencies up to 5 Hz. Using a visuomotor two-choice reaction-time task, both the spectral decomposition of Granger causality and isolated effective coherence revealed that the BOLD signal at frequency up to 3 Hz can still be used to estimate significant dominant directions of information flow consistent with results from the time-domain Granger causality analysis. We showed the specificity of estimated dominant directions of information flow at high frequencies by contrasting causality estimates using data collected during the visuomotor task and resting state. Our data suggest that hemodynamic responses carry physiological information related to inter-regional modulation at frequency higher than what has been commonly considered.

Combining parallel detection of proton echo planar spectroscopic imaging (PEPSI) measurements with a data-consistency constraint improves SNR.

One major challenge of MRSI is the poor signal-to-noise ratio (SNR), which can be improved by using a surface coil array. Here we propose to exploit the spatial sensitivity of different channels of a coil array to enforce the k-space data consistency (DC) in order to suppress noise and consequently to improve MRSI SNR. MRSI data were collected using a proton echo planar spectroscopic imaging (PEPSI) sequence at 3 T using a 32-channel coil array and were averaged with one, two and eight measurements (avg-1, avg-2 and avg-8). The DC constraint was applied using a regularization parameter λ of 1, 2, 3, 5 or 10. Metabolite concentrations were quantified using LCModel. Our results show that the suppression of noise by applying the DC constraint to PEPSI reconstruction yields up to 32% and 27% SNR gain for avg-1 and avg-2 data with λ = 5, respectively. According to the reported Cramer-Rao lower bounds, the improvement in metabolic fitting was significant (p < 0.01) when the DC constraint was applied with λ ≥ 2. Using the DC constraint with λ = 3 or 5 can minimize both root-mean-square errors and spatial variation for all subjects using the avg-8 data set as reference values. Our results suggest that MRSI reconstructed with a DC constraint can save around 70% of scanning time to obtain images and spectra with similar SNRs using λ = 5.

Enhanced spontaneous oscillations in the supplementary motor area are associated with sleep-dependent offline learning of finger-tapping motor-sequence task.

Sleep is beneficial for various types of learning and memory, including a finger-tapping motor-sequence task. However, methodological issues hinder clarification of the crucial cortical regions for sleep-dependent consolidation in motor-sequence learning. Here, to investigate the core cortical region for sleep-dependent consolidation of finger-tapping motor-sequence learning, while human subjects were asleep, we measured spontaneous cortical oscillations by magnetoencephalography together with polysomnography, and source-localized the origins of oscillations using individual anatomical brain information from MRI. First, we confirmed that performance of the task at a retest session after sleep significantly increased compared with performance at the training session before sleep. Second, spontaneous δ and fast-σ oscillations significantly increased in the supplementary motor area (SMA) during post-training compared with pretraining sleep, showing significant and high correlation with the performance increase. Third, the increased spontaneous oscillations in the SMA correlated with performance improvement were specific to slow-wave sleep. We also found that correlations of δ oscillation between the SMA and the prefrontal and between the SMA and the parietal regions tended to decrease after training. These results suggest that a core brain region for sleep-dependent consolidation of the finger-tapping motor-sequence learning resides in the SMA contralateral to the trained hand and is mediated by spontaneous δ and fast-σ oscillations, especially during slow-wave sleep. The consolidation may arise along with possible reorganization of a larger-scale cortical network that involves the SMA and cortical regions outside the motor regions, including prefrontal and parietal regions.

Improving the spatial resolution of magnetic resonance inverse imaging via the blipped-CAIPI acquisition scheme.

Using simultaneous acquisition from multiple channels of a radio-frequency (RF) coil array, magnetic resonance inverse imaging (InI) achieves functional MRI acquisitions at a rate of 100ms per whole-brain volume. InI accelerates the scan by leaving out partition encoding steps and reconstructs images by solving under-determined inverse problems using RF coil sensitivity information. Hence, the correlated spatial information available in the coil array causes spatial blurring in the InI reconstruction. Here, we propose a method that employs gradient blips in the partition encoding direction during the acquisition to provide extra spatial encoding in order to better differentiate signals from different partitions. According to our simulations, this blipped-InI (bInI) method can increase the average spatial resolution by 15.1% (1.3mm) across the whole brain and from 32.6% (4.2mm) in subcortical regions, as compared to the InI method. In a visual fMRI experiment, we demonstrate that, compared to InI, the spatial distribution of bInI BOLD response is more consistent with that of a conventional echo-planar imaging (EPI) at the level of individual subjects. With the improved spatial resolution, especially in subcortical regions, bInI can be a useful fMRI tool for obtaining high spatiotemporal information for clinical and cognitive neuroscience studies.

Mitigate B1+ inhomogeneity using spatially selective radiofrequency excitation with generalized spatial encoding magnetic fields.

PURPOSE: High-field magnetic resonance imaging (MRI) has the challenge of inhomogeneous B(1)(+), and consequently inhomogeneous flip angle distribution, which causes spatially dependent contrast and makes clinical diagnosis difficult. METHOD: We propose a two-step pulse design procedure in which (1) a combination of linear and nonlinear spatial encoding magnetic fields (SEMs) is used to remap the B(1)(+) map in order to reduce the dimensionality of the problem, (2) the locations, amplitudes, and phases of spoke pulses are estimated in one dimension. The advantage of this B(1)(+) remapping is that when the isointensity contours of a linear combination of SEMs are similar to the isointensity contours of B(1)(+), a simple pulse sequence design using time-varying SEMs can achieve a homogenous flip-angle distribution efficiently. RESULTS: We demonstrate that spatially selective radiofrequency (RF) excitation with generalized SEMs (SAGS) using both linear and quadratic SEMs in a multi-spoke k-space trajectory can mitigate the B(1)(+) inhomogeneity at 7T efficiently. Numerical simulations based on experimental data suggest that, compared with other methods, SAGS provide a formulation allowing multiple-pulse design, a similar average flip-angle distribution with less RF power, and/or a more homogeneous flip-angle distribution. CONCLUSION: Without using multiple RF coils for parallel transmission, SAGS can be used to mitigate the B(1)(+) inhomogeneity in high-field MRI experiments.

Efficient concomitant and remanence field artifact reduction in ultra-low-field MRI using a frequency-space formulation.

PURPOSE: For ultra-low-field MRI, the spatial-encoding magnetic fields generated by gradient coils can have strong concomitant fields leading to prominent image distortion. Additionally, using superconducting magnet to pre-polarize magnetization can improve the signal-to-noise ratio of ultra-low-field MRI. Yet the spatially inhomogeneous remanence field due to the permanently trapped flux inside a superconducting pre-polarizing coil modulates magnetization and causes further image distortion. METHOD: We propose a two-stage frequency-space (f-x) formulation to accurately describe the dynamics of spatially-encoded magnetization under the influence of concomitant and remanence fields, which allows for correcting image distortion due to concomitant and remanence fields. RESULTS: Our method is computationally efficient as it uses a combination of the fast Fourier transform algorithm and a linear equation solver. With sufficiently dense discretization in solving the linear equation, the performance of this f-x method was found to be stable among different choices of the regularization parameter and the regularization matrix. CONCLUSION: We present this method together with numerical simulations and experimental data to demonstrate how concomitant and remanence field artifacts in ultra-low-field MRI can be corrected efficiently.

Attention-driven auditory cortex short-term plasticity helps segregate relevant sounds from noise.

How can we concentrate on relevant sounds in noisy environments? A "gain model" suggests that auditory attention simply amplifies relevant and suppresses irrelevant afferent inputs. However, it is unclear whether this suffices when attended and ignored features overlap to stimulate the same neuronal receptive fields. A "tuning model" suggests that, in addition to gain, attention modulates feature selectivity of auditory neurons. We recorded magnetoencephalography, EEG, and functional MRI (fMRI) while subjects attended to tones delivered to one ear and ignored opposite-ear inputs. The attended ear was switched every 30 s to quantify how quickly the effects evolve. To produce overlapping inputs, the tones were presented alone vs. during white-noise masking notch-filtered ±1/6 octaves around the tone center frequencies. Amplitude modulation (39 vs. 41 Hz in opposite ears) was applied for "frequency tagging" of attention effects on maskers. Noise masking reduced early (50-150 ms; N1) auditory responses to unattended tones. In support of the tuning model, selective attention canceled out this attenuating effect but did not modulate the gain of 50-150 ms activity to nonmasked tones or steady-state responses to the maskers themselves. These tuning effects originated at nonprimary auditory cortices, purportedly occupied by neurons that, without attention, have wider frequency tuning than ±1/6 octaves. The attentional tuning evolved rapidly, during the first few seconds after attention switching, and correlated with behavioral discrimination performance. In conclusion, a simple gain model alone cannot explain auditory selective attention. In nonprimary auditory cortices, attention-driven short-term plasticity retunes neurons to segregate relevant sounds from noise.

Acoustic noise generated by TMS in typical environment and inside an MRI scanner.

BACKGROUND: The operation of a transcranial magnetic stimulation (TMS) coil produces high-intensity impulse sounds. In TMS, a magnetic field is generated by a short-duration pulse in the range of thousands of amperes in the TMS coil. When placed in a strong magnetic field, such as inside a magnetic resonance imaging (MRI) bore, the interaction of the magnetic field and the current in the TMS coil can cause strong forces on the coil casing. The strengths of these forces depend on the coil orientation in the main magnetic field (B0). Part of the energy in this process is dissipated in the form of acoustic noise. OBJECTIVE: Our objective was to measure the sound pressure levels (SPL) of TMS "click" sounds created by commercial TMS stimulators and coils in a typical environment and inside a 3-T MRI scanner and advance the knowledge of the acoustic behaviour of TMS to safely conduct TMS alone as well as concurrently with functional MRI (fMRI). METHODS: We report SPL measurements of two commercial MRI-compatible TMS systems in the 3-T B0 field of an MRI scanner and in the earth's magnetic field. Also, we present the acoustic noise measurements of four commercial TMS stimulators and three different TMS coils in a typical operational environment without the B0 field. RESULTS: The maximum peak SPL measured was 158 dB(C) inside the 3-T MRI scanner. Outside the scanner, the maximum peak SPL was 117 dB(C). Inside the scanner, the peak SPL increased by 21-45 dB(C) depending on the stimulator and the orientation of the electric field relative to the B field. CONCLUSIONS: Hearing protection is obligatory during concurrent TMS-fMRI experiments and highly recommended during any TMS experiment. The manufacturing of quieter TMS systems is encouraged to reduce the risk of hearing damage and other unwanted effects.

Omnidirectional Monolithic Marker for Intra-Operative MR-Based Positional Sensing in Closed MRI.

We present a design of an inductively coupled radio frequency (ICRF) marker for magnetic resonance (MR)-based positional tracking, enabling the robust increase of tracking signal at all scanning orientations in quadrature-excited closed MR imaging (MRI). The marker employs three curved resonant circuits fully covering a cylindrical surface that encloses the signal source. Each resonant circuit is a planar spiral inductor with parallel plate capacitors fabricated monolithically on flexible printed circuit board (FPC) and bent to achieve the curved structure. Size of the constructed marker is Ø3-mm ×5 -mm with quality factor > 22, and its tracking performance was validated with 1.5 T MRI scanner. As result, the marker remains as a high positive contrast spot under 360° rotations in 3 axes. The marker can be accurately localized with a maximum error of 0.56 mm under a displacement of 56 mm from the isocenter, along with an inherent standard deviation of 0.1-mm. Accrediting to the high image contrast, the presented marker enables automatic and real-time tracking in 3D without dependency on its orientation with respect to the MRI scanner receive coil. In combination with its small form-factor, the presented marker would facilitate robust and wireless MR-based tracking for intervention and clinical diagnosis. This method targets applications that can involve rotational changes in all axes (X-Y-Z).

Distributed source modeling of stereoencephalographic measurements of ictal activity.

OBJECTIVES: Stereoelectroencephalography (SEEG) can define the epileptogenic zone (EZ). However, SEEG is susceptible to the sampling bias, where no SEEG recording is taken within a circumscribed EZ. METHODS: Nine patients with medically refractory epilepsy underwent SEEG recording, and brain resection got positive outcomes. Ictal neuronal currents were estimated by distributed source modeling using the SEEG data and individual's anatomical magnetic resonance imaging. Using a retrospective leave-one-out data sub-sampling, we evaluated the sensitivity and specificity of the current estimates using MRI after surgical resection or radio-frequency ablation. RESULTS: The sensitivity and specificity in detecting the EZ were indistinguishable from either the data from all electrodes or the sub-sampled data (rank sum test: rank sum = 23719, p = 0.13) when at least one remaining electrode contact was no more than 20 mm away. CONCLUSIONS: The distributed neuronal current estimates of ictal SEEG data can mitigate the challenge of delineating the boundary of the EZ in cases of missing an electrode implanted within the EZ and a required second SEEG exploration. SIGNIFICANCE: Distributed source modeling can be a tool for clinicians to infer the EZ by allowing for more flexible planning of the electrode implantation route and minimizing the number of electrodes.

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

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