Thalamus segmentation based on the local diffusion direction: a group study

Robust thalamic nuclei segmentation from T1-weighted MRI using polynomial intensity transformation

Accurate segmentation of thalamic nuclei, crucial for understanding their role in healthy cognition and in pathologies, is challenging to achieve on standard T1-weighted (T1w) magnetic resonance imaging (MRI) due to poor image contrast. White-matter-nulled (WMn) MRI sequences improve intrathalamic contrast but are not part of clinical protocols or extant databases. In this study, we introduce histogram-based polynomial synthesis (HIPS), a fast preprocessing transform step that synthesizes WMn-like image contrast from standard T1w MRI using a polynomial approximation for intensity transformation. HIPS was incorporated into THalamus Optimized Multi-Atlas Segmentation (THOMAS) pipeline, a method developed and optimized for WMn MRI. HIPS-THOMAS was compared to a convolutional neural network (CNN)-based segmentation method and THOMAS modified for the use of T1w images (T1w-THOMAS). The robustness and accuracy of the three methods were tested across different image contrasts (MPRAGE, SPGR, and MP2RAGE), scanner manufacturers (PHILIPS, GE, and Siemens), and field strengths (3 T and 7 T). HIPS-transformed images improved intra-thalamic contrast and thalamic boundaries, and HIPS-THOMAS yielded significantly higher mean Dice coefficients and reduced volume errors compared to both the CNN method and T1w-THOMAS. Finally, all three methods were compared using the frequently travelling human phantom MRI dataset for inter- and intra-scanner variability, with HIPS displaying the least inter-scanner variability and performing comparably with T1w-THOMAS for intra-scanner variability. In conclusion, our findings highlight the efficacy and robustness of HIPS in enhancing thalamic nuclei segmentation from standard T1w MRI.

DiMANI: diffusion MRI for anatomical nuclei imaging-Application for the direct visualization of thalamic subnuclei

The thalamus is a centrally located and heterogeneous brain structure that plays a critical role in various sensory, motor, and cognitive processes. However, visualizing the individual subnuclei of the thalamus using conventional MRI techniques is challenging. This difficulty has posed obstacles in targeting specific subnuclei for clinical interventions such as deep brain stimulation (DBS). In this paper, we present DiMANI, a novel method for directly visualizing the thalamic subnuclei using diffusion MRI (dMRI). The DiMANI contrast is computed by averaging, voxelwise, diffusion-weighted volumes enabling the direct distinction of thalamic subnuclei in individuals. We evaluated the reproducibility of DiMANI through multiple approaches. First, we utilized a unique dataset comprising 8 scans of a single participant collected over a 3-year period. Secondly, we quantitatively assessed manual segmentations of thalamic subnuclei for both intra-rater and inter-rater reliability. Thirdly, we qualitatively correlated DiMANI imaging data from several patients with Essential Tremor with the localization of implanted DBS electrodes and clinical observations. Lastly, we demonstrated that DiMANI can provide similar features at 3T and 7T MRI, using varying numbers of diffusion directions. Our results establish that DiMANI is a reproducible and clinically relevant method to directly visualize thalamic subnuclei. This has significant implications for the development of new DBS targets and the optimization of DBS therapy.

Robust thalamic nuclei segmentation from T1-weighted MRI using polynomial intensity transformation

Accurate segmentation of thalamic nuclei, crucial for understanding their role in healthy cognition and in pathologies, is challenging to achieve on standard T1-weighted (T1w) magnetic resonance imaging (MRI) due to poor image contrast. White-matter-nulled (WMn) MRI sequences improve intrathalamic contrast but are not part of clinical protocols or extant databases. In this study, we introduce histogram-based polynomial synthesis (HIPS), a fast preprocessing transform step that synthesizes WMn-like image contrast from standard T1w MRI using a polynomial approximation for intensity transformation. HIPS was incorporated into THalamus Optimized Multi-Atlas Segmentation (THOMAS) pipeline, a method developed and optimized for WMn MRI. HIPS-THOMAS was compared to a convolutional neural network (CNN)-based segmentation method and THOMAS modified for T1w images (T1w-THOMAS). The robustness and accuracy of the three methods were tested across different image contrasts (MPRAGE, SPGR, and MP2RAGE), scanner manufacturers (PHILIPS, GE, and Siemens), and field strengths (3T and 7T). HIPS-transformed images improved intra-thalamic contrast and thalamic boundaries, and HIPS-THOMAS yielded significantly higher mean Dice coefficients and reduced volume errors compared to both the CNN method and T1w-THOMAS. Finally, all three methods were compared using the frequently travelling human phantom MRI dataset for inter- and intra-scanner variability, with HIPS displaying the least inter-scanner variability and performing comparably with T1w-THOMAS for intra-scanner variability. In conclusion, our findings highlight the efficacy and robustness of HIPS in enhancing thalamic nuclei segmentation from standard T1w MRI.

The effect of ibudilast on thalamic volume in progressive multiple sclerosis

BACKGROUND

Thalamic volume loss is known to be associated with clinical and cognitive disability in progressive multiple sclerosis (PMS).

OBJECTIVE

To investigate the treatment effect of ibudilast on thalamic atrophy more than 96 weeks in the phase 2 trial in progressive(MS Secondary and Primary Progressive Ibudilast NeuroNEXT Trial in Multiple Sclerosis [SPRINT-MS]).

METHODS

A total of 231 participants were randomized to either ibudilast (n = 114) or placebo (n = 117). Thalamic volume change was computed using Bayesian Sequence Adaptive Multimodal Segmentation tool (SAMseg) incorporating T1, fluid-attenuated inversion recovery (FLAIR), and fractional anisotropy maps and analyzed with a mixed-effects repeated-measures model.

RESULTS

