Combined imaging–histological study of cortical laminar specificity of fMRI signals

Machine Learning-Based Clustering of Layer-Resolved fMRI Activation and Functional Connectivity Within the Primary Somatosensory Cortex in Nonhuman Primates

Delineating the functional organization of mesoscale cortical columnar structure is essential for understanding brain function. We have previously demonstrated a high spatial correspondence between BOLD fMRI and LFP responses to tactile stimuli in the primary somatosensory cortex area 3b of nonhuman primates. This study aims to explore how 2D spatial profiles of the functional column vary across cortical layers (defined by three cortical depths) in both tactile stimulation and resting states using fMRI. At 9.4 T, we acquired submillimeter-resolution oblique fMRI data from cortical areas 3b and 1 of anesthetized squirrel monkeys and obtained fMRI signals from three cortical layers. In both areas 3b and 1, the tactile stimulus-evoked fMRI activation foci were fitted with point spread functions (PSFs), from which shape parameters, including full width at half maximum (FWHM), were derived. Seed-based resting-state fMRI data analysis was then performed to measure the spatial profiles of resting-state connectivity within and between areas 3b and 1. We found that the tactile-evoked fMRI response and local resting-state functional connectivity were elongated at the superficial layer, with the major axes oriented in lateral to medial (from digit 1 to digit 5) direction. This elongation was significantly reduced in the deeper (middle and bottom) layers. To assess the robustness of these spatial profiles in distinguishing cortical layers, shape parameters describing the spatial extents of activation and resting-state connectivity profiles were used to classify the layers via self-organizing maps (SOM). A minimal overall classification error (~13%) was achieved, effectively classifying the layers into two groups: the superficial layer exhibited distinct features from the two deeper layers in the rsfMRI data. Our results support distinct 2D spatial profiles for superficial versus deeper cortical layers and reveal similarities between stimulus-evoked and resting-state configurations.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Brain-wide mapping of resting-state networks in mice using high-frame rate functional ultrasound

Functional ultrasound (fUS) imaging is a method for visualizing deep brain activity based on cerebral blood volume changes coupled with neural activity, while functional MRI (fMRI) relies on the blood-oxygenation-level-dependent signal coupled with neural activity. Low-frequency fluctuations (LFF) of fMRI signals during resting-state can be measured by resting-state fMRI (rsfMRI), which allows functional imaging of the whole brain, and the distributions of resting-state network (RSN) can then be estimated from these fluctuations using independent component analysis (ICA). This procedure provides an important method for studying cognitive and psychophysiological diseases affecting specific brain networks. The distributions of RSNs in the brain-wide area has been reported primarily by rsfMRI. RSNs using rsfMRI are generally computed from the time-course of fMRI signals for more than 5 min. However, a recent dynamic functional connectivity study revealed that RSNs are still not perfectly stable even after 10 min. Importantly, fUS has a higher temporal resolution and stronger correlation with neural activity compared with fMRI. Therefore, we hypothesized that fUS applied during the resting-state for a shorter than 5 min would provide similar RSNs compared to fMRI. High temporal resolution rsfUS data were acquired at 10 Hz in awake mice. The quality of the default mode network (DMN), a well-known RSN, was evaluated using signal-noise separation (SNS) applied to different measurement durations of rsfUS. The results showed that the SNS did not change when the measurement duration was increased to more than 210 s. Next, we measured short-duration rsfUS multi-slice measurements in the brain-wide area. The results showed that rsfUS with the short duration succeeded in detecting RSNs distributed in the brain-wide area consistent with RSNs detected by 11.7-T MRI under awake conditions (medial prefrontal cortex and cingulate cortex in the anterior DMN, retrosplenial cortex and visual cortex in the posterior DMN, somatosensory and motor cortexes in the lateral cortical network, thalamus, dorsal hippocampus, and medial cerebellum), confirming the reliability of the RSNs detected by rsfUS. However, bilateral RSNs located in the secondary somatosensory cortex, ventral hippocampus, auditory cortex, and lateral cerebellum extracted from rsfUS were different from the unilateral RSNs extracted from rsfMRI. These findings indicate the potential of rsfUS as a method for analyzing functional brain networks and should encourage future research to elucidate functional brain networks and their relationships with disease model mice.

Spatially heterogeneous structure-function coupling in haemodynamic and electromagnetic brain networks

The relationship between structural and functional connectivity in the brain is a key question in connectomics. Here we quantify patterns of structure-function coupling across the neocortex, by comparing structural connectivity estimated using diffusion MRI with functional connectivity estimated using both neurophysiological (MEG-based) and haemodynamic (fMRI-based) recordings. We find that structure-function coupling is heterogeneous across brain regions and frequency bands. The link between structural and functional connectivity is generally stronger in multiple MEG frequency bands compared to resting state fMRI. Structure-function coupling is greater in slower and intermediate frequency bands compared to faster frequency bands. We also find that structure-function coupling systematically follows the archetypal sensorimotor-association hierarchy, as well as patterns of laminar differentiation, peaking in granular layer IV. Finally, structure-function coupling is better explained using structure-informed inter-regional communication metrics than using structural connectivity alone. Collectively, these results place neurophysiological and haemodynamic structure-function relationships in a common frame of reference and provide a starting point for a multi-modal understanding of structure-function coupling in the brain.

