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

Layer-specific spatiotemporal dynamics of feedforward and feedback in human visual object perception

Visual object perception is mediated by information flow between regions of the ventral visual stream along feedforward and feedback anatomical connections. However, feedforward and feedback signals during naturalistic vision are rapid and overlapping, complicating their identification and precise functional specification. Here we recorded human layer-specific fMRI responses to naturalistic object images in early visual cortex (EVC) and lateral occipital complex (LOC) to isolate feedforward and feedback information signals spatially by their cortical layer specific termination pattern. We combined these layer-specific fMRI responses with electroencephalography (EEG) responses for the same images to segregate feedforward and feedback signals in both time and space. Feedforward signals emerge early in the middle layers of EVC and LOC, followed by feedback signals in the superficial layer of both regions, and the deep layer of EVC. Comparing the identified dynamics in LOC to a visual deep neural network (DNN), revealed that early feedforward signals in LOC encode medium complexity features, whereas later feedback signals increase the representational format to high complexity features. Together this specifies the spatiotemporal dynamics and functional role of feedforward and feedback information flow mediating visual object perception.

Open-transmit and flexible receiver array for high resolution ultrahigh-field fMRI of the human sensorimotor cortex

In this study, we developed an open-transmit and 24-channel flexible receiver head coil assembly tailored for high-resolution ultrahigh-field functional magnetic resonance imaging (fMRI) of the human somatosensory and motor cortex. Leveraging the increased signal-to-noise ratio (SNR) and spatial resolution of ultrahigh field MRI, we address the technical challenges inherent in fMRI coil design. The open-birdcage transmit coil enhances patient comfort and enables visual task implementation, demonstrating superior performance in transmit efficiency and specific absorption rate distribution compared to conventional coils. Furthermore, the 24-channel flexible receiver head coil offers enhanced SNR and image quality, facilitating sub-millimeter vascular-space-occupancy imaging for precise functional mapping. These advancements provide valuable tools for unraveling the intricacies of somatosensory and motor cortex function. By enriching human brain functional studies, they contribute significantly to our understanding of the mechanisms underlying somatosensory and motor cortex function and may have valuable clinical applications in neurology and neuroscience research.

SORDINO for Silent, Sensitive, Specific, and Artifact-Resisting fMRI in awake behaving mice

Blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) has revolutionized our understanding of the brain activity landscape, bridging circuit neuroscience in animal models with noninvasive brain mapping in humans. This immensely utilized technique, however, faces challenges such as acoustic noise, electromagnetic interference, motion artifacts, magnetic-field inhomogeneity, and limitations in sensitivity and specificity. Here, we introduce Steady-state On-the-Ramp Detection of INduction-decay with Oversampling (SORDINO), a transformative fMRI technique that addresses these challenges by maintaining a constant total gradient amplitude while acquiring data during continuously changing gradient direction. When benchmarked against conventional fMRI on a 9.4T system, SORDINO is silent, sensitive, specific, and resistant to motion and susceptibility artifacts. SORDINO offers superior compatibility with multimodal experiments and carries novel contrast mechanisms distinct from BOLD. It also enables brain-wide activity and connectivity mapping in awake, behaving mice, overcoming stress- and motion-related confounds that are among the most challenging barriers in current animal fMRI studies.

Affective Visual Circuit Dysfunction in Trauma and Stress-Related Disorders

Posttraumatic stress disorder (PTSD) is widely recognized as involving disruption of core neurocircuitry that underlies processing, regulation, and response to threat. In particular, the prefrontal cortex-hippocampal-amygdala circuit is a major contributor to posttraumatic dysfunction. However, the functioning of core threat neurocircuitry is partially dependent on sensorial inputs, and previous research has demonstrated that dense, reciprocal connections exist between threat circuits and the ventral visual stream. Furthermore, emergent evidence suggests that trauma exposure and resultant PTSD symptoms are associated with altered structure and function of the ventral visual stream. In the current review, we discuss evidence that both threat and visual circuitry together are an integral part of PTSD pathogenesis. An overview of the relevance of visual processing to PTSD is discussed in the context of both basic and translational research, highlighting the impact of stress on affective visual circuitry. This review further synthesizes emergent literature to suggest potential timing-dependent effects of traumatic stress on threat and visual circuits that may contribute to PTSD development. We conclude with recommendations for future research to move the field toward a more complete understanding of PTSD neurobiology.

Reduced physiology-induced temporal instability achieved with variable-flip-angle fast low-angle excitation echo-planar technique with multishot echo planar time-resolved imaging

PURPOSE

Echo planar time-resolved imaging (EPTI) is a new imaging approach that addresses the limitations of EPI by providing high-resolution, distortion- and T2/T 2 * $$ {\mathrm{T}}_2^{\ast } $$  blurring-free imaging for functional MRI (fMRI). However, as in all multishot sequences, intershot phase variations induced by physiological processes can introduce temporal instabilities to the reconstructed time-series data. This study aims to reduce these instabilities in multishot EPTI.

THEORY AND METHODS

In conventional multishot EPTI, the time intervals between the shots comprising each slice can introduce intershot phase variations. Here, the fast low-angle excitation echo-planar technique (FLEET), in which all shots of each slice are acquired consecutively with minimal time delays, was combined with a variable flip angle (VFA) technique to improve intershot consistency and maximize signal. A recursive Shinnar-Le Roux RF pulse design algorithm was used to generate pulses for different shots to produce consistent slice profiles and signal intensities across shots. Blipped controlled aliasing in parallel imaging simultaneous multislice was also combined with the proposed VFA-FLEET EPTI to improve temporal resolution and increase spatial coverage.

RESULTS

The temporal stability of VFA-FLEET EPTI was compared with conventional EPTI at 7 T. The results demonstrated that VFA-FLEET can provide spatial-specific increase of temporal stability. We performed high-resolution task-fMRI experiments at 7 T using VFA-FLEET EPTI, and reliable BOLD responses to a visual stimulus were detected.

CONCLUSION

The intershot phase variations induced by physiological processes in multishot EPTI can manifest as specific spatial patterns of physiological noise enhancement and lead to reduced temporal stability. The VFA-FLEET technique can substantially reduce these physiology-induced instabilities in multishot EPTI acquisitions. The proposed method provides sufficient stability and sensitivity for high-resolution fMRI studies.

Enabling brain-wide mapping of layer-specific functional connectivity at 3T via layer-dependent fMRI with draining-vein suppression

