An integrative model for neuronal activity-induced signal changes for gradient and spin echo functional imaging

Layer-specific BOLD effects in gradient and spin-echo acquisitions in somatosensory cortex

PURPOSE

Previous studies have shown varied BOLD signals with gradient echo (GE) across cortical depth. To interpret these variations, and understand the effects of vascular geometry and size, the magnitudes and layer distributions of GE and spin-echo (SE) BOLD functional MRI signals were compared in the somatosensory cortex of squirrel monkeys during tactile stimulation and in a resting state at high spatial resolution and high field.

METHODS

A block-design stimulation was used to identify tactile-evoked activation signals in somatosensory Areas 3b and 1. Layer-specific connectivities were calculated using resting-state data. Signal power spectra were compared by depth and pulse sequence. The measured ratios of transverse relaxation rate changes were compared with Anderson and Weiss's model.

RESULTS

SE signals showed a 26% lower percentage signal change during tactile stimulation compared with GE, along with a slower time course. SE signals remained consistent but weaker in lower layers, whereas GE signals decreased with cortical depth. This pattern extended to resting-state power spectra. Resting-state functional connectivity indicated larger connectivity between the top layers of Area 3b and Area 1 for GE, with minimal changes for SE. Comparisons with theory suggest vessel diameters ranging from 19.4 to 9 microns are responsible for BOLD effects across cortical layers at 9.4 T.

CONCLUSION

These results provide further evidence that at high field, SE BOLD signals are relatively free of contributions from sources other than microvascular changes in response to neural activity, whereas GE signals, even in the superficial layers, are not dominated by very large vessels.

Comparison of compartmental analytical Blood-Oxygen-Level-Dependent functional Magnetic Resonance Imaging models against Monte Carlo simulations performed over cortical micro-angiograms

Blood oxygen level-dependent functional magnetic resonance imaging (BOLD fMRI) arises from a physiological and physical cascade of events taking place at the level of the cortical microvasculature which constitutes a medium with complex geometry. Several analytical models of the BOLD contrast have been developed, but these have not been compared directly against detailed bottom-up modeling methods. Using a 3D modeling method based on experimentally measured images of mice microvasculature and Monte Carlo simulations, we quantified the accuracy of two analytical models to predict the amplitude of the BOLD response from 1.5 to 7 T, for different echo time (TE) and for both gradient echo and spin echo acquisition protocols. We also showed that accounting for the tridimensional structure of the microvasculature results in more accurate prediction of the BOLD amplitude, even if the values for SO2 were averaged across individual vascular compartments. A secondary finding is that modeling the venous compartment as two individual compartments results in more accurate prediction of the BOLD amplitude compared with standard homogenous venous modeling, arising from the bimodal distribution of venous SO2 across the microvasculature in our data.

Performance of receive head arrays versus ultimate intrinsic SNR at 7 T and 10.5 T

PURPOSE

We examined magnetic field dependent SNR gains and ability to capture them with multichannel receive arrays for human head imaging in going from 7 T, the most commonly used ultrahigh magnetic field (UHF) platform at the present, to 10.5 T, which represents the emerging new frontier of >10 T in UHFs.

METHODS

Electromagnetic (EM) models of 31-channel and 63-channel multichannel arrays built for 10.5 T were developed for 10.5 T and 7 T simulations. A 7 T version of the 63-channel array with an identical coil layout was also built. Array performance was evaluated in the EM model using a phantom mimicking the size and electrical properties of the human head and a digital human head model. Experimental data was obtained at 7 T and 10.5 T with the 63-channel array. Ultimate intrinsic SNR (uiSNR) was calculated for the two field strengths using a voxelized cloud of dipoles enclosing the phantom or the digital human head model as a reference to assess the performance of the two arrays and field depended SNR gains.

RESULTS

uiSNR calculations in both the phantom and the digital human head model demonstrated SNR gains at 10.5 T relative to 7 T of 2.6 centrally, ˜2 at the location corresponding to the edge of the brain, ˜1.4 at the periphery. The EM models demonstrated that, centrally, both arrays captured ˜90% of the uiSNR at 7 T, but only ˜65% at 10.5 T, leading only to ˜2-fold gain in array SNR in going from 7 to 10.5 T. This trend was also observed experimentally with the 63-channel array capturing a larger fraction of the uiSNR at 7 T compared to 10.5 T, although the percentage of uiSNR captured were slightly lower at both field strengths compared to EM simulation results.

CONCLUSIONS

Major uiSNR gains are predicted for human head imaging in going from 7 T to 10.5 T, ranging from ˜2-fold at locations corresponding to the edge of the brain to 2.6-fold at the center, corresponding to approximately quadratic increase with the magnetic field. Realistic 31- and 63-channel receive arrays, however, approach the central uiSNR at 7 T, but fail to do so at 10.5 T, suggesting that more coils and/or different type of coils will be needed at 10.5 T and higher magnetic fields.

Imaging arterial and venous vessels using Iron Dextran enhanced multi-echo 3D gradient echo MRI at 7T

Iron Dextran is a widely used iron oxide compound to treat iron-deficiency anemia patients in the clinic. Similar to other iron oxide compounds such as Ferumoxytol, it can also be used off-label as an intravascular magnetic resonance imaging (MRI) contrast agent due to its strong iron-induced T2 and T2* shortening effects. In this study, we seek to evaluate the feasibility of using Iron Dextran enhanced multi-echo susceptibility weighted imaging (SWI) MRI at 7T to image arterial and venous blood vessels in the human brain. Phantom experiments were performed to measure the r2* relaxivity for Iron Dextran in blood, based on which the SWI sequence was optimized. Pre- and post-infusion MR images were acquired in human subjects from which maps of arteries and veins were extracted. The post-contrast SWI images showed enhanced susceptibility difference between blood and the surrounding tissue in both arteries and veins. Our results showed that the proposed Iron Dextran enhanced multi-echo SWI approach allowed the visualization of blood vessels with diameters down to ~100 μm, including small blood vessels supplying and draining small brain structures such as the hippocampus. We conclude that Iron Dextran can be an alternative iron-based MRI contrast agent for blood vessel imaging in the human brain.

Recurrent inhibition refines mental templates to optimize perceptual decisions

Translating sensory inputs to perceptual decisions relies on building internal representations of features critical for solving complex tasks. Yet, we still lack a mechanistic account of how the brain forms these mental templates of task-relevant features to optimize decision-making. Here, we provide evidence for recurrent inhibition: an experience-dependent plasticity mechanism that refines mental templates by enhancing γ-aminobutyric acid (GABA)-mediated (GABAergic) inhibition and recurrent processing in superficial visual cortex layers. We combine ultrahigh-field (7 T) functional magnetic resonance imaging at submillimeter resolution with magnetic resonance spectroscopy to investigate the fine-scale functional and neurochemical plasticity mechanisms for optimized perceptual decisions. We demonstrate that GABAergic inhibition increases following training on a visual (i.e., fine orientation) discrimination task, enhancing the discriminability of orientation representations in superficial visual cortex layers that are known to support recurrent processing. Modeling functional and neurochemical plasticity interactions reveals that recurrent inhibitory processing optimizes brain computations for perpetual decisions and adaptive behavior.

Laminar multi-contrast fMRI at 7T allows differentiation of neuronal excitation and inhibition underlying positive and negative BOLD responses

A major challenge for human neuroimaging using functional MRI is the differentiation of neuronal excitation and inhibition which may induce positive and negative BOLD responses. Here we present an innovative multi-contrast laminar functional MRI technique that offers comprehensive and quantitative imaging of neurovascular (CBF, CBV, BOLD) and metabolic (CMRO2) responses across cortical layers at 7 Tesla. This technique was first validated through a finger-tapping experiment, revealing 'double-peak' laminar activation patterns within the primary motor cortex. By employing a ring-shaped visual stimulus that elicited positive and negative BOLD responses, we further observed distinct neurovascular and metabolic responses across cortical layers and eccentricities in the primary visual cortex. This suggests potential feedback inhibition of neuronal activities in both superficial and deep cortical layers underlying the negative BOLD signals in the fovea, and also illustrates the neuronal activities in visual areas adjacent to the activated eccentricities.

7 T Lesion-Attenuated Magnetization-Prepared Gradient Echo Acquisition for Detection of Posterior Fossa Demyelinating Lesions in Multiple Sclerosis

OBJECTIVES