Does V1 response suppression initiate binocular rivalry?

During binocular rivalry (BR) only one eye's view is perceived. Neural underpinnings of BR are debated. Recent studies suggest that primary visual cortex (V1) initiates BR. One trigger might be response suppression across most V1 neurons at the onset of BR. Here, we utilize a variant of BR called binocular rivalry flash suppression (BRFS) to test this hypothesis. BRFS is identical to BR, except stimuli are shown with a ∼1s delay. If V1 response suppression was required to initiate BR, it should occur during BRFS as well. To test this, we compared V1 spiking in two macaques observing BRFS. We found that BRFS resulted in response facilitation rather than response suppression across V1 neurons. However, BRFS still reduces responses in a subset of V1 neurons due to the adaptive effects of asynchronous stimulus presentation. We argue that this selective response suppression could serve as an alternate initiator of BR.

Design and application of a multimodality-compatible 1Tx/6Rx RF coil for monkey brain MRI at 7T

OBJECTIVE

Blood-oxygen-level-dependent functional MRI allows to investigte neural activities and connectivity. While the non-human primate plays an essential role in neuroscience research, multimodal methods combining functional MRI with other neuroimaging and neuromodulation enable us to understand the brain network at multiple scales.

APPROACH

In this study, a tight-fitting helmet-shape receive array with a single transmit loop for anesthetized macaque brain MRI at 7T was fabricated with four openings constructed in the coil housing to accommodate multimodal devices, and the coil performance was quantitatively evaluated and compared to a commercial knee coil. In addition, experiments over three macaques with infrared neural stimulation (INS), focused ultrasound stimulation (FUS), and transcranial direct current stimulation (tDCS) were conducted.

MAIN RESULTS

The RF coil showed higher transmit efficiency, comparable homogeneity, improved SNR and enlarged signal coverage over the macaque brain. Infrared neural stimulation was applied to the amygdala in deep brain region, and activations in stimulation sites and connected sites were detected, with the connectivity consistent with anatomical information. Focused ultrasound stimulation was applied to the left visual cortex, and activations were acquired along the ultrasound traveling path, with all time course curves consistent with pre-designed paradigms. The existence of transcranial direct current stimulation electrodes brought no interference to the RF system, as evidenced through high-resolution MPRAGE structure images.

SIGNIFICANCE

This pilot study reveals the feasibility for brain investigation at multiple spatiotemporal scales, which may advance our understanding in dynamic brain networks.

Contrast Enhancement in MRI Using Combined Double Action Contrast Agents and Image Post-Processing in the Breast Cancer Model

Gd- and Fe-based contrast agents reduce T₁ and T₂ relaxation times, respectively, are frequently used in MRI, providing improved cancer detection. Recently, contrast agents changing both T₁/T₂ times, based on core/shell nanoparticles, have been introduced. Although advantages of the T₁/T₂ agents were shown, MR image contrast of cancerous versus normal adjacent tissue induced by these agents has not yet been analyzed in detail as authors considered changes in cancer MR signal or signal-to-noise ratio after contrast injection rather than changes in signal differences between cancer and normal adjacent tissue. Furthermore, the potential advantages of T₁/T₂ contrast agents using image manipulation such as subtraction or addition have not been yet discussed in detail. Therefore, we performed theoretical calculations of MR signal in a tumor model using T₁-weighted, T₂-weighted, and combined images for T₁-, T₂-, and T₁/T₂-targeted contrast agents. The results from the tumor model are followed by in vivo experiments using core/shell NaDyF₄/NaGdF₄ nanoparticles as T₁/T₂ non-targeted contrast agent in the animal model of triple negative breast cancer. The results show that subtraction of T₂-weighted from T₁-weighted MR images provides additional increase in the tumor contrast: over two-fold in the tumor model and 12% in the in vivo experiment.

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

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

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic-haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

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.

Improvement of sensitivity and specificity for laminar BOLD fMRI with double spin-echo EPI in humans at 7 T

Mapping mesoscopic cortical functional units such as columns or laminae is increasingly pursued by ultra-high field (UHF) functional magnetic resonance imaging (fMRI). The most popular approach for high-resolution fMRI is currently gradient-echo (GE) blood oxygenation level-dependent (BOLD) fMRI. However, its spatial accuracy is reduced due to its sensitivity to draining vessels, including pial veins, whereas spin-echo (SE) BOLD signal is expected to have higher spatial accuracy, albeit with lower sensitivity than the GE-BOLD signal. Here, we introduce a new double spin-echo (dSE) echo-planar imaging (EPI) method to improve the sensitivity of SE-BOLD contrast by averaging two spin-echoes using three radiofrequency pulses. Human fMRI experiments were performed with slices perpendicular to the central sulcus between motor and sensory cortices at 7 T during fist-clenching with touching. First, we evaluated the feasibility of single-shot dSE-EPI for BOLD fMRI with 1.5 mm isotropic resolution and found that dSE-BOLD fMRI has higher signal-to-noise ratio (SNR), temporal SNR (tSNR), and higher functional sensitivity than conventional SE-BOLD fMRI. Second, to investigate the laminar specificity of dSE-BOLD fMRI, we implemented a multi-shot approach to achieve 0.8-mm isotropic resolution with sliding-window reconstruction. Unlike GE-BOLD fMRI, the cortical profile of dSE-BOLD fMRI peaked at ~ 1.0 mm from the surface of the primary motor and sensory cortices, demonstrating an improvement of laminar specificity in humans over GE-BOLD fMRI. The proposed multi-shot dSE-EPI method is viable for high spatial resolution UHF-fMRI studies in the pursuit of resolving mesoscopic functional units.

