Search

PURPOSE: The goal of this study was to measure diffusion signals within cerebral cortex using the line-scan technique to achieve extremely high resolution in the radial direction (i.e., perpendicular to the cortical surface) and to demonstrate the utility of these measurements for investigating laminar architecture in the living human brain. METHODS: Line-scan diffusion data with 250–500 micron radial resolution were acquired at 7T on eight healthy volunteers, with each line prescribed perpendicularly to primary somatosensory cortex (S1) and primary motor cortex (M1). Apparent diffusion coefficients (ADCs), fractional anisotropy (FA) values and radiality indices were measured as a function of cortical depth. RESULTS: In the deep layers of S1, we found evidence for high anisotropy and predominantly tangential diffusion, with low anisotropy observed in superficial S1. In M1, moderate anisotropy and predominantly radial diffusion was seen at almost all cortical depths. These patterns were consistent across subjects and were conspicuous without averaging data across different locations on the cortical sheet. CONCLUSION: Our results are in accord with the myeloarchitecture of S1 and M1, known from prior histology studies: in S1, dense bands of tangential myelinated fibers run through the deep layers but not the superficial ones, and in M1, radial myelinated fibers are prominent at most cortical depths. This work therefore provides support for the idea that high-resolution diffusion signals, measured with the line-scan technique and receiving a boost in signal-to-noise ratio at 7T, may serve as a sensitive probe of in vivo laminar architecture.

Understanding how and why MR signals and their associated relaxation rates vary with cortical depth could ultimately enable the noninvasive investigation of the laminar architecture of cerebral cortex in the living human brain. However, cortical gray matter is typically only a few millimeters thick, making it challenging to sample many cortical depths with the voxel sizes commonly used in MRI studies. Line-scan techniques provide a way to overcome this challenge and here we implemented a novel line-scan GESSE pulse sequence that allowed us to measure irreversible and reversible transverse relaxation rates—R(2) and R(2)’, respectively—with extremely high resolution (250 μm) in the radial direction, perpendicular to the cortical surface. Eight healthy human subjects were scanned at 7T using this sequence, with primary visual cortex (V1) targeted in three subjects and primary motor (M1) and somatosensory cortex (S1) targeted in the other five. In all three cortical areas, a peak in R(2) values near the central depths was seen consistently across subjects—an observation that has not been made before, to our knowledge. On the other hand, no consistent pattern was apparent for R(2)’ values as a function of cortical depth. The intracortical R(2) peak reported here is unlikely to be explained by myelin content or by deoxyhemoglobin in the microvasculature; however, this peak is in accord with the laminar distribution of non-heme iron in these cortical areas, known from prior histology studies. Obtaining information about tissue microstructure via measurements of transverse relaxation (and other quantitative MR contrast mechanisms) at the extremely high radial resolutions achievable through the use of line-scan techniques could therefore bring us closer to being able to perform “in vivo histology” of the cerebral cortex.

Purpose To compare transverse relaxation rates of brain metabolites estimated from single-TE PRESS acquisitions with more conventionally derived rates estimated from multiple-TE PRESS acquisitions. Methods Single-voxel (8 mL) PRESS data within white matter from 6 subjects were acquired at five different TEs. Transverse relaxation rates R2 of N-acetylaspartate, creatine, and choline were estimated from a single TE using full versus right-side-only sampling of the echo. These R2 values were compared with R2Hahn values obtained from the multiple-TE PRESS acquisitions. Results Following baseline subtraction and RMS weighting, interindividual mean R2 values from TE = 288 ms magnitude spectra for choline, creatine, and N-acetylaspartate were highly correlated with respective R2Hahn values (r2 = 0.99). Paired individual measurements at this TE showed less correlation (r2 = 0.48), primarily due to the N-acetylaspartate resonance. Using TE = 360 ms data for N-acetylaspartate and 288 ms for choline and creatine resulted in an improved correlation coefficient (r2 = 0.80). The average absolute intra-individual differences in the estimated R2 s between single-TE and Hahn method was 9.6 ± 7.7%. Conclusion For the major brain metabolite singlets, R2Hahn values showed correlations with more fragile measurements of R2 from a single TE that are worthy of interest. Because the left side of long-TE spin echoes is available "for free" from an acquisition perspective, and although the single-TE method for estimating R2 values is associated with lower precision, the reduction in scan time may be clinically helpful.

