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

A roadmap towards standardized neuroimaging approaches for human thalamic nuclei

The thalamus has a key role in mediating cortical-subcortical interactions but is often neglected in neuroimaging studies, which mostly focus on changes in cortical structure and activity. One of the main reasons for the thalamus being overlooked is that the delineation of individual thalamic nuclei via neuroimaging remains controversial. Indeed, neuroimaging atlases vary substantially regarding which thalamic nuclei are included and how their delineations were established. Here, we review current and emerging methods for thalamic nuclei segmentation in neuroimaging data and consider the limitations of existing techniques in terms of their research and clinical applicability. We address these challenges by proposing a roadmap to improve thalamic nuclei segmentation in human neuroimaging and, in turn, harmonize research approaches and advance clinical applications. We believe that a collective effort is required to achieve this. We hope that this will ultimately lead to the thalamic nuclei being regarded as key brain regions in their own right and not (as often currently assumed) as simply a gateway between cortical and subcortical regions.

Hippocampal subfields and thalamic nuclei associations with clinical outcomes in multiple sclerosis: An ultrahigh field MRI study

BACKGROUND

Previous studies have shown that thalamic and hippocampal neurodegeneration is associated with clinical decline in Multiple Sclerosis (MS). However, contributions of the specific thalamic nuclei and hippocampal subfields require further examination.

OBJECTIVE

Using 7 Tesla (7T) magnetic resonance imaging (MRI), we investigated the cross-sectional associations between functionally grouped thalamic nuclei and hippocampal subfields volumes and T1 relaxation times (T1-RT) and subsequent clinical outcomes in MS.

METHODS

High-resolution T1-weighted and T2-weighted images were acquired at 7T (n=31), preprocessed, and segmented using the Thalamus Optimized Multi Atlas Segmentation (THOMAS, for thalamic nuclei) and the Automatic Segmentation of Hippocampal Subfields (ASHS, for hippocampal subfields) packages. We calculated Pearson correlations between hippocampal subfields and thalamic nuclei volumes and T1-RT and subsequent multi-modal rater-determined and patient-reported clinical outcomes (∼2.5 years after imaging acquisition), correcting for confounders and multiple tests.

RESULTS

Smaller volume bilaterally in the anterior thalamus region correlated with worse performance in gait function, as measured by the Patient Determined Disease Steps (PDDS). Additionally, larger volume in most functional groups of thalamic nuclei correlated with better visual information processing and cognitive function, as measured by the Symbol Digit Modalities Test (SDMT). In bilateral medial and left posterior thalamic regions, there was an inverse association between volumes and T1-RT, potentially indicating higher tissue degeneration in these regions. We also observed marginal associations between the right hippocampal subfields (both volumes and T1-RT) and subsequent clinical outcomes, though they did not survive correction for multiple testing.

CONCLUSION

Ultrahigh field MRI identified markers of structural damage in the thalamic nuclei associated with subsequently worse clinical outcomes in individuals with MS. Longitudinal studies will enable better understanding of the role of microstructural integrity in these brain regions in influencing MS outcomes.

Mapping hippocampal and thalamic atrophy in epilepsy: A 7-T magnetic resonance imaging study

OBJECTIVE

Epilepsy patients are often grouped together by clinical variables. Quantitative neuroimaging metrics can provide a data-driven alternative for grouping of patients. In this work, we leverage ultra-high-field 7-T structural magnetic resonance imaging (MRI) to characterize volumetric atrophy patterns across hippocampal subfields and thalamic nuclei in drug-resistant focal epilepsy.

METHODS

Forty-two drug-resistant epilepsy patients and 13 controls with 7-T structural neuroimaging were included in this study. We measured hippocampal subfield and thalamic nuclei volumetry, and applied an unsupervised machine learning algorithm, Latent Dirichlet Allocation (LDA), to estimate atrophy patterns across the hippocampal subfields and thalamic nuclei of patients. We studied the association between predefined clinical groups and the estimated atrophy patterns. Additionally, we used hierarchical clustering on the LDA factors to group patients in a data-driven approach.

RESULTS

