How to choose the right MR sequence for your research question at 7 T and above?

Population Neuroscience: Principles and Advances

In population neuroscience, three disciplines come together to advance our knowledge of factors that shape the human brain: neuroscience, genetics, and epidemiology (Paus, Human Brain Mapping 31:891-903, 2010). Here, I will come back to some of the background material reviewed in more detail in our previous book (Paus, Population Neuroscience, 2013), followed by a brief overview of current advances and challenges faced by this integrative approach.

Empirical transmit field bias correction of T1w/T2w myelin maps

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

7T versus 3T MRI in the presurgical evaluation of patients with drug-resistant epilepsy

BACKGROUND AND PURPOSE

MRI has a crucial role in presurgical evaluation of drug-resistant focal epilepsy patients. Whether and how much 7T MRI further improves presurgical diagnosis compared to standard of care 3T MRI remains to be established. We investigate the added value 7T MRI offers in surgical candidates with remaining clinical uncertainty after 3T MRI.

METHODS

7T brain MRI was obtained on sequential patients with drug-resistant focal epilepsy undergoing presurgical evaluation at a comprehensive epilepsy center, including patients with and without suspected lesions on standard 3T MRI. Clinical information and 3T images informed the interpretation of 7T images. Detection of a new lesion on 7T or better characterization of a suspected lesion was considered to add value to the presurgical workup.

RESULTS

Interpretable 7T MRI was acquired in 19 patients. 7T MRI identified a lesion relevant to the seizures in three of eight patients (38%) without a lesion on 3T MRI; no lesion in 7/11 patients (64%) with at least one suspected lesion on 3T MRI, contributing to the final classification of all seven as nonlesional; and confirmed and better characterized the lesion suspected at 3T MR in the remaining 4/11 patients.

CONCLUSIONS

7T MRI detected new lesions in over a third of 3T MRI nonlesional patients, confirmed and better characterized a 3T suspected lesion in one third of patients, and helped exclude a 3T suspected lesion in the remainder. Our initial experience suggests that 7T MRI adds value to surgical planning by improving detection and characterization of suspected brain lesions in drug-resistant focal epilepsy patients.

Functional magnetic resonance imaging responses during perceptual decision-making at 3 and 7 T in human cortex, striatum, and brainstem

While functional magnetic resonance imaging (fMRI) at ultra‐high field (7 T) promises a general increase in sensitivity compared to lower field strengths, the benefits may be most pronounced for specific applications. The current study aimed to evaluate the relative benefit of 7 over 3 T fMRI for the assessment of responses evoked in different brain regions by a well‐controlled cognitive task. At 3 and 7 T, the same participants made challenging perceptual decisions about visual motion combined with monetary rewards for correct choices. Previous work on this task has extensively characterized the underlying cognitive computations and single‐cell responses in cortical and subcortical structures. We quantified the evoked fMRI responses in extrastriate visual cortical areas, the striatum, and the brainstem during the decision interval and the post‐feedback interval of the task. The dependence of response amplitudes on field strength during the decision interval differed between cortical, striatal, and brainstem regions, with a generally bigger 7 versus 3 T benefit in subcortical structures. We also found stronger responses during relatively easier than harder decisions at 7 T for dopaminergic midbrain nuclei, in line with reward expectation. Our results demonstrate the potential of 7 T fMRI for illuminating the contribution of small brainstem nuclei to the orchestration of cognitive computations in the human brain.

Imaging human engrams using 7 Tesla magnetic resonance imaging

The investigation of the physical traces of memories (engrams) has made significant progress in the last decade due to optogenetics and fluorescent cell tagging applied in rodents. Engram cells were identified. The ablation of engram cells led to the loss of the associated memory, silent memories were reactivated, and artificial memories were implanted in the brain. Human engram research lags behind engram research in rodents due to methodological and ethical constraints. However, advances in multivariate analysis techniques of functional magnetic resonance imaging (fMRI) data and machine learning algorithms allowed the identification of stable engram patterns in humans. In addition, MRI scanners with an ultrahigh field strength of 7 Tesla (T) have left their prototype state and became more common around the world to assist human engram research. Although most engram research in humans is still being performed with a field strength of 3T, fMRI at 7T will push engram research. Here, we summarize the current state and findings of human engram research and discuss the advantages and disadvantages of applying 7 versus 3T fMRI to image human memory traces.