The goal of this study was to test a new formalism for extracting reversible and irreversible transverse relaxation rates from resonances within typical proton muscle spectra using only a single spin echo as acquired with routine single-voxel, point-resolved echo spectroscopy (PRESS) acquisitions. Single-voxel, non-water-suppressed PRESS acquisitions within the calf muscles of four healthy subjects were performed at 1.5 T using six echo times ranging from 30 to 576 ms. Novel transverse relaxation analyses of water, choline, creatine, and lipid resonances were performed based upon the disparate relaxation sensitivities of the left versus the right sides of spectroscopically sampled spin echoes. Irreversible and reversible transverse relaxation rates R2 and R2 ' were extracted for water, metabolites, and lipids using echo times of 288 ms and longer. The R2 values so obtained were compared with more conventional "gold standard" Hahn values, R2Hahn , evaluated from the echo-time dependence of spectral peak areas generated from right-side sampling alone. Water resonances displayed biexponential Hahn signal decays, consistent with observations of decreasing R2 values with increasing echo time via the new approach. Choline and creatine resonances displayed monoexponential echo-time decays, with R2Hahn values in reasonable agreement with R2 values obtained from the single-echo analyses at the longer echo times. Lipid methylene and methyl R2 values extracted from the new approach were also in reasonable accord with R2Hahn values. Further validation of the technique was provided through PRESS acquisitions on a water phantom to which various levels of gadolinium were added in order to manipulate transverse relaxation rates, yielding excellent agreement between water-resonance R2Hahn and single-echo R2 values. In summary, this work demonstrates the feasibility of measuring reversible and irreversible transverse relaxation rates for individual spectral peaks from single-echo PRESS acquisitions, enabling a reduction in overall scan time relative to the use of multiple acquisitions with varying echo time.

Use of spectroscopically-acquired spin echoes typically involves Fourier transformation of the right side of the echo while largely neglecting the left side. For sufficiently long echo times, the left side may have enough spectral resolution to offer some utility. Since the acquisition of this side is “free”, we deemed it worthy of attention and investigated the spectral properties and information content of this data. Theoretical expressions for left- and right-side spectra were derived assuming Lorentzian frequency distributions. For left-side spectra, three regimes were identified based upon the relative magnitudes of reversible and irreversible transverse relaxation rates, R 2′ and R 2, respectively. Point-resolved spectroscopy (PRESS) data from muscle, fat deposit and bone marrow were acquired at 1.5 T to test aspects of the theoretical expressions. For muscle water or methylene marrow resonances, left-side signals were substantially or moderately larger than right-side signals but were similar in magnitude for muscle choline and creatine resonances. Left- versus right-side spectral-peak amplitude ratios depend sensitively on the relative values of R 2 and R 2′ , which can be estimated given this ratio and a right-side linewidth measurement. Left-side spectra can be used to augment signal-to-noise and to estimate spectral R 2 and R 2′ values under some circumstances.

To investigate the heating of EEG electrodes during magnetic resonance imaging (MRI) scans and to better understand the underlying physical mechanisms with a focus on the antenna effect.Gold cup and conductive plastic electrodes were placed on small watermelons with fiberoptic probes used to measure electrode temperature changes during a variety of 1.5T and 3T MRI scans. A subset of these experiments was repeated on a healthy human volunteer.The differences between gold and plastic electrodes did not appear to be practically significant. For both electrode types, we observed heating below 4°C for straight wires whose lengths were multiples of ½ the radiofrequency (RF) wavelength and stronger heating (over 15°C) for wire lengths that were odd multiples of ¼ RF wavelength, consistent with the antenna effect.The antenna effect, which has received little attention so far in the context of EEG-MRI safety, can play as significant a role as the loop effect (from electromagnetic induction) in the heating of EEG electrodes, and therefore wire lengths that are odd multiples of ¼ RF wavelength should be avoided. These results have important implications for the design of EEG electrodes and MRI studies as they help to minimize the risk to patients undergoing MRI with EEG electrodes in place.