Layer-dependent functional magnetic resonance imaging (fMRI) is a promising yet challenging approach for investigating layer-specific functional connectivity (FC). Achieving a brain-wide mapping of layer-specific FC requires several technical advancements, including sub-millimeter spatial resolution, sufficient temporal resolution, functional sensitivity, global brain coverage, and high spatial specificity. Although gradient echo (GE)-based echo planar imaging (EPI) is commonly used for rapid fMRI acquisition, it faces significant challenges due to the draining-vein contamination. In this study, we addressed these limitations by integrating velocity-nulling (VN) gradients into a GE-BOLD fMRI sequence to suppress vascular signals from the vessels with fast-flowing velocity. The extravascular contamination from pial veins was mitigated using a GE-EPI sequence at 3T rather than 7T, combined with phase regression methods. Additionally, we incorporated advanced techniques, including simultaneous multislice (SMS) acceleration and NOise Reduction with DIstribution Corrected principal component analysis (NORDIC PCA) denoising, to improve temporal resolution, spatial coverage, and signal sensitivity. This resulted in a VN fMRI sequence with 0.9-mm isotropic spatial resolution, a repetition time (TR) of 4 seconds, and brain-wide coverage. The VN gradient strength was determined based on results from a button-pressing task. Using resting-state data, we validated layer-specific FC through seed-based analyses, identifying distinct connectivity patterns in the superficial and deep layers of the primary motor cortex (M1), with significant inter-layer differences. Further analyses with a seed in the primary sensory cortex (S1) demonstrated the reliability of the method. Brain-wide layer-dependent FC analyses yielded results consistent with prior literature, reinforcing the efficacy of VN fMRI in resolving layer-specific functional connectivity. Given the widespread availability of 3T scanners, this technical advancement has the potential for significant impact across multiple domains of neuroscience research.

Dynamic layer-specific processing in the prefrontal cortex during working memory

The dorsolateral prefrontal cortex (dlPFC) is reliably engaged in working memory (WM) and comprises different cytoarchitectonic layers, yet their functional role in human WM is unclear. Here, participants completed a delayed-match-to-sample task while undergoing functional magnetic resonance imaging (fMRI) at ultra-high resolution. We examine layer-specific activity to manipulations in WM load and motor response. Superficial layers exhibit a preferential response to WM load during the delay and retrieval periods of a WM task, indicating a lamina-specific activation of the frontoparietal network. Multivariate patterns encoding WM load in the superficial layer dynamically change across the three periods of the task. Last, superficial and deep layers are non-differentially involved in the motor response, challenging earlier findings of a preferential deep layer activation. Taken together, our results provide new insights into the functional laminar circuitry of the dlPFC during WM and support a dynamic account of dlPFC coding.

Precision fMRI and cluster-failure in the individual brain

Advances in neuroimaging acquisition protocols and denoising techniques, along with increasing magnetic field strengths, have dramatically improved the temporal signal-to-noise ratio (tSNR) in functional magnetic resonance imaging (fMRI). This permits spatial resolution with submillimeter voxel sizes and ultrahigh temporal resolution and opens a route toward performing precision fMRI in the brains of individuals. Yet ultrahigh spatial and temporal resolution comes at a cost: it reduces tSNR and, therefore, the sensitivity to the blood oxygen level-dependent (BOLD) effect and other functional contrasts across the brain. Here we investigate the potential of various smoothing filters to improve BOLD sensitivity while preserving the spatial accuracy of activated clusters in single-subject analysis. We introduce adaptive-weight smoothing with optimized metrics (AWSOM), which addresses this challenge extremely well. AWSOM employs a local inference approach that is as sensitive as cluster-corrected inference of data smoothed with large Gaussian kernels, but it preserves spatial details across multiple tSNR levels. This is essential for examining whole-brain fMRI data because tSNR varies across the entire brain, depending on the distance of a brain region from the receiver coil, the type of setup, acquisition protocol, preprocessing, and resolution. We found that cluster correction in single subjects results in inflated family-wise error and false positive rates. AWSOM effectively suppresses false positives while remaining sensitive even to small clusters of activated voxels. Furthermore, it preserves signal integrity, that is, the relative activation strength of significant voxels, making it a valuable asset for a wide range of fMRI applications. Here we demonstrate these features and make AWSOM freely available to the research community for download.

Motion compensated structured low-rank reconstruction for 3D multi-shot EPI

PURPOSE

The 3D multi-shot EPI imaging offers several benefits including higher SNR and high isotropic resolution compared to 2D single shot EPI. However, it suffers from shot-to-shot inconsistencies arising from physiologically induced phase variations and bulk motion. This work proposed a motion compensated structured low-rank (mcSLR) reconstruction method to address both issues for 3D multi-shot EPI.

METHODS

Structured low-rank reconstruction has been successfully used in previous work to deal with inter-shot phase variations for 3D multi-shot EPI imaging. It circumvents the estimation of phase variations by reconstructing an individual image for each phase state which are then sum-of-squares combined, exploiting their linear interdependency encoded in structured low-rank constraints. However, structured low-rank constraints become less effective in the presence of inter-shot motion, which corrupts image magnitude consistency and invalidates the linear relationship between shots. Thus, this work jointly models inter-shot phase variations and motion corruptions by incorporating rigid motion compensation for structured low-rank reconstruction, where motion estimates are obtained in a fully data-driven way without relying on external hardware or imaging navigators.

RESULTS

Simulation and in vivo experiments at 7T have demonstrated that the mcSLR method can effectively reduce image artifacts and improve the robustness of 3D multi-shot EPI, outperforming existing methods which only address inter-shot phase variations or motion, but not both.

CONCLUSION

The proposed mcSLR reconstruction compensates for rigid motion, and thus improves the validity of structured low-rank constraints, resulting in improved robustness of 3D multi-shot EPI to both inter-shot motion and phase variations.

Next-generation MRI scanner designed for ultra-high-resolution human brain imaging at 7 Tesla

To increase granularity in human neuroimaging science, we designed and built a next-generation 7 Tesla magnetic resonance imaging scanner to reach ultra-high resolution by implementing several advances in hardware. To improve spatial encoding and increase the image signal-to-noise ratio, we developed a head-only asymmetric gradient coil (200 mT m-1, 900 T m-1s-1) with an additional third layer of windings. We integrated a 128-channel receiver system with 64- and 96-channel receiver coil arrays to boost signal in the cerebral cortex while reducing g-factor noise to enable higher accelerations. A 16-channel transmit system reduced power deposition and improved image uniformity. The scanner routinely performs functional imaging studies at 0.35-0.45 mm isotropic spatial resolution to reveal cortical layer functional activity, achieves high angular resolution in diffusion imaging and reduces acquisition time for both functional and structural imaging.

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.

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.

