Microstructural imaging of the human brain with a 'super-scanner': 10 key advantages of ultra-strong gradients for diffusion MRI

The key component of a microstructural diffusion MRI 'super-scanner' is a dedicated high-strength gradient system that enables stronger diffusion weightings per unit time compared to conventional gradient designs. This can, in turn, drastically shorten the time needed for diffusion encoding, increase the signal-to-noise ratio, and facilitate measurements at shorter diffusion times. This review, written from the perspective of the UK National Facility for In Vivo MR Imaging of Human Tissue Microstructure, an initiative to establish a shared 300 mT/m-gradient facility amongst the microstructural imaging community, describes ten advantages of ultra-strong gradients for microstructural imaging. Specifically, we will discuss how the increase of the accessible measurement space compared to a lower-gradient systems (in terms of Δ, b-value, and TE) can accelerate developments in the areas of 1) axon diameter distribution mapping; 2) microstructural parameter estimation; 3) mapping micro-vs macroscopic anisotropy features with gradient waveforms beyond a single pair of pulsed-gradients; 4) multi-contrast experiments, e.g. diffusion-relaxometry; 5) tractography and high-resolution imaging in vivo and 6) post mortem; 7) diffusion-weighted spectroscopy of metabolites other than water; 8) tumour characterisation; 9) functional diffusion MRI; and 10) quality enhancement of images acquired on lower-gradient systems. We finally discuss practical barriers in the use of ultra-strong gradients, and provide an outlook on the next generation of 'super-scanners'.

Global signal modulation of single-trial fMRI response variability: Effect on positive vs negative BOLD response relationship

In functional magnetic resonance imaging (fMRI), the relationship between positive BOLD responses (PBRs) and negative BOLD responses (NBRs) to stimulation is potentially informative about the balance of excitatory and inhibitory brain responses in sensory cortex. In this study, we performed three separate experiments delivering visual, motor or somatosensory stimulation unilaterally, to one side of the sensory field, to induce PBR and NBR in opposite brain hemispheres. We then assessed the relationship between the evoked amplitudes of contralateral PBR and ipsilateral NBR at the level of both single-trial and average responses. We measure single-trial PBR and NBR peak amplitudes from individual time-courses, and show that they were positively correlated in all experiments. In contrast, in the average response across trials the absolute magnitudes of both PBR and NBR increased with increasing stimulus intensity, resulting in a negative correlation between mean response amplitudes. Subsequent analysis showed that the amplitude of single-trial PBR was positively correlated with the BOLD response across all grey-matter voxels and was not specifically related to the ipsilateral sensory cortical response. We demonstrate that the global component of this single-trial response modulation could be fully explained by voxel-wise vascular reactivity, the BOLD signal standard deviation measured in a separate resting-state scan (resting state fluctuation amplitude, RSFA). However, bilateral positive correlation between PBR and NBR regions remained. We further report that modulations in the global brain fMRI signal cannot fully account for this positive PBR-NBR coupling and conclude that the local sensory network response reflects a combination of superimposed vascular and neuronal signals. More detailed quantification of physiological and noise contributions to the BOLD signal is required to fully understand the trial-by-trial PBR and NBR relationship compared with that of average responses.

Evidence that the negative BOLD response is neuronal in origin: a simultaneous EEG-BOLD-CBF study in humans

Unambiguous interpretation of changes in the BOLD signal is challenging because of the complex neurovascular coupling that translates changes in neuronal activity into the subsequent haemodynamic response. In particular, the neurophysiological origin of the negative BOLD response (NBR) remains incompletely understood. Here, we simultaneously recorded BOLD, EEG and cerebral blood flow (CBF) responses to 10 s blocks of unilateral median nerve stimulation (MNS) in order to interrogate the NBR. Both negative BOLD and negative CBF responses to MNS were observed in the same region of the ipsilateral primary sensorimotor cortex (S1/M1) and calculations showed that MNS induced a decrease in the cerebral metabolic rate of oxygen consumption (CMRO2) in this NBR region. The ∆CMRO2/∆CBF coupling ratio (n) was found to be significantly larger in this ipsilateral S1/M1 region (n=0.91±0.04, M=10.45%) than in the contralateral S1/M1 (n=0.65±0.03, M=10.45%) region that exhibited a positive BOLD response (PBR) and positive CBF response, and a consequent increase in CMRO2 during MNS. The fMRI response amplitude in ipsilateral S1/M1 was negatively correlated with both the power of the 8-13 Hz EEG mu oscillation and somatosensory evoked potential amplitude. Blocks in which the largest magnitude of negative BOLD and CBF responses occurred therefore showed greatest mu power, an electrophysiological index of cortical inhibition, and largest somatosensory evoked potentials. Taken together, our results suggest that a neuronal mechanism underlies the NBR, but that the NBR may originate from a different neurovascular coupling mechanism to the PBR, suggesting that caution should be taken in assuming the NBR simply represents the neurophysiological inverse of the PBR.