7 Tesla and Beyond: Advanced Methods and Clinical Applications in Magnetic Resonance Imaging

ABSTRACT

Ultrahigh magnetic fields offer significantly higher signal-to-noise ratio, and several magnetic resonance applications additionally benefit from a higher contrast-to-noise ratio, with static magnetic field strengths of B0 ≥ 7 T currently being referred to as ultrahigh fields (UHFs). The advantages of UHF can be used to resolve structures more precisely or to visualize physiological/pathophysiological effects that would be difficult or even impossible to detect at lower field strengths. However, with these advantages also come challenges, such as inhomogeneities applying standard radiofrequency excitation techniques, higher energy deposition in the human body, and enhanced B0 field inhomogeneities. The advantages but also the challenges of UHF as well as promising advanced methodological developments and clinical applications that particularly benefit from UHF are discussed in this review article.

Cerebral white matter vasculature: still uncharted?

White matter vasculature plays a major role in the pathophysiology of permanent neurological deficits following a stroke or progressive cognitive alteration related to small vessel disease. Thus, knowledge of the complex vascularization and functional aspects of the deep white matter territories is paramount to comprehend clinical manifestations of brain ischemia. This review provides a structured presentation of the existing knowledge of the vascularization of the human cerebral white matter from seminal historical studies to the current literature. First, we revisit the highlights of prenatal development of the endoparenchymal telencephalic vascular system that are crucial for the understanding of vessel organization in the adult. Second, we reveal the tangled history of debates on the existence, clinical significance, and physiological role of leptomeningeal anastomoses. Then, we present how conceptions on white matter vascularization transitioned from the mixed ventriculopetal/ventriculofugal theory, in which a low-flow area was interposed in between concurrent arterial flows, to the purely ventriculopetal theory. The latter model explains variable white matter sensitivity to ischemia by various organizations of ventriculopetal vessel terminals having different origin/length properties and interconnection patterns. Next, arteries supplying primarily the white matter are described according to their length and overall structure. Furthermore, the known distribution territories, to date, are studied in relation to primary anatomical structures of the human cerebral white matter, emphasizing the sparsity of the "ground-truth" data available in the literature. Finally, the implications for both large vessel occlusion and chronic small vessel disease are discussed, as well as the insights from neuroimaging. All things considered, we identify the need for further research on deep white matter vascularization, especially regarding the arterial supply of white matter fiber tracts.

Can 7T MPRAGE match MP2RAGE for gray-white matter contrast?

Ultra-High Field (UHF) MRI provides a significant increase in Signal-to-Noise Ratio (SNR) and gains in contrast weighting in several functional and structural acquisitions. Unfortunately, an increase in field strength also induces non-uniformities in the transmit field (B1+) that can compromise image contrast non-uniformly. The MPRAGE is one of the most common T1 weighted (T1w) image acquisitions for structural imaging. It provides excellent contrast between gray and white matter and is widely used for brain segmentation. At 7T, the signal non-uniformities tend to complicate this and therefore, the self-bias-field corrected MP2RAGE is often used there. In both MPRAGE and MP2RAGE, more homogeneous image contrast can be achieved with adiabatic pulses, like the TR-FOCI inversion pulse, or special pulse design on parallel transmission systems, like Universal Pulses (UP). In the present study, we investigate different strategies to improve the bias-field for MPRAGE at 7T, comparing the contrast and GM/WM segmentability against MP2RAGE. The higher temporal efficiency of MPRAGE combined with the potential of the user-friendly UPs was the primary motivation for this MPRAGE-MP2RAGE comparison. We acquired MPRAGE data in six volunteers, adding a k-space shutter to reduce scan time, a kt-point UP approach for homogeneous signal excitation, and a TR-FOCI pulse for homogeneous inversion. Our results show remarkable signal contrast improvement throughout the brain, including regions of low B1⁺ such as the cerebellum. The improvements in the MPRAGE were largest following the introduction of the UPs. In addition to the CNR, both SNR and GM/WM segmentability were also assessed. Among the MPRAGEs, the combined strategy (UP + TR-FOCI) yielded highest SNR and showed highest spatial similarity between GM segments to the MP2RAGE. Interestingly, the distance between gray and white matter peaks in the intensity histograms did not increase, as better pulses and higher SNR especially benefitted the (cerebellar) gray matter. Overall, the gray-white matter contrast from MP2RAGE is higher, with higher CNR and higher intensity peak distances, even when scaled to scan time. Hence, the extra acquisition time for MP2RAGE is justified by the improved segmentability.

