Frequency preference and attention effects across cortical depths in the human primary auditory cortex

Systematic changes in neural selectivity reflect the acquired salience of category-diagnostic dimensions

Humans and other animals develop remarkable behavioral specializations for identifying, differentiating, and acting on classes of ecologically important signals. Ultimately, this expertise is flexible enough to support diverse perceptual judgments: a voice, for example, simultaneously conveys what a talker says as well as myriad cues about her identity and state. Mature perception across complex signals thus involves both discovering and learning regularities that best inform diverse perceptual judgments, and weighting this information flexibly as task demands change. Here, we test whether this flexibility may involve endogenous attentional gain to task-relevant dimensions. We use two prospective auditory category learning tasks to relate a complex, entirely novel soundscape to four classes of "alien identity" and two classes of "alien size." Identity, but not size, categorization requires discovery and learning of patterned acoustic input situated in one of two simultaneous, frequency- delimited bands. This allows us to capitalize on the coarsely segregated frequency-band- specific channels of auditory tonotopic maps using fMRI to ask whether category-relevant perceptual information is prioritized relative to simultaneous, uninformative information. Among participants expert at alien identity categorization, we observe prioritization of the diagnostic frequency band that persists even when the diagnostic information becomes irrelevant in the size categorization task. Tellingly, the neural selectivity evoked implicitly in categorization aligns closely with activation driven by explicit, sustained selective attention to other sounds presented in the same frequency band. Additionally, we observe fingerprints of individual differences in the learning trajectories taken to achieve expert-level categorization in patterns of neural activity associated with the diagnostic dimension. In all, this indicates that acquiring categories can drive the emergence of acquired attentional salience to dimensions of acoustic input.

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

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

Syllable processing is organized in discrete subregions of the human superior temporal gyrus

Modular organization at approximately 1 mm scale could be fundamental to cortical processing, but its presence in human association cortex is unknown. Using custom-built, high-density electrode arrays placed on the cortical surface of 7 patients undergoing awake craniotomy for tumor excision, we investigated receptive speech processing in the left (dominant) human posterior superior temporal gyrus. Responses to consonant-vowel syllables and noise-vocoded controls recorded with 1,024 channel micro-grids at 200 μm pitch demonstrated roughly circular domains approximately 1.7 mm in diameter, with sharp boundaries observed in 128 channel linear arrays at 50 μm pitch, possibly consistent with a columnar organization. Peak latencies to syllables in different modules were bimodally distributed centered at 252 and 386 ms. Adjacent modules were sharply delineated from each other by their distinct time courses and stimulus selectivity. We suggest that receptive language cortex may be organized in discrete processing modules.

Attention-Driven Modulation of Auditory Cortex Activity during Selective Listening in a Multispeaker Setting

Real-world listening settings often consist of multiple concurrent sound streams. To limit perceptual interference during selective listening, the auditory system segregates and filters the relevant sensory input. Previous work provided evidence that the auditory cortex is critically involved in this process and selectively gates attended input toward subsequent processing stages. We studied at which level of auditory cortex processing this filtering of attended information occurs using functional magnetic resonance imaging (fMRI) and a naturalistic selective listening task. Forty-five human listeners (of either sex) attended to one of two continuous speech streams, presented either concurrently or in isolation. Functional data were analyzed using an inter-subject analysis to assess stimulus-specific components of ongoing auditory cortex activity. Our results suggest that stimulus-related activity in the primary auditory cortex and the adjacent planum temporale are hardly affected by attention, whereas brain responses at higher stages of the auditory cortex processing hierarchy become progressively more selective for the attended input. Consistent with these findings, a complementary analysis of stimulus-driven functional connectivity further demonstrated that information on the to-be-ignored speech stream is shared between the primary auditory cortex and the planum temporale but largely fails to reach higher processing stages. Our findings suggest that the neural processing of ignored speech cannot be effectively suppressed at the level of early cortical processing of acoustic features but is gradually attenuated once the competing speech streams are fully segregated.

Intracranial Mapping of Response Latencies and Task Effects for Spoken Syllable Processing in the Human Brain

