Are brain currents detectable by means of low-field NMR? A phantom study

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.

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.

Head phantoms for electroencephalography and transcranial electric stimulation: a skull material study

Physical head phantoms allow the assessment of source reconstruction procedures in electroencephalography and electrical stimulation profiles during transcranial electric stimulation. Volume conduction in the head is strongly influenced by the skull, which represents the main conductivity barrier. Realistic modeling of its characteristics is thus important for phantom development. In the present study, we proposed plastic clay as a material for modeling the skull in phantoms. We analyzed five clay types varying in granularity and fractions of fire clay, each with firing temperatures from 550°C to 950°C. We investigated the conductivity of standardized clay samples when immersed in a 0.9% sodium chloride solution with time-resolved four-point impedance measurements. To test the reusability of the clay model, these measurements were repeated after cleaning the samples by rinsing in deionized water for 5 h. We found time-dependent impedance changes for approximately 5 min after immersion in the solution. Thereafter, the conductivities stabilized between 0.0716 S/m and 0.0224 S/m depending on clay type and firing temperatures. The reproducibility of the measurement results proved the effectiveness of the rinsing procedure. Clay provides formability, is permeable to ions, can be adjusted in conductivity value and is thus suitable for the skull modeling in phantoms.

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.

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.

Current-density imaging using ultra-low-field MRI with zero-field encoding

Electric current density can be measured noninvasively with magnetic resonance imaging (MRI). Determining all three components of the current density, however, requires physical rotation of the sample or current injection from several directions when done with conventional methods. However, the emerging technology of ultra-low-field (ULF) MRI, in which the signal encoding and acquisition is conducted at a microtesla-range magnetic field, offers new possibilities. The low applied magnetic fields can even be switched off completely within the pulse sequence, increasing the flexibility of the available sequences. In this article, we present a ULF-MRI sequence designed for obtaining all three components of a current-density pattern without the need of sample rotations. The sequence consists of three steps: prepolarization of the sample, signal encoding in the current-density-associated magnetic field without applying any MRI fields, and spatial encoding in a microtesla-range field using any standard ULF-MRI sequence. The performance of the method is evaluated by numerical simulations. The method may find applications, e.g., in noninvasive conductivity imaging of tissue.

Toward a brain functional connectivity mapping modality by simultaneous imaging of coherent brainwaves

Matching the proton-magnetic-resonance frequency to the frequency of a periodic neural oscillation (e.g., alpha or gamma band waves) by magnetic resonance imaging techniques, enables direct visualization of brain functional connectivity. Functional connectivity has been studied by analyzing the correlation between coherent neural oscillations in different areas of the brain. In electro- or magneto-encephalography, coherent source reconstruction in a source-space is very tricky due to power leaking from the correlation among the sources. For this reason, most studies have been limited to sensor-space analyses, which give doubtful results because of volume current mixing. The direct visualization of coherent brain oscillations can circumvent this problem. The feasibility of this idea was demonstrated by conducting phantom experiments with a SQUID-based, micro-Tesla NMR/MRI system. We introduce an experimental trick, an effective step-up of the measurement B-field in a pulse sequence, to decouple the magnetic resonance signal from the strong magneto-encephalographic signal at the same frequency.

An advanced phantom study assessing the feasibility of neuronal current imaging by ultra-low-field NMR

