Stimulus-induced Rotary Saturation (SIRS): a potential method for the detection of neuronal currents with MRI

Homogeneous B0 coil design method for open-access ultra-low field magnetic resonance imaging: A simulation study

A multimodal brain function measurement system integrating functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) is expected to be a tool that will provide new insights into neuroscience. To integrate fMRI and MEG, an ultra-low-field MRI (ULF-MRI) scanner that can generate a static magnetic field (B0) with an electromagnetic coil and turn off the B0 during MEG measurements is desirable. While electromagnetic B0 coil has the above advantages, it also has a trade-off between size and the broadness of the magnetic field homogeneity. In this study, we proposed a method for designing a B0 multi-stage circular coil arrangement that determines the number of coils required to maximize magnetic field homogeneity and minimize the total wiring length of the coils. The optimized multi-stage coil arrangement had an external shape of 600 mm in diameter and a maximum height of 600 mm, with an aperture of 600 mm in diameter and 300 mm in height. The magnetic field homogeneity was <100 ppm over a 210 mm diameter spherical volume (DSV). Compared to a previous two coil pairs arrangement with the same magnetic field homogeneity, the diameter was 1/1.9 times smaller, indicating that the newly designed B0 coil arrangement realized a smaller size and wider magnetic field homogeneity.

Neuronal current imaging of epileptic activity: An MRI study in patients with a first unprovoked epileptic seizure

OBJECTIVE

This study evaluates the performance of the novel MRI sequence stimulus-induced rotary saturation (SIRS) to map responses to interictal epileptic activity in the human cortex. Spin-lock pulses have been applied to indirectly detect neuronal activity through magnetic field perturbations. Following initial reports about the feasibility of the method in humans and animals with epilepsy, we aimed to investigate the diagnostic yield of spin-lock MR pulses in comparison with scalp-EEG in first seizure patients.

METHODS

We employed a novel method for measurements of neuronal activity through the detection of a resonant oscillating field, stimulus-induced rotary saturation contrast (SIRS) at spin-lock frequencies of 120 and 240 Hz acquired at a single 3T MRI system. Within a prospective observational study, we conducted SIRS experiments in 55 patients within 7 days after a suspected first unprovoked epileptic seizure and 61 healthy control subjects. In this study, we report on the analysis of data from a single 3T MRI system, encompassing 35 first seizure patients and 31 controls.

RESULTS

The SIRS method was applicable in all patients and healthy controls at frequencies of 120 and 240 Hz. We did not observe any significant age- or sex-related differences. Specificity of SIRS at 120 Hz was 90.3% and 93.5% at 240 Hz. Sensitivity was 17.1% at 120 Hz and 40.0% at 240 Hz.

SIGNIFICANCE

SIRS targets neuronal oscillating magnetic fields in patients with epilepsy. The coupling of presaturated spins to epilepsy-related magnetic field perturbations may serve as a-at this stage experimental-diagnostic test in first seizure patients to complement EEG findings as a standard screening test.

PLAIN LANGUAGE SUMMARY

Routine diagnostic tests carry several limitations when applied after a suspected first seizure. SIRS is a noninvasive MRI method to enable time-sensitive diagnosis of image correlates of epileptic activity with increased sensitivity compared to routine EEG.

Rotary excitation of non-sinusoidal pulsed magnetic fields: Towards non-invasive direct detection of cardiac conduction

PURPOSE

There is a need for high resolution non-invasive imaging methods of physiologic magnetic fields. The purpose of this work is to develop a MRI detection approach for non-sinusoidal magnetic fields based on the rotary excitation (REX) mechanism which was previously successfully applied for the detection of oscillating magnetic fields in the sub-nT range.

METHODS

The new detection concept was examined by means of Bloch simulations, evaluating the interaction effect of spin-locked magnetization and low-frequency pulsed magnetic fields. The REX detection approach was validated under controlled conditions in phantom experiments at 3 T. Gaussian and sinc-shaped stimuli were investigated. In addition, the detection of artificial fields resembling a cardiac QRS complex, which is the most prominent peak visible on a magnetocardiogram, was tested.

RESULTS

Bloch simulations demonstrated that the REX method has a high sensitivity to pulsed fields in the resonance case, which is met when the spin-lock frequency coincides with a non-zero Fourier component of the stimulus field. In the experiments, we found that magnetic stimuli of different durations and waveforms can be distinguished by their characteristic REX response spectrum. The detected REX amplitude was proportional to the stimulus peak amplitude (R2 > 0.98) and the lowest field detection was 1 nT. Furthermore, the detection of QRS-like fields with varying QRS durations yielded significant results in a phantom setup (p < 0.001).

CONCLUSION

REX detection can be transferred to non-sinusoidal pulsed magnetic fields and could provide a non-invasive, quantitative tool for spatially resolved assessment of cardiac biomagnetism. Potential applications include the direct detection and characterization of cardiac conduction.

Quantification of the rotating frame relaxation time T2ρ: Comparison of balanced spin-lock and continuous-wave Malcolm-Levitt preparations