Prior lesion, noninvasive-imaging, and intracranial-electroencephalography (iEEG) studies have documented hierarchical, parallel, and distributed characteristics of human speech processing. Yet, there have not been direct, intracranial observations of the latency with which regions outside the temporal lobe respond to speech, or how these responses are impacted by task demands. We leveraged human intracranial recordings via stereo-EEG to measure responses from diverse forebrain sites during (i) passive listening to /bi/ and /pi/ syllables, and (ii) active listening requiring /bi/-versus-/pi/ categorization. We find that neural response latency increases from a few tens of ms in Heschl's gyrus (HG) to several tens of ms in superior temporal gyrus (STG), superior temporal sulcus (STS), and early parietal areas, and hundreds of ms in later parietal areas, insula, frontal cortex, hippocampus, and amygdala. These data also suggest parallel flow of speech information dorsally and ventrally, from HG to parietal areas and from HG to STG and STS, respectively. Latency data also reveal areas in parietal cortex, frontal cortex, hippocampus, and amygdala that are not responsive to the stimuli during passive listening but are responsive during categorization. Furthermore, multiple regions-spanning auditory, parietal, frontal, and insular cortices, and hippocampus and amygdala-show greater neural response amplitudes during active versus passive listening (a task-related effect). Overall, these results are consistent with hierarchical processing of speech at a macro level and parallel streams of information flow in temporal and parietal regions. These data also reveal regions where the speech code is stimulus-faithful and those that encode task-relevant representations.

Large-scale single-neuron speech sound encoding across the depth of human cortex

Understanding the neural basis of speech perception requires that we study the human brain both at the scale of the fundamental computational unit of neurons and in their organization across the depth of cortex. Here we used high-density Neuropixels arrays1-3 to record from 685 neurons across cortical layers at nine sites in a high-level auditory region that is critical for speech, the superior temporal gyrus4,5, while participants listened to spoken sentences. Single neurons encoded a wide range of speech sound cues, including features of consonants and vowels, relative vocal pitch, onsets, amplitude envelope and sequence statistics. Neurons at each cross-laminar recording exhibited dominant tuning to a primary speech feature while also containing a substantial proportion of neurons that encoded other features contributing to heterogeneous selectivity. Spatially, neurons at similar cortical depths tended to encode similar speech features. Activity across all cortical layers was predictive of high-frequency field potentials (electrocorticography), providing a neuronal origin for macroelectrode recordings from the cortical surface. Together, these results establish single-neuron tuning across the cortical laminae as an important dimension of speech encoding in human superior temporal gyrus.

Evaluating the effect of denoising submillimeter auditory fMRI data with NORDIC

Functional magnetic resonance imaging (fMRI) has emerged as an essential tool for exploring human brain function. Submillimeter fMRI, in particular, has emerged as a tool to study mesoscopic computations. The inherently low signal-to-noise ratio (SNR) at submillimeter resolutions warrants the use of denoising approaches tailored at reducing thermal noise - the dominant contributing noise component in high resolution fMRI. NORDIC PCA is one of such approaches, and has been benchmarked against other approaches in several applications. Here, we investigate the effects that two versions of NORDIC denoising have on auditory submillimeter data. As investigating auditory functional responses poses unique challenges, we anticipated that the benefit of this technique would be especially pronounced. Our results show that NORDIC denoising improves the detection sensitivity and the reliability of estimates in submillimeter auditory fMRI data. These effects can be explained by the reduction of the noise-induced signal variability. However, we also observed a reduction in the average response amplitude (percent signal), which may suggest that a small amount of signal was also removed. We conclude that, while evaluating the effects of the signal reduction induced by NORDIC may be necessary for each application, using NORDIC in high resolution auditory fMRI studies may be advantageous because of the large reduction in variability of the estimated responses.

Modular Phoneme Processing in Human Superior Temporal Gyrus

Modular organization is fundamental to cortical processing, but its presence is human association cortex is unknown. We characterized phoneme processing with 128-1024 channel micro-arrays at 50-200µm pitch on superior temporal gyrus of 7 patients. High gamma responses were highly correlated within ~1.7mm diameter modules, sharply delineated from adjacent modules with distinct time-courses and phoneme-selectivity. We suggest that receptive language cortex may be organized in discrete processing modules.

Electrophysiological effects of smartphone notifications on cognitive control following a brief mindfulness induction

