Automatic cortical surface reconstruction of high-resolution T1 echo planar imaging data

Feasibility of interleaved multislice averaged magnetization inversion-recovery acquisitions of the spinal cord

PURPOSE

To establish an interleaved multislice variant of the averaged magnetization inversion-recovery acquisitions (AMIRA) approach for 2D spinal cord imaging with increased acquisition efficiency compared with the conventional 2D single-slice approach(es), and to determine essential prerequisites for a working interleaved multislice AMIRA approach in practice.

METHODS

The general AMIRA concept is based on an inversion recovery-prepared, segmented, and time-limited cine balanced SSFP sequence, generating images of different contrast. For AMIRA imaging of multiple, independent slices in a 2D interleaved fashion, a slice loop within the acquisition loops was programmed. The former non-selective inversions were replaced with slice-selective inversions with user-definable slice thickness.

RESULTS

The thickness of the slice-selective inversion in 2D interleaved multislice AMIRA should be doubled compared with the manufacturer's standard setting to avoid an increased sensitivity to flow and pulsation effects particularly in the CSF. However, this solution also limits its practical applicability, as slices located at directly adjacent vertebrae cannot be imaged together. Successful interleaved two-slice AMIRA imaging for a "reference" in vivo protocol with 0.50 × 0.50 mm2 in-plane resolution and 8-mm slice thickness is demonstrated, therefore halving its acquisition time per slice from 3 min down to 1.5 min.

CONCLUSION

The investigated 2D interleaved two-slice AMIRA variant facilitates spinal cord imaging that maintains similar contrast and the same resolution as the conventional 2D single-slice AMIRA approach, but does so with a halved acquisition time.

T1234: A distortion-matched structural scan solution to misregistration of high resolution fMRI data

Purpose

High-resolution fMRI at 7T is challenged by suboptimal alignment quality between functional data and structural scans. This study aims to develop a rapid acquisition method that provides distortion-matched, artifact-mitigated structural reference data.

Methods

We introduce an efficient sequence protocol termed T1234, which offers adjustable distortions. This approach involves a T1-weighted 2-inversion 3D-EPI sequence with four spatial encoding directions optimized for high-resolution fMRI. A forward Bloch model was used for T1 quantification and protocol optimization. Twenty participants were scanned at 7T using both structural and functional protocols to evaluate the utility of T1234.

Results

Results from two protocols are presented. A fast distortion-free protocol reliably produced whole-brain segmentations at 0.8mm isotropic resolution within 3:00-3:40 minutes. It demonstrates robustness across sessions, participants, and three different 7T SIEMENS scanners. For a protocol with geometric distortions that matched functional data, T1234 facilitates layer-specific fMRI signal analysis with enhanced laminar precision.

Conclusion

This structural mapping approach enables precise registration with fMRI data. T1234 has been successfully implemented, validated, and tested, and is now available to users at our center and at over 50 centers worldwide.

Simultaneous Multiparameter Mapping of the Liver in a Single Breath-Hold or Respiratory-Triggered Acquisition Using Multi-Inversion Spin and Gradient Echo MRI

BACKGROUND

Quantitative parametric mapping is an increasingly important tool for noninvasive assessment of chronic liver disease. Conventional parametric mapping techniques require multiple breath-held acquisitions and provide limited anatomic coverage.

PURPOSE

To investigate a multi-inversion spin and gradient echo (MI-SAGE) technique for simultaneous estimation of T1, T2, and T2* of the liver.

STUDY TYPE

Prospective.

SUBJECTS

Sixteen research participants, both adult and pediatric (age 17.5 ± 4.6 years, eight male), with and without known liver disease (seven asymptomatic healthy controls, two fibrotic liver disease, five steatotic liver disease, and two fibrotic and steatotic liver disease).

FIELD STRENGTH/SEQUENCE

1.5 T, single breath-hold and respiratory triggered MI-SAGE, breath-hold modified Look-Locker inversion recovery (MOLLI, T1 mapping), breath-hold gradient and spin echo (GRASE, T2 mapping), and multiple gradient echo (mGRE, T2* mapping) sequences.

ASSESSMENT

Agreement between hepatic T1, T2, and T2* estimated using MI-SAGE and conventional parametric mapping sequences was evaluated. Repeatability and reproducibility of MI-SAGE were evaluated using a same-session acquisition and second-session acquisition.

STATISTICAL TESTS

Bland-Altman analysis with bias assessment and limits of agreement (LOA) and intraclass correlation coefficients (ICC).

RESULTS

Hepatic T1, T2, and T2* estimates obtained using the MI-SAGE technique had mean biases of 72 (LOA: -22 to 166) msec, -3 (LOA: -10 to 5) msec, and 2 (LOA: -5 to 8) msec (single breath-hold) and 36 (LOA: -43 to 120) msec, -3 (LOA: -17 to 11) msec, and 4 (LOA: -3 to 11) msec (respiratory triggered), respectively, in comparison to conventional acquisitions using MOLLI, GRASE, and mGRE. All MI-SAGE estimates had strong repeatability and reproducibility (ICC > 0.72).

DATA CONCLUSION

Hepatic T1, T2, and T2* estimates obtained using an MI-SAGE technique were comparable to conventional methods, although there was a 12%/6% for breath-hold/respiratory triggered underestimation of T1 values compared to MOLLI. Both respiratory triggered and breath-hold MI-SAGE parameter maps demonstrated strong repeatability and reproducibility.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY: Stage 2.

Seven Tesla Evidence for Columnar and Rostral-Caudal Organization of the Human Periaqueductal Gray Response in the Absence of Threat: A Working Memory Study

