A dynamic causal model of the coupling between pulse stimulation and neural activity

We present a dynamic causal model that can explain context-dependent changes in neural responses, in the rat barrel cortex, to an electrical whisker stimulation at different frequencies. Neural responses were measured in terms of local field potentials. These were converted into current source density (CSD) data, and the time series of the CSD sink was extracted to provide a time series response train. The model structure consists of three layers (approximating the responses from the brain stem to the thalamus and then the barrel cortex), and the latter two layers contain nonlinearly coupled modules of linear second-order dynamic systems. The interaction of these modules forms a nonlinear regulatory system that determines the temporal structure of the neural response amplitude for the thalamic and cortical layers. The model is based on the measured population dynamics of neurons rather than the dynamics of a single neuron and was evaluated against CSD data from experiments with varying stimulation frequency (1–40 Hz), random pulse trains, and awake and anesthetized animals. The model parameters obtained by optimization for different physiological conditions (anesthetized or awake) were significantly different. Following Friston, Mechelli, Turner, and Price (2000), this work is part of a formal mathematical system currently being developed (Zheng et al., 2005) that links stimulation to the blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal through neural activity and hemodynamic variables. The importance of the model described here is that it can be used to invert the hemodynamic measurements of changes in blood flow to estimate the underlying neural activity.

Functional Magnetic Resonance Imaging Data Analysis by Autoregressive Estimator

The autoregressive (AR) estimator, a non-parametric method, is used to analyze functional magnetic resonance imaging (fMRI) data. The same method has been used, with success, in several other time series data analysis. It uses exclusively the available experimental data points to estimate the most plausible power spectra compatible with the experimental data and there is no need to make any assumption about non-measured points. The time series, obtained from fMRI block paradigm data, is analyzed by the AR method to determine the brain active regions involved in the processing of a given stimulus. This method is considerably more reliable than the fast Fourier transform or the parametric methods. The time series corresponding to each image pixel is analyzed using the AR estimator and the corresponding poles are obtained. The pole distribution gives the shape of power spectra, and the pixels with poles at the stimulation frequency are considered as the active regions. The method was applied in simulated and real data, its superiority is shown by the receiver operating characteristic curves which were obtained using the simulated data.

Muscle metabolites: functional MR spectroscopy during exercise imposed by tetanic electrical nerve stimulation

Permission from the ethics committee and informed consent were obtained. The purpose of this study was to prospectively evaluate a method developed for the noninvasive assessment of muscle metabolites during exercise. Hydrogen 1 magnetic resonance (MR) spectroscopy peaks were measured during tetanic isometric muscle contraction imposed by supramaximal repetitive nerve stimulation. The kinetics of creatine-phosphocreatine and acetylcarnitine signal changes (P < .001) could be assessed continuously before, during, and after exercise. The control peak (trimethylammonium compounds), which served as an internal reference, did not change. This technique-that is, functional MR spectroscopy-opens the possibility for noninvasive diagnostic muscle metabolite testing in a clinical setting.

A three-dimensional histological atlas of the human basal ganglia. II. Atlas deformation strategy and evaluation in deep brain stimulation for Parkinson disease

OBJECT: The localization of any given target in the brain has become a challenging issue because of the increased use of deep brain stimulation to treat Parkinson disease, dystonia, and nonmotor diseases (for example, Tourette syndrome, obsessive compulsive disorders, and depression). The aim of this study was to develop an automated method of adapting an atlas of the human basal ganglia to the brains of individual patients. METHODS: Magnetic resonance images of the brain specimen were obtained before extraction from the skull and histological processing. Adaptation of the atlas to individual patient anatomy was performed by reshaping the atlas MR images to the images obtained in the individual patient using a hierarchical registration applied to a region of interest centered on the basal ganglia, and then applying the reshaping matrix to the atlas surfaces. RESULTS: Results were evaluated by direct visual inspection of the structures visible on MR images and atlas anatomy, by comparison with electrophysiological intraoperative data, and with previous atlas studies in patients with Parkinson disease. The method was both robust and accurate, never failing to provide an anatomically reliable atlas to patient registration. The registration obtained did not exceed a 1-mm mismatch with the electrophysiological signatures in the region of the subthalamic nucleus. CONCLUSIONS: This registration method applied to the basal ganglia atlas forms a powerful and reliable method for determining deep brain stimulation targets within the basal ganglia of individual patients.

Signal changes in the spinal cord of the rat after injection of formalin into the hindpaw: Characterization using functional magnetic resonance imaging