Smartphone use is nearly ubiquitous, with 93% of adults among economically developed countries, including the United States, Canada, Israel, and South Korea owning a smartphone (Taylor & Silver, 2019). Multiple studies have demonstrated the distracting effects of smartphone notifications on behavioral measures of cognition. Fewer studies have examined the effects of notifications on neural activity underlying higher-level cognitive processes or behavioral inductions to reduce smartphone-related distraction. Using EEG spectral frequency power densities, we assessed the effects of smartphone notifications (vs. control trials) on engagement of attentional shifting processes involved in cognitive control during a Navon Letter visual oddball task. Participants were randomly assigned to a brief mindfulness induction (N = 44) or a neutral narration control condition (N = 43). Overall, participants had lower theta-band power, but higher alpha- and beta-band power densities on target letter trials preceded by smartphone notifications. Additionally, participants in the mindfulness (vs. control) condition had a larger attention shifting oddball assessed via theta power density and theta/beta ratio (TBR) values-reflecting increased engagement of cognitive control-particularly on smartphone notification (vs. control) trials. Altogether, these results provide evidence supporting the idea that smartphone notifications can decrease activity of neural correlates of cognitive control, and offer the promise of a brief mindfulness induction to buffer against the effects of smartphone notifications on cognitive control. The findings indicate a need for further research on mindfulness inductiosn as a means to reduce potential distraction caused by smartphones.

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

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

Functional neuroimaging as a catalyst for integrated neuroscience

Functional magnetic resonance imaging (fMRI) enables non-invasive access to the awake, behaving human brain. By tracking whole-brain signals across a diverse range of cognitive and behavioural states or mapping differences associated with specific traits or clinical conditions, fMRI has advanced our understanding of brain function and its links to both normal and atypical behaviour. Despite this headway, progress in human cognitive neuroscience that uses fMRI has been relatively isolated from rapid advances in other subdomains of neuroscience, which themselves are also somewhat siloed from one another. In this Perspective, we argue that fMRI is well-placed to integrate the diverse subfields of systems, cognitive, computational and clinical neuroscience. We first summarize the strengths and weaknesses of fMRI as an imaging tool, then highlight examples of studies that have successfully used fMRI in each subdomain of neuroscience. We then provide a roadmap for the future advances that will be needed to realize this integrative vision. In this way, we hope to demonstrate how fMRI can help usher in a new era of interdisciplinary coherence in neuroscience.

Ultra-High Field Imaging of Human Visual Cognition

Functional magnetic resonance imaging (fMRI), the key methodology for mapping the functions of the human brain in a noninvasive manner, is limited by low temporal and spatial resolution. Recent advances in ultra-high field (UHF) fMRI provide a mesoscopic (i.e., submillimeter resolution) tool that allows us to probe laminar and columnar circuits, distinguish bottom-up versus top-down pathways, and map small subcortical areas. We review recent work demonstrating that UHF fMRI provides a robust methodology for imaging the brain across cortical depths and columns that provides insights into the brain's organization and functions at unprecedented spatial resolution, advancing our understanding of the fine-scale computations and interareal communication that support visual cognition.

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

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

Cortical depth-dependent human fMRI of resting-state networks using EPIK

Introduction

Recent laminar-fMRI studies have substantially improved understanding of the evoked cortical responses in multiple sub-systems; in contrast, the laminar component of resting-state networks spread over the whole brain has been less studied due to technical limitations. Animal research strongly suggests that the supragranular layers of the cortex play a critical role in maintaining communication within the default mode network (DMN); however, whether this is true in this and other human cortical networks remains unclear.

Methods

Here, we used EPIK, which offers unprecedented coverage at sub-millimeter resolution, to investigate cortical broad resting-state dynamics with depth specificity in healthy volunteers.

Results

Our results suggest that human DMN connectivity is primarily supported by intermediate and superficial layers of the cortex, and furthermore, the preferred cortical depth used for communication can vary from one network to another. In addition, the laminar connectivity profile of some networks showed a tendency to change upon engagement in a motor task. In line with these connectivity changes, we observed that the amplitude of the low-frequency-fluctuations (ALFF), as well as the regional homogeneity (ReHo), exhibited a different laminar slope when subjects were either performing a task or were in a resting state (less variation among laminae, i.e., lower slope, during task performance compared to rest).

Discussion

The identification of varied laminar profiles concerning network connectivity, ALFF, and ReHo, observed across two brain states (task vs. rest) has major implications for the characterization of network-related diseases and suggests the potential diagnostic value of laminar fMRI in psychiatric disorders, e.g., to differentiate the cortical dynamics associated with disease stages linked, or not linked, to behavioral changes. The evaluation of laminar-fMRI across the brain encompasses computational challenges; nonetheless, it enables the investigation of a new dimension of the human neocortex, which may be key to understanding neurological disorders from a novel perspective.

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

INTRODUCTION

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

METHODS

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

RESULTS AND CONCLUSION

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

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

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

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

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

Layer-specific correlates of detected and undetected auditory targets during attention