In patients with mesial temporal sclerosis (MTS), we found a significant decrease in volume across all ipsilateral hippocampal subfields (false discovery rate-corrected p [pFDR] < .01) as well as in some ipsilateral (pFDR < .05) and contralateral (pFDR < .01) thalamic nuclei. In left temporal lobe epilepsy (L-TLE) we saw ipsilateral hippocampal and some bilateral thalamic atrophy (pFDR < .05), whereas in right temporal lobe epilepsy (R-TLE) extensive bilateral hippocampal and thalamic atrophy was observed (pFDR < .05). Atrophy factors demonstrated that our MTS cohort had two atrophy phenotypes: one that affected the ipsilateral hippocampus and one that affected the ipsilateral hippocampus and bilateral anterior thalamus. Atrophy factors demonstrated posterior thalamic atrophy in R-TLE, whereas an anterior thalamic atrophy pattern was more common in L-TLE. Finally, hierarchical clustering of atrophy patterns recapitulated clusters with homogeneous clinical properties.

SIGNIFICANCE

Leveraging 7-T MRI, we demonstrate widespread hippocampal and thalamic atrophy in epilepsy. Through unsupervised machine learning, we demonstrate patterns of volumetric atrophy that vary depending on disease subtype. Incorporating these atrophy patterns into clinical practice could help better stratify patients to surgical treatments and specific device implantation strategies.

The effects of P2 segment of posterior cerebral artery to thalamus blood supply pattern on gait in cerebral small vessel disease: A 7 T MRI based study

Gait disturbance is a manifestation of cerebral small vessel disease (CSVD). The posterolateral thalamus (PL), whose blood is mainly supplied by the P2 segment of posterior cerebral artery (P2-PCA), plays pivotal roles in gait regulation. We investigated the influence of the distance between P2-PCA and PL on gait with varying CSVD burden. 71 participants were divided into low and high CSVD burden groups. The distance from P2-PCA to PL was measured using 7 T TOF-MRA and categorized into an immediate or distant PCA-to-thalamus pattern. Functional connectivity (FC) and voxel-based morphometry were assessed to evaluate functional and structural alterations. In the low CSVD burden group, immediate PCA-to-thalamus supply strongly correlates with longer step length and higher wave phase time percent, and exhibited enhanced FCs in left supplementary motor area, right precentral cortex (PreCG.R). While in the high CSVD burden group, no association between PCA-to-thalamus pattern and gait was found, and we observed reduced FC in PreCG.R with immediate PCA-to-thalamus pattern. Higher CSVD burden was associated with decreased gray matter density in bilateral thalamus. However, no significant structural thalamic change was observed between the two types of PCA-to-thalamus patterns in all patients. Our study demonstrated patients with immediate PCA-to-thalamus supply exhibited better gait performance in low CSVD burden populations, which also correlated with enhanced FCs in motor-related cortex, indicating the beneficial effects of the immediate PCA-to-thalamus supply pattern. In the higher burden CSVD populations, the effects of PCA-to-thalamus pattern on gait are void, attributable to the CSVD-related thalamic destruction and impairment of thalamus-related FC.

Associations of thalamocortical networks with reduced mindfulness in alcohol use disorder

Background

Increased mindfulness is associated with reduced alcohol consumption in patients with alcohol use disorder (AUD) after residential treatment. However, the underlying neurobiological mechanism of mindfulness in AUD is unclear. Therefore, we investigate the structural and functional alterations of the thalamocortical system with a focus on the mediodorsal thalamic nucleus (MD-TN), the default mode and the salience network (DMN/SN) which has previously been associated with mindfulness in healthy subjects. We hypothesized lower mindfulness and reduced structural and functional connectivity (FC) of the thalamocortical system, particularly in the DMN/SN in AUD. We assumed that identified neurobiological alterations in AUD are associated with impairments of mindfulness.

Methods

Forty-five abstinent patients with AUD during residential treatment and 20 healthy controls (HC) were recruited. Structural and resting-state functional MRI-scans were acquired. We analysed levels of mindfulness, thalamic volumes and network centrality degree of the MD-TN using multivariate statistics. Using seed-based whole brain analyses we investigated functional connectivity (FC) of the MD-TN. We performed exploratory correlational analyses of structural and functional DMN/SN measurements with levels of mindfulness.

Results

In AUD we found significantly lower levels of mindfulness, lower bilateral thalamic and left MD-TN volumes, reduced FC between MD-TN and anterior cingulum/insula and lower network centrality degree of the left MD-TN as compared to HC. In AUD, lower mindfulness was associated with various reductions of structural and functional aspects of the MD-TN.