Spatial location and strength of BOLD activation in high-spatial-resolution fMRI of the motor cortex: a comparison of spin echo and gradient echo fMRI at 7 T

The increased blood oxygenation level-dependent contrast-to-noise ratio at ultrahigh field (7 T) has been exploited in a comparison of the spatial location and strength of activation in high-resolution (1.5  mm isotropic) gradient echo (GE) and spin echo (SE), echo planar imaging data acquired during the execution of a simple motor task in five subjects. SE data were acquired at six echo times from 30 to 55  ms. Excellent fat suppression was achieved in the SE echo planar images using slice-selective gradient reversal. Threshold-free cluster enhancement was used to define regions of interest (ROIs) containing voxels showing significant stimulus-locked signal changes from the GE and average SE data. These were used to compare the signal changes and spatial locations of activated regions in SE and GE data. T(2) and T(2)* values were measured, with means of 48.3 ± 1.1  ms and 36.5 ± 3.4  ms in the SE ROI. In addition, we identified a dark band in SE images of the motor cortex corresponding to a region in which T(2) and T(2)* were significantly lower than in the surrounding grey matter. The fractional SE signal change in the ROI was found to vary linearly as a function of TE, with a slope that was dependent on the particular ROI assessed: the mean ΔR(2) value was found to be 0.85 ± 0.11  s(-1) for the SE ROI and -0.37 ± 0.05  s(-1) for the GE ROI. The fractional signal change relative to the shortest TE revealed that the largest signal change occurred at a TE of 45  ms outside of the dark band. At this TE, the ratio of the fractional signal change in GE and SE data was found to be 0.48 ± 0.05. Phase maps produced from high-resolution GE images spanning the right motor cortex were used to identify veins. The GE ROI was found to contain 18% more voxels overlying the venous mask than the SE ROI.

Mapping human somatosensory cortex in individual subjects with 7T functional MRI

Functional magnetic resonance imaging (fMRI) is now routinely used to map the topographic organization of human visual cortex. Mapping the detailed topography of somatosensory cortex, however, has proven to be more difficult. Here we used the increased blood-oxygen-level-dependent contrast-to-noise ratio at ultra-high field (7 Tesla) to measure the topographic representation of the digits in human somatosensory cortex at 1 mm isotropic resolution in individual subjects. A "traveling wave" paradigm was used to locate regions of cortex responding to periodic tactile stimulation of each distal phalangeal digit. Tactile stimulation was applied sequentially to each digit of the left hand from thumb to little finger (and in the reverse order). In all subjects, we found an orderly map of the digits on the posterior bank of the central sulcus (postcentral gyrus). Additionally, we measured event-related responses to brief stimuli for comparison with the topographic mapping data and related the fMRI responses to anatomical images obtained with an inversion-recovery sequence. Our results have important implications for the study of human somatosensory cortex and underscore the practical utility of ultra-high field functional imaging with 1 mm isotropic resolution for neuroscience experiments. First, topographic mapping of somatosensory cortex can be achieved in 20 min, allowing time for further experiments in the same session. Second, the maps are of sufficiently high resolution to resolve the representations of all five digits and third, the measurements are robust and can be made in an individual subject. These combined advantages will allow somatotopic fMRI to be used to measure the representation of digits in patients undergoing rehabilitation or plastic changes after peripheral nerve damage as well as tracking changes in normal subjects undergoing perceptual learning.

Investigating the effect of blood susceptibility on phase contrast in the human brain

Recent work has shown a dramatic contrast between GM and WM in gradient echo phase images at high field (7 T). Although this contrast is key to the exploitation of phase in imaging normal and pathological tissue, its origin remains contentious. Several sources for this contrast have been considered including iron content, myelin, deoxy-hemoglobin, or water-macromolecule interactions. Here we quantify the contribution of intravascular dHb to the GM/WM contrast in the human brain at 7 T by modulating the susceptibility of the blood using a paramagnetic contrast agent. By carrying out high resolution, dynamic, gradient echo imaging before, during and after the injection of the contrast agent, we were able to follow the change in GM/WM phase contrast and to monitor simultaneously the susceptibility of the blood. Using these data in conjunction with the known susceptibility of venous blood we estimate the upper bound for the relative contribution of dHb in the vasculature to the measured GM/WM phase contrast to be 0.48 Hz for GM close to the pial surface, and 0.27 Hz for deeper GM. These values are up to 20% of the GM/WM phase difference observed in the human brain at 7 T. Furthermore, we found that the fractional blood volume differences required to account for the observed GM/WM phase contrast are 1.3% and 0.7% for GM close to the pial surface and for deeper GM, respectively.

Modeling and optimization of Look-Locker spin labeling for measuring perfusion and transit time changes in activation studies taking into account arterial blood volume