The periaqueductal gray (PAG) is a small midbrain structure that surrounds the cerebral aqueduct, regulates brain-body communication, and is often studied for its role in "fight-or-flight" and "freezing" responses to threat. We used ultra-high-field 7 T fMRI to resolve the PAG in humans and distinguish it from the cerebral aqueduct, examining its in vivo function during a working memory task (N = 87). Both mild and moderate cognitive demands elicited spatially similar patterns of whole-brain blood oxygenation level-dependent (BOLD) response, and moderate cognitive demand elicited widespread BOLD increases above baseline in the brainstem. Notably, these brainstem increases were not significantly greater than those in the mild demand condition, suggesting that a subthreshold brainstem BOLD increase occurred for mild cognitive demand as well. Subject-specific masks were group aligned to examine PAG response. In PAG, both mild and moderate demands elicited a well-defined response in ventrolateral PAG, a region thought to be functionally related to anticipated painful threat in humans and nonhuman animals-yet, the present task posed only the most minimal (if any) "threat," with the cognitive tasks used being approximately as challenging as remembering a phone number. These findings suggest that the PAG may play a more general role in visceromotor regulation, even in the absence of threat.

Comparison of Inter-Method Agreement and Reliability for Automatic Brain Volumetry Using Three Different Clinically Available Software Packages

Background and Objectives: No comparative study has evaluated the inter-method agreement and reliability between Heuron AD and other clinically available brain volumetric software packages. Hence, we aimed to investigate the inter-method agreement and reliability of three clinically available brain volumetric software packages: FreeSurfer (FS), NeuroQuant® (NQ), and Heuron AD (HAD). Materials and Methods: In this study, we retrospectively included 78 patients who underwent conventional three-dimensional (3D) T1-weighed imaging (T1WI) to evaluate their memory impairment, including 21 with normal objective cognitive function, 24 with mild cognitive impairment, and 33 with Alzheimer's disease (AD). All 3D T1WI scans were analyzed using three different volumetric software packages. Repeated-measures analysis of variance, intraclass correlation coefficient, effect size measurements, and Bland-Altman analysis were used to evaluate the inter-method agreement and reliability. Results: The measured volumes demonstrated substantial to almost perfect agreement for most brain regions bilaterally, except for the bilateral globi pallidi. However, the volumes measured using the three software packages showed significant mean differences for most brain regions, with consistent systematic biases and wide limits of agreement in the Bland-Altman analyses. The pallidum showed the largest effect size in the comparisons between NQ and FS (5.20-6.93) and between NQ and HAD (2.01-6.17), while the cortical gray matter showed the largest effect size in the comparisons between FS and HAD (0.79-1.91). These differences and variations between the software packages were also observed in the subset analyses of 45 patients without AD and 33 patients with AD. Conclusions: Despite their favorable reliability, the software-based brain volume measurements showed significant differences and systematic biases in most regions. Thus, these volumetric measurements should be interpreted based on the type of volumetric software used, particularly for smaller structures. Moreover, users should consider the replaceability-related limitations when using these packages in real-world practice.

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.

Simultaneous T1, T2 and T2⁎ mapping of the liver with multi-shot MI-SAGE

PURPOSE

To apply multi-shot high-resolution multi inversion spin and gradient echo (MI-SAGE) acquisition for simultaneous liver T1, T2 and T2* mapping.

METHODS

Inversion prepared spin- and gradient-echo EPI was developed with ascending slice order across measurements for efficient acquisition with T1, T2, and T2⁎ weighting. Multi-shot EPI was also implemented to minimize distortion and blurring while enabling high in-plane resolution. A dictionary-matching approach was used to fit the images to quantitative parameter maps, which were compared to T1 measured by modified Look-Locker (MOLLI), T1 measured by variable flip angle (VFA), T2 measured by multiple echo time-based Half Fourier Single-shot Turbo spin-Echo (HASTE), T2 measured by radial turbo-spin-echo (rTSE) and T2⁎ measured by multiple gradient echo (MGRE) sequences.

RESULTS

The multi-shot variant of the sequence achieved higher in-plane resolution of 1.7 × 1.7 mm2 with good image quality in 28 s. Derived quantitative maps showed comparable values to conventional mapping methods. As measured in phantom and in vivo, MOLLI, MESE and MGRE give closest values to MISAGE. VFA, HASTE and rTSE show obvious overestimation.

CONCLUSIONS

The proposed multi-shot inversion prepared spin- and gradient-echo EPI sequence allows for high-resolution quantitative T1, T2 and T2 liver tissue characterization in a single breath-hold scan.

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.

Impact of susceptibility-induced distortion correction on perfusion imaging by pCASL with a segmented 3D GRASE readout

PURPOSE

The consensus for the clinical implementation of arterial spin labeling (ASL) perfusion imaging recommends a segmented 3D Gradient and Spin-Echo (GRASE) readout for optimal signal-to-noise-ratio(SNR). The correction of the associated susceptibility-induced geometric distortions has been shown to improve diagnostic precision, but its impact on ASL data has not been systematically assessed and it is not consistently part of pre-processing pipelines. Here, we investigate the effects of susceptibility-induced distortion correction on perfusion imaging by pseudo-continuous ASL (pCASL) with a segmented 3D GRASE readout.

METHODS

Data acquired from 28 women using pCASL with 3D GRASE at 3T was analyzed using three pre-processing options: without distortion correction, with distortion correction, and with spatial smoothing (without distortion correction) matched to control for blurring effects induced by distortion correction. Maps of temporal SNR (tSNR) and relative perfusion were analyzed in eight regions-of-interest (ROIs) across the brain.

RESULTS

Distortion correction significantly affected tSNR and relative perfusion across the brain. Increases in tSNR were like those produced by matched spatial smoothing in most ROIs, indicating that they were likely due to blurring effects. However, that was not the case in the frontal and temporal lobes, where we also found increased relative perfusion with distortion correction even compared with matched spatial smoothing. These effects were found in both controls and patients, with no interactions with the participant group.

CONCLUSION

Correction of Susceptibility-induced distortions significantly impacts ASL perfusion imaging using a segmented 3D GRASE readout, and this step should therefore be considered in ASL pre-processing pipelines. This is of special importance in clinical studies, reporting perfusion across ROIs defined on relatively undistorted images and when conducting group analyses requiring the alignment of images across different subjects.