Conclusion

Our results suggest that structural and functional alterations of a network including the MD-TN and the DMN/SN underlies disturbed mindfulness in AUD.

Optimization of MP2RAGE T1 mapping with radial view-ordering for deep brain stimulation targeting at 7 T MRI

PURPOSE

Deep brain stimulation (DBS) is an effective treatment of various neurological disorders. Due to higher intrinsic signal, 7 T MRI can potentially improve delineation of DBS targets. However, the severe RF transmit field (B1+) inhomogeneity at 7 T can compromise the image contrast of traditional single-contrast sequences for DBS targeting, leading to sub-optimal target visualization. The Magnetization Prepared 2 Rapid Acquisition Gradient Echo (MP2RAGE)-based T1 mapping provides an alternative to the traditional single-contrast techniques by allowing retrospective synthesis of images at arbitrary inversion times to aid in visualization of various DBS targets. With this approach, optimization of sequence parameters to create T1 maps with low noise and low quantification bias is critical, as these characteristics directly affect the noise and uniformity of the synthetic images. In this work, we perform sequence optimization for MP2RAGE-based T1 mapping using a radial view-ordering technique to improve image quality, and demonstrate the clinical utility of T1 mapping approach for DBS targeting.

METHODS

We first introduce a systematic sequence optimization framework for 7 T MP2RAGE T1 mapping by formulating it into a constrained, multi-dimensional optimization process considering the effect of B1+ inhomogeneity on image noise, T1 quantification bias, and image blurring. With this framework, we investigate the use of radial view-order approach for T1 mapping, in lieu of the conventional linear view-ordering. Bloch's equation-based simulations were performed to compare the T1 maps generated using different approaches. Images of healthy volunteer and patients were acquired on a clinical 7 T MRI scanner for validation and to demonstrate the utility of T1 mapping for DBS targeting.

RESULTS

Numerical experiments demonstrated that the proposed framework allowed optimization of image SNR in T1 maps while controlling the quantification bias and image blurring, therefore facilitating the selection of optimal sequence parameters for visualizing DBS targets. The optimized sequence using radial view-ordering offered 40-60% noise reduction compared to the linear view-ordering. The improvement of SNR was confirmed in the in vivo examples. Clinical images showed that the synthetic images generated from the optimized T1 maps allowed clear visualization of DBS targets.

CONCLUSION

We demonstrated the optimization of MP2RAGE T1 mapping with radial view-ordering technique for DBS targeting at 7 T and showed that the optimized sequence allows retrospective generation of synthetic inversion time images commonly utilized in DBS targeting, such as fast gray matter acquisition T1 inversion recovery (FGATIR) and edge-enhancing gradient echo (EDGE) sequences.

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.

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.

Automatic segmentation of Deep Grey Nuclei using a high-resolution 7T MRI Atlas - quantification of T1 values in healthy volunteers

We present a new consensus atlas of deep grey nuclei obtained by shape-based averaging of manual segmentation of two experienced neuroradiologists and optimized from 7T MP2RAGE images acquired at (0.6mm)³ in 60 healthy subjects. A group-wise normalization method was used to build a high-contrast and high-resolution T₁ -weighted brain template (0.5mm)³ using data from 30 out of the 60 controls. Delineation of 24 deep grey nuclei per hemisphere, including the claustrum and twelve thalamic nuclei, was then performed by two expert neuroradiologists and reviewed by a third neuroradiologist according to tissue contrast and external references based on the Morel atlas. Corresponding deep grey matter structures were also extracted from the Morel and CIT168 atlases. The data-derived, Morel and CIT168 atlases were all applied at the individual level using non-linear registration to fit the subject reference and to extract absolute mean quantitative T₁ values derived from the 3D-MP2RAGE volumes, after correction for residual B₁ ⁺ biases. Three metrics (The Dice and the volumetric similarity coefficients, and a novel Hausdorff distance) were used to estimate the inter-rater agreement of manual MRI segmentation and inter-atlas variability, and these metrics were measured to quantify biases due to image registration and their impact on the measurements of the quantitative T₁ values was highlighted. This represents a fully-automated segmentation process permitting the extraction of unbiased normative T₁ values in a population of young healthy controls as a reference for characterizing subtle structural alterations of deep grey nuclei relevant to a range of neurological diseases.

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.