QSM reconstruction challenge 2.0: A realistic in silico head phantom for MRI data simulation and evaluation of susceptibility mapping procedures

PURPOSE

To create a realistic in silico head phantom for the second QSM reconstruction challenge and for future evaluations of processing algorithms for QSM.

METHODS

We created a digital whole-head tissue property phantom by segmenting and postprocessing high-resolution (0.64 mm isotropic), multiparametric MRI data acquired at 7 T from a healthy volunteer. We simulated the steady-state magnetization at 7 T using a Bloch simulator and mimicked a Cartesian sampling scheme through Fourier-based processing. Computer code for generating the phantom and performing the MR simulation was designed to facilitate flexible modifications of the phantom in the future, such as the inclusion of pathologies as well as the simulation of a wide range of acquisition protocols. Specifically, the following parameters and effects were implemented: TR and TE, voxel size, background fields, and RF phase biases. Diffusion-weighted imaging phantom data are provided, allowing future investigations of tissue-microstructure effects in phase and QSM algorithms.

RESULTS

The brain part of the phantom featured realistic morphology with spatial variations in relaxation and susceptibility values similar to the in vivo setting. We demonstrated some of the phantom's properties, including the possibility of generating phase data with nonlinear evolution over TE due to partial-volume effects or complex distributions of frequency shifts within the voxel.

CONCLUSION

The presented phantom and computer programs are publicly available and may serve as a ground truth in future assessments of the faithfulness of quantitative susceptibility reconstruction algorithms.

A plug-and-play, lightweight, single-axis gradient insert design for increasing spatiotemporal resolution in echo planar imaging-based brain imaging

The goal of this study was to introduce and evaluate the performance of a lightweight, high-performance, single-axis (z-axis) gradient insert design primarily intended for high-resolution functional magnetic resonance imaging, and aimed at providing both ease of use and a boost in spatiotemporal resolution. The optimal winding positions of the coil were obtained using a genetic algorithm with a cost function that balanced gradient performance (minimum 0.30 mT/m/A) and field linearity (≥16 cm linear region). These parameters were verified using field distribution measurements by B₀ -mapping. The correction of geometrical distortions was performed using theoretical field distribution of the coil. Simulations and measurements were performed to investigate the echo planar imaging echo-spacing reduction due to the improved gradient performance. The resulting coil featured a 16-cm linear region, a weight of 45 kg, an installation time of 15 min, and a maximum gradient strength and slew rate of 200 mT/m and 1300 T/m/s, respectively, when paired with a commercially available gradient amplifier (940 V/630 A). The field distribution measurements matched the theoretically expected field. By utilizing the theoretical field distribution, geometrical distortions were corrected to within 6% of the whole-body gradient reference image in the target region. Compared with a whole-body gradient set, a maximum reduction in echo-spacing of a factor of 2.3 was found, translating to a 344 μs echo-spacing, for a field of view of 192 mm, a receiver bandwidth of 920 kHz and a gradient amplitude of 112 mT/m. We present a lightweight, single-axis gradient insert design that can provide high gradient performance and an increase in spatiotemporal resolution with correctable geometrical distortions while also offering a short installation time of less than 15 min and minimal system modifications.

Estimation of laminar BOLD activation profiles using deconvolution with a physiological point spread function

BACKGROUND

The specificity of gradient echo (GE)-BOLD laminar fMRI activation profiles is degraded by intracortical veins that drain blood from lower to upper cortical layers, propagating activation signal in the same direction. This work describes an approach to obtain layer specific profiles by deconvolving the measured profiles with a physiological Point Spread Function (PSF).

NEW METHOD

It is shown that the PSF can be characterised by a TE-dependent peak to tail (p2t) value that is independent of cortical depth and can be estimated by simulation. An experimental estimation of individual p2t values and the sensitivity of the deconvolved profiles to variations in p2t is obtained using laminar data measured with a multi-echo 3D-FLASH sequence. These profiles are echo time dependent, but the underlying neuronal response is the same, allowing a data-based estimation of the PSF.

RESULTS

The deconvolved profiles are highly similar to the gold-standard obtained from extremely high resolution 3D-EPI data, for a range of p2t values of 5-9, which covers both the empirically determined value (6.8) and the value obtained by simulation (6.3). -Comparison with Existing Method(s) Corrected profiles show a flatter shape across the cortex and a high level of similarity with the gold-standard, defined as a subset of profiles that are unaffected by intracortical veins.