Critical factors in achieving fine-scale functional MRI: Removing sources of inadvertent spatial smoothing

Ultra‐high Field (≥7T) functional magnetic resonance imaging (UHF‐fMRI) provides opportunities to resolve fine‐scale features of functional architecture such as cerebral cortical columns and layers, in vivo. While the nominal resolution of modern fMRI acquisitions may appear to be sufficient to resolve these features, several common data preprocessing steps can introduce unwanted spatial blurring, especially those that require interpolation of the data. These resolution losses can impede the detection of the fine‐scale features of interest. To examine quantitatively and systematically the sources of spatial resolution losses occurring during preprocessing, we used synthetic fMRI data and real fMRI data from the human visual cortex—the spatially interdigitated human V2 “thin” and “thick” stripes. The pattern of these cortical columns lies along the cortical surface and thus can be best appreciated using surface‐based fMRI analysis. We used this as a testbed for evaluating strategies that can reduce spatial blurring of fMRI data. Our results show that resolution losses can be mitigated at multiple points in preprocessing pathway. We show that unwanted blur is introduced at each step of volume transformation and surface projection, and can be ameliorated by replacing multi‐step transformations with equivalent single‐step transformations. Surprisingly, the simple approaches of volume upsampling and of cortical mesh refinement also helped to reduce resolution losses caused by interpolation. Volume upsampling also serves to improve motion estimation accuracy, which helps to reduce blur. Moreover, we demonstrate that the level of spatial blurring is nonuniform over the brain—knowledge which is critical for interpreting data in high‐resolution fMRI studies. Importantly, our study provides recommendations for reducing unwanted blurring during preprocessing as well as methods that enable quantitative comparisons between preprocessing strategies. These findings highlight several underappreciated sources of a spatial blur. Individually, the factors that contribute to spatial blur may appear to be minor, but in combination, the cumulative effects can hinder the interpretation of fine‐scale fMRI and the detectability of these fine‐scale features of functional architecture.

Magnetization transfer weighted EPI facilitates cortical depth determination in native fMRI space

The increased availability of ultra-high field scanners provides an opportunity to perform fMRI at sub-millimeter spatial scales and enables in vivo probing of laminar function in the human brain. In most previous studies, the definition of cortical layers, or depths, is based on an anatomical reference image that is collected by a different acquisition sequence and exhibits different geometric distortion compared to the functional images. Here, we propose to generate the anatomical image with the fMRI acquisition technique by incorporating magnetization transfer (MT) weighted imaging. Small flip angle binomial pulse trains are used as MT preparation, with a flexible duration (several to tens of milliseconds), which can be applied before each EPI segment without constraining the acquisition length (segment or slice number). The method's feasibility was demonstrated at 7T for coverage of either a small slab or the near-whole brain at 0.8 mm isotropic resolution. Tissue contrast was found to be similar to that obtained with a state-of-art anatomical reference based on MP2RAGE. This MT-weighted EPI image allows an automatic reconstruction of the cortical surface to support laminar analysis in native fMRI space, obviating the need for distortion correction and registration.

A multi-inversion multi-echo spin and gradient echo echo planar imaging sequence with low image distortion for rapid quantitative parameter mapping and synthetic image contrasts

PURPOSE

Brain imaging exams typically take 10-20 min and involve multiple sequential acquisitions. A low-distortion whole-brain echo planar imaging (EPI)-based approach was developed to efficiently encode multiple contrasts in one acquisition, allowing for calculation of quantitative parameter maps and synthetic contrast-weighted images.

METHODS

Inversion prepared spin- and gradient-echo EPI was developed with slice-order shuffling across measurements for efficient acquisition with T1 , T2 , and T 2 ∗ weighting. A dictionary-matching approach was used to fit the images to quantitative parameter maps, which in turn were used to create synthetic weighted images with typical clinical contrasts. Dynamic slice-optimized multi-coil shimming with a B0 shim array was used to reduce B0 inhomogeneity and, therefore, image distortion by >50%. Multi-shot EPI was also implemented to minimize distortion and blurring while enabling high in-plane resolution. A low-rank reconstruction approach was used to mitigate errors from shot-to-shot phase variation.

RESULTS

The slice-optimized shimming approach was combined with in-plane parallel-imaging acceleration of 4× to enable single-shot EPI with more than eight-fold distortion reduction. The proposed sequence efficiently obtained 40 contrasts across the whole-brain in just over 1 min at 1.2 × 1.2 × 3 mm resolution. The multi-shot variant of the sequence achieved higher in-plane resolution of 1 × 1 × 4 mm with good image quality in 4 min. Derived quantitative maps showed comparable values to conventional mapping methods.

CONCLUSION

The approach allows fast whole-brain imaging with quantitative parameter maps and synthetic weighted contrasts. The slice-optimized multi-coil shimming and multi-shot reconstruction approaches result in minimal EPI distortion, giving the sequence the potential to be used in rapid screening applications.

Surface-Based Connectivity Integration: An atlas-free approach to jointly study functional and structural connectivity

There has been increasing interest in jointly studying structural connectivity (SC) and functional connectivity (FC) derived from diffusion and functional MRI. Previous connectome integration studies almost exclusively required predefined atlases. However, there are many potential atlases to choose from and this choice heavily affects all subsequent analyses. To avoid such an arbitrary choice, we propose a novel atlas-free approach, named Surface-Based Connectivity Integration (SBCI), to more accurately study the relationships between SC and FC throughout the intra-cortical gray matter. SBCI represents both SC and FC in a continuous manner on the white surface, avoiding the need for prespecified atlases. The continuous SC is represented as a probability density function and is smoothed for better facilitation of its integration with FC. To infer the relationship between SC and FC, three novel sets of SC-FC coupling (SFC) measures are derived. Using data from the Human Connectome Project, we introduce the high-quality SFC measures produced by SBCI and demonstrate the use of these measures to study sex differences in a cohort of young adults. Compared with atlas-based methods, this atlas-free framework produces more reproducible SFC features and shows greater predictive power in distinguishing biological sex. This opens promising new directions for all connectomics studies.