Acquisition and processing methods of whole-brain layer-fMRI VASO and BOLD: The Kenshu dataset

Cortical depth-dependent functional magnetic resonance image (fMRI), also known as layer-fMRI, has the potential to capture directional neural information flow of brain computations within and across large-scale cortical brain networks. E.g., layer-fMRI can differentiate feedforward and feedback cortical input in hierarchically organized brain networks. Recent advancements in 3D-EPI sampling approaches and MR contrast generation strategies have allowed proof-of-principle studies showing that layer-fMRI can provide sufficient data quality for capturing laminar changes in functional connectivity. These studies have however not shown how reliable the signal is and how repeatable the respective results are. It is especially unclear whether whole-brain layer-fMRI functional connectivity protocols are widely applicable across common neuroscience-driven analysis approaches. Moreover, there are no established preprocessing fMRI methods that are optimized to work for whole-brain layer-fMRI datasets. In this work, we aimed to serve the field of layer-fMRI and build tools for future routine whole-brain layer-fMRI in application-based neuroscience research. We have developed publicly available sequences, acquisition protocols, and processing pipelines for whole-brain layer-fMRI. These protocols are validated across 60 hours of scanning in nine participants. Specifically, we identified and exploited methodological advancements for maximizing tSNR efficiency and test-retest reliability. We are sharing an extensive multi-modal whole-brain layer-fMRI dataset (20 scan hours of movie-watching in a single participant) for the purpose of benchmarking future method developments: The Kenshu dataset. With this dataset, we are also exemplifying the usefulness of whole brain layer-fMRI for commonly applied analysis approaches in modern cognitive neuroscience fMRI studies. This includes connectivity analyses, representational similarity matrix estimations, general linear model analyses, principal component analysis clustering, etc. We believe that this work paves the road for future routine measurements of directional functional connectivity across the entire brain.

Odor-evoked layer-specific fMRI activities in the awake mouse olfactory bulb

Awake rodent fMRI is increasingly common over the use of anesthesia since it permits behavioral paradigms and does not confound normal brain function or neurovascular coupling. It is well established that adequate acclimation to the loud fMRI environment and head fixation reduces stress in the rodents and allows for whole brain imaging with little contamination from motion. However, it is unknown whether high-resolution fMRI with increased susceptibility to motion and lower sensitivity can measure small, but spatially discrete, activations in awake mice. To examine this, we used contrast-enhanced cerebral blood volume-weighted (CBVw) fMRI in the mouse olfactory bulb for its enhanced sensitivity and neural specificity. We determined that activation patterns in the glomerular layer to four different odors were spatially distinct and were consistent with previously established histological patterns. In addition, odor-evoked laminar activations were greatest in superficial layers that decreased with laminar depth, similar to previous observations. Interestingly, the fMRI response strengths in the granule cell layer were greater in awake mice than our previous anesthetized rat studies, suggesting that feedback neural activities were intact with wakefulness. We finally determined that fMRI signal changes to repeated odor exposure (i.e., olfactory adaptation) attenuated relatively more in the feedback granule cell layer compared to the input glomerular layer, which is consistent with prior observations. We, therefore, conclude that high-resolution CBVw fMRI can measure odor-specific activation patterns and distinguish changes in laminar activity of head and body restrained awake mice.

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.

Differential Laminar Activation Dissociates Encoding and Retrieval in the Human Medial and Lateral Entorhinal Cortex

The hierarchically organized structures of the medial temporal lobe are critically important for episodic memory function. Accumulating evidence suggests dissociable information processing pathways are maintained throughout these structures including in the medial and lateral entorhinal cortex. Cortical layers provide an additional dimension of dissociation as the primary input to the hippocampus derives from layer 2 neurons in the entorhinal cortex, whereas the deeper layers primarily receive output from the hippocampus. Here, novel high-resolution T2-prepared functional MRI methods were successfully used to mitigate susceptibility artifacts typically affecting MRI signals in this region providing uniform sensitivity across the medial and lateral entorhinal cortex. During the performance of a memory task, healthy human subjects (age 25-33 years, mean age 28.2 ± 3.3 years, 4 female) showed differential functional activation in the superficial and deep layers of the entorhinal cortex associated with task-related encoding and retrieval conditions, respectively. The methods provided here offer an approach to probe layer-specific activation in normal cognition and conditions contributing to memory impairment.SIGNIFICANCE STATEMENT This study provides new evidence for differential neuronal activation in the superficial versus deep layers of the entorhinal cortex associated with encoding and retrieval memory processes, respectively, in cognitively normal adults. The study further shows that this dissociation can be observed in both the medial and the lateral entorhinal cortex. The study was achieved by using a novel functional MRI method allowing us to measure robust functional MRI signals in both the medial and lateral entorhinal cortex that was not possible in previous studies. The methodology established here in healthy human subjects lays a solid foundation for subsequent studies investigating layer-specific and region-specific changes in the entorhinal cortex associated with memory impairment in various conditions such as Alzheimer's disease.

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.

Improving robustness of 3D multi-shot EPI by structured low-rank reconstruction of segmented CAIPI sampling for fMRI at 7T

Three-dimensional (3D) encoding methods are increasingly being explored as alternatives to two-dimensional (2D) multi-slice acquisitions in fMRI, particularly in cases where high isotropic resolution is needed. 3D multi-shot EPI acquisition, as the workhorse of 3D fMRI imaging, is susceptible to physiological fluctuations which can induce inter-shot phase variations, and thus reducing the achievable tSNR, negating some of the benefit of 3D encoding. This issue can be particularly problematic at ultra-high fields like 7T, which have more severe off-resonance effects. In this work, we aim to improve the temporal stability of 3D multi-shot EPI at 7T by improving its robustness to inter-shot phase variations. We presented a 3D segmented CAIPI sampling trajectory ("seg-CAIPI") and an improved reconstruction method based on Hankel structured low-rank matrix recovery. Simulation and in-vivo results demonstrate that the combination of the seg-CAIPI sampling scheme and the proposed structured low-rank reconstruction is a promising way to effectively reduce the unwanted temporal variance induced by inter-shot physiological fluctuations, and thus improve the robustness of 3D multi-shot EPI for fMRI.

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.

Macromolecular tissue volume mapping of lateral geniculate nucleus subdivisions in living human brains