One-pot synthesis of carboxymethyl-dextran coated iron oxide nanoparticles (CION) for preclinical fMRI and MRA applications

Superparamagnetic iron-oxide nanoparticles are robust contrast agents for magnetic resonance imaging (MRI) used for sensitive structural and functional mapping of the cerebral blood volume (CBV) when administered intravenously. To date, many CBV-MRI studies are conducted with Feraheme, manufactured for the clinical treatment of iron-deficiency. Unfortunately, Feraheme is currently not available outside the United States due to commercial and regulatory constraints, making CBV-MRI methods either inaccessible or very costly to achieve. To address this barrier, we developed a simple, one-pot recipe to synthesize Carboxymethyl-dextran coated Iron Oxide Nanoparticles, namely, “CION”, suitable for preclinical CBV-MRI applications. Here we disseminate a step-by-step instruction of our one-pot synthesis protocol, which allows CION to be produced in laboratories with minimal cost. We also characterized different CION-conjugations by manipulating polymer to metal stoichiometric ratio in terms of their size, surface chemistry, and chemical composition, and shifts in MR relaxivity and pharmacokinetics. We performed several proof-of-concept experiments in vivo, demonstrating the utility of CION for functional and structural MRI applications, including hypercapnic CO₂ challenge, visual stimulation, targeted optogenetic stimulation, and microangiography. We also present evidence that CION can serve as a cross-modality research platform by showing concurrent in vivo optical and MRI measurement of CBV using fluorescent-labeled CION. The simplicity and cost-effectiveness of our one-pot synthesis method should allow researchers to reproduce CION and tailor the relaxivity and pharmacokinetics according to their imaging needs. It is our hope that this work makes CBV-MRI more openly available and affordable for a variety of research applications.

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.

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.

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.

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.

Accelerated spin-echo functional MRI using multisection excitation by simultaneous spin-echo interleaving (MESSI) with complex-encoded generalized slice dithered enhanced resolution (cgSlider) simultaneous multislice echo-planar imaging

PURPOSE

Spin-echo functional MRI (SE-fMRI) has the potential to improve spatial specificity when compared with gradient-echo fMRI. However, high spatiotemporal resolution SE-fMRI with large slice-coverage is challenging as SE-fMRI requires a long echo time to generate blood oxygenation level-dependent (BOLD) contrast, leading to long repetition times. The aim of this work is to develop an acquisition method that enhances the slice-coverage of SE-fMRI at high spatiotemporal resolution.

THEORY AND METHODS

An acquisition scheme was developed entitled multisection excitation by simultaneous spin-echo interleaving (MESSI) with complex-encoded generalized slice dithered enhanced resolution (cgSlider). MESSI uses the dead-time during the long echo time by interleaving the excitation and readout of 2 slices to enable 2× slice-acceleration, while cgSlider uses the stable temporal background phase in SE-fMRI to encode/decode 2 adjacent slices simultaneously with a "phase-constrained" reconstruction method. The proposed cgSlider-MESSI was also combined with simultaneous multislice (SMS) to achieve further slice-acceleration. This combined approach was used to achieve 1.5-mm isotropic whole-brain SE-fMRI with a temporal resolution of 1.5 s and was evaluated using sensory stimulation and breath-hold tasks at 3T.

RESULTS

Compared with conventional SE-SMS, cgSlider-MESSI-SMS provides 4-fold increase in slice-coverage for the same repetition time, with comparable temporal signal-to-noise ratio. Corresponding fMRI activation from cgSlider-MESSI-SMS for both fMRI tasks were consistent with those from conventional SE-SMS. Overall, cgSlider-MESSI-SMS achieved a 32× encoding-acceleration by combining Rinplane × MB × cgSlider × MESSI = 4 × 2 × 2 × 2.

CONCLUSION

High-quality, high-resolution whole-brain SE-fMRI was acquired at a short repetition time using cgSlider-MESSI-SMS. This method should be beneficial for high spatiotemporal resolution SE-fMRI studies requiring whole-brain coverage.

Comparison of BOLD and CBV using 3D EPI and 3D GRASE for cortical layer functional MRI at 7 T

PURPOSE