Automated cerebral cortex segmentation based solely on diffusion tensor imaging for investigating cortical anisotropy

To extract Diffusion Tensor Imaging (DTI) parameters from the human cortex, the inner and outer boundaries of the cortex are usually defined on 3D-T1-weighted images and then applied to the co-registered DTI. However, this analysis requires the acquisition of an additional high-resolution structural image that may not be practical in various imaging studies. Here an automatic cortical boundary segmentation method was developed to work directly only on the native DTI images by using fractional anisotropy (FA) maps and mean diffusion weighted images (DWI), the latter with acceptable gray-white matter image contrast. This new method was compared to the conventional cortical segmentations generated from high-resolution T1 structural images in 5 participants. In addition, the proposed method was applied to 15 healthy young adults (10 cross-sectional, 5 test-retest) to measure FA, MD, and radiality of the primary eigenvector across the cortex on whole-brain 1.5 mm isotropic images acquired in 3.5 min at 3T. The proposed method generated reasonable segmentations of the cortical boundaries for all individuals and large proportions of the proposed method segmentations (more than 85%) were within ±1 mm from those generated with the conventional approach on higher resolution T1 structural images. Both FA (~0.15) and MD (~0.77 × 10^-3 mm²/s) extracted halfway between the cortical boundaries were relatively stable across the cortex, although focal regions such as the posterior bank of the central sulcus, anterior insula, and medial temporal lobe showed higher FA. The primary eigenvectors were primarily oriented radially to the middle cortical surface, but there were tangential orientations in the sulcal fundi as well as in the posterior bank of the central sulcus. The proposed method demonstrates the feasibility and accuracy of cortical analysis in native DTI space while avoiding the acquisition of other imaging contrasts like 3D T1-weighted scans.

Quantitative T1 mapping using multi-slice multi-shot inversion recovery EPI

An efficient multi-slice inversion–recovery EPI (MS-IR-EPI) sequence for fast, high spatial resolution, quantitative T₁ mapping is presented, using a segmented simultaneous multi-slice acquisition, combined with slice order shifting across multiple acquisitions. The segmented acquisition minimises the effective TE and readout duration compared to a single-shot EPI scheme, reducing geometric distortions to provide high quality T₁ maps with a narrow point-spread function. The precision and repeatability of MS-IR-EPI T₁ measurements are assessed using both T₁-calibrated and T₂-calibrated ISMRM/NIST phantom spheres at 3 and 7 T and compared with single slice IR and MP2RAGE methods. Magnetization transfer (MT) effects of the spectrally-selective fat-suppression (FS) pulses required for in vivo imaging are shown to shorten the measured in-vivo T₁ values. We model the effect of these fat suppression pulses on T₁ measurements and show that the model can remove their MT contribution from the measured T₁, thus providing accurate T₁ quantification. High spatial resolution T₁ maps of the human brain generated with MS-IR-EPI at 7 T are compared with those generated with the widely implemented MP2RAGE sequence. Our MS-IR-EPI sequence provides high SNR per unit time and sharper T₁ maps than MP2RAGE, demonstrating the potential for ultra-high resolution T₁ mapping and the improved discrimination of functionally relevant cortical areas in the human brain.

Efficient whole-brain tract-specific T1 mapping at 3T with slice-shuffled inversion-recovery diffusion-weighted imaging

PURPOSE

Most voxels in white matter contain multiple fiber populations with different orientations and levels of myelination. Conventional T₁ mapping measures 1 T₁ value per voxel, representing a weighted average of the multiple tract T₁ times. Inversion-recovery diffusion-weighted imaging (IR-DWI) allows the T₁ times of multiple tracts in a voxel to be disentangled, but the scan time is prohibitively long. Recently, slice-shuffled IR-DWI implementations have been proposed to significantly reduce scan time. In this work, we demonstrate that we can measure tract-specific T₁ values in the whole brain using simultaneous multi-slice slice-shuffled IR-DWI at 3T.

METHODS

We perform simulations to evaluate the accuracy and precision of our crossing fiber IR-DWI signal model for various fiber parameters. The proposed sequence and signal model are tested in a phantom consisting of crossing asparagus pieces doped with gadolinium to vary T₁ , and in 2 human subjects.

RESULTS

Our simulations show that tract-specific T₁ times can be estimated within 5% of the nominal fiber T₁ values. Tract-specific T₁ values were resolved in subvoxel 2 fiber crossings in the asparagus phantom. Tract-specific T₁ times were resolved in 2 different tract crossings in the human brain where myelination differences have previously been reported; the crossing of the cingulum and genu of the corpus callosum and the crossing of the corticospinal tract and pontine fibers.

CONCLUSION

Whole-brain tract-specific T₁ mapping is feasible using slice-shuffled IR-DWI at 3T. This technique has the potential to improve the microstructural characterization of specific tracts implicated in neurodevelopment, aging, and demyelinating disorders.

Scotopic Vision Is Selectively Processed in Thick-Type Columns in Human Extrastriate Cortex

In humans, visual stimuli can be perceived across an enormous range of light levels. Evidence suggests that different neural mechanisms process different subdivisions of this range. For instance, in the retina, stimuli presented at very low (scotopic) light levels activate rod photoreceptors, whereas cone photoreceptors are activated relatively more at higher (photopic) light levels. Similarly, different retinal ganglion cells are activated by scotopic versus photopic stimuli. However, in the brain, it remains unknown whether scotopic versus photopic information is: 1) processed in distinct channels, or 2) neurally merged. Using high-resolution functional magnetic resonance imaging at 7 T, we confirmed the first hypothesis. We first localized thick versus thin-type columns within areas V2, V3, and V4, based on photopic selectivity to motion versus color, respectively. Next, we found that scotopic stimuli selectively activated thick- (compared to thin-) type columns in V2 and V3 (in measurements of both overlap and amplitude) and V4 (based on overlap). Finally, we found stronger resting-state functional connections between scotopically dominated area MT with thick- (compared to thin-) type columns in areas V2, V3, and V4. We conclude that scotopic stimuli are processed in partially segregated parallel streams, emphasizing magnocellular influence, from retina through middle stages of visual cortex.

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

