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

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.

Determining laminar neuronal activity from BOLD fMRI using a generative model

Laminar fMRI using the BOLD contrast enables the non-invasive investigation of mesoscopic functional circuits in the human brain. However, the laminar neuronal activity is spatiotemporally biased in the observed cortical depth profiles of the BOLD signal. In this study, we propose a generative fMRI signal model, comprehensively covering the relationship between cortical depth-dependent changes in excitatory and inhibitory neuronal activity with the sampling of the BOLD signal with finite voxels. The generative model allowed us to investigate pertinent questions regarding the accuracy of the laminar BOLD signal relative to the neuronal activity, and we found that: a) condition differences in laminar BOLD signals may be more reflective of neuronal activity than single condition BOLD signal depth profiles; b) angular dependence of the BOLD signal induces significant signal variability, which can mask underlying activity profiles; c) even if only three neuronal depths are of interest, more BOLD signal depths should be considered in the analysis. In addition, we recommend that the laminar BOLD data should be displayed using the centroid method to appreciate its spatial distribution in the original resolution. Finally, we showed that Bayesian model inversion of the generative model can improve sensitivity and specificity of assessing depth-dependent neuronal changes both for steady-state and dynamically.

Sub-millimetre resolution laminar fMRI using Arterial Spin Labelling in humans at 7 T

Laminar fMRI at ultra-high magnetic field strength is typically carried out using the Blood Oxygenation Level-Dependent (BOLD) contrast. Despite its unrivalled sensitivity to detecting activation, the BOLD contrast is limited in its spatial specificity due to signals stemming from intra-cortical ascending and pial veins. Alternatively, regional changes in perfusion (i.e., cerebral blood flow through tissue) are colocalised to neuronal activation, which can be non-invasively measured using Arterial Spin Labelling (ASL) MRI. In addition, ASL provides a quantitative marker of neuronal activation in terms of perfusion signal, which is simultaneously acquired along with the BOLD signal. However, ASL for laminar imaging is challenging due to the lower SNR of the perfusion signal and higher RF power deposition i.e., specific absorption rate (SAR) of ASL sequences. In the present study, we present for the first time in humans, isotropic sub-millimetre spatial resolution functional perfusion images using Flow-sensitive Alternating Inversion Recovery (FAIR) ASL with a 3D-EPI readout at 7 T. We show that robust statistical activation maps can be obtained with perfusion-weighting in a single session. We observed the characteristic BOLD amplitude increase towards the superficial laminae, and, in apparent discrepancy, the relative perfusion profile shows a decrease of the amplitude and the absolute perfusion profile a much smaller increase towards the cortical surface. Considering the draining vein effect on the BOLD signal using model-based spatial “convolution”, we show that the empirically measured perfusion and BOLD profiles are, in fact, consistent with each other. This study demonstrates that laminar perfusion fMRI in humans is feasible at 7 T and that caution must be exercised when interpreting BOLD signal laminar profiles as direct representation of the cortical distribution of neuronal activity.

Large-scale functional ultrasound imaging of the spinal cord reveals in-depth spatiotemporal responses of spinal nociceptive circuits in both normal and inflammatory states

Despite a century of research on the physiology/pathophysiology of the spinal cord in chronic pain condition, the properties of the spinal cord were rarely studied at the large-scale level from a neurovascular point of view. This is mostly due to the limited spatial and/or temporal resolution of the available techniques. Functional ultrasound imaging (fUS) is an emerging neuroimaging approach that allows, through the measurement of cerebral blood volume, the study of brain functional connectivity or functional activations with excellent spatial (100 μm) and temporal (1 msec) resolutions and a high sensitivity. The aim of this study was to increase our understanding of the spinal cord physiology through the study of the properties of spinal hemodynamic response to the natural or electrical stimulation of afferent fibers. Using a combination of fUS and ultrasound localization microscopy, the first step of this study was the fine description of the vascular structures in the rat spinal cord. Then, using either natural or electrical stimulations of different categories of afferent fibers (Aβ, Aδ, and C fibers), we could define the characteristics of the typical hemodynamic response of the rat spinal cord experimentally. We showed that the responses are fiber-specific, located ipsilaterally in the dorsal horn, and that they follow the somatotopy of afferent fiber entries in the dorsal horn and that the C-fiber response is an N-methyl-D-aspartate receptor-dependent mechanism. Finally, fUS imaging of the mesoscopic hemodynamic response induced by natural tactile stimulations revealed a potentiated response in inflammatory condition, suggesting an enhanced response to allodynic stimulations.

A self-decoupled 32-channel receive array for human-brain MRI at 10.5 T

PURPOSE

Receive array layout, noise mitigation, and B₀ field strength are crucial contributors to SNR and parallel-imaging performance. Here, we investigate SNR and parallel-imaging gains at 10.5 T compared with 7 T using 32-channel receive arrays at both fields.

METHODS

A self-decoupled 32-channel receive array for human brain imaging at 10.5 T (10.5T-32Rx), consisting of 31 loops and one cloverleaf element, was co-designed and built in tandem with a 16-channel dual-row loop transmitter. Novel receive array design and self-decoupling techniques were implemented. Parallel imaging performance, in terms of SNR and noise amplification (g-factor), of the 10.5T-32Rx was compared with the performance of an industry-standard 32-channel receiver at 7 T (7T-32Rx) through experimental phantom measurements.

RESULTS

Compared with the 7T-32Rx, the 10.5T-32Rx provided 1.46 times the central SNR and 2.08 times the peripheral SNR. Minimum inverse g-factor value of the 10.5T-32Rx (min[1/g] = 0.56) was 51% higher than that of the 7T-32Rx (min[1/g] = 0.37) with R = 4 × 4 2D acceleration, resulting in significantly enhanced parallel-imaging performance at 10.5 T compared with 7 T. The g-factor values of 10.5 T-32 Rx were on par with those of a 64-channel receiver at 7 T (eg, 1.8 vs 1.9, respectively, with R = 4 × 4 axial acceleration).

CONCLUSION

Experimental measurements demonstrated effective self-decoupling of the receive array as well as substantial gains in SNR and parallel-imaging performance at 10.5 T compared with 7 T.

Early fMRI responses to somatosensory and optogenetic stimulation reflect neural information flow

fMRI has revolutionized how neuroscientists investigate human brain functions and networks. To further advance understanding of brain functions, identifying the direction of information flow, such as thalamocortical versus corticothalamic projections, is critical. Because the early hemodynamic response at microvessels near active neurons can be detected by ultrahigh field fMRI, we propose using the onset times of fMRI responses to discern the information flow. This approach was confirmed by observing the ultrahigh spatiotemporal resolution BOLD fMRI responses to bottom-up somatosensory stimulation and top-down optogenetic stimulation of the primary motor cortex in anesthetized mice. Because ultrahigh field MRI is increasingly available, ultrahigh spatiotemporal fMRI will significantly facilitate the investigation of functional circuits in humans.