In everyday life, the processing of acoustic information allows us to react to subtle changes in the auditory scene. Yet even when closely attending to sounds in the context of a task, we occasionally miss task-relevant features. The neural computations that underlie our ability to detect behavioral relevant sound changes are thought to be grounded in both feedforward and feedback processes within the auditory hierarchy. Here, we assessed the role of feedforward and feedback contributions in primary and non-primary auditory areas during behavioral detection of target sounds using submillimeter spatial resolution functional magnetic resonance imaging (fMRI) at high-fields (7 T) in humans. We demonstrate that the successful detection of subtle temporal shifts in target sounds leads to a selective increase of activation in superficial layers of primary auditory cortex (PAC). These results indicate that feedback signals reaching as far back as PAC may be relevant to the detection of targets in the auditory scene.

Selective attention sharpens population receptive fields in human auditory cortex

Selective attention enables the preferential processing of relevant stimulus aspects. Invasive animal studies have shown that attending a sound feature rapidly modifies neuronal tuning throughout the auditory cortex. Human neuroimaging studies have reported enhanced auditory cortical responses with selective attention. To date, it remains unclear how the results obtained with functional magnetic resonance imaging (fMRI) in humans relate to the electrophysiological findings in animal models. Here we aim to narrow the gap between animal and human research by combining a selective attention task similar in design to those used in animal electrophysiology with high spatial resolution ultra-high field fMRI at 7 Tesla. Specifically, human participants perform a detection task, whereas the probability of target occurrence varies with sound frequency. Contrary to previous fMRI studies, we show that selective attention resulted in population receptive field sharpening, and consequently reduced responses, at the attended sound frequencies. The difference between our results to those of previous fMRI studies supports the notion that the influence of selective attention on auditory cortex is diverse and may depend on context, stimulus, and task.

What has vision science taught us about functional MRI?

In the domain of human neuroimaging, much attention has been paid to the question of whether and how the development of functional magnetic resonance imaging (fMRI) has advanced our scientific knowledge of the human brain. However, the opposite question is also important; how has our knowledge of the brain advanced our understanding of fMRI? Here, we discuss how and why scientific knowledge about the human and animal visual system has been used to answer fundamental questions about fMRI as a brain measurement tool and how these answers have contributed to scientific discoveries beyond vision science.

The promise of layer-specific neuroimaging for testing predictive coding theories of psychosis

Predictive coding potentially provides an explanatory model for understanding the neurocognitive mechanisms of psychosis. It proposes that cognitive processes, such as perception and inference, are implemented by a hierarchical system, with the influence of each level being a function of the estimated precision of beliefs at that level. However, predictive coding models of psychosis are insufficiently constrained-any phenomenon can be explained in multiple ways by postulating different changes to precision at different levels of processing. One reason for the lack of constraint in these models is that the core processes are thought to be implemented by the function of specific cortical layers, and the technology to measure layer specific neural activity in humans has until recently been lacking. As a result, our ability to constrain the models with empirical data has been limited. In this review we provide a brief overview of predictive processing models of psychosis and then describe the potential for newly developed, layer specific neuroimaging techniques to test and thus constrain these models. We conclude by discussing the most promising avenues for this research as well as the technical and conceptual challenges which may limit its application.

Cytoarchitectonic Maps of the Human Metathalamus in 3D Space

The human metathalamus plays an important role in processing visual and auditory information. Understanding its layers and subdivisions is important to gain insights in its function as a subcortical relay station and involvement in various pathologies. Yet, detailed histological references of the microanatomy in 3D space are still missing. We therefore aim at providing cytoarchitectonic maps of the medial geniculate body (MGB) and its subdivisions in the BigBrain – a high-resolution 3D-reconstructed histological model of the human brain, as well as probabilistic cytoarchitectonic maps of the MGB and lateral geniculate body (LGB). Therefore, histological sections of ten postmortem brains were studied. Three MGB subdivisions (MGBv, MGBd, MGBm) were identified on every 5th BigBrain section, and a deep-learning based tool was applied to map them on every remaining section. The maps were 3D-reconstructed to show the shape and extent of the MGB and its subdivisions with cellular precision. The LGB and MGB were additionally identified in nine other postmortem brains. Probabilistic cytoarchitectonic maps in the MNI “Colin27” and MNI ICBM152 reference spaces were computed which reveal an overall low interindividual variability in topography and extent. The probabilistic maps were included into the Julich-Brain atlas, and are freely available. They can be linked to other 3D data of human brain organization and serve as an anatomical reference for diagnostic, prognostic and therapeutic neuroimaging studies of healthy brains and patients. Furthermore, the high-resolution MGB BigBrain maps provide a basis for data integration, brain modeling and simulation to bridge the larger scale involvement of thalamocortical and local subcortical circuits.

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.

