Robust machine learning segmentation for large-scale analysis of heterogeneous clinical brain MRI datasets.

Every year, millions of brain MRI scans are acquired in hospitals, which is a figure considerably larger than the size of any research dataset. Therefore, the ability to analyze such scans could transform neuroimaging research. Yet, their potential remains untapped since no automated algorithm is robust enough to cope with the high variability in clinical acquisitions (MR contrasts, resolutions, orientations, artifacts, and subject populations). Here, we present SynthSeg+, an AI segmentation suite that enables robust analysis of heterogeneous clinical datasets. In addition to whole-brain segmentation, SynthSeg+ also performs cortical parcellation, intracranial volume estimation, and automated detection of faulty segmentations (mainly caused by scans of very low quality). We demonstrate SynthSeg+ in seven experiments, including an aging study on 14,000 scans, where it accurately replicates atrophy patterns observed on data of much higher quality. SynthSeg+ is publicly released as a ready-to-use tool to unlock the potential of quantitative morphometry.

Diffusion-Weighted Imaging Fluid-Attenuated Inversion Recovery Mismatch on Portable, Low-Field Magnetic Resonance Imaging Among Acute Stroke Patients.

OBJECTIVE: For stroke patients with unknown time of onset, mismatch between diffusion-weighted imaging (DWI) and fluid-attenuated inversion recovery (FLAIR) magnetic resonance imaging (MRI) can guide thrombolytic intervention. However, access to MRI for hyperacute stroke is limited. Here, we sought to evaluate whether a portable, low-field (LF)-MRI scanner can identify DWI-FLAIR mismatch in acute ischemic stroke. METHODS: Eligible patients with a diagnosis of acute ischemic stroke underwent LF-MRI acquisition on a 0.064-T scanner within 24 h of last known well. Qualitative and quantitative metrics were evaluated. Two trained assessors determined the visibility of stroke lesions on LF-FLAIR. An image coregistration pipeline was developed, and the LF-FLAIR signal intensity ratio (SIR) was derived. RESULTS: The study included 71 patients aged 71 ± 14 years and a National Institutes of Health Stroke Scale of 6 (interquartile range 3-14). The interobserver agreement for identifying visible FLAIR hyperintensities was high (κ = 0.85, 95% CI 0.70-0.99). Visual DWI-FLAIR mismatch had a 60% sensitivity and 82% specificity for stroke patients <4.5 h, with a negative predictive value of 93%. LF-FLAIR SIR had a mean value of 1.18 ± 0.18 <4.5 h, 1.24 ± 0.39 4.5-6 h, and 1.40 ± 0.23 >6 h of stroke onset. The optimal cut-point for LF-FLAIR SIR was 1.15, with 85% sensitivity and 70% specificity. A cut-point of 6.6 h was established for a FLAIR SIR <1.15, with an 89% sensitivity and 62% specificity. INTERPRETATION: A 0.064-T portable LF-MRI can identify DWI-FLAIR mismatch among patients with acute ischemic stroke. Future research is needed to prospectively validate thresholds and evaluate a role of LF-MRI in guiding thrombolysis among stroke patients with uncertain time of onset. ANN NEUROL 2024.

Three-dimensional histology reveals dissociable human hippocampal long-axis gradients of Alzheimer's pathology.