Blood oxygenation level–dependent (BOLD) functional magnetic resonance imaging (fMRI) has been widely used to localize brain functions. To further advance understanding of brain functions, it is critical to understand the direction of information flow, such as thalamocortical versus corticothalamic projections. For this work, we performed ultrahigh spatiotemporal resolution fMRI at 15.2 T of the mouse somatosensory network during forepaw somatosensory stimulation and optogenetic stimulation of the primary motor cortex (M1). Somatosensory stimulation induced the earliest BOLD response in the ventral posterolateral nucleus (VPL), followed by the primary somatosensory cortex (S1) and then M1 and posterior thalamic nucleus. Optogenetic stimulation of excitatory neurons in M1 induced the earliest BOLD response in M1, followed by S1 and then VPL. Within S1, the middle cortical layers responded to somatosensory stimulation earlier than the upper or lower layers, whereas the upper cortical layers responded earlier than the other two layers to optogenetic stimulation in M1. The order of early BOLD responses was consistent with the canonical understanding of somatosensory network connections and cannot be explained by regional variabilities in the hemodynamic response functions measured using hypercapnic stimulation. Our data demonstrate that early BOLD responses reflect the information flow in the mouse somatosensory network, suggesting that high-field fMRI can be used for systems-level network analyses.

Optimizing Mental Workload Estimation by Detecting Baseline State Using Vector Phase Analysis Approach

Improper baseline return from the previous task-evoked hemodynamic response (HR) can contribute to a large variation in the subsequent HR, affecting the estimation of mental workload in brain-computer interface systems. In this study, we proposed a method using vector phase analysis to detect the baseline state as being optimal or suboptimal. We hypothesize that selecting neuronal-related HR as observed in the optimal-baseline blocks can lead to an improvement in estimating mental workload. Oxygenated and deoxygenated hemoglobin concentration changes were integrated as parts of the vector phase. The proposed method was applied to a block-design functional near-infrared spectroscopy dataset (total blocks = 1384), measured on 24 subjects performing multiple difficulty levels of mental arithmetic task. Significant differences in hemodynamic signal change were observed between the optimal- and suboptimal-baseline blocks detected using the proposed method. This supports the effectiveness of the proposed method in detecting baseline state for better estimation of mental workload. The results further highlight the need of customized recovery duration. In short, the proposed method offers a practical approach to detect task-evoked signals, without the need of extra probes.

Investigating the effects of healthy cognitive aging on brain functional connectivity using 4.7 T resting-state functional magnetic resonance imaging

Functional changes in the aging human brain have been previously reported using functional magnetic resonance imaging (fMRI). Earlier resting-state fMRI studies revealed an age-associated weakening of intra-system functional connectivity (FC) and age-associated strengthening of inter-system FC. However, the majority of such FC studies did not investigate the relationship between age and network amplitude, without which correlation-based measures of FC can be challenging to interpret. Consequently, the main aim of this study was to investigate how three primary measures of resting-state fMRI signal-network amplitude, network topography, and inter-network FC-are affected by healthy cognitive aging. We acquired resting-state fMRI data on a 4.7 T scanner for 105 healthy participants representing the entire adult lifespan (18-85 years of age). To study age differences in network structure, we combined ICA-based network decomposition with sparse graphical models. Older adults displayed lower blood-oxygen-level-dependent (BOLD) signal amplitude in all functional systems, with sensorimotor networks showing the largest age differences. Our age comparisons of network topography and inter-network FC demonstrated a substantial amount of age invariance in the brain's functional architecture. Despite architecture similarities, old adults displayed a loss of communication efficiency in our inter-network FC comparisons, driven primarily by the FC reduction in frontal and parietal association cortices. Together, our results provide a comprehensive overview of age effects on fMRI-based FC.

Cross-species neuroscience: closing the explanatory gap

Neuroscience has seen substantial development in non-invasive methods available for investigating the living human brain. However, these tools are limited to coarse macroscopic measures of neural activity that aggregate the diverse responses of thousands of cells. To access neural activity at the cellular and circuit level, researchers instead rely on invasive recordings in animals. Recent advances in invasive methods now permit large-scale recording and circuit-level manipulations with exquisite spatio-temporal precision. Yet, there has been limited progress in relating these microcircuit measures to complex cognition and behaviour observed in humans. Contemporary neuroscience thus faces an explanatory gap between macroscopic descriptions of the human brain and microscopic descriptions in animal models. To close the explanatory gap, we propose adopting a cross-species approach. Despite dramatic differences in the size of mammalian brains, this approach is broadly justified by preserved homology. Here, we outline a three-armed approach for effective cross-species investigation that highlights the need to translate different measures of neural activity into a common space. We discuss how a cross-species approach has the potential to transform basic neuroscience while also benefiting neuropsychiatric drug development where clinical translation has, to date, seen minimal success. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Validating Linear Systems Analysis for Laminar fMRI: Temporal Additivity for Stimulus Duration Manipulations

Advancements in ultra-high field (7 T and higher) magnetic resonance imaging (MRI) scanners have made it possible to investigate both the structure and function of the human brain at a sub-millimeter scale. As neuronal feedforward and feedback information arrives in different layers, sub-millimeter functional MRI has the potential to uncover information processing between cortical micro-circuits across cortical depth, i.e. laminar fMRI. For nearly all conventional fMRI analyses, the main assumption is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal adheres to the principles of linear systems theory. For laminar fMRI, however, directional blood pooling across cortical depth stemming from the anatomy of the cortical vasculature, potentially violates these linear system assumptions, thereby complicating analysis and interpretation. Here we assess whether the temporal additivity requirement of linear systems theory holds for laminar fMRI. We measured responses elicited by viewing stimuli presented for different durations and evaluated how well the responses to shorter durations predicted those elicited by longer durations. We find that BOLD response predictions are consistently good predictors for observed responses, across all cortical depths, and in all measured visual field maps (V1, V2, and V3). Our results suggest that the temporal additivity assumption for linear systems theory holds for laminar fMRI. We thus show that the temporal additivity assumption holds across cortical depth for sub-millimeter gradient-echo BOLD fMRI in early visual cortex.

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'.

Hemodynamic response varies across tactile stimuli with different temporal structures

Tactile stimuli can be distinguished based on their temporal features (e.g., duration, local frequency, and number of pulses), which are fundamental for vibrotactile frequency perception. Characterizing how the hemodynamic response changes in shape across experimental conditions is important for designing and interpreting fMRI studies on tactile information processing. In this study, we focused on periodic tactile stimuli with different temporal structures and explored the hemodynamic response function (HRF) induced by these stimuli. We found that HRFs were stimulus‐dependent in tactile‐related brain areas. Continuous stimuli induced a greater area of activation and a stronger and narrower hemodynamic response than intermittent stimuli with the same duration. The magnitude of the HRF increased with increasing stimulus duration. By normalizing the characteristics into topographic matrix, nonlinearity was obvious. These results suggested that stimulation patterns and duration within a cycle may be key characters for distinguishing different stimuli. We conclude that different temporal structures of tactile stimuli induced different HRFs, which are essential for vibrotactile perception and should be considered in fMRI experimental designs and analyses.