This work describes a new compartmental model with step-wise temporal analysis for a Look-Locker (LL)-flow-sensitive alternating inversion-recovery (FAIR) sequence, which combines the FAIR arterial spin labeling (ASL) scheme with a LL echo planar imaging (EPI) measurement, using a multireadout EPI sequence for simultaneous perfusion and T*(2) measurements. The new model highlights the importance of accounting for the transit time of blood through the arteriolar compartment, delta, in the quantification of perfusion. The signal expected is calculated in a step-wise manner to avoid discontinuities between different compartments. The optimal LL-FAIR pulse sequence timings for the measurement of perfusion with high signal-to-noise ratio (SNR), and high temporal resolution at 1.5, 3, and 7T are presented. LL-FAIR is shown to provide better SNR per unit time compared to standard FAIR. The sequence has been used experimentally for simultaneous monitoring of perfusion, transit time, and T*(2) changes in response to a visual stimulus in four subjects. It was found that perfusion increased by 83 +/- 4% on brain activation from a resting state value of 94 +/- 13 ml/100 g/min, while T*(2) increased by 3.5 +/- 0.5%.

Measurement of electric fields due to time-varying magnetic field gradients using dipole probes

The operation of dipole probes in measuring electric fields in conductive media exposed to temporally varying magnetic fields is discussed. The potential measured by the probe can be thought of as originating from two contributions to the electric field, namely the gradient of the scalar electric potential and the temporal derivative of the magnetic vector potential. Using this analysis, it is shown that the exact form of the wire paths employed when using electric field probes to measure the effects of temporally varying magnetic fields is very important and this prediction is verified via simple experiments carried out using different probe geometries in a cylindrical sample exposed to a temporally varying, uniform magnetic field. Extending this work, a dipole probe has been used to measure the electric field induced in a cylindrical sample by gradient coils as used in magnetic resonance imaging (MRI). Analytic solutions for the electric field in an infinite cylinder are verified by comparison with experimental measurements. Deviations from the analytic solutions of the electric field for the x-gradient coil due to the finite length of the sample cylinder are also demonstrated.

Magnetic-field-induced vertigo: a theoretical and experimental investigation

Vertigo-like sensations or apparent perception of movement are reported by some subjects and operators in and around high field whole body magnetic resonance body scanners. Induced currents (which modulate the firing rate of the vestibular hair cell), magneto-hydrodynamics (MDH), and tissue magnetic susceptibility differences have all been proposed as possible mechanisms for this effect. In this article, we examine the theory underlying each of these mechanisms and explore resulting predictions. Experimental evidence is summarised in the following findings: 30% of subjects display a postural sway response at a field-gradient product of 1 T(2)m(-1); a determining factor for experience of vertigo is the total unipolar integrated field change over a period greater than 1 s; the perception of dizziness is not necessarily related to a high value of the rate of change of magnetic field; eight of ten subjects reported sensations ranging from mild to severe when exposed to a magnetic field change of the order of 4.7 T in 1.9 s; no subjects reported any response when exposed to 50 ms pulses of dB/dt of 2 Ts(-1) amplitude. The experimental evidence supports the hypothesis that magnetic-field related vertigo results from both magnetic susceptibility differences between vestibular organs and surrounding fluid, and induced currents acting on the vestibular hair cells. Both mechanisms are consistent with theoretical predictions.

Theoretical optimization of multi-echo fMRI data acquisition

Simulations are used to optimize multi-echo fMRI data acquisition for detection of BOLD signal changes in this study. Optimal sequence design (echo times and sampling period (receiver bandwidth)) and the variation in sensitivity between tissues with different baseline T*(2) are investigated, taking into account the effects of physiological noise and non-exponential signal decay. In the case of a single echo, for normally distributed, uncorrelated noise, the results indicate that the sampling period should be made as long as possible (so as to produce an acceptable level of image distortion), up to a maximum sampling period of 3T*(2), (i.e. optimum TE = 1.5T*(2)). Combining the signal from multiple echoes using weighted summation improves the contrast-to-noise ratio (CNR), at a reduced optimum echo interval. If the BOLD effect causes a constant change in relaxation rate, DeltaR*(2), independent of the tissue R*(2), then a multi-echo acquisition causes considerable variation in sensitivity to BOLD signal changes with tissue T*(2), so that if the sequence is optimized for a target tissue T*(2) it will be more sensitive to BOLD signal changes in tissues with shorter T*(2) values. Fitting for DeltaR*(2) reduces the CNR, and when using this approach, the shortest echo time interval should be used, down to a limit of about 0.3T*(2), and as many echoes as possible within the constraints of TR or hardware limitations should be collected. It is also shown that the optimal sequence will remain optimum or close to optimum irrespective of whether there are physiological noise contributions.

Representations of pleasant and painful touch in the human orbitofrontal and cingulate cortices