INTRODUCTION: Three-dimensional (3D) histology analyses are essential to overcome sampling variability and understand pathological differences beyond the dissection axis. We present Path2MR, the first pipeline allowing 3D reconstruction of sparse human histology without a magnetic resonance imaging (MRI) reference. We implemented Path2MR with post-mortem hippocampal sections to explore pathology gradients in Alzheimer's disease. METHODS: Blockface photographs of brain hemisphere slices are used for 3D reconstruction, from which an MRI-like image is generated using machine learning. Histology sections are aligned to the reconstructed hemisphere and subsequently to an atlas in standard space. RESULTS: Path2MR successfully registered histological sections to their anatomic position along the hippocampal longitudinal axis. Combined with histopathology quantification, we found an expected peak of tau pathology at the anterior end of the hippocampus, whereas amyloid-beta (Aß) displayed a quadratic anterior-posterior distribution. CONCLUSION: Path2MR, which enables 3D histology using any brain bank data set, revealed significant differences along the hippocampus between tau and Aß. HIGHLIGHTS: Path2MR enables three-dimensional (3D) brain reconstruction from blockface dissection photographs. This pipeline does not require dense specimen sampling or a subject-specific magnetic resonance (MR) image. Anatomically consistent mapping of hippocampal sections was obtained with Path2MR. Our analyses revealed an anterior-posterior gradient of hippocampal tau pathology. In contrast, the peak of amyloid-beta (Aß) deposition was closer to the hippocampal body.

Detection of Intracerebral Hemorrhage Using Low-Field, Portable Magnetic Resonance Imaging in Patients With Stroke.

BACKGROUND: Neuroimaging is essential for detecting spontaneous, nontraumatic intracerebral hemorrhage (ICH). Recent data suggest ICH can be characterized using low-field magnetic resonance imaging (MRI). Our primary objective was to investigate the sensitivity and specificity of ICH on a 0.064T portable MRI (pMRI) scanner using a methodology that provided clinical information to inform rater interpretations. As a secondary aim, we investigated whether the incorporation of a deep learning (DL) reconstruction algorithm affected ICH detection. METHODS: The pMRI device was deployed at Yale New Haven Hospital to examine patients presenting with stroke symptoms from October 26, 2020 to February 21, 2022. Three raters independently evaluated pMRI examinations. Raters were provided the images alongside the patient's clinical information to simulate real-world context of use. Ground truth was the closest conventional computed tomography or 1.5/3T MRI. Sensitivity and specificity results were grouped by DL and non-DL software to investigate the effects of software advances. RESULTS: A total of 189 exams (38 ICH, 89 acute ischemic stroke, 8 subarachnoid hemorrhage, 3 primary intraventricular hemorrhage, 51 no intracranial abnormality) were evaluated. Exams were correctly classified as positive or negative for ICH in 185 of 189 cases (97.9% overall accuracy). ICH was correctly detected in 35 of 38 cases (92.1% sensitivity). Ischemic stroke and no intracranial abnormality cases were correctly identified as blood-negative in 139 of 140 cases (99.3% specificity). Non-DL scans had a sensitivity and specificity for ICH of 77.8% and 97.1%, respectively. DL scans had a sensitivity and specificity for ICH of 96.6% and 99.3%, respectively. CONCLUSIONS: These results demonstrate improvements in ICH detection accuracy on pMRI that may be attributed to the integration of clinical information in rater review and the incorporation of a DL-based algorithm. The use of pMRI holds promise in providing diagnostic neuroimaging for patients with ICH.

Accurate Bayesian segmentation of thalamic nuclei using diffusion MRI and an improved histological atlas.

The human thalamus is a highly connected brain structure, which is key for the control of numerous functions and is involved in several neurological disorders. Recently, neuroimaging studies have increasingly focused on the volume and connectivity of the specific nuclei comprising this structure, rather than looking at the thalamus as a whole. However, accurate identification of cytoarchitectonically designed histological nuclei on standard in vivo structural MRI is hampered by the lack of image contrast that can be used to distinguish nuclei from each other and from surrounding white matter tracts. While diffusion MRI may offer such contrast, it has lower resolution and lacks some boundaries visible in structural imaging. In this work, we present a Bayesian segmentation algorithm for the thalamus. This algorithm combines prior information from a probabilistic atlas with likelihood models for both structural and diffusion MRI, allowing segmentation of 25 thalamic labels per hemisphere informed by both modalities. We present an improved probabilistic atlas, incorporating thalamic nuclei identified from histology and 45 white matter tracts surrounding the thalamus identified in ultra-high gradient strength diffusion imaging. We present a family of likelihood models for diffusion tensor imaging, ensuring compatibility with the vast majority of neuroimaging datasets that include diffusion MRI data. The use of these diffusion likelihood models greatly improves identification of nuclear groups versus segmentation based solely on structural MRI. Dice comparison of 5 manually identifiable groups of nuclei to ground truth segmentations show improvements of up to 10 percentage points. Additionally, our chosen model shows a high degree of reliability, with median test-retest Dice scores above 0.85 for four out of five nuclei groups, whilst also offering improved detection of differential thalamic involvement in Alzheimer's disease (AUROC 81.98%). The probabilistic atlas and segmentation tool will be made publicly available as part of the neuroimaging package FreeSurfer (https://freesurfer.net/fswiki/ThalamicNucleiDTI).

