Physiological modeling of the BOLD signal and implications for effective connectivity: A primer.

In this primer, I provide an overview of the physiological processes that contribute to the observed BOLD signal (i.e., the generative biophysical model), including their time course properties within the framework of the physiologically-informed dynamic causal modeling (P-DCM). The BOLD signal is primarily determined by the change in paramagnetic deoxygenated hemoglobin, which results from combination of changes in oxygen metabolism, and cerebral blood flow and volume. Specifically, the physiological origin of the so-called BOLD signal "transients" will be discussed, including the initial overshoot, steady-state activation and the post-stimulus undershoot. I argue that incorrect physiological assumptions in the generative model of the BOLD signal can lead to incorrect inferences pertaining to both local neuronal activity and effective connectivity between brain regions. In addition, I introduce the recent laminar BOLD signal model, which extends P-DCM to cortical depths-resolved BOLD signals, allowing for laminar neuronal activity to be determined using high-resolution fMRI data.

Investigations of hypoxia-induced deoxyhemoglobin as a contrast agent for cerebral perfusion imaging.

The assessment of resting perfusion measures (mean transit time, cerebral blood flow, and cerebral blood volume) with magnetic resonance imaging currently requires the presence of a susceptibility contrast agent such as gadolinium. Here, we present an initial comparison between perfusion measures obtained using hypoxia-induced deoxyhemoglobin and gadolinium in healthy study participants. We hypothesize that resting cerebral perfusion measures obtained using precise changes of deoxyhemoglobin concentration will generate images comparable to those obtained using a clinical standard, gadolinium. Eight healthy study participants were recruited (6F; age 23-60). The study was performed using a 3-Tesla scanner with an eight-channel head coil. The experimental protocol consisted of a high-resolution T1-weighted scan followed by two BOLD sequence scans in which each participant underwent a controlled bolus of transient pulmonary hypoxia, and subsequently received an intravenous bolus of gadolinium. The resting perfusion measures calculated using hypoxia-induced deoxyhemoglobin and gadolinium yielded maps that looked spatially comparable. There was no statistical difference between methods in the average voxel-wise measures of mean transit time, relative cerebral blood flow and relative cerebral blood volume, in the gray matter or white matter within each participant. We conclude that perfusion measures generated with hypoxia-induced deoxyhemoglobin are spatially and quantitatively comparable to those generated from a gadolinium injection in the same healthy participant.

Magnetic resonance imaging at 9.4 T: the Maastricht journey.

The 9.4 T scanner in Maastricht is a whole-body magnet with head gradients and parallel RF transmit capability. At the time of the design, it was conceptualized to be one of the best fMRI scanners in the world, but it has also been used for anatomical and diffusion imaging. 9.4 T offers increases in sensitivity and contrast, but the technical ultra-high field (UHF) challenges, such as field inhomogeneities and constraints set by RF power deposition, are exacerbated compared to 7 T. This article reviews some of the 9.4 T work done in Maastricht. Functional imaging experiments included blood oxygenation level-dependent (BOLD) and blood-volume weighted (VASO) fMRI using different readouts. BOLD benefits from shorter T2* at 9.4 T while VASO from longer T1. We show examples of both ex vivo and in vivo anatomical imaging. For many applications, pTx and optimized coils are essential to harness the full potential of 9.4 T. Our experience shows that, while considerable effort was required compared to our 7 T scanner, we could obtain high-quality anatomical and functional data, which illustrates the potential of MR acquisitions at even higher field strengths. The practical challenges of working with a relatively unique system are also discussed.

Effects of MP2RAGE B1+ sensitivity on inter-site T1 reproducibility and hippocampal morphometry at 7T.