The cortical areas that represent affectively positive and negative aspects of touch were investigated using functional magnetic resonance imaging (fMRI) by comparing activations produced by pleasant touch, painful touch produced by a stylus, and neutral touch, to the left hand. It was found that regions of the orbitofrontal cortex were activated more by pleasant touch and by painful stimuli than by neutral touch and that different areas of the orbitofrontal cortex were activated by the pleasant and painful touches. The orbitofrontal cortex activation was related to the affective aspects of the touch, in that the somatosensory cortex (SI) was less activated by the pleasant and painful stimuli than by the neutral stimuli. This dissociation was highly significant for both the pleasant touch (P < 0.006) and for the painful stimulus (P < 0.02). Further, it was found that a rostral part of the anterior cingulate cortex was activated by the pleasant stimulus and that a more posterior and dorsal part was activated by the painful stimulus. Regions of the somatosensory cortex, including SI and part of SII in the mid-insula, were activated more by the neutral touch than by the pleasant and painful stimuli. Part of the posterior insula was activated only in the pain condition and different parts of the brainstem, including the central grey, were activated in the pain, pleasant and neutral touch conditions. The results provide evidence that different areas of the human orbitofrontal cortex are involved in representing both pleasant touch and pain, and that dissociable parts of the cingulate cortex are involved in representing pleasant touch and pain.

Optimizing the sequence parameters for double-quantum CRAZED imaging

The evolution of magnetization during repeated application of the double-quantum-(DQ)-CRAZED sequence is analyzed, with the aim of identifying sequence parameters that maximize sensitivity to signal produced by the distant dipole field (DDF). Phase cycling schemes that allow cancellation of signals following undesired coherence pathways are also described. Simulations and imaging experiments carried out at 3 T on phantoms and the human head were used to verify the analysis. The results indicate that in the absence of phase cycling, the DDF-related signal-to-noise ratio (SNR) per unit time is maximized using TR=2.05 T1, along with values of the RF flip angles (alpha approximately 90 degrees and beta approximately 60 degrees ), and echo time (TE=T2) that have previously been shown to maximize the DDF-related signal at long TR. However, with TR=2.05 T1 there can also be a significant signal contribution due to stimulated echo effects (up to 40% of the signal for water at 3 T and TE=70 ms). Using a two-step phase cycle, the stimulated echo signal is eliminated and the maximum SNR per unit time occurs for TR=2.76 T1. It is also demonstrated that sensitivity to signal changes caused by small variations in T2 is maximized by setting TE=2T2.

Mapping the absolute value of M0 using dipolar field effects

The ability to map the spatial variation of the absolute, rather than the relative value of the equilibrium magnetization could be advantageous in many areas of NMR. However, direct measurement of M(0) is usually difficult because of the multiparametric dependence of the NMR signal. Here we propose a technique for mapping the spatial variation of the absolute value of M(0), independent of relaxation weighting and flip angle calibration. This method, which works best at high field strengths, is based on the effect of the dipolar field due to the nuclear magnetization that is normally neglected in liquid-state NMR. The experimental implementation of this sequence at 3.0 T is described, and its initial application to the measurement of the water content of brain tissue is outlined.

Electric fields induced in a spherical volume conductor by temporally varying magnetic field gradients

A homogeneous spherical volume conductor is used as a model system for the purpose of calculating electric fields induced in the human head by externally applied time-varying magnetic fields. We present results for the case where magnetic field gradient coils, used in magnetic resonance imaging (MRI), form the magnetic field, and we use these data to put limits on the rates of gradient change with time needed to produce nerve stimulation. The electric field is calculated analytically for the case of ideal longitudinal and transverse linear field gradients. We also show results from computer calculations yielding the electric field maps in a sphere when the field gradients are generated by a real MRI gradient coil set. In addition, the effect of shifting the sphere within each gradient coil volume is investigated. Numerical analysis shows similar results when applied to a model human head.

Imaging the long-range dipolar field in structured liquid state samples

We describe imaging experiments in which the pattern of the dipolar field generated by spatially modulated nuclear magnetization is directly visualized in simply structured phantoms. Two types of experiment have been carried out at 11.7 T using (1)H NMR signals. In the first, the field from a single spin species is imaged via its own NMR signal. In the second, the NMR signal from one spin species is used to image the field generated by a second species. The field patterns measured in these experiments correspond well with those calculated using simple theoretical expressions for the dipolar field. The results also directly demonstrate the spatial sensitivity of the signal generated using dipolar field effects, indicating that the range of the field depends upon the inverse of the spatial frequency with which the magnetization is modulated.

Single-shot T2(*) measurement to establish optimum echo time for fMRI: studies of the visual, motor, and auditory cortices at 3.0 T

The signal change in fMRI is dependent on the echo time and the rate of decay of transverse magnetization. The latter factor may vary across regions of the brain as a result of variations in field homogeneity. Previous measurements of the signal change with echo time have generally employed relatively slow multi-echo techniques, which may be sensitive to movement and habituation effects. Here a fast T(2)(*) measurement technique, involving the generation of six low-resolution echo planar images from a single FID, is described, and its use in the evaluation of the optimum echo time for visual, motor, and auditory fMRI experiments at 3.0 T is outlined.

Cortical responses to single mechanoreceptive afferent microstimulation revealed with fMRI