Quantitative Brain Morphometry of Portable Low-Field-Strength MRI Using Super-Resolution Machine Learning.

Background Portable, low-field-strength (0.064-T) MRI has the potential to transform neuroimaging but is limited by low spatial resolution and low signal-to-noise ratio. Purpose To implement a machine learning super-resolution algorithm that synthesizes higher spatial resolution images (1-mm isotropic) from lower resolution T1-weighted and T2-weighted portable brain MRI scans, making them amenable to automated quantitative morphometry. Materials and Methods An external high-field-strength MRI data set (1-mm isotropic scans from the Open Access Series of Imaging Studies data set) and segmentations for 39 regions of interest (ROIs) in the brain were used to train a super-resolution convolutional neural network (CNN). Secondary analysis of an internal test set of 24 paired low- and high-field-strength clinical MRI scans in participants with neurologic symptoms was performed. These were part of a prospective observational study (August 2020 to December 2021) at Massachusetts General Hospital (exclusion criteria: inability to lay flat, body habitus preventing low-field-strength MRI, presence of MRI contraindications). Three well-established automated segmentation tools were applied to three sets of scans: high-field-strength (1.5-3 T, reference standard), low-field-strength (0.064 T), and synthetic high-field-strength images generated from the low-field-strength data with the CNN. Statistical significance of correlations was assessed with Student t tests. Correlation coefficients were compared with Steiger Z tests. Results Eleven participants (mean age, 50 years ± 14; seven men) had full cerebrum coverage in the images without motion artifacts or large stroke lesion with distortion from mass effect. Direct segmentation of low-field-strength MRI yielded nonsignificant correlations with volumetric measurements from high field strength for most ROIs (P > .05). Correlations largely improved when segmenting the synthetic images: P values were less than .05 for all ROIs (eg, for the hippocampus [r = 0.85; P < .001], thalamus [r = 0.84; P = .001], and whole cerebrum [r = 0.92; P < .001]). Deviations from the model (z score maps) visually correlated with pathologic abnormalities. Conclusion This work demonstrated proof-of-principle augmentation of portable MRI with a machine learning super-resolution algorithm, which yielded highly correlated brain morphometric measurements to real higher resolution images. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Ertl-Wagner amd Wagner in this issue. An earlier incorrect version appeared online. This article was corrected on February 1, 2023.

Hippocampal sclerosis of aging at post-mortem is evident on MRI more than a decade prior.

INTRODUCTION: Hippocampal sclerosis of aging (HS) is an important component of combined dementia neuropathology. However, the temporal evolution of its histologically-defined features is unknown. We investigated pre-mortem longitudinal hippocampal atrophy associated with HS, as well as with other dementia-associated pathologies. METHODS: We analyzed hippocampal volumes from magnetic resonance imaging (MRI) segmentations in 64 dementia patients with longitudinal MRI follow-up and post-mortem neuropathological evaluation, including HS assessment in the hippocampal head and body. RESULTS: Significant HS-associated hippocampal volume changes were observed throughout the evaluated timespan, up to 11.75 years before death. These changes were independent of age and Alzheimer's disease (AD) neuropathology and were driven specifically by CA1 and subiculum atrophy. AD pathology, but not HS, was associated significantly with the rate of hippocampal atrophy. DISCUSSION: HS-associated volume changes are detectable on MRI earlier than 10 years before death. Based on these findings, volumetric cutoffs could be derived for in vivo differentiation between HS and AD. HIGHLIGHTS: Hippocampal atrophy was found in HS+ patients earlier than 10 years before death. These early pre-mortem changes were driven by reduced CA1 and subiculum volumes. Rates of hippocampus and subfield volume decline were independent of HS. In contrast, steeper atrophy rates were associated with AD pathology burden. Differentiation between AD and HS could be facilitated based on these MRI findings.