CONCLUSIONS

We conclude that deconvolution is a robust approach for removing the effect of signal propagation through intracortical veins. This makes it possible to obtain profiles with high laminar specificity while benefitting from the higher efficiency of GE-BOLD sequences.

Realistic head-shaped phantom with brain-mimicking metabolites for 7 T spectroscopy and spectroscopic imaging

PURPOSE

Moving to ultra-high fields (≥7 T), the inhomogeneity of both RF (B1 ) and static (B0 ) magnetic fields increases, which further motivates us to design a realistic head-shaped phantom, especially for spectroscopic imaging. Such phantoms provide images similar to the human brain and serve as a reliable tool for developing and examining methods in MRI. This study aims to develop and characterize a realistic head-shaped phantom filled with brain-mimicking metabolites for MRS and magnetic resonance spectroscopic imaging in a 7 T MRI scanner.

METHODS

A 3D head-shaped container with three sections-mimicking brain, muscle and precranial lipid-was constructed. The phantom was designed to provide robustness to heating, mechanical damage and leakage, with easy refilling. The head's shape and the agarose mixture were optimized to provide B0 and B1 distributions and T1 /T2 relaxation values similar to those of human brain. Eight brain-tissue-mimicking metabolites were included for spectroscopy. The phantom was evaluated for localized spectroscopy, fast spectroscopic imaging and fat suppression.

RESULTS

The B0 and B1 maps showed distribution similar to that of human brain, with increased B0 inhomogeneity near the nasal and ear areas and reduced B1 in the temporal lobe and brain stem regions, as expected in vivo. The metabolites' concentrations were verified by single-voxel spectroscopy, showing an average deviation of 11%. Fast spectroscopic imaging and imaging with fat suppression were demonstrated.

CONCLUSION

A 3D head-shaped phantom for human brain imaging and spectroscopic imaging in 7 T MRI was demonstrated, making it a realistic phantom for methodology development at 7 T.

Ultra-high field fMRI reveals origins of feedforward and feedback activity within laminae of human ocular dominance columns

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

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

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

Methodological Considerations for Neuroimaging in Deep Brain Stimulation of the Subthalamic Nucleus in Parkinson's Disease Patients

Deep brain stimulation (DBS) of the subthalamic nucleus is a neurosurgical intervention for Parkinson's disease patients who no longer appropriately respond to drug treatments. A small fraction of patients will fail to respond to DBS, develop psychiatric and cognitive side-effects, or incur surgery-related complications such as infections and hemorrhagic events. In these cases, DBS may require recalibration, reimplantation, or removal. These negative responses to treatment can partly be attributed to suboptimal pre-operative planning procedures via direct targeting through low-field and low-resolution magnetic resonance imaging (MRI). One solution for increasing the success and efficacy of DBS is to optimize preoperative planning procedures via sophisticated neuroimaging techniques such as high-resolution MRI and higher field strengths to improve visualization of DBS targets and vasculature. We discuss targeting approaches, MRI acquisition, parameters, and post-acquisition analyses. Additionally, we highlight a number of approaches including the use of ultra-high field (UHF) MRI to overcome limitations of standard settings. There is a trade-off between spatial resolution, motion artifacts, and acquisition time, which could potentially be dissolved through the use of UHF-MRI. Image registration, correction, and post-processing techniques may require combined expertise of traditional radiologists, clinicians, and fundamental researchers. The optimization of pre-operative planning with MRI can therefore be best achieved through direct collaboration between researchers and clinicians.

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.

Encoding schemes in somatosensation: From micro- to meta-topography

Encoding schemes are systematic large-scale arrangements that convert incoming sensory information into a format required for further information processing. The increased spatial resolution of brain images obtained with ultra-high field magnetic resonance imaging at 7 T (7T-MRI) and above increases the granularity and precision of processing units that mediate the link between neuronal encoding and functional readouts. Here, these new developments are reviewed with a focus on human tactile encoding schemes derived from small-scale processing units (in the order of 0.5-5 mm) that are relevant for theoretical and practical concepts of somatosensory encoding and cortical plasticity. Precisely, we review recent approaches to characterize meso-scale maps, layer units, and cortical fields in the sensorimotor cortex of the living human brain and discuss their impact on theories of perception, motor control, topographic encoding, and cortical plasticity. Finally, we discuss concepts on the integration of small-scale processing units into functional networks that span multiple topographic maps and multiple cortical areas. Novel research areas are highlighted that may help to bridge the gap between cortical microstructure and meta-topographic models on brain anatomy and function.