Detection of infratentorial demyelinating lesions in multiple sclerosis (MS) presents a challenge in magnetic resonance imaging (MRI), a difficulty that is further heightened in 7 T MRI. This study aimed to assess the efficacy of a novel MRI approach, lesion-attenuated magnetization-prepared gradient echo acquisition (LAMA), for detecting demyelinating lesions within the posterior fossa and upper cervical spine on 7 T MRI and contrast its performance with conventional double-inversion recovery (DIR) and T2-weighted turbo spin echo sequences.

MATERIALS AND METHODS

We conducted a retrospective cross-sectional study in 42 patients with a confirmed diagnosis of MS. All patients had 7 T MRI that incorporated LAMA, 3D DIR, and 2D T2-weighted turbo spin echo sequences. Three readers assessed lesion count in the brainstem, cerebellum, and upper cervical spinal cord using both DIR and T2-weighted images in one session. In a separate session, LAMA was analyzed alone. Contrast-to-noise ratio was also compared between LAMA and the conventional sequences. Lesion counts between methods were assessed using nonparametric Wilcoxon signed rank test. Interrater agreement in lesion detection was estimated by intraclass correlation coefficients.

RESULTS

LAMA identified a significantly greater number of lesions than DIR + T2 (mean 6.4 vs 3.0; P < 0.001). LAMA also exhibited better interrater agreement (intraclass correlation coefficient [95% confidence interval], 0.75 [0.41-0.88] vs 0.61 [0.35-0.78]). The contrast-to-noise ratio for LAMA (3.7 ± 0.9) significantly exceeded that of DIR (1.94 ± 0.7) and T2 (1.2 ± 0.7) (all P 's < 0.001). In cases with no lesions detected using DIR + T2, at least 1 lesion was identified in 83.3% with LAMA. Across all analyzed brain regions, LAMA consistently detected more lesions than DIR + T2.

CONCLUSIONS

LAMA significantly improves the detection of infratentorial demyelinating lesions in MS patients compared with traditional methods. Integrating LAMA with standard magnetization-prepared 2 rapid acquisition gradient echo acquisition provides a valuable tool for accurately characterizing the extent of MS disease.

High Spatiotemporal Resolution Radial Encoding Single-Vessel fMRI

High-field preclinical functional MRI (fMRI) is enabled the high spatial resolution mapping of vessel-specific hemodynamic responses, that is single-vessel fMRI. In contrast to investigating the neuronal sources of the fMRI signal, single-vessel fMRI focuses on elucidating its vascular origin, which can be readily implemented to identify vascular changes relevant to vascular dementia or cognitive impairment. However, the limited spatial and temporal resolution of fMRI is hindered hemodynamic mapping of intracortical microvessels. Here, the radial encoding MRI scheme is implemented to measure BOLD signals of individual vessels penetrating the rat somatosensory cortex. Radial encoding MRI is employed to map cortical activation with a focal field of view (FOV), allowing vessel-specific functional mapping with 50 × 50 µm2 in-plane resolution at a 1 to 2 Hz sampling rate. Besides detecting refined hemodynamic responses of intracortical micro-venules, the radial encoding-based single-vessel fMRI enables the distinction of fMRI signals from vessel and peri-vessel voxels due to the different contribution of intravascular and extravascular effects.

Reliability of task-based fMRI in the dorsal horn of the human spinal cord

The application of functional magnetic resonance imaging (fMRI) to the human spinal cord is still a relatively small field of research and faces many challenges. Here we aimed to probe the limitations of task-based spinal fMRI at 3T by investigating the reliability of spinal cord blood oxygen level dependent (BOLD) responses to repeated nociceptive stimulation across two consecutive days in 40 healthy volunteers. We assessed the test-retest reliability of subjective ratings, autonomic responses, and spinal cord BOLD responses to short heat pain stimuli (1s duration) using the intraclass correlation coefficient (ICC). At the group level, we observed robust autonomic responses as well as spatially specific spinal cord BOLD responses at the expected location, but no spatial overlap in BOLD response patterns across days. While autonomic indicators of pain processing showed good-to-excellent reliability, both β-estimates and z-scores of task-related BOLD responses showed poor reliability across days in the target region (gray matter of the ipsilateral dorsal horn). When taking into account the sensitivity of gradient-echo echo planar imaging (GE-EPI) to draining vein signals by including the venous plexus in the analysis, we observed BOLD responses with fair reliability across days. Taken together, these results demonstrate that heat pain stimuli as short as one second are able to evoke a robust and spatially specific BOLD response, which is however strongly variable within participants across time, resulting in low reliability in the dorsal horn gray matter. Further improvements in data acquisition and analysis techniques are thus necessary before event-related spinal cord fMRI as used here can be reliably employed in longitudinal designs or clinical settings.

On the influence of the vascular architecture on Gradient Echo and Spin Echo BOLD fMRI signals across cortical depth: a simulation approach based on realistic 3D vascular networks

GE-BOLD contrast stands out as the predominant technique in functional MRI experiments for its high sensitivity and straightforward implementation. GE-BOLD exhibits rather similar sensitivity to vessels independent of their size at submillimeter resolution studies like those examining cortical columns and laminae. However, the presence of nonspecific macrovascular contributions poses a challenge to accurately isolate neuronal activity. SE-BOLD increases specificity towards small vessels, thereby enhancing its specificity to neuronal activity, due to the effective suppression of extravascular contributions caused by macrovessels with its refocusing pulse. However, even SE-BOLD measurements may not completely remove these macrovascular contributions. By simulating hemodynamic signals across cortical depth, we gain insights into vascular contributions to the laminar BOLD signal. In this study, we employed four realistic 3D vascular models to simulate oxygen saturation states in various vascular compartments, aiming to characterize both intravascular and extravascular contributions to GE and SE signals, and corresponding BOLD signal changes, across cortical depth at 7T. Simulations suggest that SE-BOLD cannot completely reduce the macrovascular contribution near the pial surface. Simulations also show that both the specificity and signal amplitude of BOLD signals at 7T depend on the spatial arrangement of large vessels throughout cortical depth and on the pial surface.

A fully synthetic three-dimensional human cerebrovascular model based on histological characteristics to investigate the hemodynamic fingerprint of the layer BOLD fMRI signal formation

Recent advances in functional magnetic resonance imaging (fMRI) at ultra-high field (≥7 tesla), novel hardware, and data analysis methods have enabled detailed research on neurovascular function, such as cortical layer-specific activity, in both human and nonhuman species. A widely used fMRI technique relies on the blood oxygen level-dependent (BOLD) signal. BOLD fMRI offers insights into brain function by measuring local changes in cerebral blood volume, cerebral blood flow, and oxygen metabolism induced by increased neuronal activity. Despite its potential, interpreting BOLD fMRI data is challenging as it is only an indirect measurement of neuronal activity. Computational modeling can help interpret BOLD data by simulating the BOLD signal formation. Current developments have focused on realistic 3D vascular models based on rodent data to understand the spatial and temporal BOLD characteristics. While such rodent-based vascular models highlight the impact of the angioarchitecture on the BOLD signal amplitude, anatomical differences between the rodent and human vasculature necessitate the development of human-specific models. Therefore, a computational framework integrating human cortical vasculature, hemodynamic changes, and biophysical properties is essential. Here, we present a novel computational approach: a three-dimensional VAscular MOdel based on Statistics (3D VAMOS), enabling the investigation of the hemodynamic fingerprint of the BOLD signal within a model encompassing a fully synthetic human 3D cortical vasculature and hemodynamics. Our algorithm generates microvascular and macrovascular architectures based on morphological and topological features from the literature on human cortical vasculature. By simulating specific oxygen saturation states and biophysical interactions, our framework characterizes the intravascular and extravascular signal contributions across cortical depth and voxel-wise levels for gradient-echo and spin-echo readouts. Thereby, the 3D VAMOS computational framework demonstrates that using human characteristics significantly affects the BOLD fingerprint, making it an essential step in understanding the fundamental underpinnings of layer-specific fMRI experiments.

Assessment of Small Vessel Function Using 7T MRI in Patients With Sporadic Cerebral Small Vessel Disease: The ZOOM@SVDs Study

BACKGROUND AND OBJECTIVES

Cerebral small vessel disease (cSVD) is a major cause of stroke and dementia, but little is known about disease mechanisms at the level of the small vessels. 7T-MRI allows assessing small vessel function in vivo in different vessel populations. We hypothesized that multiple aspects of small vessel function are altered in patients with cSVD and that these abnormalities relate to disease burden.

METHODS

Patients and controls participated in a prospective observational cohort study, the ZOOM@SVDs study. Small vessel function measures on 7T-MRI included perforating artery blood flow velocity and pulsatility index in the basal ganglia and centrum semiovale, vascular reactivity to visual stimulation in the occipital cortex, and reactivity to hypercapnia in the gray and white matter. Lesion load on 3T-MRI and cognitive function were used to assess disease burden.

RESULTS

Forty-six patients with sporadic cSVD (mean age ± SD 65 ± 9 years) and 22 matched controls (64 ± 7 years) participated in the ZOOM@SVDs study. Compared with controls, patients had increased pulsatility index (mean difference 0.09, p = 0.01) but similar blood flow velocity in basal ganglia perforating arteries and similar flow velocity and pulsatility index in centrum semiovale perforating arteries. The duration of the vascular response to brief visual stimulation in the occipital cortex was shorter in patients than in controls (mean difference -0.63 seconds, p = 0.02), whereas reactivity to hypercapnia was not significantly affected in the gray and total white matter. Among patients, reactivity to hypercapnia was lower in white matter hyperintensities compared with normal-appearing white matter (blood-oxygen-level dependent mean difference 0.35%, p = 0.001). Blood flow velocity and pulsatility index in basal ganglia perforating arteries and reactivity to brief visual stimulation correlated with disease burden.

DISCUSSION

We observed abnormalities in several aspects of small vessel function in patients with cSVD indicative of regionally increased arteriolar stiffness and decreased reactivity. Worse small vessel function also correlated with increased disease burden. These functional measures provide new mechanistic markers of sporadic cSVD.

Mapping cerebral perfusion in mice under various anesthesia levels using highly sensitive BOLD MRI with transient hypoxia

Cerebral perfusion is critical for the early detection of neurological diseases and for effectively monitoring disease progression and treatment responses. Mouse models are widely used in brain research, often under anesthesia, which can affect vascular physiology. However, the impact of anesthesia on regional cerebral blood volume and flow in mice has not been thoroughly investigated. In this study, we have developed a whole-brain perfusion MRI approach by using a 5-second nitrogen gas stimulus under inhalational anesthetics to induce transient BOLD dynamic susceptibility contrast (DSC). This method proved to be highly sensitive, repeatable within each imaging session, and across four weekly sessions. Relative cerebral blood volumes measured by BOLD DSC agree well with those by contrast agents. Quantitative cerebral blood volume and flow metrics were successfully measured in mice under dexmedetomidine and various isoflurane doses using both total vasculature-sensitive gradient-echo and microvasculature-sensitive spin-echo BOLD MRI. Dexmedetomidine reduces cerebral perfusion, while isoflurane increases cerebral perfusion in a dose-dependent manner.

Breath-hold BOLD fMRI without CO2 sampling enables estimation of venous cerebral blood volume: potential use in normalization of stimulus-evoked BOLD fMRI data

BOLD fMRI signal has been used in conjunction with vasodilatory stimulation as a marker of cerebrovascular reactivity (CVR): the relative change in cerebral blood flow (CBF) arising from a unit change in the vasodilatory stimulus. Using numerical simulations, we demonstrate that the variability in the relative BOLD signal change induced by vasodilation is strongly influenced by the variability in deoxyhemoglobin-containing cerebral blood volume (CBV), as this source of variability is likely to be more prominent than that of CVR. It may, therefore, be more appropriate to describe the relative BOLD signal change induced by an isometabolic vasodilation as a proxy of deoxygenated CBV (CBVdHb) rather than CVR. With this in mind, a new method was implemented to map a marker of CBVdHb, termed BOLD-CBV, based on the normalization of voxel-wise BOLD signal variation by an estimate of the intravascular venous BOLD signal from voxels filled with venous blood. The intravascular venous BOLD signal variation, recorded during repeated breath-holding, was extracted from the superior sagittal sinus in a cohort of 27 healthy volunteers and used as a regressor across the whole brain, yielding maps of BOLD-CBV. In the same cohort, we demonstrated the potential use of BOLD-CBV for the normalization of stimulus-evoked BOLD fMRI by comparing group-level BOLD fMRI responses to a visuomotor learning task with and without the inclusion of voxel-wise vascular covariates of BOLD-CBV and the BOLD signal change per mmHg variation in end-tidal carbon dioxide (BOLD-CVR). The empirical measure of BOLD-CBV accounted for more between-subject variability in the motor task-induced BOLD responses than BOLD-CVR estimated from end-tidal carbon dioxide recordings. The new method can potentially increase the power of group fMRI studies by including a measure of vascular characteristics and has the strong practical advantage of not requiring experimental measurement of end-tidal carbon dioxide, unlike traditional methods to estimate BOLD-CVR. It also more closely represents a specific physiological characteristic of brain vasculature than BOLD-CVR, namely blood volume.

The Choroid Plexus as an Alternative Locus for the Identification of the Arterial Input Function for Calculating Cerebral Perfusion Metrics Using MRI

BACKGROUND AND PURPOSE

MR imaging-based cerebral perfusion metrics can be obtained by tracing the passage of a bolus of contrast through the microvasculature of the brain parenchyma. Thus, the temporal signal pattern of the contrast agent is typically measured over a large artery such as the MCA to generate the arterial input function. The largest intracranial arteries in the brain may not always be suitable for selecting the arterial input function due to skull base susceptibility artifacts or reduced size from steno-occlusive disease. Therefore, a suitable alternative arterial input function window would be useful. The choroid plexus is a highly vascular tissue composed essentially of arterialized blood vessels and acellular stroma with low metabolic requirements relative to its blood flow and may be a suitable alternative to identify the arterial input function.

MATERIALS AND METHODS

We studied 8 healthy participants and 7 patients with gliomas who were administered a bolus of gadolinium. We selected an arterial input function from both the left and right M1 segments of the MCA and both lateral ventricles of the choroid plexus for each participant. We compared the changes in the T2* signal and the calculated resting perfusion metrics using the arterial input functions selected from the MCA and choroid plexus.

RESULTS

We found no systematic difference between resting perfusion metrics in GM and WM when calculated using an arterial input function from the MCA or choroid plexus in the same participant.

CONCLUSIONS

The choroid plexus provides an alternative location from which an arterial input function may be sampled when a suitable measure over an MCA is not available.

Multilayer Network Analysis across Cortical Depths in Resting-State 7T fMRI

In graph theory, "multilayer networks" represent systems involving several interconnected topological levels. A neuroscience example is the hierarchy of connections between different cortical depths or "lamina". This hierarchy is becoming non-invasively accessible in humans using ultra-high-resolution functional MRI (fMRI). Here, we applied multilayer graph theory to examine functional connectivity across different cortical depths in humans, using 7T fMRI (1-mm3 voxels; 30 participants). Blood oxygenation level dependent (BOLD) signals were derived from five depths between the white matter and pial surface. We then compared networks where the inter-regional connections were limited to a single cortical depth only ("layer-by-layer matrices") to those considering all possible connections between regions and cortical depths ("multilayer matrix"). We utilized global and local graph theory features that quantitatively characterize network attributes such as network composition, nodal centrality, path-based measures, and hub segregation. Detecting functional differences between cortical depths was improved using multilayer connectomics compared to the layer-by-layer versions. Superficial aspects of the cortex dominated information transfer and deeper aspects clustering. These differences were largest in frontotemporal and limbic brain regions. fMRI functional connectivity across different cortical depths may contain neurophysiologically relevant information. Multilayer connectomics could provide a methodological framework for studies on how information flows across this hierarchy.

Seven Tesla MRI in Alzheimer's disease research: State of the art and future directions: A narrative review

Seven tesla magnetic resonance imaging (7T MRI) is known to offer a superior spatial resolution and a signal-to-noise ratio relative to any other non-invasive imaging technique and provides the possibility for neuroimaging researchers to observe disease-related structural changes, which were previously only apparent on post-mortem tissue analyses. Alzheimer's disease is a natural and widely used subject for this technology since the 7T MRI allows for the anticipation of disease progression, the evaluation of secondary prevention measures thought to modify the disease trajectory, and the identification of surrogate markers for treatment outcome. In this editorial, we discuss the various neuroimaging biomarkers for Alzheimer's disease that have been studied using 7T MRI, which include morphological alterations, molecular characterization of cerebral T2*-weighted hypointensities, the evaluation of cerebral microbleeds and microinfarcts, biochemical changes studied with MR spectroscopy, as well as some other approaches. Finally, we discuss the limitations of the 7T MRI regarding imaging Alzheimer's disease and we provide our outlook for the future.

A mechanistic computational framework to investigate the hemodynamic fingerprint of the blood oxygenation level-dependent signal

Blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) is one of the most used imaging techniques to map brain activity or to obtain clinical information about human cortical vasculature, in both healthy and disease conditions. Nevertheless, BOLD fMRI is an indirect measurement of brain functioning triggered by neurovascular coupling. The origin of the BOLD signal is quite complex, and the signal formation thus depends, among other factors, on the topology of the cortical vasculature and the associated hemodynamic changes. To understand the hemodynamic evolution of the BOLD signal response in humans, it is beneficial to have a computational framework available that virtually resembles the human cortical vasculature, and simulates hemodynamic changes and corresponding MRI signal changes via interactions of intrinsic biophysical and magnetic properties of the tissues. To this end, we have developed a mechanistic computational framework that simulates the hemodynamic fingerprint of the BOLD signal based on a statistically defined, three-dimensional, vascular model that approaches the human cortical vascular architecture. The microvasculature is approximated through a Voronoi tessellation method and the macrovasculature is adapted from two-photon microscopy mice data. Using this computational framework, we simulated hemodynamic changes-cerebral blood flow, cerebral blood volume, and blood oxygen saturation-induced by virtual arterial dilation. Then we computed local magnetic field disturbances generated by the vascular topology and the corresponding blood oxygen saturation changes. This mechanistic computational framework also considers the intrinsic biophysical and magnetic properties of nearby tissue, such as water diffusion and relaxation properties, resulting in a dynamic BOLD signal response. The proposed mechanistic computational framework provides an integrated biophysical model that can offer better insights regarding the spatial and temporal properties of the BOLD signal changes.