For the quantification of rotating frame relaxation times, the T2ρ relaxation pathway plays an essential role. Nevertheless, T2ρ imaging has been studied only to a small extent compared with T1ρ, and preparation techniques for T2ρ have so far been adapted from T1ρ methods. In this work, two different preparation concepts are compared specifically for the use of T2ρ mapping. The first approach involves transferring the balanced spin-locking (B-SL) concept of T1ρ imaging. The second and newly proposed approach is a continuous-wave Malcolm-Levitt (CW-MLEV) pulse train with zero echo times and was motivated from T2 preparation strategies. The modules are tested in Bloch simulations for their intrinsic sensitivity to field inhomogeneities and validated in phantom experiments. In addition, myocardial T2ρ mapping was performed in mice as an exemplary application. Our results demonstrate that the CW-MLEV approach provides superior robustness and thus suggest that established methods of T1ρ imaging are not best suited for T2ρ experiments. In the presence of field inhomogeneities, the simulations indicated an increased banding compensation by a factor of 4.1 compared with B-SL. Quantification of left ventricular T2ρ time in mice yielded more consistent results, and values in the range of 59.2-61.1 ms (R2 = 0.986-0.992) were observed at 7 T.

Can MRI Be Used as a Sensor to Record Neural Activity?

Magnetic resonance provides exquisite anatomical images and functional MRI monitors physiological activity by recording blood oxygenation. This review attempts to answer the following question: Can MRI be used as a sensor to directly record neural behavior? It considers MRI sensing of electrical activity in the heart and in peripheral nerves before turning to the central topic: recording of brain activity. The primary hypothesis is that bioelectric current produced by a nerve or muscle creates a magnetic field that influences the magnetic resonance signal, although other mechanisms for detection are also considered. Recent studies have provided evidence that using MRI to sense neural activity is possible under ideal conditions. Whether it can be used routinely to provide functional information about brain processes in people remains an open question. The review concludes with a survey of artificial intelligence techniques that have been applied to functional MRI and may be appropriate for MRI sensing of neural activity.

Analysis of the robustness and dynamics of spin-locking preparations for the detection of oscillatory magnetic fields

Extracting quantitative information of neuronal signals by non-invasive imaging is an outstanding challenge for understanding brain function and pathology. However, state-of-the-art techniques offer low sensitivity to deep electrical sources. Stimulus induced rotary saturation is a recently proposed magnetic resonance imaging sequence that detects oscillatory magnetic fields using a spin-lock preparation. Phantom experiments and simulations proved its efficiency and sensitivity, but the susceptibility of the method to field inhomogeneities is still not well understood. In this study, we simulated and analyzed the dynamic of three spin-lock preparations and their response to field inhomogeneities in the presence of a resonant oscillating field. We show that the composite spin-lock preparation is more robust against field variations within the double resonance effect. In addition, we tested the capability of the chosen composite spin-lock preparation to recover information about the spectral components of a composite signal. This study sets the bases to move one step further towards the clinical application of MR-based neuronal current imaging.

Towards robust in vivo quantification of oscillating biomagnetic fields using Rotary Excitation based MRI

Spin-lock based functional magnetic resonance imaging (fMRI) has the potential for direct spatially-resolved detection of neuronal activity and thus may represent an important step for basic research in neuroscience. In this work, the corresponding fundamental effect of Rotary EXcitation (REX) is investigated both in simulations as well as in phantom and in vivo experiments. An empirical law for predicting optimal spin-lock pulse durations for maximum magnetic field sensitivity was found. Experimental conditions were established that allow robust detection of ultra-weak magnetic field oscillations with simultaneous compensation of static field inhomogeneities. Furthermore, this work presents a novel concept for the emulation of brain activity utilizing the built-in MRI gradient system, which allows REX sequences to be validated in vivo under controlled and reproducible conditions. Via transmission of Rotary EXcitation (tREX), we successfully detected magnetic field oscillations in the lower nano-Tesla range in brain tissue. Moreover, tREX paves the way for the quantification of biomagnetic fields.

Neuronal current imaging: An experimental method to investigate electrical currents in dogs with idiopathic epilepsy

BACKGROUND

The diagnosis of idiopathic epilepsy (IE) in dogs is based on exclusion of other potential causes of seizures. Recently, a novel magnetic resonance imaging (MRI) sequence that utilizes a variant of the rotary saturation approach has been suggested to detect weak transient magnetic field oscillations generated by neuronal currents in humans with epilepsy.

HYPOTHESIS/OBJECTIVES

Effects on the magnetic field evoked by intrinsic epileptic activity can be detected by MRI in the canine brain. As proof-of-concept, the novel MRI sequence to detect neuronal currents was applied in dogs.

ANIMALS

Twelve dogs with IE and 5 control dogs without a history of epileptic seizures were examined.

METHODS

Prospective case-control study as proof-of-concept. All dogs underwent a clinical neurological examination, scalp electroencephalography, cerebrospinal fluid analysis, and MRI. The MRI examination included a spin-locking (SL) experiment applying a low-power on-resonance radiofrequency pulse in a predefined frequency domain in the range of oscillations generated by the epileptogenic tissue.

RESULTS