There was no significant difference in thalamic volumes between treatment groups. On exploratory analysis, participants with primary progressive multiple sclerosis (PPMS) on placebo had a 0.004% greater rate of thalamic atrophy than PPMS participants on ibudilast (p = 0.058, 95% confidence interval (CI) = -0.008 to <0.001). Greater reductions in thalamic volumes at more than 96 weeks were associated with worsening multiple sclerosis functional composite (MSFC-4) scores (p = 0.002) and worsening performance on the symbol digit modality test (SDMT) (p < 0.001).

CONCLUSION

In a phase 2 trial evaluating ibudilast in PMS, no treatment effect was demonstrated in preventing thalamic atrophy. Participants with PPMS exhibited a treatment effect that trended toward significance. Longitudinal changes in thalamic volume were related to worsening of physical and cognitive disability, highlighting this outcome's clinical importance.

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

Neuroimaging at 7 Tesla: a pictorial narrative review

Neuroimaging using the 7-Tesla (7T) human magnetic resonance (MR) system is rapidly gaining popularity after being approved for clinical use in the European Union and the USA. This trend is the same for functional MR imaging (MRI). The primary advantages of 7T over lower magnetic fields are its higher signal-to-noise and contrast-to-noise ratios, which provide high-resolution acquisitions and better contrast, making it easier to detect lesions and structural changes in brain disorders. Another advantage is the capability to measure a greater number of neurochemicals by virtue of the increased spectral resolution. Many structural and functional studies using 7T have been conducted to visualize details in the white matter and layers of the cortex and hippocampus, the subnucleus or regions of the putamen, the globus pallidus, thalamus and substantia nigra, and in small structures, such as the subthalamic nucleus, habenula, perforating arteries, and the perivascular space, that are difficult to observe at lower magnetic field strengths. The target disorders for 7T neuroimaging range from tumoral diseases to vascular, neurodegenerative, and psychiatric disorders, including Alzheimer's disease, Parkinson's disease, multiple sclerosis, epilepsy, major depressive disorder, and schizophrenia. MR spectroscopy has also been used for research because of its increased chemical shift that separates overlapping peaks and resolves neurochemicals more effectively at 7T than a lower magnetic field. This paper presents a narrative review of these topics and an illustrative presentation of images obtained at 7T. We expect 7T neuroimaging to provide a new imaging biomarker of various brain disorders.

Diffusion tensor magnetic resonance imaging: is it valuable in the detection of brain microstructural changes in patients having migraine without aura?

Purpose

The aim of this study is to assess the diagnostic value of diffusion tensor magnetic resonance imaging (MRI) in the detection of brain microstructural changes in patients having migraine without aura.

Material and methods

Our prospective study included 33 patients having migraine without aura and 15 volunteers with matched age and sex, who underwent brain MRI with diffusion tensor imaging (DTI). The fractional anisotropy (FA) and mean diffusivity (MD) of selected grey and white matter regions on both sides were measured and correlated with the neurological clinical examination.

Results

Significant differences were detected in MD values in the thalamus, globus pallidus, and hippocampus head on the right side of patients versus controls. Also, significant differences of the FA values were detected in the thalamus, globus pallidus, and hippocampus head on the right side of patients versus controls. Regarding the FA values of the same regions on the left side, a significant difference in the FA value was detected only in the hippocampus head. There was a statistically significant difference in the FA values on both sides of the white matter of the frontal lobes, posterior limbs of the internal capsules, and cerebellar hemispheres in patients compared to controls. There was a statistically significant difference in MD values in the white matter of both frontal lobes, posterior limb of the right internal capsule, and both cerebellar hemispheres in patients compared to controls.

Conclusions

DTI can detect microstructural changes of the grey and white matter in patients having migraine without aura that could not be detected by conventional MRI.

In vivo high-resolution structural MRI-based atlas of human thalamic nuclei

Thalamic nuclei play critical roles in regulation of neurological functions like sleep and wakefulness. They are increasingly implicated in neurodegenerative and neurological diseases such as multiple sclerosis and essential tremor. However, segmentation of thalamic nuclei is difficult due to their poor visibility in conventional MRI scans. Sophisticated methods have been proposed which require specialized MRI acquisitions and complex post processing. There are few high spatial resolution (1 mm3 or higher) in vivo MRI thalamic atlases available currently. The goal of this work is the development of an in vivo MRI-based structural thalamic atlas at 0.7 × 0.7 × 0.5 mm resolution based on manual segmentation of 9 healthy subjects using the Morel atlas as a guide. Using data analysis from healthy subjects as well as patients with multiple-sclerosis and essential tremor and at 3T and 7T MRI, we demonstrate the utility of this atlas to provide fast and accurate segmentation of thalamic nuclei when only conventional T1 weighted images are available.

Structural assessment of thalamus morphology in brain disorders: A review and recommendation of thalamic nucleus segmentation and shape analysis

The thalamus is a central brain structure crucially involved in cognitive, emotional, sensory, and motor functions and is often reported to be involved in the pathophysiology of neurological and psychiatric disorders. The functional subdivision of the thalamus warrants morphological investigation on the level of individual subnuclei. In addition to volumetric measures, the investigation of other morphological features may give additional insights into thalamic morphology. For instance, shape features offer a higher spatial resolution by revealing small, regional differences that are left undetected in volumetric analyses. In this review, we discuss the benefits and limitations of recent advances in neuroimaging techniques to investigate thalamic morphology in vivo, leading to our proposed methodology. This methodology consists of available pipelines for volume and shape analysis, focussing on the morphological features of volume, thickness, and surface area. We demonstrate this combined approach in a Parkinson's disease cohort to illustrate their complementarity. Considering our findings, we recommend a combined methodology as it allows for more sensitive investigation of thalamic morphology in clinical populations.

Structural Changes in Thalamic Nuclei Across Prodromal and Clinical Alzheimer's Disease

BACKGROUND

Increasing evidence suggests that thalamic nuclei may atrophy in Alzheimer's disease (AD). We hypothesized that there will be significant atrophy of limbic thalamic nuclei associated with declining memory and cognition across the AD continuum.

OBJECTIVE