Functional MRI (fMRI) at the mesoscale of cortical layers and columns requires both sensitivity and specificity, the latter of which can be compromised if the imaging method is affected by vascular artifacts, particularly cortical draining veins at the pial surface. Recent studies have shown that cerebral blood volume (CBV) imaging is more specific to the actual laminar locus of neural activity than BOLD imaging using standard gradient-echo EPI sequences. Gradient and spin-echo (GRASE) BOLD imaging has also shown greater specificity when compared with standard gradient-echo EPI BOLD. Here we directly compare CBV and BOLD contrasts in high-resolution imaging of the primary motor cortex for laminar functional MRI in four combinations of signal labeling, CBV using slice-selective slab-inversion vascular space occupancy (VASO) and BOLD, each with 3D gradient-echo EPI and zoomed 3D-GRASE image readouts.

METHODS

Activations were measured using each sequence and contrast combination during a motor task. Activation profiles across cortical depth were measured to assess the sensitivity and specificity (pial bias) of each method.

RESULTS

Both CBV imaging using gradient-echo 3D-EPI and BOLD imaging using 3D-GRASE show similar specificity and sensitivity and are therefore useful tools for mesoscopic functional MRI in the human cortex. The combination of GRASE and VASO did not demonstrate high levels of sensitivity, nor show increased specificity.

CONCLUSION

Three-dimensional EPI with VASO contrast and 3D-GRASE with BOLD contrast both demonstrate sufficient sensitivity and specificity for laminar functional MRI to be used by neuroscientists in a wide range of investigations of depth-dependent neural circuitry in the human brain.

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.

Resolving multisensory and attentional influences across cortical depth in sensory cortices

In our environment, our senses are bombarded with a myriad of signals, only a subset of which is relevant for our goals. Using sub-millimeter-resolution fMRI at 7T, we resolved BOLD-response and activation patterns across cortical depth in early sensory cortices to auditory, visual and audiovisual stimuli under auditory or visual attention. In visual cortices, auditory stimulation induced widespread inhibition irrespective of attention, whereas auditory relative to visual attention suppressed mainly central visual field representations. In auditory cortices, visual stimulation suppressed activations, but amplified responses to concurrent auditory stimuli, in a patchy topography. Critically, multisensory interactions in auditory cortices were stronger in deeper laminae, while attentional influences were greatest at the surface. These distinct depth-dependent profiles suggest that multisensory and attentional mechanisms regulate sensory processing via partly distinct circuitries. Our findings are crucial for understanding how the brain regulates information flow across senses to interact with our complex multisensory world.

The cerebral angiome: High resolution MicroCT imaging of the whole brain cerebrovasculature in female and male mice

The cerebrovascular system provides crucial functions that maintain metabolic and homeostatic states of the brain. Despite its integral role of supporting cerebral viability, the topological organization of these networks remains largely uncharacterized. This void in our knowledge surmises entirely from current technological limitations that prevent the capturing of data through the entire depth of the brain. We report high-resolution reconstruction and analysis of the complete vascular network of the entire brain at the capillary level in adult female and male mice using a vascular corrosion cast procedure. Vascular network analysis of the whole brain revealed sex-related differences of vessel hierarchy. In addition, region-specific network analysis demonstrated different patterns of angioarchitecture between brain subregions and sex. Furthermore, our group is the first to provide a three-dimensional analysis of the angioarchitecture and network organization in a single reconstructed tomographic data set that encompasses all hierarchy of vessels in the brain of the adult mouse.

Myelin and Modeling: Bootstrapping Cortical Microcircuits

Histological studies of myelin-stained sectioned cadaver brain and in vivo myelin-weighted magnetic resonance imaging (MRI) show that the cerebral cortex is organized into cortical areas with generally well-defined boundaries, which have consistent internal patterns of myelination. The process of myelination is largely driven by neural experience, in which the axonal passage of action potentials stimulates neighboring oligodendrocytes to perform their task. This bootstrapping process, such that the traffic of action potentials facilitates increased traffic, suggests the hypothesis that the specific pattern of myelination (myeloarchitecture) in each cortical area reveals the principal cortical microcircuits required for the function of that area. If this idea is correct, the observable sequential maturation of specific brain areas can provide evidence for models of the stages of cognitive development.

Ultrafast fMRI of the rodent brain using simultaneous multi-slice EPI

Increasing spatial and temporal resolutions of functional MRI (fMRI) measurement has been shown to benefit the study of neural dynamics and functional interaction. However, acceleration of rodent brain fMRI using parallel and simultaneous multi-slice imaging techniques is hampered by the lack of high-density phased-array coils for the small brain. To overcome this limitation, we adapted phase-offset multiplanar and blipped-controlled aliasing echo planar imaging (EPI) to enable simultaneous multi-slice fMRI of the mouse brain using a single loop coil on a 9.4T scanner. Four slice bands of 0.3 × 0.3 × 0.5 mm³ resolution can be simultaneously acquired to cover the whole brain at a temporal resolution of 300 ms or the whole cerebrum in 150 ms. Instead of losing signal-to-noise ratio (SNR), both spatial and temporal SNR can be increased due to the increased k-space sampling compared to a standard single-band EPI. Task fMRI using a visual stimulation shows close to 80% increase of z-score and 4 times increase of activated area in the visual cortex using the multiband EPI due to the highly increased temporal samples. Resting-state fMRI shows reliable detection of bilateral connectivity by both single-band and multiband EPI, but no significant difference was found. Without the need of a dedicated hardware, we have demonstrated a practical method that can enable unparallelly fast whole-brain fMRI for preclinical studies. This technique can be used to increase sensitivity, distinguish transient response or acquire high spatiotemporal resolution fMRI.

