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.The authors would like to thank Professor Karla Miller (Oxford) for her help with the design of the ex vivo MRI acquisition; Ms. Mercedes I~niguez de Onzo~no and Mr. Francisco Romero (UCLM) for their careful technical laboratory help; and Mr. Gonzalo Artacho (UCLM) for his help with the digitization and curation of his organization of histological data. This project has received funding from the European Union's Horizon 2020 research and innovation program under the Marie Sklodowska- Curie grant agreement No 654911 (project “THALAMODEL”) and by the European Research Council (ERC) Starting Grant agreement No 677697 (“BUNGEE-TOOLS”). It was also funded by the Spanish Ministry of Economy and Competitiveness(MINECO TEC-2014-51882-P, RYC- 2014-15440, PSI2015-65696, and SEV-2015-0490), the Basque Government (PI2016-12), and UCLM Internal Research Groups grants. Support for this research was also provided in part by the National Institute of Biomedical Imaging and Bioengineering (P41EB015896, 1R01EB023281, R01EB006758, R21EB018907, R01EB019956), the National Institute on Aging (5R01AG008122, R01AG016495), the National Institute of Diabetes and Digestive and Kidney Diseases (1-R21-DK- 108277-01), the National Institute of Neurological Disorders and Stroke (R01NS0525851, R21NS072652, R01NS070963, R01NS083534, 5U01NS086625), and was made possible by the resources provided by Shared Instrumentation Grants 1S10RR023401, 1S10RR019307, and 1S- 10RR023043. Additional support was provided by the NIH Blueprint for Neuroscience Research (5U01-MH093765), part of the multiinstitutional Human Connectome Project. In addition, B.F. has a financial interest in CorticoMetrics, a company whose medical pursuits focus on brain imaging and measurement technologies. B.F.’s interests were reviewed and are managed by Massachusetts General Hospital and Partners HealthCare in accordance with their conflict of interest policies. Data collection and sharing for this project was funded by the Alzheimer's Disease Neuroimaging Initiative (National Institutes of Health Grant U01 AG024904) and DOD ADNI (DOD award number W81XWH-12-2-0012). ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: AbbVie, Alzheimers Association; Alzheimers Drug Discovery Foundation; Araclon Biotech; BioClinica, Inc.; Biogen; Bristol-Myers Squibb Company; CereSpir, Inc.; Cogstate; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann-La Roche Ltd and its affiliated company Genentech, Inc.; Fujirebio; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Lumosity; Lundbeck; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Neurotrack Technologies; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (www.fnih.org). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimers Therapeutic Research Institute at the University of Southern California. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California

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 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

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

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