The lateral geniculate nucleus (LGN) is a key thalamic nucleus in the visual system, which has an important function in relaying retinal visual input to the visual cortex. The human LGN is composed mainly of magnocellular (M) and parvocellular (P) subdivisions, each of which has different stimulus selectivity in neural response properties. Previous studies have discussed the potential relationship between LGN subdivisions and visual disorders based on psychophysical data on specific types of visual stimuli. However, these relationships remain speculative because non-invasive measurements of these subdivisions are difficult due to the small size of the LGN. Here we propose a method to identify these subdivisions by combining two structural MR measures: high-resolution proton-density weighted images and macromolecular tissue volume (MTV) maps. We defined the M and P subdivisions based on MTV fraction data and tested the validity of the definition by (1) comparing the data with that from human histological studies, (2) comparing the data with functional magnetic resonance imaging measurements on stimulus selectivity, and (3) analyzing the test-retest reliability. The findings demonstrated that the spatial organization of the M and P subdivisions was consistent across subjects and in line with LGN subdivisions observed in human histological data. Moreover, the difference in stimulus selectivity between the subdivisions identified using MTV was consistent with previous physiology literature. The definition of the subdivisions based on MTV was shown to be robust over measurements taken on different days. These results suggest that MTV mapping is a promising approach for evaluating the tissue properties of LGN subdivisions in living humans. This method potentially will enable neuroscientific and clinical hypotheses about the human LGN subdivisions to be tested.

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.

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.

On-scalp magnetocorticography with optically pumped magnetometers: Simulated performance in resolving simultaneous sources

Currently, the gold standard for high-resolution mapping of cortical electrophysiological activity is invasive electrocorticography (ECoG), a procedure that carries with it the risk of serious morbidity and mortality. Due to these risks, the use of ECoG is largely limited to pre-surgical mapping in intractable epilepsy. Nevertheless, many seminal studies in neuroscience have utilized ECoG to explore domains such as visual perception, attention, auditory processing, and sensorimotor behavior. Studies such as these, occurring in patients with epilepsy rather than healthy controls, may lack generalizability, and are limited by the placement of the electrode arrays over the presumed seizure focus. This manuscript explores the use of optically pumped magnetometers (OPMs) to create a non-invasive alternative to ECoG, which we refer to as magnetocorticography. Because prior ECoG studies reveal that most cognitive processes are driven by multiple, simultaneous independent neuronal assemblies, we characterize the ability of a theoretical 56-channel dense OPM array to resolve simultaneous independent sources, and compare it to currently available SQUID devices, as well as OPM arrays with inter-sensor spacings more typical of other systems in development. Our evaluation of this theoretical system assesses many potential sources of error, including errors of sensor calibration and position. In addition, we investigate the influence of geometrical and anatomical factors on array performance. Our simulations reveal the potential of high-density, on-scalp OPM MEG devices to localize electrophysiological brain responses at unprecedented resolution for a non-invasive device.

Automated Brain Volumetry in Patients With Memory Impairment: Comparison of Conventional and Ultrafast 3D T1-Weighted MRI Sequences Using Two Software Packages

BACKGROUND. Isotropic 3D T1-weighted imaging has long acquisition times, potentially leading to motion artifact and altered brain volume measurements. Acquisition times may be greatly shortened using an isotropic ultrafast 3D echo-planar imaging (EPI) T1-weighted sequence. OBJECTIVE. The purpose of this article was to compare automated brain volume measurements between conventional 3D T1-weighted imaging and ultrafast 3D EPI T1-weighted imaging. METHODS. This retrospective study included 36 patients (25 women, 11 men; mean age, 68.4 years) with memory impairment who underwent 3-T brain MRI. Examinations included both conventional 3D T1-weighted imaging using inversion recovery gradient-recalled echo sequence (section thickness, 1.0 mm; acquisition time, 3 minutes 4 seconds) and, in patients exhibiting motion, an isotropic ultrafast 3D EPI T1-weighted sequence (section thickness, 1.2 mm; acquisition time, 30 seconds). The 36-patient sample excluded five patients in whom severe motion artifact rendered the conventional sequence of insufficient quality for volume measurements. Automated brain volumetry was performed using NeuroQuant (version 3.0, CorTechs Laboratories) and FreeSurfer (version 7.1.1, Harvard University) software. Volume measurements were compared between sequences for nine regions in each hemisphere. RESULTS. Volumes showed substantial to almost perfect agreement between the two sequences for most regions bilaterally. However, most regions showed significant mean differences between sequences, and Bland-Altman analyses showed consistent systematic biases and wide limits of agreement (LOA). For example, for the left hemisphere using NeuroQuant, volume was significantly greater for the ultrafast sequence in four regions and significantly greater for the conventional sequence in three regions, whereas standardized effect size between sequences was moderate for four regions and large for one region. Using NeuroQuant, mean bias (ultrafast minus conventional) and 95% LOA were greatest in cortical gray matter bilaterally (-50.61 cm3 [-56.27 cm3, -44.94 cm3] for the left hemisphere; -50.02 cm3 [-54.88 cm3, -45.16 cm3] for the right hemisphere). The variation between the two sequences was observed in subset analyses of 16 patients with and 20 patients without Alzheimer disease. CONCLUSION. Brain volume measurements show significant differences and systematic biases between the conventional and ultrafast sequences. CLINICAL IMPACT. In patients in whom severe motion artifact precludes use of the conventional sequence, the ultrafast sequence may be useful to enable brain volume-try. However, the current conventional 3D T1-weighted sequence remains preferred in patients who can tolerate the standard examination.

Neuroimaging at 7 Tesla: a pictorial narrative review

Neuroimaging using the 7-Tesla (7T) human magnetic resonance (MR) system is rapidly gaining popularity after being approved for clinical use in the European Union and the USA. This trend is the same for functional MR imaging (MRI). The primary advantages of 7T over lower magnetic fields are its higher signal-to-noise and contrast-to-noise ratios, which provide high-resolution acquisitions and better contrast, making it easier to detect lesions and structural changes in brain disorders. Another advantage is the capability to measure a greater number of neurochemicals by virtue of the increased spectral resolution. Many structural and functional studies using 7T have been conducted to visualize details in the white matter and layers of the cortex and hippocampus, the subnucleus or regions of the putamen, the globus pallidus, thalamus and substantia nigra, and in small structures, such as the subthalamic nucleus, habenula, perforating arteries, and the perivascular space, that are difficult to observe at lower magnetic field strengths. The target disorders for 7T neuroimaging range from tumoral diseases to vascular, neurodegenerative, and psychiatric disorders, including Alzheimer's disease, Parkinson's disease, multiple sclerosis, epilepsy, major depressive disorder, and schizophrenia. MR spectroscopy has also been used for research because of its increased chemical shift that separates overlapping peaks and resolves neurochemicals more effectively at 7T than a lower magnetic field. This paper presents a narrative review of these topics and an illustrative presentation of images obtained at 7T. We expect 7T neuroimaging to provide a new imaging biomarker of various brain disorders.