In 11 of 12 dogs with IE, rotary saturation effects were detected by the MRI sequence. Four of 5 control dogs did not show rotary saturation effects. One control dog with a diagnosis of neuronal ceroid lipofuscinosis had SL-related effects, but did not have epileptic seizures clinically.

CONCLUSIONS AND CLINICAL IMPORTANCE

The proposed MRI method detected neuronal currents in dogs with epileptic seizures and represents a potential new line of research to investigate neuronal currents possibly related to IE in dogs.

Magnetic resonance imaging simulation with spin-lock preparations to detect tiny oscillatory magnetic fields

Spin-lock preparation was studied to detect tiny oscillatory magnetic fields such as a neural magnetic field without the blood oxygen level-dependent effect. This approach is a direct measurement and independent of static magnetic field strength. Accordingly, it is anticipated as a feasible functional magnetic resonance imaging (fMRI) in low and ultra-low-field MRI. Several reports have been published on spin-lock preparation but reports on imaging simulation are rare. Research in this area can assist in investigating magnetic resonance signal changes and, accordingly, can help to develop new spin-lock methods. In this study, we propose an imaging simulation method with an analytical solution using the Bloch equation. To demonstrate the feasibility of our proposed method, we compared simulated images with experimental results in which the number of sub-voxels and the amplitude and phase of the target oscillatory magnetic fields varied. In addition, we also applied graphics processing unit parallel computing and investigated the feasibility of avoiding an impracticable calculation time by doing so.

Dependence of stimulus-induced rotary saturation on the direction of target oscillating magnetic fields: A phantom and simulation study

Several noninvasive techniques for the direct measurement of the neuronal activity using magnetic resonance imaging (MRI) have recently been reported. As a promising candidate, we focus on a spin-lock MRI sequence (i.e., stimulus-induced rotary saturation (SIRS)) directly measuring a tiny oscillating magnetic field. Previous phantom studies on SIRS have applied the target oscillating magnetic field parallel to the direction of the static magnetic field B₀. However, in practice, the neuromagnetic fields are not always aligned in the same direction as in such a condition. This study investigates the MR signal changes during SIRS when the target magnetic field direction is not the same as that of the B₀ field through both phantom experiments and Bloch simulations. The experimental results indicate that only the target magnetic field component along the B₀ field affects the signal change, indicating that SIRS has partial sensitivity, even if the target magnetic fields are tilted from the B₀ field. Furthermore, the simulation results show good agreements with the experimental results. These results clarify the sensitivity direction of SIRS-based fMRI and lead to the possibility that the direction of the generated neuromagnetic fields can be estimated, such that we can separate directional information from the other information contained in neuromagnetic fields (e.g., phase information).

Detection of tiny oscillatory magnetic fields using low-field MRI: A combined phantom and simulation study

We demonstrated the feasibility of the spin-lock preparation sequence using low-field magnetic resonance (MR) imaging that prevents interference from blood-oxygenation-level-dependent effects. We focused on two spin-lock preparations: spin-lock Mz (SL-Mz) and stimulus-induced rotary saturation (SIRS) and analyzed the magnetization dynamics during the sequences using the Bloch equation. Next, we performed phantom experiments using a loop coil to investigate the MR signal change as a function of the target signal strength and phase. Furthermore, we performed curve fittings to consider the radio frequency, which agreed with the experimental results. Then, we investigated the detectable strength of the magnetic field, and the SL-Mz detected a signal strength of 2.34 nT. In conclusion, our experimental results showed good agreement with the results obtained using the Bloch equation.

Detection of nanotesla AC magnetic fields using steady-state SIRS and ultra-low field MRI

OBJECTIVE

Functional magnetic resonance imaging (fMRI) is commonly used to measure brain activity through the blood oxygen level dependent (BOLD) signal mechanism, but this only provides an indirect proxy signal to neuronal activity. Magnetoencephalography (MEG) provides a more direct measurement of the magnetic fields created by neuronal currents in the brain, but requires very specialized hardware and only measures these fields at the scalp. Recently, progress has been made to directly detect neuronal fields with MRI using the stimulus-induced rotary saturation (SIRS) effect, but interference from the BOLD response complicates such measurements. Here, we describe an approach to detect nanotesla-level, low-frequency alternating magnetic fields with an ultra-low field (ULF) MRI scanner, unaffected by the BOLD signal.

APPROACH

A steady-state implementation of the stimulus-induced rotary saturation (SIRS) method is developed. The method is designed to generate a strong signal at ultra-low magnetic field as well as allowing for efficient signal averaging, giving a high contrast-to-noise ratio (CNR). The method is tested in computer simulations and in phantom scans.

MAIN RESULTS

The simulations and phantom scans demonstrated the ability of the method to measure magnetic fields at different frequencies at ULF with a stronger contrast than non-steady-state approaches. Furthermore, the rapid imaging functionality of the method reduced noise efficiently. The results demonstrated sufficient CNR down to 7 nT, but the sensitivity will depend on the imaging parameters.

SIGNIFICANCE

A steady-state SIRS method is able to detect low-frequency alternating magnetic fields at ultra-low main magnetic field strengths with a large signal response and contrast-to-noise, presenting an important step in sensing biological fields with ULF MRI.