Mitigating susceptibility-induced distortions in high-resolution 3DEPI fMRI at 7T

Geometric distortion is a major limiting factor for spatial specificity in high-resolution fMRI using EPI readouts and is exacerbated at higher field strengths due to increased B0 field inhomogeneity. Prominent correction schemes are based on B0 field-mapping or acquiring reverse phase-encoded (reversed-PE) data. However, to date, comparisons of these techniques in the context of fMRI have only been performed on 2DEPI data, either at lower field or lower resolution. In this study, we investigate distortion compensation in the context of sub-millimetre 3DEPI data at 7T. B0 field-mapping and reversed-PE distortion correction techniques were applied to both partial coverage BOLD-weighted and whole brain MT-weighted 3DEPI data with matched distortion. Qualitative assessment showed overall improvement in cortical alignment for both correction techniques in both 3DEPI fMRI and whole-brain MT-3DEPI datasets. The distortion-corrected MT-3DEPI images were quantitatively evaluated by comparing cortical alignment with an anatomical reference using dice coefficient (DC) and correlation ratio (CR) measures. These showed that B0 field-mapping and reversed-PE methods both improved correspondence between the MT-3DEPI and anatomical data, with more substantial improvements consistently obtained using the reversed-PE approach. Regional analyses demonstrated that the largest benefit of distortion correction, and in particular of the reversed-PE approach, occurred in frontal and temporal regions where susceptibility-induced distortions are known to be greatest, but had not led to complete signal dropout. In conclusion, distortion correction based on reversed-PE data has shown the greater capacity for achieving faithful alignment with anatomical data in the context of high-resolution fMRI at 7T using 3DEPI.