Depth-dependent functional MRI responses to chromatic and achromatic stimuli throughout V1 and V2

In the primate visual system, form (shape, location) and color information are processed in separate but interacting pathways. Recent access to high-resolution neuroimaging has facilitated the exploration of the structure of these pathways at the mesoscopic level in the human visual cortex. We used 7T fMRI to observe selective activation of the primary visual cortex to chromatic versus achromatic stimuli in five participants across two scanning sessions. Achromatic checkerboards with low spatial frequency and high temporal frequency targeted the color-insensitive magnocellular pathway. Chromatic checkerboards with higher spatial frequency and low temporal frequency targeted the color-selective parvocellular pathway. This work resulted in three main findings. First, responses driven by chromatic stimuli had a laminar profile biased towards superficial layers of V1, as compared to responses driven by achromatic stimuli. Second, we found stronger preference for chromatic stimuli in parafoveal V1 compared with peripheral V1. Finally, we found alternating, stimulus-selective bands stemming from the V1 border into V2 and V3. Similar alternating patterns have been previously found in both NHP and human extrastriate cortex. Together, our findings confirm the utility of fMRI for revealing details of mesoscopic neural architecture in human cortex.

A Platform for Brain-wide Volumetric Functional Ultrasound Imaging and Analysis of Circuit Dynamics in Awake Mice

Imaging large-scale circuit dynamics is crucial to understanding brain function, but most techniques have a limited depth of field. Here, we describe volumetric functional ultrasound imaging (vfUSI), a platform for brain-wide vfUSI of hemodynamic activity in awake head-fixed mice. We combined a high-frequency 1,024-channel 2D-array transducer with advanced multiplexing and high-performance computing for real-time 3D power Doppler imaging at a high spatiotemporal resolution (220 × 280 × 175 μm³, up to 6 Hz). We developed a standardized software pipeline for registration, segmentation, and temporal analysis in 268 individual brain regions based on the Allen Mouse Common Coordinate Framework. We demonstrated the high sensitivity of vfUSI under multiple experimental conditions, and we successfully imaged stimulus-evoked activity when only a few trials were averaged. We also mapped neural circuits in vivo across the whole brain during optogenetic activation of specific cell types. Moreover, we identified the sequential activation of sensory-motor networks during a grasping water-droplet task.

Higher and deeper: Bringing layer fMRI to association cortex

Recent advances in fMRI have enabled non-invasive measurements of brain function in awake, behaving humans at unprecedented spatial resolutions, allowing us to separate activity in distinct cortical layers. While most layer fMRI studies to date have focused on primary cortices, we argue that the next big steps forward in our understanding of cognition will come from expanding this technology into higher-order association cortex, to characterize depth-dependent activity during increasingly sophisticated mental processes. We outline phenomena and theories ripe for investigation with layer fMRI, including perception and imagery, selective attention, and predictive coding. We discuss practical and theoretical challenges to cognitive applications of layer fMRI, including localizing regions of interest in the face of substantial anatomical heterogeneity across individuals, designing appropriate task paradigms within the confines of acquisition parameters, and generating hypotheses for higher-order brain regions where the laminar circuitry is less well understood. We consider how applying layer fMRI in association cortex may help inform computational models of brain function as well as shed light on consciousness and mental illness, and issue a call to arms to our fellow methodologists and neuroscientists to bring layer fMRI to this next frontier.

A Review of Translational Magnetic Resonance Imaging in Human and Rodent Experimental Models of Small Vessel Disease

Cerebral small vessel disease (SVD) is a major health burden, yet the pathophysiology remains poorly understood with no effective treatment. Since much of SVD develops silently and insidiously, non-invasive neuroimaging such as MRI is fundamental to detecting and understanding SVD in humans. Several relevant SVD rodent models are established for which MRI can monitor in vivo changes over time prior to histological examination. Here, we critically review the MRI methods pertaining to salient rodent models and evaluate synergies with human SVD MRI methods. We found few relevant publications, but argue there is considerable scope for greater use of MRI in rodent models, and opportunities for harmonisation of the rodent-human methods to increase the translational potential of models to understand SVD in humans. We summarise current MR techniques used in SVD research, provide recommendations and examples and highlight practicalities for use of MRI SVD imaging protocols in pre-selected, relevant rodent models.

Modeling of vascular space occupancy and BOLD functional MRI from first principles using real microvascular angiograms

PURPOSE

The vascular space occupancy (VASO) is a functional MRI technique for probing cerebral blood volume changes noninvasively, including during neuronal activation in humans. An important consideration when implementing VASO is the BOLD effect in the signal. Assessing the physical origin of this BOLD contamination and the capabilities of correction methods could improve the quantification of cerebral blood volume changes with VASO.

METHODS

Given the heterogeneity of cerebral microvascular architecture, the vascular geometry within an MRI voxel can influence both BOLD and VASO signals. To investigate this effect, 3D high-resolution images of mouse cerebral vasculature measured with two-photon microscopy were used to model BOLD and VASO signals from first principles using Monte Carlo diffusion of water protons. Quantitative plots of VASO together with intravascular and extravascular BOLD signals as a function of TE at B₀ fields 1.5 T to 14 T were obtained.

RESULTS

The BOLD contamination of the VASO response was on the order of 50% for gradient echo and 5% for spin echo at 7 T and TE = 6 ms and significantly increased with TE and B₀ . Two currently used correction schemes were shown to account for most of this contamination and recover accurate relative signal changes, with optimal correction obtained using TEs as short as possible.

CONCLUSION

These results may provide useful information for optimizing sequence parameters in VASO and BOLD functional MRI, leading the way to a wider application of these techniques in healthy and diseased brain.

Transfer functions linking neural calcium to single voxel functional ultrasound signal

Functional ultrasound imaging (fUS) is an emerging technique that detects changes of cerebral blood volume triggered by brain activation. Here, we investigate the extent to which fUS faithfully reports local neuronal activation by combining fUS and two-photon microscopy (2PM) in a co-registered single voxel brain volume. Using a machine-learning approach, we compute and validate transfer functions between dendritic calcium signals of specific neurons and vascular signals measured at both microscopic (2PM) and mesoscopic (fUS) levels. We find that transfer functions are robust across a wide range of stimulation paradigms and animals, and reveal a second vascular component of neurovascular coupling upon very strong stimulation. We propose that transfer functions can be considered as reliable quantitative reporters to follow neurovascular coupling dynamics.