Imaging faster neural dynamics with fast fMRI: A need for updated models of the hemodynamic response

Fast fMRI enables the detection of neural dynamics over timescales of hundreds of milliseconds, suggesting it may provide a new avenue for studying subsecond neural processes in the human brain. The magnitudes of these fast fMRI dynamics are far greater than predicted by canonical models of the hemodynamic response. Several studies have established nonlinear properties of the hemodynamic response that have significant implications for fast fMRI. We first review nonlinear properties of the hemodynamic response function that may underlie fast fMRI signals. We then illustrate the breakdown of canonical hemodynamic response models in the context of fast neural dynamics. We will then argue that the canonical hemodynamic response function is not likely to reflect the BOLD response to neuronal activity driven by sparse or naturalistic stimuli or perhaps to spontaneous neuronal fluctuations in the resting state. These properties suggest that fast fMRI is capable of tracking surprisingly fast neuronal dynamics, and we discuss the neuroscientific questions that could be addressed using this approach.

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

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

Ultra-High-Field Neuroimaging Reveals Fine-Scale Processing for 3D Perception

Binocular disparity provides critical information about three-dimensional (3D) structures to support perception and action. In the past decade significant progress has been made in uncovering human brain areas engaged in the processing of binocular disparity signals. Yet, the fine-scale brain processing underlying 3D perception remains unknown. Here, we use ultra-high-field (7T) functional imaging at submillimeter resolution to examine fine-scale BOLD fMRI signals involved in 3D perception. In particular, we sought to interrogate the local circuitry involved in disparity processing by sampling fMRI responses at different positions relative to the cortical surface (i.e., across cortical depths corresponding to layers). We tested for representations related to 3D perception by presenting participants (male and female, N = 8) with stimuli that enable stable stereoscopic perception [i.e., correlated random dot stereograms (RDS)] versus those that do not (i.e., anticorrelated RDS). Using multivoxel pattern analysis (MVPA), we demonstrate cortical depth-specific representations in areas V3A and V7 as indicated by stronger pattern responses for correlated than for anticorrelated stimuli in upper rather than deeper layers. Examining informational connectivity, we find higher feedforward layer-to-layer connectivity for correlated than anticorrelated stimuli between V3A and V7. Further, we observe disparity-specific feedback from V3A to V1 and from V7 to V3A. Our findings provide evidence for the role of V3A as a key nexus for disparity processing, which is implicated in feedforward and feedback signals related to the perceptual estimation of 3D structures.SIGNIFICANCE STATEMENT Binocular vision plays a significant role in supporting our interactions with the surrounding environment. The fine-scale neural mechanisms that underlie the brain's skill in extracting 3D structures from binocular signals are poorly understood. Here, we capitalize on recent advances in ultra-high-field functional imaging to interrogate human brain circuits involved in 3D perception at submillimeter resolution. We provide evidence for the role of area V3A as a key nexus for disparity processing, which is implicated in feedforward and feedback signals related to the perceptual estimation of 3D structures from binocular signals. These fine-scale measurements help bridge the gap between animal neurophysiology and human fMRI studies investigating cross-scale circuits, from micro circuits to global brain networks for 3D perception.

Structural Covariance of Depth-Dependent Intracortical Myelination in the Human Brain and Its Application to Drug-Naïve Schizophrenia: A T1w/T2w MRI Study

Aberrations in intracortical myelination are increasingly being considered as a cardinal feature in the pathophysiology of schizophrenia. We investigated the network-level distribution of intracortical myelination across various cortex depths. We enrolled 126 healthy subjects and 106 first-episode drug-naïve schizophrenia patients. We used T1w/T2w ratio as a proxy of intracortical myelination, parcellated cortex into several equivolumetric surfaces based on cortical depths and mapped T1w/T2w ratios to each surface. Non-negative matrix factorization was used to generate depth-dependent structural covariance networks (dSCNs) of intracortical myelination from 2 healthy controls datasets-one from our study and another from 100-unrelated dataset of the Human Connectome Project. For patient versus control comparisons, partial least squares approach was used; we also related myelination to clinical features of schizophrenia. We found that dSCNs were highly reproducible in 2 independent samples. Network-level myelination was reduced in prefrontal and cingulate cortex and increased in perisylvian cortex in schizophrenia. The abnormal network-level myelination had a canonical correlation with symptom burden in schizophrenia. Moreover, myelination of prefrontal cortex correlated with duration of untreated psychosis. In conclusion, we offer a feasible and sensitive framework to study depth-dependent myelination and its relationship with clinical features.