The objective of this work was to characterize volume differences in thalamic nuclei in subjects with early and late mild cognitive impairment (MCI) as well as AD when compared to healthy control (HC) subjects using a novel MRI-based thalamic segmentation technique (THOMAS).

METHODS

MPRAGE data from the ADNI database were used in this study (n = 540). Healthy control (n = 125), early MCI (n = 212), late MCI (n = 114), and AD subjects (n = 89) were selected, and their MRI data were parcellated to determine the volumes of 11 thalamic nuclei for each subject. Volumes across the different clinical subgroups were compared using ANCOVA.

RESULTS

There were significant differences in thalamic nuclei volumes between HC, late MCI, and AD subjects. The anteroventral, mediodorsal, pulvinar, medial geniculate, and centromedian nuclei were significantly smaller in subjects with late MCI and AD when compared to HC subjects. Furthermore, the mediodorsal, pulvinar, and medial geniculate nuclei were significantly smaller in early MCI when compared to HC subjects.

CONCLUSION

This work highlights nucleus specific atrophy within the thalamus in subjects with early and late MCI and AD. This is consistent with the hypothesis that memory and cognitive changes in AD are mediated by damage to a large-scale integrated neural network that extends beyond the medial temporal lobes.

Individualized Thalamic Parcellation Reveals Alterations in Shape and Microstructure of Thalamic Nuclei in Patients with Disorder of Consciousness

The thalamus plays crucial roles in consciousness generation and information processing. Previous evidence suggests that disorder of consciousness (DOC) caused by severe brain injury, is potentially related to thalamic abnormalities. However, how the morphology and microstructure change in thalamic subfields and thalamocortical fiber pathways in patients with DOC, and the relationships between these changes and the consciousness status remain unclear. Here, we generated the individual-specific thalamic parcellation in 10 DOC patients and 10 healthy controls (HC) via a novel thalamic segmentation framework based on the fiber orientation distribution (FOD) derived from 7-Tesla diffusion MRI, and investigated the shape deformation of thalamic nuclei as well as the microstructural changes associated with thalamic nuclei and thalamocortical pathways in patients with DOC. Enlargement of dorsal posterior nucleus and atrophy of anterior nucleus in the right thalamus were observed in DOC cohort relative to the HCs, and the former was closely linked to the consciousness level of the patients. We also found significant reductions of fiber density, but not fiber bundle cross-section, within several thalamic nuclei and most of the thalamocortical fiber pathways, suggesting that loss of axons might take primary responsibility for the impaired thalamocortical connections in patients with DOC rather than the change in fiber-bundle morphology. Furthermore, the individual-specific thalamic parcellation achieved 80% accuracy in classifying patients at the minimally conscious state from the vegetative state, compared with ~60% accuracy based on group-level parcellations. Our findings provide the first evidence for the shape deformation of thalamic nuclei in DOC patients and the microstructural basis of the disrupted thalamocortical connections.

Performance ramifications of abnormal functional connectivity of ventral posterior lateral thalamus with cerebellum in abstinent individuals with Alcohol Use Disorder

The extant literature supports the involvement of the thalamus in the cognitive and motor impairment associated with chronic alcohol consumption, but clear structure/function relationships remain elusive. Alcohol effects on specific nuclei rather than the entire thalamus may provide the basis for differential cognitive and motor decline in Alcohol Use Disorder (AUD). This functional MRI (fMRI) study was conducted in 23 abstinent individuals with AUD and 27 healthy controls to test the hypothesis that functional connectivity between anterior thalamus and hippocampus would be compromised in those with an AUD diagnosis and related to mnemonic deficits. Functional connectivity between 7 thalamic structures [5 thalamic nuclei: anterior ventral (AV), mediodorsal (MD), pulvinar (Pul), ventral lateral posterior (VLP), and ventral posterior lateral (VPL); ventral thalamus; the entire thalamus] and 14 "functional regions" was evaluated. Relative to controls, the AUD group exhibited different VPL-based functional connectivity: an anticorrelation between VPL and a bilateral middle temporal lobe region observed in controls became a positive correlation in the AUD group; an anticorrelation between the VPL and the cerebellum was stronger in the AUD than control group. AUD-associated altered connectivity between anterior thalamus and hippocampus as a substrate of memory compromise was not supported; instead, connectivity differences from controls selective to VPL and cerebellum demonstrated a relationship with impaired balance. These preliminary findings support substructure-level evaluation in future studies focused on discerning the role of the thalamus in AUD-associated cognitive and motor deficits.

Fast automatic segmentation of thalamic nuclei from MP2RAGE acquisition at 7 Tesla

PURPOSE

Thalamic nuclei are largely invisible in conventional MRI due to poor contrast. Thalamus Optimized Multi-Atlas Segmentation (THOMAS) provides automatic segmentation of 12 thalamic nuclei using white-matter-nulled (WMn) Magnetization Prepared Rapid Gradient Echo (MPRAGE) sequence at 7T, but increases overall scan duration. Routinely acquired, bias-corrected Magnetization Prepared 2 Rapid Gradient Echo (MP2RAGE) sequence yields superior tissue contrast and quantitative T1 maps. Application of THOMAS to MP2RAGE has been investigated in this study.

METHODS

Eight healthy volunteers and five pediatric-onset multiple sclerosis patients were recruited at Children's Hospital of Philadelphia and scanned at Siemens 7T with WMn-MPRAGE and multi-echo-MP2RAGE (ME-MP2RAGE) sequences. White-matter-nulled contrast was synthesized (MP2-SYN) from T₁ maps from ME-MP2RAGE sequence. Thalamic nuclei were segmented using THOMAS joint label fusion algorithm from WMn-MPRAGE and MP2-SYN datasets. THOMAS pipeline was modified to use majority voting to segment bias corrected T1-weighted uniform (MP2-UNI) images. Thalamic nuclei from MP2-SYN and MP2-UNI images were evaluated against corresponding nuclei obtained from WMn-MPRAGE images using dice coefficients, volume similarity indices (VSIs) and distance between centroids.