The technique of intraneural microneurography/microstimulation has been used extensively to study contributions of single, physiologically characterized mechanoreceptive afferents (MRAs) to properties of somatosensory experience in awake human subjects. Its power as a tool for sensory neurophysiology can be greatly enhanced, however, by combining it with functional neuroimaging techniques that permit simultaneous measurement of the associated CNS responses. Here we report its successful adaptation to the environment of a high-field MR scanner. Eight median-nerve MRAs were isolated and characterized in three subjects and microstimulated in conjunction with fMRI at 3.0 T. Hemodynamic responses were observed in every case, and these responses were robust, focal, and physiologically orderly. The combination of fMRI with microstimulation will enable more detailed studies of the representation of the body surface in human somatosensory cortex and further studies of the relationship of that organization to short-term plasticity in the human SI cortical response to natural tactile stimuli. It can also be used to study many additional topics in sensory neurophysiology, such as CNS responses to additional classes of afferents and the effects of stimulus patterning and unimodal/crossmodal attentional manipulations. Finally, it presents unique opportunities to investigate the basic physiology of the BOLD effect and to compare the operating characteristics of fMRI and EEG as human functional neuroimaging modalities in an unusually specific and well-characterized neurophysiological setting.

Representation of pleasant and aversive taste in the human brain

In this study, the representation of taste in the orbitofrontal cortex was investigated to determine whether or not a pleasant and an aversive taste have distinct or overlapping representations in this region. The pleasant stimulus used was sweet taste (1 M glucose), and the unpleasant stimulus was salt taste (0.1 M NaCl). We used an ON/OFF block design in a 3T fMRI scanner with a tasteless solution delivered in the OFF period to control for somatosensory or swallowing-related effects. It was found that parts of the orbitofrontal cortex were activated (P < 0.005 corrected) by glucose (in 6/7 subjects) and by salt (in 6/7 subjects). In the group analysis, separate areas of the orbitofrontal cortex were found to be activated by pleasant and aversive tastes. The involvement of the amygdala in the representation of pleasant as well as aversive tastes was also investigated. The amygdala was activated (region of interest analysis, P < 0.025 corrected) by the pleasant taste of glucose (5/7 subjects) as well as by the aversive taste of salt (4/7 subjects). Activation by both stimuli was also found in the frontal opercular/insular (primary) taste cortex. We conclude that the orbitofrontal cortex is involved in processing tastes that have both positive and negative affective valence and that different areas of the orbitofrontal cortex may be activated by pleasant and unpleasant tastes. We also conclude that the amygdala is activated not only by an affectively unpleasant taste, but also by a taste that is affectively pleasant, thus providing evidence that the amygdala is involved in effects produced by positively affective as well as by negatively affective stimuli.

Analytic calculations of the E-fields induced by time-varying magnetic fields generated by cylindrical gradient coils

Analytic expressions which allow the direct calculation of the electric field generated inside an infinite conducting cylinder by varying the current through the wires of any cylindrical coil are presented. These expressions provide some general insight into the spatial characteristics of the electric field generated inside the body by switched gradients and can be used to evaluate the locations where nerve stimulation by rapid gradient switching is likely to occur. They may also be employed at the design stage to produce gradient coils which can provide higher gradient switching rates without causing nerve stimulation. Using these expressions the electric field patterns produced by transverse and longitudinal, whole-body gradient coils were calculated. Example data are presented along with the associated magnetic field patterns. The effect on the induced electric field pattern of varying the body size and the size of the region of gradient linearity was explored.

Lip-reading ability and patterns of cortical activation studied using fMRI

Lip-reading is a complex cognitive skill with large individual differences in performance. The basis of these individual differences remains poorly understood. Functional magnetic resonance imaging (fMRI) techniques allows brain activation accompanying complex cognitive activities to be studied noninvasively. In the present paper, fMRI was used to study the patterns of cortical activation that occur during the silent lip-reading of connected speech and to investigate whether there are detectable differences in activation between subjects with widely differing lip-reading abilities. From a cohort of 26 volunteers, nine subjects who fell into three distinct lip-reading ability groups were selected. Brain activation was measured in two conditions: an experimental condition where subjects attempted to lip-read sentences; and a baseline condition where subjects passively viewed a static image of a talker's face. Relative to the baseline condition, lip-reading induced activation in several cortical areas, including the auditory cortices, despite the lack of an auditory component to the task. In comparison to the better two groups of lip-readers, subjects in the poorest group displayed significantly less activation in superior and middle temporal gyrus, but not inferior temporal gyrus. These preliminary results justify more extensive investigations of the cortical basis of individual differences in lip-reading.

Sensory-specific satiety-related olfactory activation of the human orbitofrontal cortex

When a food is eaten to satiety, its reward value decreases. This decrease is usually greater for the food eaten to satiety than for other foods, an effect termed sensory-specific satiety. In an fMRI investigation it was shown that for a region of the orbitofrontal cortex the activation produced by the odour of the food eaten to satiety decreased, whereas there was no similar decrease for the odour of a food not eaten in the meal. This effect was shown both by a voxel-wise SPM contrast (p <0.05 corrected) and an ANOVA performed on the mean percentage change in BOLD signal in the identified clusters of voxels (p <0.006). These results show that activation of a region of the human orbitofrontal cortex is related to olfactory sensory-specific satiety.

