Combination of BOLD-fMRI and VEP recordings for spin-echo MRI detection of primary magnetic effects caused by neuronal currents

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.

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.

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.

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.

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.

Detecting neuronal currents with MRI: a human study

PURPOSE

Recently developed neuronal current magnetic resonance imaging aims to directly detect neuronal currents associated with brain activity, but controversial results have been reported in different studies on human subjects. Although there is no dispute that local neuronal currents produce weak transient magnetic fields that would attenuate local MR signal intensity, there is not yet consensus as to whether this attenuation is detectable with present magnetic resonance imaging techniques. This study investigates the magnitude of neuronal current-induced signal attenuation in human visual cortex.

THEORY

A temporally well-controlled visual stimulation paradigm with a known neuronal firing pattern in monkey visual cortex provides a means of detecting and testing the magnitude of the neuronal current-induced attenuation in neuronal current magnetic resonance imaging.

METHODS

Placing a series of acquisition windows to fully cover the entire response duration enables a thorough detection of any detectable MR signal attenuation induced by the stimulus-evoked neuronal currents.

RESULTS

No significant neuronal current-induced MR signal attenuation was observed in the putative V1 in any participated subjects.

CONCLUSION

The present magnetic resonance imaging technique is not sensitive enough to detect neuronal current-induced MR signal attenuation, and the upper limit of this attenuation was found to be less than 0.07% under the study condition.

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.

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.

Octopus visual system: a functional MRI model for detecting neuronal electric currents without a blood-oxygen-level-dependent confound

PURPOSE

Despite the efforts that have been devoted to detecting the transient magnetic fields generated by neuronal firing, the conclusion that a functionally relevant signal can be measured with MRI is still controversial. For human studies of neuronal current MRI (nc-MRI), the blood-oxygen-level-dependent (BOLD) effect remains an irresolvable confound. For tissue studies where hemoglobin is removed, natural sensory stimulation is not possible. This study investigates the feasibility of detecting a physiologically induced nc-MRI signal in vivo in a BOLD-free environment.

METHODS

The cephalopod mollusc Octopus bimaculoides has vertebrate-like eyes, large optic lobes (OLs), and blood that does not contain hemoglobin. Visually evoked potentials were measured in the octopus retina and OL by electroretinogram and local field potential. nc-MRI scans were conducted at 9.4 Tesla to capture these activities.

RESULTS

Electrophysiological recording detected strong responses in the retina and OL in vivo; however, nc-MRI failed to demonstrate any statistically significant signal change with a detection threshold of 0.2° for phase and 0.2% for magnitude. Experiments in a dissected eye-OL preparation yielded similar results.

CONCLUSION

These findings in a large hemoglobin-free nervous system suggest that sensory evoked neuronal magnetic fields are too weak for direct detection with current MRI technology.

Direct detection of a single evoked action potential with MRS in Lumbricus terrestris

Functional MRI (fMRI) measures neural activity indirectly by detecting the signal change associated with the hemodynamic response following brain activation. In order to alleviate the temporal and spatial specificity problems associated with fMRI, a number of attempts have been made to detect neural magnetic fields (NMFs) with MRI directly, but have thus far provided conflicting results. In this study, we used MR to detect axonal NMFs in the median giant fiber of the earthworm, Lumbricus terrestris, by examining the free induction decay (FID) with a sampling interval of 0.32 ms. The earthworm nerve cords were isolated from the vasculature and stimulated at the threshold of action potential generation. FIDs were acquired shortly after the stimulation, and simultaneous field potential recordings identified the presence or absence of single evoked action potentials. FIDs acquired when the stimulus did not evoke an action potential were summed as background. The phase of the background-subtracted FID exhibited a systematic change, with a peak phase difference of (-1.2 ± 0.3) × 10(-5) radians occurring at a time corresponding to the timing of the action potential. In addition, we calculated the possible changes in the FID magnitude and phase caused by a simulated action potential using a volume conductor model. The measured phase difference matched the theoretical prediction well in both amplitude and temporal characteristics. This study provides the first evidence for the direct detection of a magnetic field from an evoked action potential using MR.

Phase stability in fMRI time series: effect of noise regression, off-resonance correction and spatial filtering techniques