Cortical layer-specific differences in stimulus selectivity revealed with high-field fMRI and single-vessel resolution optical imaging of the primary visual cortex

The mammalian neocortex exhibits a stereotypical laminar organization, with feedforward inputs arriving primarily into layer 4, local computations shaping response selectivity in layers 2/3, and outputs to other brain areas emanating via layers 2/3, 5 and 6. It cannot be assumed a priori that these signatures of laminar differences in neuronal circuitry are reflected in hemodynamic signals that form the basis of functional magnetic resonance imaging (fMRI). Indeed, optical imaging of single-vessel functional responses has highlighted the potential limits of using vascular signals as surrogates for mapping the selectivity of neural responses. Therefore, before fMRI can be employed as an effective tool for studying critical aspects of laminar processing, validation with single-vessel resolution is needed. The primary visual cortex (V1) in cats, with its precise neuronal functional micro-architecture, offers an ideal model system to examine laminar differences in stimulus selectivity across imaging modalities. Here we used cerebral blood volume weighted (wCBV) fMRI to examine if layer-specific orientation-selective responses could be detected in cat V1. We found orientation preference maps organized tangential to the cortical surface that typically extended across depth in a columnar fashion. We then examined arterial dilation and blood velocity responses to identical visual stimuli by using two- and three- photon optical imaging at single-vessel resolution-which provides a measure of the hemodynamic signals with the highest spatial resolution. Both fMRI and optical imaging revealed a consistent laminar response pattern in which orientation selectivity in cortical layer 4 was significantly lower compared to layer 2/3. This systematic change in selectivity across cortical layers has a clear underpinning in neural circuitry, particularly when comparing layer 4 to other cortical layers.

Comparing BOLD and VASO-CBV population receptive field estimates in human visual cortex

Vascular Space Occupancy (VASO) is an alternative fMRI approach based on changes in Cerebral Blood Volume (CBV). VASO-CBV fMRI can provide higher spatial specificity than the blood oxygenation level-dependent (BOLD) method because the CBV response is thought to be limited to smaller vessels. To investigate how this technique compares to BOLD fMRI for cognitive neuroscience applications, we compared population receptive field (pRF) mapping estimates between BOLD and VASO-CBV. We hypothesized that VASO-CBV would elicit distinct pRF properties compared to BOLD. Specifically, since pRF size estimates also depend on vascular sources, we hypothesized that reduced vascular blurring might yield narrower pRFs for VASO-CBV measurements. We used a VASO sequence with a double readout 3D EPI sequence at 7T to simultaneously measure VASO-CBV and BOLD responses in the visual cortex while participants viewed conventional pRF mapping stimuli. Both VASO-CBV and BOLD images show similar eccentricity and polar angle maps across all participants. Compared to BOLD-based measurements, VASO-CBV yielded lower tSNR and variance explained. The pRF size changed with eccentricity similarly for VASO-CBV and BOLD, and the pRF size estimates were similar for VASO-CBV and BOLD, even when we equate variance explained between VASO-CBV and BOLD. This result suggests that the vascular component of the pRF size is not dominating in either VASO-CBV or BOLD.

Using high spatial resolution fMRI to understand representation in the auditory network

Following rapid methodological advances, ultra-high field (UHF) functional and anatomical magnetic resonance imaging (MRI) has been repeatedly and successfully used for the investigation of the human auditory system in recent years. Here, we review this work and argue that UHF MRI is uniquely suited to shed light on how sounds are represented throughout the network of auditory brain regions. That is, the provided gain in spatial resolution at UHF can be used to study the functional role of the small subcortical auditory processing stages and details of cortical processing. Further, by combining high spatial resolution with the versatility of MRI contrasts, UHF MRI has the potential to localize the primary auditory cortex in individual hemispheres. This is a prerequisite to study how sound representation in higher-level auditory cortex evolves from that in early (primary) auditory cortex. Finally, the access to independent signals across auditory cortical depths, as afforded by UHF, may reveal the computations that underlie the emergence of an abstract, categorical sound representation based on low-level acoustic feature processing. Efforts on these research topics are underway. Here we discuss promises as well as challenges that come with studying these research questions using UHF MRI, and provide a future outlook.

Advances in spiral fMRI: A high-resolution study with single-shot acquisition

Spiral fMRI has been put forward as a viable alternative to rectilinear echo-planar imaging, in particular due to its enhanced average k-space speed and thus high acquisition efficiency. This renders spirals attractive for contemporary fMRI applications that require high spatiotemporal resolution, such as laminar or columnar fMRI. However, in practice, spiral fMRI is typically hampered by its reduced robustness and ensuing blurring artifacts, which arise from imperfections in both static and dynamic magnetic fields. Recently, these limitations have been overcome by the concerted application of an expanded signal model that accounts for such field imperfections, and its inversion by iterative image reconstruction. In the challenging ultra-high field environment of 7 Tesla, where field inhomogeneity effects are aggravated, both multi-shot and single-shot 2D spiral imaging at sub-millimeter resolution was demonstrated with high depiction quality and anatomical congruency. In this work, we further these advances towards a time series application of spiral readouts, namely, single-shot spiral BOLD fMRI at 0.8 mm in-plane resolution. We demonstrate that high-resolution spiral fMRI at 7 T is not only feasible, but delivers both excellent image quality, BOLD sensitivity, and spatial specificity of the activation maps, with little artifactual blurring. Furthermore, we show the versatility of the approach with a combined in/out spiral readout at a more typical resolution (1.5 mm), where the high acquisition efficiency allows to acquire two images per shot for improved sensitivity by echo combination.

Simulating Local Deformations in the Human Cortex Due to Blood Flow-Induced Changes in Mechanical Tissue Properties: Impact on Functional Magnetic Resonance Imaging