Neural magnetic field dependent fMRI toward direct functional connectivity measurements: A phantom study

Recently, the main issue in neuroscience has been the imaging of the functional connectivity in the brain. No modality that can measure functional connectivity directly, however, has been developed yet. Here, we show the novel MRI sequence, called the partial spinlock sequence toward direct measurements of functional connectivity. This study investigates a probable measurement of phase differences directly associated with functional connectivity. By employing partial spinlock imaging, the neural magnetic field might influence the magnetic resonance signals. Using simulation and phantom studies to model the neural magnetic fields, we showed that magnetic resonance signals vary depending on the phase of an externally applied oscillating magnetic field with non-right flip angles. These results suggest that the partial spinlock sequence is a promising modality for functional connectivity measurements.

Toward direct MRI of neuro-electro-magnetic oscillations in the human brain

PURPOSE

Neuroimaging techniques are widely used to investigate the function of the human brain, but none are currently able to accurately localize neuronal activity with both high spatial and temporal specificity. Here, a new in vivo MRI acquisition and analysis technique based on the spin-lock mechanism is developed to noninvasively image local magnetic field oscillations resulting from neuroelectric activity in specifiable frequency bands.

METHODS

Simulations, phantom experiments, and in vivo experiments using an eyes-open/eyes-closed task in 8 healthy volunteers were performed to demonstrate its sensitivity and specificity for detecting oscillatory neuroelectric activity in the alpha-band (8-12 Hz). A comprehensive postprocessing procedure was designed to enhance the neuroelectric signal, while minimizing any residual hemodynamic and physiological confounds.

RESULTS

The phantom results show that this technique can detect 0.06-nT magnetic field oscillations, while the in vivo results demonstrate that it can image task-based modulations of neuroelectric oscillatory activity in the alpha-band. Multiple control experiments and a comparison with conventional BOLD functional MRI suggest that the activation was likely not due to any residual hemodynamic or physiological confounds.

CONCLUSION

These initial results provide evidence suggesting that this new technique has the potential to noninvasively and directly image neuroelectric activity in the human brain in vivo. With further development, this approach offers the promise of being able to do so with a combination of spatial and temporal specificity that is beyond what can be achieved with existing neuroimaging methods, which can advance our ability to study the functions and dysfunctions of the human brain.

Dynamics of magnetization under stimulus-induced rotary saturation sequence

We studied stimulus-induced rotary-saturation preparation (which enables measurement of oscillating magnetic fields using MRI) and derived an analytical solution of the Bloch equation to understand magnetization dynamics mathematically and comprehensively and to conduct simulations without sequential-calculation techniques such as the Runge-Kutta method. We formulated the dynamics using the Bloch equation, introducing an additional rotating frame and some approximations to make it into a homogeneous differential equation. Moreover, we found that there are two modes depending on the target oscillating magnetic field. To confirm the validity of the solution, we experimentally investigated its characteristics and performed curve fitting using the analytical model. Considering the constraints on the frame, the analytical solution was found to agree with experimental data. The experimental data indicate that it is necessary to design robust sequences compensating B₀ or B_1lock spatial inhomogeneity to improve measurements. Therefore, experimenters should consider the dynamics of magnetization with RF pulses to rewind the spin phase for accurate measurements.

Magnetic Resonance Imaging technology-bridging the gap between noninvasive human imaging and optical microscopy

Technological advances in Magnetic Resonance Imaging (MRI) have provided substantial gains in the sensitivity and specificity of functional neuroimaging. Mounting evidence demonstrates that the hemodynamic changes utilized in functional MRI can be far more spatially and thus neuronally specific than previously believed. This has motivated a push toward novel, high-resolution MR imaging strategies that can match this biological resolution limit while recording from the entire human brain. Although sensitivity increases are a necessary component, new MR encoding technologies are required to convert improved sensitivity into higher resolution. These new sampling strategies improve image acquisition efficiency and enable increased image encoding in the time-frame needed to follow hemodynamic changes associated with brain activation.

Detection of fast oscillating magnetic fields using dynamic multiple TR imaging and Fourier analysis

Neuronal oscillations produce oscillating magnetic fields. There have been trials to detect neuronal oscillations using MRI, but the detectability in in vivo is still in debate. Major obstacles to detecting neuronal oscillations are (i) weak amplitudes, (ii) fast oscillations, which are faster than MRI temporal resolution, and (iii) random frequencies and on/off intervals. In this study, we proposed a new approach for direct detection of weak and fast oscillating magnetic fields. The approach consists of (i) dynamic acquisitions using multiple times to repeats (TRs) and (ii) an expanded frequency spectral analysis. Gradient echo echo-planar imaging was used to test the feasibility of the proposed approach with a phantom generating oscillating magnetic fields with various frequencies and amplitudes and random on/off intervals. The results showed that the proposed approach could precisely detect the weak and fast oscillating magnetic fields with random frequencies and on/off intervals. Complex and phase spectra showed reliable signals, while no meaningful signals were observed in magnitude spectra. A two-TR approach provided an absolute frequency spectrum above Nyquist sampling frequency pixel by pixel with no a priori target frequency information. The proposed dynamic multiple-TR imaging and Fourier analysis are promising for direct detection of neuronal oscillations and potentially applicable to any pulse sequences.