Although the majority of fMRI studies exploit magnitude changes only, there is an increasing interest regarding the potential additive information conveyed by the phase signal. This integrated part of the complex number furnished by the MR scanners can also be used for exploring direct detection of neuronal activity and for thermography. Few studies have explicitly addressed the issue of the available signal stability in the context of phase time-series, and therefore we explored the spatial pattern of frequency specific phase fluctuations, and evaluated the effect of physiological noise components (heart beat and respiration) on the phase signal. Three categories of retrospective noise reduction techniques were explored and the temporal signal stability was evaluated in terms of a physiologic noise model, for seven fMRI measurement protocols in eight healthy subjects at 3T, for segmented CSF, gray and white matter voxels. We confirmed that for most processing methods, an efficient use of the phase information is hampered by the fact that noise from physiological and instrumental sources contributes significantly more to the phase than to the magnitude instability. Noise regression based on the phase evolution of the central k-space point, RETROICOR, or an orthonormalized combination of these were able to reduce their impact, but without bringing phase stability down to levels expected from the magnitude signal. Similar results were obtained after targeted removal of scan-to-scan variations in the bulk magnetic field by the dynamic off-resonance in k-space (DORK) method and by the temporal off-resonance alignment of single-echo time series technique (TOAST). We found that spatial high-pass filtering was necessary, and in vivo a Gaussian filter width of 20mm was sufficient to suppress physiological noise and bring the phase fluctuations to magnitude levels. Stronger filters brought the fluctuations down to levels dictated by thermal noise contributions, and for 62.5mm(3) voxels the phase stability was as low as 5 mrad (0.27°). In conditions of low SNR(o) and high temporal sampling rate (short TR); we achieved an upper bound for the phase instabilities at 0.0017 ppm, which is close to the dHb contribution to the GM/WM phase contrast.

On multiple alternating steady states induced by periodic spin phase perturbation waveforms

Direct measurement of neural currents by means of MRI can potentially open a high temporal resolution (10-100 ms) window applicable for monitoring dynamics of neuronal activity without loss of the high spatial resolution afforded by MRI. Previously, we have shown that the alternating balanced steady state imaging affords high sensitivity to weak periodic currents owing to its amplification of periodic spin phase perturbations. This technique, however, requires precise synchronization of such perturbations to the radiofrequency pulses. Herein, we extend alternating balanced steady state imaging to multiple balanced alternating steady states for estimation of neural current waveforms. Simulations and phantom experiments show that the off-resonance profile of the multiple alternating steady state signal carries information about the frequency content of driving waveforms. In addition, the method is less sensitive than alternating balanced steady state to precise waveform timing relative to radiofrequency pulses. Thus, multiple alternating steady state technique is potentially applicable to MR imaging of the waveforms of periodic neuronal activity.

Modeling MR signal change induced by oxygen effect in neural tissue preparations of various geometries

Tissue preparation has recently been utilized for detection of neuronal activation in multiple non-BOLD based functional MRI studies to eliminate vascular contamination. However, undesired signal change could still occur in such systems due to the concentration change of dissolved O(2) upon tissue activation. To estimate the impact of such effects, the O(2) concentration distribution and the consequent susceptibility field in tissue-solution systems were simulated with various tissue geometries and experimental parameters. Our results indicate that substantial signal change between the resting and activated states could potentially be induced by the O(2) effect in highly oxygenated solutions, and thus caution should be taken in interpreting any signal change observed in such experiments.

Detection of neuronal current MRI in human without BOLD contamination

Controversial results regarding the detectability of neuronal current magnetic resonance imaging (ncMRI) have been reported in different studies on human subjects. In all the previous studies, the ncMRI signal was detected under a continuous and paradigm task-induced blood oxygen level dependent (BOLD) signal background. The aim of this study is to investigate the possibility of detecting ncMRI signal in human brain in the situation that task-induced BOLD background is absent or minimum. In this study, by adopting an event-related visuomotor paradigm with long interstimulus interval (=20 s), the ncMRI signal was detected when the BOLD signal fully returned to its baseline, and the potential BOLD background contamination was avoided effectively. The results showed that the visuomotor stimulation elicited BOLD activation in visual and motor cortices, but no significant ncMRI signal change (in magnitude) was detected in human brain. These experimental findings are consistent with theoretical predications.