Layer-fMRI VASO with short stimuli and event-related designs at 7 T

Layers and columns are the dominant processing units in the human (neo)cortex at the mesoscopic scale. While the blood oxygenation dependent (BOLD) signal has a high detection sensitivity, it is biased towards unwanted signals from large draining veins at the cortical surface. The additional fMRI contrast of vascular space occupancy (VASO) has the potential to augment the neuroscientific interpretability of layer-fMRI results by means of capturing complementary information of locally specific changes in cerebral blood volume (CBV). Specifically, VASO is not subject to unwanted sensitivity amplifications of large draining veins. Because of constrained sampling efficiency, it has been mainly applied in combination with efficient block task designs and long trial durations. However, to study cognitive processes in neuroscientific contexts, or probe vascular reactivity, short stimulation periods are often necessary. Here, we developed a VASO acquisition procedure with a short acquisition period and sub-millimeter resolution. During visual event-related stimulation, we show reliable responses in visual cortices within a reasonable number of trials (∼20). Furthermore, the short TR and high spatial specificity of our VASO implementation enabled us to show differences in laminar reactivity and onset times. Finally, we explore the generalizability to a different stimulus modality (somatosensation). With this, we showed that CBV-sensitive VASO provides the means to capture layer-specific haemodynamic responses with high spatio-temporal resolution and is able to be used with event-related paradigms.

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.

Patterned edge-illuminated display for clinical examination of visual evoked potentials using simultaneous magnetic resonance imaging

PURPOSE

Commonly used technologies for visual pattern stimulation cannot operate in a magnetic resonance imaging room because they can interfere with the operation of the scanner and are vulnerable to its electromagnetic and magnetic fields. The aim of this single-center prospective observational study was to introduce a novel, structurally uncomplicated, easy-to-maintain, patterned edge-illuminated display (PEID) device for visual pattern-reversal stimulation, compare it with a commonly used cathode ray tube screen, and verify the equivalence of quantitative assays.

MATERIALS AND METHODS

The left and right eyes of 36 healthy participants with undilated pupils were examined on a commercial visual evoked potential (VEP) apparatus and on the PEID device, where pattern-reversal transient VEPs were elicited by checkerboard stimuli with large (0.89°; 0.86°-0.92°) and small (0.21°; 0.20°-0.23°) checks.

RESULTS

The PEID device demonstrated the required reliability and dynamic characteristics, as well as precise time-locking required for a VEP diagnosis. The results of Deming's correlation analysis showed that both the commercial cathode ray tube monitor and the PEID device produced identical VEP results within the context of experimental uncertainty. The standard deviation of Deming's regression may indicate the uncertainty of the VEPs measured in clinical practice. The Bland-Altman analysis of the mean showed no significant difference in the amplitude and peak time of VEPs measured on the PEID device compared to that of the commercial cathode ray tube monitor.

CONCLUSIONS

The presented PEID device meets all the required standards and can be easily installed in various types of commercial magnetic resonance imaging scanners.

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.

Improved brain stimulation targeting by optimising image acquisition parameters

Functional connectivity analysis from rs-fMRI data has been used for determining cortical targets in therapeutic applications of non-invasive brain stimulation using transcranial magnetic stimulation (TMS). Reliable connectivity measures are therefore essential for every rs-fMRI-based TMS targeting approach. Here, we examine the effect of echo time (TE) on the reproducibility and spatial variability of resting-state connectivity measures. We acquired multiple runs of single-echo fMRI data with either short (TE = 30 ms) or long (TE = 38 ms) echo time to investigate inter-run spatial reproducibility of a clinically relevant functional connectivity map, i.e., originating from the sgACC. We find that connectivity maps obtained from TE = 38 ms rs-fMRI data are significantly more reliable than those obtained from TE = 30 ms data sets. Our results clearly show that optimizing sequence parameters can be beneficial for ensuring high-reliability resting-state acquisition protocols to be used for TMS targeting. The differences between reliability in connectivity measures for different TEs could inform future clinical research in optimising MR sequences.

Miniaturized MR-compatible ultrasound system for real-time monitoring of acoustic effects in mice using high-resolution MRI

Visualization of focused ultrasound in high spatial and temporal resolution is crucial for accurately and precisely targeting brain regions noninvasively. Magnetic resonance imaging (MRI) is the most widely used noninvasive tool for whole-brain imaging. However, focused ultrasound studies employing high-resolution (> 9.4 T) MRI in small animals are limited by the small size of the radiofrequency (RF) volume coil and the noise sensitivity of the image to external systems such as bulky ultrasound transducers. This technical note reports a miniaturized ultrasound transducer system packaged directly above a mouse brain for monitoring ultrasound-induced effects using high-resolution 9.4 T MRI. Our miniaturized system integrates MR-compatible materials with electromagnetic (EM) noise reduction techniques to demonstrate echo-planar imaging (EPI) signal changes in the mouse brain at various ultrasound acoustic intensities. The proposed ultrasound-MRI system will enable extensive research in the expanding field of ultrasound therapeutics.