Relationship altered between functional T1ρ and BOLD signals in bipolar disorder

INTRODUCTION

Functional neuroimaging typically relies on the blood-oxygen-level-dependent (BOLD) contrast, which is sensitive to the influx of oxygenated blood following neuronal activity. A new method, functional T1 relaxation in the rotating frame (fT1ρ) is thought to reflect changes in local brain metabolism, likely pH, and may more directly measure neuronal activity. These two methods were applied to study activation of the visual cortex in participants with bipolar disorder as compared to controls.

METHODS

Thirty-nine participants with bipolar disorder and 32 healthy controls underwent functional neuroimaging during a flashing checkerboard paradigm. Functional images were acquired in alternating blocks of BOLD and fT1ρ. Linear mixed-effect models were used to examine the relationship between these two functional imaging modalities and to test whether that relationship was altered in bipolar disorder.

RESULTS

BOLD and fT1ρ signal were strongly related in visual and cerebellar areas during the task in controls. The relationship between these two measures was reduced in bipolar disorder within the visual areas, cerebellum, striatum, and thalamus.

CONCLUSIONS

These results support a distinct mechanisms underlying BOLD and fT1ρ signals. The weakened relationship between these imaging modalities may provide a novel tool for measuring pathology in bipolar disorder and other psychiatric illnesses.

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.

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.

MR imaging of oscillatory magnetic field changes: Progressing from phantom to human

Detection of ultra-weak oscillatory magnetic field changes using MRI is of great research interest not only for neuronal current MRI of endogenous neuronal oscillations but also for direct visualization of exogenous transcranial currents or iron oxide contrast agent distribution. In this work, we present a novel oscillatory-selective detection (OSD) method that is magnitude-sensitive to the oscillatory magnetic field changes and immune to the main field inhomogeneity. In OSD, a train of 180° pulses with alternating polarity and mirror symmetry are used to refocus and accumulate magnetization changes induced by external oscillatory fields. After taking both the signal change and image signal-to-noise ratio (SNR) into account, a final 90° pulse with a phase offset of 45° is applied to store a combination of the current-induced signal change and background magnetization for the subsequent EPI acquisition. Its performance was demonstrated in phantom and human studies, both of which showed much better detection in the comparison with the recently proposed spin-lock oscillatory excitation (SLOE) method. OSD was further successfully applied in imaging transcranial alternating current stimulation (tACS) induced field changes in the human brain. These promising results suggest that OSD can overcome the limitation of field inhomogeneity impeding previous oscillatory current MRI sensitivity and be a viable tool in future tACS study.

Direct neural current imaging in an intact cerebellum with magnetic resonance imaging

The ability to detect neuronal currents with high spatiotemporal resolution using magnetic resonance imaging (MRI) is important for studying human brain function in both health and disease. While significant progress has been made, we still lack evidence showing that it is possible to measure an MR signal time-locked to neuronal currents with a temporal waveform matching concurrently recorded local field potentials (LFPs). Also lacking is evidence that such MR data can be used to image current distribution in active tissue. Since these two results are lacking even in vitro, we obtained these data in an intact isolated whole cerebellum of turtle during slow neuronal activity mediated by metabotropic glutamate receptors using a gradient-echo EPI sequence (TR=100ms) at 4.7T. Our results show that it is possible (1) to reliably detect an MR phase shift time course matching that of the concurrently measured LFP evoked by stimulation of a cerebellar peduncle, (2) to detect the signal in single voxels of 0.1mm(3), (3) to determine the spatial phase map matching the magnetic field distribution predicted by the LFP map, (4) to estimate the distribution of neuronal current in the active tissue from a group-average phase map, and (5) to provide a quantitatively accurate theoretical account of the measured phase shifts. The peak values of the detected MR phase shifts were 0.27-0.37°, corresponding to local magnetic field changes of 0.67-0.93nT (for TE=26ms). Our work provides an empirical basis for future extensions to in vivo imaging of neuronal currents.

Focal Epilepsy: MR Imaging of Nonhemodynamic Field Effects by Using a Phase-cycled Stimulus-induced Rotary Saturation Approach with Spin-Lock Preparation