Data availability: Data will be made available on request.Supplementary materials are available online at: https://www.sciencedirect.com/science/article/pii/S1053811923002744?via%3Dihub#sec0024 .Research data are available online at: https://www.sciencedirect.com/science/article/pii/S1053811923002744?via%3Dihub#ec-research-data .Copyright © 2023 The Author(s). 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.This work was primarily funded by Alzheimers Research UK (ARUK-IRG2019A003). PGs work in this area was supported by NIH NIBIB NAC P41EB015902 AYs work in this area was supported by NIH grants R01 EB021265 and R56 MH121426. DCAs work in this area was supported by EPSRC grant EP/R006032/1 and Wellcome Trust award 221915/Z/20/Z. The Dementia Research Centre is supported by Alzheimer’s Research UK, Alzheimer’s Society, Brain Research UK, and The Wolfson Foundation. This work was supported by the National Institute for Health Research (NIHR) Queen Square Dementia Biomedical Research Unit and the University College London Hospitals Biomedical Research Centre, the Leonard Wolfson Experimental Neurology Centre (LWENC) Clinical Research Facility, and the UK Dementia Research Institute, which receives its funding from UK DRI Ltd, funded by the UK Medical Research Council, Alzheimer’s Society and Alzheimer’s Research UK. This project has received funding from the European Unions Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement No. 765148, as well as from the National Institutes Of Health under project number R01NS112161. MB is supported by a Fellowship award from the Alzheimers Society, UK (AS-JF-19a-004-517). MBs work was also supported by the UK Dementia Research Institute which receives its funding from DRI Ltd, funded by the UK Medical Research Council, Alzheimers Society and Alzheimers Research UK. JDR is supported by the Miriam Marks Brain Research UK Senior Fellowship and has received funding from an MRC Clinician Scientist Fellowship (MR/M008525/1) and the NIHR Rare Disease Translational Research Collaboration (BRC149/NS/MH). JEI is supported by the European Research Council (Starting Grant 677697, project BUNGEE-TOOLS) and the NIH (1RF1MH123195-01 and 1R01AG070988-01). The collection and sharing of the ADNI data was funded by the Alzheimer’s Disease Neuroimaging Initiative (National Institutes of Health Grant U01 AG024904) and Department of Defence (W81XWH-12-2-0012). ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and the following: Alzheimer’s Association; Alzheimer’s Drug Discovery Foundation; BioClinica, Inc.; Biogen Idec Inc.; Bristol-Myers Squibb Company; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; F. Hoffmann-La Roche Ltd and affiliated company Genentech, Inc.; GE Healthcare; Innogenetics, N.V.; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Medpace, Inc.; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Synarc Inc.; and Takeda Pharmaceutical Company. The Canadian Institutes of Health Research is providing funds for ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health. The grantee is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer’s Disease Cooperative Study at the University of California, San Diego. ADNI is disseminated by the Laboratory for Neuro Imaging at the University of Southern California

Systematic validation of an automated thalamic parcellation technique using anatomical data at 3T

The thalamus is a brain region formed from functionally distinct nuclei, which contribute in important ways to various cognitive processes. Yet, much of the human neuroscience literature treats the thalamus as one homogeneous region, and consequently the unique contribution of specific nuclei to behaviour remains under-appreciated. This is likely due in part to the technical challenge of dissociating nuclei using conventional structural imaging approaches. Yet, multiple algorithms exist in the neuroimaging literature for the automated segmentation of thalamic nuclei. One recent approach developed by Iglesias and colleagues (2018) generates segmentations by applying a probabilistic atlas to subject-space anatomical images using the FreeSurfer software. Here, we systematically validate the efficacy of this segmentation approach in delineating thalamic nuclei using Human Connectome Project data. We provide several metrics quantifying the quality of segmentations relative to the Morel stereotaxic atlas, a widely accepted anatomical atlas based on cyto- and myeloarchitecture. The automated segmentation approach generated boundaries between the anterior, lateral, posterior, and medial divisions of the thalamus. Segmentation efficacy, as measured by metrics of dissimilarity (Average Hausdorff Distance) and overlap (DICE coefficient) within groups was mixed. Regions were better delineated in anterior, lateral and medial thalamus than the posterior thalamus, however all the volumes for all segmented nuclei were significantly different to the corresponding region of the Morel atlas. These mixed results suggest users should exercise care when using this approach to study the structural or functional relevance of a given thalamic nucleus

Cortical and subcortical contributions to human cognitive flexibility