Imaging the columnar functional organization of human area MT+ to axis-of-motion stimuli using VASO at 7 Tesla

Cortical columns of direction-selective neurons in the motion sensitive area (MT) have been successfully established as a microscopic feature of the neocortex in animals. The same property has been investigated at mesoscale (<1 mm) in the homologous brain area (hMT+, V5) in living humans by using ultra-high field functional magnetic resonance imaging (fMRI). Despite the reproducibility of the selective response to axis-of-motion stimuli, clear quantitative evidence for the columnar organization of hMT+ is still lacking. Using cerebral blood volume (CBV)-sensitive fMRI at 7 Tesla with submillimeter resolution and high spatial specificity to microvasculature, we investigate the columnar functional organization of hMT+ in 5 participants perceiving axis-of-motion stimuli for both blood oxygenation level dependent (BOLD) and vascular space occupancy (VASO) contrast mechanisms provided by the used slice-selective slab-inversion (SS-SI)-VASO sequence. With the development of a new searchlight algorithm for column detection, we provide the first quantitative columnarity map that characterizes the entire 3D hMT+ volume. Using voxel-wise measures of sensitivity and specificity, we demonstrate the advantage of using CBV-sensitive fMRI to detect mesoscopic cortical features by revealing higher specificity of axis-of-motion cortical columns for VASO as compared to BOLD contrast. These voxel-wise metrics also provide further insights on how to mitigate the highly debated draining veins effect. We conclude that using CBV-VASO fMRI together with voxel-wise measurements of sensitivity, specificity and columnarity offers a promising avenue to quantify the mesoscopic organization of hMT+ with respect to axis-of-motion stimuli. Furthermore, our approach and methodological developments are generalizable and applicable to other human brain areas where similar mesoscopic research questions are addressed.

Transfer function analysis assesses resting cerebral perfusion metrics using hypoxia-induced deoxyhemoglobin as a contrast agent

Introduction: Use of contrast in determining hemodynamic measures requires the deconvolution of an arterial input function (AIF) selected over a voxel in the middle cerebral artery to calculate voxel wise perfusion metrics. Transfer function analysis (TFA) offers an alternative analytic approach that does not require identifying an AIF. We hypothesised that TFA metrics Gain, Lag, and their ratio, Gain/Lag, correspond to conventional AIF resting perfusion metrics relative cerebral blood volume (rCBV), mean transit time (MTT) and relative cerebral blood flow (rCBF), respectively. Methods: 24 healthy participants (17 M) and 1 patient with steno-occlusive disease were recruited. We used non-invasive transient hypoxia-induced deoxyhemoglobin as an MRI contrast. TFA and conventional AIF analyses were used to calculate averages of whole brain and smaller regions of interest. Results: Maps of these average metrics had colour scales adjusted to enhance contrast and identify areas of high congruence. Regional gray matter/white matter (GM/WM) ratios for MTT and Lag, rCBF and Gain/Lag, and rCBV and Gain were compared. The GM/WM ratios were greater for TFA metrics compared to those from AIF analysis indicating an improved regional discrimination. Discussion: Resting perfusion measures generated by The BOLD analysis resulting from a transient hypoxia induced variations in deoxyhemoglobin analyzed by TFA are congruent with those analyzed by conventional AIF analysis.

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 T₂* at 9.4 T while VASO from longer T₁. 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.

Combining arterial blood contrast with BOLD increases fMRI intracortical contrast

BOLD fMRI is widely applied in human neuroscience but is limited in its spatial specificity due to a cortical-depth-dependent venous bias. This reduces its localization specificity with respect to neuronal responses, a disadvantage for neuroscientific research. Here, we modified a submillimeter BOLD protocol to selectively reduce venous and tissue signal and increase cerebral blood volume weighting through a pulsed saturation scheme (dubbed Arterial Blood Contrast) at 7 T. Adding Arterial Blood Contrast on top of the existing BOLD contrast modulated the intracortical contrast. Isolating the Arterial Blood Contrast showed a response free of pial-surface bias. The results suggest that Arterial Blood Contrast can modulate the typical fMRI spatial specificity, with important applications in in-vivo neuroscience.

Effect of motion, cortical orientation and spatial resolution on quantitative imaging of cortical R2* and magnetic susceptibility at 0.3 mm in-plane resolution at 7 T

MR images of the effective relaxation rate R2* and magnetic susceptibility χ derived from multi-echo T2*-weighted (T2*w) MRI can provide insight into iron and myelin distributions in the brain, with the potential of providing biomarkers for neurological disorders. Quantification of R2* and χ at submillimeter resolution in the cortex in vivo has been difficult because of challenges such as head motion, limited signal to noise ratio, long scan time, and motion related magnetic field fluctuations. This work aimed to improve the robustness for quantifying intracortical R2* and χ and analyze the effects from motion, spatial resolution, and cortical orientation. T2*w data was acquired with a spatial resolution of 0.3 × 0.3 × 0.4 mm3 at 7 T and downsampled to various lower resolutions. A combined correction for motion and B0 changes was deployed using volumetric navigators. Such correction improved the T2*w image quality rated by experienced image readers and test-retest reliability of R2* and χ quantification with reduced median inter-scan differences up to 10 s-1 and 5 ppb, respectively. R2* and χ near the line of Gennari, a cortical layer high in iron and myelin, were as much as 10 s-1 and 10 ppb higher than the region at adjacent cortical depth. In addition, a significant effect due to the cortical orientation relative to the static field (B0) was observed in χ with a peak-to-peak amplitude of about 17 ppb. In retrospectively downsampled data, the capability to distinguish different cortical depth regions based on R2* or χ contrast remained up to isotropic 0.5 mm resolution. This study highlights the unique characteristics of R2* and χ along the cortical depth at submillimeter resolution and the need for motion and B0 corrections for their robust quantification in vivo.

A vision of 14 T MR for fundamental and clinical science

OBJECTIVE

We outline our vision for a 14 Tesla MR system. This comprises a novel whole-body magnet design utilizing high temperature superconductor; a console and associated electronic equipment; an optimized radiofrequency coil setup for proton measurement in the brain, which also has a local shim capability; and a high-performance gradient set.

RESEARCH FIELDS

The 14 Tesla system can be considered a 'mesocope': a device capable of measuring on biologically relevant scales. In neuroscience the increased spatial resolution will anatomically resolve all layers of the cortex, cerebellum, subcortical structures, and inner nuclei. Spectroscopic imaging will simultaneously measure excitatory and inhibitory activity, characterizing the excitation/inhibition balance of neural circuits. In medical research (including brain disorders) we will visualize fine-grained patterns of structural abnormalities and relate these changes to functional and molecular changes. The significantly increased spectral resolution will make it possible to detect (dynamic changes in) individual metabolites associated with pathological pathways including molecular interactions and dynamic disease processes.

CONCLUSIONS

The 14 Tesla system will offer new perspectives in neuroscience and fundamental research. We anticipate that this initiative will usher in a new era of ultra-high-field MR.

Improved sensitivity and microvascular weighting of 3T laminar fMRI with GE-BOLD using NORDIC and phase regression

INTRODUCTION

Functional MRI with spatial resolution in the submillimeter domain enables measurements of activation across cortical layers in humans. This is valuable as different types of cortical computations, e.g., feedforward versus feedback related activity, take place in different cortical layers. Laminar fMRI studies have almost exclusively employed 7T scanners to overcome the reduced signal stability associated with small voxels. However, such systems are relatively rare and only a subset of those are clinically approved. In the present study, we examined if the feasibility of laminar fMRI at 3T could be improved by use of NORDIC denoising and phase regression.

METHODS

5 healthy subjects were scanned on a Siemens MAGNETOM Prisma 3T scanner. To assess across-session reliability, each subject was scanned in 3-8 sessions on 3-4 consecutive days. A 3D gradient echo EPI (GE-EPI) sequence was used for BOLD acquisitions (voxel size 0.82 mm isotopic, TR = 2.2 s) using a block design finger tapping paradigm. NORDIC denoising was applied to the magnitude and phase time series to overcome limitations in temporal signal-to-noise ratio (tSNR) and the denoised phase time series were subsequently used to correct for large vein contamination through phase regression.