Neurovascular coupling refers to changes in cerebral blood flow in response to neuronal stimulation, but to what extent this change can report neuronal activation is not known. Here the authors develop transfer functions between neural calcium signals and functional ultrasound changes in blood volume in co-registered single voxel brain volumes.

Neuroanatomical Changes in Leber's Hereditary Optic Neuropathy: Clinical Application of 7T MRI Submillimeter Morphometry

Leber’s hereditary optic neuropathy (LHON) is one of the mitochondrial diseases that causes loss of central vision, progressive impairment and subsequent degeneration of retinal ganglion cells (RGCs). In recent years, diffusion tensor imaging (DTI) studies have revealed structural abnormalities in visual white matter tracts, such as the optic tract, and optic radiation. However, it is still unclear if the disease alters only some parts of the white matter architecture or whether the changes also affect other subcortical areas of the brain. This study aimed to improve our understanding of morphometric changes in subcortical brain areas and their associations with the clinical picture in LHON by the application of a submillimeter surface-based analysis approach to the ultra-high-field 7T magnetic resonance imaging data. To meet these goals, fifteen LHON patients and fifteen age-matched healthy subjects were examined. For all individuals, quantitative analysis of the morphometric results was performed. Furthermore, morphometric characteristics which differentiated the groups were correlated with variables covering selected aspects of the LHON clinical picture. Compared to healthy controls (HC), LHON carriers showed significantly lower volume of both palladiums (left p = 0.023; right p = 0.018), the right accumbens area (p = 0.007) and the optic chiasm (p = 0.014). Additionally, LHON patients have significantly higher volume of both lateral ventricles (left p = 0.034; right p = 0.02), both temporal horns of the lateral ventricles (left p = 0.016; right p = 0.034), 3rd ventricle (p = 0.012) and 4th ventricle (p = 0.002). Correlation between volumetric results and clinical data showed that volume of both right and left lateral ventricles significantly and positively correlated with the duration of the illness (left R = 0.841, p = 0.002; right R = 0.755, p = 0.001) and the age of the LHON participants (left R = 0.656, p = 0.007; right R = 0.691, p = 0.004). The abnormalities in volume of the LHON patients’ subcortical structures indicate that the disease can cause changes not only in the white matter areas constituting visual tracts, but also in the other subcortical brain structures. Furthermore, the correlation between those results and the illness duration suggests that the disease might have a neurodegenerative nature; however, to fully confirm this observation, longitudinal studies should be conducted.

Feedback contribution to surface motion perception in the human early visual cortex

Human visual surface perception has neural correlates in early visual cortex, but the role of feedback during surface segmentation in human early visual cortex remains unknown. Feedback projections preferentially enter superficial and deep anatomical layers, which provides a hypothesis for the cortical depth distribution of fMRI activity related to feedback. Using ultra-high field fMRI, we report a depth distribution of activation in line with feedback during the (illusory) perception of surface motion. Our results fit with a signal re-entering in superficial depths of V1, followed by a feedforward sweep of the re-entered information through V2 and V3. The magnitude and sign of the BOLD response strongly depended on the presence of texture in the background, and was additionally modulated by the presence of illusory motion perception compatible with feedback. In summary, the present study demonstrates the potential of depth-resolved fMRI in tackling biomechanical questions on perception.

Direct comparison of activation maps during galvanic vestibular stimulation: A hybrid H2[15 O] PET-BOLD MRI activation study

Previous unimodal PET and fMRI studies in humans revealed a reproducible vestibular brain activation pattern, but with variations in its weighting and expansiveness. Hybrid studies minimizing methodological variations at baseline conditions are rare and still lacking for task-based designs. Thus, we applied for the first time hybrid 3T PET-MRI scanning (Siemens mMR) in healthy volunteers using galvanic vestibular stimulation (GVS) in healthy volunteers in order to directly compare H₂¹⁵O-PET and BOLD MRI responses. List mode PET acquisition started with the injection of 750 MBq H₂¹⁵O simultaneously to MRI EPI sequences. Group-level statistical parametric maps were generated for GVS vs. rest contrasts of PET, MR-onset (event-related), and MR-block. All contrasts showed a similar bilateral vestibular activation pattern with remarkable proximity of activation foci. Both BOLD contrasts gave more bilateral wide-spread activation clusters than PET; no area showed contradictory signal responses. PET still confirmed the right-hemispheric lateralization of the vestibular system, whereas BOLD-onset revealed only a tendency. The reciprocal inhibitory visual-vestibular interaction concept was confirmed by PET signal decreases in primary and secondary visual cortices, and BOLD-block decreases in secondary visual areas. In conclusion, MRI activation maps contained a mixture of CBF measured using H₂¹⁵O-PET and additional non-CBF effects, and the activation-deactivation pattern of the BOLD-block appears to be more similar to the H₂¹⁵O-PET than the BOLD-onset.

Multivoxel Pattern of Blood Oxygen Level Dependent Activity can be sensitive to stimulus specific fine scale responses

At ultra-high field, fMRI voxels can span the sub-millimeter range, allowing the recording of blood oxygenation level dependent (BOLD) responses at the level of fundamental units of neural computation, such as cortical columns and layers. This sub-millimeter resolution, however, is only nominal in nature as a number of factors limit the spatial acuity of functional voxels. Multivoxel Pattern Analysis (MVPA) may provide a means to detect information at finer spatial scales that may otherwise not be visible at the single voxel level due to limitations in sensitivity and specificity. Here, we evaluate the spatial scale of stimuli specific BOLD responses in multivoxel patterns exploited by linear Support Vector Machine, Linear Discriminant Analysis and Naïve Bayesian classifiers across cortical depths in V1. To this end, we artificially misaligned the testing relative to the training portion of the data in increasing spatial steps, then investigated the breakdown of the classifiers’ performances. A one voxel shift led to a significant decrease in decoding accuracy (p < 0.05) across all cortical depths, indicating that stimulus specific responses in a multivoxel pattern of BOLD activity exploited by multivariate decoders can be as precise as the nominal resolution of single voxels (here 0.8 mm isotropic). Our results further indicate that large draining vessels, prominently residing in proximity of the pial surface, do not, in this case, hinder the ability of MVPA to exploit fine scale patterns of BOLD signals. We argue that tailored analytical approaches can help overcoming limitations in high-resolution fMRI and permit studying the mesoscale organization of the human brain with higher sensitivities.

Linear systems analysis for laminar fMRI: Evaluating BOLD amplitude scaling for luminance contrast manipulations