Purpose To investigate whether nonhemodynamic resonant saturation effects can be detected in patients with focal epilepsy by using a phase-cycled stimulus-induced rotary saturation (PC-SIRS) approach with spin-lock (SL) preparation and whether they colocalize with the seizure onset zone and surface interictal epileptiform discharges (IED). Materials and Methods The study was approved by the local ethics committee, and all subjects gave written informed consent. Eight patients with focal epilepsy undergoing presurgical surface and intracranial electroencephalography (EEG) underwent magnetic resonance (MR) imaging at 3 T with a whole-brain PC-SIRS imaging sequence with alternating SL-on and SL-off and two-dimensional echo-planar readout. The power of the SL radiofrequency pulse was set to 120 Hz to sensitize the sequence to high gamma oscillations present in epileptogenic tissue. Phase cycling was applied to capture distributed current orientations. Voxel-wise subtraction of SL-off from SL-on images enabled the separation of T2* effects from rotary saturation effects. The topography of PC-SIRS effects was compared with the seizure onset zone at intracranial EEG and with surface IED-related potentials. Bayesian statistics were used to test whether prior PC-SIRS information could improve IED source reconstruction. Results Nonhemodynamic resonant saturation effects ipsilateral to the seizure onset zone were detected in six of eight patients (concordance rate, 0.75; 95% confidence interval: 0.40, 0.94) by means of the PC-SIRS technique. They were concordant with IED surface negativity in seven of eight patients (0.88; 95% confidence interval: 0.51, 1.00). Including PC-SIRS as prior information improved the evidence of the standard EEG source models compared with the use of uninformed reconstructions (exceedance probability, 0.77 vs 0.12; Wilcoxon test of model evidence, P < .05). Nonhemodynamic resonant saturation effects resolved in patients with favorable postsurgical outcomes, but persisted in patients with postsurgical seizure recurrence. Conclusion Nonhemodynamic resonant saturation effects are detectable during interictal periods with the PC-SIRS approach in patients with epilepsy. The method may be useful for MR imaging-based detection of neuronal currents in a clinical environment. (©) RSNA, 2016 Online supplemental material is available for this article.

A comprehensive study of sensitivity in measuring oscillatory magnetic fields using rotary saturation pulse sequences

Detecting the oscillatory currents with a specific frequency distribution may have the potential to make neuronal current MRI (ncMRI) come true. The phase shift or dephasing induced by both positive and negative episodes of oscillatory neuronal currents is likely to be canceled out over the echo time in typical BOLD-contrast fMRI experiments. Based on the contrast of rotary saturation, both of the recently developed spin-locked oscillatory excitation (SLOE) and stimulus-induced rotary saturation (SIRS) pulse sequences have been demonstrated to be able to detect weak oscillatory magnetic fields in phantoms with 3T MR scanners. In this report, through Bloch equation simulation as well as water phantom and anesthetic rats experiments, we comprehensively evaluate and compare the sensitivities of these two methods (SLOE and SIRS) in detecting the oscillatory magnetic fields for both high (100 Hz) and low (10 Hz) oscillation frequencies, while using their respective optimal imaging parameters. In agreement with the theoretical predications, both the simulated and experimental results showed that the SLOE method features a much higher detection sensitivity of weak magnetic fields than that of the SIRS method. SLOE was able to detect applied oscillatory magnetic fields as low as 0.1 nT in a water phantom and 0.5 nT in rat brains and the deteriorated noise levels in rat data may account for the reduced sensitivity in vivo. These promising results form the foundation for direct detection of in vivo neuronal currents using MRI.

Simultaneous calcium fluorescence imaging and MR of ex vivo organotypic cortical cultures: a new test bed for functional MRI

Recently, several new functional (f)MRI contrast mechanisms including diffusion, phase imaging, proton density, etc. have been proposed to measure neuronal activity more directly and accurately than blood-oxygen-level dependent (BOLD) fMRI. However, these approaches have proved difficult to reproduce, mainly because of the dearth of reliable and robust test systems to vet and validate them. Here we describe the development and testing of such a test bed for non-BOLD fMRI. Organotypic cortical cultures were used as a stable and reproducible biological model of neuronal activity that shows spontaneous activity similar to that of in vivo brain cortex without any hemodynamic confounds. An open-access, single-sided magnetic resonance (MR) "profiler" consisting of four permanent magnets with magnetic field of 0.32 T was used in this study to perform MR acquisition. A fluorescence microscope with long working distance objective was mounted on the top of a custom-designed chamber that keeps the organotypic culture vital, and the MR system was mounted on the bottom of the chamber to achieve real-time simultaneous calcium fluorescence optical imaging and MR acquisition on the same specimen. In this study, the reliability and performance of the proposed test bed were demonstrated by a conventional CPMG MR sequence acquired simultaneously with calcium imaging, which is a well-characterized measurement of neuronal activity. This experimental design will make it possible to correlate directly the other candidate functional MR signals to the optical indicia of neuronal activity in the future.

Direct detection of optogenetically evoked oscillatory neuronal electrical activity in rats using SLOE sequence

The direct detection of neuronal electrical activity is one of the most challenging goals in non-BOLD fMRI research. Previous work has demonstrated its feasibility in phantom and cell culture studies, but attempts in in vivo studies remain few and far between. Most recent in vivo studies used T2*-weighted sequences to directly detect neuronal electrical activity evoked by sensory stimulus. As neuronal electrical signal is usually comprised of a series of spectrally distributed oscillatory waveforms rather than being a direct current, it is most likely to be detected using oscillatory current sensitive sequences. In this study, we explored the potential of using the spin-lock oscillatory excitation (SLOE) sequence with spiral readout to directly detect optogenetically evoked oscillatory neuronal electrical activity, whose main spectral component can be manipulated artificially to match the resonance frequency of spin-lock RF field. In addition, experiments using the stimulus-induced rotary saturation (SIRS) sequence with spiral readout were also performed. Electrophysiological recording and MRI data acquisition were conducted on separate animals. Robust optogenetically evoked oscillatory LFP signals were observed and significant BOLD signals were acquired with the GE-EPI sequence before and after the whole SLOE and SIRS acquisitions, but no significant neuronal current MRI (ncMRI) signal changes were detected. These results indicate that the sensitivity of oscillatory current sensitive sequences needs to be further improved for direct detection of neuronal electrical activity.