Changes in metabolism and local circulation occur in the spinal cord during peripheral noxious stimulation. Evidence is presented that this stimulation also causes signal intensity alterations in functional magnetic resonance images of the spinal cord during formalin-induced pain. These results indicate the potential of functional magnetic resonance imaging in assessing noninvasively the extent and intensity of spinal cord excitation in this well characterized pain model. Therefore, the aim of this study was to establish functional magnetic resonance imaging as a noninvasive method to characterize temporal changes in the spinal cord after a single injection of 50 μl of formalin subcutaneously into the hindpaw of the anesthetized rat. This challenge produced a biphasic licking activity in the freely moving conscious animal. Images of the spinal cord were acquired within 2 min, enabling monitoring of the site and the temporal evolution of the signal changes during the development of formalin-induced hyperalgesia without the need of any surgical procedure. The time course of changes in the spinal cord functional image in the isoflurane-anesthetized animal was similar to that obtained from behavioral experiments. Also, comparable physiological data, control experiments, and the inhibition of a response through application of the local anesthetic agent lidocaine indicate that the signal changes observed after formalin injection were specifically related to excitability changes in the relevant segments of the lumbar spinal cord. This approach could be useful to characterize different models of pain and hyperalgesia and, more importantly, to evaluate effects of analgesic drugs.

Retinotopic mapping of lateral geniculate nucleus in humans using functional magnetic resonance imaging

Subcortical nuclei in the thalamus, which play an important role in many functions of the human brain, provide challenging targets for functional mapping with neuroimaging techniques because of their small sizes and deep locations. In this study, we explore the capability of high-resolution functional magnetic resonance imaging at 4 Tesla for mapping the retinotopic organization in the lateral geniculate nucleus (LGN). Our results show that the hemifield visual stimulation only activates LGN in the contralateral hemisphere, and the lower-field and upper-field visual stimulations activate the superior and inferior portion of LGN, respectively. These results reveal a similar retinotopic organization between the human and nonhuman primate LGN and between LGN and the primary visual cortex. We conclude that high-resolution functional magnetic resonance imaging is capable of functional mapping of suborganizations in small nuclei together with cortical activation. This will have an impact for studying the thalamocortical networks in the human brain.

Magnetic resonance imaging (MRI) detection of the murine brain response to light: temporal differentiation and negative functional MRI changes.

Using a 9.4 T MRI instrument, we have obtained images of the mouse brain response to photic stimulation during a period between deep anesthesia and the early stages of arousal. The large image enhancements we observe (often >30%) are consistent with literature results extrapolated to 9.4 T. However, there are also two unusual aspects to our findings. (i) The visual area of the brain responds only to changes in stimulus intensity, suggesting that we directly detect operations of the M visual system pathway. Such a channel has been observed in mice by invasive electrophysiology, and described in detail for primates. (ii) Along with the typical positive response in the area of the occipital portion of the brain containing the visual cortex, another area displays decreased signal intensity upon stimulation.

Activation of single whisker barrel in rat brain localized by functional magnetic resonance imaging.

The previously established cortical representation of rat whiskers in layer IV of the cortex contains distinct cylindrical columns of cellular aggregates, which are termed barrels and correlate in a one-to-one relation to whiskers on the contralateral rat face. In the present study, functional magnetic resonance imaging (fMRI) of the rat brain was used to map whisker barrel activation during mechanical up-down movement (+/- 2.5 mm amplitude at 8 Hz) of single/multiple whisker(s). Multislice gradient echo fMRI experiments were performed at 7 T with in-plane image resolution of 220 x 220 microns, slice thickness of 1 mm, and echo time of 16 ms. Highly significant (P < 0.001) and localized contralateral regions of activation were observed upon stimulation of single/multiple whisker(s). In all experiments (n = 10), the locations of activation relative to bregma and midline were highly correlated with the neuroanatomical position of the corresponding whisker barrels, and the results were reproducible intra- and interanimal. Our results indicate that fMRI based on blood oxygenation level-dependent image contrast has the sensitivity to depict activation of a single whisker barrel in the rat brain. This noninvasive technique will supplement existing methods in the study of rat barrel cortex and should be particularly useful for the long-term investigations of central nervous system in the same animal.

Three-dimensional functional magnetic resonance imaging of human brain on a clinical 1.5-T scanner.

Functional magnetic resonance imaging (fMRI) is a tool for mapping brain function that utilizes neuronal activity-induced changes in blood oxygenation. An efficient three-dimensional fMRI method is presented for imaging brain activity on conventional, widely available, 1.5-T scanners, without additional hardware. This approach uses large magnetic susceptibility weighting based on the echo-shifting principle combined with multiple gradient echoes per excitation. Motor stimulation, induced by self-paced finger tapping, reliably produced significant signal increase in the hand region of the contralateral primary motor cortex in every subject tested.

Illusory contours activate specific regions in human visual cortex: evidence from functional magnetic resonance imaging.