RESULTS AND CONCLUSION

NORDIC denoising resulted in tSNR values comparable to or higher than commonly observed at 7T. Layer-dependent activation profiles could thus be extracted robustly, within and across sessions, from regions of interest located in the hand knob of the primary motor cortex (M1). Phase regression led to substantially reduced superficial bias in obtained layer profiles, although residual macrovascular contribution remained. We believe the present results support an improved feasibility of laminar fMRI at 3T.

The many layers of BOLD. The effect of hypercapnic and hyperoxic stimuli on macro- and micro-vascular compartments quantified by CVR, M, and CBV across cortical depth

Ultra-high field functional magnetic resonance imaging (fMRI) offers the spatial resolution to measure neuronal activity at the scale of cortical layers. However, cortical depth dependent vascularization differences, such as a higher prevalence of macro-vascular compartments near the pial surface, have a confounding effect on depth-resolved blood-oxygen-level dependent (BOLD) fMRI signals. In the current study, we use hypercapnic and hyperoxic breathing conditions to quantify the influence of all venous vascular and micro-vascular compartments on laminar BOLD fMRI, as measured with gradient-echo (GE) and spin-echo (SE) scan sequences, respectively. We find that all venous vascular and micro-vascular compartments are capable of comparable theoretical maximum signal intensities, as represented by the M-value parameter. However, the capacity for vessel dilation, as reflected by the cerebrovascular reactivity (CVR), is approximately two and a half times larger for all venous vascular compartments combined compared to the micro-vasculature at superficial layers. Finally, there is roughly a 35% difference in estimates of CBV changes between all venous vascular and micro-vascular compartments, although this relative difference was approximately uniform across cortical depth. Thus, our results suggest that fMRI BOLD signal differences across cortical depth are likely caused by differences in dilation properties between macro- and micro-vascular compartments.

DSC MRI in the human brain using deoxyhemoglobin and gadolinium-Simulations and validations at 3T

INTRODUCTION

Dynamic susceptibility contrast (DSC) MRI allows clinicians to determine perfusion parameters in the brain, such as cerebral blood flow, cerebral blood volume, and mean transit time. To enable quantification, susceptibility changes can be induced using gadolinium (Gd) or deoxyhemoglobin (dOHb), the latter just recently introduced as a contrast agent in DSC. Previous investigations found that experimental parameters and analysis choices, such as the susceptibility amplitude and partial volume, affect perfusion quantification. However, the accuracy and precision of DSC MRI has not been systematically investigated, particularly in the lower susceptibility range.

METHODS

In this study, we compared perfusion values determined using Gd with values determined using a contrast agent with a lower susceptibility-dOHb-under different physiological conditions, such as varying the baseline blood oxygenation and/or magnitude of hypoxic bolus, by utilizing numerical simulations and conducting experiments on healthy subjects at 3T. The simulation framework we developed for DSC incorporates MRI signal contributions from intravascular and extravascular proton spins in arterial, venous, and cerebral tissue voxels. This framework allowed us to model the MRI signal in response to both Gd and dOHb.

RESULTS AND DISCUSSION

We found, both in the experimental results and simulations, that a reduced intravascular volume of the selected arterial voxel, reduced baseline oxygen saturation, greater susceptibility of applied contrast agent (Gd vs. dOHb), and/or larger magnitude of applied hypoxic bolus reduces the overestimation and increases precision of cerebral blood volume and flow. As well, we found that normalizing tissue to venous rather than arterial signal increases the accuracy of perfusion quantification across experimental paradigms. Furthermore, we found that shortening the bolus duration increases the accuracy and reduces the calculated values of mean transit time. In summary, we experimentally uncovered an array of perfusion quantification dependencies, which agreed with the simulation framework predictions, using a wider range of susceptibility values than previously investigated. We argue for caution when comparing absolute and relative perfusion values within and across subjects obtained from a standard DSC MRI analysis, particularly when employing different experimental paradigms and contrast agents.

Cerebral blood volume sensitive layer-fMRI in the human auditory cortex at 7T: Challenges and capabilities

The development of ultra high field fMRI signal readout strategies and contrasts has led to the possibility of imaging the human brain in vivo and non-invasively at increasingly higher spatial resolutions of cortical layers and columns. One emergent layer-fMRI acquisition method with increasing popularity is the cerebral blood volume sensitive sequence named vascular space occupancy (VASO). This approach has been shown to be mostly sensitive to locally-specific changes of laminar microvasculature, without unwanted biases of trans-laminar draining veins. Until now, however, VASO has not been applied in the technically challenging cortical area of the auditory cortex. Here, we describe the main challenges we encountered when developing a VASO protocol for auditory neuroscientific applications and the solutions we have adopted. With the resulting protocol, we present preliminary results of laminar responses to sounds and as a proof of concept for future investigations, we map the topographic representation of frequency preference (tonotopy) in the auditory cortex.

An optimized b-value sampling for the quantification of interstitial fluid using diffusion-weighted MRI, a genetic algorithm approach

PURPOSE

Multi-b-value diffusion-weighted MRI techniques can simultaneously measure the parenchymal diffusivity, microvascular perfusion, and a third, intermediate diffusion component. This component is related to the interstitial fluid in the brain parenchyma. However, simultaneously estimating three diffusion components from multi-b-value data is difficult and has strong dependence on SNR and chosen b-values. As the number of acquired b-values is limited due to scanning time, it is important to know which b-values are most effective to be included. Therefore, this study evaluates an optimized b-value sampling for interstitial fluid estimation.

METHOD

The optimized b-value sampling scheme is determined using a genetic algorithm. Subsequently, the performance of this optimized sampling is assessed by comparing it with a linear, logarithmic, and previously proposed sampling scheme, in terms of the RMS error (RMSE) for the intermediate component estimation. The in vivo performance of the optimized sampling is assessed using 7T data with 101 equally spaced b-values ranging from 0 to 1000 s/mm² . In this case, the RMSE was determined by comparing the fit that includes all b-values.

RESULTS

The optimized b-value sampling for estimating the intermediate component was reported to be [0, 30, 90, 210, 280, 350, 580, 620, 660, 680, 720, 760, 980, 990, 1000] s/mm² . For computer simulations, the optimized sampling had a lower RMSE, compared with the other samplings for varying levels of SNR. For the in vivo data, the voxel-wise RMSE of the optimized sampling was lower compared with other sampling schemes.

CONCLUSION

The genetic algorithm-optimized b-value scheme improves the quantification of the diffusion component related to interstitial fluid in terms of a lower RMSE.

Cortical depth profiles of auditory and visual 7 T functional MRI responses in human superior temporal areas

Invasive neurophysiological studies in nonhuman primates have shown different laminar activation profiles to auditory vs. visual stimuli in auditory cortices and adjacent polymodal areas. Means to examine the underlying feedforward vs. feedback type influences noninvasively have been limited in humans. Here, using 1-mm isotropic resolution 3D echo-planar imaging at 7 T, we studied the intracortical depth profiles of functional magnetic resonance imaging (fMRI) blood oxygenation level dependent (BOLD) signals to brief auditory (noise bursts) and visual (checkerboard) stimuli. BOLD percent-signal-changes were estimated at 11 equally spaced intracortical depths, within regions-of-interest encompassing auditory (Heschl's gyrus, Heschl's sulcus, planum temporale, and posterior superior temporal gyrus) and polymodal (middle and posterior superior temporal sulcus) areas. Effects of differing BOLD signal strengths for auditory and visual stimuli were controlled via normalization and statistical modeling. The BOLD depth profile shapes, modeled with quadratic regression, were significantly different for auditory vs. visual stimuli in auditory cortices, but not in polymodal areas. The different depth profiles could reflect sensory-specific feedforward versus cross-sensory feedback influences, previously shown in laminar recordings in nonhuman primates. The results suggest that intracortical BOLD profiles can help distinguish between feedforward and feedback type influences in the human brain. Further experimental studies are still needed to clarify how underlying signal strength influences BOLD depth profiles under different stimulus conditions.

Towards functional spin-echo BOLD line-scanning in humans at 7T

OBJECTIVE

Neurons cluster into sub-millimeter spatial structures and neural activity occurs at millisecond resolutions; hence, ultimately, high spatial and high temporal resolutions are required for functional MRI. In this work, we implemented a spin-echo line-scanning (SELINE) sequence to use in high spatial and temporal resolution fMRI.