Modeling neuronal current MRI signal with human neuron

Up to date, no consensus has been achieved regarding the possibility of detecting neuronal currents by MRI (ncMRI) in human brain. To evaluate the detectability of ncMRI, an effective way is to simulate ncMRI signal with the realistic neuronal geometry and electrophysiological processes. Unfortunately, previous realistic ncMRI models are based on rat and monkey neurons. The species difference in neuronal morphology and physiology would prevent these models from simulating the ncMRI signal accurately in human subjects. The aim of this study is to bridge this gap by establishing a realistic ncMRI model specifically for human cerebral cortex. In this model, the ncMRI signal was simulated using anatomically reconstructed human pyramidal neurons and their biophysical properties. The modeling results showed that the amplitude of ncMRI signal significantly depends on the density of synchronously firing neurons and imaging conditions such as position of imaging voxel, direction of main magnetic field (B(0) ) relative to the cortical surface and echo time. The results indicated that physiologically-evoked ncMRI signal is too weak to be detected (magnitude/phase change ≤ -1.4 × 10(-6) /0.02°), but the phase signal induced by spontaneous activity may reach a detectable level (up to 0.2°) in favorable conditions.

Direct MRI detection of the neuronal magnetic field: the effect of the dendrite branch

In recent years, neuronal current MRI (nc-MRI) was proposed as a new imaging method to directly map the magnetic field change caused by neuronal activity. Nc-MRI could offer improved spatial and temporal resolution compared to blood hemodynamics-based functional magnetic resonance imaging (fMRI). In this paper, with a finite current dipole as the model of dendrite or dendrite branch, we investigated the spatial distribution of the magnetic field generated by synchronously activated neurons to evaluate the possibility of nc-MRI. Our simulations imply that the existence of a dendrite branch may not only increase the strength of the neuronal magnetic field (NMF), but also raise the non-uniform and unsymmetry of the NMF; therefore, it can enhance the detectability of the neuronal current magnetic field by MRI directly. The results show that the signal phase shift is enlarged, but it is unstable and is still very small, <<1 radian, while the magnitude signal may be strong enough for a typical MRI voxel to be detected. We suggest making further efforts to measure the magnitude signal which may induce a large effect in an nc-MRI experiment.

Direct MRI mapping of neuronal activity evoked by electrical stimulation of the median nerve at the right wrist

Magnetic source MRI (msMRI) has being developed recently for direct detections of neuronal magnetic fields to map brain activity. However, controversial results have been reported by different research groups. In this study, more evidence was provided to demonstrate that the neuronal current signal could be detected by MRI using a rapid median nerve stimulation paradigm. The experiments were performed on six normal human participants to investigate the temporal specificity of the effect, as well as inter- and intrasubject reproducibility. Significant activation of contralateral primary sensory cortex (S1) was detected 80 ms after stimulation onset (corresponding to the P80 evoked potential peak). The 80-ms latency S1 activation was observed over three independent sessions for one subject and for all six participants. The magnitude of the signal change was 0.2-0.3%. Coinciding with our expectations, no S1 activation was found when MRI data acquisitions were targeted at the N20 and P30 peaks because of mutual cancellation of magnetic fields generated by those peaks. The results demonstrated good reproducibility of S1 activations and indicated that the S1 activations most likely originated from neuronal magnetic field rather than hemodynamic response.

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.

Dynamic MRI of small electrical activity