A fundamental assumption of nearly all functional magnetic resonance imaging (fMRI) analyses is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal can be described as following linear systems theory. With the advent of ultra-high field (7T and higher) MRI scanners, it has become possible to perform sub-millimeter resolution fMRI in humans. A novel and promising application of sub-millimeter fMRI is measuring responses across cortical depth, i.e. laminar imaging. However, the cortical vasculature and associated directional blood pooling towards the pial surface strongly influence the cortical depth-dependent BOLD signal, particularly for gradient-echo BOLD. This directional pooling may potentially affect BOLD linearity across cortical depth. Here we assess whether the amplitude scaling assumption for linear systems theory holds across cortical depth. For this, we use stimuli with different luminance contrasts to elicit different BOLD response amplitudes. We find that BOLD amplitude across cortical depth scales with luminance contrast, and that this scaling is identical across cortical depth. Although nonlinearities may be present for different stimulus configurations and acquisition protocols, our results suggest that the amplitude scaling assumption for linear systems theory across cortical depth holds for luminance contrast manipulations in sub-millimeter laminar BOLD fMRI.

Cerebral oxygenation during locomotion is modulated by respiration

In the brain, increased neural activity is correlated with increases of cerebral blood flow and tissue oxygenation. However, how cerebral oxygen dynamics are controlled in the behaving animal remains unclear. We investigated to what extent cerebral oxygenation varies during locomotion. We measured oxygen levels in the cortex of awake, head-fixed mice during locomotion using polarography, spectroscopy, and two-photon phosphorescence lifetime measurements of oxygen sensors. We find that locomotion significantly and globally increases cerebral oxygenation, specifically in areas involved in locomotion, as well as in the frontal cortex and the olfactory bulb. The oxygenation increase persists when neural activity and functional hyperemia are blocked, occurred both in the tissue and in arteries feeding the brain, and is tightly correlated with respiration rate and the phase of respiration cycle. Thus, breathing rate is a key modulator of cerebral oxygenation and should be monitored during hemodynamic imaging, such as in BOLD fMRI.

Microsleep is associated with brain activity patterns unperturbed by auditory inputs

Microsleeps are brief episodes of arousal level decrease manifested through behavioral signs. Brain activity during microsleep in the presence of external stimulus remains poorly understood. In this study, we sought to understand neural responses to auditory stimulation during microsleep. We gave participants the simple task of listening to audios of different pitches and amplitude modulation frequencies during early afternoon functional MRI scans. We found the following: 1) microsleep was associated with cortical activations in broad motor and sensory regions and deactivations in thalamus, irrespective of auditory stimulation; 2) high and low pitch audios elicited different activity patterns in the auditory cortex during awake but not microsleep state; and 3) during microsleep, spatial activity patterns in broad brain regions were similar regardless of the presence or types of auditory stimulus (i.e., stimulus invariant). These findings show that the brain is highly active during microsleep but the activity patterns across broad regions are unperturbed by auditory inputs.NEW & NOTEWORTHY During deep drowsy states, auditory inputs could induce activations in the auditory cortex, but the activation patterns lose differentiation to high/low pitch stimuli. Instead of random activations, activity patterns across the brain during microsleep appear to be structured and may reflect underlying neurophysiological processes that remain unclear.

Dynamic Contrast Optical Coherence Tomography reveals laminar microvascular hemodynamics in the mouse neocortex in vivo

Studies of flow-metabolism coupling often presume that microvessel architecture is a surrogate for blood flow. To test this assumption, we introduce an in vivo Dynamic Contrast Optical Coherence Tomography (DyC-OCT) method to quantify layer-resolved microvascular blood flow and volume across the full depth of the mouse neocortex, where the angioarchitecture has been previously described. First, we cross-validate average DyC-OCT cortical flow against conventional Doppler OCT flow. Next, with laminar DyC-OCT, we discover that layer 4 consistently exhibits the highest microvascular blood flow, approximately two-fold higher than the outer cortical layers. While flow differences between layers are well-explained by microvascular volume and density, flow differences between subjects are better explained by transit time. Finally, from layer-resolved tracer enhancement, we also infer that microvascular hematocrit increases in deep cortical layers, consistent with predictions of plasma skimming. Altogether, our results show that while the cortical blood supply derives mainly from the pial surface, laminar hemodynamics ensure that the energetic needs of individual cortical layers are met. The laminar trends reported here provide data that links predictions based on the cortical angioarchitecture to cerebrovascular physiology in vivo.

Vascular and neural basis of the BOLD signal

Neural activity in the brain is usually coupled to increases in local cerebral blood flow, leading to the increase in oxygenation that generates the BOLD fMRI signal. Recent work has begun to elucidate the vascular and neural mechanisms underlying the BOLD signal. The dilatory response is distributed throughout the vascular network. Arteries actively dilate within a second following neural activity increases, while venous distensions are passive and have a time course that last tens of seconds. Vasodilation, and thus local blood flow, is controlled by the activity of both neurons and astrocytes via multiple different pathways. The relationship between sensory-driven neural activity and the vascular dynamics in sensory areas are well-captured with a linear convolution model. However, depending on the behavioral state or brain region, the coupling between neural activity and hemodynamic signals can be weak or even inverted.

Dissociable laminar profiles of concurrent bottom-up and top-down modulation in the human visual cortex

Recent developments in human neuroimaging make it possible to non-invasively measure neural activity from different cortical layers. This can potentially reveal not only which brain areas are engaged by a task, but also how. Specifically, bottom-up and top-down responses are associated with distinct laminar profiles. Here, we measured lamina-resolved fMRI responses during a visual task designed to induce concurrent bottom-up and top-down modulations via orthogonal manipulations of stimulus contrast and feature-based attention. BOLD responses were modulated by both stimulus contrast (bottom-up) and by engaging feature-based attention (top-down). Crucially, these effects operated at different cortical depths: Bottom-up modulations were strongest in the middle cortical layer and weaker in deep and superficial layers, while top-down modulations were strongest in the superficial layers. As such, we demonstrate that laminar activity profiles can discriminate between concurrent top-down and bottom-up processing, and are diagnostic of how a brain region is activated.

A critical assessment of data quality and venous effects in sub-millimeter fMRI

Advances in hardware, pulse sequences, and reconstruction techniques have made it possible to perform functional magnetic resonance imaging (fMRI) at sub-millimeter resolution while maintaining high spatial coverage and acceptable signal-to-noise ratio. Here, we examine whether sub-millimeter fMRI can be used as a routine method for obtaining accurate measurements of fine-scale local neural activity. We conducted fMRI in human visual cortex during a simple event-related visual experiment (7 T, gradient-echo EPI, 0.8-mm isotropic voxels, 2.2-s sampling rate, 84 slices), and developed analysis and visualization tools to assess the quality of the data. Our results fall along three lines of inquiry. First, we find that the acquired fMRI images, combined with appropriate surface-based processing, provide reliable and accurate measurements of fine-scale blood oxygenation level dependent (BOLD) activity patterns. Second, we show that the highly folded structure of cortex causes substantial biases on spatial resolution and data visualization. Third, we examine the well-recognized issue of venous contributions to fMRI signals. In a systematic assessment of large sections of cortex measured at a fine scale, we show that time-averaged T2*-weighted EPI intensity is a simple, robust marker of venous effects. These venous effects are unevenly distributed across cortex, are more pronounced in gyri and outer cortical depths, and are, to a certain degree, in consistent locations across subjects relative to cortical folding. Furthermore, we show that these venous effects are strongly correlated with BOLD responses evoked by the experiment. We conclude that sub-millimeter fMRI can provide robust information about fine-scale BOLD activity patterns, but special care must be exercised in visualizing and interpreting these patterns, especially with regards to the confounding influence of the brain's vasculature. To help translate these methodological findings to neuroscience research, we provide practical suggestions for both high-resolution and standard-resolution fMRI studies.