MATERIALS AND METHODS

A line is formed by simply rotating the spin-echo refocusing gradient to a plane perpendicular to the excited slice and by removing the phase-encoding gradient. This technique promises a combination of high spatial and temporal resolution (250 μm, 500 ms) and microvascular specificity of functional responses. We compared SELINE data to a corresponding gradient-echo version (GELINE).

RESULTS

We demonstrate that SELINE showed much-improved line selection (i.e. a sharper line profile) compared to GELINE, albeit at the cost of a significant drop in functional sensitivity.

DISCUSSION

This low functional sensitivity needs to be addressed before SELINE can be applied for neuroscientific purposes.

Improved Selectivity in 7 T Digit Mapping Using VASO-CBV

Functional magnetic resonance imaging (fMRI) at Ultra-high field (UHF, ≥ 7 T) benefits from significant gains in the BOLD contrast-to-noise ratio (CNR) and temporal signal-to-noise ratio (tSNR) compared to conventional field strengths (3 T). Although these improvements enabled researchers to study the human brain to unprecedented spatial resolution, the blood pooling effect reduces the spatial specificity of the widely-used gradient-echo BOLD acquisitions. In this context, vascular space occupancy (VASO-CBV) imaging may be advantageous since it is proposed to have a higher spatial specificity than BOLD. We hypothesized that the assumed higher specificity of VASO-CBV imaging would translate to reduced overlap in fine-scale digit representation maps compared to BOLD-based digit maps. We used sub-millimeter resolution VASO fMRI at 7 T to map VASO-CBV and BOLD responses simultaneously in the motor and somatosensory cortices during individual finger movement tasks. We assessed the cortical overlap in different ways, first by calculating similarity coefficient metrics (DICE and Jaccard) and second by calculating selectivity measures. In addition, we demonstrate a consistent topographical organization of the targeted digit representations (thumb-index-little finger) in the motor areas. We show that the VASO-CBV responses yielded less overlap between the digit clusters than BOLD, and other selectivity measures were higher for VASO-CBV too. In summary, these results were consistent across metrics and participants, confirming the higher spatial specificity of VASO-CBV compared to BOLD.

CADASIL Affects Multiple Aspects of Cerebral Small Vessel Function on 7T-MRI

OBJECTIVE

Cerebral small vessel diseases (cSVDs) are a major cause of stroke and dementia. We used cutting-edge 7T-MRI techniques in patients with Cerebral Autosomal Dominant Arteriopathy with Subcortical Infarcts and Leukoencephalopathy (CADASIL), to establish which aspects of cerebral small vessel function are affected by this monogenic form of cSVD.

METHODS

We recruited 23 CADASIL patients (age 51.1 ± 10.1 years, 52% women) and 13 age- and sex-matched controls (46.1 ± 12.6, 46% women). Small vessel function measures included: basal ganglia and centrum semiovale perforating artery blood flow velocity and pulsatility, vascular reactivity to a visual stimulus in the occipital cortex and reactivity to hypercapnia in the cortex, subcortical gray matter, white matter, and white matter hyperintensities.

RESULTS

Compared with controls, CADASIL patients showed lower blood flow velocity and higher pulsatility index within perforating arteries of the centrum semiovale (mean difference - 0.09 cm/s, p = 0.03 and 0.20, p = 0.009) and basal ganglia (mean difference - 0.98 cm/s, p = 0.003 and 0.17, p = 0.06). Small vessel reactivity to a short visual stimulus was decreased (blood-oxygen-level dependent [BOLD] mean difference -0.21%, p = 0.04) in patients, while reactivity to hypercapnia was preserved in the cortex, subcortical gray matter, and normal appearing white matter. Among patients, reactivity to hypercapnia was decreased in white matter hyperintensities compared to normal appearing white matter (BOLD mean difference -0.29%, p = 0.02).

INTERPRETATION

Multiple aspects of cerebral small vessel function on 7T-MRI were abnormal in CADASIL patients, indicative of increased arteriolar stiffness and regional abnormalities in reactivity, locally also in relation to white matter injury. These observations provide novel markers of cSVD for mechanistic and intervention studies. ANN NEUROL 2023;93:29-39.

Improved laminar specificity and sensitivity by combining SE and GE BOLD signals

The most widely used gradient-echo (GE) blood oxygenation level-dependent (BOLD) contrast has high sensitivity, but low specificity due to draining vein contributions, while spin-echo (SE) BOLD approach at ultra-high magnetic fields is highly specific to neural active sites but has lower sensitivity. To obtain high specificity and sensitivity, we propose to utilize a vessel-size-sensitive filter to the GE-BOLD signal, which suppresses macrovascular contributions and to combine selectively retained microvascular GE-BOLD signals with the SE-BOLD signals. To investigate our proposed idea, fMRI with 0.8 mm isotropic resolution was performed on the primary motor and sensory cortices in humans at 7 T by implementing spin- and gradient-echo (SAGE) echo planar imaging (EPI) acquisition. Microvascular-passed sigmoidal filters were designed based upon the vessel-size-sensitive ΔR₂*/ΔR₂ value for retaining GE-BOLD signals originating from venous vessels with ≤ 45 μm and ≤ 65 μm diameter. Unlike GE-BOLD fMRI, the laminar profile of SAGE-BOLD fMRI with the vessel-size-sensitive filter peaked at ∼ 1.0 mm from the surface of the primary motor and sensory cortices, demonstrating an improvement of laminar specificity over GE-BOLD fMRI. Also, the functional sensitivity of SAGE BOLD at middle layers (0.75-1.5 mm) was improved by ∼ 80% to ∼100% when compared with SE BOLD. In summary, we showed that combined GE- and SE-BOLD fMRI with the vessel-size-sensitive filter indeed yielded improved laminar specificity and sensitivity and is therefore an excellent tool for high spatial resolution ultra-high filed (UHF)-fMRI studies for resolving mesoscopic functional units.

Magnetization transfer weighted laminar fMRI with multi-echo FLASH

Laminar functional magnetic resonance imaging (fMRI) using the gradient echo (GRE) blood oxygenation level dependent (BOLD) contrast is prone to signal changes arising from large unspecific venous vessels. Alternatives based on changes of cerebral blood volume (CBV) become more popular since it is expected that this hemodynamic response is dominant in microvasculature. One approach to sensitize the signal toward changes in CBV, and to simultaneously reduce unwanted extravascular (EV) BOLD blurring, is to selectively reduce gray matter (GM) signal via magnetization transfer (MT). In this work, we use off-resonant MT-pulses with a 3D FLASH readout to perform MT-prepared (MT-prep) laminar fMRI of the primary visual cortex (V1) at multiple echo times at 7 T. With a GRE-BOLD contrast without additional MT-weighting as reference, we investigated the influence of the MT-preparation on the shape and the echo time dependency of laminar profiles. Through numerical simulations, we optimized the sequence parameters to increase the sensitivity toward signal changes induced by changes in arterial CBV and to delineate the contributions of different compartments to the signal. We show that at 7 T, GM signals can be reduced by 30 %. Our laminar fMRI responses exhibit an increased signal change in the parenchyma at very short TE compared to a BOLD-only reference as a result of reduced EV signal intensity. By varying echo times, we could show that MT-prep results in less sensitivity toward unwanted signal changes based on changes in T₂*. We conclude that when accounting for nuclear overhauser enhancement effects in blood, off-resonant MT-prep combined with efficient short TE readouts can become a promising method to reduce unwanted EV venous contributions in GRE-BOLD and/or to allow scanning at much shorter echo times without incurring a sensitivity penalty in laminar fMRI.

Whole-brain perfusion mapping in mice by dynamic BOLD MRI with transient hypoxia