Naturalistic Stimuli: A Paradigm for Multi-Scale Functional Characterization of the Human Brain

Movies, audio stories, and virtual reality are increasingly used as stimuli for functional brain imaging. Such naturalistic paradigms are in sharp contrast to the tradition of experimental reductionism in neuroscience research. Being complex, dynamic, and diverse, naturalistic stimuli set up a more ecologically relevant condition and induce highly reproducible brain responses across a wide range of spatiotemporal scales. Here, we review recent technical advances and scientific findings on imaging the brain under naturalistic stimuli. Then we elaborate on the premise of using naturalistic paradigms for multi-scale, multi-modal, and high-throughput functional characterization of the human brain. We further highlight the growing potential of using deep learning models to infer neural information processing from brain responses to naturalistic stimuli. Lastly, we advocate large-scale collaborations to combine brain imaging and recording data across experiments, subjects, and labs that use the same set of naturalistic stimuli.

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.

Linking cortical circuit models to human cognition with laminar fMRI

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

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.

Laminar fMRI using T2-prepared multi-echo FLASH

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

The nonhuman primate neuroimaging and neuroanatomy project

Multi-modal neuroimaging projects such as the Human Connectome Project (HCP) and UK Biobank are advancing our understanding of human brain architecture, function, connectivity, and their variability across individuals using high-quality non-invasive data from many subjects. Such efforts depend upon the accuracy of non-invasive brain imaging measures. However, 'ground truth' validation of connectivity using invasive tracers is not feasible in humans. Studies using nonhuman primates (NHPs) enable comparisons between invasive and non-invasive measures, including exploration of how "functional connectivity" from fMRI and "tractographic connectivity" from diffusion MRI compare with long-distance connections measured using tract tracing. Our NonHuman Primate Neuroimaging & Neuroanatomy Project (NHP_NNP) is an international effort (6 laboratories in 5 countries) to: (i) acquire and analyze high-quality multi-modal brain imaging data of macaque and marmoset monkeys using protocols and methods adapted from the HCP; (ii) acquire quantitative invasive tract-tracing data for cortical and subcortical projections to cortical areas; and (iii) map the distributions of different brain cell types with immunocytochemical stains to better define brain areal boundaries. We are acquiring high-resolution structural, functional, and diffusion MRI data together with behavioral measures from over 100 individual macaques and marmosets in order to generate non-invasive measures of brain architecture such as myelin and cortical thickness maps, as well as functional and diffusion tractography-based connectomes. We are using classical and next-generation anatomical tracers to generate quantitative connectivity maps based on brain-wide counting of labeled cortical and subcortical neurons, providing ground truth measures of connectivity. Advanced statistical modeling techniques address the consistency of both kinds of data across individuals, allowing comparison of tracer-based and non-invasive MRI-based connectivity measures. We aim to develop improved cortical and subcortical areal atlases by combining histological and imaging methods. Finally, we are collecting genetic and sociality-associated behavioral data in all animals in an effort to understand how genetic variation shapes the connectome and behavior.

Laminar Profile of Auditory Steady-State Response in the Auditory Cortex of Awake Mice

Objective

Auditory steady-state response (ASSR) is a gamma oscillation evoked by periodic auditory stimuli, which is commonly used in clinical electroencephalographic examination to evaluate the neurological functions. Though it has been suggested that auditory cortex is the origin of ASSR, how the laminar architecture of the neocortex contributes to the ASSR recorded from the brain surface remains unclear.

Methods

We used a 16-channel silicon probe to record the local field potential and the single-unit spike activity in the different layers of the auditory cortex of unanesthetized mice. Click-trains with a repetition rate at 40-Hz were present as sound stimuli to evoke ASSR.

Results

We found that the LFPs of all cortical layers showed a stable ASSR synchronizing to the 40-Hz click stimuli, while the ASSR was strongest in the granular (thalamorecipient) layer. Furthermore, time-frequency analyses also revealed the strongest coherence between the signals recorded from the granular layer and pial surface.

Conclusion

Our results reveal that the 40-Hz ASSR primarily shows the evoked gamma oscillation of thalamorecipient layers in the neocortex, and that the ASSR may be a biomarker to detect the cognitive deficits associated with impaired thalamo-cortical connection.