RESULTS

For MP2-SYN, dice > 0.85 and VSI > 0.95 was achieved for five larger nuclei and dice > 0.6 and VSI > 0.7 was achieved for seven smaller nuclei. The dice and VSI were slightly higher, whereas the distance between centroids were smaller for MP2-SYN compared to MP2-UNI, indicating improved performance using the MP2-SYN image.

CONCLUSIONS

THOMAS algorithm can successfully segment thalamic nuclei in MP2RAGE images with essentially equivalent quality as WMn-MPRAGE, widening its applicability in studies focused on thalamic involvement in aging and disease.

Automated thalamic nuclei segmentation using multi-planar cascaded convolutional neural networks

PURPOSE

To develop a fast and accurate convolutional neural network based method for segmentation of thalamic nuclei.

METHODS

A cascaded multi-planar scheme with a modified residual U-Net architecture was used to segment thalamic nuclei on conventional and white-matter-nulled (WMn) magnetization prepared rapid gradient echo (MPRAGE) data. A single network was optimized to work with images from healthy controls and patients with multiple sclerosis (MS) and essential tremor (ET), acquired at both 3 T and 7 T field strengths. WMn-MPRAGE images were manually delineated by a trained neuroradiologist using the Morel histological atlas as a guide to generate reference ground truth labels. Dice similarity coefficient and volume similarity index (VSI) were used to evaluate performance. Clinical utility was demonstrated by applying this method to study the effect of MS on thalamic nuclei atrophy.

RESULTS

Segmentation of each thalamus into twelve nuclei was achieved in under a minute. For 7 T WMn-MPRAGE, the proposed method outperforms current state-of-the-art on patients with ET with statistically significant improvements in Dice for five nuclei (increase in the range of 0.05-0.18) and VSI for four nuclei (increase in the range of 0.05-0.19), while performing comparably for healthy and MS subjects. Dice and VSI achieved using 7 T WMn-MPRAGE data are comparable to those using 3 T WMn-MPRAGE data. For conventional MPRAGE, the proposed method shows a statistically significant Dice improvement in the range of 0.14-0.63 over FreeSurfer for all nuclei and disease types. Effect of noise on network performance shows robustness to images with SNR as low as half the baseline SNR. Atrophy of four thalamic nuclei and whole thalamus was observed for MS patients compared to healthy control subjects, after controlling for the effect of parallel imaging, intracranial volume, gender, and age (p < 0.004).

CONCLUSION

The proposed segmentation method is fast, accurate, performs well across disease types and field strengths, and shows great potential for improving our understanding of thalamic nuclei involvement in neurological diseases.

A systematic comparison of structural-, structural connectivity-, and functional connectivity-based thalamus parcellation techniques

The thalamus consists of several histologically and functionally distinct nuclei increasingly implicated in brain pathology and important for treatment, motivating the need for development of fast and accurate thalamic parcellation. The contrast between thalamic nuclei as well as between the thalamus and surrounding tissues is poor in T1- and T2-weighted magnetic resonance imaging (MRI), inhibiting efforts to date to segment the thalamus using standard clinical MRI. Automatic parcellation techniques have been developed to leverage thalamic features better captured by advanced MRI methods, including magnetization prepared rapid acquisition gradient echo (MP-RAGE), diffusion tensor imaging (DTI), and resting-state functional MRI (fMRI). Despite operating on fundamentally different image contrasts, these methods claim a high degree of agreement with the Morel stereotactic atlas of the thalamus. However, no comparison has been undertaken to compare the results of these disparate parcellation methods. We have implemented state-of-the-art structural-, diffusion-, and functional imaging-based thalamus parcellation techniques and used them on a single set of subjects. We present the first systematic qualitative and quantitative comparison of these methods. The results show that DTI parcellation agrees more with structural parcellation in the larger thalamic nuclei, while rsfMRI parcellation agrees more with structural parcellation in the smaller nuclei. Structural parcellation is the most accurate in the delineation of small structures such as the habenular, antero-ventral, and medial geniculate nuclei.

Quantifying deep grey matter atrophy using automated segmentation approaches: A systematic review of structural MRI studies

The deep grey matter (DGM) nuclei of the brain play a crucial role in learning, behaviour, cognition, movement and memory. Although automated segmentation strategies can provide insight into the impact of multiple neurological conditions affecting these structures, such as Multiple Sclerosis (MS), Huntington's disease (HD), Alzheimer's disease (AD), Parkinson's disease (PD) and Cerebral Palsy (CP), there are a number of technical challenges limiting an accurate automated segmentation of the DGM. Namely, the insufficient contrast of T1 sequences to completely identify the boundaries of these structures, as well as the presence of iso-intense white matter lesions or extensive tissue loss caused by brain injury. Therefore in this systematic review, 269 eligible studies were analysed and compared to determine the optimal approaches for addressing these technical challenges. The automated approaches used among the reviewed studies fall into three broad categories, atlas-based approaches focusing on the accurate alignment of atlas priors, algorithmic approaches which utilise intensity information to a greater extent, and learning-based approaches that require an annotated training set. Studies that utilise freely available software packages such as FIRST, FreeSurfer and LesionTOADS were also eligible, and their performance compared. Overall, deep learning approaches achieved the best overall performance, however these strategies are currently hampered by the lack of large-scale annotated data. Improving model generalisability to new datasets could be achieved in future studies with data augmentation and transfer learning. Multi-atlas approaches provided the second-best performance overall, and may be utilised to construct a "silver standard" annotated training set for deep learning. To address the technical challenges, providing robustness to injury can be improved by using multiple channels, highly elastic diffeomorphic transformations such as LDDMM, and by following atlas-based approaches with an intensity driven refinement of the segmentation, which has been done with the Expectation Maximisation (EM) and level sets methods. Accounting for potential lesions should be achieved with a separate lesion segmentation approach, as in LesionTOADS. Finally, to address the issue of limited contrast, R2*, T2* and QSM sequences could be used to better highlight the DGM due to its higher iron content. Future studies could look to additionally acquire these sequences by retaining the phase information from standard structural scans, or alternatively acquiring these sequences for only a training set, allowing models to learn the "improved" segmentation from T1-sequences alone.