Non-invasive mapping of cerebral perfusion is critical for understanding neurovascular and neurodegenerative diseases. However, perfusion MRI methods cannot be easily implemented for whole-brain studies in mice because of their small size. To overcome this issue, a transient hypoxia stimulus was applied to induce a bolus of deoxyhemoglobins as an endogenous paramagnetic contrast in blood oxygenation level-dependent (BOLD) MRI. Based on stimulus-duration-dependent studies, 5 s anoxic stimulus was chosen, which induced a decrease in arterial oxygenation to 59%. Dynamic susceptibility changes were acquired with whole-brain BOLD MRI using both all-vessel-sensitive gradient-echo and microvascular-sensitive spin-echo readouts. Cerebral blood flow (CBF) and cerebral blood volume (CBV) were quantified by modeling BOLD dynamics using a partial-volume-corrected arterial input function. In the mouse under ketamine/xylazine anesthesia, total CBF and CBV were 112.0 ± 15.0 ml/100 g/min and 3.39 ± 0.59 ml/100 g (n = 15 mice), respectively, whereas microvascular CBF and CBV were 85.8 ± 6.9 ml/100 g/min and 2.23 ± 0.27 ml/100 g (n = 7 mice), respectively. Regional total vs. microvascular perfusion metrics were highly correlated but a slight mismatch was observed in the large-vessel areas and cortical depth profiles. Overall, this non-invasive, repeatable, simple hypoxia BOLD-MRI approach is viable for perfusion mapping of rodents.

Mesoscopic in vivo human T2* dataset acquired using quantitative MRI at 7 Tesla

Mesoscopic (0.1-0.5 mm) interrogation of the living human brain is critical for advancing neuroscience and bridging the resolution gap with animal models. Despite the variety of MRI contrasts measured in recent years at the mesoscopic scale, in vivo quantitative imaging of T₂^* has not been performed. Here we provide a dataset containing empirical T₂^* measurements acquired at 0.35 × 0.35 × 0.35 mm³ voxel resolution using 7 Tesla MRI. To demonstrate unique features and high quality of this dataset, we generate flat map visualizations that reveal fine-scale cortical substructures such as layers and vessels, and we report quantitative depth-dependent T₂^* (as well as R₂^*) values in primary visual cortex and auditory cortex that are highly consistent across subjects. This dataset is freely available at https://doi.org/10.17605/OSF.IO/N5BJ7, and may prove useful for anatomical investigations of the human brain, as well as for improving our understanding of the basis of the T₂^*-weighted (f)MRI signal.

Resting-state functional connectivity in adults with 47,XXX: a 7 Tesla MRI study

Triple X syndrome is a sex chromosomal aneuploidy characterized by the presence of a supernumerary X chromosome, resulting in a karyotype of 47,XXX in affected females. It has been associated with a variable cognitive, behavioral, and psychiatric phenotype, but little is known about its effects on brain function. We therefore conducted 7 T resting-state functional magnetic resonance imaging and compared data of 19 adult individuals with 47,XXX and 21 age-matched healthy control women using independent component analysis and dual regression. Additionally, we examined potential relationships between social cognition and social functioning scores, and IQ, and mean functional connectivity values. The 47,XXX group showed significantly increased functional connectivity of the fronto-parietal resting-state network with the right postcentral gyrus. Resting-state functional connectivity (rsFC) variability was not associated with IQ and social cognition and social functioning deficits in the participants with 47,XXX. We thus observed an effect of a supernumerary X chromosome in adult women on fronto-parietal rsFC. These findings provide additional insight into the role of the X chromosome on functional connectivity of the brain. Further research is needed to understand the clinical implications of altered rsFC in 47,XXX.

Characterization of the blood oxygen level dependent hemodynamic response function in human subcortical regions with high spatiotemporal resolution

Subcortical brain regions are absolutely essential for normal human function. These phylogenetically early brain regions play critical roles in human behaviors such as the orientation of attention, arousal, and the modulation of sensory signals to cerebral cortex. Despite the critical health importance of subcortical brain regions, there has been a dearth of research on their neurovascular responses. Blood oxygen level dependent (BOLD) functional MRI (fMRI) experiments can help fill this gap in our understanding. The BOLD hemodynamic response function (HRF) evoked by brief (<4 s) neural activation is crucial for the interpretation of fMRI results because linear analysis between neural activity and the BOLD response relies on the HRF. Moreover, the HRF is a consequence of underlying local blood flow and oxygen metabolism, so characterization of the HRF enables understanding of neurovascular and neurometabolic coupling. We measured the subcortical HRF at 9.4T and 3T with high spatiotemporal resolution using protocols that enabled reliable delineation of HRFs in individual subjects. These results were compared with the HRF in visual cortex. The HRF was faster in subcortical regions than cortical regions at both field strengths. There was no significant undershoot in subcortical areas while there was a significant post-stimulus undershoot that was tightly coupled with its peak amplitude in cortex. The different BOLD temporal dynamics indicate different vascular dynamics and neurometabolic responses between cortex and subcortical nuclei.

Laminar functional magnetic resonance imaging in vision research

Magnetic resonance imaging (MRI) scanners at ultra-high magnetic fields have become available to use in humans, thus enabling researchers to investigate the human brain in detail. By increasing the spatial resolution, ultra-high field MR allows both structural and functional characterization of cortical layers. Techniques that can differentiate cortical layers, such as histological studies and electrode-based measurements have made critical contributions to the understanding of brain function, but these techniques are invasive and thus mainly available in animal models. There are likely to be differences in the organization of circuits between humans and even our closest evolutionary neighbors. Thus research on the human brain is essential. Ultra-high field MRI can observe differences between cortical layers, but is non-invasive and can be used in humans. Extensive previous literature has shown that neuronal connections between brain areas that transmit feedback and feedforward information terminate in different layers of the cortex. Layer-specific functional MRI (fMRI) allows the identification of layer-specific hemodynamic responses, distinguishing feedback and feedforward pathways. This capability has been particularly important for understanding visual processing, as it has allowed researchers to test hypotheses concerning feedback and feedforward information in visual cortical areas. In this review, we provide a general overview of successful ultra-high field MRI applications in vision research as examples of future research.

Modelling the depth-dependent VASO and BOLD responses in human primary visual cortex

Functional magnetic resonance imaging (fMRI) using a blood-oxygenation-level-dependent (BOLD) contrast is a common method for studying human brain function noninvasively. Gradient-echo (GRE) BOLD is highly sensitive to the blood oxygenation change in blood vessels; however, the spatial signal specificity can be degraded due to signal leakage from activated lower layers to superficial layers in depth-dependent (also called laminar or layer-specific) fMRI. Alternatively, physiological variables such as cerebral blood volume using the VAscular-Space-Occupancy (VASO) contrast have shown higher spatial specificity compared to BOLD. To better understand the physiological mechanisms such as blood volume and oxygenation changes and to interpret the measured depth-dependent responses, models are needed which reflect vascular properties at this scale. For this purpose, we extended and modified the "cortical vascular model" previously developed to predict layer-specific BOLD signal changes in human primary visual cortex to also predict a layer-specific VASO response. To evaluate the model, we compared the predictions with experimental results of simultaneous VASO and BOLD measurements in a group of healthy participants. Fitting the model to our experimental data provided an estimate of CBV change in different vascular compartments upon neural activity. We found that stimulus-evoked CBV change mainly occurs in small arterioles, capillaries, and intracortical arteries and that the contribution from venules and ICVs is smaller. Our results confirm that VASO is less susceptible to large vessel effects compared to BOLD, as blood volume changes in intracortical arteries did not substantially affect the resulting depth-dependent VASO profiles, whereas depth-dependent BOLD profiles showed a bias towards signal contributions from intracortical veins.

Static and dynamic BOLD fMRI components along white matter fibre tracts and their dependence on the orientation of the local diffusion tensor axis relative to the B0-field

Recent studies have reported functional MRI (fMRI) activation within cerebral white matter (WM) using blood-oxygenation-level-dependent (BOLD) contrast. Many blood vessels in WM run parallel to the fibre bundles, and other studies observed dependence of susceptibility contrast-based measures of blood volume on the local orientation of the fibre bundles relative to the magnetic field or B₀ axis. Motivated by this, we characterized the dependence of gradient-echo BOLD fMRI on fibre orientation (estimated by the local diffusion tensor) relative to the B₀ axis to test whether the alignment between bundles and vessels imparts an orientation dependence on resting-state BOLD fluctuations in the WM. We found that the baseline signal level of the T₂*-weighted data is 11% higher in voxels containing fibres parallel to B₀ than those containing perpendicular fibres, consistent with a static influence of either fibre or vessel orientation on local T₂* values. We also found that BOLD fluctuations in most bundles exhibit orientation effects expected from oxygenation changes, with larger amplitudes from voxels containing perpendicular fibres. Different magnitudes of this orientation effect were observed across the major WM bundles, with inferior fasciculus, corpus callosum and optic radiation exhibiting 14-19% higher fluctuations in voxels containing perpendicular compared to parallel fibres.