Neuronal current detection with low-field magnetic resonance: simulations and methods

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.

Feasibility of functional MRI at ultralow magnetic field via changes in cerebral blood volume

We investigate the feasibility of performing functional MRI (fMRI) at ultralow field (ULF) with a Superconducting QUantum Interference Device (SQUID), as used for detecting magnetoencephalography (MEG) signals from the human head. While there is negligible magnetic susceptibility variation to produce blood oxygenation level-dependent (BOLD) contrast at ULF, changes in cerebral blood volume (CBV) may be a sensitive mechanism for fMRI given the five-fold spread in spin-lattice relaxation time (T₁) values across the constituents of the human brain. We undertook simulations of functional signal strength for a simplified brain model involving activation of a primary cortical region in a manner consistent with a blocked task experiment. Our simulations involve measured values of T₁ at ULF and experimental parameters for the performance of an upgraded ULFMRI scanner. Under ideal experimental conditions we predict a functional signal-to-noise ratio of between 3.1 and 7.1 for an imaging time of 30 min, or between 1.5 and 3.5 for a blocked task experiment lasting 7.5 min. Our simulations suggest it may be feasible to perform fMRI using a ULFMRI system designed to perform MRI and MEG in situ.

Adaptive suppression of power line interference in ultra-low field magnetic resonance imaging in an unshielded environment

Power-line harmonic interference and fixed-frequency noise peaks may cause stripe-artifacts in ultra-low field (ULF) magnetic resonance imaging (MRI) in an unshielded environment and in a conductively shielded room. In this paper we describe an adaptive suppression method to eliminate these artifacts in MRI images. This technique utilizes spatial correlation of the interference from different positions, and is realized by subtracting the outputs of the reference channel(s) from those of the signal channel(s) using wavelet analysis and the least squares method. The adaptive suppression method is first implemented to remove the image artifacts in simulation. We then experimentally demonstrate the feasibility of this technique by adding three orthogonal superconducting quantum interference device (SQUID) magnetometers as reference channels to compensate the output of one 2nd-order gradiometer. The experimental results show great improvement in the imaging quality in both 1D and 2D MRI images at two common imaging frequencies, 1.3 kHz and 4.8 kHz. At both frequencies, the effective compensation bandwidth is as high as 2 kHz. Furthermore, we examine the longitudinal relaxation times of the same sample before and after compensation, and show that the MRI properties of the sample did not change after applying adaptive suppression. This technique can effectively increase the imaging bandwidth and be applied to ULF MRI detected by either SQUIDs or Faraday coil in both an unshielded environment and a conductively shielded room.

Modeling the effect of dendritic input location on MEG and EEG source dipoles

The cerebral sources of magneto- and electroencephalography (MEG, EEG) signals can be represented by current dipoles. We used computational modeling of realistically shaped passive-membrane dendritic trees of pyramidal cells from the human cerebral cortex to examine how the spatial distribution of the synaptic inputs affects the current dipole. The magnitude of the total dipole moment vector was found to be proportional to the vertical location of the synaptic input. The dipole moment had opposite directions for inputs above and below a reversal point located near the soma. Inclusion of shunting-type inhibition either suppressed or enhanced the current dipole, depending on whether the excitatory and inhibitory synapses were on the same or opposite side of the reversal point. Relating the properties of the macroscopic current dipoles to dendritic current distributions can help to provide means for interpreting MEG and EEG data in terms of synaptic connection patterns within cortical areas.

Current density imaging sequence for monitoring current distribution during delivery of electric pulses in irreversible electroporation

Background

Electroporation is gaining its importance in everyday clinical practice of cancer treatment. For its success it is extremely important that coverage of the target tissue, i.e. treated tumor, with electric field is within the specified range. Therefore, an efficient tool for the electric field monitoring in the tumor during delivery of electroporation pulses is needed. The electric field can be reconstructed by the magnetic resonance electric impedance tomography method from current density distribution data. In this study, the use of current density imaging with MRI for monitoring current density distribution during delivery of irreversible electroporation pulses was demonstrated.

Methods

Using a modified single-shot RARE sequence, where four 3000 V and 100 μs long pulses were included at the start, current distribution between a pair of electrodes inserted in a liver tissue sample was imaged. Two repetitions of the sequence with phases of refocusing radiofrequency pulses 90° apart were needed to acquire one current density image. For each sample in total 45 current density images were acquired to follow a standard protocol for irreversible electroporation where 90 electric pulses are delivered at 1 Hz.

Results

Acquired current density images showed that the current density in the middle of the sample increased from first to last electric pulses by 60%, i.e. from 8 kA/m² to 13 kA/m² and that direction of the current path did not change with repeated electric pulses significantly.

Conclusions

The presented single-shot RARE-based current density imaging sequence was used successfully to image current distribution during delivery of short high-voltage electric pulses. The method has a potential to enable monitoring of tumor coverage by electric field during irreversible electroporation tissue ablation.

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.

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.

MRI of the human brain at 130 microtesla