Clarifying the role of higher-level cortices in resolving perceptual ambiguity using ultra high field fMRI

The brain is organized into distinct, flexible networks. Within these networks, cognitive variables such as attention can modulate sensory representations in accordance with moment-to-moment behavioral requirements. These modulations can be studied by varying task demands; however, the tasks employed are often incongruent with the postulated functions of a sensory system, limiting the characterization of the system in relation to natural behaviors. Here we combine domain-specific task manipulations and ultra-high field fMRI to study the nature of top-down modulations. We exploited faces, a visual category underpinned by a complex cortical network, and instructed participants to perform either a stimulus-relevant/domain-specific or a stimulus-irrelevant task in the scanner. We found that 1. perceptual ambiguity (i.e. difficulty of achieving a stable percept) is encoded in top-down modulations from higher-level cortices; 2. the right inferior-temporal lobe is active under challenging conditions and uniquely encodes trial-by-trial variability in face perception.

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 field fMRI reveals origins of feedforward and feedback activity within laminae of human ocular dominance columns

Ultra-high field MRI can functionally image the cerebral cortex of human subjects at the submillimeter scale of cortical columns and laminae. Here, we investigate both in concert, by imaging ocular dominance columns (ODCs) in primary visual cortex (V1) across different cortical depths. We ensured that putative ODC patterns in V1 (a) are stable across runs, sessions, and scanners located in different continents, (b) have a width (~1.3 mm) expected from post-mortem and animal work and (c) are absent at the retinotopic location of the blind spot. We then dissociated the effects of bottom-up thalamo-cortical input and attentional feedback processes on activity in V1 across cortical depth. Importantly, the separation of bottom-up information flows into ODCs allowed us to validly compare attentional conditions while keeping the stimulus identical throughout the experiment. We find that, when correcting for draining vein effects and using both model-based and model-free approaches, the effect of monocular stimulation is largest at deep and middle cortical depths. Conversely, spatial attention influences BOLD activity exclusively near the pial surface. Our findings show that simultaneous interrogation of columnar and laminar dimensions of the cortical fold can dissociate thalamocortical inputs from top-down processing, and allow the investigation of their interactions without any stimulus manipulation.

Prior expectations evoke stimulus-specific activity in the deep layers of the primary visual cortex

The way we perceive the world is strongly influenced by our expectations. In line with this, much recent research has revealed that prior expectations strongly modulate sensory processing. However, the neural circuitry through which the brain integrates external sensory inputs with internal expectation signals remains unknown. In order to understand the computational architecture of the cortex, we need to investigate the way these signals flow through the cortical layers. This is crucial because the different cortical layers have distinct intra- and interregional connectivity patterns, and therefore determining which layers are involved in a cortical computation can inform us on the sources and targets of these signals. Here, we used ultra-high field (7T) functional magnetic resonance imaging (fMRI) to reveal that prior expectations evoke stimulus-specific activity selectively in the deep layers of the primary visual cortex (V1). These findings are in line with predictive processing theories proposing that neurons in the deep cortical layers represent perceptual hypotheses and thereby shed light on the computational architecture of cortex.

The way we perceive the world is strongly influenced by our expectations, but the neural circuitry through which the brain achieves this remains unknown. A study using ultra-high field fMRI reveals that prior expectations evoke stimulus-specific signals in the deep layers of the primary visual cortex.

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

Fine-scale computations for adaptive processing in the human brain

Adapting to the environment statistics by reducing brain responses to repetitive sensory information is key for efficient information processing. Yet, the fine-scale computations that support this adaptive processing in the human brain remain largely unknown. Here, we capitalise on the sub-millimetre resolution of ultra-high field imaging to examine functional magnetic resonance imaging signals across cortical depth and discern competing hypotheses about the brain mechanisms (feedforward vs. feedback) that mediate adaptive processing. We demonstrate layer-specific suppressive processing within visual cortex, as indicated by stronger BOLD decrease in superficial and middle than deeper layers for gratings that were repeatedly presented at the same orientation. Further, we show altered functional connectivity for adaptation: enhanced feedforward connectivity from V1 to higher visual areas, short-range feedback connectivity between V1 and V2, and long-range feedback occipito-parietal connectivity. Our findings provide evidence for a circuit of local recurrent and feedback interactions that mediate rapid brain plasticity for adaptive information processing.

Quantitative Evaluations of Geometrical Distortion Corrections in Cortical Surface-Based Analysis of High-Resolution Functional MRI Data at 7T