Thalamus Optimized Multi Atlas Segmentation (THOMAS): fast, fully automated segmentation of thalamic nuclei from structural MRI

The thalamus and its nuclei are largely indistinguishable on standard T1 or T2 weighted MRI. While diffusion tensor imaging based methods have been proposed to segment the thalamic nuclei based on the angular orientation of the principal diffusion tensor, these are based on echo planar imaging which is inherently limited in spatial resolution and suffers from distortion. We present a multi-atlas segmentation technique based on white-matter-nulled MP-RAGE imaging that segments the thalamus into 12 nuclei with computation times on the order of 10 min on a desktop PC; we call this method THOMAS (THalamus Optimized Multi Atlas Segmentation). THOMAS was rigorously evaluated on 7T MRI data acquired from healthy volunteers and patients with multiple sclerosis by comparing against manual segmentations delineated by a neuroradiologist, guided by the Morel atlas. Segmentation accuracy was very high, with uniformly high Dice indices: at least 0.85 for large nuclei like the pulvinar and mediodorsal nuclei and at least 0.7 even for small structures such as the habenular, centromedian, and lateral and medial geniculate nuclei. Volume similarity indices ranged from 0.82 for the smaller nuclei to 0.97 for the larger nuclei. Volumetry revealed that the volumes of the right anteroventral, right ventral posterior lateral, and both right and left pulvinar nuclei were significantly lower in MS patients compared to controls, after adjusting for age, sex and intracranial volume. Lastly, we evaluated the potential of this method for targeting the Vim nucleus for deep brain surgery and focused ultrasound thalamotomy by overlaying the Vim nucleus segmented from pre-operative data on post-operative data. The locations of the ablated region and active DBS contact corresponded well with the segmented Vim nucleus. Our fast, direct structural MRI based segmentation method opens the door for MRI guided intra-operative procedures like thalamotomy and asleep DBS electrode placement as well as for accurate quantification of thalamic nuclear volumes to follow progression of neurological disorders.

Comparison of MRI-based automated segmentation methods and functional neurosurgery targeting with direct visualization of the Ventro-intermediate thalamic nucleus at 7T

The ventro-intermediate nucleus (Vim), as part of the motor thalamic nuclei, is a commonly used target in functional stereotactic neurosurgery for treatment of drug-resistant tremor. As it cannot be directly visualized on routinely used magnetic resonance imaging (MRI), its clinical targeting is performed using indirect methods. Recent literature suggests that the Vim can be directly visualized on susceptibility-weighted imaging (SWI) acquired at 7 T. Our work aims to assess the distinguishable Vim on 7 T SWI in both healthy-population and patients and, using it as a reference, to compare it with: (1) The clinical targeting, (2) The automated parcellation of thalamic subparts based on 3 T diffusion MRI (dMRI), and (3) The multi-atlas segmentation techniques. In 95.2% of the data, the manual outline was adjacent to the inferior lateral border of the dMRI-based motor-nuclei group, while in 77.8% of the involved cases, its ventral part enclosed the Guiot points. Moreover, the late MRI signature in the patients was always observed in the anterior part of the manual delineation and it overlapped with the multi-atlas outline. Overall, our study provides new insight on Vim discrimination through MRI and imply novel strategies for its automated segmentation, thereby opening new perspectives for standardizing the clinical targeting.

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.

In-vivo probabilistic atlas of human thalamic nuclei based on diffusion- weighted magnetic resonance imaging

The thalamic nuclei are involved in many neurodegenerative diseases and therefore, their identification is of key importance in numerous clinical treatments. Automated segmentation of thalamic subparts is currently achieved by exploring diffusion-weighted magnetic resonance imaging (DW-MRI), but in absence of such data, atlas-based segmentation can be used as an alternative. Currently, there is a limited number of available digital atlases of the thalamus. Moreover, all atlases are created using a few subjects only, thus are prone to errors due to the inter-subject variability of the thalamic morphology. In this work, we present a probabilistic atlas of anatomical subparts of the thalamus built upon a relatively large dataset where the individual thalamic parcellation was done by employing a recently proposed automatic diffusion-based clustering method. Our analyses, comparing the segmentation performance between the atlas-based and the clustering method, demonstrate the ability of the provided atlas to substitute the automated diffusion-based subdivision in the individual space when the DW-MRI is not available.

Human Motor Thalamus Reconstructed in 3D from Continuous Sagittal Sections with Identified Subcortical Afferent Territories

Classification and delineation of the motor-related nuclei in the human thalamus have been the focus of numerous discussions for a long time. Difficulties in finding consensus have for the most part been caused by paucity of direct experimental data on connections of individual nuclear entities. Kultas-Ilinsky et al. (2011) showed that distribution of glutamic acid decarboxylase isoform 65 (GAD65), the enzyme that synthesizes inhibitory neurotransmitter γ-aminobutyric acid, is a reliable marker that allows to delineate connectionally distinct nuclei in the human motor thalamus, namely the territories innervated by nigral, pallidal, and cerebellar afferents. We compared those immunocytochemical staining patterns with underlying cytoarchitecture and used the latter to outline the three afferent territories in a continuous series of sagittal Nissl-stained sections of the human thalamus. The 3D volume reconstructed from the outlines was placed in the Talairach stereotactic coordinate system relative to the intercommissural line and sectioned in three stereotactic planes to produce color-coded nuclear maps. This 3D coordinate-based atlas was coregistered to the Montreal Neurological Institute (MNI-152) space. The current report proposes a simplified nomenclature of the motor-related thalamic nuclei, presents images of selected histological sections and stereotactic maps illustrating topographic relationships of these nuclei as well as their relationship with adjacent somatosensory afferent region. The data are useful in different applications such as functional MRI and diffusion tractography. The 3D dataset is publicly available under an open license and can also be applicable in clinical interventions in the thalamus.