Investigating human brain tissue is challenging due to the complexity and the manifold interactions between structures across different scales. Increasing evidence suggests that brain function and microstructural features including biomechanical features are related. More importantly, the relationship between tissue mechanics and its influence on brain imaging results remains poorly understood. As an important example, the study of the brain tissue response to blood flow could have important theoretical and experimental consequences for functional magnetic resonance imaging (fMRI) at high spatial resolutions. Computational simulations, using realistic mechanical models can predict and characterize the brain tissue behavior and give us insights into the consequent potential biases or limitations of in vivo, high-resolution fMRI. In this manuscript, we used a two dimensional biomechanical simulation of an exemplary human gyrus to investigate the relationship between mechanical tissue properties and the respective changes induced by focal blood flow changes. The model is based on the changes in the brain’s stiffness and volume due to the vasodilation evoked by neural activity. Modeling an exemplary gyrus from a brain atlas we assessed the influence of different potential mechanisms: (i) a local increase in tissue stiffness (at the level of a single anatomical layer), (ii) an increase in local volume, and (iii) a combination of both effects. Our simulation results showed considerable tissue displacement because of these temporary changes in mechanical properties. We found that the local volume increase causes more deformation and consequently higher displacement of the gyrus. These displacements introduced considerable artifacts in our simulated fMRI measurements. Our results underline the necessity to consider and characterize the tissue displacement which could be responsible for fMRI artifacts.

Linking cortical circuit models to human cognition with laminar fMRI

Laboratory animal research has provided significant knowledge into the function of cortical circuits at the laminar level, which has yet to be fully leveraged towards insights about human brain function on a similar spatiotemporal scale. The use of functional magnetic resonance imaging (fMRI) in conjunction with neural models provides new opportunities to gain important insights from current knowledge. During the last five years, human studies have demonstrated the value of high-resolution fMRI to study laminar-specific activity in the human brain. This is mostly performed at ultra-high-field strengths (≥ 7 T) and is known as laminar fMRI. Advancements in laminar fMRI are beginning to open new possibilities for studying questions in basic cognitive neuroscience. In this paper, we first review recent methodological advances in laminar fMRI and describe recent human laminar fMRI studies. Then, we discuss how the use of laminar fMRI can help bridge the gap between cortical circuit models and human cognition.

Age-related changes in cerebrovascular health and their effects on neural function and cognition: A comprehensive review

The process of aging includes changes in cellular biology that affect local interactions between cells and their environments and eventually propagate to systemic levels. In the brain, where neurons critically depend on an efficient and dynamic supply of oxygen and glucose, age-related changes in the complex interaction between the brain parenchyma and the cerebrovasculature have effects on health and functioning that negatively impact cognition and play a role in pathology. Thus, cerebrovascular health is considered one of the main mechanisms by which a healthy lifestyle, such as habitual cardiorespiratory exercise and a healthful diet, could lead to improved cognitive outcomes with aging. This review aims at detailing how the physiology of the cerebral vascular system changes with age and how these changes lead to differential trajectories of cognitive maintenance or decline. This provides a framework for generating specific mechanistic hypotheses about the efficacy of proposed interventions and lifestyle covariates that contribute to enhanced cognitive well-being. Finally, we discuss the methodological implications of age-related changes in the cerebral vasculature for human cognitive neuroscience research and propose directions for future experiments aimed at investigating age-related changes in the relationship between physiology and cognitive mechanisms.

ULTRAHIGH FIELD and ULTRAHIGH RESOLUTION fMRI

Functional magnetic resonance imaging (fMRI) has become one of the most powerful tools for investigating the human brain. Ultrahigh magnetic field (UHF) of 7 Tesla has played a critical role in enabling higher resolution and more accurate (relative to the neuronal activity) functional maps. However, even with these gains, the fMRI approach is challenged relative to the spatial scale over which brain function is organized. Therefore, going forward, significant advances in fMRI are still needed. Such advances will predominantly come from magnetic fields significantly higher than 7 Tesla, which is the most commonly used UHF platform today, and additional technologies that will include developments in pulse sequences, image reconstruction, noise suppression, and image analysis in order to further enhance and augment the gains than can be realized by going to higher magnetic fields.

A plug-and-play, lightweight, single-axis gradient insert design for increasing spatiotemporal resolution in echo planar imaging-based brain imaging

The goal of this study was to introduce and evaluate the performance of a lightweight, high-performance, single-axis (z-axis) gradient insert design primarily intended for high-resolution functional magnetic resonance imaging, and aimed at providing both ease of use and a boost in spatiotemporal resolution. The optimal winding positions of the coil were obtained using a genetic algorithm with a cost function that balanced gradient performance (minimum 0.30 mT/m/A) and field linearity (≥16 cm linear region). These parameters were verified using field distribution measurements by B₀ -mapping. The correction of geometrical distortions was performed using theoretical field distribution of the coil. Simulations and measurements were performed to investigate the echo planar imaging echo-spacing reduction due to the improved gradient performance. The resulting coil featured a 16-cm linear region, a weight of 45 kg, an installation time of 15 min, and a maximum gradient strength and slew rate of 200 mT/m and 1300 T/m/s, respectively, when paired with a commercially available gradient amplifier (940 V/630 A). The field distribution measurements matched the theoretically expected field. By utilizing the theoretical field distribution, geometrical distortions were corrected to within 6% of the whole-body gradient reference image in the target region. Compared with a whole-body gradient set, a maximum reduction in echo-spacing of a factor of 2.3 was found, translating to a 344 μs echo-spacing, for a field of view of 192 mm, a receiver bandwidth of 920 kHz and a gradient amplitude of 112 mT/m. We present a lightweight, single-axis gradient insert design that can provide high gradient performance and an increase in spatiotemporal resolution with correctable geometrical distortions while also offering a short installation time of less than 15 min and minimal system modifications.

Laminar fMRI using T2-prepared multi-echo FLASH