Spontaneous low-frequency fluctuations in the neural system for emotional perception in major psychiatric disorders: amplitude similarities and differences across frequency bands

BACKGROUND

Growing evidence indicates both shared and distinct features of emotional perception in schizophrenia, bipolar disorder and major depressive disorder. In these disorders, alterations in spontaneous low-frequency fluctuations have been reported in the neural system for emotional perception, but the similarities and differences in the amplitude of low-frequency fluctuation (ALFF) across the 3 disorders are unknown.

METHODS

We compared ALFF and its signal balance in the neural system for emotional perception at 2 frequency bands (slow-5 and slow-4) in 119 participants with schizophrenia, 100 with bipolar disorder, 123 with major depressive disorder and 183 healthy controls. We performed exploratory Pearson partial correlation analyses to determine the relationship between ALFF signal balance and clinical variables.

RESULTS

We observed commonalities in ALFF change patterns across the 3 disorders for emotional perception neural substrates, such as increased ALFF in the anterior cerebrum (including subcortical, limbic, paralimbic and heteromodal cortical regions) and decreased ALFF in the posterior visual cortices. Schizophrenia, bipolar disorder and major depressive disorder showed significantly decreased ALFF signal balance in the neural system for emotional perception at both slow-5 and slow-4 frequency bands, with the greatest alterations for schizophrenia, followed by bipolar disorder and major depressive disorder. We found a negative correlation between ALFF signal balance and negative/disorganized symptoms in slow-4 across the 3 disorders.

LIMITATIONS

The relatively broad age range in our sample and the cross-sectional study design may not account for our findings.

CONCLUSION

The extent of the commonalities we observed further support the concept of core neurobiological disruptions shared among the 3 disorders; ALFF signal balance could be an important neuroimaging marker for the diagnosis and treatment of schizophrenia, bipolar disorder and major depressive disorder.

Preclinical Imaging Biomarkers for Postischaemic Neurovascular Remodelling

In the pursuit of understanding the pathological alterations that underlie ischaemic injuries, such as vascular remodelling and reorganisation, there is a need for recognising the capabilities and limitations of in vivo imaging techniques. Thus, this review presents contemporary published research of imaging modalities that have been implemented to study postischaemic neurovascular changes in small animals. A comparison of the technical aspects of the various imaging tools is included to set the framework for identifying the most appropriate methods to observe postischaemic neurovascular remodelling. A systematic search of the PubMed® and Elsevier's Scopus databases identified studies that were conducted between 2008 and 2018 to explore postischaemic neurovascular remodelling in small animal models. Thirty-five relevant in vivo imaging studies are included, of which most made use of magnetic resonance imaging or positron emission tomography, whilst various optical modalities were also utilised. Notably, there is an increasing trend of using multimodal imaging to exploit the most beneficial properties of each imaging technique to elucidate different aspects of neurovascular remodelling. Nevertheless, there is still scope for further utilising noninvasive imaging tools such as contrast agents or radiotracers, which will have the ability to monitor neurovascular changes particularly during restorative therapy. This will facilitate more successful utility of the clinical imaging techniques in the interpretation of neurovascular reorganisation over time.

Functional connectivity with cortical depth assessed by resting state fMRI of subregions of S1 in squirrel monkeys

Whereas resting state blood oxygenation-level dependent (BOLD) functional MRI has been widely used to assess functional connectivity between cortical regions, the laminar specificity of such measures is poorly understood. This study aims to determine: (a) whether the resting state functional connectivity (rsFC) between two functionally related cortical regions varies with cortical depth, (b) the relationship between layer-resolved tactile stimulus-evoked activation pattern and interlayer rsFC pattern between two functionally distinct but related somatosensory areas 3b and 1, and (c) the effects of spatial resolution on rsFC measures. We examined the interlayer rsFC between areas 3b and 1 of squirrel monkeys under anesthesia using tactile stimulus-driven and resting state BOLD acquisitions at submillimeter resolution. Consistent with previous observations in the areas 3b and 1, we detected robust stimulus-evoked BOLD activations with foci were confined mainly to the upper layers (centered at 21% of the cortical depth). By carefully placing seeds in upper, middle, and lower layers of areas 3b and 1, we observed strong rsFC between upper and middle layers of these two areas. The layer-resolved activation patterns in areas 3b and 1 agree with their interlayer rsFC patterns, and are consistent with the known anatomical connections between layers. In summary, using BOLD rsFC pattern, we identified an interlayer interareal microcircuit that shows strong intrinsic functional connections between upper and middle layer areas 3b and 1. RsFC can be used as a robust invasive tool to probe interlayer corticocortical microcircuits.

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.

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.

Imaging at ultrahigh magnetic fields: History, challenges, and solutions

