Ultra-high field MRI: Advancing systems neuroscience towards mesoscopic human brain function

Nanoparticulated Bimodal Contrast Agent for Ultra-High-Field Magnetic Resonance Imaging and Spectral X-ray Computed Tomography

Bimodal medical imaging based on magnetic resonance imaging (MRI) and computed tomography (CT) is a well-known strategy to increase the diagnostic accuracy. The most recent advances in MRI and CT instrumentation are related to the use of ultra-high magnetic fields (UHF-MRI) and different working voltages (spectral CT), respectively. Such advances require the parallel development of bimodal contrast agents (CAs) that are efficient under new instrumental conditions. In this work, we have synthesized, through a precipitation reaction from a glycerol solution of the precursors, uniform barium dysprosium fluoride nanospheres with a cubic fluorite structure, whose size was found to depend on the Ba/(Ba + Dy) ratio of the starting solution. Moreover, irrespective of the starting Ba/(Ba + Dy) ratio, the experimental Ba/(Ba + Dy) values were always lower than those used in the starting solutions. This result was assigned to lower precipitation kinetics of barium fluoride compared to dysprosium fluoride, as inferred from the detailed analysis of the effect of reaction time on the chemical composition of the precipitates. A sample composed of 34 nm nanospheres with a Ba0.51Dy0.49F2.49 stoichiometry showed a transversal relaxivity (r2) value of 147.11 mM-1·s-1 at 9.4 T and gave a high negative contrast in the phantom image. Likewise, it produced high X-ray attenuation in a large range of working voltages (from 80 to 140 kVp), which can be attributed to the presence of different K-edge values and high Z elements (Ba and Dy) in the nanospheres. Finally, these nanospheres showed negligible cytotoxicity for different biocompatibility tests. Taken together, these results show that the reported nanoparticles are excellent candidates for UHF-MRI/spectral CT bimodal imaging CAs.

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

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

RF coil design strategies for improving SNR at the ultrahigh magnetic field of 10.5 Tesla

Purpose

To develop multichannel transmit and receive arrays towards capturing the ultimate-intrinsic-SNR (uiSNR) at 10.5 Tesla (T) and to demonstrate the feasibility and potential of whole-brain, high-resolution human brain imaging at this high field strength.

Methods

A dual row 16-channel self-decoupled transmit (Tx) array was converted to a 16Tx/Rx transceiver using custom transmit/receive switches. A 64-channel receive-only (64Rx) array was built to fit into the 16Tx/Rx array. Electromagnetic modeling and experiments were employed to define safe operation limits of the resulting 16Tx/80Rx array and obtain FDA approval for human use.

Results

The 64Rx array alone captured approximately 50% of the central uiSNR at 10.5T while the identical 7T 64Rx array captured ∼76% of uiSNR at this lower field strength. The 16Tx/80Rx configuration brought the fraction of uiSNR captured at 10.5T to levels comparable to the performance of the 64Rx array at 7T. SNR data obtained at the two field strengths with these arrays displayed dependent increases over a large central region. Whole-brain high resolution T 2 * and T 1 weighted anatomical and gradient-recalled echo EPI BOLD fMRI images were obtained at 10.5T for the first time with such an advanced array, illustrating the promise of >10T fields in studying the human brain.

Conclusion

We demonstrated the ability to approach the uiSNR at 10.5T over the human brain with a novel, high channel count array, achieving large SNR gains over 7T, currently the most commonly employed ultrahigh field platform, and demonstrate high resolution and high contrast anatomical and functional imaging at 10.5T.

High Spatiotemporal Resolution Radial Encoding Single-Vessel fMRI

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Amidst an amygdala renaissance in Alzheimer's disease

The amygdala was highlighted as an early site for neurofibrillary tau tangle pathology in Alzheimer's disease in the seminal 1991 article by Braak and Braak. This knowledge has, however, only received traction recently with advances in imaging and image analysis techniques. Here, we provide a cross-disciplinary overview of pathology and neuroimaging studies on the amygdala. These studies provide strong support for an early role of the amygdala in Alzheimer's disease and the utility of imaging biomarkers of the amygdala in detecting early changes and predicting decline in cognitive functions and neuropsychiatric symptoms in early stages. We summarize the animal literature on connectivity of the amygdala, demonstrating that amygdala nuclei that show the earliest and strongest accumulation of neurofibrillary tangle pathology are those that are connected to brain regions that also show early neurofibrillary tangle accumulation. Additionally, we propose an alternative pathway of neurofibrillary tangle spreading within the medial temporal lobe between the amygdala and the anterior hippocampus. The proposed existence of this pathway is strengthened by novel experimental data on human functional connectivity. Finally, we summarize the functional roles of the amygdala, highlighting the correspondence between neurofibrillary tangle accumulation and symptomatic profiles in Alzheimer's disease. In summary, these findings provide a new impetus for studying the amygdala in Alzheimer's disease and a unique perspective to guide further study on neurofibrillary tangle spreading and the occurrence of neuropsychiatric symptoms in Alzheimer's disease.

The application of functional imaging in visual field defects: a brief review

Visual field defects (VFDs) represent a prevalent complication stemming from neurological and ophthalmic conditions. A range of factors, including tumors, brain surgery, glaucoma, and other disorders, can induce varying degrees of VFDs, significantly impacting patients' quality of life. Over recent decades, functional imaging has emerged as a pivotal field, employing imaging technology to illustrate functional changes within tissues and organs. As functional imaging continues to advance, its integration into various clinical aspects of VFDs has substantially enhanced the diagnostic, therapeutic, and management capabilities of healthcare professionals. Notably, prominent imaging techniques such as DTI, OCT, and MRI have garnered widespread adoption, yet they possess unique applications and considerations. This comprehensive review aims to meticulously examine the application and evolution of functional imaging in the context of VFDs. Our objective is to furnish neurologists and ophthalmologists with a systematic and comprehensive comprehension of this critical subject matter.

The temporal specificity of BOLD fMRI is systematically related to anatomical and vascular features of the human brain

The ability to detect fast responses with functional MRI depends on the speed of hemodynamic responses to neural activity, because hemodynamic responses act as a temporal low-pass filter smoothing out rapid changes. However, hemodynamic responses (their shape and timing) are highly variable across the brain and across stimuli. This heterogeneity of responses implies that the temporal specificity of fMRI signals, or the ability of fMRI to preserve fast information, should also vary substantially across the cortex. In this work we investigated how local differences in hemodynamic response timing impact the temporal specificity of fMRI. We conducted our research using ultra-high field (7T) fMRI at high spatiotemporal resolution, using the primary visual cortex (V1) as a model area for investigation. We used visual stimuli oscillating at slow and fast frequencies to probe the temporal specificity of individual voxels. As expected, we identified substantial variability in temporal specificity, with some voxels preserving their responses to fast neural activity more effectively than others. We investigated which voxels had the highest temporal specificity and related those to anatomical and vascular features of V1. We found that low temporal specificity is only weakly explained by the presence of large veins or cerebral cortical depth. Notably, however, temporal specificity depended strongly on a voxel's position along the anterior-posterior anatomical axis of V1, with voxels within the calcarine sulcus being capable of preserving close to 25% of their amplitude as the frequency of stimulation increased from 0.05-Hz to 0.20-Hz, and voxels nearest to the occipital pole preserving less than 18%. These results indicate that detection biases in high-resolution fMRI will depend on the anatomical and vascular features of the area being imaged, and that these biases will differ depending on the timing of the underlying neuronal activity. Importantly, this spatial heterogeneity of temporal specificity suggests that it could be exploited to achieve higher specificity in some locations, and that tailored data analysis strategies may help improve the detection and interpretation of fast fMRI responses.