Neuroscience methods entailing in vivo measurements of brain activity have greatly contributed to our understanding of brain function for the past decades, from the invasive early studies in animals using single-cell electrical recordings, to the noninvasive techniques in humans of scalp-recorded electroencephalography (EEG) and magnetoencephalography (MEG), positron emission tomography (PET), and, most recently, blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI). A central objective of these techniques is to measure neuronal activities with high spatial and temporal resolution. Each of these methods, however, has substantial limitations in this regard. Single-cell recording is invasive and only typically records cellular activity in a single location; EEG/MEG cannot generally provide accurate and unambiguous delineations of neuronal activation spatially; and the most sophisticated BOLD-based fMRI methods are still fundamentally limited by their dependence on the very slow hemodynamic responses upon which they are based. Even the latest neuroimaging methodology (e.g., multimodal EEG/fMRI) does not yet unambiguously provide accurate localization of neuronal activation spatially and temporally. There is hence a need to further develop noninvasive imaging methods that can directly image neuroelectric activity and thus truly achieve a high temporal resolution and spatial specificity in humans. Here, we discuss the theory, implementation, and potential utility of an MRI technique termed Lorentz effect imaging (LEI) that can detect spatially incoherent yet temporally synchronized, minute electrical activities in the neural amplitude range (microamperes) when they occur in a strong magnetic field. Moreover, we demonstrate with our preliminary results in phantoms and in vivo, the feasibility of imaging such activities with a temporal resolution on the order of milliseconds.

The effect of physiological noise in phase functional magnetic resonance imaging: from blood oxygen level-dependent effects to direct detection of neuronal currents

Recently, the possibility to use both magnitude and phase image sets for the statistical evaluation of fMRI has been proposed, with the prospective of increasing both statistical power and the spatial specificity. In the present work, several issues that affect the spatial and temporal stability in fMRI phase time series in the presence of physiologic noise processes are reviewed, discussed and illustrated by experiments performed at 3 T. The observed phase value is a fingerprint of the underlying voxel averaged magnetic field variations. Those related to physiological processes can be considered static or dynamic in relation to the temporal scale of a 2D acquisition and will play out on different spatial scales as well: globally across the entire images slice, and locally depending on the constituents and their relative fractions inside the MRI voxel. The 'static' respiration-induced effects lead to magneto-mechanic scan-to-scan variations in the global magnetic field but may also contribute to local BOLD fluctuations due to respiration-related variations in arterial carbon dioxide. Likewise, the 'dynamic' cardiac-related effects will lead to global susceptibility effects caused by pulsatile motion of the brain as well as local blood pressure-related changes in BOLD and changes in blood flow velocity. Finally, subject motion may lead to variations in both local and global tissue susceptibility that will be especially pronounced close to air cavities. Since dissimilar manifestations of physiological processes can be expected in phase and in magnitude images, a direct relationship between phase and magnitude scan-to-scan fluctuations cannot be assumed a priori. Therefore three different models were defined for the phase stability, each dependent on the relation between phase and magnitude variations and the best will depend on the underlying noise processes. By experiments on healthy volunteers at rest, we showed that phase stability depends on the type of post-processing and can be improved by reducing the low-frequency respiration-induced mechano-magnetic effects. Although the manifestations of physiological noise were in general more pronounced in phase than in magnitude images, due to phase wraps and global Bo effects, we suggest that a phase stability similar to that found in magnitude could theoretically be achieved by adequate correction methods. Moreover, as suggested by our experimental data regarding BOLD-related phase effects, phase stability could even supersede magnitude stability in voxels covering dense microvascular networks with BOLD-related fluctuations as the dominant noise contributor. In the interest of the quality of both BOLD-based and nc-MRI methods, future studies are required to find alternative methods that can improve phase stability, designed to match the temporal and spatial scale of the underlying neuronal activity.

Microscopic investigation of the resonant mechanism for the implementation of nc-MRI at ultra-low field MRI

In this paper the possible use of the resonant mechanism between some spectral components of the neuronal activity and the spin dynamics in ultra-low field MRI experiments--for the implementation of the nc-MRI techniques and proposed by Kraus et al., 2008--is investigated by means of "realistic" simulations of the neuronal activity of a modelled neuronal network. Previously characterized digital neurons are used to reproduce neuronal currents based on biophysical details and the distribution of the local magnetic field inside a MRI cubic voxel (having a dimension of 1.2 mm) is evaluated. The properties of the water proton spin dynamics as a consequence of the neuronal field and of external applied fields are extrapolated integrating the Bloch equations. The characteristics of the expected MR signals are discussed in relation to the specifics of the NMR sequence used and to the properties of the neuronal activity. The great potentialities of the technique are provided by: a) the possible easy implementation of the technique, b) the possible cheap instrumentation required; c) the flexibility of the ultra-low field systems.

Imaging periodic currents using alternating balanced steady-state free precession