A novel histological staging of hippocampal sclerosis that is evident in gray matter loss in vivo.

INTRODUCTION: Hippocampal sclerosis of aging (HS) is defined by end-stage histological findings, strongly associated with limbic-predominant age-related TAR DNA-binding protein 43 (TDP-43) encephalopathy (LATE). We aimed to characterize features of early HS to refine the understanding of its role within combined pathology. METHODS: We studied 159 brain donations from the multimodal Vallecas Alzheimer's Center Study. A staging system (0 to IV) was developed to account for HS progression and analyzed in relation to pre-mortem cognitive and magnetic resonance imaging (MRI) data. RESULTS: Our HS staging system displayed a significant correlation with disease duration, cognitive performance, and combined neuropathologies, especially with LATE. Two-level assessment along the hippocampal longitudinal axis revealed an anterior-posterior gradient of HS severity. In vivo MRI showed focally reduced hippocampal gray matter density as a function of HS staging. DISCUSSION: The association of this staging system with clinical progression and structural differences supports its utility in the characterization and potential in vivo monitoring of HS. HIGHLIGHTS: The definition of hippocampal sclerosis of aging (HS) is currently limited to an end-stage pathological fingerprint. We characterize early HS histological features to define a complete staging system. The proposed staging displays a parallel but not identical progression to limbic-predominant age-related TAR DNA-binding protein 43 (TDP-43) encephalopathy (LATE). The proposed staging also reflects the expected demographic and cognitive differences associated with HS. In vivo magnetic resonance imaging (MRI) showed focal hippocampal gray matter loss as a function of HS staging.

High-resolution atlasing and segmentation of the subcortex: Review and perspective on challenges and opportunities created by machine learning.

This paper reviews almost three decades of work on atlasing and segmentation methods for subcortical structures in human brain MRI. In writing this survey, we have three distinct aims. First, to document the evolution of digital subcortical atlases of the human brain, from the early MRI templates published in the nineties, to the complex multi-modal atlases at the subregion level that are available today. Second, to provide a detailed record of related efforts in the automated segmentation front, from earlier atlas-based methods to modern machine learning approaches. And third, to present a perspective on the future of high-resolution atlasing and segmentation of subcortical structures in in vivo human brain MRI, including open challenges and opportunities created by recent developments in machine learning.

FreeSurfer-based segmentation of hippocampal subfields: A review of methods and applications, with a novel quality control procedure for ENIGMA studies and other collaborative efforts.

Structural hippocampal abnormalities are common in many neurological and psychiatric disorders, and variation in hippocampal measures is related to cognitive performance and other complex phenotypes such as stress sensitivity. Hippocampal subregions are increasingly studied, as automated algorithms have become available for mapping and volume quantification. In the context of the Enhancing Neuro Imaging Genetics through Meta Analysis Consortium, several Disease Working Groups are using the FreeSurfer software to analyze hippocampal subregion (subfield) volumes in patients with neurological and psychiatric conditions along with data from matched controls. In this overview, we explain the algorithm's principles, summarize measurement reliability studies, and demonstrate two additional aspects (subfield autocorrelation and volume/reliability correlation) with illustrative data. We then explain the rationale for a standardized hippocampal subfield segmentation quality control (QC) procedure for improved pipeline harmonization. To guide researchers to make optimal use of the algorithm, we discuss how global size and age effects can be modeled, how QC steps can be incorporated and how subfields may be aggregated into composite volumes. This discussion is based on a synopsis of 162 published neuroimaging studies (01/2013-12/2019) that applied the FreeSurfer hippocampal subfield segmentation in a broad range of domains including cognition and healthy aging, brain development and neurodegeneration, affective disorders, psychosis, stress regulation, neurotoxicity, epilepsy, inflammatory disease, childhood adversity and posttraumatic stress disorder, and candidate and whole genome (epi-)genetics. Finally, we highlight points where FreeSurfer-based hippocampal subfield studies may be optimized.