Following early efforts in applying nuclear magnetic resonance (NMR) spectroscopy to study biological processes in intact systems, and particularly since the introduction of 4 T human scanners circa 1990, rapid progress was made in imaging and spectroscopy studies of humans at 4 T and animal models at 9.4 T, leading to the introduction of 7 T and higher magnetic fields for human investigation at about the turn of the century. Work conducted on these platforms has provided numerous technological solutions to challenges posed at these ultrahigh fields, and demonstrated the existence of significant advantages in signal-to-noise ratio and biological information content. Primary difference from lower fields is the deviation from the near field regime at the radiofrequencies (RF) corresponding to hydrogen resonance conditions. At such ultrahigh fields, the RF is characterized by attenuated traveling waves in the human body, which leads to image non-uniformities for a given sample-coil configuration because of destructive and constructive interferences. These non-uniformities were initially considered detrimental to progress of imaging at high field strengths. However, they are advantageous for parallel imaging in signal reception and transmission, two critical technologies that account, to a large extend, for the success of ultrahigh fields. With these technologies and improvements in instrumentation and imaging methods, today ultrahigh fields have provided unprecedented gains in imaging of brain function and anatomy, and started to make inroads into investigation of the human torso and extremities. As extensive as they are, these gains still constitute a prelude to what is to come given the increasingly larger effort committed to ultrahigh field research and development of ever better instrumentation and techniques.

Evaluation of SLIce Dithered Enhanced Resolution Simultaneous MultiSlice (SLIDER-SMS) for human fMRI

High isotropic resolution fMRI is challenging primarily due to long repetition times (TR) and insufficient SNR, especially at lower field strengths. Recently, Simultaneous Multi-Slice (SMS) imaging with blipped-CAIPI has substantially reduced scan time and improved SNR efficiency of fMRI. Similarly, super-resolution techniques utilizing sub- voxel spatial shifts in the slice direction have increased both resolution and SNR efficiency. Here we demonstrate the synergistic combination of SLIce Dithered Enhanced Resolution (SLIDER) and SMS for high-resolution, high-SNR whole brain fMRI in comparison to standard resolution fMRI data as well as high-resolution data. With SLIDER-SMS, high spatial frequency information is recovered (unaliased) even in absence of super-resolution deblurring algorithms. Additionally we find that BOLD CNR (as measured by t-value in a visual checkerboard paradigm) is improved by as much as 100% relative to traditionally acquired high- resolution data. Using this gain in CNR, we are able to obtain unprecedented nominally isotropic resolutions at 3T (0.66 mm) and 7T (0.45 mm).

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

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

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.

Laminar fMRI: What can the time domain tell us?

The rapid developments in functional MRI (fMRI) acquisition methods and hardware technologies in recent years, particularly at high field (≥7 T), have enabled unparalleled visualization of functional detail at a laminar or columnar level, bringing fMRI close to the intrinsic resolution of brain function. These advances highlight the potential of high resolution fMRI to be a valuable tool to study the fundamental processing performed in cortical micro-circuits, and their interactions such as feedforward and feedback processes. Notably, because fMRI measures neuronal activity via hemodynamics, the ultimate resolution it affords depends on the spatial specificity of hemodynamics to neuronal activity at a detailed spatial scale, and by the evolution of this specificity over time. Several laminar (≤1 mm spatial resolution) fMRI studies have examined spatial characteristics of the measured hemodynamic signals across cortical depth, in light of understanding or improving the spatial specificity of laminar fMRI. Few studies have examined temporal features of the hemodynamic response across cortical depth. Temporal features of the hemodynamic response offer an additional means to improve the specificity of fMRI, and could help target neuronal processes and neurovascular coupling relationships across laminae, for example by differences in the onset times of the response across cortical depth. In this review, we discuss factors that affect the timing of neuronal and hemodynamic responses across laminae, touching on the neuronal laminar organization, and focusing on the laminar vascular organization. We provide an overview of hemodynamics across the cortical vascular tree based on optical imaging studies, and review temporal aspects of hemodynamics that have been examined across cortical depth in high spatiotemporal resolution fMRI studies. Last, we discuss the limits and potential of high spatiotemporal resolution fMRI to study laminar neurovascular coupling and neuronal processes.

Foundations of layer-specific fMRI and investigations of neurophysiological activity in the laminarized neocortex and olfactory bulb of animal models

Laminar organization of neuronal circuits is a recurring feature of how the brain processes information. For instance, different layers compartmentalize different cell types, synaptic activities, and have unique intrinsic and extrinsic connections that serve as units for specialized signal processing. Functional MRI is an invaluable tool to investigate laminar processing in the in vivo human brain, but it measures neuronal activity indirectly by way of the hemodynamic response. Therefore, the accuracy of high-resolution laminar fMRI depends on how precisely it can measure localized microvascular changes nearest to the site of evoked activity. To determine the specificity of fMRI responses to the true neurophysiological responses across layers, the flexibility to invasive procedures in animal models has been necessary. In this review, we will examine different fMRI contrasts and their appropriate uses for layer-specific fMRI, and how localized laminar processing was examined in the neocortex and olfactory bulb. Through collective efforts, it was determined that microvessels, including capillaries, are regulated within single layers and that several endogenous and contrast-enhanced fMRI contrast mechanisms can separate these neural-specific vascular changes from the nonspecific, especially cerebral blood volume-weighted fMRI with intravenous contrast agent injection. We will also propose some open questions that are relevant for the successful implementation of layer-specific fMRI and its potential future directions to study laminar processing when combined with optogenetics.