Functional magnetic resonance imaging (fMRI) using blood oxygenation level dependent (BOLD) contrast at a sub-millimeter scale is a promising technique to probe neural activity at the level of cortical layers. While gradient echo (GRE) BOLD sequences exhibit the highest sensitivity, their signal is confounded by unspecific extravascular (EV) and intravascular (IV) effects of large intracortical ascending veins and pial veins leading to a downstream blurring effect of local signal changes. In contrast, spin echo (SE) fMRI promises higher specificity towards signal changes near the microvascular compartment. However, the T₂-weighted signal is typically sampled with a gradient echo readout imposing additional T₂'-weighting. In this work, we used a T₂-prepared (T₂-prep) sequence with short GRE readouts to investigate its capability to acquire laminar fMRI data during a visual task in humans at 7 T. By varying the T₂-prep echo time (TE_prep) and acquiring multiple gradient echoes (TE_GRE) per excitation, we studied the specificity of the sequence and the influence of possible confounding contributions to the shape of laminar fMRI profiles. By fitting and extrapolating the multi-echo GRE data to a TE_GRE = 0 ms condition, we show for the first time laminar profiles free of T₂'-pollution, confined to gray matter. This finding is independent of TE_prep, except for the shortest one (31 ms) where hints of a remaining intravascular component can be seen. For TE_GRE > 0 ms a prominent peak at the pial surface is observed that increases with longer TE_GRE and dominates the shape of the profiles independent of the amount of T₂-weighting. Simulations show that the peak at the pial surface is a result of static EV dephasing around pial vessels in CSF visible in GM due to partial voluming. Additionally, another, weaker, static dephasing effect is observed throughout all layers of the cortex, which is particularly obvious in the data with shortest T₂-prep echo time. Our simulations show that this cannot be explained by intravascular dephasing but that it is likely caused by extravascular effects of the intracortical and pial veins. We conclude that even for TE_GRE as short as 2.3 ms, the T₂'-weighting added to the T₂-weighting is enough to dramatically affect the laminar specificity of the BOLD signal change. However, the bulk of this corruption stems from CSF partial volume effects which can in principle be addressed by increasing the spatial resolution of the acquisition.

Quantitative relations between BOLD responses, cortical energetics and impulse firing across cortical depth

The blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal arises as a consequence of changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen ( CMR O 2 ) that in turn are modulated by changes in neural activity. Recent advances in imaging have achieved sub-millimetre resolution and allowed investigation of the BOLD response as a function of cortical depth. Here, we adapt our previous theory relating the BOLD signal to neural activity to produce a quantitative model that incorporates venous blood draining between cortical layers. The adjustable inputs to the model are the neural activity and a parameter governing this blood draining. A three-layer version for transient neural inputs and a multi-layer version for constant or tonic neural inputs are able to account for a variety of experimental results, including negative BOLD signals.

Assessing the effective connectivity of premotor areas during real vs imagined grasping: a DCM-PEB approach

The parieto-frontal circuit underlying grasping, which requires the serial involvement of the anterior intraparietal area (aIPs) and the ventral premotor cortex (PMv), has been recently extended enlightening the role of the dorsal premotor cortex (PMd). The supplementary motor area (SMA) has been also suggested to encode grip force for grasping actions; furthermore, both PMd and SMA are known to play a crucial role in motor imagery. Here, we aimed at assessing the dynamic couplings between left aIPs, PMv, PMd, SMA and primary motor cortex (M1) by comparing executed and imagined right-hand grasping, using Dynamic Causal Modelling (DCM) and Parametrical Empirical Bayes (PEB) analyses. 24 subjects underwent an fMRI exam (3T) during which they were asked to perform or imagine a grasping movement visually cued by photographs of commonly used objects. We tested whether the two conditions a) exert a modulatory effect on both forward and feedback couplings among our areas of interest, and b) differ in terms of strength and sign of these parameters. Results of the real condition confirmed the serial involvement of aIPs, PMv and M1. PMv also exerted a positive influence on PMd and SMA, but received an inhibitory feedback only from PMd. Our results suggest that a general motor program for grasping is planned by the aIPs-PMv circuit; then, PMd and SMA encode high-level features of the movement. During imagery, the connection strength from aIPs to PMv was weaker and the information flow stopped in PMv; thus, a less complex motor program was planned. Moreover, results suggest that SMA and PMd cooperate to prevent motor execution. In conclusion, the comparison between execution and imagery reveals that during grasping premotor areas dynamically interplay in different ways, depending on task demands.

An in-vivo study of BOLD laminar responses as a function of echo time and static magnetic field strength

Layer specific functional MRI requires high spatial resolution data. To compensate the associated poor signal to noise ratio it is common to integrate the signal from voxels at a given cortical depth. If the region is sufficiently large then physiological noise will be the dominant noise source. In this work, activation profiles in response to the same visual stimulus are compared at 1.5 T, 3 T and 7 T using a multi-echo, gradient echo (GE) FLASH sequence, with a 0.75 mm isotropic voxel size and the cortical integration approach. The results show that after integrating over a cortical volume of 40, 60 and 100 mm³ (at 7 T, 3 T, and 1.5 T, respectively), the signal is in the physiological noise dominated regime. The activation profiles obtained are similar for equivalent echo times. BOLD-like noise is found to be the dominant source of physiological noise. Consequently, the functional contrast to noise ratio is not strongly echo-time or field-strength dependent. We conclude that laminar GE-BOLD fMRI at lower field strengths is feasible but that larger patches of cortex will need to be examined, and that the acquisition efficiency is reduced.

Comparing hand movement rate dependence of cerebral blood volume and BOLD responses at 7T

Functional magnetic resonance imaging (fMRI) based on the Blood Oxygenation Level Dependent (BOLD) contrast takes advantage of the coupling between neuronal activity and the hemodynamics to allow a non-invasive localisation of the neuronal activity. In general, fMRI experiments assume a linear relationship between neuronal activation and the observed hemodynamics. However, the relationship between BOLD responses, neuronal activity, and behaviour are often nonlinear. In addition, the nonlinearity between BOLD responses and behaviour may be related to neuronal process rather than a neurovascular uncoupling. Further, part of the nonlinearity may be driven by vascular nonlinearity effects in particular from large vessel contributions. fMRI based on cerebral blood volume (CBV), promises a higher microvascular specificity, potentially without vascular nonlinearity effects and reduced contamination of the large draining vessels compared to BOLD. In this study, we aimed to investigate differences in BOLD and VASO-CBV signal changes during a hand movement task over a broad range of movement rates. We used a double readout 3D-EPI sequence at 7T to simultaneously measure VASO-CBV and BOLD responses in the sensorimotor cortex. The measured BOLD and VASO-CBV responses increased very similarly in a nonlinear fashion, plateauing for movement rates larger than 1 Hz. Our findings show a tight relationship between BOLD and VASO-CBV responses, indicating that the overall interplay of CBV and BOLD responses are similar for the assessed range of movement rates. These results suggest that the observed nonlinearity of neuronal origin is already present in VASO-CBV measurements, and consequently shows relatively unchanged BOLD responses.