A biomarker of cancer aggressiveness, such as hypoxia, could substantially impact treatment decisions in the prostate, especially radiation therapy, by balancing treatment morbidity (urinary incontinence, erectile dysfunction, etc.) against mortality. R2 (*) mapping with Mono-Exponential (ME) decay modeling has shown potential for identifying areas of prostate cancer hypoxia at 1.5T. However, Gaussian deviations from ME decay have been observed in other tissues at 3T. The purpose of this study is to assess whether gradient-echo signal decays are better characterized by a standard ME decay model, or a Gaussian Augmentation of the Mono-Exponential (GAME) decay model, in the prostate at 3T. Multi-gradient-echo signals were acquired on 20 consecutive patients with a clinical suspicion of prostate cancer undergoing MR-guided prostate biopsies. Data were fitted with both ME and GAME models. The information contents of these models were compared using Akaike's information criterion (second order, AICC ), in skeletal muscle, the prostate central gland (CG), and peripheral zone (PZ) regions of interest (ROIs). The GAME model had higher information content in 30% of the prostate on average (across all patients and ROIs), covering up to 67% of cancerous PZ ROIs, and up to 100% of cancerous CG ROIs (in individual patients). The higher information content of GAME became more prominent in regions that would be assumed hypoxic using ME alone, reaching 50% of the PZ and 70% of the CG as ME R2 (*) approached 40 s(-1) . R2 (*) mapping may have important applications in MRI; however, information lost due to modeling could mask differences in parameters due to underlying tissue anatomy or physiology. The GAME model improves characterization of signal behavior in the prostate at 3T, and may increase the potential for determining correlates of fit parameters with biomarkers, for example of oxygenation status.

Purpose To assess whether R2* mapping with a standard Mono-Exponential (ME) or a Gaussian Augmentation of the Mono-Exponential (GAME) decay model better characterizes gradient-echo signal decays in gynecological cancers after external beam radiation therapy at 3T, and evaluate implications of modeling for non-invasive identification of intra-tumoral hypoxia.

The goal of this study was to measure irreversible and reversible transverse relaxation rates in the globus pallidus and putamen at 7 T, and to use these rates to make inferences about the sub-voxel structure of iron and calcification deposits. Gradient Echo Sampling of a Spin Echo (GESSE) data were acquired at 7 T on eighteen volunteers spanning a large range of ages (23-85 years), with calcifications in the globus pallidus incidentally observed in one volunteer. Maps of transverse relaxation rates were derived from the GESSE data, and the mean value of these rates in globus pallidus and putamen was estimated for each volunteer. Both irreversible and reversible transverse relaxation rates increased with the expected age-dependent iron content in these structures, except for the individual with calcifications for whom extremely large reversible relaxation rates but normal irreversible relaxation rates were found in the globus pallidus. Given the sensitivity of irreversible and reversible transverse relaxation rates to microscopic and mesoscopic field variations, respectively, our findings suggest that joint consideration of these rates may yield information not only about the amount of iron and calcification deposited in the brain, but also about the sub-voxel structure of these deposits, perhaps revealing certain aspects of their geometry and cellular distribution.