Cortical Representation of Pain and Touch: Evidence from Combined Functional Neuroimaging and Electrophysiology in Non-human Primates

Human functional MRI studies in acute and various chronic pain conditions have revolutionized how we view pain, and have led to a new theory that complex multi-dimensional pain experience (sensory-discriminative, affective/motivational, and cognitive) is represented by concurrent activity in widely-distributed brain regions (termed a network or pain matrix). Despite these breakthrough discoveries, the specific functions proposed for these regions remain elusive, because detailed electrophysiological characterizations of these regions in the primate brain are lacking. To fill in this knowledge gap, we have studied the cortical areas around the central and lateral sulci of the non-human primate brain with combined submillimeter resolution functional imaging (optical imaging and fMRI) and intracranial electrophysiological recording. In this mini-review, I summarize and present data showing that the cortical circuitry engaged in nociceptive processing is much more complex than previously recognized. Electrophysiological evidence supports the engagement of a distinct nociceptive-processing network within SI (i.e., areas 3a, 3b, 1 and 2), SII, and other areas along the lateral sulcus. Deafferentation caused by spinal cord injury profoundly alters the relationships between fMRI and electrophysiological signals. This finding has significant implications for using fMRI to study chronic pain conditions involving deafferentation in humans.

Sensitivity and specificity considerations for fMRI encoding, decoding, and mapping of auditory cortex at ultra-high field

Following rapid technological advances, ultra-high field functional MRI (fMRI) enables exploring correlates of neuronal population activity at an increasing spatial resolution. However, as the fMRI blood-oxygenation-level-dependent (BOLD) contrast is a vascular signal, the spatial specificity of fMRI data is ultimately determined by the characteristics of the underlying vasculature. At 7T, fMRI measurement parameters determine the relative contribution of the macro- and microvasculature to the acquired signal. Here we investigate how these parameters affect relevant high-end fMRI analyses such as encoding, decoding, and submillimeter mapping of voxel preferences in the human auditory cortex. Specifically, we compare a T₂* weighted fMRI dataset, obtained with 2D gradient echo (GE) EPI, to a predominantly T₂ weighted dataset obtained with 3D GRASE. We first investigated the decoding accuracy based on two encoding models that represented different hypotheses about auditory cortical processing. This encoding/decoding analysis profited from the large spatial coverage and sensitivity of the T₂* weighted acquisitions, as evidenced by a significantly higher prediction accuracy in the GE-EPI dataset compared to the 3D GRASE dataset for both encoding models. The main disadvantage of the T₂* weighted GE-EPI dataset for encoding/decoding analyses was that the prediction accuracy exhibited cortical depth dependent vascular biases. However, we propose that the comparison of prediction accuracy across the different encoding models may be used as a post processing technique to salvage the spatial interpretability of the GE-EPI cortical depth-dependent prediction accuracy. Second, we explored the mapping of voxel preferences. Large-scale maps of frequency preference (i.e., tonotopy) were similar across datasets, yet the GE-EPI dataset was preferable due to its larger spatial coverage and sensitivity. However, submillimeter tonotopy maps revealed biases in assigned frequency preference and selectivity for the GE-EPI dataset, but not for the 3D GRASE dataset. Thus, a T₂ weighted acquisition is recommended if high specificity in tonotopic maps is required. In conclusion, different fMRI acquisitions were better suited for different analyses. It is therefore critical that any sequence parameter optimization considers the eventual intended fMRI analyses and the nature of the neuroscience questions being asked.

Adaptive smoothing based on Gaussian processes regression increases the sensitivity and specificity of fMRI data

Temporal and spatial filtering of fMRI data is often used to improve statistical power. However, conventional methods, such as smoothing with fixed-width Gaussian filters, remove fine-scale structure in the data, necessitating a tradeoff between sensitivity and specificity. Specifically, smoothing may increase sensitivity (reduce noise and increase statistical power) but at the cost loss of specificity in that fine-scale structure in neural activity patterns is lost. Here, we propose an alternative smoothing method based on Gaussian processes (GP) regression for single subjects fMRI experiments. This method adapts the level of smoothing on a voxel by voxel basis according to the characteristics of the local neural activity patterns. GP-based fMRI analysis has been heretofore impractical owing to computational demands. Here, we demonstrate a new implementation of GP that makes it possible to handle the massive data dimensionality of the typical fMRI experiment. We demonstrate how GP can be used as a drop-in replacement to conventional preprocessing steps for temporal and spatial smoothing in a standard fMRI pipeline. We present simulated and experimental results that show the increased sensitivity and specificity compared to conventional smoothing strategies. Hum Brain Mapp 38:1438-1459, 2017. © 2016 Wiley Periodicals, Inc.

Tradeoffs in pushing the spatial resolution of fMRI for the 7T Human Connectome Project