PURPOSE

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

THEORY AND METHODS

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

RESULTS

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

CONCLUSION

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

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.

Layer-Specific Contributions to Imagined and Executed Hand Movements in Human Primary Motor Cortex

The human ability to imagine motor actions without executing them (i.e., motor imagery) is crucial to a number of cognitive functions, including motor planning and learning, and has been shown to improve response times and accuracy of subsequent motor actions [1, 2]. Although these behavioral findings suggest the possibility that imagined movements directly influence primary motor cortex (M1), how this might occur remains unknown [3]. Here, we use a non-blood-oxygen-level-dependent (BOLD) method for collecting fMRI data, called vascular space occupancy (VASO) [4, 5], to measure neural activations across cortical laminae in M1 while participants either tapped their thumb and forefinger together or simply imagined doing so. We report that, whereas executed movements (i.e., finger tapping) evoked neural responses in both the superficial layers of M1 that receive cortical input and the deep layers of M1 that send output to the spinal cord to support movement, imagined movements evoked responses in superficial cortical layers only. Furthermore, we found that finger tapping preceded by both imagined and executed movements showed a reduced response in the superficial layers (repetition suppression) coupled with a heightened response in the deep layers (repetition enhancement). Taken together, our results provide evidence for a mechanism whereby imagined movements can directly affect motor performance and might explain how neural repetition effects lead to improvements in behavior (e.g., repetition priming).

Highly accelerated multishot echo planar imaging through synergistic machine learning and joint reconstruction

PURPOSE

To introduce a combined machine learning (ML)- and physics-based image reconstruction framework that enables navigator-free, highly accelerated multishot echo planar imaging (msEPI) and demonstrate its application in high-resolution structural and diffusion imaging.

METHODS

Single-shot EPI is an efficient encoding technique, but does not lend itself well to high-resolution imaging because of severe distortion artifacts and blurring. Although msEPI can mitigate these artifacts, high-quality msEPI has been elusive because of phase mismatch arising from shot-to-shot variations which preclude the combination of the multiple-shot data into a single image. We utilize deep learning to obtain an interim image with minimal artifacts, which permits estimation of image phase variations attributed to shot-to-shot changes. These variations are then included in a joint virtual coil sensitivity encoding (JVC-SENSE) reconstruction to utilize data from all shots and improve upon the ML solution.

RESULTS

Our combined ML + physics approach enabled Rinplane × multiband (MB) = 8- × 2-fold acceleration using 2 EPI shots for multiecho imaging, so that whole-brain T2 and T2 * parameter maps could be derived from an 8.3-second acquisition at 1 × 1 × 3-mm3 resolution. This has also allowed high-resolution diffusion imaging with high geometrical fidelity using 5 shots at Rinplane × MB = 9- × 2-fold acceleration. To make these possible, we extended the state-of-the-art MUSSELS reconstruction technique to simultaneous multislice encoding and used it as an input to our ML network.

CONCLUSION

Combination of ML and JVC-SENSE enabled navigator-free msEPI at higher accelerations than previously possible while using fewer shots, with reduced vulnerability to poor generalizability and poor acceptance of end-to-end ML approaches.

Brainstem neuroimaging of nociception and pain circuitries

The brainstem is known to be an important brain area for nociception and pain processing, and both relaying and coordinating signaling between the cerebrum, cerebellum, and spinal cord. Although preclinical models of pain have characterized the many roles that brainstem nuclei play in nociceptive processing, the degree to which these circuitries extend to humans is not as well known. Unfortunately, the brainstem is also a very challenging region to evaluate in humans with neuroimaging. The challenges for human brainstem imaging arise from the location of this elongated brain structure, proximity to cardiorespiratory noise sources, and the size of its constituent nuclei. These challenges can require dedicated approaches to brainstem imaging, which should be adopted when study hypotheses are focused on brainstem processing of nociception or modulation of pain perception. In fact, our review will highlight many pain neuroimaging studies that have reported some brainstem involvement in nociceptive processing and chronic pain pathology. However, we note that with recent advances in neuroimaging leading to improved spatial and temporal resolution, more studies are needed that take advantage of data collection and analysis methods focused on the challenges of brainstem neuroimaging.

Functional Involvement of Human Periaqueductal Gray and Other Midbrain Nuclei in Cognitive Control

Recent theoretical advances have motivated the hypothesis that the periaqueductal gray (PAG) participates in behaviors that involve changes in the autonomic control of visceromotor activity, including during cognitively demanding tasks. We used ultra-high-field (7 tesla) fMRI to measure human brain activity at 1.1 mm resolution while participants completed a working memory task. Consistent with prior work, participants were less accurate and responded more slowly with increasing memory load-signs of increasing task difficulty. Whole-brain fMRI analysis revealed increased activity in multiple cortical areas with increasing working memory load, including frontal and parietal cortex, dorsal cingulate, supplementary motor area, and anterior insula. Several dopamine-rich midbrain nuclei, such as the substantia nigra and ventral tegmental area, also exhibited load-dependent increases in activation. To investigate PAG involvement during cognitive engagement, we developed an automated method for segmenting and spatially normalizing the PAG. Analyses using cross-validated linear support vector machines showed that the PAG discriminated high versus low working memory load conditions with 95% accuracy in individual subjects based on activity increases in lateral and ventrolateral PAG. Effect sizes in the PAG were comparable in magnitude to those in many of the cortical areas. These findings suggest that cognitive control is not only associated with cortical activity in the frontal and parietal lobes, but also with increased activity in the subcortical PAG and other midbrain regions involved in the regulation of autonomic nervous system function.SIGNIFICANCE STATEMENT Functional neuroimaging in humans has shown that cognitive control engages multiple corticostriatal networks and brainstem nuclei, but theoretical advances suggest that the periaqueductal gray (PAG) should also be engaged during cognitively demanding tasks. Recent advances in ultra-high-field fMRI provided an opportunity to obtain the first evidence that increased activation of intermediate and rostral portions of lateral and ventrolateral PAG columns in humans is modulated by cognitive load. These findings suggest that cognitive control is not solely mediated by activity in the cortex, but that midbrain structures important for autonomic regulation also play a crucial role in higher-order cognition.

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

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

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

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