Detection of subnanotesla oscillatory magnetic fields using MRI

PURPOSE

Direct mapping of neuronal currents using MRI would have fundamental impacts on brain functional imaging. Previous reports indicated that the stimulus-induced rotary saturation (SIRS) mechanism had the best potential of direct detection of neural oscillations; however, it lacked the high-sensitivity level needed. In this study, a novel strategy is proposed in an effort to improve the detection sensitivity.

METHODS

In our modified SIRS sequence, an external oscillatory magnetic field is used as the excitation pulse in place of the standard 90-degree excitation pulse. This approach could potentially lead to tens or even hundreds times of enhancement in the detection sensitivity for low field signals. It also helps to lower the physiological noise, allows for shorter pulse repetition time, and is less affected by the blood oxygen level.

RESULTS

We demonstrate that a 100-Hz oscillatory magnetic field with magnitude as low as 0.25 nanotesla generated in a current loop can be robustly detected using a 3-Tesla MRI scanner.

CONCLUSION

The modified SIRS sequence offers higher detection sensitivity as well as several additional advantages. These promising results suggest that the direct detection of neural oscillation might be within the grasp of the current MRI technology.

Selective magnetic resonance imaging of magnetic nanoparticles by acoustically induced rotary saturation

PURPOSE

The goal of this study was to introduce a new method to selectively detect iron oxide contrast agents using an acoustic wave to perturb the spin-locked water signal in the vicinity of the magnetic particles. The acoustic drive can be modulated externally to turn the effect on and off, allowing sensitive and quantitative statistical comparison and removal of confounding image background variations.

METHODS

We demonstrated the effect in spin-locking experiments using piezoelectric actuators to generate vibrational displacements of iron oxide samples. We observed a resonant behavior of the signal changes with respect to the acoustic frequency where iron oxide is present. We characterized the effect as a function of actuator displacement and contrast agent concentration.

RESULTS

The resonant effect allowed us to generate block-design "modulation response maps" indicating the contrast agent's location, as well as positive contrast images with suppressed background signal. We found that the acoustically induced rotary saturation (AIRS) effect stayed approximately constant across acoustic frequency and behaved monotonically over actuator displacement and contrast agent concentration.

CONCLUSION

AIRS is a promising method capable of using acoustic vibrations to modulate the contrast from iron oxide nanoparticles and thus perform selective detection of the contrast agents, potentially enabling more accurate visualization of contrast agents in clinical and research settings.

Magnetic resonance imaging of ionic currents in solution: the effect of magnetohydrodynamic flow

PURPOSE

Reliably detecting MRI signals in the brain that are more tightly coupled to neural activity than blood-oxygen-level-dependent fMRI signals could not only prove valuable for basic scientific research but could also enhance clinical applications such as epilepsy presurgical mapping. This endeavor will likely benefit from an improved understanding of the behavior of ionic currents, the mediators of neural activity, in the presence of the strong magnetic fields that are typical of modern-day MRI scanners.

THEORY

Of the various mechanisms that have been proposed to explain the behavior of ionic volume currents in a magnetic field, only one-magnetohydrodynamic flow-predicts a slow evolution of signals, on the order of a minute for normal saline in a typical MRI scanner.

METHODS

This prediction was tested by scanning a volume-current phantom containing normal saline with gradient-echo-planar imaging at 3 T.

RESULTS

Greater signal changes were observed in the phase of the images than in the magnitude, with the changes evolving on the order of a minute.

CONCLUSION

These results provide experimental support for the MHD flow hypothesis. Furthermore, MHD-driven cerebrospinal fluid flow could provide a novel fMRI contrast mechanism.

Physical principles for scalable neural recording

Simultaneously measuring the activities of all neurons in a mammalian brain at millisecond resolution is a challenge beyond the limits of existing techniques in neuroscience. Entirely new approaches may be required, motivating an analysis of the fundamental physical constraints on the problem. We outline the physical principles governing brain activity mapping using optical, electrical, magnetic resonance, and molecular modalities of neural recording. Focusing on the mouse brain, we analyze the scalability of each method, concentrating on the limitations imposed by spatiotemporal resolution, energy dissipation, and volume displacement. Based on this analysis, all existing approaches require orders of magnitude improvement in key parameters. Electrical recording is limited by the low multiplexing capacity of electrodes and their lack of intrinsic spatial resolution, optical methods are constrained by the scattering of visible light in brain tissue, magnetic resonance is hindered by the diffusion and relaxation timescales of water protons, and the implementation of molecular recording is complicated by the stochastic kinetics of enzymes. Understanding the physical limits of brain activity mapping may provide insight into opportunities for novel solutions. For example, unconventional methods for delivering electrodes may enable unprecedented numbers of recording sites, embedded optical devices could allow optical detectors to be placed within a few scattering lengths of the measured neurons, and new classes of molecularly engineered sensors might obviate cumbersome hardware architectures. We also study the physics of powering and communicating with microscale devices embedded in brain tissue and find that, while radio-frequency electromagnetic data transmission suffers from a severe power-bandwidth tradeoff, communication via infrared light or ultrasound may allow high data rates due to the possibility of spatial multiplexing. The use of embedded local recording and wireless data transmission would only be viable, however, given major improvements to the power efficiency of microelectronic devices.