Spatial and temporal distribution of solutes in the developing carrot taproot measured at single-cell resolution

The time-course and spatial distribution of sugars and ions in carrot (Daucus carota L.) was studied at fine resolution using single cell (SiCSA) and tissue analysis. Four phases of osmolyte accumulation in the taproot were identified: an amino acid (germination) phase, when internal sources of amino acids provide seedlings with osmotica; an ion phase, when inorganic and organic ions were the main solutes; a hexose phase, when concentrations of glucose and fructose sharply increased and reached their maximum; and a sucrose phase, when sucrose became the major solute. Spatial distribution of sugar in taproot cells showed a general trend of highest concentration on both sides of the vascular cambium (some 200 mM sucrose, 150 mM glucose) and a minimum in the pith (some 100 mM sucrose, 60 mM glucose) and in periderm. Electrolytes (e.g. potassium) followed a distribution generally reciprocal to that of sugars; minimum in the tissue adjacent to the cambium (some 10 mM) and maximum in the pith and periderm (some 60-100 mM). The cambial cells contained unexpectedly low concentrations of sugars and potassium. These spatial and temporal patterns indicate that amino acids, other electrolytes and sugars are interchangeable in the tissue osmotic balance. The nature of the solute is developmentally determined both temporally and spatially. During the accumulation of electrolytes following the initial amino acid phase, osmotic pressure to 420 mosmol kg-1 rises and then remains constant despite large changes in the concentration of individual solutes. This indicates that osmotic pressure is regulated independently of the individual concentrations of solutes.

fMRI of the responses to vibratory stimulation of digit tips

Three studies were carried out to assess the applicability of fMRI at 3.0 T to analysis of vibrotaction in humans. A novel piezoelectric device provided clean sinusoidal stimulation at 80 Hz, which was initially applied in separate runs within a scanning session to digits 2 and 5 of the left hand in eight subjects, using a birdcage RF (volume) coil. Significant clusters of activation were found in the primary somatosensory cortex (SI), the secondary somatosensory cortex (SII), subcentral gyrus, the precentral gyrus, posterior insula, posterior parietal regions (area 5), and the posterior cingulate. Digit separation in SI was possible in all subjects and the activation sites reflected the known lateral position of the representation of digit 2 relative to that of digit 5. A second study carried out in six additional subjects using a surface coil, replicated the main contralateral activation patterns detected in study one and further improved the discrimination of the digits in SI. Significant digit separation was also found in SII and in the posterior insula. A third study to investigate the frequency dependence of the response focused on the effect of an increase in vibrotactile frequency from 30 to 80 Hz, with both frequencies applied to digit 2 during the same scanning session in four new subjects. A significant increase in the number of pixels activated within both SII and the posterior insula was found, while the number of pixels activated in SI declined. No significant change in signal intensity with frequencies was found in any of the activated areas.

Sensory-specific satiety-related olfactory activation of the human orbitofrontal cortex

When a food is eaten to satiety, its reward value decreases. This decrease is usually greater for the food eaten to satiety than for other foods, an effect termed sensory-specific satiety. In an fMRI investigation it was shown that for a region of the orbitofrontal cortex the activation produced by the odour of the food eaten to satiety decreased, whereas there was no similar decrease for the odour of a food not eaten in the meal. This effect was shown both by a voxel-wise SPM contrast (p<0.05 corrected) and an ANOVA performed on the mean percentage change in BOLD signal in the identified clusters of voxels (p<0.006). These results show that activation of a region of the human orbitofrontal cortex is related to olfactory sensory-specific satiety.

Detecting activations in event-related fMRI using analysis of variance

The most common design of a functional MRI (fMRI) experiment is a block design. The use of rapid imaging, however, and carefully designed paradigms makes the separation of cognitive events possible. Such experiments make use of event-related paradigms, in which a task involving several cognitive processes is repeated. In analyzing data from such experiments, existing methods often prove inadequate, because the prediction of the exact shape or timing of the time course is difficult. Here we present an analysis of variance (ANOVA) method for analyzing fMRI data that does not require any assumptions about the shape of the activation time course. Consequently, this method can simultaneously detect brain areas showing a variety of stimulus-locked time courses in the same experiment. The utility of this technique is demonstrated by the analysis of data from two event-related paradigms in which regions of activation are detected that correspond to a variety of distinct neural processes, yielding significantly different temporal signal changes. Magn Reson Med 42:1117-1122, 1999.

Fetal brain activity demonstrated by functional magnetic resonance imaging

Functional magnetic resonance imaging was used to study fetal brain activity. This activity was in response to an auditory stimulus.

Resolution in high field echo planar microscopy