Peeking into the Black Box of Coregistration in Clinical fMRI: Which Registration Methods Are Used and How Well Do They Perform?

BACKGROUND AND PURPOSE

Interpretation of fMRI depends on accurate functional-to-structural alignment. This study explores registration methods used by FDA-approved software for clinical fMRI and aims to answer the following question: What is the degree of misalignment when registration is not performed, and how well do current registration methods perform?

MATERIALS AND METHODS

This retrospective study of presurgical fMRI for brain tumors compares nonregistered images and 5 registration cost functions: Hellinger, mutual information, normalized mutual information, correlation ratio, and local Pearson correlation. To adjudicate the accuracy of coregistration, we edge-enhanced echo-planar maps and rated them for alignment with structural anatomy. Lesion-to-activation distances were measured to evaluate the effects of different cost functions.

RESULTS

Transformation parameters were congruent among Hellinger, mutual information, normalized mutual information, and the correlation ratio but divergent from the local Pearson correlation. Edge-enhanced images validated the local Pearson correlation as the most accurate. Hellinger worsened misalignment in 59% of cases, primarily exaggerating the inferior translation; no cases were worsened by the local Pearson correlation. Three hundred twenty lesion-to-activation distances from 25 patients were analyzed among nonregistered images, Hellinger, and the local Pearson correlation. ANOVA analysis revealed significant differences in the coronal (P < .001) and sagittal (P = .04) planes. If registration is not performed, 8% of cases may have a >3-mm discrepancy and up to a 5.6-mm lesion-to-activation distance difference. If a poor registration method is used, 23% of cases may have a >3-mm discrepancy and up to a 6.9-mm difference.

CONCLUSIONS

The local Pearson correlation is a special-purpose cost function specifically designed for T2*-T1 coregistration and should be more widely incorporated into software tools as a better method for coregistration in clinical fMRI.

Dynamic 23Na MRI - A non-invasive window on neuroglial-vascular mechanisms underlying brain function

A novel magnetic resonance imaging (MRI) acquisition and reconstruction method for obtaining a series of dynamic sodium ²³Na-MRI acquisitions was designed to non-invasively assess the signal variations of brain sodium during a hand motor task in 14 healthy human volunteers on an ultra high field (7T) MR scanner. Regions undergoing activation and deactivation were identified with reference to conventional task-related BOLD functional MRI (fMRI). Activation observed in the left central regions, the supplementary motor areas and the left cerebellum induced an increase in the sodium signal observed at ultra short echo time and a decrease in the ²³Na signal observed at long echo time. Based on a simple model of two distinct sodium pools (namely, restricted and mobile sodium), the ultra short echo time measures the totality of sodium whereas the long echo time is mainly sensitive to mobile sodium. This activation pattern is consistent with previously described processes related to an influx of Na⁺ into the intracellular compartments and a moderate increase in the cerebral blood volume (CBV). In contrast, deactivation observed in the right central regions ipsilateral to the movement, the precuneus and the left cerebellum induced a slight decrease in sodium signal at ultra short echo time and an increase of sodium signal at longer echo times. This inhibitory pattern is compatible with a slight decrease in CBV and an efflux of intracellular Na⁺ to the extracellular compartments that may reflect neural dendritic spine and astrocytic shrinkage, and an increase of sodium in the extracellular fraction. In conclusion, cerebral dynamic ²³Na MRI experiments can provide access to the ionic transients following a functional task occurring within the neuro-glial-vascular ensemble. This has the potential to open up a novel non-invasive window on the mechanisms underlying brain function.

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

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

Visualizing the Human Subcortex Using Ultra-high Field Magnetic Resonance Imaging

With the recent increased availability of ultra-high field (UHF) magnetic resonance imaging (MRI), substantial progress has been made in visualizing the human brain, which can now be done in extraordinary detail. This review provides an extensive overview of the use of UHF MRI in visualizing the human subcortex for both healthy and patient populations. The high inter-subject variability in size and location of subcortical structures limits the usability of atlases in the midbrain. Fortunately, the combined results of this review indicate that a large number of subcortical areas can be visualized in individual space using UHF MRI. Current limitations and potential solutions of UHF MRI for visualizing the subcortex are also discussed.

A scalable method to improve gray matter segmentation at ultra high field MRI

High-resolution (functional) magnetic resonance imaging (MRI) at ultra high magnetic fields (7 Tesla and above) enables researchers to study how anatomical and functional properties change within the cortical ribbon, along surfaces and across cortical depths. These studies require an accurate delineation of the gray matter ribbon, which often suffers from inclusion of blood vessels, dura mater and other non-brain tissue. Residual segmentation errors are commonly corrected by browsing the data slice-by-slice and manually changing labels. This task becomes increasingly laborious and prone to error at higher resolutions since both work and error scale with the number of voxels. Here we show that many mislabeled, non-brain voxels can be corrected more efficiently and semi-automatically by representing three-dimensional anatomical images using two-dimensional histograms. We propose both a uni-modal (based on first spatial derivative) and multi-modal (based on compositional data analysis) approach to this representation and quantify the benefits in 7 Tesla MRI data of nine volunteers. We present an openly accessible Python implementation of these approaches and demonstrate that editing cortical segmentations using two-dimensional histogram representations as an additional post-processing step aids existing algorithms and yields improved gray matter borders. By making our data and corresponding expert (ground truth) segmentations openly available, we facilitate future efforts to develop and test segmentation algorithms on this challenging type of data.