Perivascular space dilation is associated with vascular amyloid-ß accumulation in the overlying cortex.

Perivascular spaces (PVS) are compartments surrounding cerebral blood vessels that become visible on MRI when enlarged. Enlarged PVS (EPVS) are commonly seen in patients with cerebral small vessel disease (CSVD) and have been suggested to reflect dysfunctional perivascular clearance of soluble waste products from the brain. In this study, we investigated histopathological correlates of EPVS and how they relate to vascular amyloid-ß (Aß) in cerebral amyloid angiopathy (CAA), a form of CSVD that commonly co-exists with Alzheimer's disease (AD) pathology. We used ex vivo MRI, semi-automatic segmentation and validated deep-learning-based models to quantify EPVS and associated histopathological abnormalities. Severity of MRI-visible PVS during life was significantly associated with severity of MRI-visible PVS on ex vivo MRI in formalin fixed intact hemispheres and corresponded with PVS enlargement on histopathology in the same areas. EPVS were located mainly around the white matter portion of perforating cortical arterioles and their burden was associated with CAA severity in the overlying cortex. Furthermore, we observed markedly reduced smooth muscle cells and increased vascular Aß accumulation, extending into the WM, in individually affected vessels with an EPVS. Overall, these findings are consistent with the notion that EPVS reflect impaired outward flow along arterioles and have implications for our understanding of perivascular clearance mechanisms, which play an important role in the pathophysiology of CAA and AD.

A deep learning toolbox for automatic segmentation of subcortical limbic structures from MRI images.

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.

Joint super-resolution and synthesis of 1 mm isotropic MP-RAGE volumes from clinical MRI exams with scans of different orientation, resolution and contrast.

Most existing algorithms for automatic 3D morphometry of human brain MRI scans are designed for data with near-isotropic voxels at approximately 1 mm resolution, and frequently have contrast constraints as well-typically requiring T1-weighted images (e.g., MP-RAGE scans). This limitation prevents the analysis of millions of MRI scans acquired with large inter-slice spacing in clinical settings every year. In turn, the inability to quantitatively analyze these scans hinders the adoption of quantitative neuro imaging in healthcare, and also precludes research studies that could attain huge sample sizes and hence greatly improve our understanding of the human brain. Recent advances in convolutional neural networks (CNNs) are producing outstanding results in super-resolution and contrast synthesis of MRI. However, these approaches are very sensitive to the specific combination of contrast, resolution and orientation of the input images, and thus do not generalize to diverse clinical acquisition protocols - even within sites. In this article, we present SynthSR, a method to train a CNN that receives one or more scans with spaced slices, acquired with different contrast, resolution and orientation, and produces an isotropic scan of canonical contrast (typically a 1 mm MP-RAGE). The presented method does not require any preprocessing, beyond rigid coregistration of the input scans. Crucially, SynthSR trains on synthetic input images generated from 3D segmentations, and can thus be used to train CNNs for any combination of contrasts, resolutions and orientations without high-resolution real images of the input contrasts. We test the images generated with SynthSR in an array of common downstream analyses, and show that they can be reliably used for subcortical segmentation and volumetry, image registration (e.g., for tensor-based morphometry), and, if some image quality requirements are met, even cortical thickness morphometry. The source code is publicly available at https://github.com/BBillot/SynthSR.

Quantification of volumetric morphometry and optical property in the cortex of human cerebellum at micrometer resolution.

