Mixed methodology in human brain research: integrating MRI and histology.

Postmortem magnetic resonance imaging (MRI) can provide a bridge between histological observations and the in vivo anatomy of the human brain. Approaches aimed at the co-registration of data derived from the two techniques are gaining interest. Optimal integration of the two research fields requires detailed knowledge of the tissue property requirements for individual research techniques, as well as a detailed understanding of the consequences of tissue fixation steps on the imaging quality outcomes for both MRI and histology. Here, we provide an overview of existing studies that bridge between state-of-the-art imaging modalities, and discuss the background knowledge incorporated into the design, execution and interpretation of postmortem studies. A subset of the discussed challenges transfer to animal studies as well. This insight can contribute to furthering our understanding of the normal and diseased human brain, and to facilitate discussions between researchers from the individual disciplines.

A unified 3D map of microscopic architecture and MRI of the human brain.

We present the first three-dimensional (3D) concordance maps of cyto- and fiber architecture of the human brain, combining histology, immunohistochemistry, and 7-T quantitative magnetic resonance imaging (MRI), in two individual specimens. These 3D maps each integrate data from approximately 800 microscopy sections per brain, showing neuronal and glial cell bodies, nerve fibers, and interneuronal populations, as well as ultrahigh-field quantitative MRI, all coaligned at the 200-µm scale to the stacked blockface images obtained during sectioning. These unprecedented 3D multimodal datasets are shared without any restrictions and provide a unique resource for the joint study of cell and fiber architecture of the brain, detailed anatomical atlasing, or modeling of the microscopic underpinnings of MRI contrasts.

A high-resolution multi-shell 3T diffusion magnetic resonance imaging dataset as part of the Amsterdam Ultra-high field adult lifespan database (AHEAD).

In order to further our understanding of brain function and the underlying networks, more advanced diffusion weighted magnetic resonance imaging (DWI MRI) data are essential. Here we present freely available high-resolution multi-shell multi-directional 3 Tesla (T) DWI MRI data as part of the 'Amsterdam Ultra-high field adult lifespan database' (AHEAD). The 3T DWI AHEAD dataset include 1.28mm isotropic whole brain DWI data of 49 healthy adult participants between 18 and 90 years old. The acquired data include DWIs at three non-zero b-values (48 directions, b-value 700 s/mm2; 56 directions, b-value 1000 s/mm2; 64 directions, b-value 1600 s/mm2) including a total of twelve volumes with a b-value of 0 s/mm2 (b0 volumes). In addition, eight b0 volumes with a reversed phase encoding direction were acquired to correct for distortions. To facilitate future use, the DWI data have been denoised, corrected for eddy currents, susceptibility-induced off-resonance field distortions, bias fields, and are skull stripped.

Charting human subcortical maturation across the adult lifespan with in vivo 7 T MRI.

The human subcortex comprises hundreds of unique structures. Subcortical functioning is crucial for behavior, and disrupted function is observed in common neurodegenerative diseases. Despite their importance, human subcortical structures continue to be difficult to study in vivo. Here we provide a detailed account of 17 prominent subcortical structures and ventricles, describing their approximate iron and myelin contents, morphometry, and their age-related changes across the normal adult lifespan. The results provide compelling insights into the heterogeneity and intricate age-related alterations of these structures. They also show that the locations of many structures shift across the lifespan, which is of direct relevance for the use of standard magnetic resonance imaging atlases. The results further our understanding of subcortical morphometry and neuroimaging properties, and of normal aging processes which ultimately can improve our understanding of neurodegeneration.

Manual delineation approaches for direct imaging of the subcortex.

The growing interest in the human subcortex is accompanied by an increasing number of parcellation procedures to identify deep brain structures in magnetic resonance imaging (MRI) contrasts. Manual procedures continue to form the gold standard for parcellating brain structures and is used for the validation of automated approaches. Performing manual parcellations is a tedious process which requires a systematic and reproducible approach. For this purpose, we created a series of protocols for the anatomical delineation of 21 individual subcortical structures. The intelligibility of the protocols was assessed by calculating Dice similarity coefficients for ten healthy volunteers. In addition, dilated Dice coefficients showed that manual parcellations created using these protocols can provide high-quality training data for automated algorithms. Here, we share the protocols, together with three example MRI datasets and the created manual delineations. The protocols can be applied to create high-quality training data for automated parcellation procedures, as well as for further validation of existing procedures and are shared without restrictions with the research community.

Structure-function similarities in deep brain stimulation targets cross-species.