Bloch simulations towards direct detection of oscillating magnetic fields using MRI with spin-lock sequence

A new MRI method using the spin-lock sequence has attracted wide attention because of its potential for detecting small oscillating magnetic fields. However, as the mechanism involved is complicated, we visualized the magnetization performance during the spin-lock sequence in order to better understand interaction of the spin-lock pulse and the externally applied oscillating magnetic fields by means of a fast-and-simple method using matrix operations to solve a time-dependent Bloch equation. To improve spin-lock imaging in the detection of small magnetic fields (in an fMRI experiment that modeled neural magnetic fields), we observed that the phenomenon decreases MR signals, which led us to investigate how spin-lock parameters cause the MR signal to decrease; based on this, we determined that MR signals decrease in oscillating magnetic fields that are resonant with the spin-lock pulse. We also determined that MR signals decrease is directly proportional to spin-lock duration. Our results suggest that MRI can feasibly detect oscillating magnetic fields directly by using of the spin-lock sequence.

Fast human brain magnetic resonance responses associated with epileptiform spikes

Neuronal currents produce local electromagnetic fields that can potentially modulate the phase of the magnetic resonance signal and thus provide a contrast mechanism tightly linked to neuronal activity. Previous work has demonstrated the feasibility of direct MRI of neuronal activity in phantoms and cell culture, but in vivo efforts have yielded inconclusive, conflicting results. The likelihood of detecting and validating such signals can be increased with (i) fast gradient-echo echo-planar imaging, with acquisition rates sufficient to resolve neuronal activity, (ii) subjects with epilepsy, who frequently experience stereotypical electromagnetic discharges between seizures, expressed as brief, localized, high-amplitude spikes (interictal discharges), and (iii) concurrent electroencephalography. This work demonstrates that both MR magnitude and phase show large-amplitude changes concurrent with electroencephalography spikes. We found a temporal derivative relationship between MR phase and scalp electroencephalography, suggesting that the MR phase changes may be tightly linked to local cerebral activity. We refer to this manner of MR acquisition, designed explicitly to track the electroencephalography, as encephalographic MRI (eMRI). Potential extension of this technique into a general purpose functional neuroimaging tool requires further study of the MR signal changes accompanying lower amplitude neuronal activity than those discussed here.

Magnetic resonance imaging of oscillating electrical currents

Functional MRI has become an important tool of researchers and clinicians who seek to understand patterns of neuronal activation that accompany sensory and cognitive processes. However, the interpretation of fMRI images rests on assumptions about the relationship between neuronal firing and hemodynamic response that are not firmly grounded in rigorous theory or experimental evidence. Further, the blood-oxygen-level-dependent effect, which correlates an MRI observable to neuronal firing, evolves over a period that is 2 orders of magnitude longer than the underlying processes that are thought to cause it. Here, we instead demonstrate experiments to directly image oscillating currents by MRI. The approach rests on a resonant interaction between an applied rf field and an oscillating magnetic field in the sample and, as such, permits quantitative, frequency-selective measurements of current density without spatial or temporal cancellation. We apply this method in a current loop phantom, mapping its magnetic field and achieving a detection sensitivity near the threshold required for the detection of neuronal currents. Because the contrast mechanism is under spectroscopic control, we are able to demonstrate how ramped and phase-modulated spin-lock radiation can enhance the sensitivity and robustness of the experiment. We further demonstrate the combination of these methods with remote detection, a technique in which the encoding and detection of an MRI experiment are separated by sample flow or translation. We illustrate that remotely detected MRI permits the measurement of currents in small volumes of flowing water with high sensitivity and spatial resolution.

Physiologically evoked neuronal current MRI in a bloodless turtle brain: detectable or not?

Contradictory reports regarding the detection of neuronal currents have left the feasibility of neuronal current MRI (ncMRI) an open question. Most previous ncMRI studies in human subjects are suspect due to their inability to separate or eliminate hemodynamic effects. In this study, we used a bloodless turtle brain preparation that eliminates hemodynamic effects, to explore the feasibility of detecting visually-evoked ncMRI signals at 9.4 T. Intact turtle brains, with eyes attached, were dissected from the cranium and placed in artificial cerebral spinal fluid. Light flashes were delivered to the eyes to evoke neuronal activity. Local field potential (LFP) and MRI signals were measured in an interleaved fashion. Robust visually-evoked LFP signals were observed in turtle brains, but no significant signal changes synchronized with neuronal currents were found in the ncMRI images. In this study, detection thresholds of 0.1% and 0.1 degrees were set for MRI magnitude and phase signal changes, respectively. The absence of significant signal changes in the MRI images suggests that visually-evoked ncMRI signals in the turtle brain are below these detectable levels.