Cognitive flexibility enables individuals to respond adaptively to an ever-changing world. Neurally, flexibility is underpinned by involvement from across the cerebrum, and there is evidence from animal and human neuroscience suggesting that integration of cortical and thalamic signals in the striatum is necessary for appropriate behavioural control. A commonly used assay of flexibility is reversal learning, an associative learning task with high inter-species translatability. Evidence from animal literature has clearly defined the importance of the striatal cholinergic system in regulating striatal activity and output from the basal ganglia, and there is nascent evidence suggesting this system operates in a similar way in humans. However, there is a need to further disentangle the role of cortical, striatal, and thalamic regions during reversal learning in humans to better understand how the system works, and whether it has heterogeneous functionality in different contexts. Furthermore, as studying these processes is not trivial, further methodological work is required to enable us to understand the system. In chapter two we systematically assess an automated parcellation technique for identifying specific thalamic nuclei. Despite generally being treated as a homologous structure in neuroimaging work, nuclei within the thalamus have dissociable roles, and have diverse contributions to cognitive functioning, including reversal learning. We found mixed efficacy for segmentations across the thalamus, with some regions being more accurately defined relative to a “gold standard” atlas than others. Crucially, we find that the centromedian and parafascicular nuclei, which have an important role in reversal learning, are clearly defined and have little overlap with contiguous regions. These results show we can use this automated parcellation technique to identify specific thalamic nuclei that are relevant for cognitive flexibility and use these parcellations to study functionally relevant processes. Recent work has demonstrated that the functional relevance of the striatal cholinergic system can be studied in vivo using magnetic resonance spectroscopy by separating the peaks of different metabolites. But this non-conventional approach has not yet been widely adopted, and work is needed to determine its reliability. Chapter three presents test-retest reliability data on the use of magnetic resonance spectroscopy to study cholinergic activity in the striatum and cortex. We find measures of choline containing compounds are highly correlated when peaks are separated and when they are not. Across time we find that choline concentrations are relatively inconsistent, and that this was due to changes in the functionally relevant metabolite choline. Conversely, metabolites that we think are not functionally relevant were stable over time. We believe these differences may underly differences in acetylcholine function over time and may explain some intra-individual behavioural variability. In chapter four we use functional magnetic resonance imaging and psychophysiological interaction analysis to study corticostriatal and thalamostriatal connectivity during serial reversal learning. Functional connectivity between the centromedian-parafascicular nuclei of the thalamus and the associative dorsal striatum, and between the lateral-orbitofrontal cortex and the associative dorsal striatum was related to processing feedback during reversal learning. Specifically, thalamostriatal connectivity was found across the task, and may reflect a general error signal used to identify potential changes in context. Conversely, corticostriatal connectivity was found to be specific to when behaviour changed and suggests this may be a mechanism for the implementing adaptive change. We also show findings from exploratory work that may explain further how the cortex supports flexibility during reversal learning. Lastly, we used magnetic resonance spectroscopy to investigate whether the state of the cholinergic system at rest is related to reversal learning performance and latent measures of behaviour using computational modelling. Choline concentrations at rest showed significant functional relevance to our measures of reversal learning. More specifically, we found that errors during reversal learning, and learning rates for positive and negative prediction errors, explained significant variance in choline. However, the relationship between choline levels and task performance presented here differ from previous work which instead used a multi-alternative reversal learning task, and suggests that the striatal cholinergic system may have dissociable roles in different contexts. Overall, we show that the striatum, its cholinergic interneuron system, and its afferent projections from the cortex and thalamus, are associated with performance during serial reversal learning. Moreover, these findings suggest that the system may operate in separable ways in different contexts which may be dependent on internal representations of task structure

Imaging the subthalamic nucleus in Parkinson’s disease

This thesis is comprised of a set of work that aims to visualize and quantify the anatomy, structural variability, and connectivity of the subthalamic nucleus (STN) with optimized neuroimaging methods. The study populations include both healthy cohorts and individuals living with Parkinson's disease (PD). PD was chosen specifically due to the involvement of the STN in the pathophysiology of the disease. Optimized neuroimaging methods were primarily obtained using ultra-high field (UHF) magnetic resonance imaging (MRI). An additional component of this thesis was to determine to what extent UHF-MRI can be used in a clinical setting, specifically for pre-operative planning of deep brain stimulation (DBS) of the STN for patients with advanced PD. The thesis collectively demonstrates that i, MRI research, and clinical applications must account for the different anatomical and structural changes that occur in the STN with both age and PD. ii, Anatomical connections involved in preparatory motor control, response inhibition, and decision-making may be compromised in PD. iii. The accuracy of visualizing and quantifying the STN strongly depends on the type of MR contrast and voxel size. iv, MRI at a field strength of 3 Tesla (T) can under certain circumstances be optimized to produce results similar to that of 7 T at the expense of increased acquisition time