Deep Brain Stimulation (DBS) is an effective neurosurgical treatment to alleviate motor symptoms of advanced Parkinson's disease. Due to its potential, DBS usage is rapidly expanding to target a large number of brain regions to treat a wide range of diseases and neuropsychiatric disorders. The identification and validation of new target regions heavily rely on the insights gained from rodent and primate models. Here we present a large-scale automatic meta-analysis in which the structure-function associations within and between species are compared for 21 DBS targets in humans. The results indicate that the structure-function association for the majority of the 21 included subcortical areas were conserved cross-species. A subset of structures showed overlapping functional association. This can potentially be attributed to shared brain networks and might explain why multiple brain areas are targeted for the same disease or neuropsychiatric disorder.

Imaging of the human subthalamic nucleus.

The human subthalamic nucleus (STN) is a small lens shaped iron rich nucleus, which has gained substantial interest as a target for deep brain stimulation surgery for a variety of movement disorders. The internal anatomy of the human STN has not been fully elucidated, and an intensive debate, discussing the level of overlap between putative limbic, associative, and motor zones within the STN is still ongoing. In this chapter, we have summarized anatomical information obtained using different neuroimaging modalities focusing on the anatomy of the STN. Additionally, we have highlighted a number of major challenges faced when using magnetic resonance imaging (MRI) approaches for the visualization of small iron rich deep brain structures such as the STN. In vivo MRI and postmortem microscopy efforts provide valuable complementary information on the internal structure of the STN, although the results are not always fully aligned. Finally, we provide an outlook on future efforts that could contribute to the development of an integrative research approach that will help with the reconciliation of seemingly divergent results across research approaches.

Measuring the iron content of dopaminergic neurons in substantia nigra with MRI relaxometry.

In Parkinson's disease, the depletion of iron-rich dopaminergic neurons in nigrosome 1 of the substantia nigra precedes motor symptoms by two decades. Methods capable of monitoring this neuronal depletion, at an early disease stage, are needed for early diagnosis and treatment monitoring. Magnetic resonance imaging (MRI) is particularly suitable for this task due to its sensitivity to tissue microstructure and in particular, to iron. However, the exact mechanisms of MRI contrast in the substantia nigra are not well understood, hindering the development of powerful biomarkers. In the present report, we illuminate the contrast mechanisms in gradient and spin echo MR images in human nigrosome 1 by combining quantitative 3D iron histology and biophysical modeling with quantitative MRI on post mortem human brain tissue. We show that the dominant contribution to the effective transverse relaxation rate (R2*) in nigrosome 1 originates from iron accumulated in the neuromelanin of dopaminergic neurons. This contribution is appropriately described by a static dephasing approximation of the MRI signal. We demonstrate that the R2* contribution from dopaminergic neurons reflects the product of cell density and cellular iron concentration. These results demonstrate that the in vivo monitoring of neuronal density and iron in nigrosome 1 may be feasible with MRI and provide directions for the development of biomarkers for an early detection of dopaminergic neuron depletion in Parkinson's disease.

Quantity and quality: Normative open-access neuroimaging databases.

The focus of this article is to compare twenty normative and open-access neuroimaging databases based on quantitative measures of image quality, namely, signal-to-noise (SNR) and contrast-to-noise ratios (CNR). We further the analysis through discussing to what extent these databases can be used for the visualization of deeper regions of the brain, such as the subcortex, as well as provide an overview of the types of inferences that can be drawn. A quantitative comparison of contrasts including T1-weighted (T1w) and T2-weighted (T2w) images are summarized, providing evidence for the benefit of ultra-high field MRI. Our analysis suggests a decline in SNR in the caudate nuclei with increasing age, in T1w, T2w, qT1 and qT2* contrasts, potentially indicative of complex structural age-dependent changes. A similar decline was found in the corpus callosum of the T1w, qT1 and qT2* contrasts, though this relationship is not as extensive as within the caudate nuclei. These declines were accompanied by a declining CNR over age in all image contrasts. A positive correlation was found between scan time and the estimated SNR as well as a negative correlation between scan time and spatial resolution. Image quality as well as the number and types of contrasts acquired by these databases are important factors to take into account when selecting structural data for reuse. This article highlights the opportunities and pitfalls associated with sampling existing databases, and provides a quantitative backing for their usage.

Multi-contrast anatomical subcortical structures parcellation.