The surface of the human cerebellar cortex is much more tightly folded than the cerebral cortex. Volumetric analysis of cerebellar morphometry in magnetic resonance imaging studies suffers from insufficient resolution, and therefore has had limited impact on disease assessment. Automatic serial polarization-sensitive optical coherence tomography (as-PSOCT) is an emerging technique that offers the advantages of microscopic resolution and volumetric reconstruction of large-scale samples. In this study, we reconstructed multiple cubic centimeters of ex vivo human cerebellum tissue using as-PSOCT. The morphometric and optical properties of the cerebellar cortex across five subjects were quantified. While the molecular and granular layers exhibited similar mean thickness in the five subjects, the thickness varied greatly in the granular layer within subjects. Layer-specific optical property remained homogenous within individual subjects but showed higher cross-subject variability than layer thickness. High-resolution volumetric morphometry and optical property maps of human cerebellar cortex revealed by as-PSOCT have great potential to advance our understanding of cerebellar function and diseases.

Reliability and sensitivity of two whole-brain segmentation approaches included in FreeSurfer - ASEG and SAMSEG.

Accurate and reliable whole-brain segmentation is critical to longitudinal neuroimaging studies. We undertake a comparative analysis of two subcortical segmentation methods, Automatic Segmentation (ASEG) and Sequence Adaptive Multimodal Segmentation (SAMSEG), recently provided in the open-source neuroimaging package FreeSurfer 7.1, with regard to reliability, bias, sensitivity to detect longitudinal change, and diagnostic sensitivity to Alzheimer's disease. First, we assess intra- and inter-scanner reliability for eight bilateral subcortical structures: amygdala, caudate, hippocampus, lateral ventricles, nucleus accumbens, pallidum, putamen and thalamus. For intra-scanner analysis we use a large sample of participants (n = 1629) distributed across the lifespan (age range = 4-93 years) and acquired on a 1.5T Siemens Avanto (n = 774) and a 3T Siemens Skyra (n = 855) scanners. For inter-scanner analysis we use a sample of 24 participants scanned on the day with three models of Siemens scanners: 1.5T Avanto, 3T Skyra and 3T Prisma. Second, we test how each method detects volumetric age change using longitudinal follow up scans (n = 491 for Avanto and n = 245 for Skyra; interscan interval = 1-10 years). Finally, we test sensitivity to clinically relevant change. We compare annual rate of hippocampal atrophy in cognitively normal older adults (n = 20), patients with mild cognitive impairment (n = 20) and Alzheimer's disease (n = 20). We find that both ASEG and SAMSEG are reliable and lead to the detection of within-person longitudinal change, although with notable differences between age-trajectories for most structures, including hippocampus and amygdala. In summary, SAMSEG yields significantly lower differences between repeated measures for intra- and inter-scanner analysis without compromising sensitivity to changes and demonstrating ability to detect clinically relevant longitudinal changes.

Infant FreeSurfer: An automated segmentation and surface extraction pipeline for T1-weighted neuroimaging data of infants 0-2 years.

The development of automated tools for brain morphometric analysis in infants has lagged significantly behind analogous tools for adults. This gap reflects the greater challenges in this domain due to: 1) a smaller-scaled region of interest, 2) increased motion corruption, 3) regional changes in geometry due to heterochronous growth, and 4) regional variations in contrast properties corresponding to ongoing myelination and other maturation processes. Nevertheless, there is a great need for automated image-processing tools to quantify differences between infant groups and other individuals, because aberrant cortical morphologic measurements (including volume, thickness, surface area, and curvature) have been associated with neuropsychiatric, neurologic, and developmental disorders in children. In this paper we present an automated segmentation and surface extraction pipeline designed to accommodate clinical MRI studies of infant brains in a population 0-2 year-olds. The algorithm relies on a single channel of T1-weighted MR images to achieve automated segmentation of cortical and subcortical brain areas, producing volumes of subcortical structures and surface models of the cerebral cortex. We evaluated the algorithm both qualitatively and quantitatively using manually labeled datasets, relevant comparator software solutions cited in the literature, and expert evaluations. The computational tools and atlases described in this paper will be distributed to the research community as part of the FreeSurfer image analysis package.