Most neuroanatomical studies are based on T1-weighted MR images, whose intensity profiles are not solely determined by the tissue's longitudinal relaxation times (T1), but also affected by varying non-T1 contributions, hampering data reproducibility. In contrast, quantitative imaging using the MP2RAGE sequence, for example, allows direct characterization of the brain based on the tissue property of interest. Combined with 7 Tesla (7T) MRI, this offers unique opportunities to obtain robust high-resolution brain data characterized by a high reproducibility, sensitivity and specificity. However, specific MP2RAGE parameter choices - e.g., to emphasize intracortical myelin-dependent contrast variations - can substantially impact image quality and cortical analyses through remnants of B1+-related intensity variations, as illustrated in our previous work. To follow up on this: we (1) validate this protocol effect using a dataset acquired with a particularly B1+ insensitive set of MP2RAGE parameters combined with parallel transmission excitation; and (2) extend our analyses to evaluate the effects on hippocampal morphometry. The latter remained unexplored initially, but can provide important insights related to generalizability and reproducibility of neurodegenerative research using 7T MRI. We confirm that B1+ inhomogeneities have a considerably variable effect on cortical T1 estimates, as well as on hippocampal morphometry depending on the MP2RAGE setup. While T1 differed substantially across datasets initially, we show the inter-site T1 comparability improves after correcting for the spatially varying B1+ field using a separately acquired Sa2RAGE B1+ map. Finally, removal of B1+ residuals affects hippocampal volumetry and boundary definitions, particularly near structures characterized by strong intensity changes (e.g. cerebral spinal fluid). Taken together, we show that the choice of MP2RAGE parameters can impact T1 comparability across sites and present evidence that hippocampal segmentation results are modulated by B1+ inhomogeneities. This calls for careful (1) consideration of sequence parameters when setting acquisition protocols, as well as (2) acquisition of a B1+ map to correct MP2RAGE data for potential B1+ variations to allow comparison across datasets.

Quantitative and simultaneous measurement of oxygen consumption rates in rat brain and skeletal muscle using 17 O MRS imaging at 16.4T.

PURPOSE: Oxygen-17 (17 O) MRS imaging, successfully used in the brain, is extended by imaging the oxygen metabolic rate in the resting skeletal muscle and used to determine the total whole-body oxygen metabolic rate in the rat. METHODS: During and after inhalations of 17 O2 gas, dynamic 17 O MRSI was performed in rats (n = 8) ventilated with N2 O or N2 at 16.4 T. Time courses of the H217 O concentration from regions of interest located in brain and muscle tissue were examined and used to fit an animal-adapted 3-phase metabolic model of oxygen consumption. CBF was determined with an independent washout method. Finally, body oxygen metabolic rate was calculated using a global steady-state approach. RESULTS: Cerebral metabolic rate of oxygen consumption was 1.97 ± 0.19 µmol/g/min on average. The resting metabolic rate of oxygen consumption in skeletal muscle was 0.32 ± 0.12 µmol/g/min and >6 times lower than cerebral metabolic rate of oxygen consumption. Global oxygen consumed by the body was 24.2 ± 3.6 mL O2 /kg body weight/min. CBF was estimated to be 0.28 ± 0.02 mL/g/min and 0.34 ± 0.06 mL/g/min for the N2 and N2 O ventilation condition, respectively. CONCLUSION: We have evaluated the feasibility of 17 O MRSI for imaging and quantifying the oxygen consumption rate in low metabolizing organs such as the skeletal muscle at rest. Additionally, we have shown that CBF is slightly increased in the case of ventilation with N2 O. We expect this study to be beneficial to the application of 17 O MRSI to a wider range of organs, although further validation is advised.

Perfusion MRI using endogenous deoxyhemoglobin as a contrast agent: Preliminary data.

PURPOSE: To demonstrate the feasibility of mapping cerebral perfusion metrics with BOLD MRI during modulation of pulmonary venous oxygen saturation. METHODS: A gas blender with a sequential gas delivery breathing circuit was used to implement rapid isocapnic changes in the partial pressure of oxygen of the arterial blood. Partial pressure of oxygen was initially lowered to a baseline of 40 mmHg. It was then rapidly raised to 95 mmHg for 20 s before rapidly returning to baseline. The induced cerebral changes in deoxyhemoglobin concentration were tracked over time using BOLD MRI in 6 healthy subjects and 1 patient with cerebral steno-occlusive disease. BOLD signal change, contrast-to-noise ratio, and time delay metrics were calculated. Perfusion metrics such as mean transit time, relative cerebral blood volume, and relative cerebral blood flow were calculated using a parametrized method with a mono-exponential residue function. An arterial input function from within the middle cerebral artery was used to scale relative cerebral blood volume and calculate absolute cerebral blood volume and cerebral blood flow. RESULTS: In normal subjects, average gray and white matter were: BOLD change = 6.3 ± 1.2% and 2.5 ± 0.6%, contrast-to-noise ratio = 4.3 ± 1.3 and 2.6 ± 0.7, time delay = 2.3 ± 0.6 s and 3.6 ± 0.7 s, mean transit time = 3.9 ± 0.6 s and 5.5 ± 0.6 s, relative cerebral blood volume = 3.7 ± 0.9 and 1.6 ± 0.4, relative cerebral blood flow = 70.1 ± 8.3 and 20.6 ± 4.0, cerebral blood flow volume = 4.1 ± 0.9 mL/100 g and 1.8 ± 0.5 mL/100 g, and cerebral blood flow = 97.2 ± 18.7 mL/100 g/min and 28.7 ± 5.9 mL/100 g/min. CONCLUSION: This study demonstrates that induced abrupt changes in deoxyhemoglobin can function as a noninvasive vascular contrast agent that may be used for cerebral perfusion imaging.