Characterisation of laminar and vascular spatiotemporal dynamics of CBV and BOLD signals using VASO and ME-GRE at 7T in humans

Interpretation of cortical laminar functional magnetic resonance imaging (fMRI) activity requires detailed knowledge of the spatiotemporal haemodynamic response across vascular compartments due to the well-known vascular biases (e.g. the draining veins). Further complications arise from the spatiotemporal hemodynamic response that differs depending on the duration of stimulation. This information is crucial for future studies using depth-dependent cerebral blood volume (CBV) measurements, which promise higher specificity for the cortical microvasculature than the blood oxygenation level dependent (BOLD) contrast. To date, direct information about CBV dynamics with respect to stimulus duration, cortical depth and vasculature is missing in humans. Therefore, we characterized the cortical depth-dependent CBV-haemodynamic responses across a wide set of stimulus durations with 0.9 mm isotropic spatial and 0.785 seconds effective temporal resolution in humans using slice-selective slab-inversion vascular space occupancy (SS-SI VASO). Additionally, we investigated signal contributions from macrovascular compartments using fine-scale vascular information from multi-echo gradient-echo (ME-GRE) data at 0.35 mm isotropic resolution. In total, this resulted in >7.5h of scanning per participant (n=5). We have three major findings: (I) While we could demonstrate that 1 second stimulation is viable using VASO, more than 12 seconds stimulation provides better CBV responses in terms of specificity to microvasculature, but durations beyond 24 seconds of stimulation may be wasteful for certain applications. (II) We observe that CBV responses show dilation patterns across the cortex. (III) While we found increasingly strong BOLD signal responses in vessel-dominated voxels with longer stimulation durations, we found increasingly strong CBV signal responses in vessel-dominated voxels only until 4 second stimulation durations. After 4 seconds, only the signal from non-vessel dominated voxels kept increasing. This might explain why CBV responses are more specific to the underlying neuronal activity for long stimulus durations.

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.

High resolution and contrast 7 tesla MR brain imaging of the neonate

Introduction

Ultra-high field MR imaging offers marked gains in signal-to-noise ratio, spatial resolution, and contrast which translate to improved pathological and anatomical sensitivity. These benefits are particularly relevant for the neonatal brain which is rapidly developing and sensitive to injury. However, experience of imaging neonates at 7T has been limited due to regulatory, safety, and practical considerations. We aimed to establish a program for safely acquiring high resolution and contrast brain images from neonates on a 7T system.

Methods

Images were acquired from 35 neonates on 44 occasions (median age 39 + 6 postmenstrual weeks, range 33 + 4 to 52 + 6; median body weight 2.93 kg, range 1.57 to 5.3 kg) over a median time of 49 mins 30 s. Peripheral body temperature and physiological measures were recorded throughout scanning. Acquired sequences included T2 weighted (TSE), Actual Flip angle Imaging (AFI), functional MRI (BOLD EPI), susceptibility weighted imaging (SWI), and MR spectroscopy (STEAM).

Results

There was no significant difference between temperature before and after scanning (p = 0.76) and image quality assessment compared favorably to state-of-the-art 3T acquisitions. Anatomical imaging demonstrated excellent sensitivity to structures which are typically hard to visualize at lower field strengths including the hippocampus, cerebellum, and vasculature. Images were also acquired with contrast mechanisms which are enhanced at ultra-high field including susceptibility weighted imaging, functional MRI, and MR spectroscopy.

Discussion

We demonstrate safety and feasibility of imaging vulnerable neonates at ultra-high field and highlight the untapped potential for providing important new insights into brain development and pathological processes during this critical phase of early life.

A 128-channel receive array for cortical brain imaging at 7 T

PURPOSE

A 128-channel receive-only array for brain imaging at 7 T was simulated, designed, constructed, and tested within a high-performance head gradient designed for high-resolution functional imaging.

METHODS

The coil used a tight-fitting helmet geometry populated with 128 loop elements and preamplifiers to fit into a 39 cm diameter space inside a built-in gradient. The signal-to-noise ratio (SNR) and parallel imaging performance (1/g) were measured in vivo and simulated using electromagnetic modeling. The histogram of 1/g factors was analyzed to assess the range of performance. The array's performance was compared to the industry-standard 32-channel receive array and a 64-channel research array.

RESULTS

It was possible to construct the 128-channel array with body noise-dominated loops producing an average noise correlation of 5.4%. Measurements showed increased sensitivity compared with the 32-channel and 64-channel array through a combination of higher intrinsic SNR and g-factor improvements. For unaccelerated imaging, the 128-channel array showed SNR gains of 17.6% and 9.3% compared to the 32-channel and 64-channel array, respectively, at the center of the brain and 42% and 18% higher SNR in the peripheral brain regions including the cortex. For R = 5 accelerated imaging, these gains were 44.2% and 24.3% at the brain center and 86.7% and 48.7% in the cortex. The 1/g-factor histograms show both an improved mean and a tighter distribution by increasing the channel count, with both effects becoming more pronounced at higher accelerations.

CONCLUSION

The experimental results confirm that increasing the channel count to 128 channels is beneficial for 7T brain imaging, both for increasing SNR in peripheral brain regions and for accelerated imaging.

A selection and targeting framework of cortical locations for line-scanning fMRI