The human subcortex is comprised of more than 450 individual nuclei which lie deep in the brain. Due to their small size and close proximity, up until now only 7% have been depicted in standard MRI atlases. Thus, the human subcortex can largely be considered as terra incognita. Here, we present a new open-source parcellation algorithm to automatically map the subcortex. The new algorithm has been tested on 17 prominent subcortical structures based on a large quantitative MRI dataset at 7 Tesla. It has been carefully validated against expert human raters and previous methods, and can easily be extended to other subcortical structures and applied to any quantitative MRI dataset. In sum, we hope this novel parcellation algorithm will facilitate functional and structural neuroimaging research into small subcortical nuclei and help to chart terra incognita.

3 versus 7 Tesla magnetic resonance imaging for parcellations of subcortical brain structures in clinical settings.

7 Tesla (7T) magnetic resonance imaging holds great promise for improved visualization of the human brain for clinical purposes. To assess whether 7T is superior regarding localization procedures of small brain structures, we compared manual parcellations of the red nucleus, subthalamic nucleus, substantia nigra, globus pallidus interna and externa. These parcellations were created on a commonly used clinical anisotropic clinical 3T with an optimized isotropic (o)3T and standard 7T scan. The clinical 3T MRI scans did not allow delineation of an anatomically plausible structure due to its limited spatial resolution. o3T and 7T parcellations were directly compared. We found that 7T outperformed the o3T MRI as reflected by higher Dice scores, which were used as a measurement of interrater agreement for manual parcellations on quantitative susceptibility maps. This increase in agreement was associated with higher contrast to noise ratios for smaller structures, but not for the larger globus pallidus segments. Additionally, control-analyses were performed to account for potential biases in manual parcellations by assessing semi-automatic parcellations. These results showed a higher consistency for structure volumes for 7T compared to optimized 3T which illustrates the importance of the use of isotropic voxels for 3D visualization of the surgical target area. Together these results indicate that 7T outperforms c3T as well as o3T given the constraints of a clinical setting.

7 Tesla MRI Followed by Histological 3D Reconstructions in Whole-Brain Specimens.

Post mortem magnetic resonance imaging (MRI) studies on the human brain are of great interest for the validation of in vivo MRI. It facilitates a link between functional and anatomical information available from MRI in vivo and neuroanatomical knowledge available from histology/immunocytochemistry. However, linking in vivo and post mortem MRI to microscopy techniques poses substantial challenges. Fixation artifacts and tissue deformation of extracted brains, as well as co registration of 2D histology to 3D MRI volumes complicate direct comparison between modalities. Moreover, post mortem brain tissue does not have the same physical properties as in vivo tissue, and therefore MRI approaches need to be adjusted accordingly. Here, we present a pipeline in which whole-brain human post mortem in situ MRI is combined with subsequent tissue processing of the whole human brain, providing a 3-dimensional reconstruction via blockface imaging. To this end, we adapted tissue processing procedures to allow both post mortem MRI and subsequent histological and immunocytochemical processing. For MRI, tissue was packed in a susceptibility matched solution, tailored to fit the dimensions of the MRI coil. Additionally, MRI sequence parameters were adjusted to accommodate T1 and T2∗ shortening, and scan time was extended, thereby benefiting the signal-to-noise-ratio that can be achieved using extensive averaging without motion artifacts. After MRI, the brain was extracted from the skull and subsequently cut while performing optimized blockface imaging, thereby allowing three-dimensional reconstructions. Tissues were processed for Nissl and silver staining, and co-registered with the blockface images. The combination of these techniques allows direct comparisons across modalities.

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.

Functional segregation and integration within the human subthalamic nucleus from a micro- and meso-level perspective.

The subthalamic nucleus (STN) is a core basal ganglia structure involved in the control of motor, cognitive, motivational and affective functions. The (challenged) tripartite subdivision hypothesis places these functions into distinct sensorimotor, cognitive/associative, and limbic subregions based on the topography of cortical projections. To a large extent, this hypothesis is used to motivate the choice of target coordinates for implantation of deep brain stimulation electrodes for treatment of neurological and psychiatric disorders. Yet, the parallel organization of basal ganglia circuits has been known to allow considerable cross-talk, which might contribute to the occurrence of neuropsychiatric side effects when stimulating the dorsolateral, putative sensorimotor, part of the STN for treatment of Parkinson's disease. Any functional segregation within the STN is expected to be reflected both at micro-level microscopy and meso-level neural population activity. As such, we review the current empirical evidence from anterograde tracing and immunocytochemistry studies and from local field potential recordings for delineating the STN into distinct subregions. The spatial distribution of immunoreactivity presents as a combination of gradients, and although neural activity in distinct frequency bands appears spatially clustered, there is substantial overlap in peak locations. We argue that regional specialization without sharply defined borders is likely most representative of the STN's functional organization.

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 R1, R2* and Quantitative Susceptibility Maps and T1-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.