Automated segmentation of the hypothalamus and associated subunits in brain MRI.

Despite the crucial role of the hypothalamus in the regulation of the human body, neuroimaging studies of this structure and its nuclei are scarce. Such scarcity partially stems from the lack of automated segmentation tools, since manual delineation suffers from scalability and reproducibility issues. Due to the small size of the hypothalamus and the lack of image contrast in its vicinity, automated segmentation is difficult and has been long neglected by widespread neuroimaging packages like FreeSurfer or FSL. Nonetheless, recent advances in deep machine learning are enabling us to tackle difficult segmentation problems with high accuracy. In this paper we present a fully automated tool based on a deep convolutional neural network, for the segmentation of the whole hypothalamus and its subregions from T1-weighted MRI scans. We use aggressive data augmentation in order to make the model robust to T1-weighted MR scans from a wide array of different sources, without any need for preprocessing. We rigorously assess the performance of the presented tool through extensive analyses, including: inter- and intra-rater variability experiments between human observers; comparison of our tool with manual segmentation; comparison with an automated method based on multi-atlas segmentation; assessment of robustness by quality control analysis of a larger, heterogeneous dataset (ADNI); and indirect evaluation with a volumetric study performed on ADNI. The presented model outperforms multi-atlas segmentation scores as well as inter-rater accuracy level, and approaches intra-rater precision. Our method does not require any preprocessing and runs in less than a second on a GPU, and approximately 10 seconds on a CPU. The source code as well as the trained model are publicly available at https://github.com/BBillot/hypothalamus_seg, and will also be distributed with FreeSurfer.

Sparsity enables estimation of both subcortical and cortical activity from MEG and EEG.

Subcortical structures play a critical role in brain function. However, options for assessing electrophysiological activity in these structures are limited. Electromagnetic fields generated by neuronal activity in subcortical structures can be recorded noninvasively, using magnetoencephalography (MEG) and electroencephalography (EEG). However, these subcortical signals are much weaker than those generated by cortical activity. In addition, we show here that it is difficult to resolve subcortical sources because distributed cortical activity can explain the MEG and EEG patterns generated by deep sources. We then demonstrate that if the cortical activity is spatially sparse, both cortical and subcortical sources can be resolved with M/EEG. Building on this insight, we develop a hierarchical sparse inverse solution for M/EEG. We assess the performance of this algorithm on realistic simulations and auditory evoked response data, and show that thalamic and brainstem sources can be correctly estimated in the presence of cortical activity. Our work provides alternative perspectives and tools for characterizing electrophysiological activity in subcortical structures in the human brain.

A probabilistic atlas of the human thalamic nuclei combining ex vivo MRI and histology.

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.

A magnetic resonance multi-atlas for the neonatal rabbit brain.

The rabbit model has become increasingly popular in neurodevelopmental studies as it is best suited to bridge the gap in translational research between small and large animals. In the context of preclinical studies, high-resolution magnetic resonance imaging (MRI) is often the best modality to investigate structural and functional variability of the brain, both in vivo and ex vivo. In most of the MRI-based studies, an important requirement to analyze the acquisitions is an accurate parcellation of the considered anatomical structures. Manual segmentation is time-consuming and typically poorly reproducible, while state-of-the-art automated segmentation algorithms rely on available atlases. In this work we introduce the first digital neonatal rabbit brain atlas consisting of 12 multi-modal acquisitions, parcellated into 89 areas according to a hierarchical taxonomy. Delineations were performed iteratively, alternating between segmentation propagation, label fusion and manual refinements, with the aim of controlling the quality while minimizing the bias introduced by the chosen sequence. Reliability and accuracy were assessed with cross-validation and intra- and inter-operator test-retests. Multi-atlas, versioned controlled segmentations repository and supplementary materials download links are available from the software repository documentation at https://github.com/gift-surg/SPOT-A-NeonatalRabbit.