Sharpness in motion corrected quantitative imaging at 7T

Sub-millimeter imaging at 7T has opened new possibilities for qualitatively and quantitatively studying brain structure as it evolves throughout the life span. However, subject motion introduces image blurring on the order of magnitude of the spatial resolution and is thus detrimental to image quality. Such motion can be corrected for, but widespread application has not yet been achieved and quantitative evaluation is lacking. This raises a need to quantitatively measure image sharpness throughout the brain. We propose a method to quantify sharpness of brain structures at sub-voxel resolution, and use it to assess to what extent limited motion is related to image sharpness. The method was evaluated in a cohort of 24 healthy volunteers with a wide and uniform age range, aiming to arrive at results that largely generalize to larger populations. Using 3D fat-excited motion navigators, quantitative R₁, R₂^* and Quantitative Susceptibility Maps and T₁-weighted images were retrospectively corrected for motion. Sharpness was quantified in all modalities for selected regions of interest (ROI) by fitting the sigmoidally shaped error function to data within locally homogeneous clusters. A strong, almost linear correlation between motion and sharpness improvement was observed, and motion correction significantly improved sharpness. Overall, the Full Width at Half Maximum reduced from 0.88 mm to 0.70 mm after motion correction, equivalent to a 2.0 times smaller voxel volume. Motion and sharpness were not found to correlate with the age of study participants. We conclude that in our data, motion correction using fat navigators is overall able to restore the measured sharpness to the imaging resolution, irrespective of the amount of motion observed during scanning.

Direct visualization and characterization of the human zona incerta and surrounding structures

The zona incerta (ZI) is a small gray matter region of the deep brain first identified in the 19th century, yet direct in vivo visualization and characterization has remained elusive. Noninvasive detection of the ZI and surrounding region could be critical to further our understanding of this widely connected but poorly understood deep brain region and could contribute to the development and optimization of neuromodulatory therapies. We demonstrate that high resolution (submillimetric) longitudinal (T1) relaxometry measurements at high magnetic field strength (7 T) can be used to delineate the ZI from surrounding white matter structures, specifically the fasciculus cerebellothalamicus, fields of Forel (fasciculus lenticularis, fasciculus thalamicus, and field H), and medial lemniscus. Using this approach, we successfully derived in vivo estimates of the size, shape, location, and tissue characteristics of substructures in the ZI region, confirming observations only previously possible through histological evaluation that this region is not just a space between structures but contains distinct morphological entities that should be considered separately. Our findings pave the way for increasingly detailed in vivo study and provide a structural foundation for precise functional and neuromodulatory investigation.

Longitudinal Automatic Segmentation of Hippocampal Subfields (LASHiS) using multi-contrast MRI

The volumetric and morphometric examination of hippocampus formation subfields in a longitudinal manner using in vivo MRI could lead to more sensitive biomarkers for neuropsychiatric disorders and diseases including Alzheimer's disease, as the anatomical subregions are functionally specialised. Longitudinal processing allows for increased sensitivity due to reduced confounds of inter-subject variability and higher effect-sensitivity than cross-sectional designs. We examined the performance of a new longitudinal pipeline (Longitudinal Automatic Segmentation of Hippocampus Subfields [LASHiS]) against three freely available, published approaches. LASHiS automatically segments hippocampus formation subfields by propagating labels from cross-sectionally labelled time point scans using joint-label fusion to a non-linearly realigned 'single subject template', where image segmentation occurs free of bias to any individual time point. Our pipeline measures tissue characteristics available in in vivo high-resolution MRI scans, at both clinical (3 ​T) and ultra-high field strength (7 ​T) and differs from previous longitudinal segmentation pipelines in that it leverages multi-contrast information in the segmentation process. LASHiS produces robust and reliable automatic multi-contrast segmentations of hippocampus formation subfields, as measured by higher volume similarity coefficients and Dice coefficients for test-retest reliability and robust longitudinal Bayesian Linear Mixed Effects results at 7 ​T, while showing sound results at 3 ​T. All code for this project including the automatic pipeline is available at https://github.com/CAIsr/LASHiS.