Existing functional brain MR imaging methods detect neuronal activity only indirectly via a surrogate signal such as deoxyhemoglobin concentration in the vascular bed of cerebral parenchyma. It has been recently proposed that neuronal currents may be measurable directly using MRI (ncMRI). However, limited success has been reported in neuronal current detection studies that used standard gradient or spin echo pulse sequences. The balanced steady-state free precession (bSSFP) pulse sequence is unique in that it can afford the highest known SNR efficiency and is exquisitely sensitive to perturbations in free precession phase. It is reported herein that when a spin phase-perturbing periodic current is locked to an RF pulse train, phase perturbations are accumulated across multiple RF excitations and the spin magnetization reaches an alternating balanced steady state (ABSS) that effectively amplifies the phase perturbations due to the current. The alternation of the ABSS signal therefore is highly sensitive to weak periodic currents. Current phantom experiments employing ABSS imaging resulted in detection of magnetic field variations as small as 0.15nT in scans lasting for 36 sec, which is more sensitive than using gradient-recalled echo imaging.

Lorentz effect imaging of ionic currents in solution

Current functional MRI techniques relying on hemodynamic modulations are inherently limited in their ability to accurately localize neural activity in space and time. To address these limitations, we previously proposed a novel technique based on the Lorentz effect and demonstrated its ability to directly image minute electrical activity with a millisecond temporal resolution in gel phantoms containing conductive wires as well as in the human median nerve in vivo. To better characterize its contrast mechanism and ultimately further improve its sensitivity for in vivo applications, we now apply this technique to image ionic currents in solution, which serve as a better model for neural conduction in biological systems than the electronic currents in conductive wires used in previous phantom studies. Our results demonstrate that ionic currents with durations and current densities on the same order of magnitude as those induced by neuroelectric activity in nerve fibers and in the brain can be detected.

Failure to direct detect magnetic field dephasing corresponding to ERP generation

Functional magnetic resonance imaging (fMRI) has become the method of choice for mapping brain activity in human subjects and detects changes in regional blood oxygenation and volume associated with local changes in neuronal activity. While imaging based on blood oxygenation level dependent (BOLD) contrast has good spatial resolution and sensitivity, the hemodynamic signal develops relatively slowly and is only indirectly related to neuronal activity. An alternative approach termed magnetic source magnetic resonance imaging (msMRI) is based on the premise that neural activity may be mapped by magnetic resonance imaging (MRI) with greater temporal resolution by detecting the local magnetic field perturbations associated with local neuronal electric currents. We used a hybrid ms/BOLD MRI method to investigate whether msMRI could detect signal changes that occur simultaneously at the time of the production of well-defined event-related potentials, the P300 and N170, in regions that previously have been identified as generators of these electrical signals. Robust BOLD activations occurred after some seconds, but we were unable to detect any significant changes in the T2*-weighted signal in these locations that correlated temporally with the timings of the evoked response potentials (ERPs).

Modeling oxygen effects in tissue-preparation neuronal-current MRI

Tissue-preparation neuronal-current MRI (ncMRI) was recently developed to directly detect neuronal activity without hemodynamic contamination. However, as a paramagnetic substance, the oxygen molecules present in the tissue may also alter the ncMRI signal through relaxivity and susceptibility effects. To study the effects of oxygen on the ncMRI signal and estimate their impact on tissue-preparation experiments, oxygen-induced MRI signal changes were formulated as a function of oxygen concentration (OC) of gas, oxygen consumption rate, and imaging parameters. Under favorable conditions of these parameters, the maximum oxygen-induced signal magnitude and phase change were estimated to be 0.32% and 3.85 degrees , respectively. Considering that the ncMRI signal changes obtained in previous tissue-preparation experiments (3-5% in magnitude, 0.8-1.7 degrees in phase) were tens or hundreds of times larger than the corresponding oxygen-induced signal changes (0.03% in magnitude, 0.03-0.07 degrees in phase), it is concluded that the oxygen had negligible effects in the previous experiments.

Realistic simulations of neuronal activity: a contribution to the debate on direct detection of neuronal currents by MRI