The Amsterdam Ultra-high field adult lifespan database (AHEAD): A freely available multimodal 7 Tesla submillimeter magnetic resonance imaging database.

Normative databases allow testing of novel hypotheses without the costly collection of magnetic resonance imaging (MRI) data. Here we present the Amsterdam Ultra-high field adult lifespan database (AHEAD). The AHEAD consists of 105 7 Tesla (T) whole-brain structural MRI scans tailored specifically to imaging of the human subcortex, including both male and female participants and covering the entire adult life span (18-80 yrs). We used these data to create probability maps for the subthalamic nucleus, substantia nigra, internal and external segment of the globus pallidus, and the red nucleus. Data was acquired at a submillimeter resolution using a multi-echo (ME) extension of the second gradient-echo image of the MP2RAGE sequence (MP2RAGEME) sequence, resulting in complete anatomical alignment of quantitative, R1-maps, R2*-maps, T1-maps, T1-weighted images, T2*-maps, and quantitative susceptibility mapping (QSM). Quantitative MRI maps, and derived probability maps of basal ganglia structures are freely available for further analyses.

Denoising High-Field Multi-Dimensional MRI With Local Complex PCA.

Modern high field and ultra high field magnetic resonance imaging (MRI) experiments routinely collect multi-dimensional data with high spatial resolution, whether multi-parametric structural, diffusion or functional MRI. While diffusion and functional imaging have benefited from recent advances in multi-dimensional signal analysis and denoising, structural MRI has remained untouched. In this work, we propose a denoising technique for multi-parametric quantitative MRI, combining a highly popular denoising method from diffusion imaging, over-complete local PCA, with a reconstruction of the complex-valued MR signal in order to define stable estimates of the noise in the decomposition. With this approach, we show signal to noise ratio (SNR) improvements in high resolution MRI without compromising the spatial accuracy or generating spurious perceptual boundaries.

The functional microscopic neuroanatomy of the human subthalamic nucleus.

The subthalamic nucleus (STN) is successfully used as a surgical target for deep brain stimulation in the treatment of movement disorders. Interestingly, the internal structure of the STN is still incompletely understood. The objective of the present study was to investigate three-dimensional (3D) immunoreactivity patterns for 12 individual protein markers for GABA-ergic, serotonergic, dopaminergic as well as glutamatergic signaling. We analyzed the immunoreactivity using optical densities and created a 3D reconstruction of seven postmortem human STNs. Quantitative modeling of the reconstructed 3D immunoreactivity patterns revealed that the applied protein markers show a gradient distribution in the STN. These gradients were predominantly organized along the ventromedial to dorsolateral axis of the STN. The results are of particular interest in view of the theoretical underpinning for surgical targeting, which is based on a tripartite distribution of cognitive, limbic and motor function in the STN.

Size and shape matter: The impact of voxel geometry on the identification of small nuclei.

How, and to what extent do size and shape of a voxel measured with magnetic resonance imaging (MRI) affect the ability to visualize small brain nuclei? Despite general consensus that voxel geometry affects volumetric properties of regions of interest, particularly those of small brain nuclei, no quantitative data on the influence of voxel size and shape on labeling accuracy is available. Using simulations, we investigated the selective influence of voxel geometry by reconstructing simulated ellipsoid structures with voxels varying in shape and size. For each reconstructed ellipsoid, we calculated differences in volume and similarity between the labeled volume and the predefined dimensions of the ellipsoid. Probability functions were derived from one or two individual raters and a simulated ground truth for reference. As expected, larger voxels (i.e., coarser resolution) and increasing anisotropy results in increased deviations of both volume and shape measures, which is of particular relevance for small brain structures. Our findings clearly illustrate the anatomical inaccuracies introduced by the application of large and/or anisotropic voxels. To ensure deviations occur within the acceptable range (Dice coefficient scores; DCS > 0.75, corresponding to < 57% volume deviation), the volume of isotropic voxels should not exceed 5% of the total volume of the region of interest. When high accuracy is required (DCS > 0.90, corresponding to a < 19% volume deviation), the volumes of isotropic voxels should not exceed 0.08%, of the total volume. Finally, when large anisotropic factors (>3) are used, and the ellipsoid is orthogonal to the slice axes, having its long axis in the imaging plane, the voxel volume should not exceed 0.005% of the total volume. This allows sufficient compensation of anisotropy effects, in order to reach accuracy in the acceptable range (DCS > 0.75, corresponding to >57% volume deviation).