Whole-brain functional magnetic resonance imaging (fMRI), in conjunction with multiband acceleration, has played an important role in mapping the functional connectivity throughout the entire brain with both high temporal and spatial resolution. Ultrahigh magnetic field strengths (7T and above) allow functional imaging with even higher functional contrast-to-noise ratios for improved spatial resolution and specificity compared to traditional field strengths (1.5T and 3T). High-resolution 7T fMRI, however, has primarily been constrained to smaller brain regions given the amount of time it takes to acquire the number of slices necessary for high resolution whole brain imaging. Here we evaluate a range of whole-brain high-resolution resting state fMRI protocols (0.9, 1.25, 1.5, 1.6 and 2mm isotropic voxels) at 7T, obtained with both in-plane and slice acceleration parallel imaging techniques to maintain the temporal resolution and brain coverage typically acquired at 3T. Using the processing pipeline developed by the Human Connectome Project, we demonstrate that high resolution images acquired at 7T provide increased functional contrast to noise ratios with significantly less partial volume effects and more distinct spatial features, potentially allowing for robust individual subject parcellations and descriptions of fine-scaled patterns, such as visuotopic organization.

Intracortical depth analyses of frequency-sensitive regions of human auditory cortex using 7TfMRI

Despite recent advances in auditory neuroscience, the exact functional organization of human auditory cortex (AC) has been difficult to investigate. Here, using reversals of tonotopic gradients as the test case, we examined whether human ACs can be more precisely mapped by avoiding signals caused by large draining vessels near the pial surface, which bias blood-oxygen level dependent (BOLD) signals away from the actual sites of neuronal activity. Using ultra-high field (7T) fMRI and cortical depth analysis techniques previously applied in visual cortices, we sampled 1mm isotropic voxels from different depths of AC during narrow-band sound stimulation with biologically relevant temporal patterns. At the group level, analyses that considered voxels from all cortical depths, but excluded those intersecting the pial surface, showed (a) the greatest statistical sensitivity in contrasts between activations to high vs. low frequency sounds and (b) the highest inter-subject consistency of phase-encoded continuous tonotopy mapping. Analyses based solely on voxels intersecting the pial surface produced the least consistent group results, even when compared to analyses based solely on voxels intersecting the white-matter surface where both signal strength and within-subject statistical power are weakest. However, no evidence was found for reduced within-subject reliability in analyses considering the pial voxels only. Our group results could, thus, reflect improved inter-subject correspondence of high and low frequency gradients after the signals from voxels near the pial surface are excluded. Using tonotopy analyses as the test case, our results demonstrate that when the major physiological and anatomical biases imparted by the vasculature are controlled, functional mapping of human ACs becomes more consistent from subject to subject than previously thought.

Systematic variation of population receptive field properties across cortical depth in human visual cortex

Receptive fields (RFs) in visual cortex are organized in antagonistic, center-surround, configurations. RF properties change systematically across eccentricity and between visual field maps. However, it is unknown how center-surround configurations are organized in human visual cortex across lamina. We use sub-millimeter resolution functional MRI at 7Tesla and population receptive field (pRF) modeling to investigate the pRF properties in primary visual cortex (V1) across cortical depth. pRF size varies according to a U-shaped function, indicating smaller pRF center size in the middle compared to superficial and deeper intra-cortical portions of V1, consistent with non-human primate neurophysiological measurements. Moreover, a similar U-shaped function is also observed for pRF surround size. However, pRF center-surround ratio remains constant across cortical depth. Simulations suggest that this pattern of results can be directly linked to the flow of signals across cortical depth, with the visual input reaching the middle of cortical depth and then spreading towards superficial and deeper layers of V1. Conversely, blood-oxygenation-level-dependent (BOLD) signal amplitude increases monotonically towards the pial surface, in line with the known vascular organization across cortical depth. Independent estimates of the haemodynamic response function (HRF) across cortical depth show that the center-surround pRF size estimates across cortical depth cannot be explained by variations in the full-width half maximum (FWHM) of the HRF.

fMRI at High Spatial Resolution: Implications for BOLD-Models

As high-resolution functional magnetic resonance imaging (fMRI) and fMRI of cortical layers become more widely used, the question how well high-resolution fMRI signals reflect the underlying neural processing, and how to interpret laminar fMRI data becomes more and more relevant. High-resolution fMRI has shown laminar differences in cerebral blood flow (CBF), volume (CBV), and neurovascular coupling. Features and processes that were previously lumped into a single voxel become spatially distinct at high resolution. These features can be vascular compartments such as veins, arteries, and capillaries, or cortical layers and columns, which can have differences in metabolism. Mesoscopic models of the blood oxygenation level dependent (BOLD) response therefore need to be expanded, for instance, to incorporate laminar differences in the coupling between neural activity, metabolism and the hemodynamic response. Here we discuss biological and methodological factors that affect the modeling and interpretation of high-resolution fMRI data. We also illustrate with examples from neuropharmacology and the negative BOLD response how combining BOLD with CBF- and CBV-based fMRI methods can provide additional information about neurovascular coupling, and can aid modeling and interpretation of high-resolution fMRI.