Many efforts have been done in order to preview the properties of the magnetic resonance (MR) signals produced by the neuronal currents using simulations. In this paper, starting with a detailed calculation of the magnetic field produced by the neuronal currents propagating over single hippocampal CA1 pyramidal neurons placed inside a cubic MR voxel of length 1.2 mm, we proceeded on the estimation of the phase and magnitude MR signals. We then extended the results to layers of parallel and synchronous similar neurons and to ensembles of layers, considering different echo times, voxel volumes and neuronal densities. The descriptions of the neurons and of their electrical activity took into account the real neuronal morphologies and the physiology of the neuronal events. Our results concern: (a) the expected time course of the MR signals produced by the neuronal currents in the brain, based on physiological and anatomical properties; (b) the different contributions of post-synaptic potentials and of action potentials to the MR signals; (c) the estimation of the equivalent current dipole and the influence of its orientation with respect to the external magnetic field on the observable MR signal variations; (d) the size of the estimated neuronal current induced phase and magnitude MR signal changes with respect to the echo time, voxel-size and neuronal density. The inclusion of realistic neuronal properties into the simulation introduces new information that can be helpful for the design of MR sequences for the direct detection of neuronal current effects and the testing of bio-electromagnetic models.

Heart beats brain: The problem of detecting alpha waves by neuronal current imaging in joint EEG–MRI experiments

It has been suggested recently that the influence of the neuro-magnetic field should make electrical brain activity directly detectable by MRI. To test this hypothesis, we performed combined EEG-MRI experiments which aim to localize the neuronal current sources of alpha waves (8-12 Hz), one of the most prominent EEG phenomena in humans. A detailed analysis of cross-spectral coherence between simultaneously recorded EEG and MRI time series revealed no sign of alpha waves. Instead the EEG-MRI approach was found to be hampered by artefacts due to cardiac pulsation, which extend into the frequency band of alpha waves. Separate brain displacement mapping experiments confirmed that not only the EEG but also the MRI signal is confounded by harmonics of the cardiac frequency even at 10 Hz and beyond. This well-known ballistocardiogram artefact cannot be avoided or eliminated entirely by available signal processing techniques. Therefore we must conclude that current EEG-MRI methodology based on correlation analysis lacks not only the sensitivity but also the specificity required for the reliable detection of alpha waves.

Modelling the magnetic signature of neuronal tissue

Neuronal communication in the brain involves electrochemical currents, which produce magnetic fields. Stimulus-evoked brain responses lead to changes in these fields and can be studied using magneto- and electro-encephalography (MEG/EEG). In this paper we model the spatiotemporal distribution of the magnetic field of a physiologically idealized but anatomically realistic neuron to assess the possibility of using magnetic resonance imaging (MRI) for directly mapping the neuronal currents in the human brain. Our results show that the magnetic field several centimeters from the centre of the neuron is well approximated by a dipole source, but the field close to the neuron is not, a finding particularly important for understanding the possible contrast mechanism underlying the use of MRI to detect and locate these currents. We discuss the importance of the spatiotemporal characteristics of the magnetic field in cortical tissue for evaluating and optimizing an experiment based on this mechanism and establish an upper bound for the expected MRI signal change due to stimulus-induced cortical response. Our simulations show that the expected change of the signal magnitude is 1.6% and its phase shift is 1 degrees . An unexpected finding of this work is that the cortical orientation with respect to the external magnetic field has little effect on the predicted MRI contrast. This encouraging result shows that magnetic resonance contrast directly based on the neuronal currents present in the cortex is theoretically a feasible imaging technique. MRI contrast generation based on neuronal currents depends on the dendritic architecture and we obtained high-resolution optical images of cortical tissue to discuss the spatial structure of the magnetic field in grey matter.

Inability to directly detect magnetic field changes associated with neuronal activity

The ability to directly detect neuronal magnetic fields by MRI would help investigators achieve the "holy grail" of neuroimaging, namely both high spatial and temporal resolution. Both positive and negative findings have been reported in the literature, with no clear consensus as to the feasibility of direct detection. The aim of this study was to replicate one of the most promising published in vivo results. A second aim was to investigate the use of steady-state visual evoked potentials (ssVEPs), which give a large evoked response and offer a well-controlled approach because the frequency of the neuronal response can be dictated by the experimenter. For both studies we used a general linear model (GLM) that included regressors for both the expected blood oxygen level-dependent (BOLD) signal and the magnetic source (MS) signal. The results showed no activity that could be attributed to the neuromagnetic signals in either study, and no frequency component corresponding to the frequency of the ssVEPs. This study demonstrates that for the particular stimuli and hardware used, the sensitivity of the magnitude MRI signal to detect evoked neuronal currents is too low to be of practical use.