A computational atlas of the hippocampal formation using ex vivo, ultra-high resolution MRI: Application to adaptive segmentation of in vivo MRI.

AbstractAutomated analysis of MRI data of the subregions of the hippocampus requires computational atlases built at a higher resolution than those that are typically used in current neuroimaging studies. Here we describe the construction of a statistical atlas of the hippocampal formation at the subregion level using ultra-high resolution, ex vivo MRI. Fifteen autopsy samples were scanned at 0.13mm isotropic resolution (on average) using customized hardware. The images were manually segmented into 13 different hippocampal substructures using a protocol specifically designed for this study; precise delineations were made possible by the extraordinary resolution of the scans. In addition to the subregions, manual annotations for neighboring structures (e.g., amygdala, cortex) were obtained from a separate dataset of in vivo, T1-weighted MRI scans of the whole brain (1mm resolution). The manual labels from the in vivo and ex vivo data were combined into a single computational atlas of the hippocampal formation with a novel atlas building algorithm based on Bayesian inference. The resulting atlas can be used to automatically segment the hippocampal subregions in structural MRI images, using an algorithm that can analyze multimodal data and adapt to variations in MRI contrast due to differences in acquisition hardware or pulse sequences. The applicability of the atlas, which we are releasing as part of FreeSurfer (version 6.0), is demonstrated with experiments on three different publicly available datasets with different types of MRI contrast. The results show that the atlas and companion segmentation method: 1) can segment T1 and T2 images, as well as their combination, 2) replicate findings on mild cognitive impairment based on high-resolution T2 data, and 3) can discriminate between Alzheimer's disease subjects and elderly controls with 88% accuracy in standard resolution (1mm) T1 data, significantly outperforming the atlas in FreeSurfer version 5.3 (86% accuracy) and classification based on whole hippocampal volume (82% accuracy)

Multimodal MRI analysis using deep learning methods

Magnetic resonance imaging (MRI) has been widely used in scientific and clinical research. It is a non-invasive medical imaging technique that reveals anatomical structures and provides useful information for investigators to explore aging and pathological processes. Different MR modalities offer different useful properties. Automatic MRI analysis algorithms have been developed to address problems in many applications such as classification, segmentation, and disease diagnosis. Segmentation and labeling algorithms applied to brain MRIs enable evaluations of the volumetric changes of specific structures in neurodegenerative diseases. Reconstruction of fiber orientations using diffusion MRI is beneficial to obtain better understanding of the underlying structures. In this thesis, we focused on development of deep learning methods for MRI analysis using different image modalities. Specifically, we applied deep learning techniques on different applications, including segmentation of brain structures and reconstruction of tongue muscle fiber orientations. For segmentation of brain structures, we developed an end-to-end deep learning algorithm for ventricle parcellation of brains with ventriculomegaly using T1-w MR images. The deep network provides robust and accurate segmentation results in subjects with high variability in ventricle shapes and sizes. We developed another deep learning method to automatically parcellate the thalamus into a set of thalamic nuclei using T1-w MRI and features from diffusion MRI. The algorithm incorporates a harmonization step to make the network adapt to input images with different contrasts. We also studied the strains associated with tongue muscles during speech production using multiple MRI modalities. To enable this study, we first developed a deep network to reconstruct crossing tongue muscle fiber orientations using diffusion MRI. The network was specifically designed for the human tongue and accounted for the orthogonality property of the tongue muscles. Next, we proposed a comprehensive pipeline to analyze the strains associated with tongue muscle fiber orientations during speech using diffusion MRI, and tagged and cine MRI. The proposed pipeline provides a solution to analyze the cooperation between muscle groups during speech production