Analysis strategies for high-resolution UHF-fMRI data

Functional MRI (fMRI) benefits from both increased sensitivity and specificity with increasing magnetic field strength, making it a key application for Ultra-High Field (UHF) MRI scanners. Most UHF-fMRI studies utilize the dramatic increases in sensitivity and specificity to acquire high-resolution data reaching sub-millimeter scales, which enable new classes of experiments to probe the functional organization of the human brain. This review article surveys advanced data analysis strategies developed for high-resolution fMRI at UHF. These include strategies designed to mitigate distortion and artifacts associated with higher fields in ways that attempt to preserve spatial resolution of the fMRI data, as well as recently introduced analysis techniques that are enabled by these extremely high-resolution data. Particular focus is placed on anatomically-informed analyses, including cortical surface-based analysis, which are powerful techniques that can guide each step of the analysis from preprocessing to statistical analysis to interpretation and visualization. New intracortical analysis techniques for laminar and columnar fMRI are also reviewed and discussed. Prospects for single-subject individualized analyses are also presented and discussed. Altogether, there are both specific challenges and opportunities presented by UHF-fMRI, and the use of proper analysis strategies can help these valuable data reach their full potential.

Inter-subject Registration of Functional Images: Do We Need Anatomical Images?

In Echo-Planar Imaging (EPI)-based Magnetic Resonance Imaging (MRI), inter-subject registration typically uses the subject's T1-weighted (T1w) anatomical image to learn deformations of the subject's brain onto a template. The estimated deformation fields are then applied to the subject's EPI scans (functional or diffusion-weighted images) to warp the latter to a template space. Historically, such indirect T1w-based registration was motivated by the lack of clear anatomical details in low-resolution EPI images: a direct registration of the EPI scans to template space would be futile. A central prerequisite in such indirect methods is that the anatomical (aka the T1w) image of each subject is well aligned with their EPI images via rigid coregistration. We provide experimental evidence that things have changed: nowadays, there is a decent amount of anatomical contrast in high-resolution EPI data. That notwithstanding, EPI distortions due to B0 inhomogeneities cannot be fully corrected. Residual uncorrected distortions induce non-rigid deformations between the EPI scans and the same subject's anatomical scan. In this manuscript, we contribute a computationally cheap pipeline that leverages the high spatial resolution of modern EPI scans for direct inter-subject matching. Our pipeline is direct and does not rely on the T1w scan to estimate the inter-subject deformation. Results on a large dataset show that this new pipeline outperforms the classical indirect T1w-based registration scheme, across a variety of post-registration quality-assessment metrics including: Normalized Mutual Information, relative variance (variance-to-mean ratio), and to a lesser extent, improved peaks of group-level General Linear Model (GLM) activation maps.

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

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

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

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

Optimized inversion-time schedules for quantitative T1 measurements based on high-resolution multi-inversion EPI

PURPOSE

Demonstrate an optimized multi-inversion echo-planar imaging technique to accelerate quantitative T₁ mapping by judicious selection of inversion times for each slice.

METHODS

Slice ordering is optimized to maximize discrimination between tissues with different T₁ values. The optimized slice orderings are tested in the International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology phantom and compared with an unoptimized 21-measurement acquisition. The utility of the method is demonstrated in a healthy subject in vivo at 3 T and validated with a gold-standard inversion-recovery sequence. The in vivo precision of our technique was tested by repeated scans of the same subject within a scan session and across scan sessions, occurring 28 days apart.

RESULTS

Phantom measurements yielded good agreement (R²  = 0.99) between the T₁ estimates from the proposed optimized protocol, reference values from the National Institute of Standards and Technology phantom and gold-standard inversion-recovery values, as well as a negligible estimation bias that was slightly lower than that from the unoptimized 21-measurement protocol (0.74 versus 19 ms). The range of values for the scan-rescan coefficient of variation was 0.86 to 0.93 (within session) and 0.83 to 0.92 (across sessions) across all scan durations tested.

CONCLUSIONS

Optimized slice orderings allow faster quantitative T₁ mapping. The optimized sequence yielded accurate and precise T₁ maps. Magn Reson Med 79:2101-2112, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Methods for cleaning the BOLD fMRI signal

Blood oxygen-level-dependent functional magnetic resonance imaging (BOLD fMRI) has rapidly become a popular technique for the investigation of brain function in healthy individuals, patients as well as in animal studies. However, the BOLD signal arises from a complex mixture of neuronal, metabolic and vascular processes, being therefore an indirect measure of neuronal activity, which is further severely corrupted by multiple non-neuronal fluctuations of instrumental, physiological or subject-specific origin. This review aims to provide a comprehensive summary of existing methods for cleaning the BOLD fMRI signal. The description is given from a methodological point of view, focusing on the operation of the different techniques in addition to pointing out the advantages and limitations in their application. Since motion-related and physiological noise fluctuations are two of the main noise components of the signal, techniques targeting their removal are primarily addressed, including both data-driven approaches and using external recordings. Data-driven approaches, which are less specific in the assumed model and can simultaneously reduce multiple noise fluctuations, are mainly based on data decomposition techniques such as principal and independent component analysis. Importantly, the usefulness of strategies that benefit from the information available in the phase component of the signal, or in multiple signal echoes is also highlighted. The use of global signal regression for denoising is also addressed. Finally, practical recommendations regarding the optimization of the preprocessing pipeline for the purpose of denoising and future venues of research are indicated. Through the review, we summarize the importance of signal denoising as an essential step in the analysis pipeline of task-based and resting state fMRI studies.

Magnetic resonance imaging of granular materials