Finding neuroelectric activity under magnetic-field oscillations (NAMO) with magnetic resonance imaging in vivo

Neuroimaging techniques are among the most important tools for investigating the function of the human nervous system and for improving the clinical diagnosis of neurological disorders. However, most commonly used techniques are limited by their invasiveness or their inability to accurately localize neural activity in space or time. Previous attempts at using MRI to directly image neuroelectric activity in vivo through the detection of magnetic field changes induced by neuronal currents have been challenging because of the extremely small signal changes and confounding factors such as hemodynamic modulations. Here we describe an MRI technique that uses oscillating magnetic field gradients to significantly amplify and detect the Lorentz effect induced by neuroelectric activity, and we demonstrate its effectiveness in imaging sensory nerve activation in vivo in the human median nerve during electrical stimulation of the wrist. This direct, real-time, and noninvasive neuroimaging technique may potentially find broad applications in neurosciences.

Challenges for detection of neuronal currents by MRI

Neuronal current MRI (nc-MRI) is an imaging method that directly maps magnetic field changes caused by neuronal currents with, at the same time, a high spatial and temporal resolution. A viable nc-MRI method would be of great benefit, both for the study of human brain function and for clinical applications in the field of epilepsy, especially for the noninvasive presurgical mapping of epileptogenic foci. A survey of fundamental issues in nc-MRI is reviewed, and challenges for future developments of the method are described within this context. Particularly, an overview of the models for signal generation is given, and the origin and physiology of different sources of neuronal currents are described. Prospects for predicting neuronal currents by electromagnetic field mapping and using this information, both a priori and a posteriori, for nc-MRI are considered. Ways of increasing specificity in nc-MRI by minimizing secondary hemodynamic and metabolic effects are described as well as means of optimizing the nc-MRI method for pushing the detection limit. Previously published works are described within these categories and future directions for nc-MRI are proposed.

Synchronized detection of minute electrical currents with MRI using Lorentz effect imaging

The blood oxygenation level-dependent (BOLD) effect is the most commonly used contrast mechanism in functional magnetic resonance imaging (fMRI), due to its relatively high spatial resolution and sensitivity. However, the ability of BOLD fMRI to accurately localize neuronal activation in space and time is limited by the inherent hemodynamic modulation. There is hence a need to develop alternative MRI methods that can directly image neuroelectric activity, thereby achieving both a high temporal resolution and spatial specificity as compared to conventional BOLD fMRI. In this paper, we extend the Lorentz effect imaging technique, which can detect spatially incoherent yet temporally synchronized minute electrical activity in a strong magnetic field, and demonstrate its feasibility for imaging randomly oriented electrical currents on the order of microamperes with a temporal resolution on the order of milliseconds in gel phantoms. This constitutes a promising step towards its application to direct imaging of neuroelectric activity in vivo, which has the same order of current density and temporal synchrony.

Frontiers of brain mapping using MRI

Over the past dozen years, the use of MRI techniques to map brain function (fMRI) has sparked a great deal of research. The ability of fMRI to image several different physiological processes concurrently (i.e., blood oxygenation, blood flow, metabolism) and noninvasively over large volumes make it the ideal choice for many different areas of neuroscience research in addition to countless applications in clinical settings. Furthermore, with the advent of high magnetic fields (and other hardware advancements, i.e., parallel imaging) for both human and animal research, spatial and temporal resolutions continue to be pushed to higher levels because of increases in the sensitivity as well as specificity of MR-detectable functional signals. fMRI methodology continues to grow and has the ability to cater to many different research applications. There seems to be no foreseeable end in sight to the advancement of fMRI techniques and its subsequent use in basic research as well as in clinical settings. In this work, fMRI techniques and the ongoing development of existing techniques are discussed with implications for the future of fMRI.