Depth-resolved functional magnetic resonance imaging (fMRI) is an emerging field growing in popularity given the potential of separating signals from different computational processes in cerebral cortex. Conventional acquisition schemes suffer from low spatial and temporal resolutions. Line-scanning methods allow depth-resolved fMRI by sacrificing spatial coverage to sample blood oxygenated level-dependent (BOLD) responses at ultra-high temporal and spatial resolution. For neuroscience applications, it is critical to be able to place the line accurately to (1) sample the right neural population and (2) target that neural population with tailored stimuli or tasks. To this end, we devised a multi-session framework where a target cortical location is selected based on anatomical and functional properties. The line is then positioned according to this information in a separate second session, and we tailor the experiment to focus on the target location. Anatomically, the precision of the line placement was confirmed by projecting a nominal representation of the acquired line back onto the surface. Functional estimates of neural selectivities in the line, as quantified by a visual population-receptive field model, resembled the target selectivities well for most subjects. This functional precision was quantified in detail by estimating the distance between the visual field location of the targeted vertex and the location in visual cortex (V1) that most closely resembled the line-scanning estimates; this distance was on average ~5.5 mm. Given the dimensions of the line, differences in acquisition, session, and stimulus design, this validates that line-scanning can be used to probe local neural sensitivities across sessions. In summary, we present an accurate framework for line-scanning MRI; we believe such a framework is required to harness the full potential of line-scanning and maximize its utility. Furthermore, this approach bridges canonical fMRI experiments with electrophysiological experiments, which in turn allows novel avenues for studying human physiology non-invasively.

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

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

Modeling the Visual Landscape: A Review on Approaches, Methods and Techniques