The neural basis for perceptual grouping operations in the human visual system, including the processes which generate illusory contours, is fundamental to understanding human vision. We have employed functional magnetic resonance imaging to investigate these processes noninvasively. Images were acquired on a GE Signa 1.5T scanner equipped for echo planar imaging with an in-plane resolution of 1.5 x 1.5 mm and slice thicknesses of 3.0 or 5.0 mm. Visual stimuli included nonaligned inducers (pacmen) that created no perceptual contours, similar inducers at the corners of a Kanizsa square that created illusory contours, and a real square formed by continuous contours. Multiple contiguous axial slices were acquired during baseline, visual stimulation, and poststimulation periods. Activated regions were identified by a multistage statistical analysis of the activation for each volume element sampled and were compared across conditions. Specific brain regions were activated in extrastriate cortex when the illusory contours were perceived but not during conditions when the illusory contours were absent. These unique regions were found primarily in the right hemisphere for all four subjects and demonstrate that specific brain regions are activated during the kind of perceptual grouping operations involved in illusory contour perception.

A distributed equivalent magnetic current based FDTD method for the calculation of E-fields induced by gradient coils

This paper evaluates a new, low-frequency finite-difference time-domain method applied to the problem of induced E-fields/eddy currents in the human body resulting from the pulsed magnetic field gradients in MRI. In this algorithm, a distributed equivalent magnetic current is proposed as the electromagnetic source and is obtained by quasistatic calculation of the empty coil's vector potential or measurements therein. This technique circumvents the discretization of complicated gradient coil geometries into a mesh of Yee cells, and thereby enables any type of gradient coil modelling or other complex low frequency sources. The proposed method has been verified against an example with an analytical solution. Results are presented showing the spatial distribution of gradient-induced electric fields in a multi-layered spherical phantom model and a complete body model. (C) 2004 Elsevier Inc. All rights reserved.

Chronotopographical analysis of the pattern onset visual evoked magnetic response (VEMR):Implications for waveform peak identification

The topography of the visual evoked magnetic response (VEMR) to a pattern onset stimulus was studied in five normal subjects using a single channel BTi magnetometer. Topographic distributions were analysed at regular intervals following stimulus onset (chronotopograpby). Two distinct field distributions were observed with half field stimulation: (1) activity corresponding to the C11 m which remains stable for an average of 34 msec and (2) activity corresponding to the C111 m which remains stable for about 50 msec. However, the full field topography of the largest peak within the first 130 msec does not have a predictable latency or topography in different subjects. The data suggest that the appearance of this peak is dependent on the amplitude, latency and duration of the half field C11 m peaks and the efficiency of half field summation. Hence, topographic mapping is essential to correctly identify the C11 m peak in a full field response as waveform morphology, peak latency and polarity are not reliable indicators. © 1993.

Topographic mapping and source localization of the pattern onset visual evoked magnetic response

The topography of the visual evoked magnetic response (VEMR) to a pattern onset stimulus was investigated using 4 check sizes and 3 contrast levels. The pattern onset response consists of three early components within the first 200ms, CIm, CIIm and CIIIm. The CIIm is usually of high amplitude and is very consistent in latency within a subject. Half field (HF) stimuli produce their strongest response over the contralateral hemisphere; the RHF stimulus exhibiting a lower positivity (outgoing field) and an upper negativity (ingoing field), rotated towards the midline. LHF stimulation produced the opposite response, a lower negative and an upper positive. Larger check sizes produce a single area of ingoing and outgoing field while smaller checks produce on area of ingoing and outgoing field over each hemisphere. Latency did not appear to vary with change in contrast but amplitudes increased with increasing contrast. A more detailed topographic study incorporating source localisation procedures suggested a source for CIIm - 4cm below the scalp, close to the midline with current flowing towards the lateral surface. Similar depth and position estimates but with opposite polarity were obtained for the pattern shift P100m previously. Hence, the P100m and the CIIm may originate in similar areas of visual cortex but reveal different aspects of visual processing. © 1992 Human Sciences Press, Inc.

Topographic mapping and source localization of the pattern reversal visual evoked magnetic response

The topography of the visual evoked magnetic response (VEMR) to pattern reversal stimulation was studied in four normal subjects using a single channel BTI magnetometer. VEMRs were recorded from 20 locations over the occipital scalp and the topographic distribution of the most consistent component (P100M) studied. A single dipole in a sphere model was fitted to the data. Topographic maps were similar when recorded two months apart on the same subject to the same stimulus. Half field (HF) stimulation elicited responses from sources on the medial surface of the calcarine fissure mainly in the contralateral hemisphere as predicted by the cruciform model. The full field (FF) responses to large checks were approximately the sum of the HF responses. However, with small checks, FF stimulation appeared to activate a different combination of sources than the two HFs. In addition, HF topography was more consistent between subjects than FF for small check sizes. Topographic studies of the VEMR may help to explain the analogous visual evoked electrical response and will be essential to define optimal recording positions for clinical applications.

Topographic study of the visual evoked magnetic response

Since the visual evoked potential to pattern reversal stimulation produces a paradoxical lateralisation of the major positive P100 component and since this paradoxical lateralisation is dependent on the stimulus parameters including check and field size, we have therefore, carried out a study of the magnetic response (VEMR) to a pattern reversal stimulus in four normal subjects using both full field and half field stimulation and two different check sizes.