Microvascular Specificity of Spin Echo BOLD fMRI: Impact of EPI Echo Train Length.

A spatially specific fMRI acquisition requires specificity to the microvasculature that serves active neuronal sites. Macrovascular contributions will reduce the microvascular specificity but can be reduced by using spin echo (SE) sequences that use a π pulse to refocus static field inhomogeneities near large veins. The microvascular specificity of a SE-echo planar imaging (SE-EPI) scan depends on the echo train length (ETL)-duration, but the dependence is not well-characterized in humans at 7T. To determine how microvascular-specific SE-EPI BOLD is in humans at 7T, we developed a Monte Carlo voxel model that computes the signal of a proton ensemble residing in a vasculature subjected to a SE-EPI pulse sequence. We characterized the ETL-duration dependence of the microvascular specificity by simulating the BOLD signal as a function of ETL, the range adhering to experimentally realistic readouts. We performed a validation experiment for our simulation observations, in which we acquired a set of SE-EPI BOLD time series with varying ETL during a hyperoxic gas challenge. Both our simulations and measurements show an increase in macrovascular contamination as a function of ETL, with an increase of 30% according to our simulation and 60% according to our validation experiment between the shortest and longest ETL durations (23.1 - 49.7 ms). We conclude that the microvascular specificity decreases heavily with increasing ETL-durations. We recommend reducing the ETL-duration as much as possible to minimize macrovascular contamination in SE-EPI BOLD experiments. We additionally recommend scanning at high resolutions to minimize partial volume effects with CSF. CSF voxels show a large BOLD response, which can be attributed to both the presence of large veins (high blood volume) and molecular oxygen-induced T1-shortening (significant in a hyperoxia experiment). The magnified BOLD signal in a GM-CSF partial volume voxel reduces the desired microvascular specificity and, therefore, will hinder the interpretation of functional MRI activation patterns.

A touch of hierarchy: population receptive fields reveal fingertip integration in Brodmann areas in human primary somatosensory cortex.

Several neuroimaging studies have shown the somatotopy of body part representations in primary somatosensory cortex (S1), but the functional hierarchy of distinct subregions in human S1 has not been adequately addressed. The current study investigates the functional hierarchy of cyto-architectonically distinct regions, Brodmann areas BA3, BA1, and BA2, in human S1. During functional MRI experiments, we presented participants with vibrotactile stimulation of the fingertips at three different vibration frequencies. Using population Receptive Field (pRF) modeling of the fMRI BOLD activity, we identified the hand region in S1 and the somatotopy of the fingertips. For each voxel, the pRF center indicates the finger that most effectively drives the BOLD signal, and the pRF size measures the spatial somatic pooling of fingertips. We find a systematic relationship of pRF sizes from lower-order areas to higher-order areas. Specifically, we found that pRF sizes are smallest in BA3, increase slightly towards BA1, and are largest in BA2, paralleling the increase in visual receptive field size as one ascends the visual hierarchy. Additionally, we find that the time-to-peak of the hemodynamic response in BA3 is roughly 0.5 s earlier compared to BA1 and BA2, further supporting the notion of a functional hierarchy of subregions in S1. These results were obtained during stimulation of different mechanoreceptors, suggesting that different afferent fibers leading up to S1 feed into the same cortical hierarchy.

Topographic representation of numerosity in the human parietal cortex.

Numerosity, the set size of a group of items, is processed by the association cortex, but certain aspects mirror the properties of primary senses. Sensory cortices contain topographic maps reflecting the structure of sensory organs. Are the cortical representation and processing of numerosity organized topographically, even though no sensory organ has a numerical structure? Using high-field functional magnetic resonance imaging (at a field strength of 7 teslas), we described neural populations tuned to small numerosities in the human parietal cortex. They are organized topographically, forming a numerosity map that is robust to changes in low-level stimulus features. The cortical surface area devoted to specific numerosities decreases with increasing numerosity, and the tuning width increases with preferred numerosity. These organizational properties extend topographic principles to the representation of higher-order abstract features in the association cortex.

Frequency specific spatial interactions in human electrocorticography: V1 alpha oscillations reflect surround suppression.

Electrical brain signals are often decomposed into frequency ranges that are implicated in different functions. Using subdural electrocorticography (ECoG, intracranial EEG) and functional magnetic resonance imaging (fMRI), we measured frequency spectra and BOLD responses in primary visual cortex (V1) and intraparietal sulcus (IPS). In V1 and IPS, 30-120 Hz (gamma, broadband) oscillations allowed population receptive field (pRF) reconstruction comparable to fMRI estimates. Lower frequencies, however, responded very differently in V1 and IPS. In V1, broadband activity extends down to 3 Hz. In the 4-7 Hz (theta) and 18-30 Hz (beta) ranges broadband activity increases power during stimulation within the pRF. However, V1 9-12 Hz (alpha) frequency oscillations showed a different time course. The broadband power here is exceeded by a frequency-specific power increase during stimulation of the area outside the pRF. As such, V1 alpha oscillations reflected surround suppression of the pRF, much like negative fMRI responses. They were consequently highly localized, depending on stimulus and pRF position, and independent between nearby electrodes. In IPS, all 3-25 Hz oscillations were strongest during baseline recording and correlated between nearby electrodes, consistent with large-scale disengagement. These findings demonstrate V1 alpha oscillations result from locally active functional processes and relate these alpha oscillations to negative fMRI signals. They highlight that similar oscillations in different areas reflect processes with different functional roles. However, both of these roles of alpha seem to reflect suppression of spiking activity.