A dynamical model of the laminar BOLD response.

High-resolution functional magnetic resonance imaging (fMRI) using blood oxygenation dependent level-dependent (BOLD) signal is an increasingly popular tool to non-invasively examine neuronal processes at the mesoscopic level. However, as the BOLD signal stems from hemodynamic changes, its temporal and spatial properties do not match those of the underlying neuronal activity. In particular, the laminar BOLD response (LBR), commonly measured with gradient-echo (GE) MRI sequence, is confounded by non-local changes in deoxygenated hemoglobin and cerebral blood volume propagated within intracortical ascending veins, leading to a unidirectional blurring of the neuronal activity distribution towards the cortical surface. Here, we present a new cortical depth-dependent model of the BOLD response based on the principle of mass conservation, which takes the effect of ascending (and pial) veins on the cortical BOLD responses explicitly into account. It can be used to dynamically model cortical depth profiles of the BOLD signal as a function of various baseline- and activity-related physiological parameters for any spatiotemporal distribution of neuronal changes. We demonstrate that the commonly observed spatial increase of LBR is mainly due to baseline blood volume increase towards the surface. In contrast, an occasionally observed local maximum in the LBR (i.e. the so-called "bump") is mainly due to spatially inhomogeneous neuronal changes rather than locally higher baseline blood volume. In addition, we show that the GE-BOLD signal laminar point-spread functions, representing the signal leakage towards the surface, depend on several physiological parameters and on the level of neuronal activity. Furthermore, even in the case of simultaneous neuronal changes at each depth, inter-laminar delays of LBR transients are present due to the ascending vein. In summary, the model provides a conceptual framework for the biophysical interpretation of common experimental observations in high-resolution fMRI data. In the future, the model will allow for deconvolution of the spatiotemporal hemodynamic bias of the LBR and provide an estimate of the underlying laminar excitatory and inhibitory neuronal activity.

Direct visualization and characterization of the human zona incerta and surrounding structures.

The zona incerta (ZI) is a small gray matter region of the deep brain first identified in the 19th century, yet direct in vivo visualization and characterization has remained elusive. Noninvasive detection of the ZI and surrounding region could be critical to further our understanding of this widely connected but poorly understood deep brain region and could contribute to the development and optimization of neuromodulatory therapies. We demonstrate that high resolution (submillimetric) longitudinal (T1) relaxometry measurements at high magnetic field strength (7 T) can be used to delineate the ZI from surrounding white matter structures, specifically the fasciculus cerebellothalamicus, fields of Forel (fasciculus lenticularis, fasciculus thalamicus, and field H), and medial lemniscus. Using this approach, we successfully derived in vivo estimates of the size, shape, location, and tissue characteristics of substructures in the ZI region, confirming observations only previously possible through histological evaluation that this region is not just a space between structures but contains distinct morphological entities that should be considered separately. Our findings pave the way for increasingly detailed in vivo study and provide a structural foundation for precise functional and neuromodulatory investigation.

Locus coeruleus imaging as a biomarker for noradrenergic dysfunction in neurodegenerative diseases.