Modeling the perception and evaluation of landscapes from the human perspective is a desirable goal for several scientific domains and applications. Human vision is the dominant sense, and human eyes are the sensors for apperceiving the environmental stimuli of our surroundings. Therefore, exploring the experimental recording and measurement of the visual landscape can reveal crucial aspects about human visual perception responses while viewing the natural or man-made landscapes. Landscape evaluation (or assessment) is another dimension that refers mainly to preferences of the visual landscape, involving human cognition as well, in ways that are often unpredictable. Yet, landscape can be approached by both egocentric (i.e., human view) and exocentric (i.e., bird's eye view) perspectives. The overarching approach of this review article lies in systematically presenting the different ways for modeling and quantifying the two 'modalities' of human perception and evaluation, under the two geometric perspectives, suggesting integrative approaches on these two 'diverging' dualities. To this end, several pertinent traditions/approaches, sensor-based experimental methods and techniques (e.g., eye tracking, fMRI, and EEG), and metrics are adduced and described. Essentially, this review article acts as a 'guide-map' for the delineation of the different activities related to landscape experience and/or management and to the valid or potentially suitable types of stimuli, sensors techniques, and metrics for each activity. Throughout our work, two main research directions are identified: (1) one that attempts to transfer the visual landscape experience/management from the one perspective to the other (and vice versa); (2) another one that aims to anticipate the visual perception of different landscapes and establish connections between perceptual processes and landscape preferences. As it appears, the research in the field is rapidly growing. In our opinion, it can be greatly advanced and enriched using integrative, interdisciplinary approaches in order to better understand the concepts and the mechanisms by which the visual landscape, as a complex set of stimuli, influences visual perception, potentially leading to more elaborate outcomes such as the anticipation of landscape preferences. As an effect, such approaches can support a rigorous, evidence-based, and socially just framework towards landscape management, protection, and decision making, based on a wide spectrum of well-suited and advanced sensor-based technologies.

Functional MRI of murine olfactory bulbs at 15.2T reveals characteristic activation patters when stimulated by different odors

Thanks to its increased sensitivity, single-shot ultrahigh field functional MRI (UHF fMRI) could lead to valuable insight about subtle brain functions such as olfaction. However, UHF fMRI experiments targeting small organs next to air voids, such as the olfactory bulb, are severely affected by field inhomogeneity problems. Spatiotemporal Encoding (SPEN) is an emerging single-shot MRI technique that could provide a route for bypassing these complications. This is here explored with single-shot fMRI studies on the olfactory bulbs of male and female mice performed at 15.2T. SPEN images collected on these organs at a 108 µm in-plane resolution yielded remarkably large and well-defined responses to olfactory cues. Under suitable T2* weightings these activation-driven changes exceeded 5% of the overall signal intensity, becoming clearly visible in the images without statistical treatment. The nature of the SPEN signal intensity changes in such experiments was unambiguously linked to olfaction, via single-nostril experiments. These experiments highlighted specific activation regions in the external plexiform region and in glomeruli in the lateral part of the bulb, when stimulated by aversive or appetitive odors, respectively. These strong signal activations were non-linear with concentration, and shed light on how chemosensory signals reaching the olfactory epithelium react in response to different cues. Second-level analyses highlighted clear differences among the appetitive, aversive and neutral odor maps; no such differences were evident upon comparing male against female olfactory activation regions.

Local and long-distance organization of prefrontal cortex circuits in the marmoset brain

The prefrontal cortex (PFC) has dramatically expanded in primates, but its organization and interactions with other brain regions are only partially understood. We performed high-resolution connectomic mapping of the marmoset PFC and found two contrasting corticocortical and corticostriatal projection patterns: "patchy" projections that formed many columns of submillimeter scale in nearby and distant regions and "diffuse" projections that spread widely across the cortex and striatum. Parcellation-free analyses revealed representations of PFC gradients in these projections' local and global distribution patterns. We also demonstrated column-scale precision of reciprocal corticocortical connectivity, suggesting that PFC contains a mosaic of discrete columns. Diffuse projections showed considerable diversity in the laminar patterns of axonal spread. Altogether, these fine-grained analyses reveal important principles of local and long-distance PFC circuits in marmosets and provide insights into the functional organization of the primate brain.

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.

High-resolution motion- and phase-corrected functional MRI at 7 T using shuttered multishot echo-planar imaging

PURPOSE

To achieve high-resolution multishot echo-planar imaging (EPI) for functional MRI (fMRI) with reduced sensitivity to in-plane motion and between-shot phase variations.

METHODS

Two-dimensional radiofrequency pulses were incorporated in a multishot EPI sequence at 7T which selectively excited a set of in-plane bands (shutters) in the phase encoding direction, which moved between shots to cover the entire slice. A phase- and motion-corrected reconstruction was implemented for the acquisition. Brain imaging experiments were performed with instructed motion to evaluate image quality for conventional multishot and shuttered EPI. Temporal stability was assessed in three subjects by quantifying temporal SNR (tSNR) and artifact levels, and fMRI activation experiments using visual stimulation were performed to assess the strength and distribution of activation, using both conventional multishot and shuttered EPI.

RESULTS

In the instructed motion experiment, ghosting was lower in shuttered EPI images without or with corrections and image quality metrics were improved with motion correction. tSNR was improved by phase correction in both conventional multishot and shuttered EPI and the acquisitions had similar tSNR without and with phase correction. However, while phase correction was necessary to maximize tSNR in conventional multishot EPI, it also increased intermittent ghosting, but did not increase intermittent ghosting in shuttered EPI. Phase correction increased activation strength in both conventional multishot and shuttered EPI, but caused increased spurious activation outside the brain and in frontal brain regions in conventional multishot EPI.

CONCLUSION

Shuttered EPI supports multishot segmented EPI acquisitions with lower sensitivity to artifacts from motion for high-resolution fMRI.

A vision of 14 T MR for fundamental and clinical science

OBJECTIVE

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

RESEARCH FIELDS

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

CONCLUSIONS

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

Multi-feed, loop-dipole combined dielectric resonator antenna arrays for human brain MRI at 7 T

OBJECTIVE

To determine whether a multi-feed, loop-dipole combined approach can be used to improve performance of rectangular dielectric resonator antenna (DRA) arrays human brain for MRI at 7 T.

MATERIALS AND METHODS

Electromagnetic field simulations in a spherical phantom and human voxel model "Duke" were conducted for different rectangular DRA geometries and dielectric constants ε_r. Three types of RF feed were investigated: loop-only, dipole-only and loop-dipole. Additionally, multi-channel array configurations up to 24-channels were simulated.

RESULTS

The loop-only coupling scheme provided the highest B₁⁺ and SAR efficiency, while the loop-dipole showed the highest SNR in the center of a spherical phantom for both single- and multi-channel configurations. For Duke, 16-channel arrays outperformed an 8-channel bow-tie array with greater B₁⁺ efficiency (1.48- to 1.54-fold), SAR efficiency (1.03- to 1.23-fold) and SNR (1.63- to 1.78). The multi-feed, loop-dipole combined approach enabled the number of channels increase to 24 with 3 channels per block.

DISCUSSION

This work provides novel insights into the rectangular DRA design for high field MRI and shows that the loop-only feed should be used instead of the dipole-only in transmit mode to achieve the highest B₁⁺ and SAR efficiency, while the loop-dipole should be the best suited in receive mode to obtain the highest SNR in spherical samples of similar size and electrical properties as the human head.

Commissioning of the Iseult CEA 11.7 T whole-body MRI: current status, gradient-magnet interaction tests and first imaging experience

OBJECTIVES

The Iseult MRI is an actively shielded whole-body magnet providing a homogeneous and stable magnetic field of 11.7 T. After nearly 20 years of research and development, the magnet successfully reached its target field strength for the first time in 2019. This article reviews its commissioning status, the gradient-magnet interaction test results and first imaging experience.

MATERIALS AND METHODS

Vibration, acoustics, power deposition in the He bath, and field monitoring measurements were carried out. Magnet safety system was tested against outer magnetic perturbations, and calibrated to define a safe operation of the gradient coil. First measurements using parallel transmission were also performed on an ex-vivo brain to mitigate the RF field inhomogeneity effect.

RESULTS

Acoustics measurements show promising results with sound pressure levels slightly above the enforced limits only at certain frequency intervals. Vibrations of the gradient coil revealed a linear trend with the B₀ field only in the worst case. Field monitoring revealed some resonances at some frequencies that are still under investigation.

DISCUSSION

Gradient-magnet interaction tests at up to 11.7 T are concluded. The scanner is now kept permanently at field and the final calibrations are on-going to pave the road towards the first acquisitions on volunteers.

Mechanisms of speed encoding in the human middle temporal cortex measured by 7T fMRI

Perception of dynamic scenes in our environment results from the evaluation of visual features such as the fundamental spatial and temporal frequency components of a moving object. The ratio between these two components represents the object's speed of motion. The human middle temporal cortex hMT+ has a crucial biological role in the direct encoding of object speed. However, the link between hMT+ speed encoding and the spatiotemporal frequency components of a moving object is still under explored. Here, we recorded high resolution 7T blood oxygen level-dependent BOLD responses to different visual motion stimuli as a function of their fundamental spatial and temporal frequency components. We fitted each hMT+ BOLD response with a 2D Gaussian model allowing for two different speed encoding mechanisms: (1) distinct and independent selectivity for the spatial and temporal frequencies of the visual motion stimuli; (2) pure tuning for the speed of motion. We show that both mechanisms occur but in different neuronal groups within hMT+, with the largest subregion of the complex showing separable tuning for the spatial and temporal frequency of the visual stimuli. Both mechanisms were highly reproducible within participants, reconciling single cell recordings from MT in animals that have showed both encoding mechanisms. Our findings confirm that a more complex process is involved in the perception of speed than initially thought and suggest that hMT+ plays a primary role in the evaluation of the spatial features of the moving visual input.

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.

BigBrain-MR: a new digital phantom with anatomically-realistic magnetic resonance properties at 100-µm resolution for magnetic resonance methods development

The benefits, opportunities and growing availability of ultra-high field magnetic resonance imaging (MRI) for humans have prompted an expansion in research and development efforts towards increasingly more advanced high-resolution imaging techniques. To maximize their effectiveness, these efforts need to be supported by powerful computational simulation platforms that can adequately reproduce the biophysical characteristics of MRI, with high spatial resolution. In this work, we have sought to address this need by developing a novel digital phantom with realistic anatomical detail up to 100-µm resolution, including multiple MRI properties that affect image generation. This phantom, termed BigBrain-MR, was generated from the publicly available BigBrain histological dataset and lower-resolution in-vivo 7T-MRI data, using a newly-developed image processing framework that allows mapping the general properties of the latter into the fine anatomical scale of the former. Overall, the mapping framework was found to be effective and robust, yielding a diverse range of realistic "in-vivo-like" MRI contrasts and maps at 100-µm resolution. BigBrain-MR was then tested in three imaging applications (motion effects and interpolation, super-resolution imaging, and parallel imaging reconstruction) to investigate its properties, value and validity as a simulation platform. The results consistently showed that BigBrain-MR can closely approximate the behavior of real in-vivo data, more realistically and with more extensive features than a more classic option such as the Shepp-Logan phantom. Its flexibility in simulating different contrast mechanisms and artifacts may also prove valuable for educational applications. BigBrain-MR is therefore deemed a favorable choice to support methodological development and demonstration in brain MRI, and has been made freely available to the community.

Multimodal 3D Mouse Brain Atlas Framework with the Skull-Derived Coordinate System

Magnetic resonance imaging (MRI) and light-sheet fluorescence microscopy (LSFM) are technologies that enable non-disruptive 3-dimensional imaging of whole mouse brains. A combination of complementary information from both modalities is desirable for studying neuroscience in general, disease progression and drug efficacy. Although both technologies rely on atlas mapping for quantitative analyses, the translation of LSFM recorded data to MRI templates has been complicated by the morphological changes inflicted by tissue clearing and the enormous size of the raw data sets. Consequently, there is an unmet need for tools that will facilitate fast and accurate translation of LSFM recorded brains to in vivo, non-distorted templates. In this study, we have developed a bidirectional multimodal atlas framework that includes brain templates based on both imaging modalities, region delineations from the Allen's Common Coordinate Framework, and a skull-derived stereotaxic coordinate system. The framework also provides algorithms for bidirectional transformation of results obtained using either MR or LSFM (iDISCO cleared) mouse brain imaging while the coordinate system enables users to easily assign in vivo coordinates across the different brain templates.

Diagnostic Performance of 0.55 T MRI for Intracranial Aneurysm Detection

OBJECTIVES

Intracranial aneurysm (IA) is the main cause of subarachnoid hemorrhages. Time-of-flight (TOF) magnetic resonance angiography (MRA) at 1.5 T or 3 T magnetic resonance imaging (MRI) is a well-established method for the diagnosis of IA. The aim of this prospective study was to evaluate the performance of a modern 0.55 T MRI in the diagnosis of IAs in comparison to digital subtraction angiography (DSA) as a standard of reference.

MATERIALS AND METHODS

Seventeen patients with suspicion of single or multiple IAs underwent TOF MRA at 0.55 T MRI 1 day before DSA. Two neuroradiologists independently measured the aneurysm neck, width, and height on 0.55 T, 1.5 T, and 3 T 3D-TOF MRA source images and 2D/3D rotational angiography. The main analysis assessed the intermodality agreement between 0.55 T TOF MRA and DSA using Bland-Altman plots, a Wilcoxon test, and the intraclass correlation coefficient (ICC). In a secondary analysis, aneurysm dimensions were compared between 0.55 T TOF MRA and 1.5/3 T TOF MRA. Interreader agreement was evaluated by ICC. A third neuroradiologist blinded to patient history screened 0.55 T TOF MRA data sets of the aforementioned 17 patients and 15 additional healthy patients for the presence and location of aneurysms.

RESULTS

A total of 19 aneurysms in 16 patients were identified in both 0.55 T MRA and DSA. Measurements of the 2 nonblinded readers showed no significant differences between 0.55 T TOF MRA and DSA in the overall aneurysm size (calculated as the mean from height/width/neck) ( P = 0.178), as well as in the mean width ( P = 0.778) and neck values ( P = 0.190). The mean height was significantly larger in 0.55 T TOF MRA in comparison to DSA ( P = 0.020). Intermodality (1.5/3 T TOF MRA) and interrater agreement were excellent (ICC > 0.94). Of the 32 data sets of patients with and without IA, the blinded reader detected all aneurysms correctly by using 0.55 T images.

CONCLUSIONS

TOF-MRA acquired with a modern 0.55 T MRI is a reliable tool for the detection and initial assessment of IAs.

Robust high spatio-temporal line-scanning fMRI in humans at 7T using multi-echo readouts, denoising and prospective motion correction

BACKGROUND

Functional magnetic resonance imaging (fMRI), typically using blood oxygenation level-dependent (BOLD) contrast weighted imaging, allows the study of brain function with millimeter spatial resolution and temporal resolution of one to a few seconds. At a mesoscopic scale, neurons in the human brain are spatially organized in structures with dimensions of hundreds of micrometers, while they communicate at the millisecond timescale. For this reason, it is important to develop an fMRI method with simultaneous high spatial and temporal resolution. Line-scanning promises to reach this goal at the cost of volume coverage.

NEW METHOD

Here, we release a comprehensive update to human line-scanning fMRI. First, we investigated multi-echo line-scanning with five different protocols varying the number of echoes and readout bandwidth while keeping the TR constant. In these, we compared different echo combination approaches in terms of BOLD activation (sensitivity) and temporal signal-to-noise ratio. Second, we implemented an adaptation of NOise reduction with DIstribution Corrected principal component analysis (NORDIC) thermal noise removal for line-scanning fMRI data. Finally, we tested three image-based navigators for motion correction and investigated different ways of performing fMRI analysis on the timecourses which were influenced by the insertion of the navigators themselves.

RESULTS

The presented improvements are relatively straightforward to implement; multi-echo readout and NORDIC denoising together, significantly improve data quality in terms of tSNR and t-statistical values, while motion correction makes line-scanning fMRI more robust.

COMPARISON WITH EXISTING METHODS

Multi-echo acquisitions and denoising have previously been applied in 3D magnetic resonance imaging. Their combination and application to 1D line-scanning is novel. The current proposed method greatly outperforms the previous line-scanning acquisitions with single-echo acquisition, in terms of tSNR (4.0 for single-echo line-scanning and 36.2 for NORDIC-denoised multi-echo) and t-statistical values (3.8 for single-echo line-scanning and 25.1 for NORDIC-denoised multi-echo line-scanning).

CONCLUSIONS

Line-scanning fMRI was advanced compared to its previous implementation in order to improve sensitivity and reliability. The improved line-scanning acquisition could be used, in the future, for neuroscientific and clinical applications.

Improved Selectivity in 7 T Digit Mapping Using VASO-CBV

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

Fluorescently-tagged magnetic protein nanoparticles for high-resolution optical and ultra-high field magnetic resonance dual-modal cerebral angiography

Extremely small paramagnetic iron oxide nanoparticles (FeMNPs) (<5 nm) can enhance positive magnetic resonance imaging (MRI) contrast by shortening the longitudinal relaxation time of water (T₁), but these nanoparticles experience rapid renal clearance. Here, magnetic protein nanoparticles (MPNPs) are synthesized from protein-conjugated citric acid coated FeMNPs (c-FeMNPs) without loss of the T₁ MRI properties and tagged with fluorescent dye (f-MPNPs) for optical cerebrovascular imaging. The c-FeMNPs shows average size 3.8 ± 0.7 nm with T₁ relaxivity (r₁) of 1.86 mM^-1 s^-1 and transverse/longitudinal relaxivity ratio (r₂/r₁) of 2.53 at 11.7 T. The f-MPNPs show a higher r₁ value of 2.18 mM^-1 s^-1 and r₂/r₁ ratio of 2.88 at 11.7 T, which generates excellent positive MRI contrast. In vivo cerebral angiography with f-MPNPs enables detailed microvascular contrast enhancement for differentiation of major blood vessels of murine brain, which corresponds well with whole brain three-dimensional time-of-flight MRI angiograms (17 min imaging time with 60 ms repetition time and 40 μm isotropic voxels). The real-time fluorescence angiography enables unambiguous detection of brain capillaries with diameter < 40 μm. Biodistribution examination revealed that f-MPNPs were safely cleared by the organs like the liver, spleen, and kidneys within a day after injection. Blood biochemical assays demonstrated no risk of iron overload in both rats and mice. With hybrid neuroimaging technologies (e.g., MRI-optical) on the rise, f-MPNPs built on this platform can generate exciting neuroscience applications.

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.

Baseline oxygen consumption decreases with cortical depth

The cerebral cortex is organized in cortical layers that differ in their cellular density, composition, and wiring. Cortical laminar architecture is also readily revealed by staining for cytochrome oxidase-the last enzyme in the respiratory electron transport chain located in the inner mitochondrial membrane. It has been hypothesized that a high-density band of cytochrome oxidase in cortical layer IV reflects higher oxygen consumption under baseline (unstimulated) conditions. Here, we tested the above hypothesis using direct measurements of the partial pressure of O2 (pO2) in cortical tissue by means of 2-photon phosphorescence lifetime microscopy (2PLM). We revisited our previously developed method for extraction of the cerebral metabolic rate of O2 (CMRO2) based on 2-photon pO2 measurements around diving arterioles and applied this method to estimate baseline CMRO2 in awake mice across cortical layers. To our surprise, our results revealed a decrease in baseline CMRO2 from layer I to layer IV. This decrease of CMRO2 with cortical depth was paralleled by an increase in tissue oxygenation. Higher baseline oxygenation and cytochrome density in layer IV may serve as an O2 reserve during surges of neuronal activity or certain metabolically active brain states rather than reflecting baseline energy needs. Our study provides to our knowledge the first quantification of microscopically resolved CMRO2 across cortical layers as a step towards better understanding of brain energy metabolism.

Reinforcement learning with associative or discriminative generalization across states and actions: fMRI at 3 T and 7 T

The model-free algorithms of "reinforcement learning" (RL) have gained clout across disciplines, but so too have model-based alternatives. The present study emphasizes other dimensions of this model space in consideration of associative or discriminative generalization across states and actions. This "generalized reinforcement learning" (GRL) model, a frugal extension of RL, parsimoniously retains the single reward-prediction error (RPE), but the scope of learning goes beyond the experienced state and action. Instead, the generalized RPE is efficiently relayed for bidirectional counterfactual updating of value estimates for other representations. Aided by structural information but as an implicit rather than explicit cognitive map, GRL provided the most precise account of human behavior and individual differences in a reversal-learning task with hierarchical structure that encouraged inverse generalization across both states and actions. Reflecting inference that could be true, false (i.e., overgeneralization), or absent (i.e., undergeneralization), state generalization distinguished those who learned well more so than action generalization. With high-resolution high-field fMRI targeting the dopaminergic midbrain, the GRL model's RPE signals (alongside value and decision signals) were localized within not only the striatum but also the substantia nigra and the ventral tegmental area, including specific effects of generalization that also extend to the hippocampus. Factoring in generalization as a multidimensional process in value-based learning, these findings shed light on complexities that, while challenging classic RL, can still be resolved within the bounds of its core computations.

RF shielding designs for birdcage coils for preclinical MRI at 9.4 T

Birdcage coils are widely used in preclinical MRI as they perform well, allow for quadrature drive, and can provide a homogeneous transmit field. Unlike in larger bore scanners, an RF shield is essential to avoid strong cross-talk with gradient coils that are in close proximity. However, gradient switching induces eddy currents that heat the shield and coil and impair the temporal signal-to-noise ratio (tSNR). The motivation of this study is to investigate the performance of different designs of RF shields on a birdcage coil used for high resolution functional MRI of small primates at 9.4 T. We found the choice of materials for RF shields significantly affected ghosting and tSNR in fMRI scans. Both ultrathin foils and a slotted pattern reduce eddy currents and improve imaging quality. Our results also demonstrate that a 9-um-thick copper foil is sufficiently thin to reduce the eddy current effects for high-resolution fMRI scans and there is no need for high-cost 4-um-thick foil. For slotted shields, our results demonstrate that the number of slots should be carefully considered, and an excessive number of slots can lead to a lower SNR and tSNR. We believe the results from this study can be used as a reference to design future RF coil shields selection for preclinical scanners.

The influence of the subcortex and brain stem on overeating: How advances in functional neuroimaging can be applied to expand neurobiological models to beyond the cortex

Functional neuroimaging has become a widely used tool in obesity and eating disorder research to explore the alterations in neurobiology that underlie overeating and binge eating behaviors. Current and traditional neurobiological models underscore the importance of impairments in brain systems supporting reward, cognitive control, attention, and emotion regulation as primary drivers for overeating. Due to the technical limitations of standard field strength functional magnetic resonance imaging (fMRI) scanners, human neuroimaging research to date has focused largely on cortical and basal ganglia effects on appetitive behaviors. The present review draws on animal and human research to highlight how neural signaling encoding energy regulation, reward-learning, and habit formation converge on hypothalamic, brainstem, thalamic, and striatal regions to contribute to overeating in humans. We also consider the role of regions such as the mediodorsal thalamus, ventral striatum, lateral hypothalamus and locus coeruleus in supporting habit formation, inhibitory control of food craving, and attentional biases. Through these discussions, we present proposals on how the neurobiology underlying these processes could be examined using functional neuroimaging and highlight how ultra-high field 7-Tesla (7 T) fMRI may be leveraged to elucidate the potential functional alterations in subcortical networks. Focus is given to how interactions of these regions with peripheral endocannabinoids and neuropeptides, such as orexin, could be explored. Technical and methodological aspects regarding the use of ultra-high field 7 T fMRI to study eating behaviors are also reviewed.

NMR and Metabolomics-A Roadmap for the Future

Metabolomics investigates global metabolic alterations associated with chemical, biological, physiological, or pathological processes. These metabolic changes are measured with various analytical platforms including liquid chromatography-mass spectrometry (LC-MS), gas chromatography-mass spectrometry (GC-MS) and nuclear magnetic resonance spectroscopy (NMR). While LC-MS methods are becoming increasingly popular in the field of metabolomics (accounting for more than 70% of published metabolomics studies to date), there are considerable benefits and advantages to NMR-based methods for metabolomic studies. In fact, according to PubMed, more than 926 papers on NMR-based metabolomics were published in 2021-the most ever published in a given year. This suggests that NMR-based metabolomics continues to grow and has plenty to offer to the scientific community. This perspective outlines the growing applications of NMR in metabolomics, highlights several recent advances in NMR technologies for metabolomics, and provides a roadmap for future advancements.

Magnetic field strength dependent SNR gain at the center of a spherical phantom and up to 11.7T

PURPOSE

The SNR at the center of a spherical phantom of known electrical properties was measured in quasi-identical experimental conditions as a function of magnetic field strength between 3 T and 11.7 T.

METHODS

The SNR was measured at the center of a spherical water saline phantom with a gradient recalled echo sequence. Measurements were performed at NeuroSpin at 3, 7, and 11.7 T. The phantom was then shipped to Maastricht University and then to the University of Minnesota for additional data points at 7, 9.4, and 10.5 T. Experiments were carried out with the exact same type of birdcage volume coil (except at 3 T, where a similar coil was used) to attempt at isolating the evolution of SNR with field strength alone. Phantom electrical properties were characterized over the corresponding frequency range.

RESULTS

Electrical properties were found to barely vary over the frequency range. Removing the influence of the flip-angle excitation inhomogeneity was crucial, as expected. After such correction, measurements revealed a gain of SNR growing as B₀ ^{1.94 ± 0.16} compared with B₀ ^2.13 according to ultimate intrinsic SNR theory.

CONCLUSIONS

By using quasi-identical experimental setups (RF volume coil, phantom, electrical properties, and protocol), this work reports experimental data between 3 T and 11.7 T, enabling the comparison with SNR theories in which conductivity and permittivity can be assumed to be constant with respect to field strength. According to ultimate SNR theory, these results can be reasonably extrapolated to the performance of receive arrays with greater than about 32 elements for central SNR in the same spherical phantom.

A novel technology for in vivo detection of cell type-specific neural connection with AQP1-encoding rAAV2-retro vector and metal-free MRI

A mammalian brain contains numerous neurons with distinct cell types for complex neural circuits. Virus-based circuit tracing tools are powerful in tracking the interaction among the different brain regions. However, detecting brain-wide neural networks in vivo remains challenging since most viral tracing systems rely on postmortem optical imaging. We developed a novel approach that enables in vivo detection of brain-wide neural connections based on metal-free magnetic resonance imaging (MRI). The recombinant adeno-associated virus (rAAV) with retrograde ability, the rAAV2-retro, encoding the human water channel aquaporin 1 (AQP1) MRI reporter gene was generated to label neural connections. The mouse was micro-injected with the virus at the Caudate Putamen (CPU) region and subjected to detection with Diffusion-weighted MRI (DWI). The prominent structure of the CPU-connected network was clearly defined. In combination with a Cre-loxP system, rAAV2-retro expressing Cre-dependent AQP1 provides a CPU-connected network of specific type neurons. Here, we established a sensitive, metal-free MRI-based strategy for in vivo detection of cell type-specific neural connections in the whole brain, which could visualize the dynamic changes of neural networks in rodents and potentially in non-human primates.

Recent Technical Advances in Accelerating the Clinical Translation of Small Animal Brain Imaging: Hybrid Imaging, Deep Learning, and Transcriptomics

Small animal models play a fundamental role in brain research by deepening the understanding of the physiological functions and mechanisms underlying brain disorders and are thus essential in the development of therapeutic and diagnostic imaging tracers targeting the central nervous system. Advances in structural, functional, and molecular imaging using MRI, PET, fluorescence imaging, and optoacoustic imaging have enabled the interrogation of the rodent brain across a large temporal and spatial resolution scale in a non-invasively manner. However, there are still several major gaps in translating from preclinical brain imaging to the clinical setting. The hindering factors include the following: (1) intrinsic differences between biological species regarding brain size, cell type, protein expression level, and metabolism level and (2) imaging technical barriers regarding the interpretation of image contrast and limited spatiotemporal resolution. To mitigate these factors, single-cell transcriptomics and measures to identify the cellular source of PET tracers have been developed. Meanwhile, hybrid imaging techniques that provide highly complementary anatomical and molecular information are emerging. Furthermore, deep learning-based image analysis has been developed to enhance the quantification and optimization of the imaging protocol. In this mini-review, we summarize the recent developments in small animal neuroimaging toward improved translational power, with a focus on technical improvement including hybrid imaging, data processing, transcriptomics, awake animal imaging, and on-chip pharmacokinetics. We also discuss outstanding challenges in standardization and considerations toward increasing translational power and propose future outlooks.

The guide to dendritic spikes of the mammalian cortex in vitro and in vivo

Half a century since their discovery by Llinás and colleagues, dendritic spikes have been observed in various neurons in different brain regions, from the neocortex and cerebellum to the basal ganglia. Dendrites exhibit a terrifically diverse but stereotypical repertoire of spikes, sometimes specific to subregions of the dendrite. Despite their prevalence, we only have a glimpse into their role in the behaving animal. This article aims to survey the full range of dendritic spikes found in excitatory and inhibitory neurons, compare them in vivo versus in vitro, and discuss new studies describing dendritic spikes in the human cortex. We focus on dendritic spikes in neocortical and hippocampal neurons and present a roadmap to identify and understand the broader role of dendritic spikes in single-cell computation.

Layer-specific activation in human primary somatosensory cortex during tactile temporal prediction error processing

The human brain continuously generates predictions of incoming sensory input and calculates corresponding prediction errors from the perceived inputs to update internal predictions. In human primary somatosensory cortex (area 3b), different cortical layers are involved in receiving the sensory input and generation of error signals. It remains unknown, however, how the layers in the human area 3b contribute to the temporal prediction error processing. To investigate prediction error representation in the area 3b across layers, we acquired layer-specific functional magnetic resonance imaging (fMRI) data at 7T from human area 3b during a task of index finger poking with no-delay, short-delay and long-delay touching sequences. We demonstrate that all three tasks increased activity in both superficial and deep layers of area 3b compared to the random sensory input. The fMRI signal was differentially modulated solely in the deep layers rather than the superficial layers of area 3b by the delay time. Compared with the no-delay stimuli, activity was greater in the deep layers of area 3b during the short-delay stimuli but lower during the long-delay stimuli. This difference activity features in the superficial and deep layers suggest distinct functional contributions of area 3b layers to tactile temporal prediction error processing. The functional segregation in area 3b across layers may reflect that the excitatory and inhibitory interplay in the sensory cortex contributions to flexible communication between cortical layers or between cortical areas.

New avenues for functional neuroimaging: ultra-high field MRI and OPM-MEG

Functional brain imaging technology has developed rapidly in recent years. On the one hand, high-field 7-Tesla magnetic resonance imaging (MRI) has excelled the limited spatial resolution of 3-Tesla MRI, allowing us to enter a new world of mesoscopic imaging from the macroscopic imaging of human brain functions. On the other hand, novel optical pumping magnetometer-magnetoencephalography (OPM-MEG) has broken down the technical barriers of traditional superconducting MEG, which brings imaging of neuronal electromagnetic signals from cortical imaging to whole-brain imaging. This article aims to present a brief introduction regarding the development of conventional MRI and MEG technology, and, more importantly, to delineate that high-field MRI and OPM-MEG complement each other and together will lead us into a new era of functional brain imaging.

A nine-channel transmit/receive array for spine imaging at 10.5 T: Introduction to a nonuniform dielectric substrate antenna

PURPOSE

The purpose of this study is to introduce a new antenna element with improved transmit performance, named the nonuniform dielectric substrate (NODES) antenna, for building transmit arrays at ultrahigh-field.

METHODS

We optimized a dipole antenna at 10.5 Tesla by maximizing the B 1 + -SAR efficiency in a phantom for a human spine target. The optimization parameters included permittivity variation in the substrate, substrate thickness, antenna length, and conductor geometry. We conducted electromagnetic simulations as well as phantom experiments to compare the transmit/receive performance of the proposed NODES antenna design with existing coil elements from the literature.

RESULTS

Single NODES element showed up to 18% and 30% higher B 1 + -SAR efficiency than the fractionated dipole and loop elements, respectively. The new element is substantially shorter than a commonly used dipole, which enables z-stacked array formation; it is additionally capable of providing a relatively uniform current distribution along its conductors. The nine-channel transmit/receive NODES array achieved 7.5% higher B 1 + homogeneity than a loop array with the same number of elements. Excitation with the NODES array resulted in 33% lower peak 10g-averaged SAR and required 34% lower input power than the loop array for the target anatomy of the spine.

CONCLUSION

In this study, we introduced a new RF coil element: the NODES antenna. NODES antenna outperformed the widely used loop and dipole elements and may provide improved transmit/receive performance for future ultrahigh field MRI applications.

The Global Configuration of Visual Stimuli Alters Co-Fluctuations of Cross-Hemispheric Human Brain Activity

We tested how a stimulus gestalt, defined by the neuronal interaction between local and global features of a stimulus, is represented within human primary visual cortex (V1). We used high-resolution fMRI, which serves as a surrogate of neuronal activation, to measure co-fluctuations within subregions of V1 as (male and female) subjects were presented with peripheral stimuli, each with different global configurations. We found stronger cross-hemisphere correlations when fine-scale V1 cortical subregions represented parts of the same object compared with different objects. This result was consistent with the vertical bias in global processing and, critically, was independent of the task and local discontinuities within objects. Thus, despite the relatively small receptive fields of neurons within V1, global stimulus configuration affects neuronal processing via correlated fluctuations between regions that represent different sectors of the visual field.SIGNIFICANCE STATEMENT We provide the first evidence for the impact of global stimulus configuration on cross-hemispheric fMRI fluctuations, measured in human primary visual cortex. Our results are consistent with changes in the level of γ-band synchrony, which has been shown to be affected by global stimulus configuration, being reflected in the level fMRI co-fluctuations. These data help narrow the gap between knowledge of global stimulus configuration encoding at the single-neuron level versus at the behavioral level.

Point-spread function of the BOLD response across columns and cortical depth in human extra-striate cortex

Columns and layers are fundamental organizational units of the brain. Well known examples of cortical columns are the ocular dominance columns (ODCs) in primary visual cortex and the column-like stripe-based arrangement in the second visual area V2. The spatial scale of columns and layers is beyond the reach of conventional neuroimaging, but the advent of high field magnetic resonance imaging (MRI) scanners (UHF, 7 Tesla and above) has opened the possibility to acquire data at this spatial scale, in-vivo and non-invasively in humans. The most prominent non-invasive technique to measure brain function is blood oxygen level dependent (BOLD) fMRI, measuring brain activity indirectly, via changes in hemodynamics. A key determinant of the ability of high-resolution BOLD fMRI to accurately resolve columns and layers is the point-spread function (PSF) of the BOLD response in relation to the spatial extent of neuronal activity. In this study we take advantage of the stripe-based arrangement present in visual area V2, coupled with sub-millimetre anatomical and gradient-echo BOLD (GE BOLD) acquisition at 7 T to obtain PSF estimates and along cortical depth in human participants. Results show that the BOLD PSF is maximal in the superficial part of the cortex (1.78 mm), and it decreases with increasing cortical depth (0.83 mm close to white matter).

Divisive normalization unifies disparate response signatures throughout the human visual hierarchy

A canonical neural computation is a mathematical operation applied by the brain in a wide variety of contexts and capable of explaining and unifying seemingly unrelated neural and perceptual phenomena. Here, we use a combination of state-of-the-art experiments (ultra-high-field functional MRI) and mathematical methods (population receptive field [pRF] modeling) to uniquely demonstrate the role of divisive normalization (DN) as the canonical neural computation underlying visuospatial responses throughout the human visual hierarchy. The DN pRF model provides a tool to investigate and interpret the computational processes underlying neural responses in human and animal recordings, but also in clinical and cognitive dimensions.

Neural processing is hypothesized to apply the same mathematical operations in a variety of contexts, implementing so-called canonical neural computations. Divisive normalization (DN) is considered a prime candidate for a canonical computation. Here, we propose a population receptive field (pRF) model based on DN and evaluate it using ultra-high-field functional MRI (fMRI). The DN model parsimoniously captures seemingly disparate response signatures with a single computation, superseding existing pRF models in both performance and biological plausibility. We observe systematic variations in specific DN model parameters across the visual hierarchy and show how they relate to differences in response modulation and visuospatial information integration. The DN model delivers a unifying framework for visuospatial responses throughout the human visual hierarchy and provides insights into its underlying information-encoding computations. These findings extend the role of DN as a canonical computation to neuronal populations throughout the human visual hierarchy.

A line through the brain: implementation of human line-scanning at 7T for ultra-high spatiotemporal resolution fMRI

Functional magnetic resonance imaging (fMRI) is a widely used tool in neuroscience to detect neurally evoked responses, e.g. the blood oxygenation level-dependent (BOLD) signal. Typically, BOLD fMRI has millimeter spatial resolution and temporal resolution of one to few seconds. To study the sub-millimeter structures and activity of the cortical gray matter, the field needs an fMRI method with high spatial and temporal resolution. Line-scanning fMRI achieves very high spatial resolution and high sampling rate, at the cost of a sacrifice in volume coverage. Here, we present a human line-scanning implementation on a 7T MRI system. First, we investigate the quality of the saturation pulses that suppress MR signal outside the line. Second, we established the best coil combination for reconstruction. Finally, we applied the line-scanning method in the occipital lobe during a visual stimulation task, showing BOLD responses along cortical depth, every 250 µm with a 200 ms repetition time (TR). We found a good correspondence of t-statistics values with 2D gradient-echo echo planar imaging (GE-EPI) BOLD fMRI data with the same temporal resolution and voxel volume (R = 0.6 ± 0.2). In summary, we demonstrate the feasibility of line-scanning in humans and this opens line-scanning fMRI for applications in cognitive and clinical neuroscience.

Laminar processing of numerosity supports a canonical cortical microcircuit in human parietal cortex

As neural signals travel through the visual hierarchy, spatial precision decreases and specificity for stimulus features increases.^1-4 A similar hierarchy has been found for laminar processing in V1, where information from the thalamus predominantly targets the central layers, while spatial precision decreases and feature specificity increases toward superficial and deeper layers.^5-17 This laminar processing scheme is proposed to represent a canonical cortical microcircuit that is similar across the cortex.^11,18-21 Here, we go beyond early visual cortex and investigate whether processing of numerosity (the set size of a group of items) across cortical depth in the parietal association cortex follows this hypothesis. Numerosity processing is implicated in many tasks such as multiple object tracking,²² mathematics,^23-25 decision making,²⁶ and dividing attention.²⁷ Neurons in the parietal association cortex are tuned to numerosity, with both a preferred numerosity tuning and tuning width (i.e., specificity).^28-30 We quantified preferred numerosity responses across cortical depth in the parietal association cortex with ultra-high field fMRI and population receptive field-based numerosity modeling.^1,28,31 We find that numerosity responses sharpen, i.e., become increasingly specific, moving away from the central layers. This suggests that the laminar processing scheme for numerosity processing in the parietal cortex is similar to primary visual cortex, providing support for the canonical cortical microcircuit hypothesis beyond primary visual cortex.