Background

Although 7T functional MRI (fMRI) provides better signal‐to‐noise ratio and higher spatial resolution than 3T fMRI, geometric distortions become more challenging because fMRI is more susceptible to distortions than structural MRI. Accurate alignment of 7T fMRI to structural MRI data is critical for precise cortical surface‐based analysis.

Purpose

To quantify the effectiveness of distortion corrections of 7T fMRI data.

Study Type

Prospective.

Subjects

Fifteen healthy individuals aged 19–26 years (mean: 21.9 years).

Field Strength/Sequence

Multiband gradient‐echo echo‐planar imaging sequence at 7T; 3D T₁/T₂‐weighted sequences (magnetization prepared rapid acquisition with gradient echo [MPRAGE] and sampling perfection with application optimized contrast using different flip angle evolution [SPACE]) at 3T.

Assessment

fMRI data at 7T were registered to cortical surfaces reconstructed from 3T structural data acquired in the same subjects. Distortions induced by B₀ inhomogeneity and gradient nonlinearity (B₀ and gradient distortions) were evaluated as cortical fallout (misregistration of noncortical areas) and displacement (misregistration along gray matter).

Statistical Tests

Repeated measures analyses of variance with post‐hoc t‐tests with Bonferroni correction.

Results

The accuracy of fully corrected fMRI images based on the intensity distribution was 89.2%. Without any corrections, 9.7% of vertices in the whole surfaces were fallout and the average displacement was 0.96 mm for the rest of the vertices. B₀ and gradient distortion corrections significantly reduced the fallout (to 2.1% and 8.7%) and displacement (to 0.29 mm and 0.86 mm). These corrections were effective even around regions with moderate distortions (the somatosensory and visual cortices for B₀ distortion, and the anterior frontal, inferior temporal, and posterior occipital cortices for gradient distortion).

Data Conclusion

B₀ distortion correction is crucial for surface‐based analysis of fine‐resolution fMRI at 7T. Gradient distortion correction should be considered when regions of interest include regions distant from the isocenter of scanners.

Evidence Level

1

Technical Efficacy Stage

1

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

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

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.

Recurrent Processing Drives Perceptual Plasticity

Learning and experience are critical for translating ambiguous sensory information from our environments to perceptual decisions. Yet evidence on how training molds the adult human brain remains controversial, as fMRI at standard resolution does not allow us to discern the finer scale mechanisms that underlie sensory plasticity. Here, we combine ultra-high-field (7T) functional imaging at sub-millimeter resolution with orientation discrimination training to interrogate experience-dependent plasticity across cortical depths that are known to support dissociable brain computations. We demonstrate that learning alters orientation-specific representations in superficial rather than middle or deeper V1 layers, consistent with recurrent plasticity mechanisms via horizontal connections. Further, learning increases feedforward rather than feedback layer-to-layer connectivity in occipito-parietal regions, suggesting that sensory plasticity gates perceptual decisions. Our findings reveal finer scale plasticity mechanisms that re-weight sensory signals to inform improved decisions, bridging the gap between micro- and macro-circuits of experience-dependent plasticity.

•

Discrimination training alters orientation representations in superficial V1 layers

•

Orientation-specific V1 plasticity is independent of task context

•

Discrimination training alters orientation representations in middle IPS layers

•

Learning enhances feedforward connectivity from visual to parietal cortex

Adaptation to feedback representation of illusory orientation produced from flash grab effect

Adaptation is a ubiquitous property of sensory systems. It is typically considered that neurons adapt to dominant energy in the ambient environment to function optimally. However, perceptual representation of the stimulus, often modulated by feedback signals, sometimes do not correspond to the input state of the stimulus, which tends to be more linked with feedforward signals. Here we investigated the relative contributions to cortical adaptation from feedforward and feedback signals, taking advantage of a visual illusion, the Flash-Grab Effect, to disassociate the feedforward and feedback representation of an adaptor. Results reveal that orientation adaptation is exclusively dependent on the perceived rather than the retinal orientation of the adaptor. Combined fMRI and EEG measurements demonstrate that the perceived orientation of the Flash-Grab Effect is indeed supported by feedback signals in the cortex. These findings highlight the important contribution of feedback signals for cortical neurons to recalibrate their sensitivity.

Feedforward-feedback signal interactions are common in the brain during sensory information processing. Here, the authors show that feedback-driven representation of perceived orientation dominates visual adaptation, despite the discrepant feedforward representation of input orientation.

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.