In ultra-low-field (ULF) NMR/MRI, a common scheme is to magnetize the sample by a polarizing field of up to hundreds of mT, after which the NMR signal, precessing in a field on the order of several μT, is detected with superconducting quantum interference devices (SQUIDs). In our ULF-NMR system, we polarize with up to 50mT and deploy a single-stage DC-SQUID current sensor with an integrated input coil which is connected to a wire-wound Nb gradiometer. We developed this system (white noise 0.50fT/√Hz) for assessing the feasibility of imaging neuronal currents by detecting their effect on the ULF-NMR signal. Magnetoencephalography investigations of evoked brain activity showed neuronal dipole moments below 50nAm. With our instrumentation, we have studied two different approaches for neuronal current imaging. In the so-called DC effect, long-lived neuronal activity shifts the Larmor frequency of the surrounding protons. An alternative strategy is to exploit fast neuronal activity as a tipping pulse. This so-called AC effect requires the proton Larmor frequency to match the frequency of the neuronal activity, which ranges from near-DC to ∼kHz. We emulated neuronal activity by means of a single dipolar source in a physical phantom, consisting of a hollow sphere filled with an aqueous solution of CuSO4 and NaCl. In these phantom studies, with physiologically relevant dipole depths, we determined resolution limits for our set-up for the AC and the DC effect of ∼10μAm and ∼50nAm, respectively. Hence, the DC effect appears to be detectable in vivo by current ULF-NMR technology.

Current-density imaging using ultra-low-field MRI with adiabatic pulses

Magnetic resonance imaging (MRI) allows measurement of electric current density in an object. The measurement is based on observing how the magnetic field of the current density affects the associated spins. However, as high-field MRI is sensitive to static magnetic field variations of only the field component along the main field direction, object rotations are typically needed to image three-dimensional current densities. Ultra-low-field (ULF) MRI, on the other hand, with B0 on the order of 10-100 μT, allows novel MRI sequences. We present a rotation-free method for imaging static magnetic fields and current densities using ULF MRI. The method utilizes prepolarization pulses with adiabatic switch-off ramps. The technique is designed to reveal complete field and current-density information without the need to rotate the object. The method may find applications, e.g., in conductivity imaging. We present simulation results showing the feasibility of the sequence.

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.

Imaging artifacts induced by electrical stimulation during conventional fMRI of the brain

Functional magnetic resonance imaging (fMRI) of brain activation during transcranial electrical stimulation is used to provide insight into the mechanisms of neuromodulation and targeting of particular brain structures. However, the passage of current through the body may interfere with the concurrent detection of blood oxygen level-dependent (BOLD) signal, which is sensitive to local magnetic fields. To test whether these currents can affect concurrent fMRI recordings we performed conventional gradient echo-planar imaging (EPI) during transcranial direct current (tDCS) and alternating current stimulation (tACS) on two post-mortem subjects. tDCS induced signals in both superficial and deep structures. The signal was specific to the electrode montage, with the strongest signal near cerebrospinal fluid (CSF) and scalp. The direction of change relative to non-stimulation reversed with tDCS stimulation polarity. For tACS there was no net effect of the MRI signal. High-resolution individualized modeling of current flow and induced static magnetic fields suggested a strong coincidence of the change EPI signal with regions of large current density and magnetic fields. These initial results indicate that (1) fMRI studies of tDCS must consider this potentially confounding interference from current flow and (2) conventional MRI imaging protocols can be potentially used to measure current flow during transcranial electrical stimulation. The optimization of current measurement and artifact correction techniques, including consideration of the underlying physics, remains to be addressed.

Magnetic resonance imaging at frequencies below 1 kHz

Within the magnetic resonance imaging (MRI) community the trend is going to higher and higher magnetic fields, ranging from 1.5 T to 7 T, corresponding to Larmor frequencies of 63.8-298 MHz. Since for high-field MRI the magnetization increases with the applied magnetic field, the signal-to-noise-ratio increases as well, thus enabling higher image resolutions. On the other hand, MRI is possible also at ultra-low magnetic fields, as was shown by different groups. The goal of our development was to reach a Larmor frequency range of the low-field MRI system corresponding to the frequency range of human brain activities ranging from near zero-frequency (near-DC) to over 1 kHz. Here, first 2D MRI images of phantoms taken at Larmor frequencies of 100 Hz and 731 Hz will be shown and discussed. These frequencies are examples of brain activity triggered by electrostimulation of the median nerve. The method will allow the magnetic fields of the brain currents to influence the magnetic resonance image, and thus lead to a direct functional imaging modality of neuronal currents.