The fine spatial scales of the structures in the human brain represent an enormous challenge to the successful integration of information from different images for both within- and between-subject ...analysis. While many algorithms to register image pairs from the same subject exist, visual inspection shows that their accuracy and robustness to be suspect, particularly when there are strong intensity gradients and/or only part of the brain is imaged. This paper introduces a new algorithm called Boundary-Based Registration, or BBR. The novelty of BBR is that it treats the two images very differently. The reference image must be of sufficient resolution and quality to extract surfaces that separate tissue types. The input image is then aligned to the reference by maximizing the intensity gradient across tissue boundaries. Several lower quality images can be aligned through their alignment with the reference. Visual inspection and fMRI results show that BBR is more accurate than correlation ratio or normalized mutual information and is considerably more robust to even strong intensity inhomogeneities. BBR also excels at aligning partial-brain images to whole-brain images, a domain in which existing registration algorithms frequently fail. Even in the limit of registering a single slice, we show the BBR results to be robust and accurate.
The false positive rates (FPR) for surface-based group analysis of cortical thickness, surface area, and volume were evaluated for parametric and non-parametric clusterwise correction for multiple ...comparisons for a range of smoothing levels and cluster-forming thresholds (CFT) using real data under group assignments that should not yield significant results. For whole cortical surface analysis, thickness showed modest inflation in parametric FPRs above the nominal level (10% versus 5%). Surface area and volume FPRs were much higher (20–30%). In the analysis of interhemispheric thickness asymmetries, FPRs were well controlled by parametric correction, but FPRs for surface area and volume asymmetries were still inflated. In all cases, non-parametric permutation adequately controlled the FPRs. It was found that inflated parametric FPRs were caused by violations in the parametric assumptions, namely a heavier-than-Gaussian spatial correlation. The non-Gaussian spatial correlation originates from anatomical features unique to individuals (e.g., a patch of cortex slightly thicker or thinner than average) and is not a by-product of scanning or processing. Thickness performed better than surface area and volume because thickness does not require a Jacobian correction.
•False positive rates (FPR) for surface-based group analysis using real MRI data.•Parametric (Gaussian-based Monte Carlo) and non-parametric (permutation).•Thickness analysis showed slightly elevated parametric FPRs.•Surface area and volume analysis showed significantly elevated parametric FPRs.•Permutation showed good control of FPRs in all cases.
The human thalamus is a brain structure that comprises numerous, highly specific nuclei. Since these nuclei are known to have different functions and to be connected to different areas of the ...cerebral cortex, it is of great interest for the neuroimaging community to study their volume, shape and connectivity in vivo with MRI. In this study, we present a probabilistic atlas of the thalamic nuclei built using ex vivo brain MRI scans and histological data, as well as the application of the atlas to in vivo MRI segmentation. The atlas was built using manual delineation of 26 thalamic nuclei on the serial histology of 12 whole thalami from six autopsy samples, combined with manual segmentations of the whole thalamus and surrounding structures (caudate, putamen, hippocampus, etc.) made on in vivo brain MR data from 39 subjects. The 3D structure of the histological data and corresponding manual segmentations was recovered using the ex vivo MRI as reference frame, and stacks of blockface photographs acquired during the sectioning as intermediate target. The atlas, which was encoded as an adaptive tetrahedral mesh, shows a good agreement with previous histological studies of the thalamus in terms of volumes of representative nuclei. When applied to segmentation of in vivo scans using Bayesian inference, the atlas shows excellent test-retest reliability, robustness to changes in input MRI contrast, and ability to detect differential thalamic effects in subjects with Alzheimer's disease. The probabilistic atlas and companion segmentation tool are publicly available as part of the neuroimaging package FreeSurfer.
Display omitted
•A probabilistic atlas of the human thalamus (26 nuclei) derived from histology.•3D histology reconstruction assisted by ex vivo MRI and blockface photographs.•Companion Bayesian method segments the nuclei from in vivo scans of any MRI contrast.•Method has excellent test-retest reliability and detects differential effects in AD.•The method is publicly available as part of FreeSurfer.
With sufficient image encoding, high-resolution fMRI studies are limited by the biological point-spread of the hemodynamic signal. The extent of this spread is determined by the local vascular ...distribution and by the spatial specificity of blood flow regulation, as well as by measurement parameters that (i) alter the relative sensitivity of the acquisition to activation-induced hemodynamic changes and (ii) determine the image contrast as a function of vessel size. In particular, large draining vessels on the cortical surface are a major contributor to both the BOLD signal change and to the spatial bias of the BOLD activation away from the site of neuronal activity. In this work, we introduce a laminar surface-based analysis method and study the relationship between spatial localization and activation strength as a function of laminar depth by acquiring 1mm isotropic, single-shot EPI at 7T and sampling the BOLD signal exclusively from the superficial, middle, or deep cortical laminae. We show that highly-accelerated EPI can limit image distortions to the point where a boundary-based registration algorithm accurately aligns the EPI data to the surface reconstruction. The spatial spread of the BOLD response tangential to the cortical surface was analyzed as a function of cortical depth using our surface-based analysis. Although sampling near the pial surface provided the highest signal strength, it also introduced the most spatial error. Thus, avoiding surface laminae improved spatial localization by about 40% at a cost of 36% in z-statistic, implying that optimal spatial resolution in functional imaging of the cortex can be achieved using anatomically-informed spatial sampling to avoid large pial vessels.
We introduce a strategy for learning image registration without acquired imaging data, producing powerful networks agnostic to contrast introduced by magnetic resonance imaging (MRI). While classical ...registration methods accurately estimate the spatial correspondence between images, they solve an optimization problem for every new image pair. Learning-based techniques are fast at test time but limited to registering images with contrasts and geometric content similar to those seen during training. We propose to remove this dependency on training data by leveraging a generative strategy for diverse synthetic label maps and images that exposes networks to a wide range of variability, forcing them to learn more invariant features. This approach results in powerful networks that accurately generalize to a broad array of MRI contrasts. We present extensive experiments with a focus on 3D neuroimaging, showing that this strategy enables robust and accurate registration of arbitrary MRI contrasts even if the target contrast is not seen by the networks during training. We demonstrate registration accuracy surpassing the state of the art both within and across contrasts, using a single model. Critically, training on arbitrary shapes synthesized from noise distributions results in competitive performance, removing the dependency on acquired data of any kind. Additionally, since anatomical label maps are often available for the anatomy of interest, we show that synthesizing images from these dramatically boosts performance, while still avoiding the need for real intensity images. Our code is available at doic https://w3id.org/synthmorph .
Gamma-aminobutyric acid (GABA) is the main inhibitory neurotransmitter in the human brain and plays a key role in several brain functions and neuropsychiatric disorders such as anxiety, epilepsy, and ...depression. For decades, several in vivo and ex vivo techniques have been used to highlight the mechanisms of the GABA system, however, no studies have currently combined the techniques to create a high-resolution multimodal view of the GABA system.
Here, we present a quantitative high-resolution in vivo atlas of the human brain benzodiazepine receptor sites (BZR) located on postsynaptic ionotropic GABAA receptors (GABAARs), generated on the basis of in vivo 11Cflumazenil Positron Emission Tomography (PET) data. Next, based on ex vivo autoradiography data, we transform the PET-generated atlas from binding values into BZR protein density. Finally, we examine the brain regional association between BZR protein density and ex vivo mRNA expression for the 19 subunits in the GABAAR, including an estimation of the minimally required expression of mRNA levels for each subunit to translate into BZR protein.
This represents the first publicly available quantitative high-resolution in vivo atlas of the spatial distribution of BZR densities in the healthy human brain. The atlas provides a unique neuroscientific tool as well as novel insights into the association between mRNA expression for individual subunits in the GABAAR and the BZR density at each location in the brain.
Despite advances in data augmentation and transfer learning, convolutional neural networks (CNNs) difficultly generalise to unseen domains. When segmenting brain scans, CNNs are highly sensitive to ...changes in resolution and contrast: even within the same MRI modality, performance can decrease across datasets. Here we introduce SynthSeg, the first segmentation CNN robust against changes in contrast and resolution. SynthSeg is trained with synthetic data sampled from a generative model conditioned on segmentations. Crucially, we adopt a domain randomisation strategy where we fully randomise the contrast and resolution of the synthetic training data. Consequently, SynthSeg can segment real scans from a wide range of target domains without retraining or fine-tuning, which enables straightforward analysis of huge amounts of heterogeneous clinical data. Because SynthSeg only requires segmentations to be trained (no images), it can learn from labels obtained by automated methods on diverse populations (e.g., ageing and diseased), thus achieving robustness to a wide range of morphological variability. We demonstrate SynthSeg on 5,000 scans of six modalities (including CT) and ten resolutions, where it exhibits unparallelled generalisation compared with supervised CNNs, state-of-the-art domain adaptation, and Bayesian segmentation. Finally, we demonstrate the generalisability of SynthSeg by applying it to cardiac MRI and CT scans.
Display omitted
•A CNN to segment brain MRI scans of any contrast and resolution without retraining.•Domain-independence is achieved by training on randomised unrealistic synthetic data.•SynthSeg sustains almost the accuracy of supervised CNNs across all tested domains.•It also outperforms state-of-the-art domain adaptation, without being retrained.•The model is implemented in FreeSurfer for easy distribution and deployment.
The serotonin (5-hydroxytryptamine, 5-HT) system modulates many important brain functions and is critically involved in many neuropsychiatric disorders. Here, we present a high-resolution, ...multidimensional, in vivo atlas of four of the human brain's 5-HT receptors (5-HT
, 5-HT
, 5-HT
, and 5-HT
) and the 5-HT transporter (5-HTT). The atlas is created from molecular and structural high-resolution neuroimaging data consisting of positron emission tomography (PET) and magnetic resonance imaging (MRI) scans acquired in a total of 210 healthy individuals. Comparison of the regional PET binding measures with postmortem human brain autoradiography outcomes showed a high correlation for the five 5-HT targets and this enabled us to transform the atlas to represent protein densities (in picomoles per milliliter). We also assessed the regional association between protein concentration and mRNA expression in the human brain by comparing the 5-HT density across the atlas with data from the Allen Human Brain atlas and identified receptor- and transporter-specific associations that show the regional relation between the two measures. Together, these data provide unparalleled insight into the serotonin system of the human brain.
We present a high-resolution positron emission tomography (PET)- and magnetic resonance imaging-based human brain atlas of important serotonin receptors and the transporter. The regional PET-derived binding measures correlate strongly with the corresponding autoradiography protein levels. The strong correlation enables the transformation of the PET-derived human brain atlas into a protein density map of the serotonin (5-hydroxytryptamine, 5-HT) system. Next, we compared the regional receptor/transporter protein densities with mRNA levels and uncovered unique associations between protein expression and density at high detail. This new in vivo neuroimaging atlas of the 5-HT system not only provides insight in the human brain's regional protein synthesis, transport, and density, but also represents a valuable source of information for the neuroscience community as a comparative instrument to assess brain disorders.
A tool was developed to automatically segment several subcortical limbic structures (nucleus accumbens, basal forebrain, septal nuclei, hypothalamus without mammillary bodies, the mammillary bodies, ...and fornix) using only a T1-weighted MRI as input. This tool fills an unmet need as there are few, if any, publicly available tools to segment these clinically relevant structures. A U-Net with spatial, intensity, contrast, and noise augmentation was trained using 39 manually labeled MRI data sets. In general, the Dice scores, true positive rates, false discovery rates, and manual-automatic volume correlation were very good relative to comparable tools for other structures. A diverse data set of 698 subjects were segmented using the tool; evaluation of the resulting labelings showed that the tool failed in less than 1% of cases. Test-retest reliability of the tool was excellent. The automatically segmented volume of all structures except mammillary bodies showed effectiveness at detecting either clinical AD effects, age effects, or both. This tool will be publicly released with FreeSurfer (surfer.nmr.mgh.harvard.edu/fswiki/ScLimbic). Together with the other cortical and subcortical limbic segmentations, this tool will allow FreeSurfer to provide a comprehensive view of the limbic system in an automated way.
Exploratory (i.e., voxelwise) spatial methods are commonly used in neuroimaging to identify areas that show an effect when a region-of-interest (ROI) analysis cannot be performed because no strong a ...priori anatomical hypothesis exists. However, noise at a single voxel is much higher than noise in a ROI making noise management critical to successful exploratory analysis. This work explores how preprocessing choices affect the bias and variability of voxelwise kinetic modeling analysis of brain positron emission tomography (PET) data. These choices include the use of volume- or cortical surface-based smoothing, level of smoothing, use of voxelwise partial volume correction (PVC), and PVC masking threshold. PVC was implemented using the Muller-Gartner method with the masking out of voxels with low gray matter (GM) partial volume fraction. Dynamic PET scans of an antagonist serotonin-4 receptor radioligand (11CSB207145) were collected on sixteen healthy subjects using a Siemens HRRT PET scanner. Kinetic modeling was used to compute maps of non-displaceable binding potential (BPND) after preprocessing. The results showed a complicated interaction between smoothing, PVC, and masking on BPND estimates. Volume-based smoothing resulted in large bias and intersubject variance because it smears signal across tissue types. In some cases, PVC with volume smoothing paradoxically caused the estimated BPND to be less than when no PVC was used at all. When applied in the absence of PVC, cortical surface-based smoothing resulted in dramatically less bias and the least variance of the methods tested for smoothing levels 5mm and higher. When used in combination with PVC, surface-based smoothing minimized the bias without significantly increasing the variance. Surface-based smoothing resulted in 2–4 times less intersubject variance than when volume smoothing was used. This translates into more than 4 times fewer subjects needed in a group analysis to achieve similarly powered statistical tests. Surface-based smoothing has less bias and variance because it respects cortical geometry by smoothing the PET data only along the cortical ribbon and so does not contaminate the GM signal with that of white matter and cerebrospinal fluid. The use of surface-based analysis in PET should result in substantial improvements in the reliability and detectability of effects in exploratory PET analysis, with or without PVC.
•This is a study of the effects of preprocessing on exploratory PET analysis.•Volumetric 3D smoothing of PET data exacerbates the partial volume effect.•Volumetric smoothing and partial volume correction combine to increase variability.•Cortical surface-based smoothing does not exacerbate the partial volume effect.•Surface-based kinetic analysis reduces intersubject variance by a factor of 4.