Pushing the limits of high-resolution functional MRI using a simple high-density multi-element coil design.

Recent studies have shown that functional MRI (fMRI) can be sensitive to the laminar and columnar organization of the cortex based on differences in the spatial and temporal characteristics of the blood oxygenation level-dependent (BOLD) signal originating from the macrovasculature and the neuronal-specific microvasculature. Human fMRI studies at this scale of the cortical architecture, however, are very rare because the high spatial/temporal resolution required to explore these properties of the BOLD signal are limited by the signal-to-noise ratio. Here, we show that it is possible to detect BOLD signal changes at an isotropic spatial resolution as high as 0.55 mm at 7 T using a high-density multi-element surface coil with minimal electronics, which allows close proximity to the head. The coil comprises of very small, 1 × 2-cm(2) , elements arranged in four flexible modules of four elements each (16-channel) that can be positioned within 1 mm from the head. As a result of this proximity, tissue losses were five-fold greater than coil losses and sufficient to exclude preamplifier decoupling. When compared with a standard 16-channel head coil, the BOLD sensitivity was approximately 2.2-fold higher for a high spatial/temporal resolution (1 mm isotropic/0.4 s), multi-slice, echo planar acquisition, and approximately three- and six-fold higher for three-dimensional echo planar images acquired with isotropic resolutions of 0.7 and 0.55 mm, respectively. Improvements in parallel imaging performance (geometry factor) were up to around 1.5-fold with increasing acceleration factor, and improvements in fMRI detectability (temporal signal-to-noise ratio) were up to around four-fold depending on the distance to the coil. Although deeper lying structures may not benefit from the design, most fMRI questions pertain to the neocortex which lies within approximately 4 cm from the surface. These results suggest that the resolution of fMRI (at 7 T) can approximate levels that are closer to the spatial/temporal scale of the fundamental functional organization of the human cortex using a simple high-density coil design for high sensitivity.

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.

Mapping the MRI voxel volume in which thermal noise matches physiological noise--implications for fMRI.

This work addresses the choice of the imaging voxel volume in blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI). Noise of physiological origin that is present in the voxel time course is a prohibitive factor in the detection of small activation-induced BOLD signal changes. If the physiological noise contribution dominates over the temporal fluctuation contribution in the imaging voxel, further increases in the voxel signal-to-noise ratio (SNR) will have diminished corresponding increases in temporal signal-to-noise (TSNR), resulting in reduced corresponding increases in the ability to detect activation induced signal changes. On the other hand, if the thermal and system noise dominate (suggesting a relatively low SNR) further decreases in SNR can prohibit detection of activation-induced signal changes. Here we have proposed and called the "suggested" voxel volume for fMRI the volume where thermal plus system-related and physiological noise variances are equal. Based on this condition we have created maps of fMRI suggested voxel volume from our experimental data at 3T, since this value will spatially vary depending on the contribution of physiologic noise in each voxel. Based on our fast EPI segmentation technique we have found that for gray matter (GM), white matter (WM), and cerebral spinal fluid (CSF) brain compartments the mean suggested cubical voxel volume is: (1.8 mm)3, (2.1 mm)3 and (1.4 mm)3, respectively. Serendipitously, (1.8 mm)3 cubical voxel volume for GM approximately matches the cortical thickness, thus optimizing BOLD contrast by minimizing partial volume averaging. The introduced suggested fMRI voxel volume can be a useful parameter for choice of imaging volume for functional studies.

Multiparametric MR assessment of pediatric brain tumors.

MR assessment of pediatric brain tumors has expanded to include physiologic information related to cellular metabolites, hemodynamic and diffusion parameters. The purpose of this study was to investigate the relationship between MR and proton MR spectroscopic imaging in children with primary brain tumors. Twenty-one patients (mean age 9 years) with histologically verified brain tumors underwent conventional MR imaging, hemodynamic MR imaging (HMRI) and proton MR spectroscopic imaging (MRSI). Fourteen patients also had diffusion-weighted MR imaging (DWMRI). Metabolic indices including choline-containing compounds (Cho), total creatine (tCr) and lipids/lactate (L) were derived by proton MRSI, relative cerebral blood volume (rCBV) by HMRI, and apparent tissue water diffusion coefficients (ADC) by DWMRI. Variables were examined by linear regression and correlation as well as by ANOVA. Cho (suggestive of tumor cellularity and proliferative activity) correlated positively with rCBV, while the relationship between Cho and ADC (suggestive of cellular density) was inverse ( P<0.001). The relationship between rCBV and ADC was also inverse ( P=0.004). Cho and lipids (suggestive of necrosis and/or apoptosis) were not significantly correlated ( P=0.51). A positive relationship was found between lipids and ADC ( P=0.002). The relationships between Cho, rCBV, ADC and lipids signify that tumor physiology is influenced by the tumor's physical and chemical environment. Normalized Cho and lipids distinguished high-grade from low-grade tumors ( P<0.05). Multiparametric MR imaging using MRSI, HMRI and DWMRI enhances assessment of brain tumors in children and improves our understanding of tumor physiology while promising to distinguish higher- from lower-malignancy tumors, a distinction that is particularly clinically important among inoperable tumors.