Functional Connectivity Parcellation of the Human Thalamus by Independent Component Analysis

As a key structure to relay and integrate information, the thalamus supports multiple cognitive and affective functions through the connectivity between its subnuclei and cortical and subcortical regions. Although extant studies have largely described thalamic regional functions in anatomical terms, evidence accumulates to suggest a more complex picture of subareal activities and connectivities of the thalamus. In this study, we aimed to parcellate the thalamus and examine whole-brain connectivity of its functional clusters. With resting state functional magnetic resonance imaging data from 96 adults, we used independent component analysis (ICA) to parcellate the thalamus into 10 components. On the basis of the independence assumption, ICA helps to identify how subclusters overlap spatially. Whole brain functional connectivity of each subdivision was computed for independent component's time course (ICtc), which is a unique time series to represent an IC. For comparison, we computed seed-region-based functional connectivity using the averaged time course across all voxels within a thalamic subdivision. The results showed that, at p < 10-6, corrected, 49% of voxels on average overlapped among subdivisions. Compared with seed-region analysis, ICtc analysis revealed patterns of connectivity that were more distinguished between thalamic clusters. ICtc analysis demonstrated thalamic connectivity to the primary motor cortex, which has eluded the analysis as well as previous studies based on averaged time series, and clarified thalamic connectivity to the hippocampus, caudate nucleus, and precuneus. The new findings elucidate functional organization of the thalamus and suggest that ICA clustering in combination with ICtc rather than seed-region analysis better distinguishes whole-brain connectivities among functional clusters of a brain region.

Robust thalamic nuclei segmentation method based on local diffusion magnetic resonance properties

The thalamus is an essential relay station in the cortical–subcortical connections. It is characterized by a complex anatomical architecture composed of numerous small nuclei, which mediate the involvement of the thalamus in a wide range of neurological functions. We present a novel framework for segmenting the thalamic nuclei, which explores the orientation distribution functions (ODFs) from diffusion magnetic resonance images at 3 T. The differentiation of the complex intra-thalamic microstructure is improved by using the spherical harmonic (SH) representation of the ODFs, which provides full angular characterization of the diffusion process in each voxel. The clustering was performed using the k-means algorithm initialized in a data-driven manner. The method was tested on 35 healthy volunteers and our results show a robust, reproducible and accurate segmentation of the thalamus in seven nuclei groups. Six of them closely matched the anatomy and were labeled as anterior, ventral anterior, medio-dorsal, ventral latero-ventral, ventral latero-dorsal and pulvinar, while the seventh cluster included the centro-lateral and the latero-posterior nuclei. Results were evaluated both qualitatively, by comparing the segmented nuclei to the histological atlas of Morel, and quantitatively, by measuring the clusters’ extent and the clusters’ spatial distribution across subjects and hemispheres. We also showed the robustness of our approach across different sequences and scanners, as well as intra-subject reproducibility of the segmented clusters using additional two scan–rescan datasets. We also observed an overlap between the path of the main long-connection tracts passing through the thalamus and the spatial distribution of the nuclei identified with our clustering algorithm. Our approach, based on SH representations of the ODFs, outperforms the one based on angular differences between the principle diffusion directions, which is considered so far as state-of-the-art method. Our findings show an anatomically reliable segmentation of the main groups of thalamic nuclei that could be of potential use in many clinical applications.

Variations in Thalamic Anatomy Affect Targeting in Deep Brain Stimulation for Epilepsy

BACKGROUND

Thalamic size and shape vary significantly across patients - with changes specific to the anterior thalamus occurring with age and in the setting of chronic epilepsy. Such ambiguity raises concerns regarding electrode position and potential implications for seizure outcomes.

METHODS

MRIs from 6 patients from a single center underwent quantitative analysis. In addition to direct measurements from postimplantation MRIs, the CRAnialVault Explorer suite was used to normalize electrode position to a common reference system. Relationships between thalamic dimensions, electrode location, and seizure outcome were analyzed.

RESULTS

Although this study group was too small to sufficiently power statistical analysis, general trends were identified. There was a trend towards smaller thalamic volumes in nonresponders. Electrode locations demonstrated more variation after normalization. There was a trend towards a more lateral, posterior, and inferior electrode position in nonresponders.

CONCLUSIONS

Variations in thalamic shape and volume necessitate direct targeting. Given that changes occur to thalamic anatomy with age and in the setting of epilepsy, improved methods for visualizing and targeting the anterior nucleus are necessary. Pronounced thalamic atrophy may preclude proper electrode placement and serve as a poor prognostic indicator. A greater understanding of thalamic anatomy and connectivity is necessary to optimize deep brain stimulation for epilepsy.

Defining thalamic nuclei and topographic connectivity gradients in vivo

The thalamus consists of multiple nuclei that have been previously defined by their chemoarchitectual and cytoarchitectual properties ex vivo. These form discrete, functionally specialized, territories with topographically arranged graduated patterns of connectivity. However, previous in vivo thalamic parcellation with MRI has been hindered by substantial inter-individual variability or discrepancies between MRI derived segmentations and histological sections. Here, we use the Euclidean distance to characterize probabilistic tractography distributions derived from diffusion MRI. We generate 12 feature maps by performing voxel-wise parameterization of the distance histograms (6 feature maps) and the distribution of three-dimensional distance transition gradients generated by applying a Sobel kernel to the distance metrics. We use these 12 feature maps to delineate individual thalamic nuclei, then extract the tractography profiles for each and calculate the voxel-wise tractography gradients. Within each thalamic nucleus, the tractography gradients were topographically arranged as distinct non-overlapping cortical networks with transitory overlapping mid-zones. This work significantly advances quantitative segmentation of the thalamus in vivo using 3T MRI. At an individual subject level, the thalamic segmentations consistently achieve a close relationship with a priori histological atlas information, and resolve in vivo topographic gradients within each thalamic nucleus for the first time. Additionally, these techniques allow individual thalamic nuclei to be closely aligned across large populations and generate measures of inter-individual variability that can be used to study both basic function and pathological processes in vivo.