Pathological alterations to the locus coeruleus, the major source of noradrenaline in the brain, are histologically evident in early stages of neurodegenerative diseases. Novel MRI approaches now provide an opportunity to quantify structural features of the locus coeruleus in vivo during disease progression. In combination with neuropathological biomarkers, in vivo locus coeruleus imaging could help to understand the contribution of locus coeruleus neurodegeneration to clinical and pathological manifestations in Alzheimer's disease, atypical neurodegenerative dementias and Parkinson's disease. Moreover, as the functional sensitivity of the noradrenergic system is likely to change with disease progression, in vivo measures of locus coeruleus integrity could provide new pathophysiological insights into cognitive and behavioural symptoms. Locus coeruleus imaging also holds the promise to stratify patients into clinical trials according to noradrenergic dysfunction. In this article, we present a consensus on how non-invasive in vivo assessment of locus coeruleus integrity can be used for clinical research in neurodegenerative diseases. We outline the next steps for in vivo, post-mortem and clinical studies that can lay the groundwork to evaluate the potential of locus coeruleus imaging as a biomarker for neurodegenerative diseases.

Non-BOLD contrast for laminar fMRI in humans: CBF, CBV, and CMRO2.

Functional magnetic resonance imaging (fMRI) using the blood oxygenation level-dependent (BOLD) contrast indirectly probes neuronal activity changes via evoked cerebral blood flow (CBF), cerebral blood volume (CBV) and cerebral metabolic rate of oxygen (CMRO2) changes. The gradient-echo BOLD signal is mostly sensitive to ascending veins in the tissue and to pial veins. Thereby, the achievable spatial specificity to neuronal activation is limited. Furthermore, the non-linear interaction of CBF, CBV and CMRO2 can hamper quantitative interpretations of the BOLD signal across cortical depths with different baseline physiology. Measuring CBF, CBV or CMRO2 directly on a depth-dependent level has the potential to overcome these limitations. Here, we review these candidates of physiologically well-defined contrasts with the particular focus on arterial spin labeling (ASL), vascular space occupancy (VASO) and calibrated fMRI. These methods are reviewed with respect to their fMRI sequence parameter space and the applicability for neuroscientific studies in humans. We show representative results of depth-dependent 'non-BOLD-fMRI' in humans and their spatiotemporal characteristics. We conclude that non-BOLD methods are promising alternatives compared to conventional fMRI as they can provide improved spatial specificity, quantifiability and, hence, physiological interpretability as a function of cortical depth. At submillimeter resolution with inherently low signal-to-noise ratio (SNR), however, their use is still challenging. Nevertheless, we believe that 'non-BOLD-fMRI' is a useful alternative for depth-dependent investigations, by providing valuable insights into neurovascular coupling models that facilitate the interpretability of fMRI for neuroscientific applications.

Linking brain vascular physiology to hemodynamic response in ultra-high field MRI.

Functional MRI using blood oxygenation level-dependent (BOLD) contrast indirectly probes neuronal activity via evoked cerebral blood volume (CBV) and oxygenation changes. Thus, its spatio-temporal characteristics are determined by vascular physiology and MRI parameters. In this paper, we focus on the spatial distribution and time course of the fMRI signal and their magnetic field strength dependence. Even though much is still unknown, the following consistent picture is emerging: a) For high spatial resolution imaging, fMRI contrast-to-noise increases supra-linearly with field strength. b) The location and spacing of penetrating arteries and ascending veins in the cortical tissue are not correlated to cortical columns, imposing limitations on achievable point-spread function (PSF) in fMRI. c) Baseline CBV distribution may vary over cortical layers biasing fMRI signal to layers with high CBV values. d) The largest CBV change is in the tissue microvasculature, less in surface arteries and even less in pial veins. e) Venous CBV changes are only relevant for longer stimuli, and oxygenation changes are largest in post-capillary blood vessels. f) The balloon effect (i.e. slow recovery of CBV to baseline) is located in the tissue, consistent with the fact that the post-stimulus undershoot has narrower spatial PSF than the positive BOLD response. g) The onset time following stimulation has been found to be shortest in middle/lower layers, both in optical imaging and high-resolution fMRI, but we argue and demonstrate with simulations that varying signal latencies can also be caused by vascular properties and, therefore, may potentially not be interpreted as neural latencies. With simulations, we illustrate the field strength dependency of fMRI signal transients, such as the adaptation during stimulation, initial dip and the post-stimulus undershoot. In sum, vascular structure and function impose limitations on the achievable PSF of fMRI and give rise to complex fMRI transients, which contain time-varying amount of excitatory and inhibitory neuronal information. Nevertheless, non-invasive fMRI at ultra-high magnetic fields not only provides high contrast-to-noise but also an unprecedented detailed view on cognitive processes in the human brain.