The application of echo planar imaging to NMR microscopy offers a temporal resolution unparalleled by other techniques. However, a major difficulty in imaging at the high field strengths used for microscopy is the effect of local field inhomogeneities caused by magnetic susceptibility effects. This can give rise to both image distortion and signal loss. In addition, the effect of diffusion in the presence of the large imaging gradients gives rise to a broadening of the point spread function and hence loss of true resolution. We compare the sensitivity of two techniques, MBEST and PEPI, to both of these effects. Analytic expressions for the signal in each echo of the two sequences are developed, and the point spread functions for the two techniques are calculated. Using PEPI, we have been able to produce images with an in-plane resolution of 50 micrometer from a single free induction decay. This technique has been extended to three dimensions allowing the generation of 64(3) images with an isotropic resolution of 80 micrometer.

Analytic approach to the design of transverse gradient coils with co-axial return paths

Transverse gradient coils with co-axial return paths offer reduced acoustic noise compared with standard cylindrical gradient coils, due to local force balancing, and can also easily be made to have a length to diameter ratio that is less than one. Analytic expressions for the magnetic field and vector potential generated by this type of coil are described here, along with a formula for calculating the coil inductance. It is shown that these expressions allow the implementation of powerful analytic methods of coil design, as well as the incorporation of active magnetic screening. It is also demonstrated how the mathematics specifies the best parameters to use when designing coils with small numbers of elements. A head gradient coil for use at 3.0 T has been designed using the analytic approach described here. The process of coil design and construction is outlined and the performance of the coil in comparison with a similar standard cylindrical coil is described.

The representation of pleasant touch in the brain and its relationship with taste and olfactory areas

Although there has been much investigation of brain pathways involved in pain, little is known about the brain mechanisms involved in processing somatosensory stimuli which feel pleasant. Employing fMRI it was shown that pleasant touch to the hand with velvet produced stronger activation of the orbitofrontal cortex than affectively neutral touch of the hand with wood. In contrast, the affectively neutral but more intense touch produced more activation of the primary somatosensory cortex than the pleasant stimulus. This indicates that part of the orbitofrontal cortex is concerned with representing the positively affective aspects of somatosensory stimuli, and in further experiments it was shown that this orbitofrontal area is different from that activated by taste and smell. The finding that three different primary or unlearned types of reinforcer (touch, taste, and smell) are represented in the orbitofrontal cortex helps to provide a firm foundation for understanding the neural basis of emotions, which can be understood in terms of states elicited by stimuli which are rewarding or punishing.

Automatic compensation of motion artifacts in MRI

Patient motion during the acquisition of a magnetic resonance image can cause blurring and ghosting artifacts in the image. This paper presents a new post-processing strategy that can reduce artifacts due to in-plane, rigid-body motion in times comparable to that required to re-scan a patient. The algorithm iteratively determines unknown patient motion such that corrections for this motion provide the best image quality, as measured by an entropy-related focus criterion. The new optimization strategy features a multi-resolution approach in the phase-encode direction, separate successive one-dimensional searches for rotations and translations, and a novel method requiring only one re-gridding calculation for each rotation angle considered. Applicability to general rigid-body in-plane rotational and translational motion and to a range of differently weighted images and k-space trajectories is demonstrated. Motion artifact reduction is observed for data from a phantom, volunteers, and patients.

Self-diffusion and molecular mobility in PVA-based dissolution-controlled systems for drug delivery

Nuclear magnetic resonance (NMR) microscopy has been used to monitor the hydration of poly(vinyl alcohol) (PVA) samples of varying molecular weight. One-dimensional profiles weighted to predominantly show the variation of water concentration were acquired every 3 min during the first 30 min of hydration and subsequently at 1 and 2 h. Diffusion-weighted profiles obtained after 30 min and 1 and 2 h were used to calculate the spatial variation of the water self-diffusion coefficient. The resulting data provide supporting evidence for the hypothesis that phenomena such as reptation are important near the glassy/rubbery interface of polymers during dissolution, while the diffusion gradually changes to Zimm type near the rubbery/solvent interface.

Strong gradients for spatially resolved diffusion measurements

A new multilayer approach to gradient coil design, which allows the production of very strong gradient coils with reasonable resistance and consequent power dissipation, has been developed. Using this approach we have designed and built a strong z-gradient coil that will accommodate vertically mounted samples contained in 5-mm nuclear magnetic resonance tubes. The coil has an efficiency of 1.73 Tm-1A-1, an inductance of 49 microH, and a resistance of 1.8 omega, with a homogeneous volume consisting of a central cylinder of 4.5-mm length and diameter. This coil has been used to monitor the diffusion of water in Nylon 6.6 at room temperature, during desorption. This system is difficult to monitor via nuclear magnetic resonance (NMR), because the diffusion coefficients are typically less than 10(-13) m2s-1, while the T2 relaxation time is less than 1 ms even when the sample is fully saturated. The resulting measurements show a strong concentration dependence of the T2 relaxation time and self-diffusion coefficient of the absorbed water. The measured concentration profiles are consistent with a Fickian diffusion process with a concentration-dependent diffusion coefficient. The measured self-diffusion values are in reasonable agreement with those inferred from the variation of the concentration profiles as a function of time, using the one-dimensional Fickian diffusion equation.