Laminar signal extraction over extended cortical areas by means of a spatial GLM

There is converging evidence that distinct neuronal processes leave distinguishable footprints in the laminar BOLD response. However, even though the achievable spatial resolution in functional MRI has much improved over the years, it is still challenging to separate signals arising from different cortical layers. In this work, we propose a new method to extract laminar signals. We use a spatial General Linear Model in combination with the equivolume principle of cortical layers to unmix laminar signals instead of interpolating through and integrating over a cortical area: thus reducing partial volume effects. Not only do we provide a mathematical framework for extracting laminar signals with a spatial GLM, we also illustrate that the best case scenarios of existing methods can be seen as special cases within the same framework. By means of simulation, we show that this approach has a sharper point spread function, providing better signal localisation. We further assess the partial volume contamination in cortical profiles from high resolution human ex vivo and in vivo structural data, and provide a full account of the benefits and potential caveats. We eschew here any attempt to validate the spatial GLM on the basis of fMRI data as a generally accepted ground-truth pattern of laminar activation does not currently exist. This approach is flexible in terms of the number of layers and their respective thickness, and naturally integrates spatial regularisation along the cortex, while preserving laminar specificity. Care must be taken, however, as this procedure of unmixing is susceptible to sources of noise in the data or inaccuracies in the laminar segmentation.

Intracortical smoothing of small-voxel fMRI data can provide increased detection power without spatial resolution losses compared to conventional large-voxel fMRI data

Continued improvement in MRI acquisition technology has made functional MRI (fMRI) with small isotropic voxel sizes down to 1 mm and below more commonly available. Although many conventional fMRI studies seek to investigate regional patterns of cortical activation for which conventional voxel sizes of 3 mm and larger provide sufficient spatial resolution, smaller voxels can help avoid contamination from adjacent white matter (WM) and cerebrospinal fluid (CSF), and thereby increase the specificity of fMRI to signal changes within the gray matter. Unfortunately, temporal signal-to-noise ratio (tSNR), a metric of fMRI sensitivity, is reduced in high-resolution acquisitions, which offsets the benefits of small voxels. Here we introduce a framework that combines small, isotropic fMRI voxels acquired at 7 T field strength with a novel anatomically-informed, surface mesh-navigated spatial smoothing that can provide both higher detection power and higher resolution than conventional voxel sizes. Our smoothing approach uses a family of intracortical surface meshes and allows for kernels of various shapes and sizes, including curved 3D kernels that adapt to and track the cortical folding pattern. Our goal is to restrict smoothing to the cortical gray matter ribbon and avoid noise contamination from CSF and signal dilution from WM via partial volume effects. We found that the intracortical kernel that maximizes tSNR does not maximize percent signal change (ΔS/S), and therefore the kernel configuration that optimizes detection power cannot be determined from tSNR considerations alone. However, several kernel configurations provided a favorable balance between boosting tSNR and ΔS/S, and allowed a 1.1-mm isotropic fMRI acquisition to have higher performance after smoothing (in terms of both detection power and spatial resolution) compared to an unsmoothed 3.0-mm isotropic fMRI acquisition. Overall, the results of this study support the strategy of acquiring voxels smaller than the cortical thickness, even for studies not requiring high spatial resolution, and smoothing them down within the cortical ribbon with a kernel of an appropriate shape to achieve the best performance-thus decoupling the choice of fMRI voxel size from the spatial resolution requirements of the particular study. The improvement of this new intracortical smoothing approach over conventional surface-based smoothing is expected to be modest for conventional resolutions, however the improvement is expected to increase with higher resolutions. This framework can also be applied to anatomically-informed intracortical smoothing of higher-resolution data (e.g. along columns and layers) in studies with prior information about the spatial structure of activation.

Resolving laminar activation in human V1 using ultra-high spatial resolution fMRI at 7T

The mesoscopic organization of the human neocortex is of great interest for cognitive neuroscience. However, fMRI in humans typically maps the functional units of cognitive processing on a macroscopic level. With the advent of ultra-high field MRI (≥7T), it has become possible to acquire fMRI data with sub-millimetre resolution, enabling probing the laminar and columnar circuitry in humans. Currently, laminar BOLD responses are not directly observed but inferred via data analysis, due to coarse spatial resolution of fMRI (e.g. 0.7-0.8 mm isotropic) relative to the extent of histological laminae. In this study, we introduce a novel approach for mapping the cortical BOLD response at the spatial scale of cortical layers and columns at 7T (an unprecedented 0.1 mm, either in the laminar or columnar direction). We demonstrate experimentally and using simulations, the superiority of the novel approach compared to standard approaches for human laminar fMRI in terms of effective spatial resolution in either laminar or columnar direction. In addition, we provide evidence that the laminar BOLD signal profile is not homogeneous even over short patches of cortex. In summary, the proposed novel approach affords the ability to directly study the mesoscopic organization of the human cortex, thus, bridging the gap between human cognitive neuroscience and invasive animal studies.

Simulations of blood as a suspension predicts a depth dependent hematocrit in the circulation throughout the cerebral cortex

Recent advances in modeling oxygen supply to cortical brain tissue have begun to elucidate the functional mechanisms of neurovascular coupling. While the principal mechanisms of blood flow regulation after neuronal firing are generally known, mechanistic hemodynamic simulations cannot yet pinpoint the exact spatial and temporal coordination between the network of arteries, arterioles, capillaries and veins for the entire brain. Because of the potential significance of blood flow and oxygen supply simulations for illuminating spatiotemporal regulation inside the cortical microanatomy, there is a need to create mathematical models of the entire cerebral circulation with realistic anatomical detail. Our hypothesis is that an anatomically accurate reconstruction of the cerebrocirculatory architecture will inform about possible regulatory mechanisms of the neurovascular interface. In this article, we introduce large-scale networks of the murine cerebral circulation spanning the Circle of Willis, main cerebral arteries connected to the pial network down to the microcirculation in the capillary bed. Several multiscale models were generated from state-of-the-art neuroimaging data. Using a vascular network construction algorithm, the entire circulation of the middle cerebral artery was synthesized. Blood flow simulations indicate a consistent trend of higher hematocrit in deeper cortical layers, while surface layers with shorter vascular path lengths seem to carry comparatively lower red blood cell (RBC) concentrations. Moreover, the variability of RBC flux decreases with cortical depth. These results support the notion that plasma skimming serves a self-regulating function for maintaining uniform oxygen perfusion to neurons irrespective of their location in the blood supply hierarchy. Our computations also demonstrate the practicality of simulating blood flow for large portions of the mouse brain with existing computer resources. The efficient simulation of blood flow throughout the entire middle cerebral artery (MCA) territory is a promising milestone towards the final aim of predicting blood flow patterns for the entire brain.