Impact of acquisition and analysis strategies on cortical depth-dependent fMRI.

Functional MRI at ultra-high magnetic fields (≥ 7T) provides the opportunity to probe columnar and laminar processing in the human brain in vivo at sub-millimeter spatial scales. However, fMRI data only indirectly reflects the neuronal laminar profile due to a bias to ascending and pial veins inherent in gradient- and spin-echo BOLD fMRI. In addition, accurate delineation of the cortical depths is difficult, due to the relatively large voxel sizes and lack of sufficient tissue contrast in the functional images. In conventional depth-dependent fMRI studies, anatomical and functional data are acquired with different image read-out modules, the fMRI data are distortion-corrected and vascular biases are accounted for by subtracting the depth-dependent activation profiles of different stimulus conditions. In this study, using high-resolution gradient-echo fMRI data (0.7 mm isotropic) of the human visual cortex, we propose instead, that depth-dependent functional information is best preserved if data analysis is performed in the original functional data space. To achieve this, we acquired anatomical images with high tissue contrast and similar distortion to the functional images using multiple inversion-recovery time EPI, thereby eliminating the need to un-distort the fMRI data. We demonstrate higher spatial accuracy for the cortical layer definitions of this approach as compared to the more conventional approach using MP2RAGE anatomy. In addition, we provide theoretical arguments and empirical evidence that vascular biases can be better accounted for using division instead of subtraction of the depth-dependent profiles. Finally, we show that the hemodynamic response of grey matter has relatively stronger post-stimulus undershoot than the pial vein voxels. In summary, we show that the choice of fMRI data acquisition and processing can impact observable differences in the cortical depth profiles and present evidence that cortical depth-dependent modulation of the BOLD signal can be resolved using gradient-echo imaging.

Techniques for blood volume fMRI with VASO: From low-resolution mapping towards sub-millimeter layer-dependent applications.

Quantitative cerebral blood volume (CBV) fMRI has the potential to overcome several specific limitations of BOLD fMRI. It provides direct physiological interpretability and promises superior localization specificity in applications of sub-millimeter resolution fMRI applications at ultra-high magnetic fields (7T and higher). Non-invasive CBV fMRI using VASO (vascular space occupancy), however, is inherently limited with respect to its data acquisition efficiency, restricting its imaging coverage and achievable spatial and temporal resolution. This limitation may be reduced with recent advanced acceleration and reconstruction strategies that allow two-dimensional acceleration, such as in simultaneous multi-slice (SMS) 2D-EPI or 3D-EPI in combination with CAIPIRINHA field-of-view shifting. In this study, we sought to determine the functional sensitivity and specificity of these readout strategies with VASO over a broad range of spatial resolutions; spanning from low spatial resolution (3mm) whole-cortex to sub-millimeter (0.75mm) slab-of-cortex (for cortical layer-dependent applications). In the thermal-noise-dominated regime of sub-millimeter resolutions, 3D-EPI-VASO provides higher temporal stability and sensitivity to detect changes in CBV compared to 2D-EPI-VASO. In this regime, 3D-EPI-VASO unveils task activation located in the cortical laminae with little contamination from surface veins, in contrast to the cortical surface weighting of GE-BOLD fMRI. In the physiological-noise-dominated regime of lower resolutions, however, 2D-SMS-VASO shows superior performance compared to 3D-EPI-VASO. Due to its superior sensitivity at a layer-dependent level, 3D-EPI VASO promises to play an important role in future neuroscientific applications of layer-dependent fMRI.

The impact of ultra-high field MRI on cognitive and computational neuroimaging.

The ability to measure functional brain responses non-invasively with ultra high field MRI (7 T and above) represents a unique opportunity in advancing our understanding of the human brain. Compared to lower fields (3 T and below), ultra high field MRI has an increased sensitivity, which can be used to acquire functional images with greater spatial resolution, and greater specificity of the blood oxygen level dependent (BOLD) signal to the underlying neuronal responses. Together, increased resolution and specificity enable investigating brain functions at a submillimeter scale, which so far could only be done with invasive techniques. At this mesoscopic spatial scale, perception, cognition and behavior can be probed at the level of fundamental units of neural computations, such as cortical columns, cortical layers, and subcortical nuclei. This represents a unique and distinctive advantage that differentiates ultra high from lower field imaging and that can foster a tighter link between fMRI and computational modeling of neural networks. So far, functional brain mapping at submillimeter scale has focused on the processing of sensory information and on well-known systems for which extensive information is available from invasive recordings in animals. It remains an open challenge to extend this methodology to uniquely human functions and, more generally, to systems for which animal models may be problematic. To succeed, the possibility to acquire high-resolution functional data with large spatial coverage, the availability of computational models of neural processing as well as accurate biophysical modeling of neurovascular coupling at mesoscopic scale all appear necessary.