We present in vivo images of the human brain acquired with an ultralow field MRI (ULFMRI) system operating at a magnetic field B0 ~ 130 μT. The system features prepolarization of the proton spins at Bp ~ 80 mT and detection of the NMR signals with a superconducting, second-derivative gradiometer inductively coupled to a superconducting quantum interference device (SQUID). We report measurements of the longitudinal relaxation time T1 of brain tissue, blood, and scalp fat at B0 and Bp, and cerebrospinal fluid at B0. We use these T1 values to construct inversion recovery sequences that we combine with Carr-Purcell-Meiboom-Gill echo trains to obtain images in which one species can be nulled and another species emphasized. In particular, we show an image in which only blood is visible. Such techniques greatly enhance the already high intrinsic T1 contrast obtainable at ULF. We further present 2D images of T1 and the transverse relaxation time T2 of the brain and show that, as expected at ULF, they exhibit similar contrast. Applications of brain ULFMRI include integration with systems for magnetoencephalography. More generally, these techniques may be applicable, for example, to the imaging of tumors without the need for a contrast agent and to modalities recently demonstrated with T1ρ contrast imaging (T1 in the rotating frame) at fields of 1.5 T and above.

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.

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.

Hybrid ultra-low-field MRI and magnetoencephalography system based on a commercial whole-head neuromagnetometer

Ultra-low-field MRI uses microtesla fields for signal encoding and sensitive superconducting quantum interference devices for signal detection. Similarly, modern magnetoencephalography (MEG) systems use arrays comprising hundreds of superconducting quantum interference device channels to measure the magnetic field generated by neuronal activity. In this article, hybrid MEG-MRI instrumentation based on a commercial whole-head MEG device is described. The combination of ultra-low-field MRI and MEG in a single device is expected to significantly reduce coregistration errors between the two modalities, to simplify MEG analysis, and to improve MEG localization accuracy. The sensor solutions, MRI coils (including a superconducting polarizing coil), an optimized pulse sequence, and a reconstruction method suitable for hybrid MEG-MRI measurements are described. The performance of the device is demonstrated by presenting ultra-low-field-MR images and MEG recordings that are compared with data obtained with a 3T scanner and a commercial MEG device.

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

A number of different methods have been developed in order to detect the spreading of neuronal currents by means of noninvasive imaging techniques. However, all of these are subjected to limitations in the temporal or spatial resolution. A new approach of neuronal current detection is based on the use of low-field nuclear magnetic resonance (LF-NMR) that records brain activity directly. In the following, we describe a phantom study in order to assess the feasibility of neuronal current detection using LF-NMR. In addition to that, necessary preliminary subject studies examining somatosensory evoked neuronal currents are presented. During the phantom study, the influences of two different neuronal time signals on (1)H-NMR signals were observed. The measurements were carried out by using a head phantom with an integrated current dipole to simulate neuronal activity. Two LF-NMR methods based on a DC and an AC (resonant) mechanism were utilized to study the feasibility of detecting both types of magnetic brain signals. Measurements were made inside an extremely magnetically shielded room by using a superconducting quantum interference device magnetometer system. The measurement principles were validated applying currents of higher intensity than those typical of the neuronal currents. Through stepwise reduction of the amplitude of the current dipole strength, the resolution limits of the two measuring procedures were found. The results indicate that it is necessary to improve the signal-to-noise ratio of the measurement system by at least a factor of 38 in order to detect typical human neuronal activity directly by means of LF-NMR. In addition to that, ways of achieving this factor are discussed.

Enhanced sensitivity current density imaging

One of the major weaknesses of current density imaging (CDI) is its poor sensitivity and therefore a need for the use of high voltage in CDI. In this work, a new CDI technique with enhanced sensitivity (ES-CDI) is presented. The ES-CDI sequence overcomes the sensitivity problem in samples with a long T(2) relaxation time that allows the use of a long current encoding period. As successful CDI detection is conditioned by a sufficiently large product of current and its application time a longer current encoding period enables the use of lower current and also lower voltage therefore significantly reducing any sample damage. In addition, the ES-CDI sequence also uses fast image signal acquisition and so enables heavy signal averaging and with it associated additional CDI sensitivity increase within the experiment time of the conventional CDI experiment. The feasibility of the ES-CDI sequence was tested on a model sample filled with physiological solution. Voltage of just 1 V and current application time of 800 ms were sufficient to detect current density of 20A/m(2) with a detection limit of 0.7A/m(2).

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.

NMR at very low fields

Although nuclear magnetic resonance in low fields around or below the Earth's magnetic field is almost as old as nuclear magnetic resonance itself, the recent years have experienced a revival of this technique that is opposed to the common trend towards higher and higher fields. The background of this development is the expectation that the low-field domain may open a new window for the study of molecular structure and dynamics. Here, we will give an overview on the specific features in the low-field domain, both from the technical and from the physical point of view. In addition, we present a short passage on the option of magnetic resonance imaging in fields of the micro-Tesla range.