Forging a path to mesoscopic imaging success with ultra-high field functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) studies with ultra-high field (UHF, 7+ Tesla) technology enable the acquisition of high-resolution images. In this work, we discuss recent achievements in UHF fMRI at the mesoscopic scale, on the order of cortical columns and layers, and examine approaches to addressing common challenges. As researchers push to smaller and smaller voxel sizes, acquisition and analysis decisions have greater potential to degrade spatial accuracy, and UHF fMRI data must be carefully interpreted. We consider the impact of acquisition decisions on the spatial specificity of the MR signal with a representative dataset with 0.8 mm isotropic resolution. We illustrate the trade-offs in contrast with noise ratio and spatial specificity of different acquisition techniques and show that acquisition blurring can increase the effective voxel size by as much as 50% in some dimensions. We further describe how different sources of degradations to spatial resolution in functional data may be characterized. Finally, we emphasize that progress in UHF fMRI depends not only on scientific discovery and technical advancement, but also on informal discussions and documentation of challenges researchers face and overcome in pursuit of their goals. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

New acquisition techniques and their prospects for the achievable resolution of fMRI

This work reviews recent advances in technologies for functional magnetic resonance imaging (fMRI) of the human brain and highlights the push for higher functional specificity based on increased spatial resolution and specific MR contrasts to reveal previously undetectable functional properties of small-scale cortical structures. We discuss how the combination of MR hardware, advanced acquisition techniques and various MR contrast mechanisms have enabled recent progress in functional neuroimaging. However, these advanced fMRI practices have only been applied to a handful of neuroscience questions to date, with the majority of the neuroscience community still using conventional imaging techniques. We thus discuss upcoming challenges and possibilities for fMRI technology development in human neuroscience. We hope that readers interested in functional brain imaging acquire an understanding of current and novel developments and potential future applications, even if they don't have a background in MR physics or engineering. We summarize the capabilities of standard fMRI acquisition schemes with pointers to relevant literature and comprehensive reviews and introduce more recent developments.

Segmented K-space blipped-controlled aliasing in parallel imaging for high spatiotemporal resolution EPI

PURPOSE

A segmented k-space blipped-controlled aliasing in parallel imaging (skipped-CAIPI) sampling strategy for EPI is proposed, which allows for a flexible choice of EPI factor and phase encode bandwidth independent of the controlled aliasing in parallel imaging (CAIPI) sampling pattern.

THEORY AND METHODS

With previously proposed approaches, exactly two EPI trajectories were possible given a specific CAIPI pattern, either with slice gradient blips (blipped-CAIPI) or following a shot-selective CAIPI approach (higher resolution). Recently, interleaved multi-shot segmentation along shot-selective CAIPI trajectories has been applied for high-resolution anatomical imaging. For more flexibility and a broader range of applications, we propose segmentation along any blipped-CAIPI trajectory. Thus, all EPI factors and phase encode bandwidths available with traditional segmented EPI can be combined with controlled aliasing.

RESULTS

Temporal SNR maps of moderate-to-high-resolution time series acquisitions at varying undersampling factors demonstrate beneficial sampling alternatives to blipped-CAIPI or shot-selective CAIPI. Rapid high-resolution scans furthermore demonstrate SNR-efficient and motion-robust structural imaging with almost arbitrary EPI factor and minimal noise penalty.

CONCLUSION

Skipped-CAIPI sampling increases protocol flexibility for high spatiotemporal resolution EPI. In terms of SNR and efficiency, high-resolution functional or structural scans benefit vastly from a free choice of the CAIPI pattern. Even at moderate resolutions, the independence of sampling pattern, TE, and image matrix size is valuable for optimized functional protocol design. Although demonstrated with 3D-EPI, skipped-CAIPI is also applicable with simultaneous multislice EPI.

Encoding schemes in somatosensation: From micro- to meta-topography

Encoding schemes are systematic large-scale arrangements that convert incoming sensory information into a format required for further information processing. The increased spatial resolution of brain images obtained with ultra-high field magnetic resonance imaging at 7 T (7T-MRI) and above increases the granularity and precision of processing units that mediate the link between neuronal encoding and functional readouts. Here, these new developments are reviewed with a focus on human tactile encoding schemes derived from small-scale processing units (in the order of 0.5-5 mm) that are relevant for theoretical and practical concepts of somatosensory encoding and cortical plasticity. Precisely, we review recent approaches to characterize meso-scale maps, layer units, and cortical fields in the sensorimotor cortex of the living human brain and discuss their impact on theories of perception, motor control, topographic encoding, and cortical plasticity. Finally, we discuss concepts on the integration of small-scale processing units into functional networks that span multiple topographic maps and multiple cortical areas. Novel research areas are highlighted that may help to bridge the gap between cortical microstructure and meta-topographic models on brain anatomy and function.

Ultra-high spatial resolution BOLD fMRI in humans using combined segmented-accelerated VFA-FLEET with a recursive RF pulse design

PURPOSE

To alleviate the spatial encoding limitations of single-shot echo-planar imaging (EPI) by developing multi-shot segmented EPI for ultra-high-resolution functional MRI (fMRI) with reduced ghosting artifacts from subject motion and respiration.

THEORY AND METHODS

Segmented EPI can reduce readout duration and reduce acceleration factors, however, the time elapsed between segment acquisitions (on the order of seconds) can result in intermittent ghosting, limiting its use for fMRI. Here, "FLEET" segment ordering, where segments are looped over before slices, was combined with a variable flip angle progression (VFA-FLEET) to improve inter-segment fidelity and maximize signal for fMRI. Scaling a sinc pulse's flip angle for each segment (VFA-FLEET-Sinc) produced inconsistent slice profiles and ghosting, therefore, a recursive Shinnar-Le Roux (SLR) radiofrequency (RF) pulse design was developed (VFA-FLEET-SLR) to generate unique pulses for every segment that together produce consistent slice profiles and signals.

RESULTS

The temporal stability of VFA-FLEET-SLR was compared against conventional-segmented EPI and VFA-FLEET-Sinc at 3T and 7T. VFA-FLEET-SLR showed reductions in both intermittent and stable ghosting compared to conventional-segmented and VFA-FLEET-Sinc, resulting in improved image quality with a minor trade-off in temporal SNR. Combining VFA-FLEET-SLR with acceleration, we achieved a 0.6-mm isotropic acquisition at 7T, without zoomed imaging or partial Fourier, demonstrating reliable detection of blood oxygenation level-dependent (BOLD) responses to a visual stimulus. To counteract the increased repetition time from segmentation, simultaneous multi-slice VFA-FLEET-SLR was demonstrated using RF-encoded controlled aliasing.

CONCLUSIONS

VFA-FLEET with a recursive RF pulse design supports acquisitions with low levels of artifact and spatial blur, enabling fMRI at previously inaccessible spatial resolutions with a "full-brain" field of view.