Multilayer Gradient Coil Design

In standard cylindrical gradient coils consisting of wires wound in a single layer, the rapid increase in coil resistance with efficiency is the limiting factor in achieving very large magnetic field gradients. This behavior results from the decrease in the maximum usable wire diameter as the number of turns is increased. By adopting a multilayer design in which the coil wires are allowed to spread out into multiple layers wound at increasing radii, a more favorable scaling of resistance with efficiency is achieved, thus allowing the design of more powerful gradient coils with acceptable resistance values. By extending the theory used to design standard cylindrical gradient coils, we have developed mathematical expressions which allow the design of multilayer coils, and the evaluation of their performance. These expressions have been used to design a four-layer, z-gradient coil of 8 mm inner diameter, which has an efficiency of 1.73 Tm-1 A-1, a resistance of 1.8 Omega, and an inductance of 50 µH. This coil produces a gradient which deviates from linearity by less than 5% within a central cylindrical region of 4.5 mm length and 4.5 mm diameter. A coil has been constructed from this design and tested in simple imaging and pulsed gradient spin echo experiments. The resulting data verify the predicted coil performance, thus demonstrating the advantages of using multilayer coils for experiments requiring very large magnetic field gradients. Copyright 1998 Academic Press.

Functional magnetic resonance imaging of single motor events reveals human presupplementary motor area

Conventional functional imaging paradigms use periods of repetitive task performance to generate sustained functional signal changes. We have developed a technique of imaging the small, transient signal changes that occur after single cognitive events. The technique uses echo-planar imaging at 3 T to generate functional images of the whole brain with a temporal resolution of 3 seconds. It uses a signal averaging technique to create time sweeps of functional activity. After a single cognitive event, widely distributed patterns of brain activation can be detected and their time course measured. This technique enables the individual cognitive tasks that constitute a paradigm to be analyzed separately and compared. We describe the application of this new technique to separate the cognitive elements in a simple "go/no-go" motor paradigm. Comparison of activation patterns during "go" and "no-go" responses reveals hierarchical subdivision of the medial premotor cortex into an anterior region (presupplementary motor area) involved in movement decision making and a posterior region (supplementary motor area proper) directly involved in motor execution.

A modified imaging sequence for accurate T2 measurements using NMR microscopy

A modified spin-echo pulse sequence is described that enables accurate T2 measurements to be made in NMR microimaging experiments. The modified sequence eliminates cumulative diffusion losses that lead to an underestimation of the T2 relaxation time using conventional spin-echo pulse sequences. The approach is theoretically justified and confirmed in comparative experiments on phantoms.

Active acoustic screening: reduction of noise in gradient coils by Lorentz force balancing

We have designed and constructed a quiet gradient set with restricted access for the combined purposes of evaluating the principles of active acoustic screening, recently introduced by Mansfield, Glover, and Bowtell, and for EPI studies of the head at 3.0 T. The design utilizes the return paths of the conductors in a closed arc loop arrangement to eliminate net Lorentz forces thereby attenuating acoustic noise especially at low frequency. This design should significantly reduce the dangers to patients of high noise levels, especially in high field magnetic resonance imaging systems.

NMR microscopy of hydrating hydrophilic matrix pharmaceutical tablets

NMR microscopy has been used to monitor the formation of the gel layer in hydrating hydrophilic polymer tablets. Such tablets are used in the controlled delivery of drugs, where it has been found that the rate and extent of the swelling of the outer gel layer critically influences the kinetics of drug release. Tablets were hydrated in distilled water at 37 degrees C and then imaged at discrete time intervals using a 500 MHz microscope. The growth of the gel layer was clearly observed in time sequences of radial and axial sections. Axial images showed some interesting dimensional changes, with the gel at the flat surface of the tablet developing a concave shape. This is probably a reflection of the occurrence of uni-axial stress relaxation as hydration proceeds. Diffusion- and T2-weighted images provided evidence that the water in the gel layer is more strongly bound close to the dry core of the tablet than at the more fully hydrated outer surface. In images of tablets containing diclofenac, disruption of the gel layer was shown to occur primarily from the flat surfaces of the tablet, whilst the distribution of particles could be seen in tablets doped with insoluble calcium phosphate.

Magnetic resonance imaging: applications of novel methods in studies of porous media

NMR imaging is finding broad applications in nonbiological areas including the study of fluid flow and fluid ingress in porous media. The porous media include at the one end mineral rocks and various building materials through various solid plastic materials to foodstuffs at the other end of the spectrum. The fluids within these various media range from crude oil and water mixtures, and water itself, to a range of organic solvents in plastic materials. This paper is concerned with the flow and ingress of water through Bentheimer sandstone and Ninian reservoir specimens, and also in solid nylon blocks.