The pulvinar nucleus and antidepressant treatment: dynamic modeling of antidepressant response and remission with ultra-high field functional MRI

Functional magnetic resonance imaging (fMRI) successfully disentangled neuronal pathophysiology of major depression (MD), but only a few fMRI studies have investigated correlates and predictors of remission. Moreover, most studies have used clinical outcome parameters from two time points, which do not optimally depict differential response times. Therefore, we aimed to detect neuronal correlates of response and remission in an antidepressant treatment study with 7 T fMRI, potentially harnessing advances in detection power and spatial specificity. Moreover, we modeled outcome parameters from multiple study visits during a 12-week antidepressant fMRI study in 26 acute (aMD) patients compared to 36 stable remitted (rMD) patients and 33 healthy control subjects (HC). During an electrical painful stimulation task, significantly higher baseline activity in aMD compared to HC and rMD in the medial thalamic nuclei of the pulvinar was detected (p = 0.004, FWE-corrected), which was reduced by treatment. Moreover, clinical response followed a sigmoid function with a plateau phase in the beginning, a rapid decline and a further plateau at treatment end. By modeling the dynamic speed of response with fMRI-data, perigenual anterior cingulate activity after treatment was significantly associated with antidepressant response (p < 0.001, FWE-corrected). Temporoparietal junction (TPJ) baseline activity significantly predicted non-remission after 2 antidepressant trials (p = 0.005, FWE-corrected). The results underline the importance of the medial thalamus, attention networks in MD and antidepressant treatment. Moreover, by using a sigmoid model, this study provides a novel method to analyze the dynamic nature of response and remission for future trials.

MP2RAGEME: T1 , T2 * , and QSM mapping in one sequence at 7 tesla

Quantitative magnetic resonance imaging generates images of meaningful physical or chemical variables measured in physical units that allow quantitative comparisons between tissue regions and among subjects scanned at the same or different sites. Here, we show that we can acquire quantitative T₁, T₂^*, and quantitative susceptibility mapping (QSM) information in a single acquisition, using a multi‐echo (ME) extension of the second gradient‐echo image of the MP2RAGE sequence. This combination is called MP2RAGE ME, or MP2RAGEME. The simultaneous acquisition results in large time savings, perfectly coregistered data, and minimal image quality differences compared to separately acquired data. Following a correction for residual transmit B₁⁺‐sensitivity, quantitative T₁, T₂^*, and QSM values were in excellent agreement with those obtained from separately acquired, also B₁⁺‐corrected, MP2RAGE data and ME gradient echo data. The quantitative values from reference regions of interests were also in very good correspondence with literature values. From the MP2RAGEME data, we further derived a multiparametric cortical parcellation, as well as a combined arterial and venous map. In sum, our MP2RAGEME sequence has the benefit in large time savings, perfectly coregistered data and minor image quality differences.

Pros and cons of ultra-high-field MRI/MRS for human application

Magnetic resonance imaging and spectroscopic techniques are widely used in humans both for clinical diagnostic applications and in basic research areas such as cognitive neuroimaging. In recent years, new human MR systems have become available operating at static magnetic fields of 7 T or higher (≥300 MHz proton frequency). Imaging human-sized objects at such high frequencies presents several challenges including non-uniform radiofrequency fields, enhanced susceptibility artifacts, and higher radiofrequency energy deposition in the tissue. On the other side of the scale are gains in signal-to-noise or contrast-to-noise ratio that allow finer structures to be visualized and smaller physiological effects to be detected. This review presents an overview of some of the latest methodological developments in human ultra-high field MRI/MRS as well as associated clinical and scientific applications. Emphasis is given to techniques that particularly benefit from the changing physical characteristics at high magnetic fields, including susceptibility-weighted imaging and phase-contrast techniques, imaging with X-nuclei, MR spectroscopy, CEST imaging, as well as functional MRI. In addition, more general methodological developments such as parallel transmission and motion correction will be discussed that are required to leverage the full potential of higher magnetic fields, and an overview of relevant physiological considerations of human high magnetic field exposure is provided.

Structure tensor informed fibre tractography at 3T