Magnetic Resonance Imaging (MRI) has become one of the most important tools to screen humans in medicine; virtually every modern hospital is equipped with a Nuclear Magnetic Resonance (NMR) tomograph. The potential of NMR in 3D imaging tasks is by far greater, but there is only "a handful" of MRI studies of particulate matter. The method is expensive, time-consuming, and requires a deep understanding of pulse sequences, signal acquisition, and processing. We give a short introduction into the physical principles of this imaging technique, describe its advantages and limitations for the screening of granular matter, and present a number of examples of different application purposes, from the exploration of granular packing, via the detection of flow and particle diffusion, to real dynamic measurements. Probably, X-ray computed tomography is preferable in most applications, but fast imaging of single slices with modern MRI techniques is unmatched, and the additional opportunity to retrieve spatially resolved flow and diffusion profiles without particle tracking is a unique feature.

HIgh b-value and high Resolution Integrated Diffusion (HIBRID) imaging

The parameter selection for diffusion MRI experiments is dominated by the "k-q tradeoff" whereby the Signal to Noise Ratio (SNR) of the images is traded for either high spatial resolution (determined by the maximum k-value collected) or high diffusion sensitivity (effected by b-value or the q vector) but usually not both. Furthermore, different brain regions (such as gray matter and white matter) likely require different tradeoffs between these parameters due to the size of the structures to be visualized or the length-scale of the microstructure being probed. In this case, it might be advantageous to combine information from two scans - a scan with high q but low k (high angular resolution in diffusion but low spatial resolution in the image domain) to provide maximal information about white matter fiber crossing, and one low q but high k (low angular resolution but high spatial resolution) for probing the cortex. In this study, we propose a method, termed HIgh b-value and high Resolution Integrated Diffusion (HIBRID) imaging, for acquiring and combining the information from these two complementary types of scan with the goal of studying diffusion in the cortex without compromising white matter fiber information. The white-gray boundary and pial surface obtained from anatomical scans are incorporated as prior information to guide the fusion. We study the complementary advantages of the fused datasets, and assess the quality of the HIBRID data compared to either alone.

High resolution data analysis strategies for mesoscale human functional MRI at 7 and 9.4T

The advent of ultra-high field functional magnetic resonance imaging (fMRI) has greatly facilitated submillimeter resolution acquisitions (voxel volume below (1 mm³)), allowing the investigation of cortical columns and cortical depth dependent (i.e. laminar) structures in the human brain. Advanced data analysis techniques are essential to exploit the information in high resolution functional measures. In this article, we use recent, exemplary 9.4 T human functional and anatomical data to review the advantages and disadvantages of (1) pooling high resolution data across regions of interest for cortical depth profile analysis, (2) pooling across cortical depths for mapping patches of cortex while discarding depth-dependent (i.e. columnar) effects, and (3) isotropic sampling without pooling to assess individual voxel’s responses. A set of cortical depth meshes may be a solution to sampling information tangentially while keeping correspondence across depths. For quantitative analysis of the spatial organization in fine-grained structures, a cortical grid approach is advantageous. We further extend this general framework by combining it with a previously introduced cortical layer volume-preserving (equi-volume) approach. This framework can readily accommodate the research questions which allow for spatial smoothing within or across layers. We demonstrate and discuss that equi-volume sampling yields a slight advantage over equidistant sampling given the current limitations of fMRI voxel size, participant motion, coregistration and segmentation. Our 9.4 T human anatomical and functional data indicate the advantage over lower fields including 7 T and demonstrate the practical applicability of T₂^* and T₂-weighted fMRI acquisitions.

•

High resolution regular cortical grids are advantageous for local applications.

•

Equi-volume sampling is slightly advantageous over equidistant sampling in-vivo.

•

Isotropic submillimeter cortical sampling without spatial pooling requires high SNR.

•

9.4 T human T2 and T2* BOLD fMRI are practically feasible and provide high SNR.

•

9.4 T T2*-weighted 0.35 mm iso. res. anatomical images for laminar contrast in vivo.

The impact of ultra-high field MRI on cognitive and computational neuroimaging

The ability to measure functional brain responses non-invasively with ultra high field MRI (7 T and above) represents a unique opportunity in advancing our understanding of the human brain. Compared to lower fields (3 T and below), ultra high field MRI has an increased sensitivity, which can be used to acquire functional images with greater spatial resolution, and greater specificity of the blood oxygen level dependent (BOLD) signal to the underlying neuronal responses. Together, increased resolution and specificity enable investigating brain functions at a submillimeter scale, which so far could only be done with invasive techniques. At this mesoscopic spatial scale, perception, cognition and behavior can be probed at the level of fundamental units of neural computations, such as cortical columns, cortical layers, and subcortical nuclei. This represents a unique and distinctive advantage that differentiates ultra high from lower field imaging and that can foster a tighter link between fMRI and computational modeling of neural networks. So far, functional brain mapping at submillimeter scale has focused on the processing of sensory information and on well-known systems for which extensive information is available from invasive recordings in animals. It remains an open challenge to extend this methodology to uniquely human functions and, more generally, to systems for which animal models may be problematic. To succeed, the possibility to acquire high-resolution functional data with large spatial coverage, the availability of computational models of neural processing as well as accurate biophysical modeling of neurovascular coupling at mesoscopic scale all appear necessary.

Challenges and opportunities for brainstem neuroimaging with ultrahigh field MRI

The human brainstem plays a central role in connecting the cerebrum, the cerebellum and the spinal cord to one another, hosting relay nuclei for afferent and efferent signaling, and providing source nuclei for several neuromodulatory systems that impact central nervous system function. While the investigation of the brainstem with functional or structural magnetic resonance imaging has been hampered for years due to this brain structure's physiological and anatomical characteristics, the field has seen significant advances in recent years thanks to the broader adoption of ultrahigh-field (UHF) MRI scanning. In the present review, we focus on the advantages offered by UHF in the context of brainstem imaging, as well as the challenges posed by the investigation of this complex brain structure in terms of data acquisition and analysis. We also illustrate how UHF MRI can shed new light on the neuroanatomy and neurophysiology underlying different brainstem-based circuitries, such as the central autonomic network and neurotransmitter/neuromodulator systems, discuss existing and foreseeable clinical applications to better understand diseases such as chronic pain and Parkinson's disease, and explore promising future directions for further improvements in brainstem imaging using UHF MRI techniques.