Disrupted Working Memory Circuitry in Schizophrenia: Disentangling fMRI Markers of Core Pathology vs Other Aspects of Impaired Performance

Working memory (WM) impairment, a core feature of schizophrenia, is often associated with aberrant dorsolateral prefrontal cortex (dlPFC) activation. Reduced resting-state connectivity within the frontoparietal control network (FPCN) has also been reported in schizophrenia. However, interpretation of WM-related dlPFC dysfunction has been limited by performance differences between patients and controls, and by uncertainty over the relevance of resting-state connectivity to network engagement during task. We contrasted brain activation in 40 schizophrenia patients and 40 controls during verbal WM performance, and evaluated underlying functional connectivity during rest and task. During correct trials, patients demonstrated normal FPCN activation, despite an inverse relationship between positive symptoms and activation. FPCN activation differed between the groups only during error trials (controls>patients). In contrast, controls demonstrated stronger deactivation of the ventromedial prefrontal cortex (vmPFC) during correct and error trials. Functional connectivity analysis indicated impaired resting-state FPCN connectivity in patients, but normal connectivity during task. However, patients showed abnormal connectivity among regions such as vmPFC, lateral orbitofrontal cortex, and parahippocampal gyrus (PHG) during both rest and task. During task, patients also exhibited altered thalamic connectivity to PHG and FPCN. Activation and connectivity patterns that were more characteristic of controls generally correlated with better performance. In summary, patients demonstrated normal FPCN activation when they remained on-task, and exhibited normal FPCN connectivity during WM, whereas vmPFC deactivation differences persisted regardless of WM performance. Our findings suggest that altered FPCN activation in patients reflects performance difference, and that limbic and thalamic dysfunction is critically involved in WM deficits in schizophrenia.

Semi-Automatic Segmentation of Optic Radiations and LGN, and Their Relationship to EEG Alpha Waves

At rest, healthy human brain activity is characterized by large electroencephalography (EEG) fluctuations in the 8-13 Hz range, commonly referred to as the alpha band. Although it is well known that EEG alpha activity varies across individuals, few studies have investigated how this may be related to underlying morphological variations in brain structure. Specifically, it is generally believed that the lateral geniculate nucleus (LGN) and its efferent fibres (optic radiation, OR) play a key role in alpha activity, yet it is unclear whether their shape or size variations contribute to its inter-subject variability. Given the widespread use of EEG alpha in basic and clinical research, addressing this is important, though difficult given the problems associated with reliably segmenting the LGN and OR. For this, we employed a multi-modal approach and combined diffusion magnetic resonance imaging (dMRI), functional magnetic resonance imaging (fMRI) and EEG in 20 healthy subjects to measure structure and function, respectively. For the former, we developed a new, semi-automated approach for segmenting the OR and LGN, from which we extracted several structural metrics such as volume, position and diffusivity. Although these measures corresponded well with known morphology based on previous post-mortem studies, we nonetheless found that their inter-subject variability was not significantly correlated to alpha power or peak frequency (p >0.05). Our results therefore suggest that alpha variability may be mediated by an alternative structural source and our proposed methodology may in general help in better understanding the influence of anatomy on function such as measured by EEG or fMRI.

Direct diffusion-based parcellation of the human thalamus

To assess stable anatomical features of the human thalamus, an unbiased diffusion tensor parcellation approach was used to segment thalamic substructures with similar spatial orientation. We determined localization, size and individual variations of 21 thalamic clusters in a group of 63 healthy human subjects (32 males/31 females). The laterality differences accounted for ± 6% and gender differences for ± 4% of the thalamic volume. Consecutively, five stable clusters in the anterior, medial, lateral and posterior thalamus were selected, which were common to 90% of all subjects and contained at least 10 voxels. These clusters could be assigned to the anteroventral nucleus (AN) group, the mediodorsal (MD) nucleus, the medial pulvinar (PuM), and the lateral nuclei group. The subcortical and cortical connectivity of these clusters revealed that: (1) the oblique cranio-caudal-oriented fibers of the AN cluster mainly connect to limbic structures, (2) the numerous dorso-frontal-oriented fibers of MD mainly project to the prefrontal cortex and the medial temporal lobe, (3) the fibers of the PuM running in parallel with the x-axis project to medio-occipital and medio-temporal areas and connect visual areas with the hippocampus and amygdala and via intrathalamic pathways with medio-frontal areas, and (4) the oblique caudo-cranial fibers of the two lateral clusters located anteriorly in the motor and posteriorly in the sensory thalamus are routing sensory-motor information from the brain stem via the internal capsule to pre- and peri-central regions of the cortex.

Contralateral cerebello-thalamo-cortical pathways with prominent involvement of associative areas in humans in vivo

In addition to motor functions, it has become clear that in humans the cerebellum plays a significant role in cognition too, through connections with associative areas in the cerebral cortex. Classical anatomy indicates that neo-cerebellar regions are connected with the contralateral cerebral cortex through the dentate nucleus, superior cerebellar peduncle, red nucleus and ventrolateral anterior nucleus of the thalamus. The anatomical existence of these connections has been demonstrated using virus retrograde transport techniques in monkeys and rats ex vivo. In this study, using advanced diffusion MRI tractography we show that it is possible to calculate streamlines to reconstruct the pathway connecting the cerebellar cortex with contralateral cerebral cortex in humans in vivo. Corresponding areas of the cerebellar and cerebral cortex encompassed similar proportion (about 80%) of the tract, suggesting that the majority of streamlines passing through the superior cerebellar peduncle connect the cerebellar hemispheres through the ventrolateral thalamus with contralateral associative areas. This result demonstrates that this kind of tractography is a useful tool to map connections between the cerebellum and the cerebral cortex and moreover could be used to support specific theories about the abnormal communication along these pathways in cognitive dysfunctions in pathologies ranging from dyslexia to autism.