Structure tensor informed fibre tractography (STIFT) based on informing tractography for diffusion-weighted images at 3T and by utilising the structure tensor obtained from gradient-recalled echo (GRE) images at 7T is able to delineate fibres when seed voxels are placed close to the fibre boundaries. However, incorporating data from two different field strengths limits the applicability of STIFT. In this study, STIFT was implemented with both diffusion-weighted images and GRE images acquired at 3T. Instead of using the magnitude GRE data directly for STIFT as in the previous work, the utility of T₂ * maps and quantitative susceptibility maps derived from complex-valued GRE data to improve fibre delineation was explored. Single-seed tractography was performed and the results show that the optic radiation reconstructed with STIFT is more distinguishable from the inferior longitudinal fasciculus/inferior fronto-occipital fasciculus complex when compared to standard diffusion-weighted imaging tractography. We further investigated the quantitative effects of STIFT in a group of five healthy volunteers and evaluated its impact on measures of structural connectivity. The framework was extended to evaluate implementations of STIFT based on T₂ *-weighted and quantitative susceptibility-weighted images in a whole-brain connectivity study. In terms of connectivity, no systematic differences were found between STIFT and diffusion-weighted imaging tractography, suggesting that local improvements in tractography are not translated to the atlas-based structural connectivity analysis. Nevertheless, the reduction in the number of statistically significant connections in the STIFT connectivity matrix suggests that STIFT can potentially reduce the false-positive connections in fibre tractography.

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.

Keeping track of sound objects in space: The contribution of early-stage auditory areas

The influential dual-stream model of auditory processing stipulates that information pertaining to the meaning and to the position of a given sound object is processed in parallel along two distinct pathways, the ventral and dorsal auditory streams. Functional independence of the two processing pathways is well documented by conscious experience of patients with focal hemispheric lesions. On the other hand there is growing evidence that the meaning and the position of a sound are combined early in the processing pathway, possibly already at the level of early-stage auditory areas. Here, we investigated how early auditory areas integrate sound object meaning and space (simulated by interaural time differences) using a repetition suppression fMRI paradigm at 7 T. Subjects listen passively to environmental sounds presented in blocks of repetitions of the same sound object (same category) or different sounds objects (different categories), perceived either in the left or right space (no change within block) or shifted left-to-right or right-to-left halfway in the block (change within block). Environmental sounds activated bilaterally the superior temporal gyrus, middle temporal gyrus, inferior frontal gyrus, and right precentral cortex. Repetitions suppression effects were measured within bilateral early-stage auditory areas in the lateral portion of the Heschl's gyrus and posterior superior temporal plane. Left lateral early-stages areas showed significant effects for position and change, interactions Category x Initial Position and Category x Change in Position, while right lateral areas showed main effect of category and interaction Category x Change in Position. The combined evidence from our study and from previous studies speaks in favour of a position-linked representation of sound objects, which is independent from semantic encoding within the ventral stream and from spatial encoding within the dorsal stream. We argue for a third auditory stream, which has its origin in lateral belt areas and tracks sound objects across space.

The impact of B1+ correction on MP2RAGE cortical T1 and apparent cortical thickness at 7T

Determination of cortical thickness using MRI has often been criticized due to the presence of various error sources. Specifically, anatomical MRI relying on T₁ contrast may be unreliable due to spatially variable image contrast between gray matter (GM), white matter (WM) and cerebrospinal fluid (CSF). Especially at ultra‐high field (≥ 7T) MRI, transmit and receive B₁‐related image inhomogeneities can hamper correct classification of tissue types. In the current paper, we demonstrate that residual B1+ (transmit) inhomogeneities in the T₁‐weighted and quantitative T₁ images using the MP2RAGE sequence at 7T lead to biases in cortical thickness measurements. As expected, post‐hoc correction for the spatially varying B1+ profile reduced the apparent T₁ values across the cortex in regions with low B1+, and slightly increased apparent T₁ in regions with high B1+. As a result, improved contrast‐to‐noise ratio both at the GM‐CSF and GM‐WM boundaries can be observed leading to more accurate surface reconstructions and cortical thickness estimates. Overall, the changes in cortical thickness ranged between a 5% decrease to a 70% increase after B1+ correction, reducing the variance of cortical thickness values across the brain dramatically and increasing the comparability with normative data. More specifically, the cortical thickness estimates increased in regions characterized by a strong decrease of apparent T₁ after B1+ correction in regions with low B1+ due to improved detection of the pial surface_. The current results suggest that cortical thickness can be more accurately determined using MP2RAGE data at 7T if B1+ inhomogeneities are accounted for.