The aim of this study was to demonstrate that gradient-echo sampling of single spin echoes can be used to isolate the signal from trabecular bone marrow, with high-quality segmentation and surface reconstructions resulting from the application of simple post-processing strategies. Theoretical expressions of the time-domain single-spin-echo signal were used to simulate signals from bone marrow, non-bone fatty deposits and muscle. These simulations were compared with and used to interpret signals obtained by the application of the gradient-echo sampling of a spin-echo sequence to image the knee and surrounding tissues at 1.5 T. Trabecular bone marrow has a much higher reversible transverse relaxation rate than surrounding non-bone fatty deposits and other musculoskeletal tissues. This observation, combined with a choice of gradient-echo spacing that accentuates Dixon-type oscillations from chemical-shift interference effects, enabled the isolation of bone marrow signal from surrounding tissues through the use of simple image subtraction and thresholding. Three-dimensional renderings of the marrow surface were then readily generated with this approach - renderings that may prove useful for bone morphology assessment, e.g. for the measurement of femoral anteversion. In conclusion, understanding the behavior of signals from bone marrow and surrounding tissue as a function of time through a spin echo facilitates the segmentation and reconstruction of bone marrow surfaces using straightforward post-processing strategies that are typically available on modern radiology workstations.

Purpose T2 relaxometry based on multiexponential fitting to a single slice multiecho sequence has been the most common MRI technique for myelin water fraction mapping, where the short T2 is associated with myelin water. However, very long acquisition times and physically unrealistic models for T2 distribution are limitations of this approach. We present a novel framework for myelin imaging which substantially increases the imaging speed and myelin water fraction estimation accuracy. Method We used the 2D multislice Carr-Purcell-Meiboom-Gill sequence to increase the volume coverage. To compensate for nonideal slice profiles, we numerically solved the Bloch equations for a range of T2 and B1 inhomogeneity scales to construct the bases for the estimation of the T2 distribution. We used a finite mixture of continuous parametric distributions to describe the complete T2 spectrum and used the constrained variable projection optimization algorithm to estimate myelin water fraction. To validate our model, synthetic, phantom, and in vivo brain experiments were conducted. Results Using the Bloch equations, we can model the slice profile and construct the forward model of the T2 curve. Our method estimated myelin water fraction with smaller error than the nonnegative least squares algorithm. Conclusions The proposed framework can be used for reliable whole brain myelin imaging with a resolution of in . Magn Reson Med, 2015. © 2015 Wiley Periodicals, Inc.

A biomarker of cancer aggressiveness, such as hypoxia, could substantially impact treatment decisions in the prostate, especially radiation therapy, by balancing treatment morbidity (urinary incontinence, erectile dysfunction, etc.) against mortality. R2 (*) mapping with Mono-Exponential (ME) decay modeling has shown potential for identifying areas of prostate cancer hypoxia at 1.5T. However, Gaussian deviations from ME decay have been observed in other tissues at 3T. The purpose of this study is to assess whether gradient-echo signal decays are better characterized by a standard ME decay model, or a Gaussian Augmentation of the Mono-Exponential (GAME) decay model, in the prostate at 3T. Multi-gradient-echo signals were acquired on 20 consecutive patients with a clinical suspicion of prostate cancer undergoing MR-guided prostate biopsies. Data were fitted with both ME and GAME models. The information contents of these models were compared using Akaike's information criterion (second order, AICC ), in skeletal muscle, the prostate central gland (CG), and peripheral zone (PZ) regions of interest (ROIs). The GAME model had higher information content in 30% of the prostate on average (across all patients and ROIs), covering up to 67% of cancerous PZ ROIs, and up to 100% of cancerous CG ROIs (in individual patients). The higher information content of GAME became more prominent in regions that would be assumed hypoxic using ME alone, reaching 50% of the PZ and 70% of the CG as ME R2 (*) approached 40 s(-1) . R2 (*) mapping may have important applications in MRI; however, information lost due to modeling could mask differences in parameters due to underlying tissue anatomy or physiology. The GAME model improves characterization of signal behavior in the prostate at 3T, and may increase the potential for determining correlates of fit parameters with biomarkers, for example of oxygenation status.