High-resolution in vivo imaging of human locus coeruleus by magnetization transfer MRI at 3T and 7T.

Locus Coeruleus (LC) is a neuromelanin-rich brainstem structure that is the source of noradrenaline in the cortex and is thought to modulate attention and memory. LC imaging in vivo is commonly performed with a 2D T1-weighted Turbo Spin Echo (TSE) MRI sequence, an approach that suffers from several drawbacks at 3T, including long acquisition times and highly anisotropic spatial resolution. In this study, we developed a high-resolution Magnetization Transfer (MT) sequence for LC imaging at both 7T and 3T and compared its performance to a TSE sequence. Results indicate that LC imaging can be achieved with an MT sequence at both 7 and 3T at higher spatial resolution than the 3T TSE. Furthermore, we investigated whether the currently disputed source of contrast in the LC region with a TSE sequence relates to MT effects or shortened T1 and T2* due to increased iron concentration. Our results suggest that the contrast in the LC area relates to MT effects. To conclude, in this study we managed to image the LC, for the first time, at 7T and at an increased resolution compared to the current state-of-the-art. Imaging the LC is highly relevant for clinical diagnostics as structural tissue properties of the LC may hold promise as a biomarker in neurodegenerative diseases.

Distortion-matched T1 maps and unbiased T1-weighted images as anatomical reference for high-resolution fMRI.

The increasing availability of ultra-high field scanners has led to a growing number of submillimetre fMRI studies in humans, typically targeting the gray matter at different cortical depths. In most analyses, the definition of surfaces at different cortical depths is based on an anatomical image with different contrast and distortions than the functional images. Here, we introduce a novel sequence providing bias-field corrected T1-weighted images and T1-maps with distortions that match those of the fMRI data, with an image acquisition time significantly shorter than standard T1-weighted anatomical imaging. For 'T1-imaging with 2 3D-EPIs', or T123DEPI, 3D-EPI volumes are acquired centred at two inversion times. These 3D-EPIs are segmented into half, quarter or smaller blocks of k-space to allow for optimisation of the inversion times. T1-weighted images and T1-maps are then generated as for MP2RAGE acquisitions. A range of T123DEPI data acquired at 7 T is shown with resolutions ranging from 0.7 mm to 1.3 mm isotropic voxels. Co-registration quality to the mean EPI of matching fMRI timecourses shows markedly less local deviations compared to co-registration of a standard MP2RAGE to the same echo planar volume. Thus, the T123DEPI T1-weighted images and T1-maps can be used to provide cortical surfaces with matched distortions to the functional data or else to facilitate co-registration between functional and undistorted anatomical data.

Individualized parcellation of the subthalamic nucleus in patients with Parkinson's disease with 7T MRI.

Deep brain stimulation of the subthalamic nucleus (STN) is a widely performed surgical treatment for patients with Parkinson's disease. The goal of the surgery is to place an electrode centered in the motor region of the STN while lowering the effects of electrical stimulation on the non-motor regions. However, distinguishing the motor region from the neighboring associative and limbic areas in individual patients using imaging modalities was until recently difficult to obtain in vivo. Here, using ultra-high field MR imaging, we have performed a dissection of the subdivisions of the STN of individual Parkinson's disease patients. We have acquired 7T diffusion-weighted images of seventeen patients with Parkinson's disease scheduled for deep brain stimulation surgery. Using a structural connectivity-based parcellation protocol, the STN's connections to the motor, limbic, and associative cortical areas were used to map the individual subdivisions of the nucleus. A reproducible patient-specific parcellation of the STN into a posterolateral motor and gradually overlapping central associative area was found in all STNs, taking up on average 55.3% and 55.6% of the total nucleus volume. The limbic area was found in the anteromedial part of the nucleus. Our results suggest that 7T MR imaging may facilitate individualized and highly specific planning of deep brain stimulation surgery of the STN.