The brain’s astonishing cognitive capacity depends on the coordination between neurons and the cerebral circulation, a system known as the neurovascular unit. The spatial and temporal coupling between these two networks is the object of intense research. However, the concise anatomical description of the cerebral circulation has so far been intractable. This paper introduces a methodology for the in silico creation of realistic models for the entire cerebral circulation. This innovation incorporates topological data from several neuroimaging modalities covering three lengths scales as input into a computer algorithm, which assembles anatomically accurate circulatory networks. When simulating blood flow as red blood cells suspended in plasma for experimental and synthetic cortical network models, we discovered that red blood cells tend to be more concentrated in deeper layers of the cortex compared to the surface. RBC fluxes are more homogenous in deeper layers. The phenomenon of depth dependent red blood cell supply supports the notion that the intricate architecture of the cortical microcirculation serves a self-regulating function to maintain uniform oxygen perfusion to neurons. We also demonstrate the practicality of predicting blood flow patterns for the entire brain with existing computer power.

Addressing challenges of high spatial resolution UHF fMRI for group analysis of higher-order cognitive tasks: An inter-sensory task directing attention between visual and somatosensory domains

Functional MRI at ultra‐high field (UHF, ≥7 T) provides significant increases in BOLD contrast‐to‐noise ratio (CNR) compared with conventional field strength (3 T), and has been exploited for reduced field‐of‐view, high spatial resolution mapping of primary sensory areas. Applying these high spatial resolution methods to investigate whole brain functional responses to higher‐order cognitive tasks leads to a number of challenges, in particular how to perform robust group‐level statistical analyses. This study addresses these challenges using an inter‐sensory cognitive task which modulates top‐down attention at graded levels between the visual and somatosensory domains. At the individual level, highly focal functional activation to the task and task difficulty (modulated by attention levels) were detectable due to the high CNR at UHF. However, to assess group level effects, both anatomical and functional variability must be considered during analysis. We demonstrate the importance of surface over volume normalisation and the requirement of no spatial smoothing when assessing highly focal activity. Using novel group analysis on anatomically parcellated brain regions, we show that in higher cognitive areas (parietal and dorsal‐lateral‐prefrontal cortex) fMRI responses to graded attention levels were modulated quadratically, whilst in visual cortex and VIP, responses were modulated linearly. These group fMRI responses were not seen clearly using conventional second‐level GLM analyses, illustrating the limitations of a conventional approach when investigating such focal responses in higher cognitive regions which are more anatomically variable. The approaches demonstrated here complement other advanced analysis methods such as multivariate pattern analysis, allowing UHF to be fully exploited in cognitive neuroscience.

Existence of Initial Dip for BCI: An Illusion or Reality

A tight coupling between the neuronal activity and the cerebral blood flow (CBF) is the motivation of many hemodynamic response (HR)-based neuroimaging modalities. The increase in neuronal activity causes the increase in CBF that is indirectly measured by HR modalities. Upon functional stimulation, the HR is mainly categorized in three durations: (i) initial dip, (ii) conventional HR (i.e., positive increase in HR caused by an increase in the CBF), and (iii) undershoot. The initial dip is a change in oxygenation prior to any subsequent increase in CBF and spatially more specific to the site of neuronal activity. Despite additional evidence from various HR modalities on the presence of initial dip in human and animal species (i.e., cat, rat, and monkey); the existence/occurrence of an initial dip in HR is still under debate. This article reviews the existence and elusive nature of the initial dip duration of HR in intrinsic signal optical imaging (ISOI), functional magnetic resonance imaging (fMRI), and functional near-infrared spectroscopy (fNIRS). The advent of initial dip and its elusiveness factors in ISOI and fMRI studies are briefly discussed. Furthermore, the detection of initial dip and its role in brain-computer interface using fNIRS is examined in detail. The best possible application for the initial dip utilization and its future implications using fNIRS are provided.

Laminar Organization of Working Memory Signals in Human Visual Cortex

The human primary visual cortex (V1) is not only activated by incoming visual information but is also engaged by top-down cognitive processes, such as visual working memory, even in the absence of visual input [1-3]. This feedback may be critical to our ability to visualize specific visual features, as higher-order regions lack the selectivity to represent such information [4]. Clearly, such internally generated signals do not trigger genuine perception of the remembered stimulus, meaning they must be organized in a manner that is different to bottom-up-driven signals. Internally generated signals may be kept separate from incoming sensory data by virtue of the laminar organization of inter-area cortical connections. Namely, bottom-up driving connections target layer 4, located in the middle of the cortical column, and feedback connections target deep and superficial layers and avoid layer 4 [5-7]. Using lamina-resolved fMRI, we simultaneously measured the activity in three early visual cortical areas (V1-V3) that are recruited to represent stimulus information during visual working memory [8]. We observed item-specific working memory signals in early visual cortex. In V1, this item-specific activity was selectively present at deep and superficial cortical depths, avoiding the middle layers, and working-memory-related activity was present at all depths in V2 and V3. These results show for the first time the laminar organization of internally generated signals during visual working memory in the human visual system and provide new insights into how bottom-up and top-down signals in visual cortex are deployed.

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

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

Biophysics of BOLD fMRI investigated with animal models

The widely-used BOLD fMRI signal depends on various anatomical, physiological, and imaging parameters. Thus, it is important to examine its biophysical and physiological source in order to optimize, model and accurately interpret fMRI. Animal models have been used to investigate these issues to take systematic measurements and combine with conventional invasive approaches. Here, we reviewed and discussed multiple issues, including the echo time-dependent intravascular contribution and extravascular contributions, gradient-echo vs. spin-echo fMRI, the physiological source of BOLD fMRI, arterial vs. venous cerebral blood volume change, cerebral oxygen consumption change, and arterial oxygen saturation change. We then discuss future directions of animal fMRI and translation to human fMRI. Systematic biophysical BOLD fMRI studies provide insight into the modeling and interpretation of BOLD fMRI in animals and humans.

Simulating laminar neuroimaging data for a visual delayed match-to-sample task