T2 relaxometry based on multiexponential fitting to a single slice multiecho sequence has been the most common MRI technique for myelin water fraction mapping, where the short T2 is associated with myelin water. However, very long acquisition times and physically unrealistic models for T2 distribution are limitations of this approach. We present a novel framework for myelin imaging which substantially increases the imaging speed and myelin water fraction estimation accuracy.We used the 2D multislice Carr-Purcell-Meiboom-Gill sequence to increase the volume coverage. To compensate for nonideal slice profiles, we numerically solved the Bloch equations for a range of T2 and B1 inhomogeneity scales to construct the bases for the estimation of the T2 distribution. We used a finite mixture of continuous parametric distributions to describe the complete T2 spectrum and used the constrained variable projection optimization algorithm to estimate myelin water fraction. To validate our model, synthetic, phantom, and in vivo brain experiments were conducted.Using the Bloch equations, we can model the slice profile and construct the forward model of the T2 curve. Our method estimated myelin water fraction with smaller error than the nonnegative least squares algorithm.The proposed framework can be used for reliable whole brain myelin imaging with a resolution of 2×2×4 mm3 in ≈17 min. Magn Reson Med 76:1301-1313, 2016. © 2015 Wiley Periodicals, Inc.

Objective The increasing spatial precision in functional brain mapping is a methodological key issue driving the progress in brain research. Mapping of previously invisible fine scale details permits the discovery of unknown functional brain properties. Currently, the available resolution of fMRI examinations is increasing rapidly. However, data analysis techniques lag behind this development, and achieve an effective resolution which is smaller as the maximum possible spatial precision. We present a data analysis technique with maximum spatial precision in group analyses of high resolution functional magnetic resonance imaging (fMRI) data. Methods The contralateral fingertip representations of thumb, index and middle finger were investigated in a group of 18 participants, for right and left hand. A block design was applied (absolute duration of 200 s) including long counting pulses as a temporal identification task for attention monitoring ( Fig. 1 ). Imaging was carried out with a 3 Tesla MRI-scanner (Verio, Siemens, Germany). FreeSurfer recommendations for the anatomy scans were followed (standard MPRAGE sequence). A standard gradient-echo EPI sequence was used for fMRI (resolution 1.5 × 1.5 × 2 mm 3 , 17 slices, field of view = 191 × 179 mm 2 , TE/TR = 3.6 ms/21 ms). Data analysis was performed with an automated protocol ( Pfannmoller et al., 2015 ), including vessel artefact removal. Cortical positions of the fingertips were computed on individual and MNI brains, using volume-based linear or surface-based non-linear normalization. An assessment of cortical distance measures, for mapping of functional properties, was carried out including straight-line distances, using the 3-D Euclidean measure, and shortest connections within the cortical surface, either with the Dijkstra or LOS-Floyd algorithms. Results In BA3b a distinct somatotopy was found for both hands. The maximum spatial fingertip spread in BA3b after linear normalization was at least two times larger than for the non-linear normalization. All properties of the BA3b fingertip somatotopy were left-right symmetric. Since the representation in BA1 did not exhibit any of those features it was neglected. The LOS-Floyd algorithm achieved an error free result for the distances within measurement accuracy. In contrast, Euclidean distances did not account for the relevant cortical properties. Dijkstra distances were sufficiently precise in the individual brains, but excessively overestimated the left-hand distances after non-linear normalization. Discussion & conclusion Since the representations in BA1 and BA3b are entirely different, the cytoarchitectonic divisions need to be accounted for. The convergence of the BA3b fingertip positions after non-linear normalization implies a removal of the anatomical variability, leaving exclusively functional variability. In the distance mapping, application of the LOS-Floyd algorithm was mandatory in the group template. The Dijkstra algorithm represents a fast alternative for distance mapping in individual brains, while the precision of Euclidean distances did not suffice for our purposes. In conclusion, our methodology improves the spatial precision of somatotopic mapping at 3 T. At larger field strengths improvements to even higher resolutions are anticipated.