Ultra-high resolution blood volume fMRI and BOLD fMRI in humans at 9.4 T: Capabilities and challenges.

Functional mapping of cerebral blood volume (CBV) changes has the potential to reveal brain activity with high localization specificity at the level of cortical layers and columns. Non-invasive CBV imaging using Vascular Space Occupancy (VASO) at ultra-high magnetic field strengths promises high spatial specificity but poses unique challenges in human applications. As such, 9.4 T B1+ and B0 inhomogeneities limit efficient blood tagging, while the specific absorption rate (SAR) constraints limit the application of VASO-specific RF pulses. Moreover, short T2* values at 9.4 T require short readout duration, and long T1 values at 9.4 T can cause blood-inflow contaminations. In this study, we investigated the applicability of layer-dependent CBV-fMRI at 9.4 T in humans. We addressed the aforementioned challenges by combining multiple technical advancements: temporally alternating pTx B1+ shimming parameters, advanced adiabatic RF-pulses, 3D-EPI signal readout, optimized GRAPPA acquisition and reconstruction, and stability-optimized RF channel combination. We found that a combination of suitable advanced methodology alleviates the challenges and potential artifacts, and that VASO fMRI provides reliable measures of CBV change across cortical layers in humans at 9.4 T. The localization specificity of CBV-fMRI, combined with the high sensitivity of 9.4 T, makes this method an important tool for future studies investigating cortical micro-circuitry in humans.

The impact of B1+ correction on MP2RAGE cortical T1 and apparent cortical thickness at 7T.

Determination of cortical thickness using MRI has often been criticized due to the presence of various error sources. Specifically, anatomical MRI relying on T1 contrast may be unreliable due to spatially variable image contrast between gray matter (GM), white matter (WM) and cerebrospinal fluid (CSF). Especially at ultra-high field (≥ 7T) MRI, transmit and receive B1 -related image inhomogeneities can hamper correct classification of tissue types. In the current paper, we demonstrate that residual B1+ (transmit) inhomogeneities in the T1 -weighted and quantitative T1 images using the MP2RAGE sequence at 7T lead to biases in cortical thickness measurements. As expected, post-hoc correction for the spatially varying B1+ profile reduced the apparent T1 values across the cortex in regions with low B1+, and slightly increased apparent T1 in regions with high B1+. As a result, improved contrast-to-noise ratio both at the GM-CSF and GM-WM boundaries can be observed leading to more accurate surface reconstructions and cortical thickness estimates. Overall, the changes in cortical thickness ranged between a 5% decrease to a 70% increase after B1+ correction, reducing the variance of cortical thickness values across the brain dramatically and increasing the comparability with normative data. More specifically, the cortical thickness estimates increased in regions characterized by a strong decrease of apparent T1 after B1+ correction in regions with low B1+ due to improved detection of the pial surface. The current results suggest that cortical thickness can be more accurately determined using MP2RAGE data at 7T if B1+ inhomogeneities are accounted for.

Cortical depth profiles of luminance contrast responses in human V1 and V2 using 7 T fMRI.

Neural activity in early visual cortex is modulated by luminance contrast. Cortical depth (i.e., laminar) contrast responses have been studied in monkey early visual cortex, but not in humans. In addition to the high spatial resolution needed and the ensuing low signal-to-noise ratio, laminar studies in humans using fMRI are hampered by the strong venous vascular weighting of the fMRI signal. In this study, we measured luminance contrast responses in human V1 and V2 with high-resolution fMRI at 7 T. To account for the effect of intracortical ascending veins, we applied a novel spatial deconvolution model to the fMRI depth profiles. Before spatial deconvolution, the contrast response in V1 showed a slight local maximum at mid cortical depth, whereas V2 exhibited a monotonic signal increase toward the cortical surface. After applying the deconvolution, both V1 and V2 showed a pronounced local maximum at mid cortical depth, with an additional peak in deep grey matter, especially in V1. Moreover, we found a difference in contrast sensitivity between V1 and V2, but no evidence for variations in contrast sensitivity as a function of cortical depth. These findings are in agreement with results obtained in nonhuman primates, but further research will be needed to validate the spatial deconvolution approach.