Invasive electrophysiological and neuroanatomical studies in nonhuman mammalian experimental preparations have helped elucidate the lamina (layer) dependence of neural computations and interregional connections. Noninvasive functional neuroimaging can, in principle, resolve cortical laminae (layers), and thus provide insight into human neural computations and interregional connections. However human neuroimaging data are noisy and difficult to interpret; biologically realistic simulations can aid experimental interpretation by relating the neuroimaging data to simulated neural activity. We illustrate the potential of laminar neuroimaging by upgrading an existing large-scale, multiregion neural model that simulates a visual delayed match-to-sample (DMS) task. The new laminar-based neural unit incorporates spiny stellate, pyramidal, and inhibitory neural populations which are divided among supragranular, granular, and infragranular laminae (layers). We simulated neural activity which is translated into local field potential-like data used to simulate conventional and laminar fMRI activity. We implemented the laminar connectivity schemes proposed by Felleman and Van Essen (Cerebral Cortex, 1991) for interregional connections. The hemodynamic model that we employ is a modified version of one due to Heinzle et al. (Neuroimage, 2016) that incorporates the effects of draining veins. We show that the laminar version of the model replicates the findings of the existing model. The laminar model shows the finer structure in fMRI activity and functional connectivity. Laminar differences in the magnitude of neural activities are a prominent finding; these are also visible in the simulated fMRI. We illustrate differences between task and control conditions in the fMRI signal, and demonstrate differences in interregional laminar functional connectivity that reflect the underlying connectivity scheme. These results indicate that multi-layer computational models can aid in interpreting layer-specific fMRI, and suggest that increased use of laminar fMRI could provide unique and fundamental insights to human neuroscience.

A guide to using functional magnetic resonance imaging to study Alzheimer's disease in animal models

Alzheimer's disease is a leading healthcare challenge facing our society today. Functional magnetic resonance imaging (fMRI) of the brain has played an important role in our efforts to understand how Alzheimer's disease alters brain function. Using fMRI in animal models of Alzheimer's disease has the potential to provide us with a more comprehensive understanding of the observations made in human clinical fMRI studies. However, using fMRI in animal models of Alzheimer's disease presents some unique challenges. Here, we highlight some of these challenges and discuss potential solutions for researchers interested in performing fMRI in animal models. First, we briefly summarize our current understanding of Alzheimer's disease from a mechanistic standpoint. We then overview the wide array of animal models available for studying this disease and how to choose the most appropriate model to study, depending on which aspects of the condition researchers seek to investigate. Finally, we discuss the contributions of fMRI to our understanding of Alzheimer's disease and the issues to consider when designing fMRI studies for animal models, such as differences in brain activity based on anesthetic choice and ways to interrogate more specific questions in rodents beyond those that can be addressed in humans. The goal of this article is to provide information on the utility of fMRI, and approaches to consider when using fMRI, for studies of Alzheimer's disease in animal models.

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

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

Anatomic & metabolic brain markers of the m.3243A>G mutation: A multi-parametric 7T MRI study

One of the most common mitochondrial DNA (mtDNA) mutations, the A to G transition at base pair 3243, has been linked to changes in the brain, in addition to commonly observed hearing problems, diabetes and myopathy. However, a detailed quantitative description of m.3243A>G patients' brains has not been provided so far. In this study, ultra-high field MRI at 7T and volume- and surface-based data analyses approaches were used to highlight morphology (i.e. atrophy)-, microstructure (i.e. myelin and iron concentration)- and metabolism (i.e. cerebral blood flow)-related differences between patients (N = 22) and healthy controls (N = 15). The use of quantitative MRI at 7T allowed us to detect subtle changes of biophysical processes in the brain with high accuracy and sensitivity, in addition to typically assessed lesions and atrophy. Furthermore, the effect of m.3243A>G mutation load in blood and urine epithelial cells on these MRI measures was assessed within the patient population and revealed that blood levels were most indicative of the brain's state and disease severity, based on MRI as well as on neuropsychological data. Morphometry MRI data showed a wide-spread reduction of cortical, subcortical and cerebellar gray matter volume, in addition to significantly enlarged ventricles. Moreover, surface-based analyses revealed brain area-specific changes in cortical thickness (e.g. of the auditory cortex), and in T₁, T₂* and cerebral blood flow as a function of mutation load, which can be linked to typically m.3243A>G-related clinical symptoms (e.g. hearing impairment). In addition, several regions linked to attentional control (e.g. middle frontal gyrus), the sensorimotor network (e.g. banks of central sulcus) and the default mode network (e.g. precuneus) were characterized by alterations in cortical thickness, T₁, T₂* and/or cerebral blood flow, which has not been described in previous MRI studies. Finally, several hypotheses, based either on vascular, metabolic or astroglial implications of the m.3243A>G mutation, are discussed that potentially explain the underlying pathobiology. To conclude, this is the first 7T and also the largest MRI study on this patient population that provides macroscopic brain correlates of the m.3243A>G mutation indicating potential MRI biomarkers of mitochondrial diseases and might guide future (longitudinal) studies to extensively track neuropathological and clinical changes.

Precise Cerebral Vascular Atlas in Stereotaxic Coordinates of Whole Mouse Brain

Understanding amazingly complex brain functions and pathologies requires a complete cerebral vascular atlas in stereotaxic coordinates. Making a precise atlas for cerebral arteries and veins has been a century-old objective in neuroscience and neuropathology. Using micro-optical sectioning tomography (MOST) with a modified Nissl staining method, we acquired five mouse brain data sets containing arteries, veins, and microvessels. Based on the brain-wide vascular spatial structures and brain regions indicated by cytoarchitecture in one and the same mouse brain, we reconstructed and annotated the vascular system atlas of both arteries and veins of the whole mouse brain for the first time. The distributing patterns of the vascular system within the brain regions were acquired and our results show that the patterns of individual vessels are different from each other. Reconstruction and statistical analysis of the microvascular network, including derivation of quantitative vascular densities, indicate significant differences mainly in vessels with diameters less than 8 μm and large than 20 μm across different brain regions. Our precise cerebral vascular atlas provides an important resource and approach for quantitative studies of brain functions and diseases.

High-Resolution CBV-fMRI Allows Mapping of Laminar Activity and Connectivity of Cortical Input and Output in Human M1

Layer-dependent fMRI allows measurements of information flow in cortical circuits, as afferent and efferent connections terminate in different cortical layers. However, it is unknown to what level human fMRI is specific and sensitive enough to reveal directional functional activity across layers. To answer this question, we developed acquisition and analysis methods for blood-oxygen-level-dependent (BOLD) and cerebral-blood-volume (CBV)-based laminar fMRI and used these to discriminate four different tasks in the human motor cortex (M1). In agreement with anatomical data from animal studies, we found evidence for somatosensory and premotor input in superficial layers of M1 and for cortico-spinal motor output in deep layers. Laminar resting-state fMRI showed directional functional connectivity of M1 with somatosensory and premotor areas. Our findings demonstrate that CBV-fMRI can be used to investigate cortical activity in humans with unprecedented detail, allowing investigations of information flow between brain regions and outperforming conventional BOLD results that are often buried under vascular biases.