In vivo characterization of the connectivity and subcomponents of the human globus pallidus

Projections from the substantia nigra and striatum traverse through the pallidum on the way to their targets. To date, in vivo characterization of these pathways remains elusive. Here we used high angular resolution diffusion imaging (N=138) to study the characteristics and structural subcompartments of the human pallidum. Our central result shows that the diffusion orientation distribution functions within the pallidum are asymmetrically oriented in a dorsal to dorsolateral direction, consistent with the orientation of underlying fiber systems. We also observed systematic differences in the diffusion signal between the two pallidal segments. Compared to the outer pallidal segment, the internal segment has more peaks in the diffusion orientation distribution and stronger anisotropy in the primary fiber direction, consistent with known cellular differences between the underlying nuclei. These differences in orientation, complexity, and degree of anisotropy are sufficiently robust to automatically segment the pallidal nuclei using diffusion properties. We characterize these patterns in one data set using diffusion spectrum imaging and replicate in a separate sample of subjects imaged using multi-shell imaging, highlighting the reliability of these diffusion patterns within pallidal nuclei. Thus the gray matter diffusion signal can be useful as an in vivo measure of the collective efferent pathways running through the human pallidum.

Ultra-high field magnetic resonance imaging of the basal ganglia and related structures

Deep brain stimulation is a treatment for Parkinson's disease and other related disorders, involving the surgical placement of electrodes in the deeply situated basal ganglia or thalamic structures. Good clinical outcome requires accurate targeting. However, due to limited visibility of the target structures on routine clinical MR images, direct targeting of structures can be challenging. Non-clinical MR scanners with ultra-high magnetic field (7T or higher) have the potential to improve the quality of these images. This technology report provides an overview of the current possibilities of visualizing deep brain stimulation targets and their related structures with the aid of ultra-high field MRI. Reviewed studies showed improved resolution, contrast- and signal-to-noise ratios at ultra-high field. Sequences sensitive to magnetic susceptibility such as T2^* and susceptibility weighted imaging and their maps in general showed the best visualization of target structures, including a separation between the subthalamic nucleus and the substantia nigra, the lamina pallidi medialis and lamina pallidi incompleta within the globus pallidus and substructures of the thalamus, including the ventral intermediate nucleus (Vim). This shows that the visibility, identification, and even subdivision of the small deep brain stimulation targets benefit from increased field strength. Although ultra-high field MR imaging is associated with increased risk of geometrical distortions, it has been shown that these distortions can be avoided or corrected to the extent where the effects are limited. The availability of ultra-high field MR scanners for humans seems to provide opportunities for a more accurate targeting for deep brain stimulation in patients with Parkinson's disease and related disorders.

Adaptive distance metric learning for diffusion tensor image segmentation

High quality segmentation of diffusion tensor images (DTI) is of key interest in biomedical research and clinical application. In previous studies, most efforts have been made to construct predefined metrics for different DTI segmentation tasks. These methods require adequate prior knowledge and tuning parameters. To overcome these disadvantages, we proposed to automatically learn an adaptive distance metric by a graph based semi-supervised learning model for DTI segmentation. An original discriminative distance vector was first formulated by combining both geometry and orientation distances derived from diffusion tensors. The kernel metric over the original distance and labels of all voxels were then simultaneously optimized in a graph based semi-supervised learning approach. Finally, the optimization task was efficiently solved with an iterative gradient descent method to achieve the optimal solution. With our approach, an adaptive distance metric could be available for each specific segmentation task. Experiments on synthetic and real brain DTI datasets were performed to demonstrate the effectiveness and robustness of the proposed distance metric learning approach. The performance of our approach was compared with three classical metrics in the graph based semi-supervised learning framework.

Visualization of intra-thalamic nuclei with optimized white-matter-nulled MPRAGE at 7T

Novel MR image acquisition strategies have been investigated to elicit contrast within the thalamus, but direct visualization of individual thalamic nuclei remains a challenge because of their small size and the low intrinsic contrast between adjacent nuclei. We present a step-by-step specific optimization of the 3D MPRAGE pulse sequence at 7T to visualize the intra-thalamic nuclei. We first measured T1 values within different sub-regions of the thalamus at 7T in 5 individuals. We used these to perform simulations and sequential experimental measurements (n=17) to tune the parameters of the MPRAGE sequence. The optimal set of parameters was used to collect high-quality data in 6 additional volunteers. Delineation of thalamic nuclei was performed twice by one rater and MR-defined nuclei were compared to the classic Morel histological atlas. T1 values within the thalamus ranged from 1400ms to 1800ms for adjacent nuclei. Using these values for theoretical evaluations combined with in vivo measurements, we showed that a short inversion time (TI) close to the white matter null regime (TI=670ms) enhanced the contrast between the thalamus and the surrounding tissues, and best revealed intra-thalamic contrast. At this particular nulling regime, lengthening the time between successive inversion pulses (TS=6000ms) increased the thalamic signal and contrast and lengthening the α pulse train time (N*TR) further increased the thalamic signal. Finally, a low flip angle during the gradient echo acquisition (α=4°) was observed to mitigate the blur induced by the evolution of the magnetization along the α pulse train. This optimized set of parameters enabled the 3D delineation of 15 substructures in all 6 individuals; these substructures corresponded well with the known anatomical structures of the thalamus based on the classic Morel atlas. The mean Euclidean distance between the centers of mass of MR- and Morel atlas-defined nuclei was 2.67mm (±1.02mm). The reproducibility of the MR-defined nuclei was excellent with intraclass correlation coefficient measured at 0.997 and a mean Euclidean distance between corresponding centers of mass found at first versus second readings of 0.69mm (±0.38mm). This 7T strategy paves the way to better identification of thalamic nuclei for neurosurgical planning and investigation of regional changes in neurological disorders.